SIMULACIONES CON ELEMENTOS FINITOS (FEM)
El método de Elementos Finitos (FEM) se ha desarrollado en un método clave, y una tecnología
indispensable en la modelación y simulación de sistemas de ingeniería avanzados en diferentes campos
como transporte, comunicaciones, etc. Los ingenieros y diseñadores analizan procesos sofisticados de
modelación, simulación, diseño, prototipo y por último fabricación. Hay que tener en cuenta que mucho
del trabajo es realizado antes de la fabricación del producto o sistema final. Este es en parte para
asegurar la trabajabilidad del producto final también como la eficiencia en los costos. El proceso es
ilustrado como un flujo en la figura siguiente. Este proceso es muchas veces de naturaleza iterativa, lo
que significa que algunos de los procedimientos son repetidos basado en los resultados obtenidos en la
etapa actual, con el fin de alcanzar un desempeño óptimo al costo más bajo para el sistema que va ha
ser construido. Por ello, las técnicas relacionadas a la modelación y simulación es una forma rápida y
efectiva que resulta en la aplicación del FEM.
El FEM fue usado por primera vez para resolver problemas de análisis de esfuerzos, y ha sido desde
entonces, aplicado a muchos otros problemas como análisis térmico, análisis de flujo de fluidos, análisis
de piezoeléctricos, y muchos otros. Básicamente, el analista, busca determinar la distribución de algúna
variable de campo como el desplazamiento en un análisis de esfuerzos, la temperatura o el flujo de calor
en un análisis térmico, la carga eléctrica en un análisis eléctrico, y así sucesivamente.
El FEM es un método numérico que busca una solución aproximada a la distribución de las variables de
campo en el dominio del problema la cuál es difícil de obtener analíticamente. Se hace mediante la
división del dominio del problema en varios elementos, tales como los presentados en la figura. Las
leyes físicas conocidas son entonces aplicadas a cada pequeño elemento, el cual usualmente tiene una
geometría simple.
Cada elemento está formado por nodos. Las incógnitas son valores discretos de las variables de campo
en los nodos. Luego sigue establecer las ecuaciones en los elementos, después de las cuáles cada
elemento es unido al otro. Este proceso guía a una serie de ecuaciones algebraicas lineales y
simultáneas para el sistema completo que pueden ser resueltas fácilmente para obtener la variable de
campo requerida.
Modelación computacional usando FEM
El comportamiento de un fenómeno en un sistema depende de la geometría o el dominio del sistema,
las propiedades de los materiales o el medio, y las fronteras. Para un sistema de ingeniería, la geometría
o el dominio puede ser muy compleja. Posteriormente, con las condiciones de frontera e iniciales
pueden ser también complicadas. Es por ello que es muy difícil resolver las ecuaciones diferenciales que
gobiernan un problema mediante la analítica. En la práctica, la mayoría de estos problemas se resuelven
mediante métodos numéricos.
El procedimiento de la modelación computacional usando el FEM consiste básicamente en cuatro pasos:
• Modelado de la geometría
• Mallado (discretización)
• Especificación de las propiedades del material
• Especificación de las condiciones iniciales y de frontera
Modelado de la geometría. Estructuras reales, componentes o dominios son muy complejos por lo que
tienen que ser reducidos a geometrías manejables. Las partes curvas de la geometría y sus fronteras
tienen que ser aproximadas a los elementos finitos. La exactitud de la representación de las partes
curvas es controlada por el número de elementos utilizados. Con mayor número de elementos, las
partes curvas se aproximan mejor con elementos que contienen líneas rectas, haciendo las bordes más
pulidos y exactos. Desafortunadamente, mientras más elementos, mayor es el tiempo computacional
necesario para resolver el problema. Entonces, debido a limitaciones de hardware y software
computacional, es necesario limitar el número de elementos, llegando a un número óptimo. Por lo que
por ejemplo, si es necesario obtener información en fronteras curveadas, se hace indispensable
aumentar el número de elementos en dichas regiones.
Referencias:
- J. Kiusalaas, Numerical methods in engineering with Matlab, Cambridge University Press (2005).
- S.C. Chapra, Applied numerical methods with Matlab, Mc Graw Hill (2008).
- C. Woodford y C. Phillips, Numerical methods with worked examples, Springer (1997).