Deformation gradient averaging regularization for third medium contact

Título: Regularización por promediado del gradiente de deformación para el contacto de tercer medio

1. Planteamiento del Problema

El método de contacto de tercer medio ha ganado popularidad como alternativa a los métodos de contacto tradicionales (como penalización o multiplicadores de Lagrange), especialmente en contextos donde la búsqueda de fronteras de contacto es compleja, como en la optimización topológica. Este método introduce un material ficticio en el espacio vacío entre cuerpos que es extremadamente compliant (baja rigidez) antes del contacto, pero se endurece drásticamente al ser comprimido a volúmenes cercanos a cero, imponiendo así las restricciones de contacto.

Sin embargo, en configuraciones de grandes deformaciones, el material de tercer medio tiende a deformarse excesivamente en la fase pre-contacto, lo que requiere técnicas de regularización para estabilizar el comportamiento. Los modelos existentes suelen basarse en la penalización de los segundos gradientes de desplazamiento (calculados localmente). Esto presenta dos desventajas principales:

1. Requiere el uso de elementos finitos de segundo orden (cuadráticos) o la introducción de grados de libertad adicionales (DOFs) en elementos de primer orden, lo que complica la implementación.
2. Penaliza modos de deformación más allá de los estrictamente necesarios, lo que puede ser restrictivo.

2. Metodología Propuesta

Los autores proponen un enfoque alternativo que evita el cálculo explícito de segundos gradientes de desplazamiento. La metodología se basa en tres componentes clave dentro de la densidad de energía de deformación del tercer medio:

1. Término de Contacto Volumétrico: Se utiliza la parte volumétrica de la ley de material hiperelástico de Neo-Hookean. Esta parte reacciona asintóticamente a la compresión extrema ($J \to 0$), impidiendo la penetración, pero se omite la parte isocórica para evitar interferencias no deseadas en la fase pre-contacto.
2. Término de Rigidez Constante (Lineal Elástica): Se introduce un comportamiento lineal elástico en la fase pre-contacto para mejorar la respuesta del material antes del contacto. Esto mantiene una rigidez constante (definida por un módulo de Young escalado respecto al material base) en lugar de depender únicamente de la respuesta hiperelástica, lo que mejora la estabilidad numérica.
3. Regularización por Promediado de Deformación (Contribución Central):
   • En lugar de penalizar los gradientes de segundo orden, se penalizan las desviaciones del gradiente de deformación ($F$) en los puntos de integración de Gauss respecto a su valor promedio en el elemento (calculado en el centroide del elemento, $\bar{F}$).
   • La energía de regularización se define como: $W_{\bar{F}} = \frac{1}{2} \kappa_{\bar{F}} \| F - \bar{F} \|^2$.
   • Este término controla la uniformidad del gradiente de deformación dentro del elemento. Si el gradiente varía mucho entre los puntos de integración, la energía penaliza esta no uniformidad.

Ventaja Metodológica: Esta formulación permite utilizar elementos finitos de primer orden (lineales) estándar sin necesidad de grados de libertad adicionales. Solo requiere elementos con múltiples puntos de integración (como cuadriláteros bilineales con integración de Lobatto o Gauss) para que exista variación en $F$ dentro del elemento.

3. Contribuciones Clave

• Simplificación de la Implementación: Se elimina la necesidad de elementos de segundo orden o DOFs adicionales, haciendo que el método sea fácil de implementar en códigos de elementos finitos existentes.
• Eliminación de Segundos Gradientes: Se reemplaza la penalización de segundos gradientes (que requiere continuidad $C^1$ o DOFs extra) por una penalización de la no uniformidad del gradiente de deformación dentro del elemento.
• Robustez en Grandes Deformaciones: La combinación de la rigidez constante pre-contacto y la regularización por promediado estabiliza las simulaciones de contacto en grandes deformaciones y optimización topológica.
• Conservación de Propiedades Energéticas: Al derivarse de un potencial de energía, el método preserva la simetría del tensor tangente de rigidez global, lo cual es crucial para la convergencia de los solucionadores Newton-Raphson.

4. Resultados Numéricos

El método fue implementado en el software de elementos finitos de código abierto OOFEM y validado mediante varios benchmarks:

• Benchmarks de C y Caja Cerrada: Se demostró que los resultados coinciden muy bien con los métodos tradicionales de contacto "nodo-a-segmento" (NTS). El método de tercer medio propuesto logra una precisión comparable en la detección del contacto y la respuesta fuerza-desplazamiento, incluso con mallas de diferentes refinamientos.
• Pandeo Restringido y Autocontacto: En un problema de inestabilidad (una viga comprimida que pandea y toca paredes), el método capturó correctamente las bifurcaciones de carga, el contacto con las paredes y el autocontacto posterior. La transferencia de fuerzas a través del tercer medio fue precisa.
• Fuerzas Configuracionales: Se validó el cálculo de fuerzas configuracionales en un problema de contacto con un indentador rígido. Los resultados numéricos coincidieron excelentemente con la solución analítica cerrada, demostrando que el método maneja correctamente la interacción normal y tangencial en grandes deformaciones y geometrías curvas.
• Estudio de Malla: Se observó que el parámetro de penalización $\kappa_{\bar{F}}$ debe escalarse con el tamaño de la malla para mantener la consistencia de los resultados.

5. Significado e Impacto

Este trabajo representa un avance significativo en la mecánica computacional de contacto, particularmente para la optimización topológica.

• Facilidad de Adopción: Al no requerir cambios estructurales en la formulación de elementos finitos (como DOFs extra o elementos de orden superior), los investigadores pueden integrar el método de tercer medio en sus códigos existentes con modificaciones mínimas.
• Versatilidad: El enfoque es universal y aplicable tanto en 2D como en 3D, y es adecuado para simulaciones multifísicas donde se necesita mapear espacios tipo "aire" o vacío.
• Estabilidad: La propuesta resuelve problemas de inestabilidad numérica comunes en simulaciones de contacto con materiales muy blandos, ofreciendo una alternativa robusta a los métodos de contacto tradicionales que a menudo fallan en topologías cambiantes.

En resumen, los autores presentan una formulación elegante y computacionalmente eficiente que democratiza el uso del método de tercer medio en problemas complejos de mecánica de sólidos no lineales.