Autor: Augusto Rodríguez, Héctor Sanzana, Dasset Gallardo, Francisco Barrios. / Centro de Ingeniería Mecánica y Diseño Industrial (CIMECDI). FIIDT.
El MEF fue al principio desarrollado en 1943 por Richard Courant, quien utilizó el método matemático de Ritz de análisis numérico y minimización de las variables de cálculo para obtener soluciones aproximadas a un sistema de vibración. Poco después, un documento publicado en 1956 por M. J. Turner, R. W. Clough, H. C. Martin, y L. J. Topp estableció una definición más amplia del análisis numérico. El documento se centró en “la rigidez y deformación de estructuras complejas.
Con la llegada de los primeros ordenadores se instaura el cálculo matricial de estructuras. Este parte de la discretización 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 como resultado de aplicar las ecuaciones de equilibrio a los nodos de la estructura. Este sistema de ecuaciones se esquematiza de la siguiente manera: F = Kd donde las incógnitas d son los desplazamientos en los nodos que se hallan a partir de las “fuerzas” o “solicitaciones” en los nodos (F) y de la rigidez de las barras (matriz de rigidez K). Conocidos dichos desplazamientos es posible determinar los esfuerzos en las barras. La solución obtenida es exacta.
Cuando se produce la llegada de los primeros equipos de cálculo en la década de los años 50, el cálculo de estructuras se encontraba en un punto en el que los métodos de cálculo predominantes consistían en método iterativos (métodos de Cross y Kani) que se realizaban de manera manual y, por tanto, resultaban bastante tediosos. El cálculo de una estructura de edificación de varios pisos, por ejemplo, podía llevar varias semanas, lo cual suponía un coste sustancial de tiempo en detrimento de la posibilidad de invertir éste en la optimización de la estructura.
La llegada de los ordenadores permitió el resurgimiento del método de los desplazamientos ya conocidos en siglos anteriores (Navier, Lagrange, Cauchy), pero que eran difíciles de aplicar dado que al final conducían a la resolución de enormes sistemas de ecuaciones inabordables desde el punto de vista manual.
Cuando las aplicaciones prácticas de elementos fi nitos crecieron en tamaño, los requerimientos de tiempo de cálculo y memoria de los ordenadores aumentaron. En ese punto el desarrollo de algoritmos más eficientes se volvió fundamental. Para la resolución de los sistemas de ecuaciones se potenció el estudio de la adaptabilidad de los algoritmos ya conocidos (Gauss, Cholesky, Crout, Gradiente conjugado, etc.). El ahorro de tiempo es fundamental y con ello el uso del método matricial se extiende.
Sin embargo, y a pesar de desarrollarse modelizaciones de elementos superficiales mediante barras (losas con emparrillados, elementos curvos mediante aproximaciones de elementos rectos, etc.), se plantean grandes dificultades ante estructuras continuas (superficies y volúmenes) y con geometrías complejas. De ahí que sea precisamente dentro del campo aeroespacial donde comiencen a desarrollarse las nuevas técnicas del MEF. Dada su generalidad el método se amplió a otros campos no estructurales como la conducción de calor, la mecánica de fluidos, etc. donde compitió con otros métodos numéricos como el de método de las diferencias finitas que, aun siendo más intuitivos, tenían de nuevo dificultades de planteamiento para geometrías complejas.
Con la llegada de los centros de cálculo y los primeros programas comerciales en los años 60, el MEF a la vez que se populariza en la industria mecánica refuerza sus bases teóricas en los centros universitarios.
En los años 70 se produce un gran crecimiento de la bibliografía científica (el primer libro sobre elementos finitos propiamente dicho fue publicado en 1967 por Zienkiewicz y Cheung) así como su extensión del método a otros problemas como los no lineales. En esta década, el MEF estaba limitado a caros ordenadores centrales generalmente en propiedad de las industrias aeronáuticas, de automoción, de defensa y nucleares. Se estudian nuevos tipos de elementos y se sientan las bases matemáticas rigurosas del método, que había aparecido antes más como técnica de la ingeniería que como método numérico de la matemática.
Por último, a partir de la década de los 80, con la generalización de los ordenadores personales, se extiende el uso de los programas comerciales que se especializan en los diversos campos, instaurándose el uso de pre y postprocesadores gráficos que realizan el mallado y la representación gráfica de los resultados. Se continúa en el estudio de la aplicación del método a nuevos modelos de comportamiento (plasticidad, fractura, daño continuo, etc.) y en el análisis de los errores.
En la actualidad, dentro del campo estructural, el MEF comparte protagonismo con el método matricial, siendo muchos los programas que mezclan el análisis por ambos métodos, debido sobre todo a la mayor necesidad de memoria que requiere el análisis por elementos fi nitos frente a un cálculo de barras
MEF (Método de Elemento Finito)
La simulación por el método de los elementos finitos es un recurso muy extendido en el ámbito de la ingeniería y consiste en el análisis numérico de un determinado sistema físico (análisis estructural, dinámica de fluidos, transferencia de calo entre otros), donde es necesario definir un dominio o espacio de interés y unas condiciones de contorno.
En el diseño mecánico existen algunas variables cuyo análisis se hace imprescindible para poder lanzar un producto al mercado o incluso para poder empezar su construcción de una pieza equipo u otro. Me refiero al estudio de esfuerzos y de formaciones que se presentan en una pieza o ensamblaje sujeto a cargas y restricciones. La obtención de estos resultados nos dará la información necesaria para saber si el material que utilizamos es el adecuado, para validar la seguridad de nuestro diseño o para considerar reducir la cantidad de material y seguir trabajando con un buen factor de seguridad, además de predecir las características del sistema bajo condiciones de servicio. La ingeniería mecánica, plenamente consolidada en fundamentos físicos y matemáticos, ofrece métodos para resolver estas variables en un análisis estático y dinámico del modelo que es objeto de estudio, sin embargo la resolución de problemas a “lápiz y papel” demandan demasiado tiempo, y posiblemente limitación en exactitud y optimización de resultados. El método de elemento finito, es una discretización de una pieza en una cantidad finita de elementos para los cuales se resuelven simultáneamente ecuaciones diferenciales que describen el desplazamiento de los nodos, o puntos de intersección de los elementos. Dada la imposibilidad técnica para dar solución a una inmensa matriz de ecuaciones diferenciales simultaneas, este método se clasifica como un método numérico, ya que es necesario el empleo de una computadora con excelente capacidad de procesamiento, de tal manera que el análisis de un sistema complejo ameritaría un buen servidor para que las posibilidades de un óptimo análisis con características muy aproximadas a la realidad pueden ser realizadas en tiempo récord. No obstante, para que los resultados nos proporcionen información de comportamientos reales, es necesario de un muy buen modelo de partida y con características y condiciones de frontera correctamente establecidas.
El método de los elementos finitos permite obtener una solución numérica aproximada sobre un cuerpo, estructura o dominio (cuando es un medio continuo) sobre el que están definidas ciertas ecuaciones diferenciales en forma débil o integral que caracterizan el comportamiento físico del problema , dividiéndolo en un número elevado de subdominios no-intersectantes entre sí denominados (elementos finitos).
El conjunto de elementos finitos forma una subdivisión del dominio también denominada discretización. Dentro de cada elemento se distinguen una serie de puntos representativos llamados (nodos). Dos nodos son adyacentes si pertenecen al mismo elemento finito; además, 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 cálculos se realizan sobre una malla de puntos (llamados nodos), que sirven a su vez de base para discretización del dominio en elementos finitos. La generación de la malla se realiza usualmente con programas especiales llamados generadores de mallas, en una etapa previa a los cálculos que se denomina pre-proceso.
Generalmente, la modelización de un sistema complejo requiere trabajo y experiencia, y su simulación puede tardar horas o incluso días dependiendo de la potencia computacional a disposición y del nivel de complejidad del sistema. También se utiliza en casos donde la experimentación empírica es difícil, costosa o peligrosa.
Seguidamente toda la información del modelo se la transmiten entre los diferentes elementos mediante los nodos. Es decir, si los nodos no son comunes o no están relacionados de alguna manera, aunque ocupen la misma posición no habrá transferencia de información.
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. 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 y al número de grados de libertad de cada nodo.
En definitiva, con esta metodología en lugar de obtener la solución exacta a una variable que cumple un sistema de ecuaciones diferenciales, se obtiene el valor aproximado de dicha variable en unos puntos determinados (nodos) estimándose posteriormente el resto de las variables asociadas.
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 transmisión de calor, de mecánica de fluidos para calcular campos de velocidades y presiones (mecánica de fluidos computacional, CFD) 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.
Las aplicaciones del método son muy amplias, como ejemplo en: Mecanismos especiales, sistemas mecánicos, procesos de manufactura, transferencia de calor y energía, la aeroespacial, la robótica, biomedicina, cardiología, la mecánica del deporte, etcétera. Pueden presentar las oportunidades de aplicación según se trate de comprender algún comportamiento.
En resumen, las etapas a seguir podemos agruparlo en los siguientes pasos.
Pre-procesador.
Donde se prepara el modela para el cálculo, en él se realizas las operaciones de:
- Dibujo del modelo o importación si se ha generado por un CAD
- Selección de elementos a emplear
- Selección de materiales a emplear
- Asignación de Elementos y propiedades de los materiales
- Mallado de los componentes del modelo
- Aplicación de las cargas exteriores (puntuales , lineales o de superficie )
- Aplicación de las condiciones de contorno
Calculador
Es la parte del programa donde se realiza todo el cálculo por el MEF y genera las soluciones, para ello se tiene:
- Selección de tipo de cálculo a realizar, por ejemplo si es un tipo de cálculo transitorio en régimen armónico, estático, etc.
- Configuración de los parámetros de cálculo, selección de intervalos de tiempo, números de interacciones entre otros.
- Inicio del cálculo el programa empieza transfiriendo las cargas al modelo, genera las matrices de rigidez, realiza la triangulación de la matriz, resuelve el sistema de ecuaciones y genera la solución.
Post procesador
Es la herramienta que permite la representación gráfica, así como resultados indirectos que se pueden obtener utilizando las soluciones del modelo. Lo cual el profesional en el área respectiva analizara dichos resultados y tomara las decisiones correspondientes.
Bibliografía
- J. ABAFFY AND E. SPEDICATO, ABS Projection Methods, Halstead Press, 1989.
- L. M. ADAMS AND H. JORDAN, Is SOR color-blind?. SIAM Journal on Scientific and Statistical Computing.
- BUDYNAS, R.G. Advanced Strength and Applied Stress Analysis McGraw-Hill, New York, 1977.Google Scholar
- FONG, H.H. Interactive graphics and commercial finite element codes. Mechanical Engineering 106 (6) June 1984.
- TURNER, M.J.; CLOUGH, R.W.; Martin, H.C.; TOPP, L.J. Stiffness and deflection analysis of complex structures. J. Aero. Sci. 23, 1956.
- https://ingeoexpert.com/articulo/fundamento-y-origen-del-metodo-de-los-elementos-finitos-mef/?v=a99877f71bd9
- TOCHER, James L., SAE Transactions Vol. 89, Section 1: 800004–800252 (1980), pp. 24-28 Published by: SAE International, https://www.jstor.org/stable/44632404
Contacto: orlandoaco2@gmail.com
