Universidad De Oriente

Ncleo de Anzotegui
Escuela de Ingeniera y Ciencias Aplicadas
Departamento de Ingeniera Civil
Ctedra: Anlisis estructural II

Mtodo De
Elementos Finitos

BNA, 03/05/201

INDICE
INTRODUCCIN..............................................................................................................................4
1.- CONCEPTOS BSICOS DEL MTODO DE ELEMENTOS FINITOS (MEF)........................5
1.1- Resumen..................................................................................................................................5
1.2. Desarrollo................................................................................................................................5
1.3.-Objetivos del mtodo...............................................................................................................6
2.- ETAPAS BSICAS EN LA UTILIZACIN DEL MTODO DE LOS ELEMENTOS
FINITOS.............................................................................................................................................9
2.1.- Definicin del problema y su dominio.................................................................................10
2.2- Discretizacin del dominio....................................................................................................10
2.3.- Identificacin de la(s) estado................................................................................................11
2.4.- Formulacin del problema...................................................................................................12
2.5.- establecimiento de los sistemas de referencia......................................................................12
2.6.- construccin de las funciones de aproximacin de los elementos.......................................14
2.7.- determinacin de las ecuaciones a nivel de cada elemento..................................................15
2.8.- transformacin de coordenadas............................................................................................15
2.9.- ensamblaje de las ecuaciones de los elementos....................................................................15
2.10.- introduccin de las condiciones de contorno.....................................................................16
3.11. Solucin del sistema de ecuaciones resultantes...................................................................16
2.11.1. Mtodos directos...........................................................................................................16
2.11.2.- mtodo iterativo...........................................................................................................17
2.12. Interpretacin de los resultados...........................................................................................17
3.- ESQUEMA DE METODOLOGA DE CLCULO..................................................................18
4. EJERCICIO..................................................................................................................................18
CONCLUSION................................................................................................................................22
BIBLIOGRAFIA..............................................................................................................................23

INTRODUCCIN
El Mtodo de Elementos Finitos (MEF) ha evolucionado desde su invencin en los
aos sesenta hasta convertirse en una herramienta indispensable en cualquier mbito ingenieril.
Actualmente hay abiertas muchas lneas de investigacin para mejorar y depurar las tcnicas y
teoras utilizadas, entre las cuales se encuentra el refinamiento de malla adaptativo.
Como toda herramienta, debemos entender que el uso del MEF requiere una persona
cualificada y con experiencia en el campo. El refinamiento de malla adaptativo busca facilitar el
uso del MEF realizando tareas esenciales en todo modelo, tales como la estimacin del error
cometido y la correccin de la malla en zonas que lo requieran. Todava no existen programas
comerciales con mtodos de refinamiento de malla adaptativos, de modo que sigue siendo
imperativa la presencia de un profesional experimentado para utilizar los programas en desarrollo.

1.- CONCEPTOS BSICOS DEL MTODO DE ELEMENTOS FINITOS (MEF)

1.1- Resumen
Muchos problemas de importancia prctica que frecuentemente aparecen en la
ingeniera, resultan de una complejidad matemtica tal que, aunque la deduccin de las ecuaciones
diferenciales que gobiernan tales problemas no resulta muy difcil, su solucin por mtodos
exactos de anlisis, aun despus de introducir algunas hiptesis simplificadoras, no se logra si no
para ciertos problemas de geometra, condiciones de contorno y/o sistemas de carga muy
particulares. Por esto, aunque este tipo de solucin es la que mas informacin proporciona sobre el
comportamiento de las variables involucradas en un problema dado, se debe recurrir a los mtodos
numricos, los cuales permiten elaborar un anlisis y diseos con alto grado de sofistificacin y
precisin
1.2. Desarrollo
El Mtodo de Elementos Finitos fue al principio desarrollado en 1943 por Richard
Courant, quien utiliz el mtodo de Ritz de anlisis numrico y minimizacin de las variables de
clculo para obtener soluciones aproximadas a un sistema de vibracin. Poco despus, un
documento publicado en 1956 por M. J. Turner, R. W. Clough, H. C. Martin, y L. J. Topp
estableci una definicin ms amplia del anlisis numrico. El documento se centr en la rigidez
y deformacin de estructuras complejas. Con la llegada de los primeros ordenadores instaura el
clculo matricial de estructuras. ste parte de la discretizacin de la estructura en elementos
lineales tipo barra de los que se conoce su rigidez frente a los desplazamientos de sus nodos. Se
plantea entonces un sistema de ecuaciones resultado de aplicar las ecuaciones de equilibrio a los
nodos de la estructura. Este sistema de ecuaciones se esquematiza de la siguiente manera:

Donde las incgnitas son los desplazamientos en los nodos (vector u) que se hallan
a partir de las "fuerzas" o "solicitaciones" en los nodos (vector ) y de la rigidez de las barras
(matriz de rigidez

). Conocidos dichos desplazamientos es posible determinar los esfuerzos en

las barras. La solucin obtenida es exacta

1.3.-Objetivos del mtodo
El Mtodo de Elementos Finitos permite obtener una solucin
numrica aproximada sobre un cuerpo, estructura o dominio sobre el que estn definidas
ciertas ecuaciones diferenciales en forma dbil o integral que caracterizan el comportamiento
fsico del problema dividindolo en un nmero elevado de subdominios no intersectantes entre s
denominados elementos finitos. El conjunto de elementos finitos forma una particin del
dominio tambin denominada discretizacin. Dentro de cada elemento se distinguen una serie de
puntos representativos llamados nodos. Dos nodos son adyacentes si pertenecen al mismo
elemento finito; adems, un nodo sobre la frontera de un elemento finito puede pertenecer a varios
elementos. El conjunto de nodos considerando sus relaciones de adyacencia se llama malla.
Los clculos se realizan sobre una malla de puntos (llamados nodos), que sirven a
su vez de base para discretizacin del dominio en elementos finitos. La generacin de la malla se
realiza usualmente con programas especiales llamados generadores de mallas, en una etapa previa
a los clculos que se denomina pre-proceso. De acuerdo con estas relaciones de adyacencia o
conectividad se relaciona el valor de un conjunto de variables incgnitas definidas en cada nodo y
denominadas grados de libertad. El conjunto de relaciones entre el valor de una determinada
variable entre los nodos se puede escribir en forma de sistema de ecuaciones lineales (o

linealizadas). La matriz de dicho sistema de ecuaciones se llama matriz de rigidez del sistema. El
nmero de ecuaciones de dicho sistema es proporcional al nmero de nodos.
Tpicamente el anlisis de los elementos finitos se programa computacionalmente
para calcular el campo de desplazamientos y, posteriormente, a travs de relaciones cinemticas y
constitutivas las deformaciones y tensiones respectivamente, cuando se trata de un problema de
mecnica de slidos deformables o ms generalmente un problema de mecnica de medios
continuos. El mtodo de los elementos finitos es muy usado debido a su generalidad y a la
facilidad de introducir dominios de clculo complejos (en dos o tres dimensiones). Adems el
mtodo es fcilmente adaptable a problemas de transmisin de calor, de mecnica de fluidos para
calcular campos de velocidades y presiones (mecnica de fluidos computacional, CFD) o
de campo electromagntico. Dada la imposibilidad prctica de encontrar la solucin analtica de
estos problemas, con frecuencia en la prctica ingenieril los mtodos numricos y, en particular,
los elementos finitos, se convierten en la nica alternativa prctica de clculo.
Hace aproximadamente cuarenta aos esta tecina se ha destacado, tanto por su
capacidad para modelar dominios irregulares, condiciones de contorno, no-linealidades, como por
la facilidad en la seleccin del mecanismo de aproximacin de las variables involucradas en un
problema especifico. Visto globalmente, la solucin de un problema dado se puede esquematizar
tal como se muestra en la fig. 1.1.

Fig.1.1 diagrama esquemtico del modelaje matemtico de un problema.

___________ proceso deductivo / analtico.
----------- fuentes de discrepancia: realidad / modelo matemtico

2.- ETAPAS BSICAS EN LA UTILIZACIN DEL MTODO DE LOS

ELEMENTOS FINITOS.
Independientemente de la naturaleza fsica del problema, el anlisis del mismo
mediante MEF sigue los siguientes pasos:

1.- Definicin del problema y su dominio.

2.- Discretizacin del dominio.
3.- Identificacin de la(s) variables de estado.
4.- Formulacin del problema.
5.- Establecimiento de los sistemas de referencia.
6.- construccin de las funciones de aproximacin de los elementos.
7.- determinacin de las ecuaciones a nivel de cada elemento.
8.- transformacin de coordenadas.
9.- ensamblaje de las ecuaciones de los elementos.
10.- introduccin de las condiciones de contorno.

11. Solucin del sistema de ecuaciones resultantes.

12. Interpretacin de los resultados.

2.1.- Definicin del problema y su dominio.

El anlisis de un problema dado va el mef, tiene implcito tres tipos de
aproximacin. La primera se relaciona con la definicin del dominio (fsica y geomtrica) del
problema, las otras dos estn asociadas a la discretizacin da las ecuaciones gobernantes, y a los
algoritmos empleados en la solucin del sistema de ecuaciones algebraicas simultaneas resultante.
Las aproximaciones usadas en la definicin de las caractersticas fsicas de las
diferentes regiones del dominio, depende fundamentalmente del tipo de problema a resolver. Sin
embargo, la definicin geomtrica del dominio, requiere el establecimiento de ejes coordenadas
globales en referencia a los cuales se describen las coordenadas de ciertos puntos (nodos), los
cuales, a su vez, definen las ecuaciones de las lneas, superficies y/o volumen de elementos. Este
sistema no necesita ser rectangular y cartesiano, para algunos problemas especficos, resulta mas
adecuado utiliza algn tipo de sistema coordenado curvilneo.
2.2- Discretizacin del dominio
Puesto que usualmente el problema est definido sobre un dominio continuo, las
ecuaciones gobernantes de un problema, con excepcin de las condiciones de contorno, son

10

validas tanto en todo el dominio como en cualquier parte del. Esto permite idealizar el dominio a
travs de regiones de tamao finito (elementos), interconectados de diferente forma y tamao, tal
como se muestra en la fig. 1.2. Esta forma de discretizacin introduce ciertas aproximaciones. Sin
embargo, colocando un nmero suficiente de elementos (o elementos de orden superior), se podr
reproducir el dominio tan aproximadamente cuando queramos
Aun cuando es cierto que, en general, reduciendo el tamao de los elementos se
obtienen mejores resultado, tambin es cierto que un refinamiento excesivo conduce a grandes
sistemas de ecuaciones, lo cual puede tornarse imprctico desde un punto de vista computacional.
De tal modo que, inicialmente, se debe entonces responder las siguientes pregunta: Cul es el
tipo de elemento ms eficiente y cul ser el tamao adecuado?

Fig. 1.2 distribucin del dominio con diferentes elementos finitos

2.3.- Identificacin de la(s) estado

11

Hasta el momento no se ha hecho referencia a la naturaleza fsica del problema ya

que las etapas anteriores son comunes a cualquier tipo de problema, ya sea este de transferencia de
calor, de la mecnica de los fluidos, de la mecnica de los slidos, etc. A continuacin, y para cada
problema en particular, la descripcin matemtica del fenmeno fsico conducir al
correspondiente problema de valor de contorno, el cual contendr las variables de estado asociadas
al mismo. Estas variables se relacionan entre s a travs de las ecuaciones constitutivas, las cuales
representan una expresin matemtica de una ley fsica en particular. La tabla 1.1 muestra varios
problemas con las variables de estado asociadas, y las correspondientes ecuaciones constitutivas,
muchos problemas reales involucran el anlisis de dos o ms tipos de problemas mostrados en
dicha tabla, de modo conjunto y simultaneo.
2.4.- Formulacin del problema
Frecuentemente, un problema fsico esta formulado a travs de un conjunto de
ecuaciones diferenciales con sus correspondientes condiciones de contorno, o mediante una
ecuacin integral (un funcional) sujeto a un requerimiento estacionario (mximo o mnimo). En el
primer caso se dice que el problema fsico est referido a su forma diferencial y en el segundo, a su
forma variaciones. En ambos casos se llega al mismo resultado. En este texto se presentaran las
dos formulaciones como forma de establecer las ecuaciones de los elementos.
2.5.- establecimiento de los sistemas de referencia
Adems de los ejes globales de referencia del sistema completo, existen dos
importantes razones para seleccionar, adicionalmente, un sistema de referencia local para los
elementos: la facilidad con la que se construyen las llamadas funciones de forma de los elementos
y la facilidad con la que se integra en el interior de los mismos, con respecto al sistema local de

12

cada elemento en particular. Sin embargo, puesto que los elementos se ensamblan en el sistema
global de referencia, este paso introduce a una transformacin de coordenadas.
Tabla 1.1 Descripcin matemtica de varios problemas de valor de contorno

A pesar que todos los clculos en el mef se pueden realizar directamente en el

sistema global, este procedimiento es muy complicado para cualquier problema de inters practico
y, puesto que la transformacin de coordenadas entre cualesquiera dos sistemas de coordenados
est bien definida y es una operacin matemticamente sencilla, se deben deducir las ecuaciones
de los elementos con relacin a su sistema local de referencia el cual puede ser cartesiano o
curvilneo, dependiendo de la forma de un elemento dado. En la fig.1.3 se muestra un elemento
bidimensional y los sistemas, global y local de referencia

13

Fig. 1.3 Sistema de referencias usadas en el mtodo de los elementos finitos

(s) Sistema local de referencia; (b) sistema global de referencia

2.6.- construccin de las funciones de aproximacin de los elementos

Una vez que se han seleccionado el sistema de coordenado local y la(s) variables(s)
de estado, estas pueden ser aproximadas de diferentes formas. En el mef, la aproximacin tanto del
dominio del problema como de las variables involucradas en el mismo, se realiza mediante
funciones algebraicas. Si el elemento es plano o de lados rectos, las coordenadas de los nodos
primarios, definirn la forma exacta del mismo. Debido a esto, la discretizacin del dominio
muchas veces se realiza mediante elementos de lado recto. Sin embargo, para algunos problemas
los elementos, pueden producir errores inaceptables y la discretizacin debe ser realizada con
elementos de orden superior como lo muestra la fig. 1.4

14

Nodos para definir la forma del elemento (geometra).

o Nodos para aproximar las variables de estado (fsica).

Fig.1.4 ejemplo de elementos finitos subparametricos, isoparametricos,
superparametrico

2.7.- determinacin de las ecuaciones a nivel de cada elemento

A esta altura el modelaje del problema, se ha completado. Usando algn modelo
matemtico (mtodo de residuos pesados, trabajo virtual, mtodo de energa, etc.), se debe
establecer, a continuacin, sobre cada elemento, las ecuaciones discretas del problema continuo.
Este paso involucra determinar la llamada matriz de rigidez de cada elemento con respecto a su
sistema local de referencia.
2.8.- transformacin de coordenadas
Una vez determinadas las matrices de rigidez de todos los elementos que conforman
la discretizacin del dominio del problema, y antes de proceder al ensamblaje de todas las
matrices, para as obtener el comportamiento de todo el sistema, es necesario realizar la
transformacin de coordenadas, que permita transformar matrices de rigidez de los elementos,
desde sus respectivos ejes coordenados locales, al sistema global de referencia.

15

2.9.- ensamblaje de las ecuaciones de los elementos

El ensamblaje de las matrices de las ecuaciones de los elementos, se realiza de
acuerdo con la configuracin topolgica de los mismos, despus que estas han sido transformadas
al sistema global de referencia. Dicha configuracin se obtiene a travs del establecimiento de una
relacin entre la numeracin local de los nodos de los elementos, y la numeracin global de los
mismos. El ensamblaje se efecta considerando nicamente los nodos de las interfaces, los cuales
son comunes a los elementos adyacentes. La matriz resultante se denomina matriz global del
sistema.
2.10.- introduccin de las condiciones de contorno
En este paso se introducen las condiciones de contorno en la matriz global del
sistema, con lo cual esta matriz se podr reducir o condensar a su forma final, aun cuando en
algunos casos se prefiere, para no aadir nuevos algoritmos a la solucin del problema, dejar e
sistema global con su tamao inicial.
3.11. Solucin del sistema de ecuaciones resultantes
El paso final en la solucin de un problema por la mef, lo constituye la resolucin del
sistema de ecuaciones simultneas resultante. Debido a la naturaleza determinantica del mef, los
procedimientos de solucin de dicho sistema se pueden clasificar en dos grupos:
2.11.1. Mtodos directos
Tales como: mtodo de Gauss y factorizacin de Cholesky, los cuales son los ms
utilizados para sistema de ecuaciones pequeos.

16

El mef se desarrollo, al inicio de los 50 a partir del llamado mtodo directo asociado
al clculo estructural, el cual fue ampliamente usado en la solucin de diversos problemas
estructurales relacionados con la industria aeronutica. Las relaciones entre los desplazamientos y
las fuerzas que los originan, se expresan mediante un conjunto de ecuaciones, dando el origen a lo
que se dio en llamar matriz de rigidez de cada elemento estructural, y se desarrollaron tcnicas
para realizar el ensamblaje de estas matrices en una matriz global, que expresara el
comportamiento de todas las estructuras en estudio. Desafortunadamente este mtodo es difcil de
aplicar en problemas bidimensionales y tridimensionales, los cuales son, precisamente, los casos
donde el mef es ms til.
2.11.2.- mtodo iterativo
Tales como el mtodo de Gauss-Seidel y el de Jacobi, los cuales son ms apropiados
para sistemas de grandes ordenes. En estos mtodos, el tiempo de la solucin es considerablemente
menor que en los mtodos directos.
2.12. Interpretacin de los resultados
Con la resolucin del sistema se obtienen los valores aproximados de la(s)
variable(s) en los puntos discretos (nodos) del dominio. Generalmente estos valores son
interpretados y usados en el clculo de otras cantidades fsicas, tales como los esfuerzos,
deformaciones, el flujo del calor, etc.

17

3.- ESQUEMA DE METODOLOGA DE CLCULO

4. EJERCICIO.
Para este caso usaremos el mtodo directo, este mtodo es usado la mayora de las
veces para anlisis esttico lineal, lo que quiere decir que tenemos la zona elstica y la zona
plstica, en este caso solo trabajaremos en la zona elstica, donde el esfuerzo es igual al modulo de
elasticidad por la deformacin unitaria, y la deformacin unitaria es igual a la deformacin total
sobre la longitud original de la pieza como se muestra en la fig. 1.5

18

Fig.1.5 zona platica y elstica

El primer elemento con el que trataremos ser un elemento finito ms simple, en

este caso ser el elemento barra, el cual tendr una longitud L y un rea de seccin A.
Lo siguiente ser comenzar la formulacin del elemento finito utilizando el mtodo directo. Antes
de hacer la formulacin utilizaremos una simplificacin muy simple por ejemplo de nuestro
elemento barra, como se muestra en la fig.1.6, al cual le hemos aplicado una fuerza lo
simplificaremos por un resorte, con rigidez k el cual al aplicar la fuerza tendr un
desplazamiento u de manera que la fuerza para el resorte es igual a K por U, donde K es
igual al modulo de elasticidad (E) por el rea de seccin (A) sobre la longitud de la pieza.

19

Fig.1.6 simplificado por un resorte

Iniciamos la formulacin, e imaginemos que tenemos un elemento al cual le aplicamos una
fuerza F1 y una fuerza F2. La fuerza F1 aplicada en el nodo 1 provocara un desplazamiento
U1, y la fuerza en el nodo 2, generara un desplazamiento U2. Fig.1.7

Fig.1.7 Formulacin
Lo primero que haremos, ser analizar este sistema haciendo quel desplazamiento en U1
sea igual 0 , para que esto sea posible, si mantenemos fijo el nodo 1, la fuerza F2
generara un desplazamiento U2 pero tendremos que tener una reaccin F1 en el sentido

20

opuesto. De manera que podemos escribir a la fuerza F2 como K por U2 y por equilibrio
diramos que F1 es igual (-F2) que a su vez es igual a (-K*U2). fig.1.8

FIG.1.8
Lo siguiente que haremos ser fijar al nodo 2,y procederemos igual que en el caso
anterior, nodo 1. Fig.1.9

Fig.1.9

21

Una vez hecho estos dos anlisis, podemos utilizar el mtodo de superposicin para unirlos
y quedara de la manera en que se muestra en la siguiente figura. Fig.1.10
Escribiendo esto es forma matricial tendramos lo siguiente. fig. 1.10

fig. 1.10

22

CONCLUSION

Los mtodos de los elementos finitos, de diferencias finitas, de volumen de control (bien sea
basado en diferencias finitas o elementos finitos) y de contorno, son apenas algunos, entre una
gran gama de mtodos numricos que se han venido desarrollando y usando exitosamente, en la
solucin de muchos problemas en distintas reas de la ciencia . aun cuando todos estos
mtodos constituyen una muy poderosa herramienta matemtica , no dejan de ser mtodos
aproximados , debindose tener por lo tanto un especial cuidado en su utilizacin , ya que la
calidad de las soluciones que se obtengan depende de varios factores , entre los cuales se
pueden destacar la distribucin de la discretizacin espacial de la regin en estudio , el tipo de
discretizacin en el tiempo en los problemas no permanentes , la aplicacin apropiada de las
condiciones de contorno , las correctas inclusiones en el modelo de las propiedades fsicas de
los materiales que intervienen en el problema etc. .el correcto posicionamiento de estos aspectos
requieren del sentido comn y alguna experiencia del anlisis , independientemente del mtodo
seleccionado .

23

El mtodo de los elementos finitos (MEF en castellano o FEM en ingls) es

un mtodo numrico general para la aproximacin de soluciones de ecuaciones diferenciales
parciales muy utilizado en diversos problemas de ingeniera y fsica.
El MEF est pensado para ser usado en computadoras y permite resolver
ecuaciones diferenciales asociadas a un problema fsico sobre geometras complicadas. El MEF se
usa en el diseo y mejora de productos y aplicaciones industriales, as como en la simulacin de
sistemas fsicos y biolgicos complejos. La variedad de problemas a los que puede aplicarse ha
crecido enormemente, siendo el requisito bsico que las ecuaciones constitutivas y ecuaciones de
evolucin temporal del problema a considerar sean conocidas de antemano.

BIBLIOGRAFIA