Composite Failure Criteria Simulator Back
Structural Analysis / Composites

Composite Material Failure Criteria Simulator

Calculate ply stresses in CFRP/GFRP laminates and evaluate failure using Tsai-Wu, Hashin, and Max Stress criteria. Fiber/matrix failure modes and failure envelopes displayed in real time.

Parameters
Fiber Angle θ
°
In-plane Load Nₓ
N/mm
In-plane Load Nᵧ
N/mm
In-plane Load Nₓᵧ
N/mm
Material
Failure Criterion
Number of Plies
Results
Failure index FI
Strength ratio R = 1/√FI
σ₁ (MPa)
σ₂ (MPa)
τ₁₂ (MPa)
FailureMode
Failure index bar (FI)SAFE
Env
Sensitivity Sweep
Theory & Key Formulas
$F_1\sigma_1 + F_2\sigma_2 + F_{11}\sigma_1^2$
$+ F_{22}\sigma_2^2 + F_{66}\tau_{12}^2 + 2F_{12}\sigma_1\sigma_2 \leq 1$

F₁=1/Xt−1/Xc, F₁₁=1/(XtXc)
Xt, Xc: fiber tensile/compressive strength
Yt, Yc: matrix tensile/compressive strength

What is a Composite Failure Criterion?

🙋
What exactly is a "failure criterion" for composites? Why can't we just use the maximum stress like we do for steel?
🎓
Great question! Basically, composites like CFRP are anisotropic—their strength depends on direction. A single stress value doesn't tell the whole story. In practice, failure depends on the interaction between stresses along the fiber, across the fiber, and in shear. Try the simulator: select "Max Stress" and apply a load. You'll see it fails only when one stress hits its limit, ignoring interactions.
🙋
Wait, really? So the other criteria account for interactions? What's the difference between Hashin and Tsai-Wu?
🎓
Exactly. Hashin's criterion separates failure modes—like fiber rupture vs. matrix cracking. It's more physical. Tsai-Wu is a single, smooth mathematical surface that blends all interactions. A common case is a biaxial load: Tsai-Wu might predict failure at a lower load because it accounts for stress coupling. Switch between them in the simulator and watch how the failure envelope changes shape.
🙋
So which one should I trust? And what's this "Failure Index" number the simulator shows?
🎓
In practice, engineers use both! Hashin is great for diagnosing how a ply fails (useful for damage tolerance). Tsai-Wu is simpler for initial sizing. The Failure Index (FI) is key: FI = 1 means you're on the failure boundary. FI < 1 is safe. Move the stress sliders and watch the FI update in real-time. For design, you'd aim for a Strength Ratio (R) > 1.5, meaning your stresses are only about 2/3 of the failure level.

Physical Model & Key Equations

The Tsai-Wu criterion is a generalized quadratic failure model that accounts for stress interactions in anisotropic materials. It creates a single, smooth failure envelope in stress space.

$$F_1\sigma_1 + F_2\sigma_2 + F_{11}\sigma_1^2 + F_{22}\sigma_2^2 + F_{66}\tau_{12}^2 + 2F_{12}\sigma_1\sigma_2 = 1$$

Where $\sigma_1$ is stress along fibers, $\sigma_2$ is transverse stress, and $\tau_{12}$ is in-plane shear. The coefficients F are derived from basic strengths: $F_1 = 1/X_t - 1/X_c$, $F_{11}= 1/(X_t X_c)$, where $X_t, X_c$ are tensile and compressive strengths in the fiber direction. The interaction term $F_{12}$ is often determined experimentally.

Hashin's criterion separates failure into distinct, physically meaningful modes. This allows engineers to understand whether failure is driven by fiber breakage, matrix cracking, or shear.

$$ \begin{align*}\text{Fiber Tensile}(\sigma_1 \ge 0): & \quad \left(\frac{\sigma_1}{X_t}\right)^2 + \left(\frac{\tau_{12}}{S}\right)^2 = 1 \\[4pt] \text{Matrix Tensile}(\sigma_2 \ge 0): & \quad \left(\frac{\sigma_2}{Y_t}\right)^2 + \left(\frac{\tau_{12}}{S}\right)^2 = 1 \end{align*} $$

Here, $Y_t, Y_c$ are transverse strengths, and $S$ is the shear strength. The simulator calculates each mode separately. The one that reaches 1 first dictates the failure mode, which is crucial for post-failure analysis.

Frequently Asked Questions

The Tsai-Wu criterion uses a quadratic interaction formula of stress components to determine only the presence or absence of failure, whereas the Hashin criterion evaluates fiber failure and matrix failure separately. This tool calculates both simultaneously and displays in real time which mode causes failure.
The inside of the envelope is the safe region, and the outside is the failure region. It is displayed as an ellipse or polygon on the σ1-σ2 plane, and is asymmetric because tensile and compressive strengths differ. By moving the stress point with the mouse, the failure margin at that position can be confirmed numerically.
Values for typical CFRP/GFRP are listed in material databases or manufacturers' technical documents. If experimental values are not available, you can estimate them from the fiber volume fraction or refer to the default values in this tool (e.g., T300/5208).
Input the in-plane stresses (Nx, Ny, Nxy) acting on the entire laminate or the strain of each ply. This tool transforms the stresses to the material principal axis direction of each ply based on orthotropic elasticity and performs failure evaluation.

Real-World Applications

Aerospace Wing & Fuselage Design: CFRP skins on aircraft wings experience complex biaxial tension and shear during flight. Engineers use these criteria to optimize ply angles (try changing the "Ply Angle" selector) to ensure the laminate doesn't fail under gust loads or during pressurization cycles, targeting a safe Strength Ratio.

Wind Turbine Blade Design: GFRP blades undergo massive bending and torsional loads. Hashin's criterion is vital here to distinguish between compressive fiber buckling on one side and matrix cracking on the other, informing where to add extra material or different ply sequences.

Automotive Chassis & Body Panels: In high-performance cars, carbon fiber tubs must absorb crash energy. Simulating failure under multi-axial impact loads (by adjusting the $\sigma_1$ and $\sigma_2$ sliders together) helps design a controlled, predictable crushing failure rather than sudden brittle fracture.

Sports Equipment (Bikes, Rackets): For a bicycle frame, the failure envelope under combined pedaling forces (axial) and cornering (bending/torsion) determines the layup. Using the simulator to find the "weakest" criterion for a given load case ensures the product is both lightweight and durable.

Common Misconceptions and Points to Note

When you start using this simulator, there are several pitfalls that beginners in CAE often encounter. The first is the misconception that if the Failure Index (FI) exceeds 1, the entire part immediately shatters. In actual composite materials, especially laminates, "progressive failure" occurs where other plies take up the load even if one ply fails. An FI=1.2 is a sign that "failure begins in that ply," not the ultimate strength of the entire component. In practice, it's common to set a safety factor that keeps the FI around 0.5~0.8.

Next is input errors for material constants. For example, if you mistakenly swap the tensile strength (Xt) and compressive strength (Xc), the results become meaningless. For CFRP, the compressive strength in the fiber direction is often about 60-70% of the tensile strength. For instance, if Xt=1500MPa, then Xc=900~1000MPa is a good guideline. If this relationship is reversed, it's a sign you should review your inputs.

Finally, the tendency to search for the "strongest failure criterion". Tsai-Wu is convenient but doesn't tell you the failure mode. On the other hand, Hashin indicates the mode but may not fully account for fiber-matrix interactions. The key is "at which stage of design, and what do you want to know?" A practical approach is to use Tsai-Wu for a rough FI during conceptual design and then use Hashin to identify weak failure modes during detailed design.

How to Use

  1. Input lamina orientation angle (theta) in degrees relative to load direction; typical aerospace values range 0° to 90°
  2. Enter in-plane stresses: normal stress nx (MPa) along fiber direction and transverse stress ny (MPa); example CFRP laminate: nx=500 MPa, ny=20 MPa
  3. Select failure criterion from dropdown: Tsai-Wu (interactive quadratic), Hashin (fiber/matrix separation), or Maximum Stress (linear safety factors)
  4. Run simulation to compute ply stresses via coordinate transformation, then evaluate against material allowables
  5. Review failure mode output: fiber tension, fiber compression, matrix cracking, or matrix crushing for CFRP with Xt=1500 MPa, Xc=1200 MPa, Yt=50 MPa, Yc=200 MPa

Worked Example

Carbon fiber epoxy (AS4/3501-6) lamina at theta=30° under biaxial loading: nx=400 MPa, ny=15 MPa. Transform to material axes using rotation matrix. Ply stresses: sigma_1≈365 MPa, sigma_2≈32 MPa, tau_12≈-166 MPa. Tsai-Wu criterion with F1=1/1500, F2=1/50, F11=1/(1500×1200), F22=1/(50×200), F66=1/300² yields interaction coefficient ≈0.68 (safe, reserve factor ≈1.47). Hashin checks fiber tension (365/1500=0.24, safe) and matrix cracking (32/50=0.64, approaching limit). Maximum Stress method shows tau_12/Sxy exceeds 0.55 safety ratio.

Practical Notes

  1. Angle sensitivity: Quasi-isotropic [0/±45/90]s laminates reduce peak ply stress by 40% versus unidirectional at critical theta=45° compared to theta=0°
  2. Hashin criterion detects matrix failure modes earlier (lower safety margin) in low-fiber-angle plies under transverse loading; use for pressure vessels, pipes
  3. Tsai-Wu includes interaction term F12; critical for off-axis loading in wind turbine blades where combined bending/torsion occurs at theta=±30°
  4. Temperature effects: reduce Yt by ~8% per 10°C above Tg for epoxy; critical for aerospace high-altitude missions