Redes de Flujo
Redes de Flujo
Redes de Flujo
INTRODUCCIÓN.
Las filtraciones dependen de la carga hidráulica que origina la trayectoria del agua,
dependiendo también de las características físicas del material. El flujo puede llegar a
erosionar y socavar los terraplenes de la presa o sus fundaciones, ocasionando
asentamientos, desprendimientos, etc.
ASPECTOS TEÓRICOS.
Para calcular el flujo de agua a través del suelo se requiere la aplicación directa de la
ley de Darcy. El flujo de agua a través de un suelo no solo es en una dirección, y no es
uniforme sobre toda el área perpendicular al flujo. En tales casos, el flujo del agua
subterránea se calcula generalmente usando gráficas llamadas redes de flujo. El concepto
de red de flujo se basa en la ecuación de continuidad de Laplace, que gobierna la condición
de flujo permanente para un punto dado de la masa de suelo.
2.2.1 Ley de Darcy. - En 1856, Henri Philibert Gaspar Darcy publicó una simple ecuación
empírica para la velocidad de descarga del agua a través de suelos saturados; se basa
principalmente en observaciones relativas al flujo de agua a través de arenas limpias y se
expresa por:
v=k*i (2.5)
Donde:
v = Velocidad de descarga
k = Coeficiente de permeabilidad [cm /s]
i = Gradiente Hidraulico
∂v x ∂v y
+ =
0
∂x ∂y
∂ϕ * ⎛ ∂h⎞
vx =k x ⎜ − ⎟ =
⎝ ∂x ⎠ ∂x
h ∂
v y = k y * ⎛⎜ −∂ ⎞ ⎟ = ϕ (2.10)
⎝ ∂y ⎠ ∂y
Donde:
kx y ky , son permeabilidades en las direcciones horizontales y verticales.
ϕ , es función potencial tal que:
ϕ = −k * h
En la figura 2.5 se muestra la familia de las curvas ortogonales con sus respectivos
nombres.
Para completar la construcción gráfica de una red de flujo, se deben dibujar las
líneas de flujo y equipotenciales de manera que las equipotenciales crucen a las de flujo
según ángulos rectos y que los elementos de flujo formados sean aproximadamente
cuadrados.
Las líneas de flujo (Nf) forman canales de flujo y dos líneas equipotenciales (Nd)
forman las caídas equipotenciales.
Nf
Q=k*H* (2.17)
Nd
Redes de flujo en el cuerpo de la presa. - El flujo de agua a través de presas de tierra
constituye uno de los casos de mayor importancia en la aplicación de la teoría del flujo
de agua a los problemas prácticos, ello se debe tanto a la importancia que la presa tiene
en si como estructura, como al hecho de que en este caso el problema del flujo
presenta características especiales a las que ha de dedicarse un estudio especial para poder
llegar a soluciones apropiadas.
El flujo de agua a través de una presa de tierra produce los siguientes efectos:
◘ Uno directo, de pérdida de agua, que suele ser menos importante y más fácil
de controlar o subsanar.
◘ Un estado de presiones internas con componentes opuestas al efecto
estabilizador del peso. Además, al estar mojados los materiales, disminuye
su cohesión y su resistencia al rozamiento, añadiéndose estos efectos al de
sus componentes desestabilizadoras de las presiones internas.
◘ El paso del agua a través de las zonas con materiales finos tiende a arrastrar
esas partículas, con el consiguiente peligro de erosión interna progresiva.
Este fenómeno se llama sifonamiento (piping).
La presa de tierra es, en muchos sentidos una región de flujo como otra cualquiera;
trazando su red de flujo para las condiciones de frontera que se tengan se podrá calcular el
caudal de filtración, los gradientes hidráulicos, las velocidades del agua en cualquier punto,
las presiones hidrodinámicas, las fuerzas de filtración, etc.
La particularidad del problema del flujo a través de presas radica en el hecho de que
son necesarios métodos especiales para lograr trazar la red de flujo, en la región de flujo
que es la presa de tierra, no se conoce a priori una de las fronteras, de modo que no se
satisface el pre-requisito básico para resolver el problema, el conocimiento de todas las
fronteras de la región para el trazo de la red de flujo.
Fig.2.9: Condiciones de frontera en el problema del flujo de agua a través de una presa de tierra
(Fuente: Mecánica de suelos III - Juárez Badillo)
cuando el talud aguas arriba de la presa está invertido ( α 〉 90º ), según se muestra en la
figura 2.10.
Fig.2.10: Condiciones de entrada de la línea de corriente superior en la presa de tierra
(Fuente: Mecánica de suelos III Juárez Badillo)
Las condiciones de salida de la línea de corriente superior, o sea el ángulo con que dicha
línea intercepta al talud aguas abajo de la presa en el punto 4 (figura 2.9), dependen del
ángulo que dicho talud forma con la horizontal, cuando el ángulo es menor o igual que 90º
como sucede con el ángulo α de la figura 2.9, la línea de corriente superior debe salir
tangente al talud aguas abajo, siendo el punto de tangencia el punto 4.
Fig. 2.11: Estudio de la condición de salida de la línea de corriente superior para α ≤ 90º
(Fuente: Mecánica de suelos III Juárez Badillo)
como sigue:
mn a
m4 = = (2.19)
cos β cos β
Y por lo tanto:
a
Δ h2 = m4 cos (90 − α ) = (2.20)
senα
cos β
Debe tenerse que:
Δ h1 = Δ h2 (2.21)
Por lo tanto:
senα
sen(α − β ) = (2.22)
cos β
senα
〉 senα (2.24)
cos β
La línea de corriente superior existe en todas las presas de materiales sueltos en las
que hay flujo. Solo en casos particulares el flujo puede producirse exclusivamente a través
de la cimentación.
Es de gran interés determinar la línea de corriente superior, pues con ella se define el
contorno de la red y se facilita la obtención de las líneas restantes.
2.3.1 Método de Albert Casagrande para la línea de corriente superior en una presa
de tierra α=180º.- Este método es sugerido por Albert Casagrande, en 1937, aplica una
solución para el caso ilustrado de la figura 2.13 común en presas de tierra y en el que el
ángulo α tiene el valor de 180º.
Fig. 2.13: Presa de tierra con drenaje horizontal α=180º
(Fuente: Adaptado de texto docente)
Este método consiste en trazar una parábola geométrica para corregirla con respecto
a los extremos de entrada y salida, como ya se describió para la determinación de la línea
de corriente superior.
◘ Ecuación de la parábola
Distancia foco a punto:
d 1 = x2 + y 2 (2.25)
De la gráfica 2.15:
d1 = d 2 (2.27)
Reemplazando los valores de d1 y de d2 en la igualdad anterior (ecuación
2.27), se tiene la ecuación de la parábola base:
2 2
x + y = x+s (2.28)
s = d 2 + h2 − d (2.29)
2 2 2 2
x + y = x +s +2*x* (2.30)
s
despejando:
y = s2 + 2 * x * (2.31)
s
δy s
= 2 (2.32)
δx s +2*x*s
cuando x = 0 :
dy dy
=1≅ (2.33)
dx ds
dy
q=k*y* (2.34)
dx
dy
Si x = 0 ; y=s ; =1
dx
q=k*S (2.35)
21
Tratado básico de presas - Eugenio Vallarino
Fig. 2.17: Presa de tierra sin drenaje cuando α ∠ 90º
(Fuente: adaptado de texto docente)
El caudal infiltrado
dy
q=k*y* (2.36)
ds
Integrando la ecuación anterior
∫ q * ds = ∫ k * y * dy
2
y
q*s=k* + Ctte
2
Reemplazando el punto B conocido de coordenadas s = s 0 , y = h (longitud DB )
en la anterior ecuación se tiene:
k*h
Ctte = q * s0 − 2 (2.37)
2
dy
En el punto de tangencia = sin α ; y = a * sin α
ds
2
q = k * a * sin α * sin α = k * a * sen (2.38)
α
2
q = k * a * sen α (2.39)
k 2
k * a * sen2 α = *(y + (2.40)
2
h)
2 * (s − s0)
Con la siguiente ecuación se puede plotear la línea de flujo, s hasta cualquier punto
podría nuevamente ser tomado aproximadamente como la distancia entre el pie D.
2
(y −h )
2
2
a * sen α = (2.41)
2 * (s − s 0)
Cuando s = a , y = a * sen α
a sen α − h
2 2
2
2
a sen α = (2.42)
2(a − s 0)
Despejando a :
h2
CF = a = s0 − s02 − 2 (2.43)
sen α
S es la suma de la longitud del arco de parábola BK más el tramo KF del paramento libre.
Para α ∠ 60º puede suponerse con suficiente aproximación que esa longitud es la
hipotenusa del triángulo F B Bº (figura 2.19):
s ≈ h2 + d 2 (2.44)
y resulta:
2 2 2
a = h2 + d 2 − d − h cot α (2.45)
línea de flujo pasa tangencialmente por el punto C, que esta sobre el talud de aguas abajo a
una distancia a del foco medida sobre el talud. Según la gráfica: FC = a
Esta construcción gráfica vale para α < 60º, lo que se cumple en todas las presas
homogéneas. Pero en las heterogéneas hay núcleos verticales o inclinados y drenes
interiores; la salida del agua al medio drenante puede dar α > 60º e incluso α > 90º. Puede
darse este caso en una presa homogénea con drenajes horizontales profundos o drenes
chimenea. Para estos casos Casagrande da un gráfico (figura 2.20) para determinar el
desplazamiento Δ a = KC del punto de salida respecto a la intersección de la parábola,
22
como es lógico, puesto que sería el caso de un dren horizontal en una presa homogénea.
22
Tratado básico de presas - Eugenio Vallarino
Fig.2.20: Valores de Δa en función de α
a + Δa
2.3.3 Método de Casagrande – Línea de corriente en una presa de tierra con pie de
roca 90º < α<180º.- Para dibujar la línea de corriente superior para el caso que se muestra
en la figura 2.21, primero debemos dibujar la parábola base que sigue el procedimiento
23
explicado anteriormente.
Fig.2.22: Parábola base para una presa de tierra con pie de roca cuando 90º ∠ α ∠180º
(Fuente: Adaptado de texto docente)
23
Adaptado de texto docente
La parábola base corta al pie de roca en el punto K , sirve de guía para la
construcción de la línea de corriente. La distancia del punto de corte K al foco F es
a+ Δ a . Es decir: FK = a +Δa .
En la figura 2.23 se muestra un detalle de los puntos sobre la cara aguas arriba del
pie de roca y de la caída de la línea de corriente superior. El punto C se encuentra sobre el
talud de aguas arriba del pie de roca, la línea de corriente cae perpendicular al punto C .
1. Que para pequeñas inclinaciones de la línea de corriente superior las líneas de flujo
pueden considerarse horizontales y consecuentemente, las líneas equipotenciales
como verticales.
2. Que el gradiente hidráulico es igual a la pendiente de la línea de corriente superior
en el punto de que se trate y es constante en cualquier punto de la vertical que se
dy dy
trace por aquél ( en lugar de ).
dx ds
Este método toma como consideración hasta cierto punto las condiciones de entrada y
salida de la línea de corriente superior (figuras: 2.10- 2.12).
q=k*y*
dy
dx
∫ q * dx = ∫ k * y * dy
y2
q*x=k* + Ctte (2.46)
La constante de integración (Ctte) puede obtenerse manejando las coordenadas del
punto conocido, por el que pasa la parábola. Reemplazando el punto B conocido de
coordenadas x = d, y = h en la ecuación 2.46 se tiene:
2
k*
Ctte = q * d − (2.47)
h2
y k* 2
q * x = k *2 + q * d − h (2.48)
2
2
dy
i= = tan α
dx
k 2
k * a * sin α * tan α = *(y −H (2.51)
2
)
2*(x−d)
Despejando a :
d 2
d − h
2
a= − 2 2 (2.53)
cos α cos α sin α
25
2.3.5 Método gráfico para determinar “a”.- L. Casagrande propuso en este caso una
solución gráfica de la ecuación 2.43 que se describe a continuación:
26
2.3.6 Método de Albert Casagrande.- Existe otro método para dibujar la parábola, que
es el cálculo del flujo laminar a través de los medios porosos de las presas de materiales
sueltos.
a = d 2 + h2 − d 2 − h2 cot 2 α
Fig. 2.27: Esquema de la línea de corriente superior
(Fuente: Presas de corrección de torrenteras y retención de sedimentos Luís Suárez Villar)
26
Presas de corrección de torrenteras y retención de sedimentos Luís Suárez Villar
◘ Se determina la distancia “a” utilizando la ecuación 2.46 quedando así
determinado el punto S.
◘ El segmento SR se divide en el mismo número de partes iguales en que se
dividió el segmento BR (4 partes en el caso de la figura) y los puntos se
numeran secuencialmente, siendo 0 el punto S y 4 el punto R.
◘ Se unen mediante rectas los numerados del segmento BR con el punto S.
Pendiente dy dy
tan α ( ) sin α ( )
dx ds
30º 0.577 0.5
60º 1.732 0.866
90º ∞ 1.00
Tabla.2.2: Errores de los métodos para diferentes ángulos
(Fuente: Adaptado de texto docente)
2.5 EJEMPLOS.
Ejemplo 5.1
La presa del dibujo se asienta sobre materiales cuya conductividad hidráulica es 0,3 m/día.
Bajo dichos materiales se encuentra un sustrato impermeable. Se pide:
Estrato Impermeable
Solución:
Se han dibujado tres líneas de flujo, pero eso no es regla, la red puede ser más densa
con la cual se obtendría una solución mas precisa. Pero sería más difícil o imposible de
realizar a mano.
Hemos de suponer que una gota procedente del infinito (a la izquierda del dibujo)
circula pegada al fondo y finalmente asciende fuera del dibujo por la derecha. Esto nos hace
considerar que además de todos los cuadros dibujados también aparecen dos "cuadros
abiertos", a la izquierda y a la derecha de la red dibujada. Por el contrario, el agua a ambos
lados de la presa es una línea equipotencial: todos los puntos del fondo del lago tienen el
mismo potencial. Por tanto, las líneas de flujo nacen y terminan perpendicularmente.
Aunque parezca que cumple las normas (huecos cuadrados, cortes perpendiculares)
hemos dibujado círculos inscritos en algunos huecos, observando que algunos pueden ser
aceptables (en verde, punteados), pero la mayor parte de ellos no son cuadrados (los de
color rojo, interior sin puntos).
Es muy difícil dibujar a mano (como es el caso) una red de flujo perfecta, pero para
mostrar cualitativamente el flujo y para el cálculo que realizaremos ahora, la precisión es
suficiente.
Y dentro de este tubo aplicaremos la Ley de Darcy a la sección de una presa de las
superficies equipotenciales, cuya ancho sería a.
Consideremos ese tubo aislado y apreciamos que la sección es igual a:
Área = a*1
La distancia entre esa equipotencial y la siguiente es Δx y la diferencia de potencial
entre ellas será:
H 6
Δh = = = 0,4metros
Nd 15
h
Q=k* * Nf
Nd
Q = 0.48 m3 día
Este seria el caudal por cada metro.
H 6
Δh = = = 0,4 metros
Nd 15
Solución:
GE = 67,5 m
1
BE = GE = 22,5 m
3
d = 111 m
1
tan α =
2,5
Obteniendo el valor de a con la ecuación 2.45:
2 2
d + h + d − h cot α ; a = 26,122 m
2 2 2
a=
−4
b) k h=9k , Esta condición nos dice que es un suelo anisotrópico k h = 9 x 10 y
v donde
−4
k v = 10
BE = 7.5 m
O B1 = d = 37 m
a = d 2 + h2 − d 2 − h2 cot 2 α
a = 16,55 m
Continuando con el procedimiento determinamos las coordenadas de los puntos B, 1, 2, S
B = (15 , 27)
1 = (27 , 22.2)
2 = (35.8 , 14.5)
S = (41.4 , 12.7)
Teniendo las coordenadas se vuelve a la escala original, para obtener la verdadera línea freática,
las coordenadas de los puntos por donde pasa la línea freática se mantendrá en su coordenada
vertical pero escalando en la coordenada horizontal.
Entonces las coordenadas transformadas a escala original
son:
B = (45 , 27)
1 = (81 , 22.2)
2 = (106 , 14.5)
S = (124.2 , 12.7)
Con estas nuevas coordenadas se procede a dibujar la nueva línea de flujo en las dimensiones
originales de la presa, se puede tomar más puntos de referencia para tener más exactitud en el
recorrido de la línea de flujo.
q=S*K
Se tiene:
-4
K= k h * k v = 3 X 10 cm seg
-6
K = 3 X 10 m
seg
2
área = S = 810 m
3
Entonces el caudal es: q = 2.43 x 10 – m seg
3