6.1 Introduction: Beyond Linear Elastic Fracture Mechanics

Chapter 5 established LEFM for situations where plastic deformation at the crack tip is confined to a small zone. When the plastic zone is large relative to crack size, ligament, or specimen thickness, LEFM loses validity. This occurs in:

• Tough, ductile metals (low-strength steels, aluminum) where extensive blunting precedes fracture.
• Thin sheets where plane-stress conditions allow large-scale yielding.
• Crack growth under cyclic loading where the plastic zone grows with each cycle (ratcheting).
• High-temperature creep where time-dependent deformation enlarges the process zone.

Elastic-Plastic Fracture Mechanics (EPFM) extends fracture analysis through:

1. Paris law — empirical crack growth under cyclic loading (still in the LEFM $K$ regime).
2. $J$-integral — energy-based parameter valid for nonlinear materials.
3. Crack Tip Opening Displacement (CTOD) — physical measure of crack blunting.

This chapter covers Paris law in depth (the primary tool for damage-tolerant design) and introduces $J$ and CTOD for large-scale yielding.

---

6.2 Fatigue Crack Growth: The Three Stages

When a cyclically loaded structure contains a crack, the crack grows incrementally per cycle. On a log–log plot of crack growth rate $da/dN$ vs stress intensity range $\Delta K$, three distinct regions appear:

6.2.1 Region I — Threshold

Below a threshold $\Delta K_{th}$, cracks do not propagate (or grow at imperceptible rates). $\Delta K_{th}$ depends on $R$-ratio, environment, and material microstructure.

6.2.2 Region II — Paris (steady-state) regime

Power-law relationship between growth rate and $\Delta K$:

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

This is the Paris law (1963), valid over several orders of magnitude in $da/dN$ for most engineering alloys in air.

6.2.3 Region III — Fast fracture

As $K_{\max}$ approaches $K_{Ic}$, growth accelerates rapidly to unstable fracture. This region is governed by $K_{\max}$, not $\Delta K$ alone.

---

6.3 Paris Law: Parameters and Physical Meaning

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

where:

• $da/dN$ = crack growth per cycle (m/cycle or mm/cycle).
• $\Delta K = K_{\max} - K_{\min}$ = stress intensity range.
• $C$ = material/environment constant (units depend on $m$).
• $m$ = Paris exponent (typically 2–4 for metals).

6.3.1 Typical Paris parameters

Units caution: $C$ depends on the units of $da/dN$ and $\Delta K$. Always verify consistent SI or US customary units before integration.

6.3.2 Physical interpretation

The Paris exponent $m$ reflects the sensitivity of crack growth to stress intensity:

• $m \approx 2$: mild sensitivity (some aluminum alloys).
• $m \approx 3$: typical for steels.
• $m \approx 4$–$6$: high sensitivity (brittle materials, corrosive environments).

Higher $m$ means small increases in $\Delta K$ cause dramatically faster growth — demanding tighter inspection intervals.

---

6.4 Integration of Paris Law for Remaining Life

Given an initial crack size $a_i$ and a final (critical) size $a_f$, integrate Paris law to find the number of cycles to grow from $a_i$ to $a_f$:

$$\boxed{N = \int_{a_i}^{a_f} \frac{da}{C(\Delta K)^m}}$$

6.4.1 Constant $Y$ and constant $\Delta\sigma$

If $\Delta K = Y \Delta\sigma \sqrt{\pi a}$ and $Y$, $\Delta\sigma$ are constant:

$$N = \int_{a_i}^{a_f} \frac{da}{C(Y\Delta\sigma\sqrt{\pi})^m a^{m/2}}$$

For $m \ne 2$:

$$\boxed{N = \frac{a_f^{1-m/2} - a_i^{1-m/2}}{C(Y\Delta\sigma\sqrt{\pi})^m (1 - m/2)}}$$

For $m = 2$ (special case):

$$N = \frac{1}{C(Y\Delta\sigma\sqrt{\pi})^2} \ln\frac{a_f}{a_i}$$

6.4.2 Variable $Y(a)$

When $Y$ changes with crack length (finite width), integrate numerically (Simpson's rule, or the Live Solver in this chapter).

---

6.5 $R$-Ratio Effects on Crack Growth

The load ratio $R = K_{\min}/K_{\max} = \sigma_{\min}/\sigma_{\max}$ strongly affects crack growth:

6.5.1 Threshold dependence on $R$

$\Delta K_{th}$ decreases as $R$ increases (more tensile mean load). Empirical form (Walker, Forman):

$$\Delta K_{th}(R) = \Delta K_{th,0}(1-R)^{\gamma}$$

where $\gamma \approx 1$–$2$ and $\Delta K_{th,0}$ is the threshold at $R = 0$.

6.5.2 Growth rate dependence on $R$

At constant $\Delta K$, higher $R$ (higher $K_{\max}$) increases $da/dN$ because:

• Crack faces are less likely to close (reduced roughness-induced closure).
• $K_{\max}$ approaches $K_{Ic}$, entering Region III acceleration.

6.5.3 Effective stress intensity range

Elber's crack closure model introduces an effective $\Delta K$:

$$\Delta K_{\mathrm{eff}} = K_{\max} - K_{\mathrm{op}}$$

where $K_{\mathrm{op}}$ is the stress intensity at which the crack faces first open during the loading cycle. For $R \ge 0$: $\Delta K_{\mathrm{eff}} \approx \Delta K$. For $R < 0$ (compression in part of cycle): $\Delta K_{\mathrm{eff}} < \Delta K$.

---

6.6 Threshold and Crack Growth Threshold Design

Below $\Delta K_{th}$, cracks do not grow (or grow at $< 10^{-10}$ m/cycle). In damage-tolerant design:

1. Compute $\Delta K$ at the inspection limit crack size $a_i$.
2. If $\Delta K < \Delta K_{th}$: the crack is non-propagating — infinite life (in that environment).
3. If $\Delta K > \Delta K_{th}$: integrate Paris law for remaining cycles.

6.6.1 Non-propagating crack concept

Small cracks at notches may arrest if $\Delta K$ falls below threshold as the crack grows into the compressive residual stress field (shot-peened surfaces). This is the basis of damage tolerance with arrest features.

---

6.7 Introduction to the $J$-Integral

When LEFM is invalid (large plastic zone), the $J$-integral (Rice, 1968) provides an energy-based fracture parameter that remains path-independent for nonlinear elastic (or deformation theory plastic) materials:

$$\boxed{J = \oint_\Gamma \left(W\, dy - T_i \frac{\partial u_i}{\partial x}\, ds\right)}$$

where:

• $W$ = strain energy density.
• $T_i$ = traction vector on the contour $\Gamma$.
• $u_i$ = displacement field.
• $ds$ = arc length element.

6.7.1 Physical meaning

$J$ equals the energy release rate $G$ for linear elasticity and generalizes it to nonlinear materials. At initiation of crack growth:

$$J \ge J_{Ic}$$

where $J_{Ic}$ is the critical $J$ value (fracture toughness in EPFM).

6.7.2 Relation to $K_{Ic}$

In the linear elastic limit:

$$J = \frac{K_I^2}{E'} \quad\Rightarrow\quad K_{Jc} = \sqrt{J_{Ic} E'}$$

For a material where LEFM is borderline valid, $J_{Ic}$ can be measured on a smaller specimen than $K_{Ic}$.

6.7.3 $J$-resistance curve ($J_R$)

For ductile tearing, $J$ increases with crack extension $\Delta a$ (tearing resistance). The $J_R$ curve characterizes stable crack growth before final instability — essential for pressure vessel and pipeline design.

---

6.8 Crack Tip Opening Displacement (CTOD)

CTOD ($\delta$) is the physical opening displacement at the crack tip, measured at the intersection of the crack faces with a $90°$ vertex (or at $0.2$ mm offset from the tip).

6.8.1 Irwin relation (LEFM limit)

$$\delta = \frac{K_I^2}{E' S_y} = \frac{G}{S_y}$$

6.8.2 Wells' CTOD criterion

For large-scale yielding, fracture occurs when:

$$\delta \ge \delta_c$$

where $\delta_c$ is the critical CTOD (material property, typically 0.1–0.5 mm for structural steels).

6.8.3 Dugdale model

The strip-yield model gives:

$$\delta = \frac{8 S_y a}{\pi E} \ln\sec\left(\frac{\pi K_I}{2 S_y}\right)$$

For small $K_I/S_y$: $\delta \approx K_I^2/(E S_y)$ (Irwin relation).

6.8.4 Design use

CTOD is the primary fracture criterion in offshore structural steel design (BS 7910, API 579). It directly measures the blunting that precedes ductile fracture.

---

6.9 Damage-Tolerant Design Procedure

A complete damage-tolerant design cycle:

1. Assume initial flaw $a_i$ = NDT detection limit (e.g., 1.27 mm for dye penetrant).
2. Compute $\Delta K$ at $a_i$ under service loading spectrum.
3. Check threshold: if $\Delta K < \Delta K_{th}$, non-propagating — done.
4. Integrate Paris law from $a_i$ to critical $a_f$ (where $K_{\max} \to K_{Ic}$).
5. Compare $N$ to required service life with safety factor (typically 2× on life).
6. Set inspection interval: inspect at $N/2$ or when crack reaches $a_i$ again after repair.
7. If $N$ insufficient: reduce $\Delta\sigma$ (redesign), improve material ($K_{Ic}$, lower $C$), or add arrest features.

---

6.10 Environmental and Spectrum Effects

6.10.1 Corrosion fatigue

In corrosive environments (seawater, H₂S), $C$ increases by 10–100× and $\Delta K_{th}$ decreases dramatically. Use environment-specific Paris data.

6.10.2 Variable amplitude loading

Real spectra (aircraft gusts, vehicle road loads) require cycle counting (rainflow) and Miner's rule on crack growth increments, or sophisticated retardation models (Willenborg, Wheeler).

6.10.3 Temperature

Crack growth rates generally increase with temperature. Use temperature-corrected $C$ and $m$ from databases.

---

6.11 Summary

Design takeaway: Paris law integration is the primary tool for setting inspection intervals in aerospace, pressure vessels, and rotating machinery. Always use environment-specific $C$, $m$, and $\Delta K_{th}$. When LEFM validity is marginal, supplement with $J$ or CTOD analysis.