3. Introduccin al mtodo del Elemento finito.

Dos de los objetivos ms importantes que se busca obtener mediante el anlisis

ingenieril es la capacidad de identificar los principios fsicos bsicos que rigen el
comportamiento de un sistema y adems transformar esos principios en modelos
matemticos compuestos de una o varias ecuaciones que puedan ser resueltas, con el fin
de predecir el comportamiento cuantitativo y cualitativo del sistema, teniendo presente
siempre que esta prediccin debe ser precisa.
El modelo matemtico resultante esta compuesto de una ecuacin o de un sistema de
ecuaciones, cuya solucin debe ser consistente y debe ser representativa de las bases y
principios fsicos del sistema.
En situaciones donde el sistema es relativamente simple, es posible analizar el problema
mediante el uso de sistemas tradicionales. Como son los mtodos aprendidos en cursos
elementales como ecuaciones diferenciales. Sin embargo, los sistemas actuales tienden
cada ves a ser ms complejos, por lo que la solucin de los sistemas de ecuaciones
diferenciales principales o la regin donde se puede localizar la solucin demandan el
uso de un mtodo de aproximacin o un mtodo numrico, para as extraer la
informacin relacionada con el comportamiento del sistema.
Hasta hace poco tiempo, un mtodo muy conocido y muy empleado por ingenieros,
fsicos y matemticos para analizar problemas complejos que involucran sistemas de
ecuaciones diferenciales, ordinales y parciales era el mtodo de diferencia finita o
diferencial finito (finite difference method) utilizado en una de sus mltiples formas.
Durante el ltimo cuarto de siglo, el mtodo del elemento finito (FEM, por sus siglas en
ingls) se ha convertido en una alternativa viable que en muchas aplicaciones tienes
ventajas sobre el mtodo antes mencionado.
El mtodo del elemento finito es una tcnica numrica que soluciona o se aproxima a
una solucin a un sistema de ecuaciones diferenciales relacionadas con un problema de
carcter fsico o ingenieril.
Como en el uso del mtodo de diferencia finita. El mtodo del elemento finito requiere
que el problema se encuentre definido en un espacio geomtrico, o dominio, para as ser
subdividido en un nmero finito de regiones pequeas, formando una especie de red o
malla (mesh). Sin embargo, existe una diferencia importante entre estos dos mtodos.
En el mtodo de diferencia finita, la malla consiste de filas y columnas de lneas
ortogonales, mientras que en el mtodo de elemento finito, cada divisin es nica y no
necesariamente ortogonal. (figura 11).


Por ejemplo, es posible utilizar tringulos o cuadrilteros en dos dimensiones, y
tetraedros o hexaedros en tres dimensiones. Sobre cada elemento finito, las variables
desconocidas (velocidad, temperatura, etc.) son aproximadas mediante el uso de
funciones, funciones que pueden ser lineales o de orden ms alto como polinomial.
Dicho orden depende de la localizacin de los puntos geomtricos o vrtices que
definen la forma del elemento finito (nodos).
En contraste con el procedimiento del mtodo diferencial, consistente en
discretizaciones convencionales, las ecuaciones principales en el mtodo del elemento
finito son integradas para cada uno de los elementos que componen la malla y la
solucin se suma alrededor de todo el dominio. Como consecuencia de estas
operaciones, se obtiene un conjunto de ecuaciones lineales finitas, que estn en funcin
de un conjunto de parmetros desconocidos para cada elemento.
La solucin de estas ecuaciones se logra mediante mtodos algebraicos.
La historia del mtodo del elemento finito es particularmente interesante,
principalmente porque dicho mtodo surgi hace poco tiempo.
Las ideas bsicas del mtodo surgieron a partir de los avances en el rea de anlisis
estructural relacionado con la aviacin. En 1941, Hrenikoff present una solucin a un
problema de elasticidad usando lo que se conoca como mtodo de estructura de trabajo
(frame work method). Un mtodo alterno apareci en 1943, su creador fue Courant, y
usaba interpolaciones polinomiales sobre regiones triangulares para la solucin de
problemas relacionados con la torsin.
A mediados de los aos 50as comenzaron a aparecer los resultados producto del
esfuerzo dedicado a resolver problemas continuos de elasticidad con el uso de pequeos
elementos para describir el comportamiento de una barra elstica.
El progreso de estos mtodos continuo a lo largo de esa dcada, y en 1960 surgi el
nombre de mtodo de elemento finito formalmente en una publicacin de Clough.
En su edad temprana, el uso de los elementos finitos estaba restringido a la aplicacin
de estas tcnicas para problemas estructurales y para cuestiones relacionadas. Sin
embargo, la versatilidad del mtodo y sus fuertes bases matemticas inherentes a dicho
mtodo fueron reconocidas y aplicadas a otras reas, siendo Zienkiewicz y Cheung los
primeros en aplicar, o al menos en publicar el mtodo del elemento finito aplicado a
problemas de campo (conductores de calor, fluidos irrotacionales, etc.).
El desarrollo posterior se dio en forma sbita, en tan diversas aplicaciones como son
problemas de transferencia de calor, como por ejemplo lo discutido en Huebner (1975).
Han sido publicados numerosos artculos y textos, como revisiones y descripciones del
mtodo, Bickford (1990), o Zienkiewicz y Taylor (1989) por mencionar algunos.
Para ahondarse en el mtodo del elemento finito es necesario retomar y ms que nada
reafirmar algunos conceptos bsicos necesarios en el desarrollo de ste.
1. Esfuerzos y equilibrio.
Un cuerpo tridimensional que ocupa un volumen, denotado por V, y que tiene una
superficie S se muestra en la siguiente figura (figura 12).

Cualquier punto en este cuerpo est definido por las coordenadas x, y, z. La superficie o
frontera est restringida a una regin, donde el desplazamiento est especificado. Sobre
una parte de la superficie se aplica una carga distribuida por unidad de rea, esta fuerza
de traccin se denomina T. Bajo esta carga, el cuerpo se deforma. La deformacin de un
punto [ ] ) , , ( x
T
z y x est dada por los tres componentes de su desplazamiento:
[ ]
T
, , u w v u
La fuerza distribuida por unidad de volumen, por ejemplo, el peso por unidad de
volumen es el vector f definido por
[ ]
T
, , f z y x f f f
La fuerza aplicada que acta sobre el volumen elemental dV se observa en la figura. La
fuerza superficial T puede definirse por los valores de sus componentes sobre puntos
definidos en la superficie.
[ ]
T
, , T z y x T T T
Ejemplos de traccin son fuerzas distribuidas de contacto y accin de presin. Una
carga P que acta en un punto i se representa por sus tres componentes
[ ]
T
, , P
i
z y x i P P P
Los esfuerzos que actan sobre el volumen dV se muestran en la siguiente figura (figura
13).

Donde el volumen dV se reduce a un punto, el esfuerzo tensor se representa mediante la
colocacin de sus componentes en una matriz simtrica (3 x 3). Sin embargo,
representamos el esfuerzo por los seis componentes independientes de la siguiente
forma:
[ ]
T
, , , , , xy xz yz z y x
donde z y x , , son esfuerzos normales y , , , xy xz yz son esfuerzos cortantes. Ahora bien,
si se considera que el volumen mostrado en la figura se encuentra en equilibrio es
posible obtener las ecuaciones de equilibrio para el elemento. Primero se obtienen las
fuerzas por la multiplicacin de los esfuerzos por sus reas correspondientes.
Escribiendo

0 x F ,

0 y F y

0 z F , y estableciendo dV = dx dy dz.
Las ecuaciones sern...
0 +

x
xz xy x
f
z y x


0 +

y
yz y xy
f
z y x


0 +

z
z yz xz
f
z y x


2. Condiciones de frontera.
En el anlisis de la figura 2, es posible observar que hay condiciones de frontera para el
desplazamiento y condiciones de carga en la superficie. Si u est definida en parte de la
frontera denotada por u S , se tiene
u S en 0 u
tambin es posible considerar las condiciones de frontera como u = a, donde a es un
desplazamiento previamente definido.
Considrese el equilibrio de un tetraedro ABCD, que se muestra en la siguiente figura.
(figura 14)

Donde DA, DB, y DC son paralelos a los ejes x, y y z respectivamente, y el rea ABC,
denotada por dA se halla en la superficie. Si [ ]
T
, , n z y x n n n es la unidad normal a dA,
entonces el rea dA n BDC x , rea dA n ADC y , y el rea . dA n ADB z
La consideracin de equilibrio a lo largo de los tres ejes da como resultado...
z z z y yz x xz
y z yz y y x xy
x z xz y xy x x
T n n n
T n n n
T n n n
+ +
+ +
+ +




Estas condiciones deben de cumplirse en la frontera T S , donde las fuerzas de traccin
son aplicadas. En esta descripcin, las cargas puntuales deben ser tratadas como cargas
distribuidas sobre reas pequeas pero finitas.
3. Relacin deformacin-desplazamiento.
Las deformaciones del elemento se representan en un vector que corresponde a los
esfuerzos.
[ ]
T
, , , , , xy xz yz z y x
donde z y x y , , son las deformaciones normales y xy xz y y , , son las deformaciones por
esfuerzo cortante.
4. Relacin esfuerzo-deformacin.
Para materiales que presentan un comportamiento elstico lineal, la relacin esfuerzo-
deformacin viene de la generalizacin de la regla de Hooke. Para materiales isotpicos,
las dos propiedades de los materiales son el mdulo de Young o mdulo de elasticidad
E y la razn de Poisson v. Considerando un cubo elemental dentro del cuerpo, la ley de
Hooke da como resultado.
E E
v
E
v
E
v
E E
v
E
v
E
v
E
z y x
z
z y x
y
z y x
x

+
+


G
G
G
xy
xy
xz
xz
yz
yz

Donde G es el mdulo de rigidez, y est dado por.
) 1 ( 2 v
E
G
+

de la ley de Hooke, se nota que
) (
) 2 1 (
z y x z y x
E
v
+ +

+ +
si se substituye por ) ( z y + en la primera ecuacin producto de la regla de Hooke, se
obtienen las relaciones inversas D
D representa la matriz del material (6 x 6) dada por:
1
1
1
1
1
1
1
1
]
1

v
v
v
v v v
v v v
v v v
v v
E
5 . 0 0 0 0 0 0
0 5 . 0 0 0 0 0
0 0 5 . 0 0 0 0
0 0 0 1
0 0 0 1
0 0 0 1
) 2 1 )( 1 (
D
Dentro de la relacin esfuerzo-deformacin existen adems casos especiales que se
describen brevemente a continuacin.
Unidimensional. En una dimensin, hay esfuerzo normal a lo largo de x y la
deformacin correspondiente. La relacin esfuerzo-deformacin es simplemente
E
Bidimensional. En dos dimensiones, los problemas son modelados como esfuerzos
planos y deformaciones en el plano.
Esfuerzos planos. Un cuerpo plano delgado, que est sometido a una carga en el
plano en el borde de su superficie, se dice est sometido a esfuerzo plano. Un anillo
sujetado a presin en una flecha es un ejemplo de este tipo de carga. (Figura 15).

En este caso los esfuerzos 0 y , , yz xz z , por lo que la ley de Hooke da
) (
) 1 ( 2
y x z
xy xy
y x
y
y x
x
E
v
E
v
E E
v
E
v
E

+
+

+


las relaciones inversas son

'

1
1
1
1
]
1

'

xy
y
x
xy
y
x
v
v
v
v
E

2
1
0 0
0 1
0 1
1
2
que es usada como D
Deformacin planar. Si un cuerpo largo de seccin transversal constante est
sujeto a una carga transversal a lo largo de su longitud, en pequeo espesor en la zona
con carga aplicada, como lo muestra la figura. (figura 16)

Puede ser tratada como sujeta a deformacin en el plano. En este caso yz xz z , , se
asumen como 0. El esfuerzo z puede ser un valor diferente de cero. La relacin
esfuerzo-deformacin puede obtenerse directamente de la matriz del material D.

'

1
1
1
1
]
1

'

xy
y
x
xy
y
x
v
v v
v v
v v
E

2
1
0 0
0 1
0 1
) 2 1 )( 1 (

en este caso D es una matriz (3 x 3), ya que representa tres esfuerzos y tres
deformaciones.
Cuerpos anisotrpicos, con orientacin uniforme, pueden considerarse usando la matriz
D apropiada para el material.
Existen otros muchos conceptos bsicos que son parte fundamental en el mtodo del
elemento finito, sin embargo, los mencionados antes son probablemente los ms
importantes. A continuacin se mencionan otros, sin profundizar en su definicin.
5. Efecto de la temperatura. Si existe un incremento en la temperatura ) , , ( z y x T con
respecto a la temperatura original, la deformacin asociada se puede considerar
fcilmente. Para los materiales isotpicos, el incremento de temperatura da como
resultado una deformacin uniforme, la cual depende del coeficiente de expansin lineal
del elemento. Dicho coeficiente representa el cambio de longitud de temperatura. Se
asume como constante dentro del rango de variacin de la temperatura. La deformacin
producto de este cambio de temperatura no provoca ningn esfuerzo cuando el cuerpo
tiene libertad de deformarse.
6. Energa potencial y equilibrio. Mtodo de Rayleigh-Ritz. En mecnica de slidos, el
problema es determinar el desplazamiento u de un cuerpo tridimensional con volumen
V, mediante una solucin que satisfaga las ecuaciones de equilibrio del elemento.
Ntese que los esfuerzos estn relacionados con las deformaciones, que a su vez, estn
relacionadas con los desplazamientos. Esto lleva a la necesidad de resolver ecuaciones
diferenciales de segundo orden. La solucin de este conjunto de ecuaciones se refiere
generalmente a una solucin exacta. Dichas soluciones exactas son disponibles para
geometras simples y condiciones de carga. Pero para geometras complejas y en
general para las condiciones de carga y de frontera, la obtencin de estas soluciones es
una labor casi imposible. Algunos mtodos de aproximacin utilizan la energa
potencial o mtodos variacionales, en las cuales las condiciones en las funciones son
ms flexibles.
La energa potencial es la suma de la energa total de deformacin y el trabajo potencial.
Mtodo de Rayleigh-Ritz. La energa potencial puede utilizarse para encontrar una
solucin aproximada, el mtodo Rayleigh-Ritz involucra la construccin de un supuesto
campo de desplazamiento, donde dicho desplazamiento debe ser cinematicamente
admisible. Esto es, que las componentes del desplazamiento deben satisfacer
condiciones de frontera especficas.
7. Mtodo de Galerkin. El mtodo de Galerkin usa un conjunto de ecuaciones en el
desarrollo de una forma integral. Usualmente se presenta como uno de los mtodos
residuales de mayor importancia. Se considera una representacin de una ecuacin
general de una ecuacin rectora en una regin V.
Lu = P
Para una barra unidimensional por ejemplo, la ecuacin principal es la ecuacin
diferencial
0 ) (
dx
du
EA
dx
d

si se considera L como el operador
) (
dx
d
EA
dx
d

operando en u.
La solucin exacta necesita ser vlida para cada punto x. Si se busca una solucin
aproximada u
~
, introduce un error (x) que recibe el nombre de residual:
P u L x
~
) (
los mtodos aproximados se mueven alrededor determinando el residual relativo a una
funcin W, a cero.

V
i dV p u L W 0 )
~
( i = 1 a n
las funciones i W son elegidas de acuerdo a las funciones base utilizadas para la
construccin de u
~
. Si u
~
est representado por

n
i
i iG Q u
1
~

Donde n, a 1 , i Gi son funciones base (usualmente funciones polinomiales de x, y, z).
Se escogen las funciones de ms peso a ser una combinacin de funciones base i G .
Especficamente, se considera una funcin arbitraria, dada por

n
i
i iG
1

8. Principio de Saint Venant. En ocasiones es necesario hacer aproximaciones en la
definicin de las condiciones de frontera para representar una interfase soporte-
estructura.
Considrese una viga en cantiliver, libre en uno de sus extremos y sujeta en el otro a
una columna mediante remaches. Algunas cuestiones surgen acerca de la rigidez de la
unin remachada, y acerca de la igualdad en las condiciones de frontera de cada uno de
los elementos situados en la seccin transversal del extremo fijo. Saint Venant
consider el efecto que diferentes aproximaciones tienen sobre la solucin total del
problema. Este principio establece que mientras las aproximaciones realizadas sean
estticamente equivalentes, las soluciones resultantes sern vlidas siempre y cuando el
enfoque sea en una regin suficientemente alejada del soporte. Esto quiere decir que, las
soluciones pueden variar significativamente solo dentro de la proximidad del soporte.
Hasta ahora, se han discutido las bases necesarias para el mtodo del elemento finito,
por lo que se observarn algunas aplicaciones de este mtodo. Que por su versatilidad es
aplicado cada vez a ms reas. Adems se incluyen breves ejemplo para algunas de
stas aplicaciones.

3.1 Aplicaciones del mtodo del elemento finito en ingeniera.
1. Mecnica de slidos. Estructuras.
La industria aerospacial usa el mtodo del elemento finito para determinar la respuesta
esttica y dinmica de los aviones y las naves espaciales ante la gran variedad de
ambientes y condiciones que encuentran durante su funcionamiento.
Para realizar un anlisis los elementos se dividen en reiteradas ocasiones para el anlisis
de los componentes.
2. Mecnica de fluidos.
Un ejemplo de una aplicacin del mtodo del elemento finito en la mecnica de fluidos
es el problema del flujo de aire alrededor del ala de un avin. Como se muestra en la
imagen de la seccin transversal de un ala. (Figura 17).


Lo ms interesante de esto son las fuerzas de flotacin y las fuerzas que actan en
contra del avance del cuerpo (drag force). La regin de flujo cercana al ala se divide en
elementos como se indica. La solucin del modelo de elemento finito permite el
cmputo de las fuerzas previamente mencionadas.
3. Transferencia de calor.
Un ejemplo de la aplicacin del mtodo del elemento finito a esta rea es en torno a un
motor de turbina de gas. La temperatura de funcionamiento alcanza niveles muy altos en
el motor y ciertos componentes deben ser enfriados para que soporten y puedan seguir
funcionando correctamente. Las aspas en el rotor generalmente contienen cavidades
para que permitan el paso de una corriente de aire refrigerante hacia el interior de la
aspa para cumplir con el objetivo de evitar el sobrecalentamiento.
Un modelo de elemento finito para un aspa se muestra en la siguiente figura (figura 18).

Este modelo se usa para determinar el nmero apropiado de cavidades, su tamao y su
localizacin para un enfriamiento apropiado del aspa y el motor.
4. Filtracin de agua subterrnea.
Una importante aplicacin en el rea geofsica es el problema de la filtracin de agua
subterrnea. Una tpica situacin de este problema se muestra en la siguiente imagen,
(figura 19).

Donde el agua es retenida detrs de un dique impenetrable. La tarea es determinar la
cantidad de agua que se pierde debido a la filtracin bajo el dique en la tierra. La malla
para el modelo de elemento finito del suelo se indica tambin.
El mtodo del elemento finito es aplicado fcilmente y comnmente a una enorme
cantidad de problemas fsicos y matemticos como los que se ha indicado hasta ahora.
Actualmente, existen muchos software que operan con el mtodo del elemento finito.
Estos software realizan todas las operaciones matemticas y todo el anlisis de los
modelos en cuestin de minutos o incluso segundos.
A continuacin se mencionan algunas empresas o casos que mediante el uso de un
software de elemento finito, han desarrollado productos o servicios que ofrecen al
mercado.
Actualmente existen cientos de industrias dedicadas al desarrollo de productos que han
sido desarrollados sobre la base del anlisis del elemento finito. As como empresas que
han desarrollado estos mismos programas, que cada vez ofrecen anlisis ms completos,
ms precisos y ms rpidos.
1. Una de las reas donde se ha aplicado sta tecnologa con excelentes resultados es en
la industria automotriz, ya que al permitir el modelado y anlisis de componentes que
van a estar sometidos a diversos tipos de carga, se pueden hacer mejoras y ahorrar
dinero que se gastara en mltiples pruebas. (Figura 20)

2. La industria aerospacial tambin ha aprovechado la tecnologa del FEM, para el
modelado de aviones, y sus partes, as como para el diseo de los componentes internos.
(figura 21).

3. en la industria manufacturera, la aplicacin del modelado por mtodo de elemento
finito, permite hacer un anlisis para manufacturabilidad, as como para disear
componentes que cumplan ciertos requisitos, y la principal ventaja radica en que este
anlisis se puede hacer antes de haber construido el elemento o siquiera un prototipo.
Por ejemplo, en la siguiente imagen. (Figura 22) se observa un anlisis de diseo para
fatiga, sobre un componente.

4. Con una gran cantidad de modelos fsicos y eficientes mtodos de solucin, FIDAP
es la herramienta de modelacin de flujos a elegir para aplicaciones como el
procesamiento de polmeros, recubrimientos de capas delgadas, flujos en biomedicina,
crecimiento de cristales semiconductores, metalurgia, procesamiento de vidrio y otras
reas relacionadas.
Basado en el mtodo de elemento finito, FIDAP ofrece una flexibilidad total en cuanto a
la definicin de mallas as como soluciones robustas y eficientes. Est disponible para
sistemas UNIX y Windows/NT y est optimizado para procesamiento en paralelo.
5. En el rea de Fisicoqumica de Fluidos se est haciendo investigacin sobre el tema
de las ecuaciones integrales mencionadas antes, desde varios aspectos:
Desarrollo de nuevos mtodos numricos confiables y eficientes que permitan resolver
estos problemas en tiempos de cmputo que sean accesibles. Dentro de este esfuerzo, se
desarroll un mtodo numrico que, basado en las tcnicas numricas de Elemento
Finito, ha demostrado ya en varios casos ser lo suficientemente poderoso para competir
con ventajas con los mtodos numricos tradicionales. Aplicaciones del mtodo
mencionado en a) en distintas situaciones de inters fsico. Ejemplo de esto es la
determinacin del comportamiento crtico de diferentes sistemas modelo, como el del
potencial de pozo cuadrado y el de Yukawa. Generalizacin del mtodo numrico
mencionado en a) para el tratamiento de sistemas no uniformes y en particular para el
tratamiento de la doble capa elctrica segn el modelo primitivo de un electrolito.
Bsqueda de soluciones analticas de las ecuaciones integrales de la fsica de lquidos.
Como se mencion antes, las soluciones analticas son muy contadas. En este campo se
hizo un avance importante al encontrar la solucin analtica de la ecuacin integral de la
aproximacin esfrica media para el modelo de potencial intermolecular de Yukawa.
6. Aplicaciones del Mtodo del elemento Finito a la biomecnica ortopdica.
As como estos, podemos encontrar cientos de casos ms, sin embargo, lo importante es
notar como la aplicacin de este mtodo, principalmente utilizando programas que
trabajan bajo dicho mtodo, representa una mejora cada da en campos ms variados.
Para este proyecto, el software a utilizar es Algor , el cul se describe en el captulo 4.