EL MÉTODO DE LOS ELEMENTOS FINITOS
METODOLOGIA GENERAL DE LOS ELEMENTOS FINITOS
El método de los elementos finitos, es un procedimiento basado en técnicas computacionales, que
puede ser usado para analizar estructuras y diferentes sistemas continuos. Es un método numérico
versátil, y que es ampliamente aplicado para resolver problemas que cubren casi todo el espectro de
análisis ingenieriles. Sus aplicaciones comunes, incluyen el comportamiento de sistemas estáticos,
dinámicos y térmicos. Los avances en el hardware, han facilitado y aumentado la eficiencia del
software de elementos finitos, para la solución de sistemas complejos de ingeniería sobre
computadores personales.
Los resultados obtenidos con el análisis de elementos finitos, son raramente exactos. Sin embargo,
una solución adecuada puede ser obtenida, si se usa un modelo apropiado de elementos finitos.
Si el objetivo del ingeniero, es el desarrollo de código de elementos finitos, entonces una profunda
comprensión de la teoría de elementos finitos es esencial. Si el objetivo, es el uso del código de
elementos finitos, entonces para el análisis es necesario tener:
1. Básica comprensión de los conceptos fundamentales del método de los elementos fi n i to s
2. Práctica incluyendo el conocimiento de las capacidades y limitaciones en el programa
computacional que va a ser usado
El documento presenta algunos conceptos básicos de la teoría de elementos finitos. Esta discusión
ayuda a entender cómo funciona el código de los elementos finitos y como se debe ser usado. Se
tratarán conceptos, referentes a la aplicación del método en sólidos mecánicos.
ORIGENES DEL METODO DEL ELEMENTO FINITO
El uso de métodos clásicos, es probablemente la mejor forma de analizar estructuras simples; no
obstante, su uso es poco aconsejable cuando el sistema es complejo. En estos casos la mejor
alternativa, es usualmente una solución obtenida con el método de los elementos finitos.
La primera diferencia entre los métodos clásicos y los elementos finitos son la forma de ver la
estructura y el consiguiente procedimiento de solución. Los métodos clásicos consideran la
estructura como continúo, cuyo comportamiento es gobernado por ecuaciones diferenciales
parciales ordinarias. El método de elementos finitos considera la estructura como el ensamble de un
número finito de partículas pequeñas. El comportamiento de las partículas, y de toda la estructura,
es obtenida por la formulación de un sistema algebraico de ecuaciones que puede ser
solucionado por medio de un computador.
Las partículas de tamaño finito, son llamadas elementos finitos. Los puntos donde los
elementos finitos son interconectados, son conocidos como nodos, y el procedimiento de selección
de nodos es llamado discretización o modelización,
Figura 1 Tanque cilíndrico modelado con elementos finitos
Normalmente, el análisis de elementos finitos involucra siete pasos. Pasos 1,2,4,5 y 7 requieren
decisiones realizadas por el usuario del programa de elementos finitos. El resto de los pasos son
realizados automáticamente por el programa de computadora.
PASOS PARA EL ANÁLISIS DE ELEMENTOS FINITOS:
1. Discretización o modelado de la estructura: La estructura es dividida en una cantidad finita de
elementos, con ayuda de un preprocesador. Este paso es uno de los más cruciales para obtener
una solución exacta del problema, de esta forma, determinar el tamaño o la cantidad de elementos
en cierta área o volumen del elemento a analizar representa una ventaja del método, pero a la vez
implica que el usuario debe estar muy conciente de esto para no generar cálculos innecesarios o
soluciones erróneas.
2. Definir las propiedades del elemento: En este paso el usuario debe definir las propiedades del
elemento.
3. Ensamblar las matrices de rigidez de los elementos: La matriz de rigidez de un elemento,
consiste de coeficientes los cuales pueden ser derivados del equilibrio, residuos ponderados o
métodos de energía. La matriz de rigidez del elemento se refiere a los desplazamientos nodales al
ser aplicadas fuerzas en los nodos (K*F = U). El ensamble de las matrices de rigidez, implica la
aplicación de equilibrio para toda la estructura.
4. Aplicación de las cargas: Fuerzas externas concentradas o fuerzas uniformes y momentos son
especificados en este paso.
5. Definir las condiciones de frontera: Las condiciones de apoyo deben ser dadas, por ejemplo, si
el desplazamiento de ciertos nodos es conocido. Usando los elementos de la frontera se pueden
determinar las reacciones en los mismos.
6. Solucionar el sistema de ecuaciones algebraicas lineales: La secuencial aplicación de los pasos
descritos, conduce a un sistema de ecuaciones algebraicas simultáneas, donde los
desplazamientos nodales son desconocidos.
7. Calcular los esfuerzos: El usuario puede entonces calcular los esfuerzos, reacciones,
deformaciones u otra información relevante. El post-procesador ayuda a visualizar la salida en forma
gráfica.
Los elementos de malla pueden ser elementos 1D, 2D o 3D en función del modelo a simular. El modelo FEA es
siempre un modelo tridimensional, por lo que si es 1D o 2D se utilizan elementos, es necesario definir 1 o 2
dimensiones en para ser analizado.
Estas dimensiones están definidas por las propiedades del elemento.
TIPOS DE MALLADO
El mallado es uno de los componentes clave para obtener resultados precisos de un modelo FEA. Los elementos
de la malla deben tener en cuenta muchos aspectos para poder discretizar los gradientes de tensión con precisión.
Normalmente, cuanto menor es el tamaño y mayor el número de elementos de la malla, más precisa es la solución,
ya que los diseños se muestrean mejor en los dominios físicos. La compensación es que cuanto mayor es la
precisión, más grandes se vuelven las simulaciones y, por lo tanto, los tiempos de resolución se extienden.
Estas piezas más pequeñas se denominan elementos (finitos). Los elementos están conectados entre sí en los
nodos. Estos elementos y nodos también se llaman mallas. El ensamblaje de elementos y nodos es llamado
modelo de elementos finitos.
PROPIEDADES DE LOS ELEMENTOS
Los elementos de malla pueden ser elementos 1D, 2D o 3D en función del modelo a simular. El modelo FEA es
siempre un modelo tridimensional, por lo que si es 1D o 2D se utilizan elementos, es necesario definir 1 o 2
dimensiones en para ser analizado.
Estas dimensiones están definidas por las propiedades del elemento
ELEMENTO 1D
Los elementos 1D se utilizan cuando la longitud de la estructura es 8 ~ 10 veces mayor que su ancho y espesor y
cuando las cargas externas se aplican solo a las juntas.
Hay 3 tipos de elementos 1D:
• Elemento de varilla: solo realizar compresión
• Elemento de celosía (barra): realiza tensión y compresión
• Elemento de viga: asumir tensión, compresión, cortantey momento
Las propiedades de sección y material deben asignarse para
elementos
Elementos 2D
Los elementos 2D se utilizan cuando la longitud y el ancho de la estructura es significativamente mayor que su
espesor.
Hay diferentes tipos de elementos 2D:
• Placa: elemento 2D general
• Membrana: no se dobla
• Superficie: sin espesor
• Deformación plana: la deformación normal a la superficie es 0
• Carcasa compuesta: elementos de placa 2D con capas compuestas
Elementos 3D
Es necesario asignar propiedades de espesor y material para elementos
Las propiedades de los elementos 3D solo están presentes para asignar materiales y sistemas de coordenadas de
materiales a mallas 3D.
Hay 2 tipos principales de elementos 3D, elementos tetraédricos (4 caras, 4 nodos) y elementos hexaédricos (6
caras, 8 nodos).
Los elementos de segundo orden, también llamados elementos cuadráticos, agregan 1 nodo medio a cada borde.
El cálculo es más preciso para elementos hexaédricos que para elementos tetraédricos. Y es más precisa para
elementos de segundo orden porque la función de aproximación será cuadrática en lugar de lineal.
Tipos de Elementos Finitos
Esta sección describe muchas características sobresalientes de los elementos más utilizados;
denominados, truss, beam, plane stress, plane strain, axisymmetric, membrane, plate, shell, solid ó
brick, tetrahedral, hexahedral, boundary, y gap . Los programas comerciales de elementos finitos
poseen una gran cantidad de elementos en sus librerías. Sin embargo, la mayoría de las
estructuras y aplicaciones mecánicas pueden ser solucionadas con los elementos básicos ya
mencionados.
Dependiendo la dimensión, los elementos básicos se pueden dividir en tres categorías: elemento
de línea, área y volumen. Truss, beam y los elementos de restricción, son de línea. Plane stress,
plain strain, axisymmetric, membrane, plate y shell son elementos de área. Solid ó brick, tetrahedral
y hexahedral son elementos de volumen. Los criterios para las selecciones del elemento apropido
para cada aplicación se verán más adelante.
Elementos ‘Truss’
El elemento truss, es un elemento caracterizado básicamente porque solo puede comportarse
como un miembro sometido a dos fuerzas (se sabe por tanto que estas cargas deben estar
dirigidas a lo largo del eje longitudinal del elemento).
Una estructura los elementos se pueden modelar como un elemento Truss si cumplen estos tres
requerimientos:
a. Su longitud es mucho mayor que su alto o ancho (entre 8 y 10 veces);
b. Esta es conectada con el resto de la estructura con pasadores que no transfieren momentos.;
y
c. Las cargas externas solo son aplicadas en el extremo de los elementos, y son paralelas al
mismo (Carga Axial).
Los elementos Truss solo pueden ser sometidos a tracción o compresión. De esta forma, la única
propiedad de la sección que se debe especificar es el área axial del elemento. La figura 2.2.1
muestra la geometria y las fuerzas nodales en un elemento truss tridimencional. Como se muestra
en la figura, un elemento truss tridimensional posee tres grados de libertad por nodo, esto es tres
desplazamientos sobre los ejes globales X, Y y Z.
Figura 2.2.1 Elemento Truss tridimencional
Elementos ‘Beam’
El elemento Beam, es probablemente el más usado. Además de sus aplicaciones obvias en estructuras,
muchos otros sistemas, como uniones mecánicas, sistemas de conductos, tuberías y vigas en puentes
pueden ser modeladas con el elemento ‘beam’.
Para miembros estructurales para ser modelados con elementos ‘Beam’, una de sus dimensiones debe
ser mucho mayor, por lo menos 10 veces más grande que las otras dos. Contrario al elemento truss, el
elemento beam puede estar sometido a cargas transversales y/o momentos flectores en adición a la
tracción y compresión. La geometría y los desplazamientos/rotación son mostrados en la figura 2.2.2.
Note que el elemento beam tridimensional posee seis grados de libertad por nodo, esto es, tres
desplazamientos y tres rotaciones sobre los ejes globales X, Y y Z.
Figura 2.2.2 Elemento beam tridimencional.
.
Los perfiles comunes de elementos beam, son la sección I, sección en T, caja, circular y
canales. Dentro de las propiedades de la sección, se debe especificar el área axial, la
resistencia a la torsión y el momento de inercia.
Elementos Elásticos bidimensionales
Hay tres tipos de elementos bidimensionales:
1. Plane Stress Elements (Esfuerzo plano)
2. Plane Strain Elements (Deformación plana)
3. Axisymmetric Elements (Elementos Axisimétricos).
Elementos sometidos a Esfuerzo Plano y Deformación Unitaria Plana
La explicación sobre la diferencia entre los casos de esfuerzo plano y deformación unitaria plana ya
fue definida en el primer capítulo del curso.
Para el caso de análisis plano existen principalmente dos tipos de elementos: Triangular y
Cuadrilátero. Dependiendo el tipo del tipo de esfuerzo al que está sometido el elemento, este se
debe modelar como esfuerzo plano o deformación unitaria plana.
Figura 2.2.3 Elemento bidimensional en esfuerzo plano: (a) Triangular; (b) Cuadrilátero.
Como regla, se prefieren los elementos cuadriláteros a los triangulares por razones de isotropía
geométrica. Sin embargo, se sugiere el uso de elementos triangulares cuando se presentan
irregularidades en la geometría del elemento a modelar, como se muestra en la figura 2.2.4.
Figura 2.2.4 Elementos cuadrilateros y triangulares.
Elementos Axisimetricos
Tanques de acero y concreto, rotores, conchas, toberas y contenedores son algunos ejemplos
representativos de estructuras axisimétricas. De forma similar a las estructuras tridimensionales
que están bajo condición de esfuerzo plano o deformación plana, las estructuras axisimétricas
sometidas a cargas axisimétricas, pueden ser analizadas en un modelo bidimensional. Para
analizar una estructura axisimétrica, como un cilindro de pared delgada t, sujeta a una presión
constante p, el modelo es la intersección del cilindro con el plano YZ como se muestra en la figura
2.2.5.
La carga p, es aplicada al modelo de elementos finitos como se muestra en la figura 2.2.5(b).
Cuadrilateros y triangulos axisimétricos poseen dos grados de libertad en cada nodo, figuras.
3.2.5(c) y (d).
Figura 2.2.5 (a) Cilindro bajo presión interna; (b) Modelo Axisimétrico; (c) Elemento Cuadrilatero; (d)
Elemento Triangulart.
Sólidos elásticos tridimensionales ó elementos ‘Brick’
Los elementos sólidos son elementos tridimencionales con tres grados de libertad translacional por
nodo, ver figura 2.2.6. Los nodos son usualmente introducidos en la intersección de los tres planos,
o la mitad de la intersección de dos planos Un elemento brick de 8 nodos, con sus respectivos
grados de libertad se puede apreciar en la figura 3 .2 .6 .
Figura 2.2.6 Elemento brick y sus grados de libertad por nodo (8 nodos).
El elemento brick puede proveer información acerca de la variación tridimensional de los
esfuerzos y deformaciones del elemento.
Elementos ‘Tetrahedral’ and ‘Hexahedral’
Así como los elementos brick, los elementos ‘tetrahedral’ y ‘hexahedral’ pueden ser usados para
modelar estructuras tridimensionales. El tetraedro puede ser visto como un triángulo en tercera
dimensión, como se ve en a figura 2.2.7, mientras que el hexaedro puede ser visto como un
cuadrilátero extendido en la tercera dimensión. Se puede apreciar entonces que el hexaedro tiene
la misma geometría del elemento brick de 8 nodos. La diferencia e n t r e
estos dos, es la formulación y precisión computacional. Por lo general los elementos tetraedro y el
hexaedro poseen solo tres grados de libertad por nodo, y la precisión de estos elementos se puede
incrementar colocando nodos en la mitad de los lados.
Figura 2.2.7 (a) Tetraedro de 4 nodos; (b) tetraedro de 10 nodos
Selección del tipo de Elementos
Antes de seleccionar el tipo de elemento para el modelo de una estructura, se debe primero
dibujar un bosquejo del sistema físico indicando su geometría, condiciones de frontera, cargas y
discontinuidades geométricas o de material.
El bosquejo debe además incluir un sistema de coordenadas globales y las dimensiones de la
estructura.
Después, se debe examinar si el modelo puede ser reducido, o simplificado. Esfuerzo plano,
deformación plana y modelos axisimétricos, permiten la reducción de problemas tridimensionales
a bidimensionales. Además, la presencia de planos de simetría permite modelar sólo una parte de
la estructura. El uso apropiado de la simetría en modelos estáticos ya fue discutido.
El bosquejo de un sistema físico, puede ayudar en la selección del elemento apropiado. Por
ejemplo, para modelar cargas transversales o axiales en elementos mecánicos, eléctricos y
estructuras civiles, se pueden usar elementos beam o truss. Elementos de esfuerzos planos son
apropiados para modelar en el plano de acción, placas y vigas cortas. Elementos de deformación
plana son usualmente utilizados para modelar paredes de contención y largos diques. Los
elementos axisimétricos son usados para modelar estructuras que son rotacionalmente simétricas
sobre uno de los ejes y cargado simétrica o antisimétricamente sobre el mismo eje, como los
cilindros sometidos a presión interna.
La selección del tipo de elemento o elementos también depende del tipo de resultados esperados.
Por ejemplo, el cilindro mostrado en la figura 2.3.1. con los extremos empotrados y sujeto a una
carga puntual, puede ser modelado de diferentes formas:
a. Un modelo de viga puede ser usado está interesado en una aproximación de la deflexión del
elemento.
b. Un modelo usando elementos shell/plate puede ser usado si el objetivo es obtener un cálculo
de los esfuerzos.
Figura 2.3.1 (a) Cilindro con extremos empotrados y carga punctual en la mitad; (b) Modelo de viga;
(c) Modelo con elementos Plate/shell.
Ciertos problemas pueden ser solucionados con más de un tipo de análisis. Por ejemplo, el cilindro
bajo presión interna mostrado en la figura 2.3.2 puede ser analizado de tres formas diferentes:
a. Un problema de deformación plana. En este caso la geometría del modelo consiste en una
sección circular obtenida por la intersección del cilindro con el plano YZ. Por el uso de elementos de
deformación plana, se pueden obtener resultados que son válidos para una sección del cilindro que
está lejos de los extremos empotrados.
b. Un problema axisimétrico. El modelo usando elementos axisimétricos, es mostrado en la figura
2.2.5. En este modelo todos los nodos sobre los lados AB y CD están restringidos, y la presión p es
aplicada en el lado AC. Las deformaciones y esfuerzos obtenidos del análisis del modelo
axisimétrico es válido para toda la estructura. Note que se puede analizar el cilindro con este
modelo, porque la carga también es axisimétrica.
Una estructura tridimensional. El cilindro es modelado usando elementos plate como se muestra en
la figura 2.3.1 (c). Este modelo es el más elaborado, pero menos eficiente, puede ser una opción, si
las cargas no fueran axisimétricas.
Figura 2.3.2 Cilindro bajo presión interna.
APLICACIONES PRACTICAS
1. Defina los nodos sobre o cerca de los puntos de aplicación de cargas, dixcontinuidades
geométricas, soportes y todas aquellas regiones que requieran información acerca de los esfuerzos
o desplazamientos.
Una malla de elementos finitos debe ser uniforme donde sea práctico. Sin embrago, la no
uniformidad es requerida para obtener resultados apropiados en cambios bruscos de la geometría y
carga. Solo estas partes, donde la geometría, las cargas o los esfuerzos, cambian dramáticamente,
se debe refinar la malla, ver figura 2.4.1
Figura 2.4.1 Placa con concentrador de esfuerzo.
2. Prefiera el uso de cuadriláteros, elmentos sólidos de 6 lados y hexagonos, excepto donde los
elementos triangulares y tetraedros son necesarios para acomodar irregularidades geométricas y
cargas.
3. Una malla más refinada es requerida para obtener con precisión esfuerzos y desplazamientos.
Si es necesario realizar un estudio de convergencia, por ejemplo, empiece con un número
relativamente pequeño de elementos y progresivamente genere modelos más refinados. La figure
2.4.2 muestra el procedimiento con dos modelos de esfuerzo plano en una viga corta en cantilever
modelada en el plano YZ.
Figure 3.4.2 Modelo de 4 y 6 elementos de una viga corta en cantilever.
Proporción
La proporción de los elementos es definida por la relación entre la dimensión más grande y más
pequeña del elemento. La figura 2.4.3 muestra la configuración de tres elementos bidimensionales
clasificados en ‘buenos’, ‘pobres’ e ‘ilegales’. Los buenos elementos se caracterizan por que su
proporción es cercana a la unidad y los ángulos están cerca a los 90 grados. Lo elementos pobres
deben ser evitados, pues estos pueden generar resultados inexactos. Los elementos ilegales son
inaceptables, y no deben ser usados para modelos de elementos finitos. Cuando elementos
ilegales son dibujados inadvertidamente, estos generan modelos de elementos finitos inválidos.
Muchos post-procesadores de programas de elementos finitos comerciales, permiten la
identificación de éstos errores.
Figura 2.4.3 Elementos (a)
Buenos; (b) Pobres; (c) legales.
4. Si el campo de esfuerzos tiene similar gradiente en todas las direcciones, trate de mantener
la proporción en 1. Elementos con alta proporción pueden ser inválidos, ver figura
2.4.4 (a). Sin embargo, los elementos catalogados como pobres, pueden ser utilizados en el análisis
de regiones de la estructura donde el gradiente de esfuerzos varía gradualmente sobre una gran
dimensión de los elementos. Los elementos ilegales, debe ser siempre disperse
Figure 3.4.4 (a) Proporción muy alta; (b) Proporción en un rango aceptable.
7. Para regiones con pequeñas variaciones de esfuerzos, el factor de proporción puede ser
alto, de 40 a 1, y aun así obtener buenos resultados. Sin embargo, como regla general, se
recomienda mantener el factor de proporción debajo de 10 para análisis de deformación, y bajo 3
para análisis de esfuerzo, ver figura 2.4.4 (b). Se puede notar que, en el análisis de deformación, el
énfasis se da en la precisión para calcular los desplazamientos de cada uno de los nodos, mientras
que, en el análisis de esfuerzos, se está interesado en la precisión del cálculo de ambos, esfuerzos
y desplazamientos.
8. El cambio repentino en el tamaño de los objetos debe ser evitado. Cuando el uso de
elementos de diferentes tamaños es necesario, la relación de dimensión de los objetos adyacentes
debe ser menor a dos.
CONCLUSIONES
- La implementación de la metodología para la construcción del modelo tridimensional del pozo
es fundamental para el análisis geo mecánico del mismo porque relaciona cada material que se está
perforando con su inmediatamente anterior y así establece un modelo heterogéneo a partir de
materiales homogéneos en sí mismos.
- Debido a la creación de bucles en el software, lo cual permite aprovechar las relaciones de
enumeración de las líneas, se puede construir el modelo tridimensional de forma automática, así
como también el enmallado, la aplicación de la presión de lodo y las restricciones.
- El método de los elementos finitos se presenta como otra alternativa de análisis de modelamiento
numérico para establecer el estado de esfuerzos y deformaciones de un pozo de petról
