Paris-Erdogan Fatigue Crack Growth & Fracture Calculator

Comprehensive Linear Elastic Fracture Mechanics (LEFM) Fatigue Guide

Fatigue failure accounts for over 80% of all catastrophic structural collapses in engineering service, including aircraft fuselages, railway axles, offshore drilling rigs, pressure piping, and bridge trusses. Linear Elastic Fracture Mechanics (LEFM) treats structural safety not by assuming pristine, uncracked metal, but by acknowledging the inevitable presence of microscopic manufacturing flaws and calculating the exact cyclic rate at which these cracks grow toward critical rupture under service fatigue.

1. The Stress Intensity Factor ($K$) & Fracture Toughness ($K_{Ic}$)

In linear elasticity, George R. Irwin (1957) proved that the asymptotic stress singularity near a sharp crack tip in an isotropic medium is uniquely scaled by the **Stress Intensity Factor ($K$)** with units of $text{MPa}sqrt{text{m}}$:

For an opening-mode (Mode I) crack of length $a$ in a tensile stress field $sigma$:

where $Y$ is a dimensionless geometry correction factor dependent on component boundary conditions and crack aspect ratio ($a/W$):

• Center-Cracked Plate in Tension (CCT): $Y = sqrt{secleft(frac{pi a}{2W}right)} approx 1.0$ for wide plates.
• Single Edge Notch Tension (SENT): $Y approx 1.12 - 0.231(a/W) + 10.55(a/W)^2 - 21.72(a/W)^3 + 30.39(a/W)^4 approx 1.12$ for shallow edge cracks.
• Semi-Elliptical Surface Flaw: $Y approx frac{1.12}{Phi} approx 0.70 dots 1.10$.

Catastrophic, rapid, brittle fracture occurs instantaneously when the maximum stress intensity reaches the material's plane-strain **fracture toughness ($K_{Ic}$)**: $K_{max} = K_{Ic}$.

2. The Three Regimes of Fatigue Crack Growth

When cyclic tensile loading is applied, the crack length $a$ advances by an increment $Delta a$ each cycle. Plotting the crack propagation rate $frac{da}{dN}$ against the cyclic stress intensity factor range $Delta K = K_{max} - K_{min} = Y Deltasigma sqrt{pi a}$ on logarithmic axes reveals three distinct behavioral regimes:

1. Regime I (Threshold Near $Delta K_{th}$): At extremely low stress intensities ($frac{da}{dN} < 10^{-9}text{ m/cycle}$), crack growth is governed by crystallographic microstructural barriers (grain boundaries). If $Delta K < Delta K_{th}$, crack propagation arrests entirely, defining the infinite-life threshold.
2. Regime II (Paris Steady-State Striation Growth): The intermediate regime where crack extension advances via continuum plastic blunting and resharpening, leaving microscopic fatigue striations on the fracture surface. This power-law region is governed by the classical **Paris-Erdogan Law**.
3. Regime III (High $Delta K$ Acceleration Toward $K_{Ic}$): As $K_{max}$ approaches $K_{Ic}$, local static fracture mechanisms (microvoid coalescence, cleavage bursts) intermingle with fatigue striations, causing acceleration toward infinite growth rate ($da/dN to infty$).

3. Paris-Erdogan Equation & Mean Stress Ratio ($R$)

In 1963, Paul Paris and Fazil Erdogan established the empirical power law governing Regime II:

where $C$ and $m$ are empirical constants determined from ASTM E647 compact tension (CT) testing. For standard structural steels, $m approx 3.0$; for aerospace aluminum alloys, $m approx 3.2 dots 4.0$.

Mean Stress Ratio ($R$) Influence & The Walker Model

The stress ratio $R = sigma_{min} / sigma_{max}$ alters the effective crack tip opening stress due to plasticity-induced crack closure (Elber's closure mechanism). Higher mean stress ($R > 0.5$) elevates growth rates. The **Walker Equation** accounts for this shift:

where $gamma approx 0.5 dots 0.7$ is the Walker material sensitivity exponent.

4. Analytical Lifespan Integration ($N_f$) and Critical Crack Size ($a_c$)

Rearranging the Paris law and integrating from an initial detectable flaw size $a_0$ to the critical catastrophic failure crack size $a_c$:

Assuming geometry factor $Y$ remains approximately constant during crack propagation, the integral yields the exact closed-form solution for $m ne 2$:

The **critical crack length ($a_c$)** is governed by fracture toughness:

Notice that because $m approx 3$, the term $a_0^{1 - m/2} = a_0^{-0.5} = frac{1}{sqrt{a_0}}$ heavily dominates the numerator compared to $a_c^{-0.5} = frac{1}{sqrt{a_c}}$. Over **85% to 95% of a component's total fatigue life is consumed while the crack is tiny** (growing from $a_0$ to $2 a_0$), leaving very little time once the crack becomes large!

5. Damage-Tolerant Inspection Intervals

Under modern aerospace and nuclear damage-tolerance design paradigms (e.g., USAF Damage Tolerance Handbook, ASME Section XI), components are designed to tolerate flaws throughout service. Periodic Non-Destructive Testing (NDT) inspections (eddy current, ultrasonic, X-ray) are mandated at intervals:

This ensures at least two opportunities to detect and repair an advancing flaw before catastrophic fracture.

6. Step-by-Step Worked Fatigue Crack Growth Calculation

Assess an edge-notched aircraft landing gear component forged from AISI 4340 high-strength steel:

• Initial NDT detectable flaw size $a_0 = 2.0text{ mm} = 0.002text{ m}$
• Edge crack geometry factor $Y = 1.12$
• Cyclic stress: $sigma_{max} = 220text{ MPa}$, $sigma_{min} = 20text{ MPa} implies Deltasigma = 200text{ MPa}$, $R = 0.091$
• Fracture toughness $K_{Ic} = 65text{ MPa}sqrt{text{m}}$
• Paris parameters: $C = 1.5 times 10^{-11}text{ m/(cycle}cdot(text{MPa}sqrt{text{m}})^3)$, $m = 3.0$

Step 1: Calculate Critical Crack Size ($a_c$)

Step 2: Calculate Initial Stress Intensity Range ($Delta K_0$)

Because $Delta K_0 = 17.76text{ MPa}sqrt{text{m}} > Delta K_{th} approx 3.5text{ MPa}sqrt{text{m}}$, the flaw is active and propagates immediately.

Step 3: Initial Crack Growth Rate

Step 4: Integrate Total Cycles to Failure ($N_f$)

The component sustains approximately $33,300$ flights before catastrophic separation. The maximum allowable inspection interval is $N_f / 2 approx 16,600$ flight cycles.

7. Standard Paris Law Material Constants Reference Library