Crack Propagation & Critical Crack Length

The Crack Propagation & Critical Crack Length Simulator is an interactive Linear Elastic Fracture Mechanics (LEFM) lab. It computes the Mode I stress intensity factor K_I = Y · σ · √(πa), finds the critical crack length a_c at which fast brittle fracture occurs, and animates fatigue crack growth using the Paris law da/dN = C(ΔK)^m.

It targets Bachelor-level mechanical and materials-science students: design of safety-critical components, fracture toughness selection, fatigue life prediction, and damage-tolerance analysis.

The simulator opens in Simulate → Static mode with a Centre Crack in 4340 steel. The top canvas shows the cracked plate with a colour-mapped σ_yy stress field (with the r^−½ singularity around the crack tip), the crack itself, and the Irwin plastic zone. The chart below plots K_I vs crack length a with the K_IC horizontal line and the critical crack length marker.

Pick a geometry, pick a material (12 presets covering ductile metals, high-strength alloys, ceramics, and polymers), then drag the a, σ, and W sliders. The status badge turns SAFE → CAUTION → FRACTURE as K_I approaches and exceeds K_IC.

In Static sub-tab, the simulator evaluates K_I at the current crack length and compares it to K_IC. The 4 geometries available:

• Centre Crack — through-thickness 2a centre crack in a plate of width W. Y = sec(πa/W)^½ (Feddersen secant).
• Single Edge Crack — edge-through crack of length a in a plate of width W under tension. Y polynomial in a/W (Tada et al.).
• Double Edge Crack — symmetric cracks of length a at both edges.
• Semi-elliptical Surface Crack — a/c = 0.4, deepest-point K_I (Newman–Raju approximation).

Switch the Analysis sub-tab to Paris Law. Now the simulator integrates da/dN = C(ΔK)^m starting from an initial crack a₀ until a reaches the critical crack length a_c. Adjust Δσ (cyclic stress range) and stress ratio R. Press Play to animate the crack growing visibly in the canvas; the sparkline in the Learning Panel logs a vs N.

Material constants C and m are taken from the selected material’s database entry where available, otherwise generic ferritic-steel values are used and flagged.

Five categories: Fundamentals (Griffith, LEFM, K-field), Modes I/II/III (opening, sliding, tearing), Materials & K_IC (database with toughness ranges), Fatigue (Paris) (three regions of da/dN−ΔK, threshold ΔK_th), and Design & Standards (damage tolerance, ASTM E399, E647).

Practice generates random LEFM problems — compute K_I, solve for a_c, find cycles to failure from Paris law. Submit a numeric answer; click Show Solution for a full worked walk-through.

Quiz is a 5-question, mixed-format session (multiple choice + numeric) with star rating at the end.

The SI / Imperial toggle switches displayed units: stress (MPa↔ksi), length (mm↔in), toughness (MPa·√m ↔ ksi·√in). Internal computation always uses SI.

Right-click the canvas for Export PNG / Copy K_I / Reset. The action bar exposes CSV export of the K vs a curve and PNG export with watermark. Keyboard: arrow keys nudge a focused slider; Escape closes the calculation modal.

Geometry factors Y(a/W) are taken from standard handbook fits:

• Centre crack — Feddersen secant, Y = √(sec(π·a/W)); a is the half-crack length, valid for 2a/W < 1.
• Single edge — Tada/Paris/Irwin universal trigonometric form (accurate to 0.5% for all a/W).
• Double edge — Anderson Table 2.4 polynomial, F(α) = 1.12 + 0.43α − 4.79α² + 15.46α³ with α = a/W (W = full plate width).
• Surface (semi-elliptical, a/c = 0.4) — Newman-Raju with M₁ = 1.094, Q = 1.316 and a simple finite-width correction; assumes plate thickness ≈ width.

Fatigue (Paris) physics: Walker effective ΔK with γ = 0.5, ΔK_eff = ΔK/√(1−R) for R ≥ 0; for R < 0, the compressive half is suppressed (crack closure) and ΔK_eff = ΔK. Threshold scales as ΔK_th(R) = ΔK_th0·√(1−R).

Fast-fracture criterion: K_max = σ_max·Y·√(πa) ≥ K_IC, with σ_max = Δσ/(1−R). The critical crack length a_c is solved iteratively (bisection) because Y depends on a/W.

Out of scope (not modelled): Mode II/III mixed loading; stress-corrosion cracking; hold-time / dwell effects; plasticity-induced closure beyond simple Walker; environmental effects; residual stress; small-crack regime (a < ~10 grain diameters); FOREMAN ramp-up to K_IC; load-sequence effects (overload retardation).

• LEFM small-scale yielding requirement: the plastic zone r_p must be small compared to a and the remaining ligament (W−a). The simulator warns when this assumption is violated.
• Plane stress vs plane strain: thin plates yield in plane stress (larger r_p); thick sections are in plane strain (smaller r_p, lower apparent toughness). K_IC in the database is the plane-strain value.
• Damage tolerance design: always assume an initial flaw of detectable size (often 0.5–1 mm by NDT) and ensure cycle life exceeds inspection interval × safety factor.
• Brittle materials (ceramics, glass) have very low K_IC (< 5 MPa·√m) — even hairline cracks are critical.
• Threshold ΔK_th — below this value, fatigue cracks do not grow. Useful for infinite-life design.