A Finite Element Process Modeling Framework for Predicting Residual Stresses and Microcracking in Ceramic Matrix Composites

Resumen Técnico: Un marco de modelado de procesos mediante elementos finitos para la predicción de tensiones residuales y microfisuración en compuestos de matriz cerámica

Planteamiento del problema
Los compuestos de matriz cerámica (CMC) son críticos para aplicaciones hipersónicas de próxima generación debido a su excepcional resistencia térmica e integridad estructural a temperaturas superiores a 1200–1600 °C. Sin embargo, el método de fabricación principal, la infiltración y pirólisis de polímero (PIP), sufre una formación significativa de defectos, específicamente microfisuras y huecos. Estos defectos requieren múltiples y costosos ciclos de reinfiltración para lograr la densificación completa. Las causas fundamentales son complejas, involucrando reacciones termoquímicas acopladas y la evolución simultánea de las propiedades mecánicas dentro del polímero precerámico (PCP) durante la pirólisis. Estos procesos generan grandes tensiones residuales que son difíciles de cuantificar experimentalmente debido a las temperaturas extremas de procesamiento y la ocurrencia simultánea de la cinética de reacción, pérdida de masa, densificación, contracción volumétrica y transporte de gases. Por consiguiente, existe una necesidad crítica de marcos computacionales basados en la física para simular estos mecanismos acoplados y predecir la formación de defectos.

Metodología
Este trabajo introduce un marco de elementos finitos (FE) basado en la física diseñado para la fabricación digital de CMC, centrándose específicamente en la simulación a microescala. El enfoque utiliza un novedoso proceso de dos pasos que vincula el curado (etapa de C) directamente con la pirólisis posterior.

• Física gobernante: El marco resuelve tres ecuaciones acopladas: una ley constitutiva elástica lineal instantánea, una ecuación de calor transitoria con un término de fuente exotérmica y un modelo cinético para el curado y la pirólisis. La deformación total se descompone en componentes mecánicas, térmicas y de contracción. La contracción se modela como una función lineal del grado de curado para la etapa de curado y como una función del rendimiento volumétrico (derivado de la pérdida de masa y la densificación) para la etapa de pirólisis.
• Modelado de daño: Se incorpora un modelo de daño progresivo basado en la teoría de banda de grieta de Bažant y Oh para capturar la microfisuración in-situ. La iniciación del daño sigue un criterio de tensión principal máxima, y el comportamiento post-pico se modela mediante tracción-separación lineal con ablandamiento por deformación. Un parámetro de daño escalar ($D$) degrada el tensor de cumplimiento.
• Implementación: Las ecuaciones gobernantes se implementan en el entorno de elementos finitos de Abaqus utilizando subrutinas definidas por el usuario (UMAT y UMATHT). El marco emplea un análisis de dos etapas: un análisis de curado establece el estado inicial (tensiones residuales, daño y propiedades del material), que luego se importa como la condición inicial para el análisis de pirólisis.
• Modelado de materiales: Debido a los limitados datos experimentales para el sistema específico (matriz SMP-10 con fibras de carbono AS4), el estudio construye una descripción de material de sustitución. Las propiedades de curado se derivan de sistemas análogos (fenólico SC-1008 y epoxi RIM R135-H1366), mientras que la cinética de pirólisis y la contracción se basan en modelos establecidos para SMP-10 por Key et al. Se asume que las propiedades del material durante la pirólisis evolucionan como una función de la densidad cerámica cambiante.
• Geometría: El marco se demuestra en un Elemento de Volumen Representativo (RVE) que contiene 20 fibras de carbono AS4 dispersas aleatoriamente (fracción de volumen del 45%) dentro de una matriz de SMP-10, sujeto a condiciones de contorno periódicas.

Contribuciones clave

1. Marco unificado de dos pasos: El artículo presenta un vínculo directo entre el curado y la pirólisis, permitiendo que el estado de tensión residual y el daño acumulado durante el curado sirvan como las condiciones iniciales para la simulación de la pirólisis. Esto contrasta con trabajos previos que a menudo trataban estas etapas por separado o se centraban en mesoescalas donde las fibras no se modelaban explícitamente.
2. Integración de daño progresivo: La inclusión de un modelo de daño progresivo permite la simulación de la microfisuración in-situ a través de ambas etapas de procesamiento, capturando la interacción entre la evolución de la tensión y la degradación del material.
3. Supuestos físicamente consistentes: Ante la ausencia de datos completos del material, los autores desarrollaron un modelo de aproximación robusto para las propiedades de SMP-10, asegurando la consistencia física al derivar la evolución de las propiedades a partir de los cambios de densidad y modelos cinéticos validados.
4. Resolución a microescala: El modelo resuelve explícitamente la interfaz fibra-matriz a la microescala, lo que permite observar concentradores de tensión y localización de daños que a menudo se homogeneizan en modelos de mayor escala.

Resultados
Las simulaciones del microestructura del CMC revelaron lo siguiente:

• Etapa de curado: Las tensiones residuales generadas durante el curado fueron relativamente bajas (alcanzando un pico de ~7 MPa de tensión durante la rampa y ~65 MPa a temperatura ambiente tras el enfriamiento). La iniciación de daño fue mínima, con el parámetro de daño promedio por volumen ($D$) manteniéndose por debajo de 0.015 al final del curado. Se observaron concentraciones de tensión principalmente donde las fibras estaban en proximidad cercana.
• Etapa de pirólisis: La microfisuración significativa de la matriz fue impulsada por la contracción inducida por la pirólisis. El daño se inició tempranamente en la rampa térmica (alrededor de 285 °C, coincidiendo con el inicio de la descomposición rápida). El parámetro de daño promedio por volumen aumentó abruptamente, superando el 0.7 antes de que la temperatura alcanzara la mitad del objetivo final de 1200 °C.
• Tensión y contracción: Las tensiones de tracción máximas ocurrieron temprano en la rampa térmica (aproximadamente 300–400 °C), alcanzando ~55 MPa antes de que el daño progresivo limitara la acumulación de tensión adicional. La simulación predijo una contracción volumétrica final de aproximadamente el 8%.
• Evolución de propiedades: El módulo elástico final y la resistencia a la tracción de la matriz se degradaron significativamente (~160 GPa y 1.4 GPa, respectivamente) en comparación con las propiedades de volumen de SiC completamente denso, debido al daño acumulado. La microestructura resultante exhibió morfologías de fragmentación consistentes con las observaciones de la literatura.

Significancia y Reivindicaciones
Los autores posicionan este trabajo como una prueba de concepto para un marco computacional que vincula el curado y la pirólisis para predecir la formación de defectos en los CMC. El estudio afirma establecer una base para un "modelado predictivo informado por el proceso" que pueda investigar sistemáticamente cómo las condiciones de procesamiento (ciclos térmicos, parámetros de PIP) influyen en la tensión residual, la contracción y la microfisuración.

El artículo señala modestamente que, debido a los limitados datos de materiales, los resultados actuales dependen de descripciones y supuestos de aproximación. Sin embargo, afirma que el marco proporciona una plataforma necesaria para comprender los mecanismos acoplados de formación de defectos. Los autores concluyen que, a medida que haya más datos de caracterización de materiales precisos y se incorporen físicas adicionales (como el transporte de gas de pirólisis), este enfoque tiene el potencial de guiar el diseño de estrategias de fabricación para mejorar el rendimiento as-procesado de los CMC. El trabajo no pretende resolver todos los problemas de optimización de inmediato, sino que proporciona el conjunto de herramientas computacionales esenciales para comenzar tales investigaciones.