⚡ electrical engineering

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

This paper presents a physics-based finite element framework that links curing and pyrolysis stages to predict residual stresses, shrinkage, and microcracking in ceramic matrix composites, demonstrating its capability to simulate defect formation during polymer infiltration and pyrolysis processing.

Original authors: Michael N. Olaya, Marianna Maiaru

Published 2026-06-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Michael N. Olaya, Marianna Maiaru

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to bake a very special, super-strong ceramic cake. But instead of flour and eggs, you are using a liquid plastic that turns into rock when heated. This is how scientists make Ceramic Matrix Composites (CMCs), materials used in things like hypersonic jets that need to survive extreme heat.

The problem is that this "baking" process is tricky. As the plastic turns into rock, it shrinks and changes chemically. This creates invisible internal tensions (like a rubber band being pulled too tight) that cause tiny cracks and holes (voids) to form inside the material. If these cracks form, the final product is weak and useless.

This paper introduces a computer simulation framework—a digital "virtual oven"—that helps engineers predict exactly where and when these cracks will happen before they ever build a real part.

Here is how the authors' "virtual oven" works, broken down into simple steps:

1. The Two-Step "Baking" Process

The authors realized you can't just jump straight to the high heat. You have to simulate the process in two distinct stages:

Stage 1: The "Setting" Phase (Curing): Imagine the liquid plastic is like gelatin. First, you have to let it set into a solid jelly. The computer simulates this "setting" process (called C-staging). During this time, the material hardens, but it doesn't shrink much yet. The computer checks if any tiny cracks form just from the material hardening.
Stage 2: The "Rocking" Phase (Pyrolysis): Once the jelly is set, you turn up the heat to turn it into rock. This is where the magic (and the danger) happens. The material loses weight (gas escapes), shrinks significantly, and hardens into a ceramic. The computer tracks how the material shrinks and how the internal pressure builds up during this intense heat.

2. The "Digital Twin" Microstructure

To make the simulation realistic, the authors didn't just look at a block of material. They built a digital twin of the material's tiny internal structure.

The Fibers: Think of these as the "skeleton" or the "reinforcing steel bars" inside the cake. They used a model of carbon fibers (AS4) that stay strong and don't change much.
The Matrix: This is the "filling" or the "cake batter" (SMP-10) that surrounds the fibers. This is the part that shrinks and cracks.
The Simulation: They created a tiny, repeating digital box containing 20 fibers surrounded by this "batter" and watched what happened as they heated it up in the computer.

3. The "Crack Detector"

The most clever part of their framework is a progressive damage model.

Imagine the material has a "breaking point." As the computer simulates the heating, it constantly checks the stress on the material.
If the stress gets too high (like pulling a rubber band until it snaps), the computer "breaks" that tiny spot in the simulation.
Once a spot breaks, the computer updates the material's properties, making it weaker in that spot, which changes how the stress spreads to the neighbors. This allows the simulation to show a network of tiny cracks forming, just like real life.

4. What They Found (The Results)

Since they didn't have all the perfect real-world data for this specific material, they used the best available information from other studies and made logical guesses (like using a similar material as a "stand-in").

Here is what their digital oven revealed:

The Big Shrink: When the material turned from plastic to ceramic, it shrank by about 8%. That's a lot for a solid object!
The Stress Spike: The biggest danger wasn't at the very end of the baking. The highest stress (and the start of the most cracking) happened early on, when the temperature was just starting to climb (around 300–400°C).
The Cracks: Because the "batter" shrank so much while the "skeleton" (fibers) stayed the same size, the batter got pulled apart. By the end of the simulation, about 85% of the matrix material had some level of cracking.
The Outcome: The final material was much weaker and less stiff than a perfect, uncracked piece of ceramic would be.

Why This Matters

The authors aren't saying this computer model is perfect yet; they admit they had to fill in some gaps with assumptions because real-world data is hard to get. However, they have built a foundation.

Think of this framework as a flight simulator for manufacturing. Just as pilots use simulators to learn how to handle a plane in a storm without crashing a real one, engineers can now use this tool to test different heating recipes. They can ask, "What if we heat it slower?" or "What if we change the fiber arrangement?" to see if they can stop the cracks from forming, all inside the computer before spending money on expensive real-world experiments.

In short: They built a digital tool that simulates the "baking" of super-strong ceramics, showing that the process causes the material to shrink and crack significantly, and proving that we need to carefully control the heating process to fix these defects.

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

Problem Statement
Ceramic Matrix Composites (CMCs) are critical for next-generation hypersonic applications due to their exceptional thermal resistance and structural integrity at temperatures exceeding 1200–1600 °C. However, the primary manufacturing method, Polymer Infiltration and Pyrolysis (PIP), suffers from significant defect formation, specifically microcracks and voids. These defects necessitate multiple, costly re-infiltration cycles to achieve full densification. The root causes are complex, coupled thermochemical reactions and the concurrent evolution of mechanical properties within the pre-ceramic polymer (PCP) during pyrolysis. These processes generate large residual stresses that are difficult to quantify experimentally due to extreme processing temperatures and the simultaneous occurrence of reaction kinetics, mass loss, densification, volumetric shrinkage, and gas transport. Consequently, there is a critical need for physics-based computational frameworks to simulate these coupled mechanisms and predict defect formation.

Methodology
This work introduces a physics-based finite element (FE) framework designed for the digital manufacturing of CMCs, specifically focusing on microscale simulation. The approach utilizes a novel two-step process linking curing (C-staging) directly to subsequent pyrolysis.

Governing Physics: The framework solves three coupled equations: an instantaneous linear elastic constitutive law, a transient heat equation with an exothermic source term, and a kinetic model for curing and pyrolysis. The total strain is decomposed into mechanical, thermal, and shrinkage components. Shrinkage is modeled as a linear function of the degree of cure for the curing step and as a function of volumetric yield (derived from mass loss and densification) for the pyrolysis step.
Damage Modeling: A progressive damage model based on Bažant and Oh's crack band theory is incorporated to capture in-situ microcracking. Damage initiation follows a maximum principal stress criterion, and post-peak behavior is modeled via linear traction-separation with strain softening. A scalar damage parameter ( $D$ ) degrades the compliance tensor.
Implementation: The governing equations are implemented in the Abaqus FE environment using user-defined subroutines (UMAT and UMATHT). The framework employs a two-stage analysis: a curing analysis establishes the initial state (residual stresses, damage, and material properties), which is then imported as the initial condition for the pyrolysis analysis.
Material Modeling: Due to limited experimental data for the specific system (SMP-10 matrix with AS4 carbon fibers), the study constructs a stand-in material description. Curing properties are derived from analogous systems (phenolic SC-1008 and epoxy RIM R135-H1366), while pyrolysis kinetics and shrinkage are based on established models for SMP-10 by Key et al. Material properties during pyrolysis are assumed to evolve as a function of the changing ceramic density.
Geometry: The framework is demonstrated on a Representative Volume Element (RVE) containing 20 randomly dispersed AS4 carbon fibers (45% volume fraction) within an SMP-10 matrix, subject to periodic boundary conditions.

Key Contributions

Unified Two-Step Framework: The paper presents a direct link between curing and pyrolysis, allowing the residual stress state and damage accumulated during curing to serve as the initial conditions for the pyrolysis simulation. This contrasts with previous works that often treated these stages separately or focused on mesoscales where fibers were not explicitly modeled.
Progressive Damage Integration: The inclusion of a progressive damage model allows for the simulation of in-situ microcracking throughout both processing stages, capturing the interaction between stress evolution and material degradation.
Physically Consistent Assumptions: In the absence of complete material data, the authors developed a robust proxy model for SMP-10 properties, ensuring physical consistency by deriving property evolution from density changes and validated kinetic models.
Microscale Resolution: The model explicitly resolves the fiber-matrix interface at the microscale, enabling the observation of stress risers and damage localization that are often homogenized in larger-scale models.

Results
Simulations of the CMC microstructure revealed the following:

Curing Stage: Residual stresses generated during curing were relatively low (peaking at ~7 MPa tension during the ramp and ~65 MPa at room temperature after cooldown). Damage initiation was minimal, with the volume-averaged damage parameter ( $D$ ) remaining below 0.015 at the end of the cure. Stress concentrations were observed primarily where fibers were in close proximity.
Pyrolysis Stage: Significant matrix microcracking was driven by pyrolysis-induced shrinkage. Damage initiated early in the thermal ramp (around 285°C, coinciding with the onset of rapid decomposition). The volume-averaged damage parameter rose steeply, exceeding 0.7 before the temperature reached the halfway point to the final 1200°C target.
Stress and Shrinkage: Peak tensile stresses occurred early in the thermal ramp (approx. 300–400°C), reaching ~55 MPa before progressive damage limited further stress accumulation. The simulation predicted a final volumetric shrinkage of approximately 8%.
Property Evolution: The final elastic modulus and tensile strength of the matrix were significantly degraded (~160 GPa and 1.4 GPa, respectively) compared to the theoretical bulk properties of fully dense SiC, due to the accumulated damage. The resulting microstructure exhibited fragmentation morphologies consistent with literature observations.

Significance and Claims
The authors position this work as a proof-of-concept for a computational framework that links curing and pyrolysis to predict defect formation in CMCs. The study claims to establish a foundation for "predictive, process-informed modeling" that can systematically investigate how processing conditions (thermal cycles, PIP parameters) influence residual stress, shrinkage, and microcracking.

The paper modestly notes that due to limited material data, the current results rely on proxy descriptions and assumptions. However, it asserts that the framework provides a necessary platform for understanding the coupled mechanisms of defect formation. The authors conclude that as more accurate material characterization data becomes available and additional physics (such as pyrolysis gas transport) are incorporated, this approach has the potential to guide the design of manufacturing strategies to improve the as-processed performance of CMCs. The work does not claim to solve all optimization problems immediately but rather provides the essential computational toolset to begin such investigations.