Métodos de Elementos Finitos

Los métodos de elementos finitos reformulan una Ecuación Diferencial Parcial (EDP) en forma débil (variacional) y aproximan su solución mediante funciones polinómicas a trozos sobre una malla de elementos simples, lo que permite obtener soluciones precisas en geometrías complejas.

Definition

El método de elementos finitos es una técnica numérica que aproxima la solución de una EDP proyectando su formulación débil sobre un espacio de funciones polinómicas a trozos de dimensión finita, definidas sobre una malla, reduciendo el problema a un sistema de ecuaciones algebraicas.

Scope

Este tema abarca las formulaciones débiles y los entornos de espacios de Sobolev, el método de Galerkin y el lema de Cea, la construcción de espacios de elementos finitos en triangulaciones, el ensamblaje de la matriz de rigidez y el vector de carga, las estimaciones de error a priori y las estimaciones a posteriori que impulsan el refinamiento adaptativo de la malla.

Core questions

• ¿Cómo amplía la formulación débil la clase de soluciones admisibles y sustenta el método?
• ¿Cómo relaciona la proyección de Galerkin, a través del lema de Cea, el error discreto con la mejor aproximación?
• ¿Cómo se construyen los espacios de elementos finitos y se ensambla el sistema global a partir de las contribuciones de los elementos locales?
• ¿Cómo cuantifican la precisión y guían la adaptación de la malla las estimaciones de error a priori y a posteriori?

Key theories

Formulación débil y Lax-Milgram

Multiplicar la EDP por funciones de prueba e integrar la reformula como un problema variacional en un espacio de Sobolev; el teorema de Lax-Milgram garantiza una solución débil única cuando la forma bilineal asociada es acotada y coercitiva, proporcionando la base rigurosa para el método.

Ortogonalidad de Galerkin y lema de Cea

La solución de elementos finitos satisface la ortogonalidad de Galerkin, y el lema de Cea acota su error por una constante multiplicada por el error de la mejor aproximación en el espacio de elementos finitos, reduciendo el análisis de convergencia al poder de aproximación de los elementos elegidos.

Estimación a posteriori y adaptabilidad

Los estimadores de error a posteriori computables acotan el error real utilizando solo la solución discreta y los datos, lo que permite algoritmos adaptativos que refinan la malla donde el error es mayor para lograr una precisión objetivo de manera eficiente.

Mechanisms

El dominio se divide en elementos (triángulos, tetraedros o cuadriláteros), y en cada elemento la solución se representa mediante funciones base polinómicas cuyos soportes se superponen solo en las caras compartidas, lo que da lugar a funciones base globales con soporte local. La sustitución de estas en la forma débil produce un sistema lineal disperso: la matriz de rigidez de la forma bilineal y el vector de carga de los datos, ambos ensamblados elemento por elemento. La resolución del sistema produce los coeficientes de la solución aproximada. Las estimaciones a priori relacionan el error con el tamaño de la malla y el grado polinómico, mientras que los estimadores a posteriori dirigen el refinamiento adaptativo.

Clinical relevance

El método de elementos finitos es la tecnología de simulación dominante en mecánica estructural y de sólidos, transferencia de calor, electromagnetismo y biomecánica, y se utiliza ampliamente en dinámica de fluidos; su capacidad para manejar geometrías complejas, propiedades de materiales variadas y refinamiento adaptativo lo convierte en la columna vertebral de la mayoría del software comercial de análisis de ingeniería.

History

El método surgió de la ingeniería estructural en la década de 1950 y se le dio una base matemática variacional basándose en el trabajo anterior de Courant; la teoría de aproximación rigurosa fue desarrollada por Ciarlet, Babuska y otros en la década de 1970, convirtiendo el método de elementos finitos tanto en una herramienta práctica como en un área profunda del análisis numérico.

Key figures

• Richard Courant
• Olgierd Zienkiewicz
• Philippe Ciarlet
• Susanne C. Brenner

Seminal works

• brenner2008
• ern2004

Frequently asked questions

¿Por qué reformular la EDP en forma débil?

La forma débil reduce la diferenciabilidad requerida de la solución y plantea el problema en un entorno de espacio de Hilbert donde la existencia, unicidad y aproximación pueden analizarse rigurosamente, y se adapta naturalmente a las aproximaciones polinómicas a trozos en mallas complejas.

¿Qué hace que los elementos finitos sean buenos para geometrías complejas?

Debido a que el dominio se divide en elementos pequeños y de forma simple que pueden dimensionarse y orientarse para ajustarse al contorno, las mallas de elementos finitos se adaptan a formas intrincadas mucho más fácilmente que las cuadrículas regulares que requieren los métodos de diferencias finitas.