Efficient strain-space hyperreduction in large-deformation solid mechanics

Resumen Técnico: Hiperreducción Eficiente en el Espacio de Deformaciones para la Mecánica de Sólidos de Grandes Deformaciones

Planteamiento del Problema
Las técnicas de Reducción de Orden de Modelo (MOR) que operan en el espacio de deformaciones han demostrado recientemente un rendimiento superior en la homogeneización computacional, particularmente para Elementos de Volumen Representativo (RVE). Sin embargo, su aplicación se ha limitado en gran medida a este contexto específico. Extender estas técnicas a problemas generales de mecánica de sólidos de grandes deformaciones presenta desafíos significativos. A diferencia de las formulaciones basadas en desplazamientos, los métodos en el espacio de deformaciones no pueden imponer directamente condiciones de contorno (BC) de Dirichlet, ya que las deformaciones de los puntos de integración en la frontera no se conocen explícitamente, y los valores de contorno se prescriben típicamente como desplazamientos, no como deformaciones. Además, los enfoques estándar de espacio de deformaciones suelen depender de condiciones de contorno periódicas, las cuales no son aplicables a problemas de macroescala con valores de contorno arbitrarios.

Metodología
Los autores proponen un marco generalizado para aplicar MOR en el espacio de deformaciones a problemas con condiciones de contorno de Dirichlet parametrizadas y arbitrarias. La metodología central consta de tres componentes principales:

1. Manejo de Contornos mediante Levantamiento (Lifting): Para satisfacer las condiciones de contorno de Dirichlet no nulas, los autores emplean una estrategia de levantamiento. Se computa un campo de referencia consistente con el contorno offline para valores de desplazamiento unitarios en cada frontera independiente. El campo solución real se descompone en este campo de referencia (escalado por los valores de contorno reales) y un "campo de fluctuación". El campo de fluctuación se anula en todos los contornos de Dirichlet por construcción, lo que permite que la base reducida (derivada de la Descomposición de Órdenes Propios, POD) satisfaga las condiciones de contorno a priori. Este enfoque transforma las condiciones de contorno inhomogéneas en homogéneas para el resolvedor reducido.

2. Formulaciones en el Espacio de Deformaciones: El artículo generaliza tres técnicas específicas de hiperreducción a este nuevo entorno:

   • ECM (Método de Cubatura Empírica) en el Espacio de Deformaciones: Reformulado para realizar la integración reducida directamente en los puntos de Gauss, evitando operaciones a nivel de elemento.
   • E3C (Cubatura de Clúster Corregida Empíricamente) en el Espacio de Deformaciones: Generalizada para problemas mecánicos arbitrarios. E3C optimiza los puntos de integración y sus pesos en el espacio de deformaciones, permitiendo que estos sean distintos de los puntos de Gauss de la malla original.
   • EMSL (Muestreo de Material y Linealización Empírica) en el Espacio de Deformaciones: Extendida más allá de los RVE. EMSL linealiza el comportamiento del material alrededor de las deformaciones promedio por clúster. A diferencia de otros métodos que integran aproximadamente una ley de material exacta, EMSL integra exactamente una aproximación lineal por clúster, lo que resulta en un problema afín por paso de carga que no requiere un resolvedor iterativo de Newton-Raphson durante la fase online.

3. Línea Base de Comparación: Estos métodos de espacio de deformaciones se contrastan con el método estándar de Muestreo y Ponderación de Conservación de Energía (ECSW) en el espacio de desplazamientos, el cual sirve como una línea base robusta y preservadora de la estructura.

Contribuciones Clave

• Generalización de la MOR en el Espacio de Deformaciones: El trabajo extiende con éxito las técnicas de MOR en el espacio de deformaciones (ECM, E3C, EMSL) desde la homogeneización computacional hacia la mecánica de sólidos de grandes deformaciones con condiciones de contorno de Dirichlet arbitrarias.
• Tratamiento de Contorno Novedoso: La introducción de una formulación de campo de fluctuación basada en levantamiento permite la imposición de condiciones de contorno arbitrarias en el espacio de deformaciones sin requerir datos de deformación en la frontera.
• Benchmarking Exhaustivo: Los autores proporcionan una comparación rigurosa de los métodos de espacio de deformaciones frente a la reducción de orden en el espacio de desplazamientos (ECSW) utilizando dos ejemplos hiperelásticos de estimación de parámetros inversos, con parámetros de material variables y trayectorias de carga complejas (tensión, cortante y mixta).

Resultados
Los métodos fueron probados en dos problemas de placas hiperelásticas con agujeros elípticos:

1. Ejemplo 1 (Deformación Moderada): Una placa bajo carga de tracción con hasta un 20% de deformación.
2. Ejemplo 2 (Grandes Deformaciones y Carga Compleja): La misma geometría sometida a tensión (hasta 50%), cortante (hasta 25%) y carga mixta, con una malla más fina y mayor costo computacional.

Hallazgos de Rendimiento:

• Precisión vs. Tiempo de Ejecución: En ambos ejemplos, los métodos de espacio de deformaciones superaron consistentemente al método ECSW en el espacio de desplazamientos en la relación entre precisión y tiempo de ejecución.
• Aceleraciones (Speedups):
  • E3C: Logró aceleraciones de aproximadamente 10,000 veces (Ejemplo 1) y 30,000 veces (Ejemplo 2) en comparación con el modelo de orden completo, manteniendo una alta precisión (por ejemplo, <1% de error). E3C generalmente proporcionó la mayor precisión para un presupuesto computacional dado entre los métodos de espacio de deformaciones cuando se utilizaron suficientes puntos de integración.
  • EMSL: Demostró una eficiencia de tiempo de ejecución excepcional, logrando aceleraciones de hasta 100,000 veces (Ejemplo 1) y 300,000 veces (Ejemplo 2). EMSL destacó cuando los presupuestos online y offline eran muy limitados, superando a menudo a E3C en velocidad, aunque a veces con una precisión ligeramente menor en regímenes altamente no lineales en comparación con E3C con muchos puntos de integración.
  • ECM: Funcionó mejor que ECSW, pero generalmente se quedó rezagado respecto a E3C y EMSL en las métricas específicas probadas.
• Costos Offline: Aunque E3C tiene una fase de entrenamiento offline más costosa debido a la optimización iterativa, los autores señalan que esto es manejable y significativamente menor que el tiempo requerido para generar los datos de entrenamiento. EMSL posee el costo offline más bajo.

Significancia y Reivindicaciones
El artículo sostiene que la metodología propuesta justifica la investigación adicional en las técnicas de MOR en el espacio de deformaciones más allá de la homogeneización computacional. Los autores destacan que estas técnicas ofrecen:

• Eficiencia Superior: Reducciones drásticas en el tiempo de ejecución online (órdenes de magnitud) manteniendo una alta precisión.
• No Intrusividad: Los métodos requieren únicamente datos de desplazamiento, deformación y esfuerzo (fácilmente extraídos de resolvedores de caja negra) y una rutina de material. Evitan las operaciones a nivel de elemento y no requieren conocimiento de los detalles de la discretización (por ejemplo, tipos de elementos) o de los vectores de residuo.
• Aplicabilidad: El manejo exitoso de condiciones de contorno de Dirichlet arbitrarias elimina una barrera importante para la aplicación de estas técnicas eficientes en problemas estructurales de macroescala.

Los autores concluyen que, si bien persisten desafíos para problemas altamente no lineales (como la elastoplasticidad) y la estimación de deformaciones por clúster en EMSL, el rendimiento demostrado sugiere que estos métodos podrían habilitar nuevas aplicaciones de ingeniería en áreas como el control, la optimización y el ajuste de modelos de materiales, donde el costo computacional es actualmente prohibitivo.