Dynamic brittle fracture using Lip-field approach in an explicit dynamics context

Resumen Técnico: Fractura Frágil Dinámica mediante el Enfoque de Campo Lip en un Contexto de Dinámica Explícita

Planteamiento del Problema
Este artículo aborda la simulación numérica de la fractura frágil dinámica en materiales sometidos a cargas de alta tasa de deformación, como impactos o explosiones. Mientras que la mecánica de fractura clásica (por ejemplo, el modelo de Griffith) presenta dificultades con la iniciación y ramificación de grietas, se han desarrollado enfoques modernos como los Modelos de Zona Cohesiva (CZM) y los modelos de daño difuso. Sin embargo, los modelos de daño difuso locales sufren de una sensibilidad de malla patológica, lo que requiere una regularización para introducir una escala de longitud y asegurar resultados independientes de la malla. Las técnicas de regularización existentes, como los métodos de Campo de Fase (Phase-Field) y de Conjunto de Nivel Grueso (Thick Level Set - TLS), suelen implicar funcionales de energía no locales que pueden ser computacionalmente costosos. El desafío específico abordado aquí es la extensión del enfoque de campo Lip —un método variacional que introduce escalas de longitud mediante restricciones de Lipschitz en el campo de daño— desde sus aplicaciones previas en cuasiestática y dinámica unidimensional hacia la fractura dinámica bidimensional dentro de un marco de integración temporal explícita.

Metodología
Los autores proponen una formulación variacional para la mecánica de daño dinámica basada en la minimización de un potencial incremental. El núcleo de la metodología consiste en:

1. Formulación Variacional: El problema mecánico se define mediante la minimización de un potencial incremental $F$ que incluye la energía cinética, la energía libre y el potencial de disipación. Las variables de estado son el desplazamiento ($u$) y una variable escalar de daño ($d$). El modelo constitutivo utiliza materiales de estado estándar generalizado con funciones específicas de degradación de energía ($g(d)$) y ablandamiento ($h(d)$).
2. Regularización de Campo Lip: A diferencia de los modelos de Campo de Fase que modifican el funcional de energía con un término de gradiente, el enfoque de campo Lip mantiene el funcional de energía local (no regularizado). En su lugar, impone restricciones de Lipschitz al campo de daño, requiriendo que el gradiente del daño esté acotado por un valor crítico ($1/l$). Esto restringe el campo de daño a un espacio de Lipschitz específico $\mathcal{L}$, introduciendo efectivamente una escala de longitud $l$ sin alterar la densidad de energía local.
3. Esquema Escalonado Explícito: El artículo emplea un esquema de integración temporal de Newmark explícito (método de diferencias centrales). Se utiliza un enfoque escalonado donde los campos de desplazamiento y daño se resuelven secuencialmente. Para un paso de tiempo dado, el desplazamiento se actualiza primero, seguido de la actualización del daño.
4. Solucionador de Daño Eficiente: Una contribución algorítmica clave es el uso de la estimación de límites (derivada de trabajos previos) para optimizar el cálculo del daño. El proceso involucra:
   • Resolver un problema de minimización puramente local para encontrar un campo de daño local ($d_{loc}$).
   • Proyectar $d_{loc}$ en límites superiores e inferiores basados en la restricción de Lipschitz.
   • Identificar las regiones donde estos límites difieren.
   • Resolver el problema de minimización no local (con restricción de Lipschitz) únicamente en la pequeña porción del dominio donde los límites no son iguales. En las regiones donde los límites coinciden, se acepta la solución local. Esto reduce significativamente el costo computacional en comparación con resolver el problema no local en todo el dominio.
5. Discretización: La discretización espacial utiliza elementos triangulares lineales para el campo de desplazamiento. Se utiliza una "malla Lip" separada (construida a partir de los centroides de los elementos) para almacenar y restringir el campo de daño, asegurando que el número de restricciones de Lipschitz se minimice.

Contribuciones Clave

• Extensión a la Dinámica 2D: El artículo extiende con éxito el enfoque de campo Lip, validado previamente en fragmentación dinámica y casos cuasiestáticos en 1D, hacia problemas de fractura dinámica bidimensional.
• Integración de Dinámica Explícita: Describe una implementación específica del enfoque de campo Lip utilizando un esquema de integración temporal explícita, lo cual es computacionalmente eficiente para problemas de impacto de alta velocidad.
• Eficiencia Computacional mediante Límites: Los autores demuestran que el uso del estimador de límites permite que la costosa minimización no local se restrinja a subdominios pequeños, haciendo que el enfoque de campo Lip sea computacionalmente competitivo para simulaciones dinámicas.
• Distinción Variacional: El trabajo clarifica la distinción entre los enfoques de campo Lip y de Campo de Fase: el campo Lip mantiene un funcional de energía local con restricciones no locales, mientras que el Campo de Fase utiliza un funcional de energía no local.

Resultados de la Simulación
El modelo propuesto fue validado mediante ejemplos numéricos bidimensionales:

1. Tensión de Borde Simple (SENT): El modelo simuló con éxito la propagación y ramificación de grietas dinámicas. Los resultados mostraron convergencia de malla en la energía potencial y disipada. El fenómeno de ramificación de grietas se capturó sin criterios ad-hoc, y las velocidades de la punta de la grieta generalmente se mantuvieron por debajo de la velocidad de la onda Rayleigh, excepto durante picos transitorios durante la ramificación. Los resultados fueron comparables con las simulaciones de Campo de Fase y TLS encontradas en la literatura.
2. Influencia de la Intensidad de la Carga: Al variar la tracción aplicada, el modelo demostró la capacidad de capturar múltiples ramas de grietas, ocurriendo la ramificación más cerca de la muesca bajo cargas más altas.
3. Prueba de Kalthoff-Winkler: El modelo simuló la iniciación y propagación de grietas bajo carga de corte dinámica.
   • A bajas velocidades de impacto ($33$ m/s), predijo un ángulo de propagación de grieta Modo-I de $\approx 62.4^\circ$, consistente con las observaciones experimentales de $\approx 70^\circ$.
   • A altas velocidades de impacto ($100$ m/s), capturó la ramificación sucesiva de grietas y múltiples grietas.
   • Limitación: Los autores señalan que, si bien los experimentos muestran una transición hacia la falla dúctil (fractura por corte a $\approx -10^\circ$), el modelo actual de daño elástico no captura este comportamiento dúctil, ya que carece de plasticidad.

Significancia y Reivindicaciones
El artículo afirma que el enfoque de campo Lip ofrece un marco variacional robusto para la fractura dinámica que preserva la relevancia matemática del problema mientras introduce las escalas de longitud necesarias. Su principal significancia radica en su eficiencia computacional lograda a través del solucionador basado en límites, que evita resolver ecuaciones no locales globalmente. Los autores aseguran que el enfoque maneja eficazmente fenómenos de fractura complejos como la iniciación, ramificación y coalescencia de grietas sin criterios ad-hoc.

El artículo concluye modestamente que, aunque el enfoque de campo Lip es efectivo para la fractura frágil y la fragmentación, está limitado al comportamiento de daño elástico. Reconoce que simular la falla dúctil a altas tasas de deformación requeriría una extensión a modelos de daño elástico-plástico, lo cual está fuera del alcance del presente trabajo. El método se presenta como una alternativa viable a los métodos de Campo de Fase y TLS, particularmente donde la eficiencia computacional en la dinámica explícita es crítica.