5.- cÁlculo de flujo en la presa heterogÉnea
TRANSCRIPT
70
5.- CÁLCULO DE FLUJO EN LA PRESA HETEROGÉNEA 5.1.- GEOMETRÍA DE LA PRESA
El cálculo del problema de flujo se ha realizado a partir del caso B2, es decir,
para el caso de presa heterogénea calculada con la teoría de Mohr-Coulomb. A
los resultados obtenidos en este caso, se han añadido dos fases más
correspondientes al llenado de la presa y un desembalse rápido de la misma.
En ambos casos se ha partido de desplazamientos nulos, sin acumular los que
anteriormente se habían calculado para el caso B2. Por tanto, la geometría de
la presa y la malla de cálculo siguen siendo idénticas a las del caso B2.
5.2.- PROPIEDADES DE LOS MATERIALES
Las permeabilidades que se han tomado para el cálculo del flujo son las que se
resumen en la tabla 14:
FLAC PLAXIS
K (cm/s) K (m3·s/kg) K (m/s)
Material 1 Espladones 1.00E-03
1.00E-09
1.00E-05
Material 2 Núcleo 1.00E-05
1.00E-11
1.00E-07
Material 3 Roca 1.00E-05
1.00E-11
1.00E-07
Tabla 14. Permeabilidades de los materiales
En FLAC se ha de introducir los datos de la permeabilidad como una
permeabilidad dinámica, cuya equivalencia con la permeabilidad del material
es:
Permeabilidad dinámica (FLAC) ≡ Permeabilidad (cm/s) · 10-6
71
Para evitar errores de cálculo, se han escogido permeabilidades idénticas para
la roca y para el núcleo central, ya que permeabilidades muy diferentes dan
lugar a errores en la convergencia del cálculo del flujo. En FLAC, el efecto de
introducir permeabilidades muy diferentes (a partir de dos órdenes de magnitud
de diferencia), implica un aumento en el tiempo de cálculo considerable debido
a que, al ser un método de solución explícita, necesita realizar muchas más
iteraciones para llegar a la convergencia. En PLAXIS la tolerancia frente a
valores muy diferentes de la permeabilidad, aunque también puede llegar a
tener problemas en el cálculo cuando las permeabilidades son muy diferentes
(más de 6 ordenes de magnitud), pero a diferencia de FLAC, PLAXIS o bien no
converge, o bien no realiza el cálculo, pero en cualquier caso el tiempo que
necesita es muy inferior al que es empleado por FLAC.
Casos analizados
Caso C1: Cálculo de flujo para un embalse de altura de agua H = 28 m.
Caso C2: Cálculo de flujo para un desembalse rápido desde una altura de
agua H = 28 m, hasta H = 5 m.
5.3.- CASO C1. CÁLCULO DE FLUJO. EMBALSE DE LA PRESA
5.3.1.- CÁLCULO REAL DEL PROBLEMA DE FLUJO ACOPLADO. MÉTODO DE MOHR-COULOMB
En este caso se ha procedido a añadir una etapa más de cálculo a las 8 etapas
anteriores realizadas en el caso B2, tanto en FLAC como en PLAXIS. Para ello,
se han definido en ambos programas la altura de la superficie del agua en el
contorno del talud de aguas arriba de la presa, y se ha dejado la salida libre en
el contorno del talud de aguas debajo de la presa (ver figuras 30 y 31).
72
Figura 30. Condiciones de contorno para la determinación de la red de
flujo en PLAXIS.
Figura 31. Condiciones de contorno para la determinación de la red de
flujo en FLAC.
5.3.2.- RESUMEN DE LOS RESULTADOS OBTENIDOS
En este caso analizaremos los resultados de la superficie piezométrica dentro
de la presa, la distribución de las presiones de agua y el efecto sobre las
tensiones totales y los desplazamientos que produce el embalse de la presa.
73
Superficie piezométrica En las figuras 32 y 33 se presentan los resultados obtenidos por FLAC y por
PLAXIS en el cálculo de la superficie piezométrica.
Figura 32a. Vectores de flujo en FLAC 2 0 .00 0 30 .0 0 0 4 0 .0 00 5 0.0 0 0 60 .0 0 0 7 0 .00 0 80 .0 0 0 9 0 .0 00 1 00 .0 0 0 11 0 .0 00 1 2 0.0 0 0 13 0 .0 00
-2 0 .0 00
-1 0 .0 00
0 .0 00
1 0 .0 00
2 0 .0 00
3 0 .0 00
4 0 .0 00
F low Fie ldE xtr em e vel oc ity 5 .5 5 *1 0-6 m /s To tal d is ch ar g e 4 .8 4* 10 -6 m 3 /s/m
Figura 32b. Vectores de flujo en PLAXIS
Figura 32. Vectores de flujo obtenidos por ambos programas. El valor máximo del vector de flujo es mayor en PLAXIS que en FLAC.
74
Figura 33a. Red de flujo obtenida por FLAC.
CDEFGHIJKLMNO
P
20.000 30.000 40.000 50.000 60.000 70.000 80.000 90.000 100.000 110.000 120.000
-20.000
-10.000
0.000
10.000
20.000
30.000
40.000
50.000
Active groundwater headExtreme groundwater head 28.00 m
[ m]
A : -2.000
B : 0.000
C : 2.000
D : 4.000
E : 6.000
F : 8.000
G : 10.000
H : 12.000
I : 14.000
J : 16.000
K : 18.000
L : 20.000
M : 22.000
N : 24.000
O : 26.000
P : 28.000
Q : 30.000
Figura 33b. Líneas equipotenciales en PLAXIS.
Figura 33. Las líneas equipotenciales del núcleo son muy parecidas en los dos programas. La cota de entrada del agua al núcleo es de H = 28 m, y la de salida del núcleo es de H = 6.5 m, aproximadamente.
75
Presiones de agua
La distribución de las presiones de agua en la presa se muestran en la figura
34, y en la tabla 15 se resumen los valores de la presión de agua tomados en
el eje central del núcleo de la presa, en diferentes alturas.
Tabla 15. Presiones de agua en el eje de la presa
Figura 34a. Presiones de agua en FLAC
Pw(KN/m2)
H(m)
FLAC PLAXIS 28 0.00 0.00 24 12.60 16.48 20 36.20 43.08 16 57.40 65.56 12 78.30 87.69 8 101.60 112.96 4 129.50 143.92 0 162.00 179.15 -3 217.00 208.90 -6 241.20 236.99
76
D E FG H I J K L MNOP QRS
20.000 30.000 40.000 50.000 60.000 70.000 80.000 90.000 100.000 110.000 120.000
-20.000
-10.000
0.000
10.000
20.000
30.000
40.000
50.000
Active pore pressuresExtremeactivepore pressure -340.00 kN/m2
(pressure= negativ
e
)
[ kN/m2]
A : -360.0B : -340.0C : -320.0D : -300.0E : -280.0F : -260.0G : -240.0H : -220.0I : -200.0J : -180.0K : -160.0L : -140.0M : -120.0N : -100.0O : -80.00P : -60.00Q : -40.00R : -20.00S : 0.000T : 20.00
Figura 34b. Presiones de agua en PALXIS
Figura 34. Distribución de las presiones de agua. Aparentemente, las presiones son
idénticas en ambos programas, pero comparando los valores en los puntos del eje central del núcleo, vemos que las presiones de FLAC son menores que las de PLAXIS. En cambio, en la roca las presiones aumentan en FLAC y superan a las de PLAXIS.
Tabla 16. Resumen de las presiones de agua, tensiones totales verticales y tensiones
Efectivas medidas en el eje central de la presa.
Pw(KN/m2)
σyy(KN/m2)
σ'yy(KN/m2)
H(m)
FLAC PLAXIS FLAC PLAXIS FLAC PLAXIS 28 0.00 0.00 -55.80 -65.18 -55.80 -65.18 24 12.60 16.48 -122.10 -125.83 -109.50 -109.35 20 36.20 43.08 -186.00 -177.93 -149.80 -134.85 16 57.40 65.56 -247.40 -241.69 -190.00 -176.13 12 78.30 87.69 -309.80 -299.07 -231.50 -211.38 8 101.60 112.96 -375.90 -357.06 -274.30 -244.10 4 129.50 143.92 -455.10 -429.04 -325.60 -285.12 0 162.00 179.15 -542.70 -498.71 -380.70 -319.56 -3 217.00 208.90 -628.60 -533.87 -411.60 -324.97 -6 241.20 236.99 -687.30 -566.92 -446.10 -329.93
77
Tensiones
La introducción del flujo de agua en la presa provoca una modificación de las
tensiones totales con respecto al caso B2 (ver figura 36). En la tabla 16 se
resumen los resultados de las tensiones verticales tomadas en el eje central del
núcleo (x = 64 m), así como de las tensiones efectivas. Las curvas que
relacionan tensiones con altura de la presa se presentan en la figura 35.
En cuanto a los valores máximos y mínimos obtenidos se resumen a
continuación, en la tabla 17.
Tabla 17. Valores máximos y mínimos de la tensión total vertical. Las presiones de agua se
han tomado en estos puntos obteniéndose así la tensión efectiva.
MAXσyy (KN/m2)
pw (KN/m2)
σ'yy (KN/m2)
FLAC PLAXIS FLAC PLAXIS FLAC PLAXIS
(x,y) (56.12;0.00) (57.00;0.00) (56.12;0.00) (57.00;0.00) (56.12;0.00) (57.00;0.00)
tensión -740.70 -723.03 263.30 276.60 -477.40 -446.43
MINσyy (KN/m2)
pw (KN/m2)
σ'yy (KN/m2)
FLAC PLAXIS FLAC PLAXIS FLAC PLAXIS
(x,y) (71.00;0.00) (72.00;0.00) (70.00;0.00) (72.00;0.00) (70.00;0.00) (72.00;0.00)
tensión -229.30 -115.74 70.31 63.53 -158.99 -52.21
78
Tensiones totalesFLACvsPLAXIS
-6-4-202468
101214161820222426283032
-800-700-600-500-400-300-200-1000
Tensión (KN/m2)
Altu
ra p
resa
H (m
)
syy PLAXISsxx PLAXISsyy FLACsxx FLAC
Tensiones efectivasFLACvsPLAXIS
-6-4-202468101214161820222426283032
-600 -500 -400 -300 -200 -100 0
Presión agua (KN/m2)
Altu
ra p
resa
H (m
)
pw PLAXISpw FLACs'yy PLAXISs'yy FLAC
Figura 35. Curvas de tensión vertical y horizontal total, curvas de presión de agua y curvas de tensión vertical específica tomadas en el núcleo central de la presa.
79
Figura 36a. Distribución de las tensiones totales verticales en FLAC.
B C D E
FF F
G H I J K L M N
10 .000 20.000 30.000 40.000 50 .000 60 .000 70 .000 80.000 90.000 100 .000
-20 .000
-10 .000
0 .000
10.000
20.000
30.000
40.000
50.000
M ean s tressesExtrem e m ean stress -581 .07 kN /m2
[ kN /m2]
A : -600 .00
B : -560 .00
C : -520 .00
D : -480 .00
E : -440 .00
F : -400 .00
G : -360 .00
H : -320 .00
I : -280 .00
J : -240 .00
K : -200 .00
L : -160 .00
M : -120 .00
N : -80.000
O : -40.000
P : 0 .000
Q : 40 .000
Figura 36b. Distribución de las tensiones principales totales en PALXIS,
Figura 36. Distribución de las tensiones totales en FLAC y en PLAXIS.
80
Desplazamientos En la figura 37 podemos ver la distribución de los desplazamientos en la presa.
Como se puede ver, el efecto del agua produce una descarga en el espaldón
de aguas arriba de la presa debido a las presiones de agua, provocando
desplazamientos verticales hacia arriba, mientras que en el resto de la presa
los movimientos son horizontales y verticales hacia abajo. En las figuras 38 y
39 aparecen los movimientos verticales y horizontales respectivamente, donde
se observa mejor las diferencias entre los resultados de ambos programas.
El valor del desplazamiento máximo es de 13.68 cm en FLAC y se produce en
el punto x = 55 m , y = 22 m, es decir, en el espaldón de aguas arriba de la
presa, cerca del contacto con el núcleo central. En PLAXIS el valor del
desplazamiento máximo es de 10.30 cm y se produce en el núcleo de la presa,
a una altura de la base de y = 9 m, y en x = 64 m.
Figura 37a. Vectores de desplazamiento. El valor máximo es de 13.68 cm y se da en el
punto (54;23).
81
35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000 95.000
-10.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
To ta l displacementsE xtrem e tota l d isplacem ent 102.99*10-3 m
Figura 37b. Desplazamientos en PLAXIS
Figura 37. Desplazamientos provocados por el flujo del agua. Hay diferencias entre los
resultados obtenidos por FLAC que por PLAXIS. FLAC da movimientos máximos del orden de 13 cm, situados en el talud de aguas arriba de la presa mientras
que PLAXIS da un valor máximo de 10.99 cm, y corresponde a un punto situado en el núcleo (x = 64 m; y = 9 m).
Figura 38a. Distribución de los desplazamientos verticales en FLAC
82
BC
DEFGH
IJ
K
35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
Vertical displacem entsExtreme vertical displacement 56.07*10-3 m
[*10-3 m]
A : -50.000
B : -40.000
C : -30.000
D : -20.000
E : -10.000
F : 0.000
G : 10.000
H : 20.000
I : 30.000
J : 40.000
K : 50.000
L : 60.000
Figura 38b. Distribución de los desplazamientos verticales en PLAXIS
Figura 38. Desplazamientos verticales en PLAXISy en FLAC. El valor máximo en PLAXIS es
de 5.61 cm y se da en el punto (50.0;18.8), hacia arriba. El vector de desplazamiento vertical máximo hacia abajo se produce en el punto (56.35;13.20),
y su valor es de 4.40 cm. El valor máximo en FLAC es de 13.45 cm, y se da en el punto (53.38;23.64).
5.3.3.- INTERPRETACIÓN DE LOS RESULTADOS
En este nuevo caso se ha resuelto el problema de flujo acoplado en los dos
programas, donde se analiza los resultados correspondientes a la red de flujo,
a las presiones de agua, a las tensiones totales y efectivas, y a los
desplazamientos.
En cuanto a la red de flujo, tanto FLAC como PLAXIS dan resultados muy
parecidos. Al final del cálculo, ambos programas dan el vector de caudal
máximo por unidad de área, en el que FLAC da un valor inferior a PLAXIS
(qFLAC = 3.23E-06 m/s , qPLAXIS = 5.55E-06 m/s). Esta diferencia es muy
pequeña y no se tiene en cuenta ya que los dos valores están dentro del mismo
orden de magnitud. Tanto las líneas de corriente como las equipotenciales
83
Figura 39a. Desplazamientos horizontales en FLAC
B
CD
EF
G
H
IJ
K
L
35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000 80.000 85.000 90.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
35.000
Horizontal displacementsExtreme horizontal displacement 94.13*10-3 m
[*10-3 m]
A : -20.000
B : -10.000
C : 0.000
D : 10.000
E : 20.000
F : 30.000
G : 40.000
H : 50.000
I : 60.000
J : 70.000
K : 80.000
L : 90.000
M : 100.00
Figura 39b. Desplazamientos horizontales en PALXIS
Figura 39. Desplazamientos horizontales. En PLAXIS el valor máximo es de 9.41 cm, y se
da en el punto (63.50;8.00). El valor máximo es de 11.43 cm, y se da en el punto (63.13;7.86).
84
presentan una distribución prácticamente idéntica. Hay que señalar que en
FLAC se producen pequeños errores en los vectores de flujo de agua en
algunos nodos, tal y como se puede ver en la figura 32a), que son comparables
a los caudales reales, ya que las permeabilidades de trabajo son de por sí
valores muy pequeños.
En las curvas de presión de agua-altura de la presa en el eje del núcleo central
se observa como de nuevo PLAXIS da valores superiores de presión en el
núcleo de la presa, mientras que al llegar a la roca, las presiones de agua
aumentan en FLAC y superan las del PLAXIS. Este salto en los valores de las
presiones de agua también se produce para el caso de las tensiones totales
calculadas en los casos B1 y B2 y puede ser debido a la dificultad que
encuentra el programa FLAC en el contacto entre dos materiales con
diferencias grandes entre sus módulos elásticos.
También se han analizado las tensiones totales y efectivas en el núcleo. En
este caso la tendencia de las tensiones totales cambia con respecto a los
casos anteriores ya que las tensiones verticales obtenidas por FLAC superan a
las del PLAXIS a medida que nos acercamos a la base de la presa.
En cambio, las tensiones horizontales no varían de un programa a otro. Los
resultados obtenidos dibujan una curva de tensión horizontal-altura de la presa
con numerosos saltos cuantitativos que se pueden deber a errores de tipo
numérico y se dan en los dos programas, aunque con más gravedad en
PLAXIS.
Estas variaciones entre las presiones de agua y las tensiones totales hacen
que aumente la diferencia entre las tensiones efectivas en el núcleo de la
presa, llegando a ser la diferencia de un 16% entre las tensiones verticales
efectivas superior en FLAC que en PLAXIS (380.70 KN/m2 en FLAC frente a
los 319.56 KN/m2 en PLAXIS).
En general, la distribución de las tensiones en todo el cuerpo de la presa es
muy parecida en los dos programas, produciéndose una variación del efecto
85
arco producida por el flujo de agua. En el contacto entre espaldón de aguas
arriba con el núcleo aumentan las tensiones con respecto al caso B2, mientras
que en el contacto entre núcleo y espaldón de aguas debajo de la presa, se
mantiene el efecto arco, e incluso su efecto aumenta. Las tensiones mínimas
se producen en el núcleo, pero no en su eje (como ocurría en el caso de
construcción de la presa), sino que éstas se desplazan hacia el espaldón de
aguas abajo debido fundamentalmente a la carga de agua introducida sobre el
espaldón de aguas arriba de la presa. Las diferencias más notables entre FLAC
y PLAXIS se dan en la tensión vertical total mínima (en FLAC es un 50% mayor
que en PLAXIS), y por tanto, en la tensión vertical efectiva mínima (en FLAC es
un 67% superior a la de PLAXIS).
Curiosamente, no ocurre lo mismo con la tensión vertical máxima, donde la
diferencia entre los resultados de ambos programas es de tan solo un 2.5% a
favor de FLAC.
Estas diferencias en los valores de las tensiones efectivas se traduce también
en la aparición de más puntos en estado plástico en PLAXIS que en FLAC (ver
figura 40).
Si comparamos estos resultados con los del caso anterior, vemos como el
efecto del embalse de la presa produce una �isminución en las tensiones
totales máximas y en las tensiones totales mínimas (ver tabla _ donde se
comparan los resultados de los casos B2 y C1).
En cuanto a los desplazamientos producidos por el embalse de la presa, en
líneas generales el movimiento global que se da en la presa es igual en los dos
programas. El valor del desplazamiento máximo es superior en FLAC en un
2%, aunque el punto donde se da este valor varía entre un programa y otro.
FLAC obtiene el desplazamiento máximo en el espaldón de aguas arriba, a una
altura de 22 m, mientras que en PLAXIS el valor máximo se da en el núcleo, a
una altura media de la presa. Estas diferencias se aprecian mejor en la
86
Figura 40a. Puntos en estado plástico en FLAC 20.000 30.000 40.000 50.000 60.000 70.000 80.000 90.000 100.000 110.000 120.000 130.0
-20.000
-10.000
0.000
10.000
20.000
30.000
40.000
50.000
Plastic Points
Plastic Mohr-Coulomb point Tension cut-off point Figura 40b. Puntos en estado plástico en PLAXIS
Figura 40. Puntos en estado plástico en FLAC y en PLAXIS. Vemos como en PLAXIS
plastifican los puntos del espaldón donde las tensiones verticales disminuyen. En FLAC estos puntos no plastifican.
87
comparación entre los desplazamientos verticales y horizontales, donde FLAC
obtiene valores del orden de los 10 cm, mientras que PLAXIS son del orden de
los 5 cm. Esta diferencia en los desplazamientos se explica porque se obtienen
a partir de las tensiones, y en este caso FLAC obtiene tensiones efectivas más
elevadas que PLAXIS.
CASO B2
CASO C1
FLAC
PLAXIS
FLAC
PLAXIS
X
(m) Y
(m) σyy
(KN/m2) X
(m) Y
(m)σyy
(KN/m2) X
(m)Y
(m)σyy
(KN/m2) X
(m) Y
(m) σyy
(KN/m2)
NÚCLEO 64 0 364 64 0 377 64 0 543 64 0 499
X
(m) Y
(m) MIN σyy (KN/m2)
X (m)
Y (m)
MIN σyy(KN/m2)
X (m)
Y (m)
MIN σyy (KN/m2)
X (m)
Y (m)
MIN σyy (KN/m2)
71 0 169 71.6 0.33 149 71 0 229 72 0 115
ESPALDÓN X
(m) Y
(m) MAX σyy (KN/m2)
X (m)
Y (m)
MAX σyy(KN/m2)
X (m)
Y (m)
MAX σyy (KN/m2)
X (m)
Y (m)
MAX σyy (KN/m2)
72 0 755 57 0.5 815 57 0 741 57 0 723
Tabla 18. Resumen y comparación de los resultados de las tensiones totales máximas,
mínimas y la que se producen en el eje de la presa para H = 0 m, en las casos B2 y C1.
5.4.- CASO C2: CÁLCULO DE FLUJO. DESEMBALSE DE LA PRESA 5.4.1.- CÁLCULO REAL DEL PROBLEMA DE FLUJO ACOPLADO. MÉTODO DE MOHR-COULOMB En este caso se realiza el cálculo de flujo acoplado añadiendo una nueva etapa
de cálculo a las 9 anteriores (8 para construcción de la presa y 1 para el cálculo
del embalse de la presa), y que corresponde al caso de desembalse rápido de
la presa.
88
Las condiciones de contorno aplicadas corresponden a una variación en el
nivel de altura del agua en el espaldón de aguas arriba de la presa, pasando de
una altura de H = 28 m, a una de H = 5 m, y su correspondiente carga de agua
(ver figuras 41 y 42).
Figura 41. Condiciones de contorno definidas en PLAXIS.
Figura 42. Condiciones de contorno definidas en FLAC.
89
En PLAXIS el usuario tiene la opción de elegir el tipo de problema a resolver. El
caso de desembalse rápido de la presa correspondería a un caso no drenado
ya que las presiones de agua tardan en disiparse. En cambio, con FLAC no
puedes elegir entre el caso drenado y no drenado debido a que el propio
programa, al utilizar el tiempo como variable de cálculo, es el usuario el que ha
de interpretar si el problema al final es drenado o no drenado.
8.3.2.- RESUMEN DE LOS RESULTADOS OBTENIDOS
Al igual que en el caso anterior, analizaremos los resultados de la superficie
piezométrica dentro de la presa, la distribución de presiones de agua y el efecto
sobre las tensiones totales y los desplazamientos que se generan por el flujo
de agua en la presa.
Superficie piezométrica
En la figura 43 se presenta la red de flujo que se obtiene con FLAC y las líneas
equipotenciales obtenidas con PLAXIS.
Figura 43a. Red de flujo obtenida por FLAC.
90
IJ
KK
K
LL
L
MMM
NOPQRS
10.000 15.000 20.000 25.000 30.000 35.000 40.000 45.000 50.000 55.000 60.000 65.000
-10.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
35.000
Active groundwater headExtreme groundwater head 43.47 m
[ m]
A : -2.500
B : 0.000
C : 2.500
D : 5.000
E : 7.500
F : 10.000
G : 12.500
H : 15.000
I : 17.500
J : 20.000
K : 22.500
L : 25.000
M : 27.500
N : 30.000
O : 32.500
P : 35.000
Q : 37.500
R : 40.000
S : 42.500
T : 45.000
Figura 43b. Líneas equipotenciales obtenidas por PLAXIS.
Figura 43. Red de flujo calculada por FLAC (fig. 43a), y por PLAXIS (fig. 43-b).
Presiones de agua
La distribución de las presiones de agua se muestran en la figura 44. En la
tabla 19 se resumen los valores obtenidos de la presión de agua en el núcleo
de la presa. En este caso las diferencias entre los resultados obtenidos por los
dos programas aumentan, siendo PLAXIS el que da los valores más altos de
presión de agua.
91
pw(KN/m2)
H(m)
FLAC PLAXIS 28 0.00 -13.98 24 5.70 2.43 20 29.80 42.42 16 51.30 73.18 12 73.10 102.68 8 98.00 138.51 4 126.00 170.74 0 148.10 208.31 -3 204.50 221.62 -6 239.60 253.23
Tabla 19. Presiones de agua en el núcleo.
En la figura 44 se observa la distribución de las presiones de agua obtenida por
los programas.
Figura 44a. Presiones de agua en FLAC.
92
BC
C
C
CC
CC
DD
D
DD
E
EE
FF
FG
5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000
-10.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
35.000
40.000
Active pore pressuresExtreme active pore pressure -297.23 kN/m2
(pressure = negative)
[ kN/m2]
A : -300.00
B : -250.00
C : -200.00
D : -150.00
E : -100.00
F : -50.000
G : 0.000
Figura 44b. Presiones de agua en PLAXIS.
Figura 44. Presiones de agua obtenidas por FLAC y por PLAXIS. En general, las presiones de agua obtenidas por PLAXIS son mayores que las obtenidas por FLAC. La diferencia puede hallarse en que PLAXIS se ha realizado un cálculo no drenado, mientras que FLAC ha invertido más tiempo para la resolución del problema, con lo que el flujo está más refinado. Tensiones
En este caso vuelve a producirse una redistribución de las tensiones debido a
la descarga del agua en el espaldón de aguas arriba de la presa. Se han
analizado las tensiones verticales en el eje central de la presa (x = 64 m), con
tal de poder comparar los dos programas así como la variación entre este caso
y el anterior. Como se puede ver en la figura 45, las curvas de tensión vertical y
horizontal total, son muy parecidas en este caso en ambos programas. La
mayor diferencia está en el resultado de las tensiones efectivas, ya que éstas
dependen directamente del valor de las presiones de agua (ver tabla 20).
93
σyy(KN/m2)
σ'yy(KN/m2)
H(m)
FLAC PLAXIS FLAC PLAXIS28 -67.50 -85.58 -67.50 -99.5624 -156.90 -156.64 -151.20 -154.2020 -231.40 -221.79 -201.60 -179.3716 -285.60 -277.60 -234.30 -204.4212 -336.80 -329.73 -263.70 -227.058 -390.20 -390.79 -292.20 -252.284 -458.60 -459.31 -332.60 -288.580 -539.00 -529.02 -390.90 -320.70-3 -576.10 -568.29 -371.60 -346.67-6 -628.90 -583.45 -389.30 -330.22
Tabla 20. Tensiones verticales totales y efectivas.
Tensiones totalesFLACvsPLAXIS
-6-4-202468
101214161820222426283032
-800-700-600-500-400-300-200-1000100
Tensión (KN/m2)
Altu
ra p
resa
H (m
)
syy PLAXISsxx PLAXISsyy FLACsxx FLAC
Figura 45a. Tensiones verticales y horizontales totales medidas en el núcleo.
94
Tensiones efectivasFLACvsPLAXIS
-6-4-202468101214161820222426283032
-600 -500 -400 -300 -200 -100 0
Presión agua (KN/m2)
Altu
ra p
resa
H (m
)pw PLAXISpw FLACs'yy PLAXISs'yy FLAC
Figura 45b. Presiones de agua y tensiones efectivas verticales medidas en el núcleo.
Figura 45. Comparación entre FLAC y PLAXIS de las tensiones medidas en el eje
central del núcleo de la presa. En general, la distribución de las tensiones totales en el cuerpo de la presa no
presentan grandes variaciones entre un programa y otro (ver figura 45). No
obstante, sí que existen variaciones en los valores de tensión máxima y
mínima, cuyas diferencias se mantienen con respecto a las ya obtenidas en el
caso C1 (tabla 20).
Tabla 20. Tensiones máximas y mínimas obtenidas en el cálculo por FLAC y PLAXIS.
MAX σyy (KN/m2)
pw (KN/m2)
σ'yy
(KN/m2)
FLAC PLAXIS FLAC PLAXIS FLAC PLAXIS (x,y) (83.15;0.00) (57.00;0.00) (83.15;0.00) (57.00;0.00) (83.15;0.00) (57.00;0.00)
tensión -711.70 -647.51 73.23 240.11 -638.47 -407.40
MIN σyy (KN/m2)
pw (KN/m2)
σ'yy
(KN/m2)
FLAC PLAXIS FLAC PLAXIS FLAC PLAXIS (x,y) (71.00;0.00) (72.00;0.00) (71.00;0.00) (72.00;0.00) (71.00;0.00) (72.00;0.00)
tensión -244.90 -137.84 80.42 80.65 -164.48 -57.19
95
Desplazamientos
Los desplazamientos que se producen en este caso no son comparables en
ambos programas en términos cuantitativos, aunque si lo son en términos
cualitativos ya que en ambos casos se produce una rotura del talud de aguas
arriba de la presa.
En esta fase de cálculo PLAXIS deja de iterar cuando llega la estructura a
colapsar, obteniéndose un desplazamiento máximo de 0.50 m, mientras que
FLAC sigue iterando a pesar de que el sistema no sea estable, con lo que los
desplazamientos máximos que se producen son del orden de varios metros
(ver figura 46), aunque éstos corresponden a puntos del contorno del talud de
aguas arriba y los valores no son continuos sino que son muy dispersos.
Figura 46a. Vectores de desplazamiento obtenidos por FLAC. Después de casi 500.000
iteraciones, el programa no acaba de llegar a converger y continua deformándose, obteniendo movimientos muy elevados en el contorno.
96
0.000 10.000 20.000 30.000 40.000 50.000 60.000 70.000
-10.000
0.000
10.000
20.000
30.000
Total displacementsE xtreme total displacement 498.26*10-3 m
Figura 46b. Desplazamientos obtenidos por PLAXIS. El programa deja de iterar cuando la
estructura colapsa, obteniendo un valor máximo del desplazamiento de 0.50 m.
Figura 46. Desplazamientos calculados por FLAC y por PLAXIS.
5.4.3.- INTERPRETACIÓN DE LOS RESULTADOS
En este último caso analizado aparecen más diferencias entre los dos
programas en los resultados de la red de flujo. PLAXIS obtiene una red de flujo
más heterogénea que FLAC, debido a las condiciones no drenadas impuestas
y a las pocas iteraciones empleadas para el cálculo ya que éste queda
interrumpido cuando la estructura llega al colapso. En cambio, en FLAC se
invierten más iteraciones para el cálculo del flujo acoplado, con lo que al final
obtiene una red más homogénea, aunque aparecen zonas en el espaldón de
aguas arriba donde las líneas de flujo y las equipotenciales quedan cortadas.
Lo mismo ocurre con las presiones de agua, donde FLAC da valores inferiores
a PLAXIS tanto en el núcleo como en el espaldón de aguas arriba. En las
97
curvas de presión de agua-altura de la presa tomadas en el eje del núcleo, se
observa un aumento en la diferencia entre FLAC y PLAXIS con respecto al
caso anterior. En este caso, la diferencia entre la presión de agua medida en la
base de la presa (H = 0 m), es de un 29% superior en PLAXIS, mientras que en
caso C1 la diferencia era de tan sólo un 10%. Estas diferencias aumentan en la
zona del espaldón de aguas arriba, donde la permeabilidad del material es
mayor que en la del núcleo y, por tanto, FLAC disipa las presiones de agua de
forma más rápida que en PLAXIS, donde se le ha impuesto la condición de
material no drenado.
Estas diferencias en la presión del agua no se aprecian en las tensiones
totales, donde la diferencia entre los valores de ambos programas medidos en
la base del núcleo (H = 0 m), es de un 2%. En cambio, sí se siguen
manteniendo las diferencias entre las tensiones máximas y mínimas, donde
FLAC, tal y como ocurrió en el caso C1, da valores superiores a PLAXIS.
De todas formas, aunque la tendencia es que FLAC de valores de tensión
efectiva superiores a PLAXIS, tanto en uno como en otro programa el talud de
aguas arriba de la presa entra en rotura debido a la presencia de un alto
número de puntos en estado plástico, ya que disminuyen las tensiones por el
desembalse del agua pero las presiones de agua se mantienen. El número de
puntos en estado plástico sigue siendo superior en PLAXIS que en FLAC (ver
figura 47).
Los desplazamientos en este caso no se pueden valorar cuantitativamente, ya
que PLAXIS al dejar de iterar da un valor de desplazamiento máximo
correspondiente al último cálculo realizado, mientras que FLAC al continuar
iterando obtiene valores muy alejados de la realidad. Cualitativamente, los
resultados obtenidos son similares en ambos programas, donde los
desplazamientos máximos se producen en el espaldón de aguas arriba de la
presa, produciendo una superficie de rotura circular en el talud.
98
Figura 47a. Puntos en estado plástico en FLAC.
5.000 10.000 15.000 20.000 25.000 30.000 35.000 40.000 45.000 50.000 55.000 60.000 65.000 70.000 75.000
-10.000
-5.000
0.000
5.000
10.000
15.000
20.000
25.000
30.000
35.000
Plastic Points
Plastic Mohr-Coulomb point Tens ion cut-off point
Figura 47b. Puntos en estado plástico en PLAXIS.
Figura 47. En PLAXIS prácticamente todo el espaldón de aguas arriba de la presa está dentro del campo de la plasticidad.