A multi-phase-field model for fiber-reinforced composite laminates based on puck failure theory

Este artículo propone un modelo multifase basado en la teoría de falla de Puck, que utiliza variables de campo de fase independientes y un método de superposición de mallas para predecir con precisión los modos de daño en láminas de composites reforzados con fibra bajo diversas condiciones de carga.

Lo esencial

Imagina que los materiales compuestos, como los que se usan en las alas de un avión o en los coches de carreras, son como las capas de una lasaña gigante. Cada capa (o "placa") tiene fibras que la hacen fuerte en una dirección específica, como espaguetis alineados. Cuando empujas o tiras de esta lasaña, a veces se rompe de formas muy complicadas: las fibras se rompen, la "salsa" (la resina) que las une se agrieta, o las capas se separan entre sí.

Predecir exactamente cómo y cuándo se romperá esta lasaña es un pesadilla para los ingenieros. Si fallan, el avión o el coche podrían tener un problema grave.

Aquí es donde entra este artículo. Los autores han creado un nuevo "simulador de videojuego" matemático (un modelo computacional) para predecir esas roturas con mucha precisión. Vamos a desglosarlo con analogías sencillas:

1. El Problema: ¿Cómo se rompe la lasaña?

Antes, los modelos eran como intentar adivinar si una lasaña se romperá mirando solo la superficie. No entendían bien la diferencia entre romper los espaguetis (las fibras) y romper la salsa (la matriz). Además, si tenías muchas capas con diferentes direcciones, los modelos antiguos se volvían locos o eran demasiado lentos.

2. La Solución: El Modelo de "Dos Focos" (Multi-Phase-Field)

Los autores proponen un modelo que usa dos "focos" o sensores independientes para vigilar la salud de la lasaña:

• Foco 1 (Fibras): Vigila si los espaguetis se están rompiendo.
• Foco 2 (Matriz): Vigila si la salsa se está agrietando.

Esto es genial porque en la vida real, a veces la salsa se rompe primero, y otras veces las fibras ceden. Al tener dos sensores separados, el modelo puede ver cómo un daño afecta al otro. Es como tener dos cámaras de seguridad: una para el techo y otra para el suelo.

3. La Regla del Juego: La Teoría de Puck

Para saber cuándo se rompe algo, el modelo sigue un manual de instrucciones muy famoso llamado Teoría de Puck. Imagina que es como un árbitro de fútbol muy estricto.

• Si la tensión en las fibras es demasiado alta, el árbitro pita: "¡Fallo de fibra!".
• Si la presión en la salsa es demasiado fuerte, pita: "¡Fallo de matriz!".
  El modelo usa estas reglas para decidir exactamente en qué momento empieza la grieta.

4. El Truco de Magia: La "Superposición de Mallas" (Mesh Overlay)

Aquí está la parte más creativa. Imagina que tienes una hoja de papel transparente (la primera capa de la lasaña) y la pones sobre una mesa. Luego, tomas otra hoja transparente (la segunda capa) con un dibujo diferente y la pones exactamente encima de la primera, sin moverla.

• El truco: Todas las capas comparten los mismos puntos de anclaje (nodos), pero cada una tiene su propia "ropa" interna (dirección de las fibras).
• El resultado: En lugar de tener que construir un modelo 3D gigante y pesado (como una lasaña real en 3D), el modelo calcula todo en 2D, pero "apilando" las capas virtualmente. Esto hace que el cálculo sea muy rápido y eficiente, como usar una calculadora en lugar de construir un superordenador.

5. ¿Qué probaron? (Los Experimentos)

Para ver si su "simulador" funciona, lo pusieron a prueba contra la realidad en cuatro escenarios clásicos:

1. Tirar y empujar tiras de material: Como estirar una goma elástica hasta que se rompe.
2. Agujero en el medio (Open Hole): Imagina una tabla de madera con un agujero en el centro. Al tirar de ella, las grietas suelen salir desde el agujero. El modelo predijo exactamente cómo se expandían esas grietas.
3. Tensión compacta: Un tipo de prueba donde se tira de un material con una muesca para ver cuánto aguanta antes de partirse.
4. Doble muesca: Como tirar de una hoja de papel que tiene dos cortes en los bordes.

6. El Resultado Final

El modelo funcionó increíblemente bien.

• Calidad: Podía "ver" las grietas formándose de la misma manera que en los experimentos reales (caminando a lo largo de las fibras o atravesándolas).
• Cantidad: Podía predecir casi exactamente cuánta fuerza soportaría el material antes de romperse.

En Resumen

Este artículo presenta una herramienta digital inteligente que permite a los ingenieros diseñar aviones y coches más seguros y ligeros sin tener que destruir tantos prototipos reales.

En lugar de adivinar, ahora pueden usar este modelo para ver "dentro" de las capas del material, entender si se romperán las fibras o la resina, y hacerlo todo de forma rápida gracias a su truco de superponer las capas como si fueran transparencias. Es un gran paso para hacer que las estructuras compuestas sean más predecibles y seguras.

Resumen técnico

Resumen Técnico: Modelo Multi-Campo de Fase para Laminados Compuestos Reforzados con Fibra Basado en la Teoría de Fallo de Puck

1. Planteamiento del Problema

Los materiales compuestos reforzados con fibra (FRC) son fundamentales en industrias como la aeroespacial y automotriz debido a su alta relación resistencia-peso. Sin embargo, predecir su fallo es extremadamente complejo debido a la naturaleza multifacética de sus mecanismos de rotura, que incluyen:

• Agrietamiento de la matriz intralaminar (transversal y longitudinal).
• Fallo de las fibras (tracción y compresión).
• Delaminación (fallo interlaminar).

Los métodos numéricos existentes (como elementos cohesivos, XFEM o peridinámica) a menudo enfrentan dificultades para capturar la nucleación y propagación de grietas en laminados con secuencias de apilamiento complejas, especialmente en configuraciones bidimensionales (in-plane). Además, muchos modelos de campo de fase (Phase-Field, PF) existentes no logran predecir correctamente la iniciación del fallo bajo cargas de compresión o requieren ajustes ad hoc para calibrar las escalas de longitud. Existe una necesidad crítica de un modelo que integre criterios de fallo fenomenológicos robustos (como la teoría de Puck) dentro de un marco de campo de fase para simular la interacción entre el fallo de la fibra y el de la matriz en laminados multidireccionales.

2. Metodología

Los autores proponen un modelo bidimensional multi-campo de fase que combina la teoría de variacional de fractura con la teoría de fallo de Puck. Los componentes clave de la metodología son:

• Variables de Campo de Fase Independientes: Se introducen dos variables de campo de fase escalares, $d_f$ (fallo de fibra) y $d_m$ (fallo inter-fibra/matriz), para modelar separadamente los modos de fallo dominados por la fibra y por la matriz. Esto permite capturar la interrelación entre ambos tipos de daño.
• Integración de la Teoría de Puck: A diferencia de modelos anteriores que dependen solo de la energía elástica, este modelo utiliza la teoría de Puck para determinar la iniciación del daño. Los factores de exposición de Puck ($F_f$ para fibra y $F_m$ para matriz) se utilizan para calcular la energía de conducción del daño ($H_f, H_m$), asegurando que el fallo comience bajo las condiciones de tensión/compresión y cizalladura correctas definidas por Puck.
• Método de Superposición de Mallas (Mesh Overlay): Para simular laminados con múltiples orientaciones de fibra sin aumentar drásticamente la complejidad computacional (evitando modelos 3D completos), se utiliza una técnica de superposición de mallas.
  • Se utiliza una malla base compartida para todas las capas.
  • Cada capa (ply) tiene su propio conjunto de elementos "dummy" con una orientación de fibra específica ($\theta_i$).
  • Los desplazamientos están acoplados entre capas, pero las tensiones y el daño se calculan de forma desacoplada según la orientación local de cada capa.
  • Limitación: Este enfoque 2D omite la delaminación (fallo interlaminar), asumiendo interfaces perfectas.
• Formulación Variacional: El modelo se basa en la minimización de un funcional de energía total que incluye la energía de deformación degradada y la energía de superficie difusa, garantizando la consistencia termodinámica mediante la desigualdad de Clausius-Duhem y las condiciones de Karush-Kuhn-Tucker (KKT) para la irreversibilidad del daño.
• Implementación: El modelo se implementó en el código comercial Abaqus mediante un elemento de usuario (UEL) que aprovecha el solucionador Newton-Raphson integrado.

3. Contribuciones Clave

1. Modelo Multi-Campo de Fase para Laminados: Es uno de los primeros modelos de campo de fase que considera más de dos capas de secuencia de laminado para estudiar la propagación de grietas en el plano (in-plane), superando las limitaciones de estudios previos restringidos a unidireccionales o laminados cruzados simples.
2. Acoplamiento con Puck: La integración exitosa de los criterios de fallo de Puck dentro del marco de campo de fase resuelve el problema de la predicción de la nucleación de grietas, algo que los modelos puramente basados en energía a menudo fallan en cargas de compresión.
3. Eficiencia Computacional: El uso del método de superposición de mallas permite simular laminados complejos con diferentes secuencias de apilamiento utilizando una malla 2D, ofreciendo una alternativa eficiente a los costosos modelos 3D.
4. Validación Exhaustiva: El modelo se valida contra cuatro casos de prueba de referencia (benchmarks) con datos experimentales de la literatura, cubriendo tanto cargas de tracción como de compresión.

4. Resultados

El modelo se evaluó mediante cuatro ejemplos de referencia:

• Pruebas de Probeta (Coupon Tests): Se simularon probetas con secuencias $[(\theta/-\theta)_8]_s$ bajo tracción y compresión. El modelo reprodujo con precisión los patrones de agrietamiento de la matriz y las curvas fuerza-desplazamiento, mostrando un buen acuerdo con los datos experimentales de Flatscher et al., aunque con una ligera sobreestimación de la carga máxima debido a la ausencia de plasticidad en el modelo.
• Tracción con Agujero (Open Hole Tension - OHT):
  • Para la secuencia $[0^4/90^4]_s$, el modelo capturó correctamente la propagación simétrica de grietas de la matriz paralelas a la carga, seguida del fallo de las fibras.
  • Para la secuencia $[45^4/-45^4]_s$, el modelo reprodujo el fenómeno de "blindaje de grietas" (crack shielding) y la propagación asimétrica, coincidiendo cualitativamente con los experimentos.
• Compacidad de Tensión (Compact Tension): Se analizaron laminados bloqueados $[0^2/90^2]_{2s}$ y $[0^4/90^4]_{2s}$. El modelo predijo bien el comportamiento lineal inicial y la propagación de grietas. Las discrepancias en la carga máxima se atribuyeron a la falta de modelado de la delaminación, que en los experimentos disipa energía.
• Tracción con Doble Entalla (Double-Notched Tension): Se estudiaron secuencias cruzadas, cuasi-isotrópicas y mixtas. El modelo demostró su capacidad para simular la competencia entre el fallo de la fibra y de la matriz, así como la nucleación de grietas subcríticas, mostrando una correlación excelente con los resultados experimentales de [30].

En todos los casos, el modelo logró predecir tanto cualitativamente (patrones de grietas) como cuantitativamente (cargas máximas y curvas fuerza-desplazamiento) el comportamiento del material.

5. Significancia

Este trabajo representa un avance significativo en la mecánica computacional de materiales compuestos al:

• Cerrar la brecha entre la teoría de fallo fenomenológica y los métodos de campo de fase, permitiendo una predicción más física de la iniciación del daño.
• Proveer una herramienta eficiente para el diseño de laminados complejos, donde la interacción entre diferentes orientaciones de fibra es crítica.
• Establecer un estándar de validación mediante la comparación directa con múltiples experimentos de la literatura.
• Promover la transparencia científica, ya que los autores ponen a disposición el código fuente y los datos de los ejemplos, facilitando la investigación futura y la replicabilidad.

Aunque el modelo tiene limitaciones (principalmente la omisión de la delaminación y el comportamiento plástico), ofrece una base sólida y robusta para la evaluación cualitativa y cuantitativa de la fractura en laminados compuestos bajo cargas in-plane.