Introducción a elementos finitos

Aproxima campos continuos mediante funciones locales sobre una malla de elementos.

Contexto conceptual

El método de elementos finitos transforma un dominio continuo en una malla de regiones pequeñas llamadas elementos. En cada elemento, el campo físico se aproxima mediante funciones locales sencillas. Los nodos comparten información y permiten construir una solución global que satisface equilibrio, continuidad y condiciones de contorno. La idea central no consiste en copiar exactamente el campo, sino en representar sus variaciones con un espacio aproximante controlado. Al reducir $h$, aparecen más grados de libertad y el $u_h$ puede seguir mejor la solución física. La calidad se evalúa mediante $E_h$ y su ritmo de disminución, descrito por $p$.

🟢 Nivel esencial

Un elemento finito es una pequeña pieza de la geometría donde el campo se representa con una regla simple.

En un elemento lineal, $N_1$ y $N_2$ mezclan $u_1$ y $u_2$ para construir $u_h$. Muchas piezas conectadas forman una malla; cada una aporta una relación local y el ensamblaje enlaza desplazamientos con fuerzas en todo el dominio. Refinar significa disminuir $h$, lo que suele mejorar la representación de curvas y gradientes. La comparación debe conservar las mismas cargas, propiedades y restricciones para que el cambio observado mida realmente el efecto de la discretización. También conviene revisar visualmente que los elementos no estén deformados y que las regiones con variación rápida reciban suficiente resolución. El beneficio, sin embargo, depende de la suavidad de la solución, la calidad de los elementos y las condiciones de contorno. Una malla más densa no arregla automáticamente un modelo físico incorrecto: geometría, propiedades y restricciones deben seguir representando el problema real.

🔵 Nivel formal

Fórmulas básicas

Dentro de un elemento lineal, el campo se interpola con los valores nodales y las funciones de forma:

Interpolación lineal en un elemento

Las funciones asociadas a los dos nodos se definen sobre $\xi$ y forman una partición de unidad:

Funciones de forma lineales unidimensionales

Para una barra axial lineal, $E$, $A$ y $L_e$ determinan la escala de $k_e$:

Rigidez elemental de una barra unidimensional

Preguntas reales del alumno

¿Por qué usar elementos si ya existen ecuaciones diferenciales? Porque muchas geometrías y condiciones no admiten una solución analítica manejable.

¿Los valores nodales son la solución exacta? Son coeficientes del $u_h$. Su precisión debe verificarse.

¿Más elementos siempre mejoran? Solo si la formulación es consistente, la malla tiene calidad y otros errores no dominan.

¿Por qué la matriz puede ser singular? Puede quedar un movimiento rígido libre o existir un error de conectividad.

¿Qué diferencia hay entre refinar h y aumentar p? Refinar h reduce $h$; aumentar p eleva el orden polinómico del elemento.

¿Qué orden debo observar? Depende del elemento, la norma y la suavidad, no solo del grado polinómico.

Síntesis final

Los elementos finitos aproximan campos mediante piezas locales coordinadas globalmente. La solución gana precisión al reducir $h$, pero la convergencia depende de regularidad, calidad de malla, formulación y fronteras. El criterio decisivo es demostrar que $E_h$ disminuye con un $p$ coherente.