The semi-explicit nonsmooth Newmark time integrator for robust unilateral contact in dynamic fragmentation simulations

Este artículo introduce y valida un esquema de integración temporal de Newmark semiexplícito no suave (NSN) que maneja de manera robusta el contacto unilateral en simulaciones de fragmentación dinámica mediante la imposición estricta de restricciones, logrando así una estabilidad y precisión superiores sobre los métodos basados en penalización, al tiempo que revela que la disipación por contacto puede aumentar paradójicamente el recuento de fragmentos al mejorar la localización del daño.

Lo esencial

Imagina que estás intentando simular qué sucede cuando un objeto sólido, como un plato de cerámica o una roca, es golpeado con tanta fuerza que se hace añicos en miles de pedazos diminutos. Esto no es solo una ruptura simple; es una explosión caótica donde los fragmentos vuelan, chocan entre sí, rebotan en las paredes y se muelen unos contra otros.

El artículo presenta un nuevo programa informático (un "integrador de tiempo") diseñado para simular este caos sin que la computadora falle o dé resultados absurdos. Así es como funciona, desglosado en conceptos simples:

1. El Problema: La trampa del "Resorte"

Para simular la rotura, los científicos suelen utilizar un método donde fingen que el material está hecho de pequeños resortes. Cuando el material se rompe, los resortes se sueltan. Cuando las piezas chocan entre sí, utilizan "resortes de penalización" para empujarlas y separarlas para que no se atraviesen unas a otras.

La Analogía: Imagina intentar detener una bola de boliche con una banda elástica.

• Si la banda es demasiado floja (baja rigidez), la bola la atravesará (no físico).
• Si la banda es súper tensa (alta rigidez) para detener la bola perfectamente, actúa como una pared rígida. Pero si la haces demasiado tensa, la computadora tiene que dar pasos diminutos, minúsculos, para calcular el rebote, haciendo que la simulación tarde una eternidad.
• La Reclamación del Artículo: La forma antigua (usando estos resortes tensos) es inestable. Provoca que la computadora se desvíe, pierda energía o colapse, especialmente cuando ocurren millones de colisiones al mismo tiempo.

2. La Solución: El "Agente de Tránsito" (Nonsmooth Newmark)

Los autores crearon un nuevo método llamado Nonsmooth Newmark (NSN). En lugar de usar bandas elásticas para empujar las piezas, este método actúa como un estricto agente de tránsito en una intersección concurrida.

La Analogía:

• El Cuerpo (El Coche): El cuerpo principal del objeto se mueve libre y suavemente. La computadora predice hacia dónde iría el coche si no hubiera obstáculos. Esta parte se calcula de forma muy rápida (explícitamente).
• El Contacto (La Intersección): Si el coche choca con una pared o con otro coche, el "agente de tránsito" interviene. En lugar de empujar el coche hacia atrás con un resorte, el agente dice instantáneamente: "¡Alto! No puedes ir ahí". Impone una regla estricta: Prohibido atravesar.
• La Magia: Este método trata la regla de "no atravesar" como una ley dura de la física en lugar de un resorte blando. Permite que la computadora tome pasos de tiempo mucho más grandes porque no tiene que preocuparse de que la banda elástica se tense demasiado.

3. El enfoque de la "Personalidad Dividida"

El artículo describe este método como "semi-explícito". Piensa en esto como una danza de dos pasos:

1. Paso A (La Predicción): La computadora adivina dónde estará todo en el siguiente momento, ignorando las colisiones.
2. Paso B (La Corrección): Si la predicción muestra que dos piezas se superponen, la computadora corrige instantáneamente su velocidad y posición para arreglar la superposición, tal como una bola de billar golpeando a otra y cambiando de dirección instantáneamente.

Esto permite que la simulación sea rápida (como la predicción) pero estable y precisa (como la corrección).

4. Lo que Encontraron (Los Experimentos)

Los autores probaron este nuevo método de "agente de tránsito" contra los antiguos métodos de "banda elástica" usando tres escenarios:

• La Bola que Rebota: Una bola simple rebotando en el suelo. El nuevo método fue tan preciso como los mejores métodos existentes, pero manejó los rebotes sin perder energía o volverse errático.
• La Barra en Impacto: Una barra de metal golpeando una pared. Los métodos antiguos tuvieron dificultades con la velocidad del impacto, pero el nuevo método manejó el "crujido" perfectamente, manteniendo correctos los cálculos de energía.
• La Barra que se Fragmenta: Una barra que ya tiene grietas y luego golpea una pared. Los métodos antiguos requerían pasos de tiempo tan diminutos para mantenerse estables que eran increíblemente lentos. El nuevo método pudo tomar pasos enormes, corriendo 27 veces más rápido siendo más preciso.

5. El Descubrimiento Sorprendente: Fragmentación Confinada

La parte más interesante del artículo involucra un experimento "confinado". Imagina una barra rompiéndose dentro de una caja pequeña en lugar de en un espacio abierto.

• La Vieja Intuición: Podrías pensar que si las piezas rebotan contra las paredes y pierden energía (disipación), habría menos energía restante para romper el material, lo que resultaría en trozos más grandes y menos numerosos.
• El Hallazgo del Artículo: Sucedió lo contrario. Cuando las piezas rebotaban contra las paredes y perdían un poco de energía (disipación por contacto), el material en realidad se rompió en más piezas y más pequeñas.
• ¿Por qué? Los autores explican que el "rebote" actúa como un filtro. En un mundo perfectamente elástico (con mucho rebote), las ondas de tensión rebotan salvajemente, causando que el material se "confunda" y desarrolle muchas grietas diminutas y débiles que no llegan a separarse del todo. Cuando las paredes absorben algo de esa energía, las ondas se calman. Esto permite que la tensión se concentre en puntos específicos, impulsando las grietas a través de todo el material para crear fragmentos limpios y separados.

Resumen

El artículo presenta una nueva herramienta matemática para simular la rotura de objetos tratando las colisiones como reglas duras e instantáneas en lugar de resortes blandos. Esto hace que la simulación por computadora sea:

1. Más Estable: No colapsa ni se desvía.
2. Más Rápida: Puede tomar pasos de tiempo más grandes.
3. Más Precisa: Predice correctamente en cuántas piezas se romperá un objeto.

Los autores concluyen que esta herramienta está lista para ser utilizada en simulaciones 3D complejas, como entender cómo se desprenden los escombros de una nave espacial o cómo se fragmentan las rocas en una avalancha, proporcionando una forma robusta de manejar la danza caótica de millones de fragmentos en colisión.

Resumen técnico

Resumen Técnico: El integrador de tiempo Newmark semi-explícito no suave para el contacto unilateral robusto en simulaciones de fragmentación dinámica

Planteamiento del problema
Las simulaciones numéricas de fragmentación dinámica implican un desafío dual: capturar la rápida evolución topológica de un continuo sólido que se desintegra en fragmentos discretos y resolver los contactos densos y simultáneos que ocurren entre estos fragmentos y las caras de las grietas. Aunque el Método de los Elementos Finitos (FEM) con modelos de zona cohesiva (CZM) extrínsecos gestiona eficazmente la fractura, los tratamientos de contacto convencionales plantean dificultades significativas. Los métodos basados en penalización, frecuentemente utilizados por su compatibilidad con la integración de tiempo explícita, sufren un compromiso fundamental: una baja rigidez permite una interpenetración no física, mientras que una alta rigidez restringe el paso de tiempo estable a niveles computacionalmente prohibitivos. Además, la interacción entre la fractura y el contacto regularizado introduce una falta de suavidad (nonsmoothness) que conduce a una deriva de energía numérica persistente e inestabilidad, enmascarando la física subyacente.

Metodología
Los autores proponen y validan un nuevo esquema de integración de tiempo semi-explícito: el método Newmark no suave (NSN). Este enfoque adapta el marco de la Dinámica de Contacto No Suave (NSCD) a la fragmentación dinámica mediante la desvinculación del tratamiento de la dinámica de volumen y las interacciones de contacto:

1. Formulación semi-explícita: El esquema divide la dinámica en una parte "suave" (deformación de volumen, carga externa y fractura cohesiva) y una parte "no suave" (reacciones de contacto).
   • La dinámica suave se integra explícitamente utilizando un esquema Newmark-$\beta$ de segundo orden ($\beta=0, \gamma=1/2$), manteniendo una alta precisión durante las fases de vuelo libre.
   • La dinámica no suave (contacto unilateral) se trata implícitamente al nivel de la velocidad utilizando un enfoque de Euler implícito de primer orden. Esto impone estrictamente las restricciones de brecha de Signorini sin una complacencia artificial.
2. Implementación algorítmica: El problema de contacto se formula como un Problema de Complementariedad Lineal (LCP), que se resuelve como un problema de Programación Cuadrática (QP) convexa restringido al conjunto de contacto activo. Esto permite soluciones eficientes y paralelizables que escalan con el número de contactos activos en lugar del tamaño total del sistema.
3. Ley de Tracción-Separación (TSL) Modificada: Para abordar el desafío numérico de la rigidez cohesiva infinita a medida que el daño se aproxima a cero en esquemas explícitos, los autores introducen un límite de rigidez ($\tilde{k}$). Esta regularización evita que el paso de tiempo tienda a cero mientras preserva la convexidad del problema de contacto.
4. Validación y Aplicación: El marco se implementa dentro de la biblioteca FEM de código abierto Akantu. Se valida mediante soluciones analíticas y métodos no suaves establecidos (Moreau-Jean y CD-Lagrange) utilizando pruebas de referencia como la bola que rebota y la barra de impacto. Posteriormente, se aplica a la fragmentación dinámica en 1D bajo expansión libre y confinada.

Contribuciones Clave y Resultados

• Precisión y Estabilidad: Las pruebas de referencia demuestran que el esquema NSN logra una precisión comparable a los métodos no suaves establecidos (Moreau-Jean y CD-Lagrange) y supera significativamente a los enfoques basados en penalización. En la prueba de la barra de impacto, el NSN mantiene una precisión de segundo orden para el desplazamiento y una precisión de velocidad superior comparado con CD-Lagrange, particularmente en regímenes elásticos.
• Eficiencia Computacional: Aunque el método NSN incurre en un mayor costo computacional por paso de tiempo (aproximadamente 5 veces más que un esquema de penalización puramente explícito debido a la resolución del QP), su estabilidad robusta permite pasos de tiempo significativamente mayores. Para las pruebas de referencia en 1D, esto resulta en un tiempo total de solución comparable o incluso mejor que los enfoques puramente explícitos, con errores de conservación de energía reducidos en órdenes de magnitud (cercanos a la precisión de máquina).
• Percepciones Físicas en la Fragmentación Confinada:
  • Cambio en el Presupuesto de Energía: En escenarios de expansión confinada (donde los fragmentos colisionan con paredes rígidas), el presupuesto de energía para la fractura se desplaza de la energía cinética local de los fragmentos individuales hacia la energía cinética global de todo el sistema. Esto conduce a una reducción drástica en el tamaño medio de los fragmentos y un aumento en la disipación de la energía total de fractura en comparación con la expansión libre.
  • Papel Contraintuitivo de la Disipación de Contacto: El estudio revela que introducir disipación de contacto (colisiones inelásticas, $e < 1$) reduce la energía total disipada por la fractura, pero aumenta el recuento final de fragmentos. Esto ocurre porque la disipación de contacto actúa como un filtro de paso bajo, suprimi de las oscilaciones de ondas de estrés de alta frecuencia. Este filtrado promueve una propagación de ondas de estrés más limpia y una mejor localización del daño, impulsando las grietas hacia una separación completa en lugar de permitir que el daño permanezca distribuido en zonas parcialmente dañadas ($d < 1$).

Significado
El artículo establece el esquema NSN como una herramienta robusta para generar estadísticas de fragmentación de alta fidelidad en escenarios de contacto denso. Al tratar rigurosamente el contacto como no suave mientras retiene la eficiencia explícita para la dinámica de volumen, el método supera las inestabilidades numéricas inherentes a las formulaciones de penalización. Los autores afirman que este marco es esencial para simular interacciones complejas de múltiples cuerpos, como las encontradas en escenarios de escombros orbitales, donde las frecuentes interacciones fragmento-fragmento impulsan la fragmentación secundaria. El trabajo constituye un paso crítico hacia la extensión del análisis de fragmentación de alta fidelidad a dimensiones superiores, aprovechando la capacidad del método para preservar la fidelidad física sin el costo prohibitivo de las resoluciones implícitas globales.