TRABAJO FINAL ENERO 2019
Informe
Ejemplo de aplicación Poisson
2D – Flujo bajo una presa
impermeable
Universidad Técnica Federico Santa María
15 de enero de 2019
Nombres Apellidos Rol
Cristóbal Hernán Acevedo Farías 201211001-7
MÉTODOS DE ELEMENTOS FINITOS
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 2
MÉTODOS DE ELEMENTOS FINITOS
I. RESUMEN
En esta memoria de cálculo, se aplica un algoritmo de resolución el cual permite resolver,
en un dominio definido, la ecuación de Poisson haciendo uso del método de elementos
finitos. Específicamente, se analiza el flujo bajo una presa impermeable
Para ello, se define previamente la malla en el dominio, las condiciones de borde,
propiedades de los suelos bajo la presa y orden de la cuadratura Gaussiana a utilizar.
Específicamente, se calcula el caudal que escurre bajo la presa y se dibujan las líneas
equipotenciales correspondientes. Se consideran suelos sin y con estratificación y se
estudia cómo afectan a los resultados la anisotropía en las capas de suelo.
Se concluye de este trabajo que el orden en que se presentan las capas de suelo influye
en el caudal de salida; a pesar de que una estratificación esté conformada por los mismos
suelos, la configuración que presente en la superficie el suelo más impermeable,
presentará un menor caudal de salida. De ahí la importancia de colocar un relleno
impermeable en la superficie si se quiere retener el fluido, como en un embalse para
riego, por ejemplo.
Por otra parte, la anisotropía afecta la red de flujo; cuando el flujo va en la dirección que
posee la menor permeabilidad, se genera una caída de carga muy rápida. Cuando la
permeabilidad es menor en la dirección vertical respecto de la horizontal, el flujo posee un
cambio de dirección más fuerte y, por tanto, pasa más cerca de la presa.
MÉTODOS DE ELEMENTOS FINITOS ENERO
2019
MEMORIA DE CÁLCULO 1
MÉTODOS DE ELEMENTOS FINITOS
ÍNDICE DE CONTENIDO
I. Resumen...................................................................................................................1
1. Introducción...............................................................................................................4
2. Marco teórico.............................................................................................................5
3. problema....................................................................................................................7
3.1 enunciado..................................................................................................................7
3.2 solución.....................................................................................................................9
4. Conclusiones...........................................................................................................10
ÍNDICE DE FIGURAS
Figura 1 Elemento cilíndrico - Problema...............................................................................7
Figura 2 Condición de borde 1..............................................................................................8
Figura 3 Condición de borde 2..............................................................................................8
Figura 4 Distribución de temperatura en una barra a lo largo del tiempo.............................9
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 2
MÉTODOS DE ELEMENTOS FINITOS
ÍNDICE DE TABLAS
NO SE ENCUENTRAN ELEMENTOS DE TABLA DE ILUSTRACIONES.
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 3
MÉTODOS DE ELEMENTOS FINITOS
1. INTRODUCCIÓN
En nuestro país es extendido el uso de presas, ya sea para riego, como tranques de
relave o para la generación de energía hidroeléctrica, por ejemplo. El uso de presas
impermeables origina filtraciones bajo ésta, las cuales puede ocasionar problemas de
estabilidad, pérdida de agua o una descarga de sustancias contaminantes al ambiente.
Por ello, es de interés evaluar estos caudales asociados, en distintos escenarios respecto
a las propiedades hidrogeológicas de los suelos sobre los que se fundan estas presas.
Existen modelos donde se puede evaluar el caudal de infiltración mediante una fórmula
analítica, los cuales consideran condiciones ideales como medios homogéneos,
isotrópicos y régimen permanente. Para cambiar alguna de estas condiciones es
necesario recurrir a programas como MODFLOW o SEEP/W. La similitud entre ambos es
la utilización de modelos de elementos finitos, para así ampliar el espacio de problemas
que pueden ser resueltos.
En este trabajo, es usado un programa llamado POISSON el cual permite resolver
problemas que involucren ecuaciones tipo Poisson. Este programa permite la definición
de más de un estrato de suelo y, además, considerar anisotropía
En la sección 2 se da un marco teórico respecto al problema que se va a tratar. En la
sección 3 se plantea el problema específico que se pretende resolver. En la sección 4 se
desarrolla el problema, se obtienen resultados y se dan comentarios respecto al mismo.
Se finaliza en la sección 5 con las conclusiones del trabajo.
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 4
MÉTODOS DE ELEMENTOS FINITOS
2. MARCO TEÓRICO
La ecuación de flujo en dos dimensiones está basada en la Ley de Conservación de la
Masa y en la Ley de Darcy.
El primer principio establece que el flujo que entra a una unidad de volumen debe ser
igual al flujo que sale. En esta afirmación existen algunas suposiciones:
El flujo se da en dos dimensiones (x-y).
El fluido es incompresible.
No existe acumulación dentro de los volúmenes.
Se desprecian las pérdidas.
El esquema descrito se muestra en la Figura 1.
Figura 1 Elemento de volumen unitario con caudales
De acuerdo a la Figura 1, se escribe la ecuación de continuidad:
∂vx ∂vy
v x dydz+ v y dxdz=(v ¿¿ x + dx)dydz +(v ¿¿ y + dy)dxdz ¿ ¿ (1)
∂x ∂y
El término v i es la componente de la velocidad del flujo en la dirección i . Se debe notar
que la ecuación (1) es un balance de caudales. El caudal Q se define como:
Q=v∗A (2)
, donde v es la velocidad del flujo y A es el área cuya componente normal coincide con la
dirección del flujo.
Eliminando términos de la ecuación (1), se obtiene:
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 5
MÉTODOS DE ELEMENTOS FINITOS
∂vx ∂vy
+ =0 (2)
∂x ∂ y
El segundo principio se refiere a la Ley de Darcy, la cual se expresa matemáticamente en
la siguiente ecuación:
dH
v i=k i , i=x , y (3)
di
En la ecuación precedente, k i se refiere a la permeabilidad del suelo en la dirección i, y H
es la carga hidráulica. Esta se compone de:
2
u v
H= + z+ (4)
γ 2g
2
u v
, donde es la carga de presión de poros, z es la carga de elevación y es la carga de
γ 2g
velocidad. Éste último término es muy bajo en los suelos, dado que las velocidades de
infiltración son muy bajas.
Combinando las ecuaciones (2) y (3) se obtiene la ecuación de continuidad de Laplace.
∂
∂x(kx
dH
+
∂
dx ∂ y) (
ky
dH
dy
=0 ) (5)
A partir de la ecuación (5), se obtienen líneas de flujo y líneas equipotenciales. Las
primeras representan la trayectoria a lo largo de la cual una partícula de agua viaja. Por
otra parte, dentro de una línea equipotencial la carga hidráuica total (potencial) es
constante. Ambas líneas son ortogonales entre sí. La representación gráfica en conjunto
de ambos tipos de líneas es conocida como red de flujo.
3. PLANTEAMIENTO DEL PROBLEMA
En la Figura 2 se muestra un esquema típico de una presa impermeable.
Figura 2 Presa impermeable
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 6
MÉTODOS DE ELEMENTOS FINITOS
Es posible ver que existen dos niveles de agua, uno a cada lado de la presa de alturas h1
y h2 . Debido a la permeabilidad del suelo y la diferencia energética entre los niveles (
h1−h2), existe un flujo desde izquierda a derecha.
En el problema en particular, se va a considerar que h1=15[m] y h2 =5[m]. De acuerdo
con esto, se definen las siguientes condiciones de borde para el dominio llamado SUELO
en la Figura 3.
Figura 3 Condiciones de borde
La condición de altura H están dadas por las alturas de agua existentes a cada lado de la
presa. Existen flujos nulos en los contacto entre el suelo y la presa, dado que esta se
definió como impermeable. En la parte inferior del dominio también existe flujo 0 dado que
se considera la existencia de una superficie rocosa impermeable en ese nivel. Por último,
a los costados del dominio también se definió flujo nulo, dado que se considera sólo
aporte de flujo por parte de los niveles de agua definidos en la Figura 2. En otras
palabras, estos límites en los costados se deben tomar a una distancia lo suficientemente
grande de la presa de manera que los flujos que aportan las columnas de agua
consideradas sean verticales y lo suficientemente pequeña para que no considere aportes
de flujo de otras fuentes.
En la Figura 4 se ven las dimensiones del problema. Además, en la Figura 5 se ve una
estratificación con dos tipos de suelo con permeabilidades diferentes, caso que también
va a ser evaluado.
Por otro lado, es común que los suelos posean la permeabilidad como una propiedad
anisotrópica; debido a fenómenos de sedimentación y compactación, es común que la
permeabilidad vertical sea menor a la horizontal. Es más, la permeabilidad horizontal
puede ser más de 10 veces mayor según mediciones en terreno. Este efecto también se
considera en este estudio.
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 7
MÉTODOS DE ELEMENTOS FINITOS
Figura 4 Dimensiones y configuración del problema (medidas en [m])
Figura 5 Configuración de dos estratos (medidas en [m])
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 8
MÉTODOS DE ELEMENTOS FINITOS
Se consideran los siguientes casos:
Tabla 1 Casos de evaluación
Casos Capa Permeabilidad [m/día ]
1 Suelo k x =k y =2
2 Suelo k x =k y =1
Suelo 1 k x =k y =2
3
Suelo 2 k x =k y =1
Suelo 1 k x =k y =1
4
Suelo 2 k x =k y =2
5 Suelo k x =2∨¿ k y =0.2
6 Suelo k x =1∨¿ k y =0. 1
Suelo 1 k x =2∨¿ k y =0.2
7
Suelo 2 k x =1∨¿ k y =0.1
Suelo 1 k x =1∨¿ k y =0. 1
8
Suelo 2 k x =2∨¿ k y =0.2
Los casos del 1 al 4 contemplan suelos isotrópicos, mientras que del 5 al 8 son suelos
anisotrópicos. Los casos 1,2,5 y 6 presentan 1 estrato, mientras que los casos 3,4,7 y 8
presentan 2 estratos.
Las permeabilidades elegidas pueden corresponder a suelos tipo arena muy fina o
mezclas limo-arenosas con arcilla.
4. RESOLUCIÓN
Se aplica el algoritmo de resolución POISSON. De acuerdo a la (5), el potencial de la
ecuación de Poisson es la carga hidráulica total H .
Para aplicar el algoritmo, se debe definir una malla. Esta será definida con un paso
Δ x =Δ y=5 [m]. Se intentó definir con una malla el doble de fina, pero el programa no
corre por la gran cantidad de información.
5. PROBLEMA
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 9
MÉTODOS DE ELEMENTOS FINITOS
5.1 PROBLEMA 1
5.1.1 ENUNCIADO
Definir el sistema de ecuaciones de 6x6 (matriz de coeficientes y vector del lado derecho)
para el problema plano de Poisson utilizando la malla mostrada en la Figura 6.
Figura 6 Problema Poisson 2D - Problema 1
5.1.2 SOLUCIÓN
Se deben calcular los valores necesarios para obtener las expresiones de la ecuación
Error: Reference source not found. Estas hacen referencia a los elementos que forman
parte de la malla de la Figura 6. Como los elementos cuadrados poseen las mismas
dimensiones, basta con calcular las matrices necesarias una sola vez para luego asociar
los elementos de cada matriz a cada nodo de la malla. Lo mismo ocurre tanto para los
elementos triangulares del lado derecho como para los del lado izquierdo.
Se escoge como elemento base para los elementos cuadrados el mostrado en la Figura 7,
con el origen local de coordenadas mostrado.
De esta manera, se obtiene que:
L −L −L L
+y −y +y −y
2 2 2 2
2 2 2 2
L L L L
( )
L L −L −L
+x −x +x −x
2 2 2 2
2 2 2 2
L L L L
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 10
MÉTODOS DE ELEMENTOS FINITOS
Figura 7 Elemento cuadrado
5.2 PROBLEMA 2
5.2.1 ENUNCIADO
Derivar las expresiones para la matriz Jacobiana [J ] para las siguientes transformaciones:
Figura 8 Transformación N°1 - Problema 2
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 11
MÉTODOS DE ELEMENTOS FINITOS
Figura 9 Transformación N°2 - Problema 2
5.2.2 SOLUCIÓN
5.3 PROBLEMA 3
5.3.1 ENUNCIADO
Considerar la siguiente transformación de un elemento maestro de tres nodos.
Figura 10 Transformación 1D - Problema 3
5.3.2 SOLUCIÓN
6. CONCLUSIONES
¿?????????????????????
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 12
MÉTODOS DE ELEMENTOS FINITOS
Para que la transformación de elementos sea 1 a 1, el desplazamiento de los nodos
respecto al elemento maestro u original no debe ser exagerado.
MÉTODOS DE ELEMENTOS FINITOS ENERO 2019
MEMORIA DE CÁLCULO 13
