🩹 Microcracks (Paris Law)
Interactive fatigue crack growth simulator using the Paris law da/dN = C(ΔK)^m. Watch microcracks propagate cycle by cycle until fracture toughness K_IC is reached.
About this simulation
This simulation models fatigue crack growth in a metal specimen under cyclic loading, using the Paris law da/dN = C·(ΔK)^m, the standard engineering equation for stable crack propagation. Each cycle the crack tip stress-intensity range ΔK = Y·Δσ·√(πa) is recomputed from the current crack length a, and the crack advances by da/dN millimetres per cycle. Growth continues until the peak stress intensity K_max = Y·σ_max·√(πa) reaches the material's fracture toughness K_IC, at which point the specimen fractures catastrophically in a single frame rather than continuing to grow slowly.
🔬 What it shows
A rectangular specimen pulled between σ_max and σ_min every cycle, with a crack growing from an initial length a₀ down through the material. The live plot tracks crack length a against cycles N on a log scale, marking the critical crack length a_c where K_max = K_IC, and the 3D view splits the specimen and shows a fracture surface once failure occurs.
🎮 How to use
Choose a material preset from the Material dropdown (Steel A36, Aluminium 7075, Titanium Ti-6Al-4V, Cast Iron) or dial in your own C, m and K_IC values. Adjust Initial crack a₀, σ_max, σ_min, Specimen width W and Y geometry factor sliders, set Cycles per frame to speed up or slow down the simulation, and use Pause/Play and Reset to control the run.
💡 Did you know?
Because da/dN scales with ΔK raised to the power m, doubling the applied stress range Δσ does not just double the crack growth rate — for a typical steel with m≈3 it multiplies it by roughly 2³ = 8, which is why small increases in cyclic load can slash a component's fatigue life dramatically.
Frequently asked questions
What is the Paris law and what do C and m mean?
The Paris law is an empirical equation, da/dN = C·(ΔK)^m, describing how fast a fatigue crack grows per load cycle once it is in the stable propagation regime. C and m are material constants fitted from test data: C sets the overall growth rate and m controls how sensitive that rate is to the stress-intensity range ΔK. In this simulator C is entered in units of ×10⁻¹¹ and m typically ranges from about 2.9 for titanium to 4.5 for cast iron.
How is the stress-intensity factor calculated?
The simulation uses the standard linear-elastic fracture mechanics formula K = Y·σ·√(πa), where σ is the applied stress, a is the current crack length and Y is a dimensionless geometry factor that accounts for specimen shape (adjustable via the Y geometry factor slider, default 1.12 for an edge crack). ΔK uses the stress range Δσ = σ_max − σ_min, while K_max uses σ_max directly; failure is triggered as soon as K_max reaches K_IC.
What determines when the specimen fractures?
Fracture occurs when the peak stress intensity K_max = Y·σ_max·√(πa) reaches the material's fracture toughness K_IC, entered in MPa·√m via the K_IC slider. As the crack a grows, K_max rises even though the applied stresses stay constant, so every material and load combination has a critical crack length a_c (shown as a dashed line on the plot) beyond which the crack becomes unstable and the specimen snaps rather than growing slowly.
Is this a full physical crack simulation or a simplified model?
It is a simplified engineering model, not a finite-element or full physical simulation. The crack path itself is a randomly jittered line for visual effect, and the specimen geometry is a basic rectangular bar rather than a real fracture-test coupon. The growth physics, however, is the genuine Paris law integrated numerically cycle by cycle with a real stress-intensity formula, so the trends in crack length, growth rate and cycles-to-failure reflect real fatigue behaviour.
Why do the different material presets fail after such different numbers of cycles?
Each preset (Steel A36, Aluminium 7075, Titanium Ti-6Al-4V, Cast Iron) has its own C, m and K_IC values drawn from typical published fatigue data, and the estimated Cycles to failure statistic is computed by integrating the Paris law from the current crack length to the critical length a_c. Materials with low K_IC (like cast iron at 20 MPa·√m) reach their critical crack length sooner, while a high Paris exponent m makes growth accelerate much faster once cracking begins, so small differences in material properties translate into large differences in fatigue life.
Cyclic stress slowly grows microscopic cracks per cycle by da/dN = C·(ΔK)^m. Watch the crack creep along its meandering path — then snap catastrophically when K_max reaches the fracture toughness K_IC.
3D · Three.js / WebGL renderer · 60 FPS target · runs fully client-side, no install