Puntes Hidraulica Fluvial 2013-I PDF
Puntes Hidraulica Fluvial 2013-I PDF
Puntes Hidraulica Fluvial 2013-I PDF
1.1 INTRODUCCIÓN
Desde sus orígenes, el hombre ha dependido del agua no sólo como elemento vital, sino también como
promotor de su desarrollo. Las más importantes civilizaciones se asentaron en las riberas de grandes
ríos: Mesopotamia en una región muy fértil alimentada por los ríos Tigris y Eufrates; Egipto debe su vida
al Nilo; China, a Yangtze; la India, al Indo y al Ganges, las civilizaciones azteca e Inca en México y Perú,
respectivamente, entre otras culturas.
Tales civilizaciones comprendieron que el agua permitía la vida y el desarrollo; sin embargo también
observaron que junto con tales beneficios existían .importantes riesgos, como la destrucción de sus
parcelas debido a inundaciones causadas por precipitaciones intensas.
Dentro de este binomio beneficio-perjuicio, dada su particular forma de vida, eran muy superiores los
beneficios en relación con los daños. Asimismo, tenían conocimiento de la variación espacial y temporal
del recurso, lo que, aunado a la necesidad de contar con el agua y protegerse de ella, motivó el
planteamiento de importantes obras hidráulicas, cuyos vestigios sorprenden al advertir su concepción y
funcionalidad.
Cuatro mil años antes de Cristo, estos antiguos pobladores construían presas de almacenamiento,
canales para riego agrícola y acueductos, con los que hacían llegar el agua de la fuente a la ciudad para
usarla de manera doméstica.
En la antigüedad el control de los ríos ya tenía una importancia fundamental; su estudio y aplicación de
acertadas técnicas llevó a un ingeniero hidráulico a ser emperador de China en el año 2278 A.C. El
ingeniero-emperador llamado Yu hizo la regulación de nueve ríos de acuerdo con sus características y
propiedades particulares, entre los que destacan el Hwang Ho y el Yangtze Kiang. Estos ríos
permanecieron estables en sus cauces por casi 1,700 años.
La hidráulica fluvial fue iniciada por Guglielmini, quien en 1690 publicó su libro Aquarum fluentium
mensura nova methodo inquisita, que consistía en un método para medir el flujo de agua mediante una
esfera suspendida por la corriente. Las contribuciones de Guglielmini a la hidráulica fluvial se dieron
principalmente a partir de observaciones de campo.
La ingeniería de ríos comprende el control y utilización de los ríos en beneficio del hombre.
En un sentido amplio sus objetivos pueden incluir el manejo de ríos, diseño de canales, control de
avenidas, abastecimiento de agua, mejoramiento de la navegación, diseño de estructuras hidráulicas,
mitigación de riesgos y comportamiento ambiental.
Aunque el dominio de la ingeniería de los ríos está limitado a los cauces y la escala de interés de la
ingeniería es el tiempo, el comportamiento del río es afectado por el sistema fluvial en su conjunto,
siendo importante en este aspecto la escala geológica. Por tanto la ingeniería de los ríos debe estar
basada en la comprensión del sistema fluvial y de conceptos de geomorfología.
a) El sistema fluvial
Un río esta dentro del dominio del sistema fluvial, que consiste de la cuenca de drenaje y de los
reservorios, lagos u océanos ubicados aguas abajo. Schumrn (1977) ha dividido el sistema fluvial en tres
partes, como se muestra en la figura 1.1.
La parte superior o Zona 1, es la porción de cuenca donde se origina la mayor cantidad de agua y de
sedimentos. Quebradas pequeñas en esta zona se caracterizan por ser inestables y a menudo por
encontrarse trenzados, Debido a la inestabilidad de los cauces, el estudio de la geomorfología del río
solo puede ser analizado sobre la base de asunciones gruesas más no detalladas.
La parte media o Zona 2, es el tramo en el cual el río es más estable y donde su configuración está
mejor definida. Ríos grandes presentan tramos largos en esta Zona 2, pero esta zona puede estar
ausente en ríos pequeños. Este es el tramo en el cual se realizan los mayores estudios, modelamientos
y obras de control.
La Zona 3 está cerca de la boca de salida donde el río aluvial está bajo la influencia de las variaciones
de las mareas. Los ríos en esta zona a menudo se encuentran trenzados.
o Canales artificiales y ríos tienen en común que transporte agua en lámina libre. Sin embargo en un
canal artificial hay determinaciones previas que responden a las siguientes preguntas:
o En un río no hay determinaciones previas, las respuestas son estudio de la hidrología, geomorfología o
la hidráulica fluvial.
o Un canal artificial es prismático y definido por una sección tipo. Un río no.
o La rugosidad en un canal es un parámetro bien definido y determinante de su capacidad.
o En un río, el caudal y la altura de agua están relacionados de una manera más compleja.
Existe una resistencia al flujo por el tamaño del grano del material de fondo y otra añadida por las formas
del fondo granular (dunas, etc).
UNIDAD DIDACTICA Nº 2
2.1 RIOS
a) Conceptos Generales
Los ríos son cursos de aguas permanentes o temporarios que recorren la superficie terrestre siguiendo la
pendiente del suelo y desembocan en el mar, otros ríos, lagos o se pierden en el interior de los
continentes.
Se llama cuenca hidrográfica de un río a la superficie terrestre que recorren todos los ríos y afluentes
que desaguan un mismo río principal.
El origen de los mismos se debe a la acumulación de agua de lluvia o al deshielo de las cumbres
montañosas. De acuerdo a su origen, el caudal de los mismos, no es igual en todas las épocas del año,
pudiendo presentar variaciones en su régimen.
Elementos de un río
Régimen de un río
Se llama caudal de un río a la cantidad de agua que transportan por unidad de tiempo. De acuerdo a su
origen, pueden distinguirse tres regímenes diferentes:
El período de menor caudal recibe la denominación de estiaje, mientras que el de mayor caudal se
denomina creciente. Existen ríos que presentan regímenes regulares y mantienen su nivel de caudal
durante todo el año.
CAUDAL MEDIO
- Diario: Volumen de agua que pasa a través de una sección transversal del río durante el día dividido
por el número de segundos del día.
- Mensual: Media aritmética de los caudales medios diarios del mes
- Anual: Media aritmética de los caudales medios diarios del año
MODULO ANUAL:
Es la media aritmética de los caudales medios anuales observados en un período de tiempo
determinado, mínimo 10 años.
Qr 5 Baja
5 Qr 15 Media
15 Qr 40 Alta
Qr 40 Muy alta
IRREGULARIDAD DE LOS RIOS: Variaciones entre los valores extremos del caudal
Se expresa por medio del coeficiente de irregularidad, que se obtiene dividiendo el caudal máximo por el
mínimo.
Q m áx
Ci
Qmín
Ci 3 Régimen irregular
3 Ci 7 Régimen de transición
Ci 7 Régimen irregular
b) Morfología de un Río
La predicción de los cambios de forma de un río es importante por muchas razones. Por ejemplo, para la
selección de la ubicación de un puente o de una obra de toma es muy importante conocer cómo podría
ser el comportamiento del río en el futuro, de manera que se puedan proteger las estructuras. Es
importante predecir las formas de migración de un río en el caso de estructuras ya existentes o
poblaciones cercanas. Es necesario también evaluar la respuesta de un río a diferentes métodos de
protección de márgenes. En un río que migra es valioso evaluar la erosión de las márgenes y el
desarrollo del meandreo. Es importante estimar el ancho total ocupado por un río meandreante. En el
caso de la construcción de presas de almacenamiento o de derivación es importante evaluar el cambio
de régimen del río y, a su vez, el meandreo aguas abajo de la presa.
Las características morfológicas de un río desde el punto de vista de las características geométricas,
incluyen una revisión del concepto de:
o Régimen
o Descarga formativa
o Perfil longitudinal del río
o Clasificación de los ríos
o Geometría hidráulica
o Formación de meandros
o El análisis geomorfológico de la respuesta de los ríos.
c) Concepto de Régimen
El concepto de régimen generalmente es considerado como sinónimo de equilibrio. El concepto se
originó del estudio de canales aluviales estables, que con fondo móvil y bancos de tierra, no presentaban
socavación ni sedimentación durante un ciclo operativo. Sin embargo un canal aluvial utilizado para
irrigación, es usualmente operado bajo una descarga constante.
Debido a la variación natural de las descargas, el régimen real o dinámico de un río natural puede nunca
ser alcanzado, aunque cada río por su cuenta se ajusta constantemente hacia esa dirección.
Mackin (1948) definió el término "gradual" como una condición de equilibrio en ríos que actúan como
agente de transporte. Definió río graduado, como aquel que en un periodo de años, la pendiente se
ajusta suavemente para proveer, con la descarga disponible y las condiciones prevalentes de la
geometría del canal, la velocidad requerida para transportar la carga de sedimentos proporcionada por la
cuenca.
Esta simplificación aproximada se justifica en el hecho de que descargas bajas, con movimiento mínimo
de sedimentos, contribuyen menos a la formación del cauce. También, el incremento en la descarga por
encima del nivel de los bancos llenos, es mayormente absorbido por el ancho de la planicie de
inundación y por lo tanto generalmente tiene menos efecto en la forma del canal.
1. Método hidrológico
Los análisis hidrológicos basados en el régimen hidrológico de avenidas consisten simplemente en el
criterio de atribuir al caudal dominante un periodo de retorno (T) comprendido entre 1.5 y 7 años para el
caso de una hidrología irregular, a más irregularidad en la hidrología más alto es el periodo.
El trabajo puede considerarse el más completo y afinado en cuanto a las leyes de frecuencia de los
caudales máximos del río en estudio, ya que, además de usar una serie de “n” valores de caudales
instantáneos máximos anuales, incorpora al cálculo datos históricos. Se calculan caudales máximos, en
lo que se refiere a los bajos periodos de retorno entre 1.5 y 7 años, es decir Q T=2 años, QT=5 años y
QT= 7 años.
2. Método hidráulico
Para la aplicación del análisis hidráulico se toman algunas secciones características del río en estudio.
Este debe ser un lugar en el que el río muestra todavía hoy una morfología natural, donde exista una
gran barra de gravas y arenas en la parte interior de la curva y una zona profunda, cauce de las aguas
permanentes, en la parte exterior. Las secciones en esta región, consideramos poco alterada, son
preferibles a las secciones del río en otros lugares dentro del tramo de estudio, probablemente más
alteradas.
En las secciones transversales se determina el nivel de cauce lleno estudiando dónde el cociente entre
la anchura y la profundidad tiene un cambio fuerte de valor. El cauce determinado mediante esta
operación geométrica se identifica con el cauce principal en sentido geomorfológico modelado por el
caudal dominante. El valor de este caudal se obtiene utilizando la fórmula de Manning con un coeficiente
de rugosidad de 0.035 ó 0.040 (se justifica el uso de fórmulas de régimen uniforme en lugar de cálculos
en régimen gradualmente variado dado el carácter aproximado de esta determinación).
e) Tipos de Ríos
Diferentes factores influyen en una corriente para tomar una u otra forma, entre estos, los parámetros
hidráulicos, propiedades del fluido y características del flujo, características del material de fondo y de los
bordos (tamaño graduación, forma, etcétera), geometría y estructura de los bordes (altura, pendiente
cohesión, estratificación o tipo rocoso), características biológicas (tipo de vegetación tales como pasto,
HIDRAULICA FLUVIAL 6 ING. M. Sc. CESAR MILLA VERGARA
arbustos y árboles), factores humanos como agricultura, urbanización, drenaje, desarrollo de las llanuras
de inundación, y bordos de protección.
- Ríos rectos: Siguen una alineación recta. Existen en planicies que son inadecuadas para permitir
velocidades erosivas, o en pendientes pronunciadas donde se pueden alcanzar altas velocidades. En el
último caso es posible que la alineación se deba a la alta concentración de sedimento, la cual requiere
para su transporte toda la capacidad que proporciona una pendiente fuerte. Sin embargo, esta condición
no es un factor determinante, ya que bajo esas condiciones una corriente natural con gasto fluctuante
tiende a formar un cauce trenzado.
Ríos realmente rectos son raros en la naturaleza. Un río recto, en esta clasificación, se refiere como
aquel que no tiene comportamiento meándrico esto es, su sinuosidad es menor de 1.5.
- Ríos trenzados: Son aquellos formados por canales interconectados aleatoriamente, separados por
barras, que presentan la apariencia de una trenza. Este tipo de ríos se encuentra raramente en
1/4
pendientes relativamente fuertes, según Lane (1957), si S= 0.10 Q ó mayores, y según Leopold y
-O.44
Wolman (1957), si S= 0.06 Q o mayores. En el primer caso, Q es el gasto promedio en pies cúbicos
por segundo, y en el segundo caso es el gasto total contenido por el bordo. El término S es la pendiente
en pies por miles de pies. En ríos con formaciones de barras por el depósito de sedimentos puede ocurrir
el trenzado para diferentes pendientes, pero en ríos sedimentarios balanceados se cree que el trenzado
resulta por los depósitos de materiales transportados durante avenidas grandes en cantidades o
tamaños bastante considerables. Estos depósitos frecuentemente forman barras en donde florece la
vegetación.
- Ríos con meandro: Un río meándrico es aquel que posee una serie de curvas consecutivas, visto en
planta, forma de S.
En general se puede establecer, como punto de referencia, que un río se considera meándrico cuando
tiene su sinuosidad (S) mayor de 1.5, definiendo la sinuosidad como la relación existente entre la
longitud del cauce principal (Lt) y la longitud del valle (Lv) por el cual discurre dicho cauce (S=Lt/Lv).
Una de las clasificaciones más utilizadas por los geomorfólogos y ampliamente aceptadas en ingeniería
es la de jóvenes, maduros y viejos.
Los cauces se desarrollan primeramente mediante la erosión del agua sobre la superficie del terreno; por
ello, los ríos jóvenes tienen generalmente valles irregulares en forme de V. Sus cauces adoptan la
misma forma y están constituidos por materiales fracturados que pueden o no ser erosionables. Casi
todos los ríos de montaña y sus tributarios son ejemplos de ríos jóvenes.
En los ríos maduros, el valle fluvial se ensancha, las pendientes longitudinales son suaves y la erosión
lateral es mucho más significativa que la erosión en el fondo. El lecho del río alcanza una condición de
equilibrio, es decir, la pendiente y la energía del río son justamente suficientes para transportar el
material aportado al rio. En un cauce maduro hay llanuras de inundación angostas y se tiene al
desarrollo de los meandros.
Los cauces clasificados como ríos viejos, son extensiones en el tiempo de los maduros. A medida que
la erosión continúa, los valles fluviales se desarrollan hasta que sus características pasan a ser las de
mayor ancho y menor pendiente.
f) Ríos meandreantes
De acuerdo con Leopold y otros investigadores, el meandro consiste de un par de curvas opuestas. Las
pendientes de un río con meandreo son relativamente bajas en comparación con las pendientes de los
ríos rectos o trenzados; las pendientes suaves están relacionadas con velocidades bajas y números de
Froude pequeños.
Cuatro tipos principales de ríos se presentan en la figura siguiente. Estos pueden ser: ríos de canal
sinuoso, ríos sinuosos con barras, ríos sinuosos trenzados y ríos no sinuosos trenzados.
Los ríos de canal sinuoso tiene una pendiente plana, un ancho uniforme característico y moderada a alta
sinuosidad. El canal es relativamente angosto y profundo, con alto contenido de barro arcilloso y gran
estabilidad en los taludes.
Los ríos sinuosos trenzados, para la misma descarga, son más escarpados y anchos que los ríos
sinuosos con barras caracterizados por tasas rápidas de migración lateral y rápidos cambios en la
posición del talweg. Estos ríos tienen bastante transporte de material grueso en el fondo pero menos
contenido de barro arcilloso. Las barras son más irregulares a medida que el trenzado se incrementa.
Los ríos no sinuosos trenzados sin barras existen en pendientes escarpadas con transporte de material
grueso en el fondo y bajo contenido de barro arcilloso. Estos ríos son altamente trenzados y tienen
moderadas tasas de migración lateral en lugares donde las trenzas chocan contra los taludes.
La relación pendiente-descarga obtenida por Leopold y Wolman para el límite que separa ríos
meándricos y ríos escarpados trenzados es:
S 0.0125Q0.44
donde:
3
Q = caudal formativo en pies /s
Esta ecuación fue obtenida utilizando datos de ríos con fondo de arena y con fondos de grava.
Henderson (1961) probó refinar este criterio incluyendo el efecto del tamaño de la partícula del fondo
obteniendo la siguiente ecuación para ríos con fondo de grava.
S 0.0002d 1.15Q0.46
donde:
d = tamaño medio de la partícula en mm.
El conocimiento de la geomorfología es útil para el análisis cualitativo de la respuesta de los ríos. Una
relación muy útil para predecir la respuesta de un cauce fluvial fue la propuesta por Lane y la dispuso en
la analogía de la balanza, como se muestra en la figura.
qS q s D
donde:
q = caudal líquido
qs = caudal sólido de fondo
S = pendiente
D = tamaño del sedimento
Las características de los sedimentos reflejan los procesos de erosión, el transporte y el depósito del
material erosionado de la superficie de la cuenca. La densidad, el tamaño y la forma de las partículas
influyen determinantemente en las varias etapas del transporte de sedimentos que es de gran interés
para el diseño de las obras hidráulicas (como obras de protección en los ríos y para estimar los
volúmenes de azolve en vasos).
Es de interés el estudio de las propiedades de los sedimentos por su relación con las formas de fondo y
por consiguiente, con la resistencia y los regímenes de flujo. Asimismo, los procedimientos para calcular
el transporte de sedimentos en cauces utilizan diámetros representativos de las partículas del fondo
obtenidos de muestras del mismo.
Otro de los propósitos para determinar los parámetros de los sedimentos es comparar los análisis de las
muestras para correlacionar el tipo de sedimento con su ambiente. Además, la normalización de los
parámetros permite las comparaciones con resultados de otros estudios y, en algunos casos, inferir
condiciones a partir de los datos y resultados de otros cauces.
Las propiedades que se describen en este acápite son la densidad, el tamaño, la forma, la textura, la
composición mineral, la velocidad de caída y el ángulo de reposo, las cuales resultan ser las de mayor
influencia en los procesos de hidráulica fluvial y las más utilizadas para realizar el cálculo de los
parámetros de diseño en las obras de protección y control de cauces.
Por otra parte, el sedimento no está constituido por partículas idénticas, sino que éstas varían entre sí en
peso, tamaño y forma; por ello es necesario determinar o seleccionar parámetros que permitan
representar el comportamiento del conjunto de tamaños.
Con base en lo expuesto anteriormente, y de acuerdo con lo que señala la mayoría de los
investigadores, la determinación del tamaño de una partícula consiste en considerar algunas de las
dimensiones siguientes:
a) Diámetro de cribado: se define como la abertura mínima de la malla por la cual pasa la partícula.
b) Diámetro de sedimentación: es el diámetro de una esfera con la misma densidad y velocidad de caída
a la de una partícula al caer en un líquido a una temperatura constante.
c) Diámetro nominal: diámetro de una esfera con igual volumen que la partícula
d) Diámetro de caída estándar: diámetro de una esfera con densidad relativa de 2.65, que tiene la misma
velocidad de caída de la partícula cuando ambas se precipitan en agua destilada a 24° C.
e) Dimensiones triaxiales: longitud máxima y mínima, y la que resulte, medidas en la dirección de tres
ejes mutuamente perpendiculares.
Z
a Y
C a b c
X b
Z
Z
a
Y
Los diámetros de cribado y sedimentación son los que más se utilizan, sólo en ocasiones se usan las
dimensiones triaxiales. Cuando un material se ha cribado adecuadamente, sin forzar el paso de las
partículas del suelo a través de las mallas, el diámetro de cribado corresponde aproximadamente al de
sedimentación.
Las dimensiones triaxiales o la media de ellas se utiliza para fijar el tamaño de cantos rodados o boleos.
Por otra parte, el diámetro de cribado sirve para definir el tamaño de las arenas y gravas. Por último, el
diámetro de sedimentación se emplea para determinar el tamaño de partículas más finas como limos y
arcillas.
En hidráulica fluvial, la clasificación propuesta por la American Geophysical Unión es la más usada, la
cual se muestra en la siguiente tabla.
ii) Forma
Se refiere a la forma geométrica completa de una partícula sin considerar su tamaño o composición. Las
partículas de formas geométricas muy diferentes, pero con el mismo volumen y densidad pueden
comportarse en forma similar en un fluido, por lo que la forma puede definirse en términos de su
comportamiento dinámico.
La forma de las partículas puede ser muy variada y afecta su comportamiento dentro del agua. Las
partículas pueden tener formas tendientes a esferas, discos, láminas, elipsoides o ser completamente
irregulares, como se muestra en la siguiente figura:
bc
3
a2
Donde:
a, b y c son las dimensiones triaxiales de acuerdo con la definición siguiente:
Dnom
a
Donde:
Dnom = diámetro nominal
Con este valor y una vez medida a, b y c, puede aplicarse el diagrama de la siguiente figura, para la
clasificación de las partículas.
En la siguiente tabla, se presentan los valores de redondez de las partículas y una descripción de sus
características físicas.
Redondez Descripción SF
Partículas con superficie de fractura recientes,
1. Muy angular múltiples aristas y bordes cortantes 0.12 0.17
Presentan esquinas y bordes ásperos, lastiman, no son
2. Angular tan cortantes o afiladas sino de formas más bien 0.17 0.25
prismáticas
3. Subangular Se distinguen bordes y aristas, pero están
despuntadas o ligeramente redondeadas 0.25 0.35
4. Subredondeada Se distinguen, pero no se sienten las aristas ni los
bordes. Se aprecia que estos han sido redondeados 0.35 0.49
Tiende a equidimensional. No se diferencian esquinas
5. Redondeada o puntas al rotarlo con los dedos 0.49 0.70
6. Bien redondeadas Partículas con superficie llana o pulida.
Aproximadamente esféricas o elipsoidales 0.70 1.00
c) Velocidad de caída
Se define como la velocidad máxima que adquiere la partícula al caer dentro del agua; se alcanza
cuando el peso sumergido de la partícula se equilibra con el empuje del agua.
También se le llama velocidad terminal.
Su cálculo o medición es fundamental para obtener los parámetros que se utilizan en el análisis del
transporte de sedimentos y en algunas relaciones para predecir los regímenes de flujo.
Los criterios para el cálculo de la velocidad caída se basan en deducciones hechas para esferas y en
correcciones que han sido obtenidas de experimentos, según sea la forma de las partículas.
La velocidad de caída de una partícula esférica se determina cuando la fuerza de arrastre F D se equilibra
con el peso sumergido FS.
FD FS
C D 2 A
FD
2
S gD 3
FS
6
Donde:
CD : coeficiente de arrastre
: densidad del fluido
: velocidad de caída terminal de una esfera dentro de un fluido en reposo
A : área de la partícula esférica, D 2 / 4
S : densidad de la partícula
D : diámetro de la esfera
g : aceleración de la gravedad
Donde
1. Se supone un valor de CD
2. Se calcula la velocidad de caída con la ecuación anterior
3. Se obtienen el número de Reynolds de la partícula
4. Se determina el nuevo valor de CD
5. Se repiten los pasos hasta que el coeficiente de arrastre calculado es igual al supuesto
Además:
24 gD 2
Si para Re 1, entonces C D
Re 18
0.5
24 3 4g
Si para 0.1< Re 5, entonces C D 1 Re D
Re 16 3C D
0.5
4g
Si para 5< Re 800, entonces C D
24
1 0.15 Re 0.687 D
Re 3C D
0.5
4g
Si para 1000 < Re 100000, entonces CD toma un valor constante de 0.4 D
3C D
Ecuación de Rubey
0.5
2 36 2 6
gD 2 para d < 0.1 mm
3 D D
10 0.01gD 3
0.5
D
1 1 para 0.1 < d 1 mm
2
0.806(gD)0.5
En una muestra grande de suelo no cohesivo, es de interés conocer la forma en la que están distribuidos
los tamaños de las partículas más que el tamaño correspondiente a una sola. Esto se obtiene mediante
la curva granulométrica que representa la distribución de las partículas en una muestra de sedimento.
La curva granulométrica se obtiene midiendo la distribución de tamaños de las partículas en una muestra
representativa. Esta medición puede hacerse en arenas y gravas con una separación por tamizado y
para la fracción fina por sedimentación.
Por ejemplo, si d75 = 0.524 mm, significa que el 75 por ciento, en peso, del sedimento está constituido
por partículas cuyos tamaños son menores o iguales a 0.524 mm.
El d75 es utilizado por Lane en su análisis concerniente al inicio del arrastre de sedimentos.
- En el estudio realizado por Meyer – Meter y Muller sobre el arrastre de los sedimentos se considera el
diámetro medio, como el calculado con la relación.
1
dm
100
Pi d i
HIDRAULICA FLUVIAL 17 ING. M. Sc. CESAR MILLA VERGARA
Donde:
Pi : porcentaje del peso total de la muestra correspondiente
di : diámetro de la partícula para el cual el i% en peso de la muestra es menor o igual que ese
tamaño
El análisis de la distribución granulométrica de los sedimentos naturales muestra que están constituidos
por partículas cuyos tamaños no se distribuyen siguiendo una función de probabilidad determinada, ya
que su distribución se encuentra influenciada por diversos factores, donde destaca el mecanismo de
transporte y la pendiente sobre la cual los sedimentos son arrastrados.
Distribución Normal
Las arenas finas y los limos que se encuentran en el cauce de los ríos de planicie frecuentemente
presentan una distribución granulométrica que se aproxima a una distribución normal.
Si resulta que los puntos graficados sobre un papel de distribución normal se alinean sobre una recta, los
diámetros de las partículas del suelo se ajustan a una función de probabilidad normal o Gaussiana, y la
distribución granulométrica puede describirse con la siguiente ecuación:
d i d 50 Z n
Donde:
di : diámetro cuyo porcentaje en peso es menor o igual
d50 : diámetro que corresponde a la mediana de la granulometría para la cual 50% de las partículas
tiene un tamaño menor y el otro 50% es mayor
Zn : es la variable aleatoria estandarizada, que tiene una distribución normal, con media igual a cero
y desviación estándar igual a uno.
: desviación estándar que representa una medida de dispersión que indica qué tan alejados
están los datos respecto a un valor central y se define como:
d 84 d 50
Si se dibuja el logaritmo del tamaño (log d) contra el porcentaje acumulado que pasa y se obtiene una
recta como curva de ajuste, se dice que la muestra se distribuye de acuerdo con una distribución log –
normal.
La distribución granulométrica del tipo log – normal puede describirse mediante la expresión siguiente:
Donde:
g : desviación estándar geométrica definida por:
d 84
g
d 50
Si la granulometría tiene una distribución log – normal, cualquiera de las igualdades siguientes se
cumple:
d 84 d 50 d 84
g
d 50 d16 d16
d 50 d 84 / d16
La distribución log – normal no es simétrica, ya que la mediana d50 y la media dm no son iguales.
1
d m d 50 exp (log g ) 2
2
Distribución circular
Cuando los sedimentos transportados por el agua se depositan en lugares donde escurren corrientes
sobre pendientes pronunciadas, se observa gran abundancia de cantos rodados y boleos. Si la
granulometría del fondo sigue una distribución circular, los tamaños de las partículas se distribuyen de
acuerdo con la expresión siguiente:
p
2
d i d máx 1 1
100
Donde:
di : diámetro cuyo porcentaje en peso es menor o igual
dmáx : diámetro máximo de partícula
p : por ciento que pasa
Para analizar si una muestra de sedimento sigue una ley circular es necesario graficar la curva
granulométrica en un espacio donde el eje de las abscisas corresponde al diámetro de las partículas en
escala aritmética y el eje de las ordenadas corresponde al por ciento que pasa utilizando la escala que
resulta de emplear la expresión de la ley circular. Si la distribución granulométrica sigue la ley circular, se
observa que los puntos de la curva quedan alineados sobre una recta.
Si se adoptan en papel aritmético escalas tales que las distancias representativas del diámetro máximo y
de 100% sean iguales, se obtiene un cuarto de circunferencia con radio igual a D máx.
Distribución logarítmica
Si al dibujar los puntos de la curva granulométrica en papel semilogarítmico, éstos quedan alineados
sobre una recta, la distribución de los tamaños de las partículas es logarítmica y puede describirse
mediante la ecuación:
d n d 50 gpn
d n d 50 gqn
Donde:
qn : variable que depende del porcentaje n correspondiente al diámetro que interesa conocer. Se
calcula como:
n
q n 4.43835Log
50
Es este caso se verifica que:
0.45531 0.31286
d d d
g 84 50 84
d 50 d16 d16
d 50 d16 d 84 0.68714
0.31286
f) Densidad (ρ)
La densidad se define como la masa por unidad de volumen, es decir:
masa
volumen
3
Las unidades en el sistema internacional son kg/m . La densidad varía con la temperatura. A 4° C el
3
agua tiene una densidad de 1000 kg/m .
W
V
3 2 2
Las unidades en el sistema internacional son N/m , o bien kg/(m s ).
g
9810 N / m3 9.81kN / m3
HIDRAULICA FLUVIAL 20 ING. M. Sc. CESAR MILLA VERGARA
h) Peso volumétrico
El peso volumétrico, V de un conjunto de partículas es el peso de la materia entre el volumen total que
ocupa el conjunto, contenidos dentro de él los huecos o vacíos. Sus unidades y dimensiones son las
mismas que las indicadas para el peso específico.
En el peso volumétrico de una muestra de sedimento habrá que distinguir el peso volumétrico seco que
se tiene cuando todos los vacíos están ocupados por aire, el parcialmente saturado y el sumergido
cuando parte o la totalidad de ellos están ocupados por agua.
El peso volumétrico es útil para estimar la vida útil y los períodos con los que deberán programarse
acciones de desazolve en embalses.
Las relaciones de mayor interés que se utilizan al evaluar el volumen real ocupado por los sedimentos al
depositarse, se muestra en la figura:
Volúmenes Pesos
Vv
VT VW Fase líquida
WW WT
VS Fase sólida
Ws
WS
V
VT
WS : peso de sólidos
VT : volumen total
Si se define la porosidad n como la relación entre el volumen de vacíos al volumen total como:
VV
n
VT
Donde:
VV : volumen de vacíos
VT : volumen total
Se cumple:
V S (1 n)
Por otra parte, la relación de vacíos e se define como el volumen de vacíos entre el volumen de sólidos,
por tanto, se pueden escribir las ecuaciones siguientes:
HIDRAULICA FLUVIAL 21 ING. M. Sc. CESAR MILLA VERGARA
VV n
e
VS 1 n
Donde:
VV : volumen de vacíos
VS : volumen de sólidos
e
n
1 e
S
V
1 e
i) Peso específico sumergido (γS)
Cuando las partículas están dentro del agua, el empuje hidrostático influye en el peso específico.
S' S
En hidráulica fluvial es común utilizar el peso específico relativo de los sólidos sumergidos ( ), el cual
se expresa como:
s' s 1
S S
s' s 1
j) Viscosidad
Es la propiedad de los fluidos, de resistencia a un movimiento interno o a su deformación angular.
Newton obtuvo que en un fluido en movimiento, la fuerza interna de frotamiento por unidad de área o
esfuerzo tangencial t es proporcional al gradiente transversal de velocidades du / dy .
du
dy
Donde:
: es el coeficiente de proporcionalidad, que se denomina viscosidad dinámica o absoluta, y es
característico de un fluido.
Si la relación entre las velocidades de deformación y los esfuerzos que las producen es lineal, el fluido se
denomina Newtoniano. El aire y el agua son un ejemplo de los fluidos Newtonianos.
La viscosidad varía con la temperatura del fluido, en García Flores y Maza (1998) pueden obtenerse
valores de la viscosidad y de otras propiedades del agua para diferentes temperaturas. Para una
3 2
temperatura de 20° C, la viscosidad dinámica del agua es de 1x10 N.s/m = 1 centi-poise = 0.001 Pa.s =
2
0.01 Poise = 0.001 N.s/m .
m w Cv( S w )
W
m
1
Cs1
S W
S
m W 1
Cs2
S W
S
m W
Cs3
S W
S
Esta propiedad tiene importancia en el diseño de canales sin revestimiento; la pendiente de los taludes
de los canales se relaciona con el ángulo de reposo, que se incrementa al disminuir la redondez y al
aumentar la rugosidad de la textura.
Para estimar el ángulo de reposo, θ, por debajo del nivel del agua para un rango de
0.16 d 0.59 mm y 0.90 s 1.63 , Gibson propuso la expresión siguiente, Garde y Ranga Raju
(1985):
Para el diseño de la estabilidad de canales, Van Rajin (1993) propuso, de acuerdo con la agudeza de las
aristas de las partículas, los valores del ángulo de reposo que se indican en la tabla siguiente.
La variación del ángulo de reposo fue estudiada por Simona y Senturk (1976), para materiales uniformes
de lignita, bakelita, piedra pómez, arenas y gravas.
Para cada tipo de material, el ángulo de reposo decrece con el incremento del tamaño de la partícula
hasta alcanzar un mínimo luego del cual el ángulo de reposo aumenta nuevamente. Para arenas y
gravas el valor mínimo de θ ocurre con un valor de d = 2.4 mm, Ver figura siguiente.
Los datos del ángulo de reposo de un enrrocado (riprap) están contenidos en dos conjuntos de curvas
publicadas por Lane (1955) y Simona (1957), esta última se muestra en la figura siguiente.
El ángulo de reposo está dado como una función directa del tamaño medio y la angulosidad de la piedra.
La curva de Lane y Simona, difieren más por tamaños de piedra pequeñas que para las grandes.
piedra quebrada
grados
40
muy angular
35 muy redondeado
ANGULO DE REPOSO
30
25
20
2 4 7 2 4 7 2 4 7 2 4
10-3 10-2 10-1 100
TAMAÑO MEDIO DE PIEDRA (d50),pies
Un importante aspecto de los procesos fluviales en ríos, es que la hidráulica del flujo, el transporte de
sedimentos, la formación del cauce y los cambios migratorios están interrelacionados. Como el flujo en
ríos es un flujo en canal abierto, los principios y ecuaciones que se aplican a canales no prismáticos con
flujo en canales abiertos, se aplican también a la hidráulica de ríos.
Debido a que el flujo en canal abierto es un flujo con superficie libre, una importante tarea en la hidráulica
de ríos, es determinar la superficie libre, y de esta manera determinar los parámetros hidráulicos, tales
como radio hidráulico, velocidad, gradiente de energía.
a) Esfuerzo Cortante
Cuando el agua fluye en un canal, se desarrolla una fuerza que actúa sobre el fondo en dirección del
flujo, esta fuerza que es el empuje del agua sobre el área mojada, se conoce como fuerza tractiva o
esfuerzo cortante. El esfuerzo cortante es la fuerza por unidad de área en la dirección del flujo. Si
consideramos flujo uniforme permanente bidimensional en un canal como se muestra en la Fig. 1. el
esfuerzo interno 't en un nivel z por encima del fondo del canal puede ser obtenido considerando un
volumen de control con su superficie definida por ABCD. y un ancho unitario perpendicular a la
superficie. El eje coordenado x es en la dirección del flujo a lo largo de la pendiente del canal S. y el eje
coordenado z es perpendicular al flujo.
Las fuerzas que actúan sobre el elemento con una longitud dx en la dirección x son las siguientes:
Las fuerzas de presión hidrostática sobre AB y CD, respectivamente.
El esfuerzo cortante T actuando sobre BC
El peso W x del fluido en el volumen de control en la componente x
Bajo las condiciones de flujo uniforme, las fuerzas hidrostáticas se encuentran en equilibrio y por lo tanto
Wx debe ser balanceada por el esfuerzo cortante, esto es:
WX dx 0 (1)
WS dx 0 (2)
Reemplazando el peso W por γ(D -z)dx, donde γ es el peso específico del fluido y D es el tirante de
flujo, en esta ecuación, se tiene:
D z S (3)
Si el canal es tridimensional, el esfuerzo cortante promedio τ0 sobre la superficie mojada del canal, se
puede obtener siguiendo el mismo procedimiento.
Consideremos el volumen de control AEFD, con una longitud de tramo dx y una sección transversal del
canal con área A. El peso del volumen de fluido es γAdx. La componente x del peso W se equilibra con
el esfuerzo cortante sobre la superficie mojada, esto es:
0 RS (5)
i) Fórmula de Manning
La fórmula de Manning, es quizás la más ampliamente utilizada. En el sistema internacional tiene la
siguiente expresión:
1 2 / 3 1/ 2
U R S (6)
n
Donde
U= Q/A. es la velocidad media de la sección en m/s
n = coeficiente de rugosidad
Valores de n para diferentes condiciones de canal, se encuentran en tablas.
Aquellos valores de n para ríos aluviales deben ser considerados como promedios puesto que el flujo
induce tipos de rugosidades de acuerdo a la condición del flujo.
La vegetación genera mayor resistencia al flujo a velocidades bajas, mientras que causa resistencia muy
baja a flujos grandes.
U C RS (7)
Para evaluar el coeficiente de C, se puede utilizar la siguiente expresión obtenida a partir de la teoría de
Prandt y de las experiencias de Nikuradse para secciones rectangulares muy amplias:
11.1R
C 18Log (8)
kS
Donde kS es el diámetro de las partículas si el fondo es plano en m. En cauces naturales se recomienda
tomar ks= 2dm (donde dm es el diámetro medio).
Donde τ0= es el esfuerzo cortante medio en las fronteras y ρ es la densidad de masa del fluido. Si el
esfuerzo cortante es reemplazado por γRS, la ecuación resulta:
1/ 2
8 gRS U 8
U ó
U * f
(10)
f
uz
R (12)
Donde:
u = velocidad a una distancia z desde la frontera sólida
υ = es la velocidad cinemática
Dentro de la sub capa laminar la velocidad puede ser asumirse que varía linealmente con respecto a z, y
así el esfuerzo cortante τ0 está dado por la ley de viscosidad de Newton.
du 0z
0 ó u (13)
dz
Donde:
μ = viscosidad dinámica
En el final de la sub capa laminar, z = δ (espesor de la sub capa laminar) y la ecuación (13) da la
velocidad como:
0
u (14)
El número de Reynolds en el cual el flujo laminar comienza su transición a flujo turbulento ha sido
uz
determinado experimentalmente como un valor cercano a 11.6, esto es, = 11.6. Luego
el espesor de la subcapa laminar es:
11.6 (15)
U*
Aunque la sub capa laminar existe sobre superficies lisas, este puede ser disturbado en superficies
rugosas, si los elementos rugosos sobresalen hacia la zona turbulenta. Sobre la base del tamaño de la
k SU *
5 (16)
Donde ks= diámetro medio de granos de arena utilizado por Nikuradse y es conocido como rugosidad de
arenas de Nikuradse.
Una lista de valores equivalentes de ks para diferentes tipos de superficies se muestra en la tabla 3.1
Tabla 3.1: Valores de ks en pies
Material ks (pies)
Cemento 0.0013 – 0.0040
Concreto 0.0015 – 0.0100
Lecho de río natural 0.1000 – 3.0000
Para régimen de transición, el número de; Reynolds de corte está en el rango de 5 - 70. Los elementos
rugosos penetran parcialmente la zona de la sub capa laminar para afectar el flujo turbulento.
Para régimen hidráulicamente rugoso, el número de Reynolds de corte es mayor de 70. Los elementos
rugosos penetran a través de la sub capa laminar hacia el fondo de la zona turbulenta.
l t (17)
Donde Tl es el esfuerzo de corte causado por la viscosidad dinámi.ca μ y τ t es el esfuerzo de corte
causado por los remolinos turbulentos. La cantidad τ l está dada por la ley de Viscosidad de Newton, en
donde el esfuerzo de corte es directamente proporcional al gradiente de velocidad duJdz.
El esfuerzo cortante turbulento está dado por:
Donde u' w' es el tiempo promedio del producto de las velocidades fluctuantes u' y w' en las direcciones
x y z, respectivamente. Definiciones de tiempo-velocidad promedio u y su velocidad fluctuante u', están
dados en la Figura 3. La velocidad instantánea es la suma de estas dos cantidades.
Velocity at a poitn
U
Time
Fig. 3: Diagrama esquemático de la velocidad en flujo turbulento
Prandt sugirió una relación para el esfuerzo turbulento en términos del tiempo-velocidad promedio u,
como:
2
du
l l
2
(19)
dz
Donde 1 es la longitud de mezcla propuesta por Prandtl. Como se muestra en la Figura 2, una partícula
es llevada a través de una cierta distancia vertical l como resultado de un positivo w'. Concomitante a
este movimiento, la velocidad de la partícula cambia por una cantidad del orden de -1 du/dz. El signo
negativo indica a una magnitud menor que el tiempo promedio. Se asume que u' se debe a esta
diferencia de velocidad y que w' tiene el mismo orden de magnitud como u', entonces:
du du
u ' l ; w' l (20)
dz dz
du du
t l 2 (21)
dz dz
l z (22)
Donde K es una constante conocida como constante de von Karman. Sustituyendo l=Kz en la ecuación
(19), resulta:
2
du
t K z
2 2
(23)
dz
Esta ecuación para flujo turbulento es válida solo en el núcleo turbulento afuera de la sub capa laminar.
Esta ecuación puede ser escrita como:
du
Kz U* (24)
dz
Aquí existe una sub capa laminar sobre una superficie lisa. De la ecuación (13), para z=δ, la velocidad u
es τ0 δ/μ. Sustituyendo estos valores en la ecuación (25), se obtiene el valor de C, obteniendo así el perfil
de velocidad para la región turbulenta:
100Z
V ( Z ) 5.75V* Log (26)
33Z
V ( Z ) 5.75V* Log (27)
S
K
33Z
V ( Z ) 5.75V* Log (28)
K S 0.3
e) Resistencia al Flujo en Cauces de Fondo Fijo
La resistencia al flujo puede asumirse que consiste de dos componentes: la rugosidad de la partícula y la
rugosidad de la forma.
Cualquier de estas rugosidades puede ser asumida en términos del coeficiente n de Manning, C de
Chezy o f de Darcy Weisbach.
En canales cubiertos con arenas o gravas, la resistencia al flujo, en ausencia de la forma del lecho,
puede ser considerada que es causada principalmente por la rugosidad de las partículas.
La fórmula de Strickler (1923) que define la rugosidad de Manning como una función del tamaño de la
partícula, es:
d 1/ 6
n (29)
21.1
Donde d representa el diámetro de arena uniforme en metros. Sustituyendo la fórmula de Strickler para n
en la ecuación de Manning, se tiene la llamada fórmula de Manning-Strickler.
1/ 6
U R
6.74 (30)
U* d
HIDRAULICA FLUVIAL 31 ING. M. Sc. CESAR MILLA VERGARA
Donde U * gRS , es la velocidad de fricción. Las cantidades R y d tienen las mismas unidades.
Muchas otras fórmulas caen dentro de la misma forma que la fórmula de Strickler. Así por ejemplo,
Meyer-Peter y Muller (1948) desarrollaron la siguiente relación para arenas mixtas.
n
d 90 1/ 6 (31)
26
Donde d90 está en metros. En el caso de canales revestidos con grava) Lane y Carlson (1953) sugirieron
la siguiente fórmula basada en un estudio que realizaron en el Canal San Luis en Colorado:
n
d 75 1/ 6 (32)
39
Donde d75 está en pulgadas.
U R
6 2.5 ln (33)
U* kS
1/ 6
ng 1 / 2 d
0.225 50 n 0.0395d 50
1/ 6
ó (34)
D
1/ 6
D
Donde d84 es usado como la altura de rugosidad representativa del material gravoso no uniforme. El
coeficiente a es usado para definir el efecto de la geometría de la sección transversal sobre la resistencia
al flujo. Su valor se encuentra en el rango de 11.1 a 13.46.
La ecuación de Rey fue desarrollada con datos de descargas desde menores de 1 hasta 444 pies
cúbicos por segundo.
Usando datos de 67 ríos en Alberta, Canada, Bray (1979) obtuvo coeficientes de mejor ajuste para una
ecuación logarítmica de resistencia:
1 D
1/ 2
0.248 2.36 log (36)
f d 50
Donde:
H = es la carga de energía total
Z = elevación de la superficie de agua
α = coeficiente de energía.
La altura de velocidad se multiplica por α para considerar el flujo no uniforme en la sección transversal,
cuyo valor está dado por:
u dA u
3 3
Ai
i
(39)
U 3A U A 3
Donde ui es la velocidad media del subarea Ai. El valor de α se incrementa cuando es mayor la no
uniformidad de flujo. Para canales naturales sin cubierta de hielo, el valor de α varía en el rango de 1.15
a 1.50, con un promedio de 1.30.
Diferenciando la ecuación (38) con respecto a la distancia a lo largo del canal, se tiene:
dH d U2
Z S
2 g
(40)
dx dx
dH
Donde S es el gradiente de energía que es igual a por definición. Escribiendo en forma de
dx
diferencias finitas, la ecuación anterior se transforma en:
U 2j U 2j 1
Zj j Z j 1 j 1 Sx (41)
2 g 2 g
Donde j es el sub índice de la sección transversal, contabilizado de aguas arriba hacia aguas abajo, y ∆x
es la distancia entre las secciones j y j+1. El gradiente de energía S es para el tramo ∆x, y usualmente
es tomado como la media geométrica de dos secciones adyacentes:
S S j S j 1
1/ 2
(42)
Donde:
H’l = es la pérdida de carga debido a la resistencia del contorno
He = es la pérdida de carga debido a las contracciones y expansiones de flujo
”
H l = es la pérdida de carga debido a las corrientes secundarias
Si se quiere tomar ambos, flujos sub crítico y supercrítico. Es necesario aplicar las relaciones para
tirantes alternos y conjugados.
El método del paso estándar utiliza la ecuación de energía para el cálculo del perfil hidráulico, tanto para
canales artificiales como naturales. La geometría del canal está definida en la sección transversal
discretizada. Los parámetros hidráulicos pertinentes como área de flujo, perímetro mojado y rugosidad
deben ser expresados en términos del estado para cada sección transversal. La elevación de la
superficie de agua es calculada de sección a sección siguiendo la dirección apropiada. Avanzando de la
sección j+ 1 a la sección j, por ejemplo, el estado de la sección j denominada Zj es calculada del
conocimiento de la etapa Zj+l. Puesto que Zj está implícito en la ecuación (58) este se determina
siguiendo una aproximación ensayo-error, como sigue:
a) Generalidades
El inicio de movimiento de las partículas que componen el lecho ocurre cuando los esfuerzos
hidrodinámicos actuantes superan los esfuerzos de resistencia. Este inicio de movimiento no es
instantáneo para todas las partículas de un determinado tamaño que cubren el lecho. Solo una parte de
estas partículas entran en movimiento, mientras que otra parte permanece en reposo. Esto se debe a la
naturaleza turbulenta del flujo, que determina la fuerza tractiva sobre la partícula. La condición crítica de
inicio de transporte es definida como el estado en que una parte representativa del material del lecho
empieza a moverse. Esta condición es determinada a través de observaciones y tiene un carácter
subjetivo.
Los esfuerzos de resistencia, al movimiento de las partículas, dependen del tamaño y la composición
granulométrica de los sedimentos. Los sedimentos muy finos, que contienen una cantidad apreciable de
limos y arcillas, resisten al movimiento a través de esfuerzos de cohesión. La complejidad del fenómeno
de cohesión entre las partículas, quizás sea el factor principal para explicar la existencia de poquísimos
trabajos relativos a este tipo de material. Los sedimentos no cohesivos, constituidos por las arenas,
gravas y piedras resisten al movimiento, principalmente debido al peso de los granos.
Existen básicamente dos tipos de enfoques del problema, que comprenden casi la totalidad de los
trabajos existentes. El primero, es el criterio de utilización de la velocidad crítica y el segundo el criterio
de utilización de la fuerza tractiva crítica.
La Figura 1 (a) y (b) muestran las fuerzas actuante sobre una partícula de sedimento sumergida para el
caso de flujo laminar y turbulento respectivamente.
Donde:
FD 0 C2 D 2 es la fuerza de arrastre
U2
FL CLC3 D 2 la fuerza de sustentación
2
En la Figura 1 y las ecuaciones anteriores: θ es la pendiente del canal, Φ el ángulo de reposo del
material, Cl el coeficiente de uniformidad del material e igual a π/6 para partículas esféricas; C 2 y C3
coeficientes de forma y D el diámetro característico de la partícula.
i) Análisis para Flujo Laminar: Cuando el flujo es laminar, el efecto viscoso predomina y el flujo bordea
la partícula y la resultante de la fuerza de arrastre actúa por encima del punto C.
White estudió el equilibrio de una partícula en flujo laminar y determinó la siguiente ecuación para definir
el esfuerzo de corte crítico:
Desde el punto de vista de ingeniería θ es muy pequeño por lo que: tan 0 y cos 1:
C1a1
C D S tan (4)
C2 a2
C1a1
k 0.18
C2 a2
ii) Análisis para Flujo Turbulento: De la Figura 1b se puede ver que a1=a2, por lo que
C1
k y la ecuación (4) se transforma en:
C2
C
k . tan (5)
D S
El primer miembro de la ecuación (5) es un adimensional que relaciona la fuerza de arrastre y la fuerza
gravitacional y es muy usado en la definición de inicio del movimiento de la partícula.
c) Análisis de Shields
Shields, fue el primero en estudiar el inicio del movimiento de la partícula considerando las fuerzas que
actúan sobre ella y aplicando los principios de similitud. Este análisis se presenta en detalle a
continuación.
La fuerza, FS requerida para mover una partícula de tamaño D y de específico γ S desplazado dentro de
un fluido de peso específico γ es: F C1 S D 3 ; donde C1 es un coeficiente que depende solo de
las características del sedimento como forma de la partícula.
HIDRAULICA FLUVIAL 36 ING. M. Sc. CESAR MILLA VERGARA
De manera similar, la fuerza hidrodinámica ejercida por el fluido sobre la partícula es:
u a2
Fl C D C 2 D 2 , donde ua es la velocidad característica, CD es el coeficiente de arrastre de la
2
3
partícula al número de Reynolds correspondiente a ua y C2 es un coeficiente tal que C2D da el área
proyectada de la partícula. La velocidad característica ua, puede ser tomado como la velocidad en el tope
de la partícula, ud. Empleando la ecuación de Karman – Prandtl para la distribución de velocidades, ud
ud u D 0
puede ser expresado como: f1 * , donde u* es la velocidad de corte dada por . Como
u*
u D u D
CD f 2 d puede asumirse que: C D f 3 * .
u* D
Denotando Re * , F1 puede ser expresado como: F1 f 3 (Re* ) u*2 f12 (Re* )C 2 D 2
2
Igualando F y F1 para una condición de movimiento incipiente e introduciendo la letra “c” para denotar la
condición crítica, se obtiene:
u*2c 2C
C1 S D 3 f 3 Re *c u*2 c f12 Re *c C2 D 3 de donde: 1 f Re *C siendo:
2 S D C2
u D
Re *C *C y u*C la velocidad de corte para un movimiento incipiente. Así para partículas de una
forma dada la condición crítica de movimiento incipiente es:
C
f Re *C (6)
S D
C
La figura 2 muestra la variación de con Re *c obtenido por Shields, basado en datos
S D
experimentales. La correlación obtenida en la Figura es significativa y se puede afirmar que Re *C es
D
proporcional a la relación entre el tamaño de la partícula y espesor de la sub capa límite laminar: . De
esta forma el gráfico representado por la Figura es similar al gráfico del factor de fricción para tuberías en
la zona de transición. La porción de línea recta del lado izquierdo del gráfico representa el caso donde
' D y la partícula esta sumergida completamente en la sub capa laminar. Cuando el espesor de la
sub capa laminar es del mismo orden de magnitud que D, la turbulencia disturba la sub capa laminar y
afecta el flujo alrededor de la partícula. La depresión en la curva de la Figura para 0.25 Re *C 40
C
representa este caso. Cuando D , la sub capa laminar es destruida y para el
S D
D
movimiento incipiente llega a ser independiente de . Así para materiales muy gruesos en el
'
C
movimiento incipiente se da 0.06
S D
Este criterio considera que el movimiento ocurre debido a la acción del impacto del flujo sobre la
partícula. La velocidad de referencia, que puede ser una velocidad en las proximidades del lecho, o
velocidad media, es relacionada con el diámetro de la partícula. Las primeras observaciones de la
condición de velocidad crítica para el inicio de movimiento fueron hechas por DuBuat, y otro ejemplo
clásico es el trabajo de Fortier y Scobey, que fijaron las velocidades permisibles en canales que fueron
usados en diseño de canales por muchos años. El trabajo de Hjulstrom, luego de un análisis extenso de
datos obtenidos por diversos autores, dio como resultado una relación entre la velocidad media del flujo
en el inicio del movimiento y el tamaño de los sedimentos. Las curvas fueron determinadas para flujos
con profundidad mínima de 1,0 m.
La crítica a este método es que la velocidad no es suficiente para proveer informaciones sobre el inicio
de movimiento de las partículas. Se sabe que dos flujos con la misma fuerza tractiva en el fondo,
granulometrías idénticas y las mismas distribuciones de velocidades, pueden tener velocidades medias
diferentes si las profundidades fuesen diferentes. Por esta razón, es recomendable que se emplee el
criterio del esfuerzo crítico de corte siempre que sea posible.
Agua
Coeficiente de Agua clara Agua transportando transportando
Tipo de Material
Rugosidad limo coloidal limo, arena y grava
(n) (m/s) (m/s) (m/s)
Arena fina coloidal 0.020 0.46 0.76 0.46
Greda arenosa no coloidal 0.020 0.53 0.76 0.61
Greda limosa no coloidal 0.020 0.61 0.91 0.61
Limo aluvial 0.020 0.61 1.07 0.61
Greda firme 0.020 0.76 1.07 0.69
Ceniza volcánica 0.020 0.76 1.07 0.61
Grava fina 0.020 0.76 1.52 1.14
Arcilla dura coloidal 0.025 1.14 1.52 0.91
Greda graduada a guijarro 0.030 1.14 1.52 1.52
Limo aluvial coloidal 0.025 1.15 1.52 0.91
Limo graduado a guijarro 0.030 1.22 1.68 1.52
Grava gruesa 0.025 1.22 1.83 1.98
Guijarro y ripio 0.035 1.52 1.68 1.98
Capas dura 0.025 1.83 1.83 1.52
i) Ecuación de Kramer:
Obtenida en canales experimentales de 14 m de longitud, 0.81 m de ancho y 0.30 m de profundidad
usando partículas de cuarzo de densidad relativa de 2.70:
10 4
C S D (7)
6 M
Donde:
C Fuerza tractiva crítica en N/m 2
S Peso específico del sedimento en N/m3
Peso específico del agua en N/ m3
HIDRAULICA FLUVIAL 38 ING. M. Sc. CESAR MILLA VERGARA
D = Diámetro medio en mm (Debe de variar entre 0.24 mm y 6.52 mm.)
M = Coeficiente de uniformidad de Kramer. (Debe variar entre 0.265 y 1.0)
1/ 2
D
C 0.285 S (8)
M
Donde:
C Fuerza tractiva crítica en N/m 2
S Densidad del sedimento en kg/m 3
Densidad del agua en kg/m3
D = Diámetro medio en mm (Debe de variar entre 0.205 mm y 4.077 mm.)
M = Coeficiente de uniformidad de Kramer. (Debe variar entre 0.28 y 0.643)
S
C 0.754 D (10)
Donde;
C Fuerza tractiva crítica en N/m 2
S Densidad del sedimento en kg/m 3
Densidad del agua en kg/m3
D = Diámetro medio en mm
v) Fórumula de Indri:
S 1
Cuando D 1.0 mm; entonces C 0.130D 0.12
m
(11)
S 1
Cuando D 1.0 mm; entonces C 0.538D 0.73
m
Donde;
a) Generalidades
El régimen de flujo en lecho móvil, es definido, por Garde y Albertson, de la siguiente manera: La
naturaleza de la configuración del lecho y de la superficie líquida varían de acuerdo con las
características del sedimento, del flujo y/o características del fluido. Estos tipos de configuraciones tanto
del lecho como de la superficie líquida son clasificados de acuerdo con sus características y
denominados "Regímenes de Flujo".
Se debe tomar cuidado de no confundir esta definición con otras nomenclaturas semejantes de la
hidráulica de canales.
Los diferentes regímenes de flujo fueron observados en canales naturales y descritos por Albertson,
Simons y Richardson; a partir del reposo y por sucesión de ocurrencias conforme la velocidad del flujo
aumente.
i) Fondo Plano: Hasta el momento en que los sedimentos no alcanzan las condiciones limites para el
inicio del movimiento, el lecho se mantiene en reposo.
ii) Rizos: Cuando el sedimento inicia el movimiento, ocurren pequeñas deformaciones cuyo corte
longitudinal se asemeja a los dientes de una sierra. En general el talud de aguas arriba es bastante
suave y el de aguas abajo mas inclinado, alcanzando el ángulo de reposo natural del sedimento. Si el
material de fondo fuera fino, los rizos se forman rápidamente luego del inicio del movimiento. Los
materiales groseros con diámetros del orden de 1.0 mm o mayor, no producen este tipo de formación y el
lecho permanece plano por mas tiempo hasta la aparición de las dunas.
iii) Dunas: Cuando la velocidad aumenta, aparecen conformaciones periódicas mayores, con una forma
semejante al de los rizos y con la superficie más irregular. Las dunas pueden alcanzar grandes
proporciones, que a veces reciben el nombre de bancos. Si el material del lecho fuera relativamente fino
puede ocurrir la formación de rizos en el dorso de las dunas, que a su vez pueden ser barridas a medida
que la velocidad aumenta.
iv) Transición: El régimen de transición se caracteriza por una situación bastante inestable, donde
pueden ocurrir cambios rápidos en la forma de la superficie libre y del lecho con solo pequeños cambios
de las condiciones de flujo. Generalmente ocurre cuando el número de Froude es del orden de 0.8. Con
el aumento progresivo de la velocidad, las dunas se van alargando y disminuyendo en amplitud y si el
material fuera relativamente fino, el lecho puede pasar a la forma plana.
vi) Régimen de Rápidos: En este régimen ocurre una sucesión de regímenes rápidos y lentos separados
por resaltos hidráulicos. Ocurren en los estados avanzados del flujo.
Un resumen de un estudio exhaustivo realizado por el U.S. Geologycal Survey, Colorado State
University, sobre deformaciones del lecho, fue preparado por Simons y otros, la misma que se muestra
en la Tabla siguiente.
a) Generalidades:
Cuando el esfuerzo de corte promedio en el fondo excede la fuerza tractiva crítica para el material,
estadísticamente las partículas del fondo empiezan a moverse en la dirección del flujo. Las partículas se
mueven de diferentes formas dependiendo de las condiciones del flujo, tamaño y peso específico de las
partículas. Una forma de movimiento de las partículas es por el rodamiento o deslizamiento a lo largo del
lecho. Tal tipo de movimiento es generalmente discontinuo; la partícula puede deslizarse o rodar por
algún tiempo y quedar estacionario por otro tiempo y nuevamente empezar el movimiento por algún otro
tiempo. El sedimento transportado de esta forma es conocido como transporte por contacto. Una
segunda forma de movimiento del sedimento es conocida como transporte por saltación. Saltación es un
modo importante de transporte en caso de materiales no cohesivos de velocidades de caída
relativamente altas. El tercer modo de transporte es el transporte en suspensión; en este caso las
partículas de sedimento son continuamente soportados por la turbulencia del flujo.
b) Ecuaciones de Cálculo
i) Fórmula de Du Boys
La primera fórmula de transporte sólido de fondo fue propuesta por Du Boys en 1879, que imaginó el
movimiento de los sedimentos de fondo en capas paralelas. La velocidad de cada capa varía desde un
máximo en la superficie del fondo a cero en la capa inferior.
qB K . 0 0 C
0.17 3
K (m /kg/seg)
D3/ 4
El esfuerzo crítico con la expresión:
0
*
S Dm
0 RS
hS
*
Dm
Mínimo Máximo
- Tamaño de las partículas (D y Dm) en m 0.0004 0.030
3
- Pesos específicos, en Kg/m 1250 4200
- Pendientes 0.0004 0.020
- Tirantes, en m 0.01 1.20
2 36 2 36 2
F1
3 gDm3 S gDm3 S
Misri, Garde y Ranga Raju, analizaron una gran cantidad de datos de transporte de sedimentos de
cuarzo uniformes. Basado en el argumento de que el esfuerzo de corte del grano es el responsable del
3
n 2
movimiento de las partículas, ellos establecieron una relación entre: * S *
'
y
n
1/ 2 1/ 2
q 1
B
3
S S
.
gD
v) Ecuación de Einstein
a) CALCULO DE R’ Y R”
1. Asumir un valor de R’, conociendo el valor de R
2. Calcular: V*' g R 'S
ν
3. Calcular: δ' 11.6 '
V*
K d 65
Asumimos:
'
K d 65 g R 'S
Reemplazando se tiene:
δ 11.6 ν
K
4. Con , se obtiene “X” de la lámina 24.
5. calcular:
R 'X
V 5.75V*' Log12.27
0
K d
0 65
X X
R ' X2
V 5.75V*' Log12.27
d 65
6. Calcular:
γs γ w d
' 35
γw R S
V
7. De la lamina 25, se obtiene , que está en función de
V*"
V
8. Calcular V*" , se conoce y V (Paso 5)
V*"
g S
HIDRAULICA FLUVIAL 47 ING. M. Sc. CESAR MILLA VERGARA
10. Verificar si se cumple:
Rc R' R"
R Rc
4. Calcular X :
X
β x Log10.6
0
β Log10.6
β2
5. Calcular * ξY 2
β*
6. Calcular * , conocido * a partir de la lámina 27
γw
TS : transporte de fondo kg/m*s
d : diámetro de las partículas sólidas en m
1.6
TRAMO 1
1
1
X 0.69359
, [2.238 , 100]
K
0.99006 1.1035 * exp 0.55304 *
1.4
TRAMO 2
1.2
X
0.8
0.6
0.1 1 10 100
K/δ
ORIGINAL AJUSTE TOTAL 1° TRAMO 2° TRAMO
ψ en función de (V/V'* )
100
V/V+ '
10
V 0.568
41.959 * exp 0.001533 * 35.407
'
V*
1
0.1 1 10 100
ψ
ORIGINAL AJUSTE
3.6585
K
Y 0.023971* exp 0.20809 * 0.39556, [1.5 , 5]
K
Y
TRAMO 2
TRAMO 1
K K K
Y 1 1.8458 1.3522 * 1.677 0.010291* Ln , [0.3 , 1.5]
0.1
0.1 1 10
K/δ
ORIGINAL AJUSTE 1° TRAMO AJUSTE 2° TRAMO JUSTE
100
75.114
74.035
1 1.4615 * exp 7.7776 *
d
X
10
ξ
0.1
0.1 1 10
d/X
ORIGINAL AJUSTADA
0.1
0.010135
Ф*
0.01
Ln* 1.38482 * 857.148 * exp 10.898 * exp 0.022142 * *
* 14.508 * exp
0.8213
TRAMO 2
0.001 TRAMO 1
0.0001
0.00001
0 5 10 15 20 25 30 35
ψ*
a) Mecanismo de Suspensión
Uno de los problemas de mayor interés en la mecánica de suspensión es el estudio de un método exacto
de cómo las partículas de sedimentos es transportado en suspensión. Es ampliamente conocido que la
turbulencia del flujo es el responsable por la suspensión de las partículas en un curso de agua. Las
partículas en suspensión están sujetas a la acción de la gravedad que causa la sedimentación de las
partículas que tienen mayor peso específico que el agua. Si consideramos, en un primer momento, que
la concentración del material en suspensión es constante a lo largo de la vertical y que el flujo de agua
hacia arriba y hacia abajo son iguales; una mayor cantidad de sedimentos serán transportados hacia
abajo que hacia arriba debido a que la sedimentación ayuda el flujo de sedimentos hacia abajo. Por lo
tanto se ha creado una gradiente de concentraciones debido al transporte de sedimentos y estos
entrarán eventualmente en equilibrio. Bajo condiciones de equilibrio se puede establecer la ecuación de
difusión:
C
C s 0 (1)
y
C
Donde C es el transporte de sedimentos neto hacia abajo y s el transporte de sedimentos neto
y
hacia arriba; C es la concentración del sedimento, y la altura media desde el fondo, s es el coeficiente
de difusión del sedimento y la velocidad de sedimentación.
b) Integración de la Ecuación de Distribución de Sedimentos
La ecuación (1), es la ecuación diferencial de la distribución del material suspendido en la vertical. Esta
C y
ecuación puede ser ordenada de la siguiente manera e integrada entre los limites a e y
C s
perfil de perfil de
velocidad concentración
u C
h
Ca
y
a
fondo
Figura Nº 1: Distribución de la Velocidad y Concentración en la Vertical
Conociendo que el coeficiente de difusión s varía en función de y; la ecuación (1) fue integrada por
Rouse:
Z
C h y a
C a y h a
(3)
El valor de Z es el exponente de la ecuación de distribución de sedimentos y k es la constante de
u* k
Karman y para agua clara es igual a 0.4.
c) Método de Einstein
h
q S C.u.dy (4)
a
a
Z
1 Y Z 30.2hx 1 Y
Z
q S 5.75u*C a h LogY .dY Log
dY (7)
ha A Y D65 A Y
I 1 0.216
1 AZ A Y dY (8 a)
A Z 1 1 Y
Z
I 2 0.216
1 AZ A Y LnY .dY (8 b)
30.2hx
q S 11.6u*C a a 2.3Log I 1 I 2 (9)
D65
Según Einstein a 2D . Las dos integrales I 1 e I 2 fueron resueltos numéricamente por Einstein y
graficados en función de A y Z (figuras 31 y 32). Considerando que la descarga sólida de fondo puede
ser representada por: i B q B 11.6u*Ca a la ecuación de transporte sólido en suspensión puede ser
'
30.2hx
iS q S i B q B 2.3Log I 1 I 2 i B q B P.I 1 I 2 (10)
D65
Garde y Pande obtuvieron a partir de datos de campo una relación entre el caudal sólido en suspensión
y el caudal líquido específico:
4
qS u
0.000051 * (11)
q
e) Método de Samaga
S 30 *6 (12)
Donde:
1/ 2 1/ 2
RS 0 q 1
* S S
3
S S
; (13)
S D S D gD
1.0 1000
s
m/
nc
ne
cio
/s
cm
ta
en
en
10-1 n
cio 100
dim
ta
en
W cm/s
W cm/s
m
di
se
se
de
de
ad
cid
ad
lo
ve
cid
lo
ve
10-2 10
1.0
0.01 0.1 1.0 10 100 1000
Tamaño de grano en mm
LAMINA N° 30
FUNCION I1 EN TERMINOS DE
z=
0
1000
0
.2
100
0
.4
0.6
10
0.8
1.0
1.0
1.2
1.5
2.0
2.5
3.0
0.1
3.5
4.0
5.0
0.1
0.01
0.00001 0.0001 0.001 0.01 LAMINA Nº 31
Z=
0
1000
0.
2
0.
4
100
0.
6
I2 0.8
10
1.0
1.2
1.5
2.0
1.0 3.0
4.0
5.0
0.1
0.01 0.1
0.00001 0.0001 0.001 0.01
LAMINA Nº 32
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.02
0.01
0.2
1/6
n/h = 0.03
0.1
0.09
0.08
0.07
0.06
0.05
0.01
0.03
0.02
0.01
1 2 1 4 5 6 7 8 9 10 2 3 4 5 6 7 8 9 10
w/(t0/p) (1/2)
LAMINA N° 33
4.1 SOCAVACION
a) Generalidades
En los últimos años se han presentado grandes desastres en nuestro país debido a problemas de
erosión y flujo de Iodos y material grueso arrastrado por avenidas extraordinarias. En muchos de los
casos el problema es atribuible a una mala planeación, diseño o construcción de obras, y asentamientos
cercanos y dentro del área de influencia del cauce de los ríos. Para un buen diseño de dichas obras es
necesario evaluar la posible socavación que pueda afectar sus estructuras. Es indispensable que en la
solución de este problema se considere la experiencia y el juicio de los ingenieros especialistas en
hidráulica, así como el análisis extenso de los casos de falla que se han presentado, ya que de ellos
pueden obtenerse lecciones provechosas para conocer mejor este complejo problema.
Uno de los objetivos principales del diseño hidráulico de una obra de cruce debe ser la reducción de las
turbulencias provocadas por la obstrucción que representan las estructuras al libre curso de la corriente.
Se recomienda por ejemplo, evitar en lo posible cualquier obstrucción del cauce si esto no es posible se
deben diseñar transiciones graduales cuidando mantener la capacidad del cauce, que en el caso de un
puente puede resultar en claros considerablemente grandes.
Generalmente el cauce de los ríos se forma con el efecto de avenidas ordinarias de duración
considerable, no precisamente las máximas, de aquí que el criterio de selección del caudal de diseño
sea fundamental para el proyecto adecuado de la obra.
Cuando por una obra o por un proceso natural el ancho del cauce se reduce, el fondo tenderá a
erosionarse. Si se incrementa la pendiente al rectificar un cauce, se puede generar una socavación local.
Si los cambios de dirección se fijan en algún punto, es posible que se aceleren en otro. En menor o
mayor grado, todas las secciones y tramos de un río se encuentran sujetas a un proceso de erosión y
sedimentación.
Las necesidades de la actividad humana modifican los procesos naturales, en algunas ocasiones
acelerándolos o retardándolos mediante obras y construcciones que afectan el cauce, las curvas de los
ríos o los almacenamientos. Es importante tener presente que la existencia de una construcción en el río
tiene un efecto sobre el flujo en consecuencia, la alteración de su equilibrio natural queda sujeta a este
efecto. Estas alteraciones pueden ser locales, como el caso de una pila de puente ubicada al centro del
cauce o puede incluir varios kilómetros, como es el caso de las presas donde, se afecta tanto el tramo
aguas arriba, el almacenamiento y el tramo aguas abajo de la cortina.
b) Tipos de Socavación
Los tipos de socavación más importantes que se consideran en los estudios son los siguientes:
Socavación general
Socavación transversal
Socavación en curvas
Socavación local
Socavación aguas abajo de grandes presas
Socavación al pie de obras de excedencias y deflectores de descarga
Socavación bajo tuberías
i) Socavación general
Consiste en una disminución generalizada del nivel del fondo como consecuencia del incremento de la
capacidad del flujo para transportar material durante las avenidas. También se debe al desequilibrio
ocasionado por la diferencia de capacidad de transporte entre dos secciones consecutivas que definen a
HIDRAULICA FLUVIAL 62 ING. M. Sc. CESAR MILLA VERGARA
un tramo; aquí las erosiones tienen lugar cuando el volumen de sedimento que sale por la sección de
aguas abajo en un cierto periodo es mayor que el que entra al tramo por la sección aguas arriba en el
mismo periodo. Este fenómeno es un proceso natural que puede ocurrir a todo lo largo del río y no es
provocado por factores humanos.
Generalmente la profundidad máxima se observa en el tramo de la mitad de la curva hacia aguas abajo;
sin embargo y por su posición, depende de la dirección del flujo a la entrada y su trayectoria a lo largo de
la curva. Para problemas prácticos se considera que puede existir en cualquier sección de la curva.
Se considera una distribución homogénea del material en el subsuelo cuando existe sólo un tipo de
material en toda la profundidad de socavación, y una distribución heterogénea cuando dos o más capas
de diferente tipo de material pueden ser erosionadas.
Ue Ur (1)
Donde:
Ue : velocidad media de la corriente necesaria para erosionar al material del fondo (inicio de
arrastre), en m/s.
i) Cálculo de Ur
La hipótesis fundamental consiste en suponer que el caudal que pasa por cualquier franja de ancho fijo
de la sección permanece constante durante el proceso de erosión.
La variación de la velocidad media del flujo Ur, como función de la profundidad en un punto de la sección
puede ser obtenida mediante el análisis de una franja vertical de ancho B (ver figura 1).
La ecuación de Manning aplicada a una franja vertical, de ancho B , para el caudal correspondiente
Q es:
1 1/ 2 5 / 3
Q UA S d 0 B (2)
n
Donde:
S : gradiente hidráulico
do : profundidad inicial de la franja analizada, medida sobre su eje vertical entre el nivel del agua
correspondiente al caudal de diseño y el nivel del fondo correspondiente al estiaje, en m
B : ancho total de la sección. Para encontrar B se traza una línea perpendicular a las líneas de
corriente proyectada sobre la superficie libre del agua. Así se toma en cuenta además el
esviajamiento de la corriente.
a=0
b1 b2
Planta
B
d0
ds
Perfil antes de
la socavacion
Perfil de equilibrio al P
final de la
socavacuion
Perfil
Figura 1: Sección transversal.- Variables utilizadas para determinar Ur y Ue
Q d 05 / 3 B (3)
El valor de expresado en forma general como una función de la profundidad media dm del flujo en la
sección antes de la socavación, la velocidad media en la sección Ur, y el caudal de diseño Qd:
1 1/ 2 5 / 3
Qd S d m Be (4)
n
donde:
Be : ancho efectivo en la sección, descontados todos los obstáculos, en metros. Para encontrar Be
se traza una línea perpendicular a las líneas de corriente, sobre esta línea se proyectan todos
los obstáculos y Be es la suma de todos los espacios libres. Así se toma en cuenta el
esviajamiento de la corriente (ver figura 2). Es muy importante conocer Be cuando existe una
obra en la sección de estudio. Cuando no existen obstáculos, Be = B.
dm : profundidad media de la sección que resulta de dividir su área hidráulica por el ancho efectivo
Be en metros: A/Be.
0
b1 b2
B Be3
Be1 Be 2
0
Be= Bei=Be1 +Be2+Be3
Planta
Figura 2: Determinación del ancho efectivo Be
Considerando la turbulencia que el flujo produce cerca de las pilas y estribos, es necesario aplicar en la
fórmula el coeficiente de contracción :
Qd S 1 / 2 d m5 / 3 Be (5)
n
Donde es un coeficiente que considera el efecto de contracción producido por las pilas y los estribos.
Se puede calcular como función de la velocidad media del flujo y el claro entre las pilas y estribos L
utilizando la siguiente expresión:
0.387U 0 Q
1 ; U0 d (5a)
L A
Al sustituir en la ecuación (5):
Qd d m5 / 3 Be (6)
despejando :
Qd
(7)
d m5 / 3 Be
Para una franja de ancho B , cuando d0 aumenta a cualquier valor de profundidad de socavación, ds,
la velocidad decrece a un nuevo valor Ur y Q se puede expresar como función de la velocidad y la
profundidad:
Q U r d s B (8)
Donde ds es la profundidad después de producirse la socavación del fondo, medida desde la superficie
libre del agua al pasar la avenida hasta el fondo de la sección erosionada. Sustituyendo en la ecuación
(3):
U r d s B d 05 / 3 B (9)
d 05 / 3
Ur (10)
ds
El valor propuesto en la ecuación (10) es válido mientras el ancho de la sección permanezca constante
durante el tiempo de la avenida, y el fondo descienda uniformemente en todo el ancho de la sección; la
primera condición generalmente se cumple mientras que la segunda ocurre sólo en los cauces de ríos
donde todo el ancho del cauce principal se encuentra cubierto de agua durante la época de sequía.
HIDRAULICA FLUVIAL 66 ING. M. Sc. CESAR MILLA VERGARA
Esta es una de las deficiencias del método ya que para cumplir completamente la condición establecida,
todos los tramos de la sección deben-presentar una resistencia similar a la erosión; incluso para todos
los estratos descubiertos. Se pueden mencionar dos ejemplos sencillos donde esta condición no se
cumple totalmente:
El primero es cuando el canal principal está cubierto de vegetación, excepto en la zona donde fluye el
agua con caudales bajos. Al presentarse la avenida se alcanzan mayores velocidades en la zona
desprotegida que en la zona con vegetación, por lo tanto será socavado con mayor rapidez
incrementándose también el caudal en esta zona.
El segundo ejemplo es cuando una sección está compuesta por una parte de arcilla y otra de arena.
Durante la avenida, la socavación será mayor en la parte arenosa ya que se requiere de mayor tiempo
para erosionar la parte arcillosa. Sin embargo, cuando sucede algo similar, al menos es posible
conocer que se tendrá una socavación mayor en las zonas más débiles, aún cuando no sea posible
evaluarla exactamente.
Una particularidad del primer ejemplo se presenta en los ríos donde el ancho del cauce durante la época
de sequía es menor que en la época de lluvias; aparecen islas y zonas descubiertas sobre las que crece
vegetación, lo cual incrementa la rugosidad y disminuye la velocidad. Como resultado se tiene que en la
época de avenidas la corriente principal tiende a fluir sobre el cauce de la época de sequía; cuando esto
sucede, la velocidad del flujo se incrementa en esa parte de la sección. En estas zonas, con velocidades
más altas, se produce también mayor socavación que la esperada.
ii) Cálculo de Ue
La velocidad mínima necesaria para arrastrar los materiales erosionados, Ue, depende de la naturaleza
de los mismos:
Si 0.05mm Dm 2.6mm
U e 4.7D84
0.03
0.28 0.322/ Dm
dS (11a)
Si 2.6mm Dm 182mm
U e 4.7D84
0.092
0.28 0.223/ Dm
dS (11b)
Si 182mm Dm 1000mm
U e 4.7D84
0.187
0.28 0.191/ Dm
dS (11c)
Donde:
: Coeficiente que toma en cuenta el periodo de retorno T, correspondiente al caudal de diseño.
Su valor se puede determinar con la fórmula siguiente:
0.8416 0.03342LnT
La velocidad Ue, calculada con la ecuación (11) es muy similar a la que se requiere para el inicio del
arrastre de las partículas.
Conocido el tipo de suelo que existe en el sitio y suponiendo que la rugosidad es constante en toda la
sección, la profundidad hasta la que llegará la socavación general se obtiene al igualar el valor de la
velocidad media Ur con el valor de la velocidad mínima necesaria para arrastrar los materiales
erosionados Ue.
d 05 / 3
4.7 D84
0.092
0.28 0.223/ Dm
dS (13b)
ds
d 05 / 3
4.7 D84
0.187
0.28 0.191/ Dm
dS (13 c)
ds
Si 0.05mm Dm 2.6mm
0.03 / 0.322 Dm
0.03
Dm
d 5/3
ds 0
(15 a)
4.7 D
0.28
m
Si 2.6mm Dm 182mm
0.092/ 0.223 Dm
0.092
Dm
d 5/3
ds 0
(15 b)
4.7 D
0.28
m
Si 182mm Dm 1000mm
0.187 / 0.191 Dm
0.187
Dm
d 5/3
ds 0
(15 c)
4.7 D
0.28
m
HIDRAULICA FLUVIAL 68 ING. M. Sc. CESAR MILLA VERGARA
Para suelos cohesivos:
0.725/ 66.28 m
0.725
m
5780d 5/3
ds 0
(16)
1.18
m
Se considera que existe rugosidad uniforme cuando la misma sección es utilizada por el flujo en las
épocas de estiaje y de avenidas por no haber cauces de inundación o bien cuando estos últimos son
pequeños y no crecen en ellos plantas o arbustos.
Las ecuaciones (15 y 16) se aplican en las franjas verticales en que se divide la sección (ver figura 3); en
cada una se obtiene una profundidad de socavación d S función de la profundidad inicial d0. Al unir todos
los puntos obtenidos se obtiene el perfil teórico de la sección erosionada.
P6
Sección transversal
P1
antes de la socavación P5
P2
R1
P4
R6
R2 P3
R5
Sección transversal
Para el cálculo de la socavación general se presupone una condición inicial de equilibrio, por lo que la
socavación calculada ocurre cuando se presenta el caudal de diseño.
Así también, si el tramo de río está sujeto a procesos de sedimentación y erosión el cálculo de la erosión
será válido sólo para el periodo de avenidas que siguen al periodo seco durante el que se realizó el
levantamiento topográfico de la sección transversal utilizada en el cálculo.
Un suelo heterogéneo está integrado por dos o más tipos de materiales muy diferentes entre sí,
conformados en capas. Puede consistir de un perfil estratigráfico de capas intercaladas de material
cohesivo y granular con diferentes densidades, distribuciones de diámetro, etc. No importa el tipo de
estratificación, la profundidad de equilibrio puede ser calculada aplicando las ecuaciones (15 y 16) en
forma analítica por el método de prueba y error o utilizando el método semigráfico.
De la distribución estratigráfica de los materiales bajo una vertical (ver Figura 4), se escoge el estrato
superior y de acuerdo con la naturaleza del material se aplica la fórmula correspondiente para determinar
ds. Si la profundidad obtenida resulta bajo el nivel del límite inferior del estrato, se escoge el segundo
estrato y se repite en forma similar el cálculo anterior, considerando la fórmula correspondiente al tipo de
suelo del segundo estrato. En el tanteo en el que la profundidad ds calculada esté dentro del espesor del
estrato en estudio, se habrá obtenido la ds buscada. El cálculo se debe hacer en forma ordenada,
iniciando por el estrato superior y continuando hacia los más profundos. Cuando la ds calculada se
HIDRAULICA FLUVIAL 69 ING. M. Sc. CESAR MILLA VERGARA
encuentra a un nivel superior al del estrato analizado, significa que ese material es muy resistente a la
erosión. En ese caso, la profundidad de erosión corresponde al nivel de la frontera superior del estrato
considerado. En todos los tanteos se utiliza el mismo valor de d0.
Sección transversal
antes de la socavación P6
P1 P5
P3 P4
P2 R6
R1 Arena
R5 Arena con grava
R2
Limo R4
R3 Sección transversal
después de la socavación
Pi:Puntos bajo estudio antes de la socavación
Ri:Puntos teóricos estimados de socavación
Sección transversal
Figura 4: Cálculo de la socavación para suelos heterogéneos- Método de prueba y error
- Método semigráfico
Para cada franja vertical se calculan las velocidades Ur y Ue a diferentes profundidades, por ejemplo las
correspondientes a los niveles de las fronteras superiores e inferiores de cada estrato. Los resultados
obtenidos se grafican en un sistema coordenado, con la velocidad sobre el eje horizontal y las
profundidades sobre el eje vertical. El punto de intersección de ambas curva determina la profundidad de
equilibrio de la socavación y la velocidad media correspondiente.
m Velocidad
P1 P6
Arena
Ue
P5
P4
P2 P3 ds O1
Ur
0
A II A2
A1 Ue
Ue
B3 B2
C3
m- 1 Profundidad
Graficas de Ue y Ur contra ds
Seccion transversal
para el punto P3
Figura 5 Cálculo de la socavación para suelos heterogéneos.- Método semigráfico
En el caso más general, cuando existen dos o más zonas con diferente rugosidad a lo ancho de la
sección en estudio, el procedimiento de cálculo es semejante, con la diferencia de trabajar en forma
aislada con cada franja y que para cada una debe calcularse el coeficiente a, correspondiente:
Qdi
i (17)
i d m5 / 3 Bei
Generalmente este problema se presenta cuando se tiene un cauce de avenidas extenso y cubierto de
vegetación, aunque también se presenta en caso de suelos heterogéneos. En la ecuación (17), el sub
índice i de las variables que intervienen en ella se refiere al valor para cada franja en que se divide la
sección en estudio.
ei i di
i 1
d i1 / 6
Ci (19)
ni
Donde:
3
Qdi : caudal que pasa por cada franja de la sección en m /s
Aei : área hidráulica efectiva antes de la erosión de la sección en estudio, es decir, el área hidráulica
total menos el área proyectada por los obstáculos en el plano perpendicular al flujo
3
Qd : caudal total de diseño en m /s
Ci : coeficiente de rugosidad de Chézy de cada franja
Di : profundidad media en cada franja en m
Ni : coeficiente de rugosidad de Manning de cada franja
vi) Efecto en la socavación cuando el flujo tiene una alta concentración de sedimento en suspensión
Cuando las condiciones aguas arriba favorecen el transporte de materiales finos en suspensión, existe
una reducción en la profundidad de socavación calculada para una determinada velocidad media. Esto
se puede atribuir a que se requiere cierta turbulencia para levantar las partículas del fondo. La
turbulencia, a su vez, es función de la velocidad del flujo dividida por la viscosidad cinemática del fluido.
Por lo tanto, para un cierto flujo con sólidos en suspensión es necesario incrementar la velocidad media
con el fin de obtener el mismo grado de socavación que con agua limpia. Lo anterior se logra
introduciendo el coeficiente como producto en las ecuaciones (11 y 12) este coeficiente depende del
valor del peso específico del agua con sedimentos suspendidos m y se puede calcular con la siguiente
expresión.
2
0.38 m (20)
1272
La profundidad de socavación general para material granular se obtiene utilizando las ecuaciones
siguientes:
Si 0.05mm Dm 2.6mm
0.03 / 0.322 Dm
0.03
Dm
d 5/3
ds 0
0.28 (21a)
4.7 Dm
Si 2.6mm Dm 182mm
0.092/ 0.223 Dm
0.092
Dm
d 5/3
ds 0
0.28 (21b)
4.7 Dm
Si 182mm Dm 1000mm
d 5/3
ds 0
0.28 (21c)
4.7 Dm
0.725/ 66.28 m
0.725
m
5780d 5/3
ds 0
(22)
1.18
m
Se puede calcular con el método para la socavación general, ya que se toman en cuenta las reducciones
producidas dentro del valor del área efectiva, Ae, y el ancho efectivo Be. Sin embargo, para tener una
idea aproximada del valor de la socavación transversal se puede utilizar la fórmula de Straub.
Partiendo de que el caudal sólido total por cualquier sección del tramo de río en estudio debe ser el
mismo, es decir:
GB 0 GB1 (23)
Donde GB es el caudal sólido transportado por el fondo y los sub índices 0 y 1 corresponden a las
secciones normal y reducida, respectivamente.
Al aplicar en la ecuación (11) para material granular, la fórmula de Meyer, Peter y Müller para calcular el
transporte de sedimentos, y se consideran las condiciones hidráulicas durante la ocurrencia de la
avenida, se obtiene la siguiente ecuación:
2/3
B
d mi d m 0 0 (24)
B1
donde:
dmi : profundidad media de la sección analizada
Bi : ancho efectivo de la sección analizada
El sub índice 1 corresponde a la sección reducida
El sub índice 0 corresponde a los valores en una sección inalterada localizada aguas arriba.
En el caso de que debido a las condiciones hidráulicas en ambas secciones las pendientes hidráulicas S,
sean diferentes, la fórmula anterior se modifica como:
2/3
S0 B0
d mi d m 0 (25)
S1 B1
La fórmula de Straub debe ser aplicada en fondos arenosos con una distribución homogénea del
material. Dada su simplicidad, es conveniente utilizarla para estimaciones preliminares. Es útil para
calcular la erosión en una sección reducida de ancho fijo, ya que se supone que el ancho de la superficie
libre permanece constante; es decir, que no se tiene erosión lateral.
El fundamento del método de Straub es suponer que el transporte de sedimentos entre una sección no
alterada aguas arriba y la sección de ancho reducido es constante.
Como acciones preventivas de una erosión intensa localizada en las secciones de los puentes, es
conveniente mantener el cauce limpio de vegetación en una longitud de al menos igual a una o dos
veces el ancho del cauce (2Be) hacia aguas arriba y a la mitad (Be/2) hacia aguas abajo.
HIDRAULICA FLUVIAL 72 ING. M. Sc. CESAR MILLA VERGARA
e) Estimación de la Socavación en Curvas
Existen varias formas de resolver este problema. Si se dispone de la sección transversal en estiaje se
puede aplicar el método de la socavación general, ya que éste indicará mayores profundidades en el
lado exterior de las curvas, además, su aplicación presenta la ventaja de que permite evaluar la
profundidad máxima y la forma aproximada de la sección transversal ya socavada.
Si no se tiene la sección de estiaje, las profundidades media y máxima se calculan a partir de las
características de la curva en planta, tales como el radio de curvatura r, medido al centro del cauce, y el
ancho de la superficie libre B. Maza Álvarez propone la siguiente fórmula basada en datos
proporcionados por Altunin:
2.24
B
d e 1.9 dm (26)
r
Donde:
dm : profundidad media en el tramo recto situado aguas arriba de la curva, en metros
dc : profundidad media a lo largo de la zona erosionada de la curva, en metros
Además, en un momento dado y en una de las secciones de la curva, la profundidad máxima, d cmáx.,
puede ser determinada con:
r
Si 2.0 5.91 :
B
0.38r
d c max 3.73 d m (27)
B
r
Si 5.91 :
B
B
d c max 1.27 1.26 d m (28)
r
Existen otros métodos para evaluar las profundidades medias en curvas y meandros. La máxima
profundidad puede ocurrir en cualquier zona a lo largo de la curva, pero no en toda al mismo tiempo. Con
mayor frecuencia aparecen aguas abajo del final de la curva. El valor obtenido con la ecuación (26) toma
en cuenta la socavación general probable.
Existen dos casos de socavación local que son de interés: cuando se produce al pie de obstáculos
rodeados por la corriente y por obstáculos pegados a una margen y que su función es desviar la
corriente. En el primer caso se pueden incluir las pilas de puente y en el segundo, los espigones y
estribos.
Este caso ha sido estudiado por muchos investigadores, entre los que se encuentran: Jarocki (1961),
Maza y Sánchez (1963), Breusers (1967), Shen-Shneider y Karaki (1969) y Melvilla (1976) entre otros.
Maza y Sánchez obtuvieron los tres diagramas mostrados en las figuras (6), (7) y (8), mediante los
cuales es posible determinar la socavación local como una función de la relación entre el ancho de la pila
y la profundidad del agua y el cuadrado del número de Froude.
b
U
b a
b1
U
b1=b;fc=1 Direcciones
U Donde:
S1
5.5
8
5.0
7 4.5
(d+s)/b1=(ST)/b1
4.0
6
3.5
d/b1
3.0
U
2.5
4 2.0
1.5
3
1.0
2
0.5
0
0.1 0.15 0.2 0.25 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
F fc F
Para =0º y F>0.25
Arena negra diam. = 0.17mm. Arena cafe diam. =0.56mm. Arena rosa, diam. =1.30mm.
Angulo de ataque = 15º Angulo de ataque =30º
b1
U
b1=b ; fc=1
U
d
Sr
Donde:
b1 es la proyeccion de la pila sobre el
plano normal al flujo.
S
99
8 5.5
d/b1
5.0
5.5
7
4.5
5.0
(d*s)/b1=Sr/b1
6 4.0
4.5
3.5
d/b1
4.0
5 3.0
3.5
2.5
4
3.0
2
2.5
3 1.5
2.0
1.0
2
1.5
0.5
1.0
1
0.5
0
0.1 0.15 0.2 0.25 0.3 0.4 0.5 0.6 0.7 0.8 0.91.0
f fc F
para 0 =0° y F >0.25
Arena negra diam =0.17mm o Arena cafe diam =0.56mm Arena rosa diam =1.30mm
Angulo de ataque de 0 =15.7 Angulo de ataque de 0 =30°
U
d
St
8 5.5
5.5
5.0
7 5.0 4.5
4.0
(d*s)/b 1=Sr/b 1
4.5
6
4.0
d/b1
3.5
3.5 3.0
5
3.0 2.5
4
2.5 2
1.5
3 2.0
1.5
1.0
2
1.0
0.5
1
1.5
0
0.1 0.15 0.2 0.25 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
d : tirante de la corriente frente a la pila en una zona no afectada por la erosión local, obtenido
después de calcular la socavación general, transversal o en curvas. Es decir, las erosiones locales se
calculan a partir de las condiciones de la corriente, una vez que se han producido las erosiones que
afectan al fondo del cauce.
La erosión local que se obtiene con ayuda de estas figuras es la máxima que puede presentarse.
Este método es útil solamente para fondos compuestos de arena y grava. Los diagramas han sido
desarrollados para pilas con secciones de formas rectangulares, rectangulares con frentes redondeados,
y circulares. Son también útiles cuando el eje de las pilas forma un cierto ángulo con la dirección del
flujo. Para pilas esviajadas, el número de Froude debe corregirse por el factor fc, si éste es mayor que
0.06; el valor de fc se indica en las figuras 6 y 7 como una función de (ángulo que forman los ejes de
las pilas con la dirección del flujo).
Para aplicar el método no es necesario conocer el mecanismo de socavación desarrollado, sin embargo
es conveniente conocer el proceso con la finalidad de tener una idea clara de la formación y el desarrollo
de la socavación al pie de las pilas de un puente.
Las siguientes observaciones fueron realizadas en un canal del laboratorio de hidráulica del Instituto de
Ingeniería de la Universidad Nacional Autónoma de México (UNAM). Los parámetros que influyen en la
profundidad de socavación al pie de una pila de puente se pueden agrupar, de acuerdo con su
naturaleza, en cinco diferentes clases:
a) Parámetros hidráulicos
1. Velocidad media del flujo.
2. Profundidad del flujo en la cara de la pila.
3. Distribución de velocidades.
4. Dirección del flujo con respecto a los ejes de las pilas.
e) Parámetros de tiempo
17. Duración del pico de la avenida.
18. Tiempo requerido para remover el material y alcanzar una condición estable.
Debido a la forma en que se realizaron los estudios, se despreciaron los factores 10, 14, 16, 17 y 18. Así
también se desarrollaron experimentos con tres tipos de material granular, cuatro tipos de pilas y
diferentes ángulos de incidencia; en cada caso, para un determinado material, tipo de pila y su
localización, sólo se varió la profundidad y la Velocidad del flujo.
Cuando el flujo alcanza una cierta velocidad media, la erosión se inicia en las dos esquinas de la cara de
aguas arriba y el material erosionado se deposita a los lados. Conforme se incrementa la velocidad, la
erosión en las esquinas continúa hasta alcanzar la misma profundidad a lo ancho de la cara de aguas
arriba (ver Figura 9). Al aumentar aún más la velocidad se inicia el transporte del fondo aguas arriba y se
deposita en la zona de erosión. Si el tirante del flujo se mantiene constante y se incrementa
paulatinamente la velocidad, la erosión aumenta hasta que se alcanza la estabilidad de socavación. El
material erosionado que al principio se deposita a los lados, es desplazado aguas abajo debido al
incremento de la velocidad; este desplazamiento persiste hasta que alcanza el otro extremo de la pila,
desde donde el material es transportado por el flujo. Con velocidades altas siempre habrá suspensión de
material en ese extremo, y en algunos casos otro depósito a poca distancia aguas debajo de la pila.
Conforme la profundidad de erosión varía, la estabilización del proceso es interferida por la formación de
rizos y dunas. Cuando la cresta de una duna se aproxima a la zona de erosión, la profundidad se
incrementa, alcanzando un máximo cuando la cresta se encuentra en la frontera de la zona de
socavación. Este valor se mantiene constante hasta que media duna ha sido desplazada aguas abajo,
decreciendo después hasta que pasa la siguiente duna. Los experimentos realizados en el Instituto de
Ingeniería mostraron que se requerían de dos a cinco horas para obtener el equilibrio de socavación.
Zonas de deposito
La erosión se inicia en la esquina exterior de la cara de aguas arriba ya medida que aumenta la
velocidad del flujo se hace uniforme frente a esa cara. Para velocidades altas con mucho transporte de
sedimentos, la máxima socavación ocurre en la esquina de aguas abajo opuesta a aquella en la que se
inició, es decir, en la esquina de aguas abajo expuesta al flujo y no protegida por el cuerpo de la pila. En
una pila enviajada, las magnitudes de las erosiones pueden ser dos o más veces mayores que para pilas
alineadas, por lo que es conveniente evitar esta condición. Por otro lado, como la máxima erosión se
presenta aguas abajo de la pila, se observa que cuando en grandes avenidas las pilas esviajadas fallan,
se inclinan hacia aguas abajo; esto se interpreta erróneamente, al suponer que el empuje de la corriente
es lo que las hace fallar.
En estos casos la erosión se desarrolla en forma similar a la de las pilas rectangulares, pero con ciertas
peculiaridades: la erosión se inicia en dos regiones aproximadamente a 65° a cada lado del eje de la pila
(ver Figuras 10 y 11). El resto del proceso es el mismo, excepto por las diferencias registradas de las
velocidades a las cuales los dos conos de erosión llegan a juntarse, ya las que el mismo grado de
erosión ocurre en la pila y lleva a la máxima erosión. Aguas abajo de las pilas circulares, las velocidades
se incrementan en un grado mayor que en las rectangulares, formando un vórtice con gran cantidad de
material en suspensión, Cuando no se puede garantizar que la dirección del flujo se mantendrá fija,
conviene utilizar pilas circulares.
Cuando se trata de pilas alargadas y esviajadas, la máxima socavación que puede ocurrir casi no
depende de la forma de sus extremos, por lo que la erosión obtenida con la figura 9 para pilas esviajadas
se aplica cualquiera que sea la forma de las mismas.
Zonas de deposito
E zona de A E
E deposito
D D D
C C C
zona de
zona de deposito deposito
Transporte
La misma profundidad Transporte
es alcanzada en toda la Esta zona aparece al
zona entre 65° medidos aumentar el flujo y esta
sobre cada lado del eje cargada de material en
Transporte
suspension
Transporte
zona de
deposito
zona de
18° deposito
Transporte zona de
deposito
Transporte
La socavacion maxima aparece en
la esquina C y es
considerablemente mayor que la
A A del frente de una pila alineada E
B E
65° B
D D E
C C E
Transporte
v) Comentarios generales
La erosión local al pie de obras rodeadas por la corriente se incrementa notablemente cuando la
estructura está esviajada, por lo que debe evitarse esta situación, y cuando se conozca de antemano
que la dirección del flujo pueda variar durante la vida útil de la obra, se deben seleccionar las pilas de
sección circular.
Si el fondo del río contiene boleos, la erosión local al pie de pilas es menor que la teórica que se obtiene
con el procedimiento descrito. No se conoce ningún método confiable para tener en cuenta esa
reducción.
Se desconoce la forma de evaluar la erosión local cuando el material del fondo es cohesivo.
Aunque este tipo de erosión ha sido también estudiado por varios autores, entre quienes pueden citarse
a Laursen, Shen, Veiga Da Cundha, es la más difícil de determinar y estudiar. El problema se da por la
cantidad de parámetros que es necesario tomar en cuenta y porque varían notablemente de "m punto a
otro, como es el caso del ancho total del estribo con respecto al primer claro del puente, el ángulo de
ataque de la corriente, el caudal teórico interceptado por el estribo, sus taludes y la forma e inclinación
de las esquinas cuando la cara del estribo es vertical, además de considerar la velocidad y el tirante de
la corriente y las propiedades del material del fondo. Para valuar esta erosión, los autores proponen las
relaciones que a continuación se indican y que fueron obtenidas a partir de datos de Artamonov y Veiga
Da Cundha, principalmente. La socavación al pie de un estribo, medida desde la superficie libre del
agua, está dada por:
d e K K k K q d a (28a)
Donde:
de : tirante del flujo en el extremo del estribo considerando la erosión local.
da :tirante del flujo aguas arriba del estribo en una zona no afectada por la erosión del estribo, pero
que incluye la erosión general, transversal, en curvas y otro tipo de erosión que afecte el fondo.
K :coeficiente que depende del ángulo ; su valor se obtiene con la relación:
K 0.782e 0.027
: ángulo medido aguas abajo del eje del puente y formado entre ese eje y la dirección del flujo.
Kk : coeficiente que depende del talud, k, del extremo del estribo; se obtiene con ayuda de la
expresión:
Si 0.001 k 1.3
0.85
Kk
k 0.0237
Si 1.5 k 3
k
Kk
2.8k 2.34
Kq : coeficiente que depende del cociente de Q1 y Q, donde Q1 es el caudal teórico que podría pasar
a través del área ocupada por el estribo si éste no existiera, y Q es el caudal total en el río.
d e 0.75K K k K q d a (28 b)
Aun cuando no existe forma de determinar la socavación cuando el espigón está cubierto de agua, y ya
que el método propuesto da los valores máximos posibles, es aconsejable considerar Q 1 como el flujo
teórico máximo interceptado por el espigón hasta su cresta.
El error máximo entre los resultados de las ecuaciones propuestas y los datos de Artamonov es menor
que 1.8 %.
HIDRAULICA FLUVIAL 81 ING. M. Sc. CESAR MILLA VERGARA
Pendiente
del estribo.
k1 k2
Q2
2)
1-Q
Q1 (Q
Q-
Flujo teorico interceptado Flujo teorico interceptado
por el estribo izquierdo. por el estribo derecho.
Q
Episto Total.
Figura 12: Esquema de la disposición de estribos. Fuente: Maza Alvares, México 1987
Para una compuerta colocada conforme a la figura 13, en el que la descarga es libre y el salto hidráulico
se forma inmediatamente después de la salida sobre un fondo erosionable, se puede determinar la
profundidad de la erosión utilizando el método propuesto por Valenti. Sin embargo, este método no
determina la distancia de la compuerta a la que ocurre la máxima erosión. La figura mencionada incluye
los valores de todos los parámetros necesarios para trazarla. Fr, es el número de Froude asociado con la
sección 1 y D90, es el diámetro representativo del material erosionado.
La erosión producida por los saltos de ski y deflectores deben ser estudiados mediante modelos
hidráulicos de fondo móvil, lo anterior es necesario debido a que el chorro es ocasionalmente continuo y
su forma y el aire mezclado en él, pueden variar considerablemente de una estructura a otra. El modelo
debe cumplir con la ley de similitud de Froude, ya que la escala de velocidades de caída de la partícula
es igual a la escala de velocidades del flujo. El problema en este tipo de análisis es reproducir la
cementación y el proceso de disgregación del material que constituye el área sobre la que golpea el
agua en su descarga. Si la obra existe y se tiene información de la erosión producida por una descarga
conocida, es posible calibrar el fondo móvil probando con diferentes grados de cementación hasta lograr
la representación del evento observado. Una vez que se logra la calibración es posible probar para
diferentes condiciones.
(a)
(b)
dk : Profundidad de socavacion causada por la
descarga de una compuerta en metros.
8
Fr1=U1/((gd1)^1/2)
0
0 2 4 6 8 10 12
(dk/d1)((D90/D1)^0.55)
(c)
Para cruzar un río, las tuberías de agua, oleoductos y gasoductos deben pasar bajo el fondo del cauce.
Al producirse la erosión general durante una avenida puede suceder que la tubería quede parcialmente
descubierta, lo cual induce una socavación local bajo la tubería cuyo valor puede estimarse con ayuda
de la figura 14 propuesto por Maza. Los parámetros adimensionales que intervienen en la figura son:
a/D, F, y S/D, donde “a” es la distancia del fondo del cauce a la parte inferior de la tubería (puede ser
negativa), en metros; D es el diámetro de la tubería en metros; Fr es igual a U / gd , donde U y d son
la velocidad y tirante medios de la corriente, en metros por segundo y metros, respectivamente, y S es la
socavación bajo la tubería medida desde el nivel original del fondo en metros.
Las curvas de la figura 14 dan el valor de la máxima erosión bajo la tubería. La socavación ocurre
cuando más de la mitad del diámetro de la tubería queda descubierta y el número de Froude es mayor a
0.1. Para evitar que la tubería falle, se debe emplazarla por debajo del nivel de socavación determinado
por el método de socavación general.
0.300
0.6
0.250
0.4
0.2
0.200
F1=0.100 0.150
0
-0.5 0 0.15 2 3 4
a/D
a.1) Generalidades
• Los diques en estudio son estructuras que evitan que el agua salga del cauce e inunde extensas áreas.
• Se tienen de tierra, enrocados con núcleo impermeable
• Los diques de tierra y enrocado tienen la ventaja de que son económicos y se aprovechan los
materiales del lugar
• Un dique de tierra de cualquier altura requiere de un ancho de base grande, y algunas veces no se
puede construir en ciudades, por el espacio que ocupan, se recurre en estos casos a los diques de
concreto
3
Descarga de Diseño (m /s) Borde Libre (m)
Menos de 200 0.6
200 – 500 0.8
500 – 2000 1.0
* PROYECTO DE LA CIMENTACION
Las cimentaciones se agrupan en tres clases principales de acuerdo con sus características
predominantes:
- Cimentaciones de roca
Este tipo de cimentación no presenta ningún problema de resistencia para los diques. Las
principales consideraciones son las peligrosas filtraciones erosivas y el excesivo flujo de agua por las
puntas, fisuras, hendiduras, estratos permeables, y a lo largo de los planos de falla.
Si en algún tramo se encuentra este tipo de problema, y en base a la importancia y características del
dique, se debe tomar en cuenta la inyección de lechada a presión para tapar hendiduras, juntas u
otras aberturas de la roca fija hasta una profundidad igual a la del máximo tirante de agua. Las
inyecciones pueden ser hechas con cemento puro y agua, en una relación 1:5. Se le añade arena o
arcilla, si se encuentran grandes huecos.
• Si no se toman las medidas adecuadas para controlar los asentamientos excesivos, puede ocurrir la
falla del dique:
- Por asentamientos diferenciales que producen la ruptura del núcleo impermeable del terraplén
- Por asentamiento del terraplén, por lo que se reduce el borde libre
Talud de Aguas
Clasificación de los Suelos Arriba Abajo
GW, GP, SW, SP No adecuado No adecuado
GC, CM, SC, SM 1V:2.5H 1V:2H
CL, ML 1V:3H 1V:2.5H
CH, MH 1V:3.5H 1V:2.5H
En la zona de saturación el material esta en suspensión debido a la presión hidrostática, afectado por las
fuerzas de filtración (presión hridrodinámica), que tiende a desplazar a las partículas pequeñas del
suelo en dirección hacia aguas abajo, originando el fenómeno de tubificación.
- Levantamiento o reventón
Es la separación o desplazamiento de un determinado volumen de suelo hacia la superficie, por acción
de las fuerzas de filtración.
- Filtros
a. Abatir la presión neutral en el agua que se infiltra en la cortina, con lo que se tendrá un
mejoramiento de la resistencia al esfuerzo cortante del material y de la estabilidad del dique
b. Un control del agua que se infiltra a través del dique, a la que se impide arrastrar el material fino
del mismo.
El dren puede estar compuesto de material granular (filtro), el cual debe satisfacer las siguientes
condiciones:
Análisis de Estabilidad
Mediante el análisis de estabilidad se verifica los taludes asumidos en el dique, en base a las
recomendaciones dadas en función del tipo de suelo usado
Un método muy usado para los análisis de estabilidad es el de Bishop. Este método considera una
superficie de falla circular, la masa de falla se divide en una serie de dovelas verticales y se
considera el equilibrio de cada una de ellas.
- Revestimiento de protección.
. Debido a que los diques se construyen generalmente de grava, arena, arcilla, es necesario protegerlos
de la erosión y de las filtraciones de agua.
. Se pueden usara, entre otros materiales, los enrocados de protección, o los colchones de gabiones.
. Se debe proteger también el pie del dique contra los fenómenos de socavación.
b.1) Generalidades
Los espigones son obras transversales que avanzan desde la orilla existente hasta la nueva línea de
orilla, para reducir las anchuras excesivas del lecho, provocando la sedimentación de la zona limitada
por ellos.
Las ventajas de los espigones sobre las obras longitudinales estriban en que permiten modificar
posteriormente con facilidad el ancho del lecho.
• Pero en cambio ofrecen el inconveniente de no formar un lecho regular, ya que constituyen
solamente guías a relativamente grandes distancias, y entre ellas el lecho se forma libremente
• La razón principal de que se empleen a pesar de este inconveniente, radica en que implican
menos desembolsos
• En ríos de montaña con considerable transporte de sedimentos los espigones no dan buen resultado
• Estas estructuras han sido muy efectivas en ríos pequeños y medianos.
• Los espigones se usan comúnmente en protección de tramos en curva
Los espigones se construyen con plataformas sumergidas, con cestones, con enrocados, con gaviones,
con cajones de piedras, de tierra con taludes protegidos, de concreto, etc.
Los espigones impermeables son construídos de piedra, grava y roca, como diques de tierra con
protección en los taludes, o como estructuras de concreto.
• Alrededor del morro ocurren fenómenos de socavación local, por lo tanto se debe hacer mucho más
seguro la protección tanto en el talud, como al pie de este.
• Los espigones sumergidos necesitan también tener protección en la cresta.
• El talud en el morro es mayor que los taludes del frente y la espalda del espigón
a) Espigones no sumergibles
Son los más efectivos causan mayores profundidades en los fenómenos de socavación local
• Al proteger una sola curva o un tramo completo, los primeros tres espigones aguas arriba deben tener
longitud variable: el primero será el de menor longitud posible (igual al tirante) y los otros dos
aumentarán uniformemente, de tal manera que el cuarto ya tenga la longitud del proyecto. La
pendiente longitudinal de la corona debe ser uniforme en todos ellos
La longitud total de un espigón se divide en una longitud de anclaje o empotramiento (L A), y en una
longitud de trabajo (LT)
Los valores de B, y del tirante deben ser los correspondientes al gasto formativo, aquel gasto que de
permanecer constante a lo largo del año, transportará la misma cantidad de material de fondo que el
hidrograma anual.
b) Otras recomendaciones
Sp = 4LT a 4.5LT
Sp = B a 2B
Separación en curvas
- Si la curva es regular, la separación recomendable es:
Sp = 2.5LT a 4LT
b) Espigones sumergibles
En este tipo de espigones la cresta es ascendente hacia arriba. Si el morro del espigón queda
por debajo del nivel de estiaje, se le denomina espigón bañado, y si toda la cresta se encuentra
por debajo de dicho nivel, se dice que el espigón es de tipo sumergido.
Se inician a la elevación de la margen, o a la elevación de la superficie libre al ocurrir el gasto
formativo.
El extremo dentro del cauce puede tener alturas máximas de 0.5 m. Sobre el fondo actual
La pendiente de la cresta puede estar entre 0.05 a 0.25
Se recomienda una pendiente S uniforme hasta el fondo. El piso C de los espigones debe
construirse primero para evitar erosiones locales durante la construcción.
b.10) Taludes
b) Talud de Morro:
El talud del morro es mayor que los taludes laterales, varían entre 1:2.5 a 1:5
- El revestimiento tiene por objetivo proteger el material de las riberas contra los efectos de la erosión.
- El revestimiento se apoya sobre un talud estable.
- Se realiza con los gabiones tipo “Colchón”, también conocidos como “Reno”
Procedimiento de Cálculo:
1.- Predimensionamiento de colchones.
2.- Verificación de arrastre de colchones por velocidades.
3.- Verificación de arrastre de colchones por esfuerzos cortantes.
4.- Verificación de deformaciones en el colchón.
5.- Verificación de estabilidad de material base por flujo en espacio colchón-filtro.
6.- Diseño de protección de pie de talud.
Velocidad Crítica: Velocidad que puede soportarse sin inicio de movimiento de piedras en colchón.
Velocidad Límite: Velocidad que puede soportarse admitiendo modestas deformaciones debido al
movimiento de piedras en colchón
B.- Verificación:
Si V > Vc habrá arrastre de colchones,
Se debe incrementar espesor del colchón.
Ta K . .h.S
donde:
Ta: Esfuerzo cortante máximo
:Peso específico del Agua
h: Tirante máximo
K: Coeficiente que depende del talud.
donde:
2
Tc : Esfuerzo Cortante Crítico en la orilla (kg/m )
C’* : Coeficiente de Shields (0.10)
s
3
: Peso Específico de la piedra del colchón (2500 kg/m )
w
3
: Peso Específico del agua(1000 kg/m )
Dm : Diámetro Medio de la piedra del colchón (m)
K = (1-sen /sen )
2 2 1/2
donde:
: Angulo que forma el talud con la horizontal
: Angulo de reposo de la piedra del colchón (41°)
C.- Verificación.
Si Ta > Tc el revestimiento será inestable
Se debe incrementar espesor del colchón
C.- Verificación:
Si mín <dm, habrá peligro en material base
Se debe incrementar espesor de colchón
Espesor:
2
S> df50[1-(Ve/Va) ]/0.05 y S <0.15-0.20 m>
Tipo 2:
Para grandes profundidades de erosión general z >1.5 m, extender horizontalmente protección de pie
de talud en una longitud x > 1.5 z.
- Normalmente los taludes de apoyo son perfilados con pendiente de 1:1.5 y 1:2 según el tipo de
terreno.
- 1:1.5 - para terrenos vegetales de mediana consistencia.
- 1:2 para terrenos arenosos finos incoherentes.
CONSIDERACIONES BASICAS
Los enrocados de protección de riberas y de diques son una parte importante en los trabajos de
tratamiento de ríos, y sirven a los siguientes propósitos:
• Tratamiento del río para mantener el alineamiento de la ribera
• Protección de los terrenos adyacentes contra la erosión
HIDRAULICA FLUVIAL 98 ING. M. Sc. CESAR MILLA VERGARA
• Protección de los diques de defensa contra inundaciones
• Protección de estructuras, como puentes, barrajes, presas, etc
3
Densidad de diferentes tipos de materiales en Kgmasa/m
Material S ( S W ) / W
Arena, grava 2650 1.65
Concreto 2200 1.2
Concreto armado 2400 1.4
Concreto asfáltico 2300 a 2400 1.3 a 1.4
Granito 2500 a 3100 1.5 a 2.1
Basalto 3000 a 4000 2a3
i) Fórmula de Maynord:
d50 V
C1F 3 F C2
y gy
Donde: d 50 es el diámetro medio de las rocas, y los valores recomendados de C1 y C2 se muestran a
continuación:
V 8.8 y
0.75Log
gd d
Donde:
V = velocidad media del flujo
d = diámetro mínimo de las rocas
y = tirante medio
0.001V 3
I
d50 Sistema Inglés
y 0.5k11.5
K1 1 Sen 2 / Sen 2 0.5
2.12 FS
1.5
C0 CsgCsf Csg Csf
DR 11.5 1.2
Donde:
= es el ángulo de inclinación del talud
= es el ángulo de reposo del enrocado
Dr = densidad relativa
FS = factor de seguridad
En el siguiente cuadro se muestra valores del factor de seguridad FS:
Diámetros mínimos de las piedras de protección (cm), para un tirante igual a 1m (*)
1
V V1 y a donde a
2 y
Conocidos V e “y”, se despeja V1, y se regresa al cuadro para conocer el diámetro de las piedras.
Se tiene una ola, de altura H, longitud L, que impacta sobre un enrocado, avanzando hacia arriba, luego
el agua regresa a lo largo del talud.
Se tienen las siguientes fórmulas para el cálculo de la masa mínima de las rocas, de tal manera que
puedan resistir el ataque de las olas:
- Fórmula de Iribarren:
fH 3 r
M
3 Cos Sen
3
FILTROS:
• La estabilidad del revestimiento, en una ribera o en un dique, depende no solamente del tipo y
construcción del enrocado, sino que depende en gran medida del tipo y composición del filtro.
• El filtro protege al suelo de la erosión debido a la corriente de agua, ataque de olas, y de flujos
subterráneos; y evita el movimiento de las partículas finas que conforma el suelo protegido.
• Las curvas granulométricas del filtro y del material del suelo deberán ser más o menos paralelas
en el rango de los diámetros pequeños
• En cuanto a la granulometría del material del filtro se han hecho muchas investigaciones entre las
que tenemos:
- Simons y Senturk recomiendan que la granulometría de los filtros debe cumplir con las
siguientes ecuaciones:
Se sugiere que el espesor mínimo del filtro de grava sea la mitad del espesor del enrocado.
- Terzaghi recomienda que la granulometría de los filtros debe cumplir con las siguientes
ecuaciones:
b. Filtros de Geotextil
b1) Ventajas y desventajas
Los filtros sintéticos son otra alternativa con respecto a los filtros granulares.
Ventajas:
- La instalación es generalmente rápida y eficiente
- Son consistentes y tienen una calidad de material más confiable
- Son capaces de deformarse con el enrocado y permanecen continuos
Desventajas:
- Puede haber dificultad para colocarlos debajo el agua
- El desarrollo de bacterias dentro del suelo, o sobre el filtro puede alterar el comportamiento
hidráulico definido en las especificaciones de fábrica
Los geotextiles que se colocan debajo de los enrocados de protección se recomiendan que cumplan con
las siguientes especificaciones mínimas:
• Un buen contacto entre el geotextil y el suelo es esencial. Por esta razón la superficie de la ribera o
del dique debe ser una superficie lisa, libre de protuberancias, depresiones y lentes de material
suelto.
• Debe ser colocado suavemente, sin pliegues, de arriba hacia abajo
• Se debe tener mucho cuidado al colocar el enrocado, pues puede romper el geotextil.
• Si las rocas tienen aristas filudas se debe colocar una subcapa granular entre el enrocado y el
geotextil
• Despues de colocar el enrocado, el geotextil debe ser asegurado al pie de este, tal como se indica en
la figura, y anclado en la parte alta de la ribera o dique
• Por lo tanto se debe proteger la base del talud con enrocado. En la siguiente figura se muestra un
esquema de protección.
• Los extremos del tramo protegido deben tener un tratamiento adecuado, de tal manera que no haya
peligro de que los flujos de avenida tiendan también a discurrir por debajo de los enrocados,
haciéndoles fallar.
• En cada extremo debe haber un anclaje hacia el interior de la ribera. En la siguiente figura se muestra
una instalación típica de enrocado en los extremos del tramo protegido.