Fracture Mechanics — CAE Glossary

Category: Glossary | 2026-03-28
CAE visualization for fracture mechanics - technical simulation diagram

What is Fracture Mechanics

🧑‍🎓

What is fracture mechanics? How does it differ from ordinary strength calculation — like checking that stress stays below the yield stress?


🎓

Good question. Ordinary strength calculation assumes the material is sound — meaning cracks don't exist — and checks whether stress remains below an allowable value. But real structures always contain some form of discontinuity: welding defects, fatigue cracks, corrosion pits. That's where fracture mechanics comes in.


🧑‍🎓

So it's a discipline that evaluates structures assuming cracks already exist?


🎓

Exactly. Fracture Mechanics is the science of quantitatively assessing whether an existing crack will remain stable or lead to unstable fracture (sudden rupture). For example, if engineers discover a small crack in an aircraft fuselage, fracture mechanics lets them predict whether the crack will grow to a critical size before the next inspection. This underpins damage-tolerant design.


Stress Intensity Factor K (SIF: Stress Intensity Factor)

🧑‍🎓

I keep hearing about the stress intensity factor K, but isn't it completely different from the stress concentration factor Kt?


🎓

Completely different. The stress concentration factor $K_t$ is a dimensionless ratio showing how many times the nominal stress amplifies near a notch or hole. The stress intensity factor $K$, on the other hand, has units of $\mathrm{MPa}\sqrt{\mathrm{m}}$ and quantifies the strength of the stress field at the crack tip itself. Its definition is:

$$K_I = \sigma \sqrt{\pi a} \cdot Y(a/W)$$

Here $\sigma$ is the remote stress, $a$ is the crack length, and $Y(a/W)$ is a geometry correction factor. For an infinite plate with a central crack, $Y=1$.


🧑‍🎓

What does the subscript "I" in $K_I$ mean?


🎓

It refers to the crack deformation mode. Perfect timing to discuss the different crack modes.


Crack Modes I / II / III

🧑‍🎓

There are different ways for cracks to open?


🎓

Yes, three of them. These are critically important in practice, so learn them well.

Each mode has its own stress intensity factor: $K_I$, $K_{II}$, and $K_{III}$. In real structures, multiple modes often occur together ("mixed-mode"), but Mode I is the most design-critical. The fracture toughness value $K_{Ic}$ is also defined for Mode I.


🧑‍🎓

So if $K_I < K_{Ic}$ the structure is safe, and if $K_I \geq K_{Ic}$ it fractures?


🎓

Basically yes. That's the failure criterion for LEFM (Linear Elastic Fracture Mechanics). But there's an important precondition. Let me explain the difference between LEFM and EPFM.


LEFM vs EPFM

🧑‍🎓

What do LEFM and EPFM stand for? When do we use each?


🎓

LEFM is Linear Elastic Fracture Mechanics, and EPFM is Elastic-Plastic Fracture Mechanics. Roughly speaking, you pick based on whether the plastic zone at the crack tip is small or large.


🧑‍🎓

How small is "small"?


🎓

LEFM is valid under Small-Scale Yielding (SSY) conditions, when the plastic zone size $r_p$ is much smaller than the crack length $a$ or specimen dimensions. Irwin's plastic zone estimate gives:

$$r_p = \frac{1}{2\pi}\left(\frac{K_I}{\sigma_Y}\right)^2 \quad \text{(plane stress)}$$

If $r_p$ is less than about 1/50 of the crack length, LEFM is sufficient. But ductile materials — aluminum alloys, austenitic stainless steels — develop large plastic zones. Then $K$ alone is inadequate, and we need EPFM parameters like the J-integral or CTOD (Crack Tip Opening Displacement).


🧑‍🎓

So high-strength steel uses LEFM and soft materials use EPFM?


🎓

Roughly, but it's not that simple. Even a ductile material can exhibit brittle fracture at low temperature or in thin sections (plane strain). Conversely, a high-strength steel in thick plates can show plasticity. Industry codes like ASME BPVC Section XI and BS 7910 use the FAD (Failure Assessment Diagram) to handle both LEFM and EPFM in a unified framework.


J-Integral

🧑‍🎓

What exactly is the J-integral? The mathematics is a bit intimidating…


🎓

Intuitively, it's the energy released when a crack extends by a unit area. Defined by Rice (1968) as a path integral around the crack tip:

$$J = \int_\Gamma \left( W \, dy - \mathbf{T} \cdot \frac{\partial \mathbf{u}}{\partial x} \, ds \right)$$

Here $W$ is strain energy density, $\mathbf{T}$ is the surface traction vector, and $\mathbf{u}$ is displacement.


🧑‍🎓

Does J give the same value regardless of where I draw the path?


🎓

That's the magic — J is path-independent. You get the same value whether you compute it very close to the crack tip or far from it. This is huge for FEM: even if your crack-tip mesh isn't perfectly refined, you can compute J accurately on outer contours. Abaqus and other solvers automatically compute J over multiple contours so you can verify convergence.


🧑‍🎓

How are J and K related?


🎓

Under small-scale yielding, there's a direct relationship:

$$J = \frac{K_I^2}{E'} \quad \text{where} \quad E' = \begin{cases} E & \text{(plane stress)} \\ \frac{E}{1-\nu^2} & \text{(plane strain)} \end{cases}$$

So within LEFM's domain, $J$ and $K_I$ are equivalent. But once large-scale yielding sets in, $K$ breaks down conceptually, and only $J$ remains physically meaningful. That's why EPFM uses $J_{Ic}$ (critical J-value) as the fracture toughness metric.


XFEM (Extended Finite Element Method)

🧑‍🎓

XFEM is mentioned everywhere lately. How does it differ from traditional FEM for crack analysis?


🎓

In traditional FEM, you must align mesh element boundaries along the crack surface. The crack tip has a $1/\sqrt{r}$ stress singularity, so you place quarter-point singular elements (with mid-side nodes at 1/4 position) there. Every time the crack grows, you remesh. In 3D, this becomes incredibly tedious.


🧑‍🎓

XFEM eliminates remeshing?


🎓

Exactly. XFEM augments standard FEM shape functions with enrichment functions:

$$\mathbf{u}^h(\mathbf{x}) = \sum_i N_i(\mathbf{x}) \, \mathbf{u}_i + \sum_j N_j(\mathbf{x}) \, H(\mathbf{x}) \, \mathbf{a}_j + \sum_k N_k(\mathbf{x}) \sum_{\alpha=1}^{4} F_\alpha(\mathbf{x}) \, \mathbf{b}_k^\alpha$$

The first term is standard FEM, the second uses the Heaviside step function $H$ to model the discontinuity across the crack, and the third introduces branch functions $F_\alpha$ for the crack-tip singularity. The crack geometry becomes independent of the mesh — you define it via level-set functions.


🧑‍🎓

Which commercial solvers support it?


🎓

Abaqus 6.9 and later have XFEM as a standard feature, making crack insertion and growth relatively straightforward. ANSYS Mechanical offers Smart Crack Growth. Open-source Code_Aster supports X-FEM too. But XFEM isn't a panacea — complex 3D crack branching and contact-aided fracture still pose challenges. For those cases, Cohesive Zone Model (CZM) or phase-field methods are worth considering.


Fracture Mechanics in CAE Practice

🧑‍🎓

When I run fracture mechanics analyses, what's the most critical thing to watch?


🎓

The mesh at the crack tip. For traditional FEM, create a concentric "spiderweb" mesh with singular quarter-point elements and collapse the crack-tip nodes to a single point. The innermost element size should be 1/20 of the crack length or smaller. Use 5–10 contours for J-integral computation and verify convergence on the outer contours.


🧑‍🎓

How do you analyze fatigue crack growth?


🎓

First, compute the stress intensity factor range $\Delta K = K_{\max} - K_{\min}$ at each load cycle. Then apply the Paris law to get the crack growth rate:

$$\frac{da}{dN} = C (\Delta K)^m$$

$C$ and $m$ are material constants (for steel, $m \approx 3$; aluminum, $m \approx 3$–$4$). Numerically integrate from initial crack size $a_0$ to critical size $a_c$ (where $K_{I\max} = K_{Ic}$) to find remaining life in cycles. The workflow is typically: compute K using Abaqus's *CONTOUR INTEGRAL, then feed results to tools like fe-safe for propagation analysis.


🧑‍🎓

I'm drowning in choices: LEFM, EPFM, XFEM, CZM… How do I pick?


🎓

Here's a quick guide:

Start by assessing plastic zone size and crack geometry. If unsure, begin with LEFM, then transition to J-integral if $r_p$ is too large. This stepwise approach works well.


Precise understanding of CAE terminology is the foundation of team communication. — NovaSolver Project supports engineering learning.

Share your fracture mechanics challenges

NovaSolver addresses real CAE engineer pain points — setup complexity, computational cost, result interpretation. Your practical experience can drive better tool development.

Contact Us (Coming Soon)
Article Evaluation
Thank you for your response!
Helpful
More
Detail
Report
Error
Helpful
0
More Detail
0
Report Error
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — Sitemap
View Profile