A locking free mixed FEM based on a pure pseudostress based formulation for the elasticity eigenproblem

Resumen Técnico: Un Método de Elementos Finitos Mixto Libre de Bloqueo Basado en una Formulación de Pseudoesfuerzo Puro para el Problema de Autovalores de Elasticidad

Planteamiento del Problema
El artículo aborda la aproximación numérica del problema de autovalores de elasticidad lineal, que gobierna la vibración de estructuras elásticas. El principal desafío que se aborda es el fenómeno de "bloqueo" (locking) que ocurre en los métodos numéricos estándar cuando el coeficiente de Poisson ($\nu$) se aproxima a $1/2$ (el régimen casi incompresible o perfectamente incompresible). En este régimen, el coeficiente de Lamé $\lambda$ tiende al infinito, lo que provoca que las formulaciones estándar fallen o pierdan precisión. Además, los autores pretenden desarrollar una formulación que evite la imposición de la simetría del esfuerzo (ya sea débil o fuerte), un requisito común en los métodos de elasticidad mixta que complica la implementación y el análisis.

Metodología
Los autores proponen un nuevo método de elementos finitos mixtos (FEM) basado exclusivamente en el tensor de pseudoesfuerzo ($\rho$). El campo de desplazamiento ($u$) no es una incógnita independiente en el sistema primario, sino que se recupera mediante un paso de post-procesamiento.

1. Formulación de Pseudoesfuerzo Puro:
   El problema continuo se reformula para depender únicamente del tensor de pseudoesfuerzo $\rho$ y del desplazamiento $u$. Al utilizar la identidad que relaciona $\rho$ con $\nabla u$, el sistema se reduce a un sistema de ecuaciones diferenciales parciales de pseudoesfuerzo puro. La formulación variacional se define en el espacio $X = H(\text{div}, \Omega)$, específicamente en un subespacio $X_0$ donde la traza del tensor integra cero.

2. Análisis Matemático:

   • Teoría de Operadores No Compactos: A diferencia de las formulaciones de elasticidad estándar donde el operador de solución es compacto, el operador de solución $T^\lambda$ asociado con esta formulación de pseudoesfuerzo es no compacto. En consecuencia, los autores analizan la convergencia y las estimaciones de error utilizando la teoría de operadores no compactos (referenciando los trabajos [7, 8]).
   • Caracterización Espectral: El espectro del operador de solución se descompone en tres partes: un autovalor de multiplicidad infinita en 1 (asociado a tensores divergencia-cero), un autovalor en 0, y una secuencia de autovalores de multiplicidad finita que convergen a 0.
   • Análisis de Límite: El artículo demuestra rigurosamente que, a medida que $\lambda \to \infty$ (es decir, $\nu \to 1/2$), el espectro del problema casi incompresible converge al espectro del caso perfectamente incompresible, el cual coincide con el sistema de autovalores de Stokes. Esta convergencia se establece sin bloqueo numérico.

3. Discretización:
   El método emplea elementos finitos de Raviart–Thomas (RT), los cuales son naturalmente adecuados para los espacios de conformidad $H(\text{div})$. Los espacios discretos se construyen para satisfacer la condición de traza nula requerida para $X_0$.

4. Análisis de Error:

   • Estimaciones A Priori: Los autores derivan estimaciones de error de orden óptimo tanto para los autovalores como para los autofunciones. Notablemente, el error del autovalor converge con orden $O(h^{2\min\{r,k\}})$, donde $r$ es la regularidad de la autofunción y $k$ es el grado polinómico.
   • Estimaciones A Posteriori: Se introduce un estimador de error a posteriori basado en residuos. Los autores demuestran su fiabilidad (límite superior) y su eficiencia (límite inferior) tanto para los regímenes casi incompresibles como para los perfectamente incompressibles. El estimador está diseñado para ser robusto con respecto al coeficiente de Poisson.

Contribuciones Clave

• Formulación Libre de Bloqueo: El método propuesto evita el bloqueo numérico cuando $\nu \to 1/2$, manteniendo tasas de convergencia óptimas tanto en regímenes casi como perfectamente incompresibles.
• Sin Imposición de Simetría: La formulación no requiere la imposición explícita de la simetría del esfuerzo, lo que simplifica la implementación de los elementos finitos.
• Enfoque de Pseudoesfuerzo Puro: El trabajo extiende el concepto de pseudoesfuerzo puro (previamente aplicado a problemas de Stokes) al problema de autovalores de elasticidad, permitiendo la recuperación del desplazamiento y el esfuerzo mediante post-procesamiento.
• Análisis de Operadores No Compactos: El artículo aplica con éxito la teoría de operadores no compactos al problema de autovalores de elasticidad, un paso necesario debido a la estructura específica de la formulación de pseudoesfuerzo.
• Estimador A Posteriori Robusto: Se propone un indicador de error novedoso que permanece eficiente y fiable incluso en el límite de la incompresibilidad, una característica que no se encuentra comúnmente en los estimadores existentes para esta formulación específica.

Resultados Numéricos
Los hallazgos teóricos se validan mediante experimentos numéricos en 2D y 3D utilizando el entorno FEniCSx.

• Convergencia: Las pruebas en un dominio cuadrado y un cubo unitario demuestran que el método alcanza una convergencia de segundo orden ($O(h^2)$) para los autovalores a través de varios coeficientes de Poisson, incluyendo $\nu = 0.5$.
• Refinamiento Adaptativo: Los experimentos en dominios no convexos (un sector 2D con una esquina entrante y un cubo con muescas) muestran que la estrategia adaptativa impulsada por el estimador propuesto logra recuperar tasas de convergencia óptimas ($O(\text{dof}^{-1}$ en 2D y $O(\text{dof}^{-2/3}$ en 3D) a pesar de la reducida regularidad de las autofunciones cerca de las singularidades.
• Robustez: Los índices de eficiencia permanecen acotados y estables a medida que $\nu$ se aproxima a $0.5$, confirmando la naturaleza libre de bloqueo del método y la robustez del estimador de error.

Significado
El artículo afirma que este trabajo proporciona un marco numérico robusto y libre de bloqueo para analizar la vibración de estructuras elásticas, particularmente en regímenes donde los materiales son casi o perfectamente incompresibles. Al eliminar la necesidad de imponer la simetría y utilizar una formulación de pseudoesfuerzo puro, el método ofrece una alternativa simplificada pero matemáticamente rigurosa a los métodos mixtos tradicionales. La capacidad de capturar con precisión el espectro en el límite incompresible sin autovalores espurios o pérdida del orden de convergencia representa un avance significativo en la mecánica computacional de los problemas de autovalores de elasticidad.