An Investigation of Stabilization Scaling in Finite-Strain Virtual Element Methods for Hyperelasticity

Este trabajo propone un método de estabilización libre de submalla para elementos virtuales en hiperelasticidad de gran deformación que desacopla los canales desviador y volumétrico, eliminando la sensibilidad a la teselación interna y garantizando una robustez uniforme en el régimen casi incompresible mediante un escalado espectral consistente con el tensor tangente de Newton.

Lo esencial

Imagina que estás construyendo un puente gigante, pero en lugar de usar ladrillos cuadrados perfectos, usas piedras de formas extrañas y desiguales (triángulos, pentágonos, hexágonos). Quieres que el puente sea tan fuerte y flexible como la realidad, capaz de soportar terremotos o cargas pesadas sin romperse.

En el mundo de la ingeniería computacional, esto se hace con un método llamado Método de Elemento Virtual (VEM). Es como un "super-ingeniero" que puede calcular cómo se deforman esas piedras extrañas bajo presión.

Sin embargo, hay un problema. Cuando las piedras se deforman mucho (como cuando estiras una goma elástica hasta el límite), el cálculo principal del ingeniero a veces se pierde en detalles que no puede ver. Es como si el ingeniero solo mirara la forma general de la piedra, pero ignorara las pequeñas vibraciones o "temblores" internos que ocurren en las esquinas. Si no controlamos esos temblores, el puente podría parecer más rígido de lo que es en realidad, o incluso colapsar en la simulación.

Aquí es donde entra el estabilizador. Es como un "seguro" o un "amortiguador" que añadimos para controlar esos temblores invisibles.

El Problema: El "Amortiguador" que se vuelve un "Bloqueo"

En los métodos actuales, este seguro se calcula de una manera complicada:

1. Dividen la piedra en triángulos pequeños (una sub-triangulación) para hacer los cálculos. Esto es como si tuvieras que dibujar líneas imaginarias dentro de cada piedra solo para medir su estabilidad.
2. Usan una "receta" de materiales que mezcla dos cosas: la resistencia a ser estirado (cizalladura) y la resistencia a ser comprimido (volumen).

El problema es que, cuando el material es casi incompresible (como el agua o la goma muy dura, donde no puedes reducir su volumen), esta receta se descontrola. El seguro empieza a aplicar una fuerza de compresión gigante a los temblores que deberían ser solo de estiramiento.

La analogía: Imagina que tienes un resorte que debe controlar solo el movimiento lateral de una puerta. Pero, por error, el resorte también se conecta a la cerradura. Cuando intentas abrir la puerta (cambiar el volumen), el resorte se pone tan tenso que la puerta se traba y no se mueve. En ingeniería, esto se llama "bloqueo" (locking). El puente parece más rígido de lo que es y no se deforma correctamente.

La Solución: Un Seguro "Desacoplado" y "Sin Subdivisiones"

Los autores de este paper (Paulo y Rodrigo) proponen una nueva forma de diseñar este seguro. Su idea es genialmente simple y elegante:

1. Sin subdivisiones (Submesh-free): En lugar de dibujar triángulos imaginarios dentro de la piedra, miran directamente a los vértices (las esquinas) de la piedra. Es como medir la estabilidad de un edificio solo revisando cómo se mueven sus pilares principales, sin necesidad de inspeccionar cada ladrillo interior.
2. Desacoplamiento (Decoupling): Separan completamente el seguro en dos canales independientes:
   • Canal de Cizalladura (Estiramiento): Controla los temblores laterales. Se escala solo con la "dureza al estirar" del material.
   • Canal Volumétrico (Compresión): Controla los cambios de volumen. Se escala con la "dureza a comprimir".

La analogía: Imagina que tienes dos amortiguadores separados para tu coche. Uno controla las baches laterales (cizalladura) y el otro controla si el coche se hunde o sube (volumen). En los métodos viejos, si el coche se hundía mucho (casi incompresible), el amortiguador de baches se activaba por error y hacía que el coche se moviera muy duro. En el nuevo método, si el coche se hunde, el amortiguador de baches no se activa. Solo se activa el de volumen. Así, el coche sigue manejando suavemente en los baches, sin importar cuánto se hunda.

¿Por qué es importante esto?

El paper demuestra matemáticamente y con pruebas de computadora que:

• Es más robusto: Funciona bien incluso si las piedras tienen formas muy extrañas o deformadas.
• Evita el bloqueo: En materiales casi incompresibles (como goma o tejidos biológicos), el nuevo método permite que la simulación se deforme de manera realista, mientras que los métodos antiguos se quedaban "trabados" y daban resultados erróneos.
• Es más rápido y limpio: Al no necesitar subdivisiones internas y tener fórmulas más directas, es computacionalmente más eficiente.

En resumen

Este trabajo es como inventar un nuevo tipo de seguro de coche que sabe exactamente cuándo actuar. Si el coche se mueve de lado, el seguro lateral actúa. Si el coche se aplasta, el seguro vertical actúa. Pero si el coche se aplasta, el seguro lateral no se pone nervioso.

Gracias a esto, los ingenieros pueden simular materiales muy complejos (como el corazón humano, neumáticos o gomas industriales) con formas irregulares, sabiendo que sus cálculos no están "trabados" por un error en la forma en que calculan la estabilidad. Es un avance que hace que las simulaciones sean más precisas, rápidas y confiables.

Resumen técnico

Resumen Técnico: Escalamiento de Estabilización en Métodos de Elemento Virtual para Hiperelasticidad de Gran Deformación

1. Planteamiento del Problema

El Método de Elemento Virtual (VEM) es un marco de Galerkin conforme que permite el uso de mallas poligonales y poliedrales generales. En su versión de orden bajo (k=1), las contribuciones consistentes de la energía se calculan mediante proyecciones polinómicas, mientras que los modos no resueltos (el núcleo del proyección, $\ker(\Pi^\nabla_E)$) deben controlarse mediante un término de estabilización.

En la práctica actual de VEM para hiperelasticidad de gran deformación, los esquemas de estabilización existentes presentan dos deficiencias críticas, especialmente en el régimen casi incompresible ($\nu \to 1/2$):

1. Dependencia de Sub-mallas: Suelen evaluar energías de sustitución (surrogate) sobre una triangulación interna auxiliar arbitraria, lo que introduce sensibilidad a la calidad y orientación de dicha sub-malla.
2. Contaminación Volumétrica en Modos Isocóricos: Las escalas de estabilización a menudo utilizan parámetros de Lamé modificados o factores de compresibilidad que inyectan "proxies" volumétricos en las penalizaciones de tipo cortante. Esto provoca que los modos isocóricos (que no cambian de volumen) en el núcleo del proyección reciban una rigidez artificialmente alta dependiente del módulo volumétrico ($\lambda$), lo que conduce a un bloqueo (locking) severo y a una respuesta deformada incorrecta en materiales casi incompresibles.

El objetivo es diseñar una estabilización que sea libre de sub-mallas, desacople los canales desviador y volumétrico, y escale consistentemente con la energía tangente de Newton actual, sin degradarse cuando $\nu \to 1/2$.

2. Metodología

Los autores proponen un nuevo esquema de estabilización basado en los siguientes pilares metodológicos:

• Enfoque Espectral: Se reformula el requisito de estabilidad del VEM como una equivalencia espectral entre la forma bilineal de estabilización y la energía tangente de Newton restringida al núcleo del proyección ($\ker(\Pi^\nabla_E)$). Esto se caracteriza mediante cocientes de Rayleigh generalizados y límites de autovalores, proporcionando una métrica independiente de la base para evaluar cómo se asigna la energía a los modos no resueltos.
• Estabilización Solo-Núcleo (Kernel-Only): Se propone un término de estabilización que depende exclusivamente del residuo del proyección en los grados de libertad (nodos). Esto elimina la necesidad de triangulaciones internas auxiliares.
• Desacoplamiento Desviador/Volumétrico:
  • Término Desviador: Se escala únicamente por una medida de corte ($\mu_E$) y se enriquece con pesos direccionales basados en la geometría (ejes principales de inercia del elemento) para manejar anisotropías. No incluye dependencias de $\lambda$.
  • Término Volumétrico: Se escala por una medida volumétrica independiente ($\kappa_E$) y se aplica como una penalización en los bordes sobre las componentes normales del residuo. Este término puede ser suprimido o acotado cuando $\nu \to 1/2$ para evitar rigidez residual en modos isocóricos.
• Análisis Teórico: Se demuestran teoremas que garantizan que la estabilización desviadora propuesta es equivalente a $\mu_E |\cdot|_{1,E}^2$ en el núcleo, con constantes uniformes independientes del tamaño de la malla ($h_E$) y del coeficiente de Poisson ($\nu$).

3. Contribuciones Clave

1. Marco Espectral para VEM No Lineal: Se establece un marco teórico donde la estabilidad se define mediante la equivalencia espectral en el núcleo del proyección, ofreciendo una herramienta diagnóstica precisa para evaluar la asignación de energía a modos de deformación no resueltos.
2. Nuevo Esquema de Estabilización Desacoplado: Se introduce una formulación libre de sub-mallas que separa estrictamente los canales desviador y volumétrico. La parte desviadora escala solo con el módulo de corte y utiliza pesos geométricos acotados para elementos anisotrópicos.
3. Resultados de Escalamiento Uniforme en $\nu$: Se demuestra teóricamente que la estabilización desviadora propuesta mantiene constantes de equivalencia independientes de $\nu$, evitando el sobrecalentamiento de modos isocóricos en el límite incompresible.
4. Diagnósticos a Nivel de Elemento: Se desarrollan pruebas específicas (modo isocórico de un solo elemento y análisis modal restringido al núcleo) para detectar la contaminación volumétrica en las estabilizaciones clásicas.
5. Validación Numérica Robusta: Se demuestra mediante el benchmark de la membrana de Cook que el nuevo método elimina el bloqueo volumétrico observado en métodos clásicos, manteniendo una convergencia coherente en familias de mallas distorsionadas y poligonales irregulares.

4. Resultados

Los experimentos numéricos validan la superioridad del método propuesto frente a las formulaciones clásicas (basadas en energías de sustitución y sub-triangulaciones):

• Diagnóstico de Modo Isocórico: En un elemento unitario, al aumentar $\nu$ hacia 0.5, la estabilización clásica asigna una energía creciente a modos isocóricos (debido a la dependencia de $\lambda$), mientras que la propuesta permanece constante y escalada solo por $\mu$.
• Análisis Modal: El espectro de autovalores restringido al núcleo confirma que la nueva estabilización:
  • Anula exactamente los modos polinómicos.
  • Redistribuye la rigidez en elementos anisotrópicos según los pesos direccionales diseñados.
  • Mantiene una separación clara entre modos desviadores (constantes al variar $\kappa_E$) y volumétricos.
• Membrana de Cook (Gran Deformación):
  • Regímen Casi Incompresible ($\nu = 0.499$): La formulación clásica muestra un bloqueo severo, subestimando drásticamente el desplazamiento en mallas gruesas y convergiendo lentamente.
  • Propuesta: La formulación desacoplada produce una respuesta libre de bloqueo, con una convergencia rápida y consistente hacia un valor de referencia en todas las familias de mallas (cuadriláteros regulares, distorsionados y teselaciones de Voronoi), independientemente de la distorsión geométrica.

5. Significado e Impacto

Este trabajo aborda una limitación fundamental en la aplicación del VEM a problemas de mecánica de sólidos no lineales con materiales casi incompresibles.

• Robustez Geométrica: Al eliminar la dependencia de triangulaciones internas y utilizar pesos direccionales adaptativos, el método es robusto frente a mallas poligonales distorsionadas o anisotrópicas.
• Eliminación del Bloqueo Volumétrico: La separación explícita de los canales desviador y volumétrico, junto con la eliminación de proxies volumétricos en la escala de corte, resuelve el problema de bloqueo que ha limitado la precisión del VEM en régimen incompresible.
• Eficiencia Computacional: Al ser libre de sub-mallas y tener derivadas analíticas cerradas (tangentes consistentes), el método es computacionalmente eficiente y compatible con métodos de Newton-Raphson estándar.
• Fundamento Teórico: Proporciona una justificación teórica rigurosa basada en la equivalencia espectral, guiando el diseño futuro de estabilizaciones que deben ser "conscientes de los modos" (mode-aware) y físicamente consistentes.

En conclusión, la propuesta representa un avance significativo hacia la implementación robusta y precisa del VEM en problemas de hiperelasticidad complejos, ofreciendo una alternativa superior a las estrategias de estabilización basadas en sub-mallas y escalas mixtas tradicionales.