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

Resumo Técnico: Uma Estrutura de Modelagem de Processo por Elementos Finitos para Prever Tensões Residuais e Microfissuração em Compósitos de Matriz Cerâmica

Definição do Problema
Os Compósitos de Matriz Cerâmica (CMCs) são críticos para aplicações hipersônicas de próxima geração devido à sua excepcional resistência térmica e integridade estrutural em temperaturas superiores a 1200–1600 °C. No entanto, o principal método de fabricação, a Infiltração de Polímero e Pirólise (PIP), sofre com a formação significativa de defeitos, especificamente microfissuras e vazios. Esses defeitos necessitam de múltiplos e dispendiosos ciclos de reinfiltração para alcançar a densificação total. As causas raízes são complexas, envolvendo reações termoquímicas acopladas e a evolução simultânea das propriedades mecânicas do polímero pré-cerâmico (PCP) durante a pirólise. Esses processos geram grandes tensões residuais que são difíceis de quantificar experimentalmente devido às temperaturas extremas de processamento e à ocorrência simultânea de cinética de reação, perda de massa, densificação, contração volumétrica e transporte de gases. Consequentemente, há uma necessidade crítica de estruturas computacionais baseadas em física para simular esses mecanismos acoplados e prever a formação de defeitos.

Metodologia
Este trabalho apresenta uma estrutura de elementos finitos (FE) baseada em física, projetada para a manufatura digital de CMCs, focando especificamente na simulação em microescala. A abordagem utiliza um novo processo de duas etapas que vincula a cura (estágio de cura) diretamente à pirólise subsequente.

• Física Regente: A estrutura resolve três equações acopladas: uma lei constitutiva elástica linear instantânea, uma equação de calor transiente com um termo de fonte exotérmica e um modelo cinético de cura e pirólise. A deformação total é decomposta em componentes mecânicos, térmicos e de contração. A contração é modelada como uma função linear do grau de cura para a etapa de cura e como uma função do rendimento volumétrico (derivado da perda de massa e densificação) para a etapa de pirólise.
• Modelagem de Dano: Um modelo de dano progressivo baseado na teoria de banda de fissura de Bažant e Oh é incorporado para capturar a microfissuração in-situ. A iniciação do dano segue um critério de tensão principal máxima, e o comportamento pós-pico é modelado via separação-tração linear com amolecimento de deformação. Um parâmetro de dano escalar ($D$) degrada o tensor de complacência.
• Implementação: As equações regentes são implementadas no ambiente de FE do Abaqus usando sub-rotinas definidas pelo usuário (UMAT e UMATHT). A estrutura emprega uma análise de dois estágios: uma análise de cura estabelece o estado inicial (tensões residuais, dano e propriedades do material), que é então importada como a condição inicial para a análise de pirólise.
• Modelagem de Material: Devido aos dados experimentais limitados para o sistema específico (matriz SMP-10 com fibras de carbono AS4), o estudo constrói uma descrição de material substituta. As propriedades de cura são derivadas de sistemas análogos (fenólicos SC-1008 e epóxi RIM R135-H1366), enquanto a cinética de pirólise e a contração baseiam-se em modelos estabelecidos para SMP-10 por Key et al. As propriedades do material durante a pirólise são assumidas como evoluindo em função da mudança de densidade cerâmica.
• Geometria: A estrutura é demonstrada em um Elemento de Volume Representativo (RVE) contendo 20 fibras de carbono AS4 dispersas aleatoriamente (fração de volume de 45%) dentro de uma matriz SMP-10, sujeita a condições de contorno periódicas.

Principais Contribuições

1. Estrutura Unificada de Duas Etapas: O artigo apresenta um elo direto entre a cura e a pirólise, permitendo que o estado de tensão residual e o dano acumulado durante a cura sirvam como as condições iniciais para a simulação de pirólise. Isso contrasta com trabalhos anteriores que frequentemente tratavam essas etapas separadamente ou focavam em mesoescalas onde as fibras não eram explicitamente modeladas.
2. Integração de Dano Progressivo: A inclusão de um modelo de dano progressivo permite a simulação de microfissuração in-situ ao longo de ambas as etapas de processamento, capturando a interação entre a evolução da tensão e a degradação do material.
3. Premissas Fisicamente Consistentes: Na ausência de dados completos de material, os autores desenvolveram um modelo de proxy robusto para as propriedades do SMP-10, garantindo consistência física ao derivar a evolução das propriedades a partir das mudanças de densidade e modelos cinéticos validados.
4. Resolução de Microescala: O modelo resolve explicitamente a interface fibra-matriz na microescala, permitindo a observação de concentradores de tensão e localização de danos que são frequentemente homogeneizados em modelos de maior escala.

Resultados
As simulações da microestrutura do CMC revelaram o seguinte:

• Estágio de Cura: As tensões residuais geradas durante a cura foram relativamente baixas (atingindo pico de ~7 MPa de tensão durante a rampa e ~65 MPa em temperatura ambiente após o resfriamento). O início do dano foi mínimo, com o parâmetro de dano médio volumétrico ($D$) permanecendo abaixo de 0,015 ao final da cura. Concentrações de tensão foram observadas principalmente onde as fibras estavam em proximidade estreita.
• Estágio de Pirólise: A microfissuração significativa da matriz foi impulsionada pela contração induzida pela pirólise. O dano iniciou-se cedo na rampa térmica (por volta de 285 °C, coincidindo com o início da decomposição rápida). O parâmetro de dano médio volumétrico subiu acentuadamente, excedendo 0,7 antes que a temperatura atingisse a metade do alvo final de 1200 °C.
• Tensão e Contração: As tensões de tração de pico ocorreram cedo na rampa térmica (aprox. 300–400 °C), atingindo ~55 MPa antes que o dano progressivo limitasse o acúmulo adicional de tensão. A simulação previu uma contração volumétrica final de aproximadamente 8%.
• Evolução de Propriedades: O módulo elástico final e a resistência à tração da matriz foram significativamente degradados (~160 GPa e 1,4 GPa, respectivamente) em comparação com as propriedades de bulk teóricas do SiC totalmente denso, devido ao dano acumulado. A microestrutura resultante exibiu morfologias de fragmentação consistentes com observações da literatura.

Significância e Alegações
Os autores posicionam este trabalho como uma prova de conceito para uma estrutura computacional que vincula a cura e a pirólise para prever a formação de defeitos em CMCs. O estudo alega estabelecer uma base para "modelagem preditiva informada pelo processo" que pode investigar sistematicamente como as condições de processamento (ciclos térmicos, parâmetros de PIP) influenciam a tensão residual, a contração e a microfissuração.

O artigo nota modestamente que, devido aos dados de material limitados, os resultados atuais dependem de descrições e suposições de proxy. No entanto, afirma que a estrutura fornece uma plataforma necessária para compreender os mecanismos acoplados de formação de defeitos. Os autores concluem que, à medida que dados de caracterização de material mais precisos se tornarem disponíveis e outras físicas (como o transporte de gás de pirólise) forem incorporadas, esta abordagem tem o potencial de guiar o design de estratégias de manufatura para melhorar o desempenho pós-processamento dos CMCs. O trabalho não alega resolver todos os problemas de otimização imediatamente, mas sim fornecer o conjunto de ferramentas computacionais essencial para iniciar tais investigações.