A comprehensive anisotropic phase-field fracture formulation with a novel positive/negative projection scheme

Resumen Técnico: Una Formulación Integral de Fractura de Campo de Fase Anisotrópico con un Novedoso Esquema de Proyección Positivo/Negativo

Planteamiento del Problema
La teoría de campo de fase (PF, por sus siglas en inglés) ha surgido como un enfoque eficiente para simular patrones de grietas complejos sin trayectorias predefinidas. Sin embargo, extender estos modelos a materiales frágiles anisotrópicos presenta desafíos significativos. Mientras que los modelos de PF isotrópicos suelen depender de descomposiciones de la energía de deformación (como las divisiones espectrales o las divisiones volumétrico-desviatóricas) para evitar el crecimiento de grietas no físico bajo compresión, estos métodos frecuentemente fallan en mantener la consistencia termodinámica o la ortogonalidad de la energía cuando se aplican a la elasticidad anisotrópica. En consecuencia, existe una falta de un marco rigurosamente consistente que maneje simultáneamente la elasticidad anisotrópica, asegure el equilibrio termodinámico y suprima el daño bajo estados de compresión sin recurrir a aproximaciones heurísticas.

Metodología
Los autores desarrollaron un marco de campo de fase integral para materiales frágiles anisotrópicos, implementado dentro del software comercial de elementos finitos Abaqus mediante subrutinas de elemento de usuario (UEL) y de material de usuario (UMAT). El núcleo de la metodología involucra tres componentes principales:

1. Marco Variacionalmente Consistente: La formulación se adhiere a las leyes fundamentales de equilibrio y a la primera y segunda leyes de la termodinámica. Utiliza una estrategia de solución monolítica para resolver simultáneamente las ecuaciones acopladas de desplazamiento y de campo de fase, asegurando una convergencia robusta.
2. Elasticidad Anisotrópica: El modelo incorpora la simetría de material trigonal (representativa de cristales tipo cuarzo) mediante la rotación del tensor de rigidez desde el sistema de coordenadas del cristal hacia el marco global de elementos finitos. Esto permite capturar la respuesta mecánica y la rigidez dependientes de la orientación.
3. Novedoso Esquema de Proyección Positivo/Negativo Anisotrópico (PNP): Para abordar la no coaxialidad de la tensión y la deformación en materiales anisotrópicos, los autores proponen un novedoso esquema PNP basado en la norma de la energía complementaria.
   • Estrategia de Mapeo: El método mapea los tensores de tensión y deformación anisotrópicos reales a un espacio isotrópico auxiliar y ficticio utilizando tensores de transformación derivados de la equivalencia de energía.
   • Proyección: En este espacio isotrópico auxiliar, donde la tensión y la deformación son coaxes, se resuelve un Problema de Complementariedad Lineal (LCP) para descomponer la deformación en partes positiva (tracción) y negativa (compresión). Esto asegura que la descomposición sea ortogonal en la norma de energía y termodinámicamente consistente.
   • Retransformación: Las cantidades proyectadas se mapean de vuelta al espacio anisotrópico real. La función de degradación del campo de fase se aplica únicamente a la contribución de tracción (positiva), mientras que la contribución de compresión (negativa) permanece intacta, suprimiendo eficazmente el crecimiento de grietas bajo compresión.

Contribuciones Clave

• División Anisotrópica Termodinámicamente Consistente: El artículo introduce una extensión rigurosa del esquema PNP de Wu y Cervera a la elasticidad anisotrópica. Al utilizar un mapeo de equivalencia de energía a un espacio isotrópico auxiliar, la formulación garantiza que la energía que impulsa la grieta esté asociada consistentemente con la deformación por tracción, incluso cuando la tensión y la deformación no son coaxes.
• Implementación Monolítica en Abaqus: Los autores proporcionan una implementación práctica utilizando un UEL para el sistema acoplado de desplazamiento-daño y un UMAT para la visualización del post-procesamiento. Este enfoque aprovecha la robustez del solver de Abaqus permitiendo al mismo tiempo un comportamiento constitutivo anisotrópico complejo.
• Verificación Exhaustiva: El marco se verifica mediante una jerarquía de pruebas, incluyendo pruebas de un solo elemento bajo tracción, compresión y cortante; pruebas de parche con mallas distorsionadas para asegurar la robustez de la formulación del elemento; y múltiples hitos de cuasi-estática de fractura (tensión con entalla de un solo borde, compresión, cortante y flexión en tres puntos).

Resultos
Los resultados numéricos demuestran la eficacia de la formulación propuesta:

• Benchmarks Isotrópicos: El modelo reproduce soluciones analíticas y resultados de la literatura para estándares benchmarks isotrópicos (por ejemplo, tensión con entalla de un solo borde y flexión en tres puntos) con alta precisión. Suprime con éxito el daño bajo compresión pura y exhibe respuestas simétricas para cargas de cortante positivas y negativas.
• Comportamiento Anisotrópico: Para materiales trigonales (cristales de arena de sílice), el modelo captura respuestas dependientes de la orientación distintivas.
  • Rigidez y Resistencia: Las simulaciones muestran cargas pico más altas y respuestas más rígidas cuando la carga se alinea con la dirección "fuerte" cristalográfica (eje c) en comparación con la dirección "débil" (eje a).
  • Trayectorias de Grieta: El modelo predice trayectorias de grietas dependientes de la orientación. En la orientación fuerte, las grietas se desvían de la dirección de la entalla inicial debido a términos de acoplamiento anisotrópico (por ejemplo, $C_{14}$), curvándose hacia una dirección energéticamente preferida antes de reorientarse cerca de los bordes libres. En la orientación débil, las grietas se propagan de manera más lineal.
  • Compresión y Cortante: La división basada en LCP logra prevenir la evolución del daño bajo carga de compresión para ambas orientaciones. Bajo cortante, el modelo captura la asimetría en la respuesta entre cargas de cortante positivas y negativas para la orientación fuerte, una consecuencia directa de la simetría trigonal y el acoplamiento constitutivo, manteniendo la simetría para la orientación débil.
• Estudios Paramétricos: El estudio confirma que el refinamiento de la malla, la longitud de escala de regularización y la tasa de liberación de energía crítica ($G_c$) influyen en la carga pico y el comportamiento de ablandamiento post-pico, de acuerdo con las expectativas de la teoría de campo de fase.

Significancia
El artículo afirma que este trabajo proporciona un "marco integral y novedoso" para cuantificar el comportamiento de fractura de materiales frágiles anisotrópicos. Su principal significancia radica en cerrar la brecha entre la consistencia termodinámica rigurosa y el modelado práctico de la fractura dependiente de la dirección. Al extender el esquema PNP a la elasticidad anisotrópica mediante un mapeo de equivalencia de energía, la formulación supera las limitaciones de las divisiones heurísticas previas que asumían la coaxialidad tensión-deformación. Los autores sostienen que el modelo es capaz de representar fenómenos de fractura en sólidos frágiles con dependencia direccional, capturando no solo la rigidez y las cargas pico dependientes de la orientación, sino también trayectorias de grietas complejas impulsadas por la anisotropía. El trabajo se presenta como un paso hacia simulaciones más precisas de materiales como el cuarzo y otras estructuras cristalinas, aunque los autores señalan las limitaciones actuales respecto a la simetría trigonal, geometrías 2D y la energía de fractura independiente de la orientación, sugiriendo estos como objetivos para futuras extensiones.