Universidad autónoma de Chiapas

Facultad de Ingeniería

"Método del elemento finito"

Materia:
Métodos numéricos

Docente:
MTRO. Saul Chanona Vázquez

Alumno:
Francisco Eduardo Domínguez Recillas

Segundo semestre
Grupo "C"

P á g i n a 1 | 14
 Índice

Portada____________________________________________________________01

Índice_____________________________________________________________02

Introducción________________________________________________________03

Método del elemento finito

Generalidades______________________________________________________04

Historia____________________________________________________________05

Delimitación del método_______________________________________________07

Aplicaciones________________________________________________________11

Conclusión_________________________________________________________13

Bibliografía_________________________________________________________14

P á g i n a 2 | 14
 Introducción

En el presente proyecto se hablará sobre el método del elemento finito que no es

mas que un método numérico para la resolución de ecuaciones diferenciales, que se
presentan en varios ámbitos, generalmente en problemas de la ingeniería y en la
física. Se dará a conocer sobre sus generalidades, algo de su historia, algunas
delimitaciones que tiene y aplicaciones, en lo que son las aplicaciones se dará un
ejemplo y links para saber mas sobre el tema con muchos mas ejemplos.

P á g i n a 3 | 14
 Método del elemento finito

Generalidades

El método del elemento finito (MEF en español o FEM en inglés) es un método

numérico para la resolución de ecuaciones diferenciales, utilizado en diversos
problemas de ingeniería y física.

El método se basa en dividir el cuerpo, estructura o dominio (medio continuo) sobre

el que están definidas ciertas ecuaciones integrales que caracterizan el
comportamiento físico del problema (figura 8.1), en una serie de subdominios no
intersectantes entre sí denominados elementos finitos. El conjunto de elementos
finitos forma una partición del dominio también llamada discretización.

Dentro de cada elemento se distinguen una serie de puntos representativos llamados

nodos. Dos nodos son adyacentes sí pertenecen al mismo elemento finito, además,
un nodo sobre la frontera de un elemento finito puede pertenecer a varios elementos
(figura 8.2). El conjunto de nodos considerando sus relaciones de adyacencia se
conoce como MODELADO DE PROCESOS DE MANUFACTURA 234 malla. Los
cálculos se realizan sobre una malla o discretización creada a partir del dominio con
programas generadores de mallas, en una etapa previa a los cálculos que se
denomina pre-proceso. De acuerdo con estas relaciones de adyacencia o
conectividad se relaciona el valor de un conjunto de variables incógnitas 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 número de ecuaciones de dicho sistema es
proporcional al número de nodos.
P á g i n a 4 | 14
Típicamente, el método del elemento finito se programa computacionalmente para
calcular el campo de desplazamientos y, posteriormente, a través de relaciones
cinemáticas y constitutivas, las deformaciones y tensiones respectivamente, cuando
se trata de un problema de mecánica de sólidos deformables o más generalmente un
problema de mecánica del medio continuo. El método de los elementos finitos es
muy usado debido a su generalidad y a la facilidad de introducir dominios de cálculo
complejos (en dos o tres dimensiones).

Además, el método es fácilmente adaptable a problemas de difusión del calor, de

mecánica de fluidos para calcular campos de velocidades y presiones, o de campo
electromagnético. Dada la imposibilidad práctica de encontrar la solución analítica de
estos problemas, con frecuencia, en la práctica ingenieril, los métodos numéricos y,
en particular, los elementos finitos se convierten en la única alternativa práctica de
cálculo.

Una importante propiedad del método es la convergencia, si se consideran

particiones de elementos finitos sucesivamente más finas, la solución numérica
calculada converge rápidamente hacia la solución exacta del sistema de ecuaciones

Historia

Diversos autores han considerado que Arquímedes (figura 8.3) utilizó un método
semejante al del elemento finito para determinar el volumen de algunos sólidos.
Aunque él calculó áreas, longitudes y volúmenes de objetos geométricos,
dividiéndolos en otros más sencillos y luego sumando sus contribuciones, el
concepto de aproximación variacional no se observa por ningún lado. La relación con
la definición de MEF es muy pobre. Se puede argumentar que la medida del volumen
P á g i n a 5 | 14
(área, longitud) de un objeto es una función escalar de su geometría. Cambiando
“medida” por energía y “objetos” por elementos en las líneas anteriores, la
descripción se aproxima a lo establecido por el MEF “la energía del sistema es igual
a la suma de la energía de cada elemento”. Sin embargo, Arquímedes necesitaba las
definiciones de derivada para realizar sus cálculos de energía y el Cálculo no fue
inventado sino hasta 20 siglos después.

En 1941, Hrenikoff presentó una solución para problemas elásticos usando el

“método de trabajo del marco”. En un artículo publicado en 1943, Courant usó
interpolación polinomial por partes sobre subregiones triangulares para modelar
problemas de torsión. Las ideas básicas del método del elemento finito se originaron
en el análisis estructural de las aeronaves. En el periodo de 1950-1962, Turner
trabajando para Boeing formula y perfecciona el Método por Rigidez Directo. Turner y
otros investigadores obtuvieron matrices de rigidez para armaduras, vigas y otros
elementos y presentaron sus resultados en 1956. Clough fue el primero en acuñar y
emplear el término elemento finito en 1960.

En los primeros años de la década de 1960, los ingenieros usaron el método para
obtener soluciones aproximadas en problemas de análisis de esfuerzos, flujo de
fluidos, transferencia de calor y otras áreas. Un libro de Argyris, publicado en 1955,
sobre teoremas de energía y métodos matriciales, cimentó métodos adicionales en
los estudios de elemento finito. El primer libro sobre elementos finitos por Zienkiewicz
y Cheng fue publicado en 1967. A finales de la década de 1960 y principios de la
siguiente, el análisis por elemento finito se aplicó a problemas no lineales y de
grandes deformaciones. El libro de Oden sobre continuos no lineales apareció en
1972.

P á g i n a 6 | 14
Delimitación del método

El campo de la mecánica puede ser subdividido en tres áreas (figura 8.4). La

mecánica teórica estudia las leyes y principios fundamentales de la mecánica, pero
considerando solo su valor científico. La mecánica aplicada transfiere el
conocimiento teórico hacia aplicaciones científicas y de ingeniería, atendiendo
principalmente la construcción de modelos matemáticos que representen fenómenos
físicos. La mecánica computacional resuelve problemas específicos aplicando
métodos numéricos implementados en computadoras digitales (simulación).

Mecánica computacional

Varias ramas de la mecánica computacional pueden distinguirse de acuerdo a la

escala física en que se enfocan (figura 8.5)

La nano mecánica trata con fenómenos de la materia a nivel molecular y atómico.

Por esta razón se encuentra ligada a la física y química de las partículas. La
micromecánica, por su parte, trabaja a nivel cristalográfico y granular de la materia.

P á g i n a 7 | 14
Su principal aplicación tecnológica es el diseño y fabricación de materiales y micro
dispositivos.

La mecánica del continuo estudia los cuerpos a nivel macroscópico, utilizando

modelos continuos en los cuales la microestructura es considerada como
homogénea gracias a promedios cualitativos. Las áreas tradicionales de aplicación
son la mecánica de sólidos y de fluidos. Se incluyen las estructuras debido a que,
obviamente, son fabricadas con sólidos. La mecánica computacional de sólidos usa
aproximaciones de ciencias aplicadas, mientras que la mecánica estructural
computacional hace énfasis en aplicaciones tecnológicas para el análisis y diseño de
estructuras.

La mecánica computacional de fluidos estudia problemas que involucran el

movimiento y equilibrio de líquidos y gases. Las áreas beneficiadas con sus
desarrollos tecnológicos son la hidrodinámica, aerodinámica, acústica, física de la
atmósfera, combustión y propulsión.

Finalmente, los sistemas identifican objetos mecánicos, naturales o artificiales, que

realizan una función fácilmente reconocible. Ejemplos de sistemas artificiales son
todos aquellos dispositivos hechos por el hombre. Los demás sistemas en los que el
hombre no tiene participación se consideran como naturales.

Estática vs. Dinámica

Los problemas del continuo pueden ser subdivididos de acuerdo a si sus efectos
inerciales son considerados o no (figura 8.6). En dinámica, la dependencia del tiempo
actual debe ser considerada explícitamente, porque el cálculo de fuerzas inerciales
requiere tomar derivadas con respecto al tiempo actual.

P á g i n a 8 | 14
Los problemas en estática también pueden depender del tiempo actual, pero ignoran
o desprecian las fuerzas inerciales. Por consiguiente, los problemas estáticos pueden
ser clasificados en estrictamente estáticos y cuasi estáticos. En los primeros, el
tiempo no necesita ser considerado explícitamente; cualquier tiempo registrado como
respuesta, se usa como parámetro de orden, si es que alguno se necesita. En los
problemas cuasi estáticos, tales como la edificación de cimientos, fluencia en metal,
plasticidad o fatiga cíclica, una medida real es requerida pero las fuerzas inerciales
aún son despreciadas.

Lineal vs. no lineal

Una clasificación de los problemas estáticos, particularmente relevante para estas

notas, se muestra en la figura 8.7. El análisis estático lineal involucra problemas
estáticos en los cuales la respuesta es lineal en el sentido causa y efecto. Por
ejemplo, si las fuerzas aplicadas se duplican, los desplazamientos y esfuerzos
internos también se duplican. Los problemas que caen fuera de este dominio son
clasificados como no lineales.

Métodos de discretización

Finalmente, se pueden clasificar los análisis estáticos basándose en el método de

discretización (figura 8.8) por medio del cual el modelo matemático del continuo es
discretizado en el espacio, es decir, convertirlo en un modelo discreto con un número
finito de grados de libertad:

P á g i n a 9 | 14
En los problemas lineales, el método del elemento finito domina actualmente la
escena en lo relativo a discretización espacial. Los métodos por fronteras finitas
constituyen una segunda alternativa en áreas de aplicación específicas. Para
problemas no lineales, el dominio de los métodos por elemento finito es enorme.

Los métodos por diferencias finitas aplicados en sólidos y mecánica de estructuras

han desaparecido virtualmente debido a que son poco prácticos. Sin embargo, para
mecánica de fluidos estos métodos son aún importantes. Los métodos por
volúmenes finitos, que se relacionan directamente con la discretización de las leyes
de conservación, son importantes en problemas difíciles de la mecánica de fluidos;
por ejemplo, en dinámica de gases con número de Reynolds elevado.

Los métodos espectrales se basan en correspondencias que transforman

dimensiones espaciales y/o temporales a espacios (por ejemplo, el dominio de
frecuencias), donde el problema es mucho más sencillo de resolver.

De reciente aparición se tienen los métodos por mallado automático. Estos combinan
técnicas y herramientas del elemento finito, por ejemplo, la formulación variacional e
interpolación, con características de diferencias finitas como el soporte no localizado.

Variantes del método por elemento finito

Actualmente el término MEF identifica a una amplia gama de técnicas (figura 8.9)
que comparten características descritas en los puntos anteriores. Dos
subclasificaciones que se ajustan a la mecánica de estructuras son:

P á g i n a 10 | 14
Utilizando la clasificación anterior, podemos establecer el tema principal de estas
notas: Análisis computacional por elemento finito de problemas estáticos lineales. De
las variantes mostradas, se pondrá especial atención a la formulación por
desplazamiento y a la solución por rigidez. Esta combinación recibe el nombre de
Método directo por rigidez o DSM.

Aplicaciones

El objeto de estas notas es estudiar el Método del elemento finito, pero ¿qué es un
elemento finito? El concepto básico puede ser parcialmente ilustrado a través de un
antiguo problema: encontrar el perímetro 𝐿 de un círculo cuyo diámetro es 𝑑. Como
𝐿 = 𝜋𝑑, esto equivale a obtener un valor numérico para 𝜋. Se dibuja un círculo de
radio 𝑟 y diámetro 𝑑 = 2𝑟 como se muestra en la figura 8.10a. Se inscribe un
polígono regular de 𝑛 lados, donde 𝑛 = 8 en la figura 8.10b. Se renombran los lados
del polígono como elementos y los vértices como nodos. Las etiquetas de los nodos
son enteros que van de 1 a 8. Considérese un elemento típico, el que une los nodos
4-5, como se muestra en la figura 8.10c. Este es un caso del elemento genérico 𝑖 − 𝑗
𝜋
mostrado en la figura 8.10d. La longitud del elemento es 𝐿𝑖𝑗 = 2𝑟𝑠𝑒𝑛 (𝑛) =
180°
2𝑟𝑠𝑒𝑛( ). Como todos los elementos tienen la misma longitud, el perímetro del
𝑛
𝐿𝑛 180°
polígono es 𝐿𝑛 = 𝑛𝐿𝑖𝑗 , de donde la aproximación para 𝜋 resulta 𝜋𝑛 = = 𝑛𝑠𝑒𝑛( ).
𝑑 𝑛

P á g i n a 11 | 14
Para facilitar la explicación se ha implementado un sencillo algoritmo en
Mathematica®, mismo que se anexa en las siguientes páginas (Archimedes_nb).

Algunas ideas del MEF, pueden identificarse gracias al ejemplo anterior. El círculo,
visto como un objeto matemático, es reemplazado por polígonos. Estos constituyen
la aproximación discreta del círculo. Los lados, renombrados como elementos, están
completamente identificados por los nodos en sus extremos. Los elementos pueden
separarse desconectando sus nodos, un proceso llamado desensamble en el MEF.
Gracias a este proceso, un elemento genérico puede ser definido,
independientemente del círculo original, por el segmento que conecta dos nodos 𝑖 y
𝑗. La propiedad relevante del elemento, en este ejemplo, es la longitud de su lado 𝐿𝑖𝑗 ,
misma que puede ser calculada en el elemento genérico independientemente de los
otros, una propiedad conocida como soporte local en el MEF. La propiedad objetivo:
el perímetro del polígono, es obtenido al reconectar n elementos y sumando su
longitud; los pasos correspondientes en el MEF son el ensamble y la solución,
respectivamente. Por supuesto que no existe nada particular en el problema del
círculo, pues la misma técnica se puede utilizar para obtener la longitud de curvas
suaves.

Para ver mas aplicaciones visitar los siguientes links:

http://eprints.uanl.mx/6095/1/1080087879.PDF

http://oa.upm.es/33144/1/SAMARTIN_058.pdf

http://www.simytec.com/docs/pdf_2.pdf

https://cimat.repositorioinstitucional.mx/jspui/bitstream/1008/202/2/TE%20322.pdf

P á g i n a 12 | 14
 Conclusión

Como algo personal a mí me gustaría desarrollarme cono profesor, no

específicamente a un área de la ingeniería, como plan a futuro tengo planeado
estudiar un poco algunas áreas básicas sobre la ingeniería, para ser un profesor que
pueda guiar a sus alumnos en algo que les guste, por lo cual considero bueno
conocer este método para poder enseñarles a las siguientes generaciones de el y las
facilidades que trae este método, y no solo este, si no como muchos métodos
pueden facilitar el camino de la ingeniería.

P á g i n a 13 | 14
 Bibliografía

http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/2548/011
-MPM-Cap8-Final.pdf?sequence=11

http://eprints.uanl.mx/6095/1/1080087879.PDF

http://oa.upm.es/33144/1/SAMARTIN_058.pdf

http://www.simytec.com/docs/pdf_2.pdf

https://cimat.repositorioinstitucional.mx/jspui/bitstream/1008/202/2/TE%20322.pdf

P á g i n a 14 | 14