FLUJO BIDIMENSIONAL

MECÁNICA DE
SUELOS
PRIMER SEMESTRE 2023
Prof. Catalina Lozada

1
FLUJO BIDIMENSIONAL

FLUJO BI-DIMENSIONAL

• Ecuación de continuidad LAPLACE

• Redes de flujo
• Calculo de flujo bajo estructuras hidráulicas
• Fuerzas de subpresión bajo una estructura
• Análisis numérico de flujo
• Fuerzas de infiltración por unidad de volumen

2
FLUJO BIDIMENSIONAL

Ecuación de continuidad LAPLACE

3
 FLUJOBIDIMENSIONAL
FLUJO BIDIMENSIONAL
Ecuación de continuidad LAPLACE
𝑞𝑖𝑛 = 𝑞𝑜𝑢𝑡

𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 = 𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 + d𝑞𝑥 + 𝑑𝑞𝑦 + 𝑑𝑞𝑧

𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 = 𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 + d𝑞𝑥 + 𝑑𝑞𝑦 + 𝑑𝑞𝑧

0 = d𝑞𝑥 + 𝑑𝑞𝑦 + 𝑑𝑞𝑧

0 = d𝑣𝑥 𝑑𝐴𝑥 + d𝑣𝑦 𝑑𝐴𝑦 + d𝑣𝑧 𝑑𝐴𝑧

𝜕2 ℎ 𝜕2 ℎ 𝜕2 ℎ
0= 𝑘𝑥 𝜕𝑥 𝜕𝑥𝜕𝑦𝜕𝑧 + 𝑘𝑦 𝜕𝑦 𝜕𝑦𝜕𝑥𝜕𝑧 + 𝑘𝑧 𝜕𝑧 𝜕𝑧𝜕𝑥𝜕𝑦

𝜕2 ℎ 𝜕2 ℎ 𝜕2 ℎ
0= 𝑘𝑥 𝜕𝑥 𝜕𝑥𝜕𝑦𝜕𝑧 + 𝑘𝑦 𝜕𝑦 𝜕𝑦𝜕𝑥𝜕𝑧 + 𝑘𝑧 𝜕𝑧 𝜕𝑧𝜕𝑥𝜕𝑦

𝜕2 ℎ 𝜕2 ℎ 𝜕2 ℎ
0= 𝑘𝑥 + 𝑘𝑦 + 𝑘𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

4
FLUJO BIDIMENSIONAL

𝜕2 ℎ 𝜕2 ℎ 𝜕2 ℎ
0= 𝑘𝑥 𝜕𝑥 + 𝑘𝑦 𝜕𝑦 + 𝑘𝑧 𝜕𝑧

𝑠𝑖 𝑘𝑥 = 𝑘𝑦 = 𝑘𝑧

𝜕2 ℎ 𝜕2 ℎ 𝜕2 ℎ
0= + + Ecuación de continuidad LAPLACE
𝜕𝑥 𝜕𝑦 𝜕𝑧
3D

𝜕2 ℎ 𝜕2 ℎ
0= + 2D
𝜕𝑥 𝜕𝑦

𝜕2 ℎ
0= 1D
𝜕𝑥

5
FLUJO BIDIMENSIONAL

Ecuación de continuidad LAPLACE para flujo 2D

Condiciones:
• Flujo permanente
• Agua incompresible y no hay cambio de volumen
• Suelo isotrópico con respecto a la permeabilidad

𝜕2ℎ 𝜕2ℎ
+ =0
𝜕𝑥 2 𝜕𝑧 2

Para flujo en 2D y kx = ky

6
FLUJO BIDIMENSIONAL

Ejemplo para un problema de flujo 1d

𝜕2ℎ
=0 ℎ = 𝐶1𝑥 + 𝐶2
𝜕𝑧 2
Solución analítica

• Solución exacta de la ecuación de Laplace para problemas simplificados.

Geometría y condiciones de frontera simples.

• El flujo es únicamente en una dirección.

• Encontrar las ecuaciones por medio de la determinación de las

condiciones de borde.

7
FLUJO BIDIMENSIONAL

Ejercicio No 1. Ecuación de continuidad LAPLACE flujo 1D

• El flujo ocurre únicamente en una
dirección. En el eje x.
• La longitud y los coeficientes de
permeabilidad de las dos capas son
conocidas .
• Las alturas de las cabezas totales en las
secciones 1 y 3 son conocidas.

Se debe calcular:

a) Cabeza total en 2 h2
b) Una función para la cabeza total entre x
=0 y x = La
c) Una función para la cabeza total entre
x= La y x = La + Lb

8
FLUJO BIDIMENSIONAL

Ejercicio No 1. Ecuación de continuidad LAPLACE flujo 1D

solución

a) Cabeza total en 2 h2

Altura total h2 𝑘𝐴 ℎ1
ℎ2 =
𝑘 𝑘
𝐿𝐴 (𝐿𝐴 + 𝐿𝐵 )
𝐴 𝐵

9
FLUJO BIDIMENSIONAL

Ejercicio No 1. Ecuación de continuidad LAPLACE flujo 1D

solución
a) Una función para la cabeza total entre x =0 y x = La
b) Una función para la cabeza total entre x= La y x = La + Lb

Para flujo unidimensional: 𝜕2ℎ

=0
𝜕𝑥 2

Integrando dos veces: ℎ = 𝐶2 𝑥 + 𝐶1

Condiciones de frontera para el flujo a través del suelo A

1) X = 0 ; h=h1
2) X = L ; h = h2
Condiciones de frontera para el flujo a través del suelo B
1) X = LA ; h = h2
2) X = LA + LB ; h = 0
10
FLUJO BIDIMENSIONAL

Ecuación de continuidad LAPLACE

Solución gráfica: Redes de flujo

Solución numérica: Diferencias finitas (Matlab, Excel)

Solución numérica: Elementos finitos (Plaxis, Abaqus, Rocsience)

11
FLUJO BIDIMENSIONAL

Ecuación de continuidad LAPLACE Solución gráfica: Redes de flujo

https://www.youtube.com/watch?v=rM38JiyXDU8

12
FLUJO BIDIMENSIONAL

Redes de flujo– Líneas de flujo y equipotenciales

La solución de la ecuación de
continuidad representa dos familias
de curvas ortogonales: líneas de
flujo y líneas equipotenciales.

Líneas de flujo: Línea a través la

cual el agua viaja de un lado aguas
arriba a uno aguas abajo en un
medio permeable

Líneas equipotenciales: Es una

línea la cual la cabeza de potencial
en todos los puntos es la misma. Si
varios piezómetros son ubicados en
las líneas equipotenciales todos
marcaran la misma altura.

13
FLUJO BIDIMENSIONAL

Redes de flujo

Combinación de varias líneas

de flujo y líneas equipotenciales.

Construidos para calcular el

flujo de agua subterránea y las
cabezas en el medio.

Para graficar estas líneas se

considera:

1. Las líneas equipotenciales

intersectan las líneas de flujo en
ángulos rectos.

2. Los elementos de flujo

formados son aproximadamente
rectángulos o cuadrados.

14
FLUJO BIDIMENSIONAL

Redes de flujo
• Nf es el numero de canales de flujo • ∆H es la pérdida de cabeza total en el
• Nd es el numero de caídas de sistema H1 - H2
potencial.

Condiciones de frontera:

• Condición 1: Las superficies

permeables aguas arriba y aguas NR
abajo son líneas equipotenciales b-a
y d-e 6
• Condición 2: Debido a que estas
1
líneas son equipotenciales. Las líneas
2 3 5
de flujo las interceptan en ángulos 4
rectos.
• Condición 3: El borde de la capa
impermeable: línea f-g es una línea
de flujo.
• Condición 4: Las líneas
equipotenciales intersectan acd y fg
en ángulos rectos.

15
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

El canal entre dos líneas de flujo

adyacentes se llama canal de
flujo.

h1, h2, h3, h4 son los niveles

piezometricos correspondientes a
las línea equipotenciales.
El caudal a través de cada canal
de flujo se puede calcular como:

16
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

De la ley de Darcy q =kiA. De la

ecuación anterior:

17
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

Si los elementos de flujo se dibujan como aproximadamente cuadrados, la caída

de la altura piezométrica de dos líneas equipotenciales adyacentes se llama
caída de potencial.

= ∆h

18
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

Si el numero de canales de flujo es: Nf

El caudal total en todos los canales esta dado por:

19
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

Si los elementos de flujo son

rectangulares.

• Cuál será el caudal para cada

canal de flujo?

• Cuál será el caudal total para

todos los canales?

20
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

Flujo confinado Flujo inconfinado

21
FLUJO BIDIMENSIONAL
• A partir de la red de flujo mostrada en la siguiente figura calcular:
• La cabeza total en los puntos a, b, c y d.
• Los esfuerzos efectivos y las presiones de poros en los mismos puntos. Asumir que el peso
unitario saturado para la arcilla es γsat= 21 kN/m3.
• Distribución de presiones en ambos lados de la tablestaca.
• El caudal por unidad de longitud en el canal de Flujo II si b/l=1.
• El gradiente hidráulico a la salida iexit y el gradiente hidráulico crítico.

22
FLUJO BIDIMENSIONAL

23
FLUJO BIDIMENSIONAL

Factor de seguridad contra tubificación - piping

Gradiente hidráulico crítico

∆ℎ 𝛾 ′
=
𝐿 𝛾𝑤

Fs = icr/ iexit

24
FLUJO BIDIMENSIONAL
Flujo confinado
Fuerzas de subpresión bajo una estructura

El calculo de la red de flujo puede utilizarse para determinar las subpresiones en

la base de una estructura hidráulica
Asumir kx = kz = k Determinar el diagrama de distribución de presiones en la base
de la presa.

25
FLUJO BIDIMENSIONAL
Flujo confinado
Ejemplo No. 2: Calcular la fuerza de subpresión bajo la presa

1. Calcular la cabeza de posición, total y de presión en a-f.

2. La fuerza es el área del diagrama de presiones.

26
FLUJO BIDIMENSIONAL
Flujo confinado

Nd = 7
∆H = 7m
Δh = 7 / 7 = 1m

27
FLUJO BIDIMENSIONAL
Flujo confinado
Ejemplo No 3

a) Calcular las cabezas total, posición y presión en los puntos A-E

b) Calcular la presión del agua en el punto C

28
FLUJO BIDIMENSIONAL
Ejemplo No 4:
Flujo confinado

a) Determinar el caudal total

si k = 20x10-4 cm/s.
Longitud de la presa de
120 m.
b) Gradiente de salida en x
c) Distribución de presiones
0.5∆h 0.4∆h en la base (3.5 – 6.9)

29
FLUJO BIDIMENSIONAL
Flujo confinado
Ejemplo No 5:

30
FLUJO BIDIMENSIONAL

Calculo de Redes de flujo

Flujo inconfinado

hp = 0
ht = hzb
Línea freática

0
ht = hp+hz
hzb ht = hz

Filtro para permitir

El drenaje

31
FLUJO BIDIMENSIONAL

PARA ESTUDIAR…

 Ejercicios en Moodle

Y Libros para estudiar:

- Lambe, T. William, and Robert V. Whitman. Mecánica de suelos. No. 624.151
LAM. 1972.
- Budhu, Muni. SOIL MECHANICS AND FOUNDATIONS. John Wiley & Sons,
2008.

32
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

 La solución se evalúa en un numero finitos de puntos (malla o grilla).

 Las derivadas se aproximan como expresiones algebraicas con los valores

de la variable dependiente en el numero de puntos seleccionados.

 El valor de los puntos seleccionados se convierten en incógnitas. El sistema

de ecuaciones algebraicas debe ser resuelto.

33
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

Problema: Solución de la ecuación de Laplace para flujo bidimensional en un

suelo isotrópico con respecto a la permeabilidad.

𝜕2ℎ 𝜕2ℎ
+ =0
𝜕𝑥 2 𝜕𝑧 2

En el caso de problemas de flujo confinados, la cabeza total o la velocidad del

fluido esta indicada en las fronteras.
En términos matemáticos:
Las condiciones de frontera se dan en cabeza total o gradientes de cabeza
total.

34
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

Discretización de las derivadas.

Una función continua puede
definirse en términos de
valores discretos fi
correspondientes a los valores
xi en el eje x.
Si la función es diferenciable
se puede expandir utilizando
las series de Taylor.

35
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

Discretización de las derivadas. Expansión de Taylor

Si x = xi

36
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

La ecuación de primer orden se puede obtener restando (4) y (5):

La derivada de segundo orden se puede obtener sumando (4) y (5)

Aproximaciones de la primera y segunda derivada. Los errores entre la

solución exacta y aproximada son proporcionales a Δx2

37
FLUJO BIDIMENSIONAL

Flujo bidimensional – Diferencias finitas

La ecuación de primer orden se puede obtener restando (4) y (5):

La derivada de segundo orden se puede obtener sumando (4) y (5)

Aproximaciones de la primera y segunda derivada. Los errores entre la

solución exacta y aproximada son proporcionales a Δx2

38
FLUJO BIDIMENSIONAL

Discretización de problemas en dos dimensiones.

Las anteriores ecuaciones aplican para flujo bidimensional. El espacio de dos

dimensiones se discretiza con una malla y unos nodos. Las coordenadas se
denominan i, j.
Δx es la separación en el eje x y Δy es la separación en el eje y

39
FLUJO BIDIMENSIONAL

Discretización de problemas en dos dimensiones.

𝜕2ℎ 𝜕2ℎ
La ecuación de Laplace discretizada en el punto i, j es: + =0
𝜕𝑥 2 𝜕𝑧 2

Como únicamente los nodos alrededor de i,j contribuyen en esta ecuación. Si

Δx = Δy:

Donde: α = kx / ky. Si kx = ky : α = 1

40
FLUJO BIDIMENSIONAL

Discretización de problemas en dos dimensiones.

Condiciones de Borde:
Para considerar fronteras de flujo se consideran únicamente las
fronteras impermeables. La velocidad de infiltración en la frontera es
cero.

De la solución discretizada de esta ecuación se obtienen diferentes

ecuaciones para las fronteras.

41
FLUJO BIDIMENSIONAL

Discretización de problemas en dos dimensiones.

Condiciones de Borde:

42
FLUJO BIDIMENSIONAL

PARA ESTUDIAR diferencias finitas…

Y Libros para estudiar:

- Budhu, Muni. SOIL MECHANICS AND FOUNDATIONS. John Wiley & Sons,
2008.
- Bardet, Jean-Pierre. Experimental soil mechanics. Prentice Hall, 1997.

43