Resumen Técnico: SoliDualSPHysics

Planteamiento del Problema

La Hidrodinámica de Partículas Suavizadas (SPH) se ha consolidado como un método sin malla (meshfree) poderoso para la dinámica de fluidos y la mecánica de sólidos de grandes deformaciones debido a su naturaleza lagrangiana. Sin embargo, el panorama actual de los marcos de trabajo SPH de alto rendimiento y código abierto está fuertemente sesgado hacia la dinámica de fluidos computacional (CFD) y la interacción fluido-estructura (FSI) que involucra cuerpos rígidos o elasticidad de deformación pequeña. Existe una brecha significativa en los marcos unificados capaces de manejar la mecánica de sólidos deformables autónoma con comportamientos constitutivos avanzados, específicamente la plasticidad de deformación finita y la fractura frágil físicamente consistente.

Las formulaciones de fractura SPH existentes suelen depender de modelos de daño empíricos ad hoc, enfoques de partículas de grieta o modelos de zona cohesiva, los cuales sufren de dependencia de la malla, resultados no convergentes o la necesidad de un seguimiento de grietas explícito o enriquecimiento local. Además, integrar estas capacidades avanzadas de mecánica de sólidos en arquitecturas escalables CPU/GPU adecuadas para simulaciones a gran escala sigue siendo un desafío. La necesidad específica identificada es una formulación SPH unificada que soporte hiperelasticidad, plasticidad de deformación finita y fractura frágil sin seguimiento de grietas explícito, manteniendo la eficiencia computacional mediante la paralelización.

Metodología

Los autores presentan SoliDualSPHysics, una extensión de código abierto y acelerada por GPU del marco de trabajo DualSPHysics. La metodología se basa en una formulación de SPH de Lagrangiano Total (TLSPH), la cual permite la aplicación directa de cargas externas y condiciones de contorno en la configuración no deformada.

Marco Teórico

1. Balance de Momento: El solver emplea una discretización TLSPH de forma fuerte para la ecuación de balance de momento. Utiliza una corrección de gradiente de kernel de primer orden para restaurar la consistencia e incluye un término de viscosidad artificial para mitigar las inestabilidades numéricas y los modos de energía cero espurios.
2. Modelado Constitutivo:
   • Hiperelasticidad: El marco soporta los modelos de St. Venant–Kirchhoff y neo–Hooke compresible. Estos se acoplan con una formulación de campo de fase para simular la fractura frágil.
   • Plasticidad de Deformación Finita: Se implementa una extensión a la plasticidad $J_2$ independiente de la tasa con endurecimiento isotrópico. Esta utiliza una descomposición multiplicativa del gradiente de deformación ($F = F^e F^p$) y un algoritmo de retorno de mapa exponencial para asegurar la objetividad de marco y la consistencia termodinámica bajo grandes rotaciones.
   • Nota de Acoplamiento: En la implementación actual, la plasticidad y la fractura de campo de fase no están acopladas; el modelo de plasticidad describe la respuesta elastoplástica sin fractura, mientras que el modelo de campo de fase gestiona la fractura frágil en materiales hiperelásticos.
3. Fractura Frágil (Campo de Fase): La fractura se modela mediante una formulación de campo de fase hiperbólica. A diferencia de las formulaciones elípticas que requieren resolver sistemas lineales globales (incompatibles con el paso de tiempo explícito), este enfoque utiliza una ecuación de evolución hiperbólica para la variable de daño $s$. Esto permite que la iniciación, propagación, ramificación y coalescencia de grietas emerjan naturalmente sin seguimiento o enriquecimiento explícito. El modelo impone la irreversibilidad mediante un funcional de historia y restringe la propagación de grietas a regiones de tracción mediante la descomposición de la energía de deformación elástica.
4. Integración Temporal: El sistema se avanza utilizando esquemas de integración temporal explícitos (Verlet o Simpléctico/Leapfrog), con pasos de tiempo restringidos por condiciones de tipo CFL basadas en la propagación de ondas y la aceleración.

Arquitectura de Software e Implementación

• Plataforma: El código está escrito en C++ y CUDA, aprovechando la arquitectura paralela nativa de CPU (OpenMP) y GPU (CUDA) de DualSPHysics.
• Entrada/Salida: La configuración se gestiona mediante archivos de entrada XML. La extensión introduce nuevos bloques (<deformstrucs>, <mathexpressions>) para definir cuerpos deformables, propiedades de materiales, condiciones de contorno y expresiones definidas por el usuario.
• Expresiones Definidas por el Usuario: Una característica clave es la inclusión de un analizador de expresiones matemáticas que permite a los usuarios prescribir cantidades dependientes del tiempo y el espacio (por ejemplo, campos de velocidad, cargas de tracción) utilizando variables integradas (posición, tiempo, desplazamiento) y operadores matemáticos/lógicos estándar.
• Condiciones de Contorno: Se pueden aplicar condiciones de contorno de velocidad y tracción flexibles de forma global o a grupos de superficie específicos definidos por geometrías auxiliares.

Contribuciones Clave

1. Marco Unificado: SoliDualSPHysics proporciona un entorno único de código abierto para simular hiperelasticidad, plasticidad de deformación finita y fractura frágil dentro del paradigma SPH.
2. Fractura de Campo de Fase en SPH: Implementa una formulación de campo de fase hiperbólica para la fractura frágil que evita el costo computacional de los solvers elípticos, permitiendo una dinámica explícita eficiente para la ramificación y coalescencia de grietas sin criterios ad hoc.
3. Plasticidad de Deformación Finita: La integración de un modelo de plasticidad $J_2$ multiplicativa con un mapeo de retorno de mapa exponencial extiende la aplicabilidad del código a problemas plásticos de gran deformación.
4. Escalabilidad: Al utilizar la infraestructura de CPU/GPU de DualSPHysics, el software logra una aceleración sustancial para simulaciones a gran escala que involucran millones de partículas.
5. Flexibilidad: La inclusión de un analizador de expresiones definido por el usuario y definiciones de condiciones de contorno flexibles mejora la adaptabilidad del código para escenarios de carga complejos y aplicaciones de investigación.
6. Disponibilidad de Código Abierto: El código fuente completo, la documentación y los casos de prueba se ponen a disposición de forma gratuita a través de un repositorio público en GitHub para garantizar la reproducibilidad.

Resultados y Validación

El software fue verificado y validado contra una serie de problemas numéricos de referencia y datos experimentales:

• Elasto-Dinámica: La oscilación libre de vigas y placas (2D y 3D) demostró un excelente acuerdo con las soluciones analíticas y los resultados de referencia del Método de Elementos Finitos (FEM). Las simulaciones mostraron un intercambio de energía casi conservativo y una disipación numérica mínima.
• Gran Deformación: Una viga en voladizo 3D sometida a una gran tracción superficial distribuida confirmó la capacidad de manejar la no linealidad geométrica y condiciones de contorno complejas sin oscilaciones espurias.
• Dinámica Torsional: Un caso de columna en torsión 3D validó la capacidad del marco para resolver la propagación de ondas torsionales y la recuperación elástica, con resultados que coinciden con estudios de referencia.
• Fractura: Un benchmark de ramificación de grieta dinámica (placa con muesca previa bajo tensión) demostró la capacidad del modelo para capturar patrones de fractura complejos, incluyendo la ramificación, de manera consistente con la literatura establecida de campo de fase.
• Rendimiento: Se demostró que la implementación en GPU ofrece una aceleración significativa en comparación con las ejecuciones solo en CPU, soportando simulaciones escalables.

Significado y Alcance

El artículo posiciona a SoliDualSPHysics como un paso significativo para cerrar la brecha entre los solvers de fluidos SPH de alto rendimiento y las capacidades avanzadas de la mecánica de sólidos. Su principal significancia radica en proporcionar una herramienta unificada, de código abierto y acelerada por GPU que elimina la necesidad de seguimiento de grietas ad hoc en simulaciones de fractura frágil y soporta la plasticidad de deformación finita dentro de un marco sin malla.

Los autores señalan explícitamente las limitaciones de la versión actual:

• Está restringida a la respuesta de material isotrópica (hiperelasticidad isotrópica, plasticidad $J_2$ isotrópica).
• La plasticidad y la fractura no están acopladas; no se incluyen mecanismos de fractura dúctil.
• Los modelos constitutivos anisotrópicos, los materiales compuestos y los operadores dedicados de control de hourglass aún no se han implementado.

A pesar de estas limitaciones, este trabajo establece una base sólida para futuros desarrollos en mecánica de sólidos y modelado de fractura dentro de la comunidad SPH, ofreciendo una plataforma reproducible y extensible para investigadores.