PracticeVideosBlogFlashcards

5.2 Static Failure Theories for Brittle Materials & Fracture Mechanics

Key Takeaways

  • Brittle materials exhibit percent elongation %EL<5%\%EL < 5\% and significant strength asymmetry (Suc≫SutS_{uc} \gg S_{ut}, often 3 to 4 times) because internal micro-flaws open under tension but close under compression.

  • The Maximum Normal Stress (MNS) theory predicts fracture when σ1≥Sut\sigma_1 \ge S_{ut} or σ3≤−Suc\sigma_3 \le -S_{uc}; while valid in biaxial tension (Quadrant 1), it is dangerously non-conservative in Quadrant 4 where compressive stress accelerates tensile fracture.

  • The Modified Mohr (MM / Dowling) theory eliminates the unwarranted conservatism of the Brittle Coulomb-Mohr (BCM) theory in Quadrant 4 when ∣σB/σA∣≤1|\sigma_B/\sigma_A| \le 1 (n=Sut/σAn = S_{ut}/\sigma_A), accurately tracking empirical biaxial test data.

  • Linear Elastic Fracture Mechanics (LEFM) governs crack propagation via the stress intensity factor KI=YσπaK_I = Y \sigma \sqrt{\pi a}, predicting brittle fracture when KI≥KIcK_I \ge K_{Ic} where KIcK_{Ic} is the plane-strain toughness measured on ASTM E399 specimens sized by the 2.5(KIc/Sy)22.5(K_{Ic}/S_y)^2 rule.

  • The Leak-Before-Break (LBB) criterion requires that the critical crack size exceed the wall thickness (ac>ta_c > t), ensuring through-wall penetration produces detectable fluid leakage before catastrophic fast fracture can initiate.

Last updated: October 2026

5.2 Static Failure Theories for Brittle Materials & Fracture Mechanics

Quick Summary: Brittle materials sustain negligible plastic deformation prior to catastrophic fracture (%EL<5%\%EL < 5\%) and exhibit markedly asymmetric strength, with ultimate compressive strengths exceeding tensile capacities by factors of three to four. Classical stress concentration formulations break down at sharp flaws where theoretical notch radii approach zero. This section covers the governing brittle failure criteria—Maximum Normal Stress, Brittle Coulomb-Mohr, and the Modified Mohr (Dowling) theory—and introduces Linear Elastic Fracture Mechanics (LEFM), including stress intensity (KIK_I), ASTM E399 plane strain fracture toughness (KIcK_{Ic}), valid thickness constraints, NDE crack detection limits, and the Leak-Before-Break (LBB) design criterion.


1. Mechanical Behavior and Asymmetric Strength of Brittle Materials

Engineering materials exhibiting percent elongation at fracture %EL<5%\%EL < 5\% are classified as brittle. Common engineering brittle materials include gray cast iron, white cast iron, hardened tool steels (quenched without tempering), structural ceramics (silicon nitride, alumina, zirconia), glass, concrete, and powder metallurgy components.

Cleavage Fracture vs. Plastic Slip

Unlike ductile metals that deform via dislocation glide along close-packed slip planes, brittle solids undergo cleavage fracture—the rapid, direct separation of atomic crystallographic planes normal to the maximum tensile stress direction. Because cleavage requires very little energy dissipation compared to plastic shear slip, brittle fracture is sudden, catastrophic, and unaccompanied by gross dimensional warning (no necking or yield plateau).

The Physical Origin of Compressive-to-Tensile Strength Asymmetry

The defining mechanical characteristic of brittle materials is their pronounced strength asymmetry:

Suc≫SutS_{uc} \gg S_{ut}

For standard grades of gray cast iron (ASTM A48 Class 20 to Class 60), the ultimate compressive strength is typically three to four times greater than the ultimate tensile strength (Suc≈3–4SutS_{uc} \approx 3\text{--}4 S_{ut}):

Material SpecificationTensile Strength SutS_{ut}Compressive Strength SucS_{uc}Ratio Suc/SutS_{uc} / S_{ut}
Gray Cast Iron Class 20150 MPa150\text{ MPa} (22 ksi22\text{ ksi})570 MPa570\text{ MPa} (83 ksi83\text{ ksi})3.83.8
Gray Cast Iron Class 30214 MPa214\text{ MPa} (31 ksi31\text{ ksi})750 MPa750\text{ MPa} (109 ksi109\text{ ksi})3.53.5
Gray Cast Iron Class 40293 MPa293\text{ MPa} (42.5 ksi42.5\text{ ksi})965 MPa965\text{ MPa} (140 ksi140\text{ ksi})3.33.3
Structural Concrete3–5 MPa3\text{--}5\text{ MPa} (0.4–0.7 ksi0.4\text{--}0.7\text{ ksi})30–40 MPa30\text{--}40\text{ MPa} (4.3–5.8 ksi4.3\text{--}5.8\text{ ksi})8–108\text{--}10
Sintered Alumina Ceramic250 MPa250\text{ MPa} (36 ksi36\text{ ksi})2500 MPa2500\text{ MPa} (360 ksi360\text{ ksi})10.010.0

This extreme asymmetry arises from internal microstructural defects. Gray cast iron contains an interconnected network of sharp, brittle graphite flakes dispersed throughout a pearlitic or ferritic iron matrix:

  1. Under Tensile Stress: The sharp tip radii of graphite flakes act as severe internal stress concentrators. Tensile stresses open micro-cracks at flake tips, which link together and propagate across the cross section at low nominal tensile loads (SutS_{ut}).
  2. Under Compressive Stress: Compressive forces push the surfaces of the graphite flakes and micro-cracks into mechanical contact. Flaws are squeezed shut, allowing the solid metallic iron matrix to transfer compressive load across the entire gross cross-sectional area. Failure can only occur when the iron matrix itself crushes or shears along oblique planes at very high stress levels (SucS_{uc}).

Fracture Morphology in Tension vs. Torsion

Because brittle fracture is governed by normal tensile cleavage:

  • Uniaxial Tension: Fractures along a flat plane oriented precisely perpendicular (90°) to the tensile loading axis.
  • Pure Torsion: Under pure shear stress τ\tau, Mohr's circle indicates that the principal stresses are σ1=+τ\sigma_1 = +\tau and σ3=−τ\sigma_3 = -\tau, acting on planes inclined at 45° to the shaft longitudinal axis. Consequently, a brittle cast iron shaft tested in pure torsion fractures along a characteristic 45° helical cleavage plane.

2. Maximum Normal Stress (MNS) Theory

Formulated by Scottish engineer William Rankine, the Maximum Normal Stress (MNS) theory is the simplest static failure model for brittle materials.

The MNS Postulate

Brittle fracture occurs when the maximum principal tensile stress in a multiaxial stress state reaches the ultimate tensile strength SutS_{ut}, or when the minimum principal stress reaches the ultimate compressive strength −Suc-S_{uc}.

Mathematical Formulation

Ordering principal stresses algebraically such that σ1≥σ2≥σ3\sigma_1 \ge \sigma_2 \ge \sigma_3, the MNS failure criterion states that fracture occurs if:

σ1≥Sutorσ3≤−Suc\sigma_1 \ge S_{ut} \quad \text{or} \quad \sigma_3 \le -S_{uc}

The design factor of safety (nn) is evaluated separately for tension and compression:

nt=Sutσ1(if σ1>0),nc=Suc∣σ3∣(if σ3<0)n_t = \frac{S_{ut}}{\sigma_1} \quad (\text{if } \sigma_1 > 0), \qquad n_c = \frac{S_{uc}}{|\sigma_3|} \quad (\text{if } \sigma_3 < 0)

The governing factor of safety is n=min⁡(nt,nc)n = \min(n_t, n_c).

Evaluation and Severe Limitations of MNS

In plane stress (σ3=0\sigma_3 = 0), the MNS criterion plots as a simple rectangle on the σ1−σ2\sigma_1 - \sigma_2 plane bounded by σ1=Sut\sigma_1 = S_{ut}, σ1=−Suc\sigma_1 = -S_{uc}, σ2=Sut\sigma_2 = S_{ut}, and σ2=−Suc\sigma_2 = -S_{uc}:

  • Quadrant 1 (Biaxial Tension: σ1>0,σ2>0\sigma_1 > 0, \sigma_2 > 0): MNS exhibits excellent agreement with experimental fracture data. Superimposed transverse tensile stress does not significantly alter the tensile fracture capacity.
  • Quadrant 4 (Combined Tension and Compression: σ1>0,σ2<0\sigma_1 > 0, \sigma_2 < 0): Extensive experimental testing (notably by Coffin, Grassi, and Cornet on thin-walled gray iron tubes) reveals that MNS is dangerously non-conservative. When compressive stress is applied simultaneously with tensile stress, the tensile stress required to trigger brittle fracture drops well below SutS_{ut}. Designing with MNS in Quadrants 2 and 4 can lead to unexpected catastrophic rupture.

3. Coulomb-Mohr and Modified Mohr (MM) Theories

To account for strength asymmetry (Suc>SutS_{uc} > S_{ut}) and overcome the non-conservatism of MNS in mixed tension-compression stress states, two Mohr-based failure criteria are used in machine design.

Brittle Coulomb-Mohr (BCM) Theory

The Brittle Coulomb-Mohr (BCM) theory (also known as the internal friction theory) constructs a straight-line failure envelope tangent to the uniaxial tensile test circle (diameter SutS_{ut}) and the uniaxial compressive test circle (diameter SucS_{uc}) on a Mohr's stress plane.

For 2D plane stress with ordered in-plane principal stresses σA≥σB\sigma_A \ge \sigma_B:

  1. Quadrant 1 (σA≥σB≥0\sigma_A \ge \sigma_B \ge 0): n=SutσAn = \frac{S_{ut}}{\sigma_A}
  2. Quadrants 2 & 4 (σA≥0≥σB\sigma_A \ge 0 \ge \sigma_B): σASut−σBSuc=1n  ⟹  n=1σASut−σBSuc\frac{\sigma_A}{S_{ut}} - \frac{\sigma_B}{S_{uc}} = \frac{1}{n} \implies n = \frac{1}{\frac{\sigma_A}{S_{ut}} - \frac{\sigma_B}{S_{uc}}}
  3. Quadrant 3 (0≥σA≥σB0 \ge \sigma_A \ge \sigma_B): n=−SucσB=Suc∣σB∣n = -\frac{S_{uc}}{\sigma_B} = \frac{S_{uc}}{|\sigma_B|}

While BCM prevents the dangerous non-conservatism of MNS in Quadrant 4, extensive testing shows that BCM is excessively conservative when compressive stress is small (∣σB∣≤Sut|\sigma_B| \le S_{ut}).

Modified Mohr (MM) Theory (Dowling Equations)

The Modified Mohr (MM) theory—standardized by N.E. Dowling and adopted in Shigley's Mechanical Engineering Design—modifies the BCM envelope to match experimental biaxial test data with exceptional precision. It recognizes that a minor compressive stress does not degrade the tensile fracture capacity until the compressive stress exceeds SutS_{ut}.

For ordered in-plane principal stresses σA≥σB\sigma_A \ge \sigma_B:

                               σB
                                ^
                                |   Quadrant 1
                            Sut + - - - - - - - .
                                |               |
                                |   Modified    |
                                |     Mohr      |
               -Suc             |    Envelope   |       Sut
     <-----------+--------------+---------------+--------+------> σA
                 |              |               |
                 |              |               |
                 |              |   Quadrant 4  |
                 |              * - - - - - - - '
                 |             /   (Sut, -Sut)
                 |            /
                 |           /
                 |          /
            -Suc . - - - - '
                   Quadrant 3

The factor of safety (nn) is evaluated using Dowling's four distinct piecewise conditions:

Condition 1: Biaxial Tension (Quadrant 1: σA≥σB≥0\sigma_A \ge \sigma_B \ge 0)

Both principal stresses are tensile. Failure is governed strictly by tensile cleavage:

n=SutσAn = \frac{S_{ut}}{\sigma_A}

Condition 2: Tension-Dominant Mixed Stress (Quadrant 4: σA≥0≥σB\sigma_A \ge 0 \ge \sigma_B with ∣σB/σA∣≤1|\sigma_B / \sigma_A| \le 1)

Here, the magnitude of the compressive stress is less than or equal to the tensile stress (meaning σB≥−σA\sigma_B \ge -\sigma_A, which keeps the stress state above the −45∘-45^\circ diagonal line). Experimental data confirm that tensile fracture remains completely uninhibited by the transverse compression:

n=SutσAn = \frac{S_{ut}}{\sigma_A}

Condition 3: Compression-Dominant Mixed Stress (Quadrant 4: σA≥0≥σB\sigma_A \ge 0 \ge \sigma_B with ∣σB/σA∣>1|\sigma_B / \sigma_A| > 1)

When compressive stress exceeds tensile stress (∣σB∣>σA|\sigma_B| > \sigma_A), the compressive stress significantly accelerates crack propagation. Dowling's linear interaction equation governs:

(Suc−Sut)σASucSut−σBSuc=1n\frac{(S_{uc} - S_{ut})\sigma_A}{S_{uc} S_{ut}} - \frac{\sigma_B}{S_{uc}} = \frac{1}{n}

Solving explicitly for the factor of safety nn:

n=SucSut(Suc−Sut)σA−SutσBn = \frac{S_{uc} S_{ut}}{(S_{uc} - S_{ut})\sigma_A - S_{ut} \sigma_B}

Condition 4: Biaxial Compression (Quadrant 3: 0≥σA≥σB0 \ge \sigma_A \ge \sigma_B)

Both principal stresses are compressive. Failure is governed strictly by compressive crushing or shear splitting along the minimum principal stress axis:

n=−SucσB=Suc∣σB∣n = -\frac{S_{uc}}{\sigma_B} = \frac{S_{uc}}{|\sigma_B|}

Quadrant-by-Quadrant Theory Comparison

Stress QuadrantStress StateMaximum Normal Stress (MNS)Brittle Coulomb-Mohr (BCM)Modified Mohr (MM)
Quadrant 1σA≥σB≥0\sigma_A \ge \sigma_B \ge 0n=SutσAn = \frac{S_{ut}}{\sigma_A}n=SutσAn = \frac{S_{ut}}{\sigma_A}n=SutσAn = \frac{S_{ut}}{\sigma_A}
Quadrant 4 (∣σB/σA∣≤1\lvert\sigma_B/\sigma_A\rvert \le 1)σA>0>σB\sigma_A > 0 > \sigma_Bn=SutσAn = \frac{S_{ut}}{\sigma_A} (Dangerously Non-conservative)n=1σASut−σBSucn = \frac{1}{\frac{\sigma_A}{S_{ut}} - \frac{\sigma_B}{S_{uc}}} (Overly Conservative)n=SutσAn = \frac{S_{ut}}{\sigma_A} (Most Accurate)
Quadrant 4 (∣σB/σA∣>1\lvert\sigma_B/\sigma_A\rvert > 1)σA>0>σB\sigma_A > 0 > \sigma_Bn=min⁡(SutσA,Suc∣σB∣)n = \min(\frac{S_{ut}}{\sigma_A}, \frac{S_{uc}}{\lvert\sigma_B\rvert})n=1σASut−σBSucn = \frac{1}{\frac{\sigma_A}{S_{ut}} - \frac{\sigma_B}{S_{uc}}}n=SucSut(Suc−Sut)σA−SutσBn = \frac{S_{uc} S_{ut}}{(S_{uc} - S_{ut})\sigma_A - S_{ut}\sigma_B}
Quadrant 30≥σA≥σB0 \ge \sigma_A \ge \sigma_Bn=Suc∣σB∣n = \frac{S_{uc}}{\lvert\sigma_B\rvert}n=Suc∣σB∣n = \frac{S_{uc}}{\lvert\sigma_B\rvert}n=Suc∣σB∣n = \frac{S_{uc}}{\lvert\sigma_B\rvert}

4. Linear Elastic Fracture Mechanics (LEFM)

The Breakdown of Classical Stress Concentration

Classical elasticity analyzes geometric notches (fillets, holes, keyways) using stress concentration factors (Kt=σmax⁡/σnomK_t = \sigma_{\max} / \sigma_{\text{nom}}). For an elliptical hole with major axis 2a2a and notch root radius ρ\rho oriented perpendicular to a tensile stress, the Inglis solution yields:

σmax⁡=σ(1+2aρ)\sigma_{\max} = \sigma\left(1 + 2\sqrt{\frac{a}{\rho}}\right)

In actual manufactured hardware, physical cracks, fatigue tears, and weld lack-of-fusion defects have atomic-scale root radii approaching zero (ρ→0\rho \to 0). Under classical elasticity, as ρ→0\rho \to 0, Kt→∞K_t \to \infty and σmax⁡→∞\sigma_{\max} \to \infty. Infinite stress predictions make classical stress analysis incapable of assessing whether a cracked component can safely carry load. Linear Elastic Fracture Mechanics (LEFM) resolves this mathematical singularity by treating the crack tip not as an infinite point stress, but through energy dissipation and stress field intensity.

Griffith Energy Criterion (1921)

A.A. Griffith established that a crack propagates catastrophically when the rate of elastic strain energy released by infinitesimal crack advance (dU/dadU/da) equals or exceeds the rate of surface energy required to create two new fracture surfaces:

G=−dUda≥Gc=2γsG = -\frac{dU}{da} \ge G_c = 2\gamma_s

For an ideal brittle material, the critical fracture stress is:

σf=2Eγsπa\sigma_f = \sqrt{\frac{2E\gamma_s}{\pi a}}

In engineering metals, E. Orowan and G.R. Irwin showed that localized dislocation plasticity always occurs in a microscopic zone at the crack tip, absorbing thousands of times more energy than surface tension alone (γp≫γs\gamma_p \gg \gamma_s). Irwin modified the criterion to:

Gc=2(γs+γp)≈2γp  ⟹  σf=EGcπaG_c = 2(\gamma_s + \gamma_p) \approx 2\gamma_p \implies \sigma_f = \sqrt{\frac{E G_c}{\pi a}}

The Three Modes of Crack Surface Displacement

Crack tip stress fields are categorized into three fundamental kinematic modes:

  1. Mode I (Opening Mode): Tensile normal stresses act perpendicular to the crack plane, pulling crack faces directly apart (yy-direction). Mode I is the predominant failure mode in mechanical machine design, pressure vessels, and rotating shafts.
  2. Mode II (Sliding / In-Plane Shear Mode): Shear stresses act parallel to the crack plane and perpendicular to the crack front (xx-direction), causing crack faces to slide over one another.
  3. Mode III (Tearing / Anti-Plane Shear Mode): Out-of-plane shear stresses act parallel to the crack front (zz-direction), causing crack faces to twist and tear apart.

Mode I Stress Intensity Factor (KIK_I)

G.R. Irwin demonstrated that the asymptotic elastic stress field at distance rr and angle θ\theta from a Mode I crack tip is governed by a single scalar parameter:

σij(r,θ)=KI2πrfij(θ)+O(r0)\sigma_{ij}(r, \theta) = \frac{K_I}{\sqrt{2\pi r}} f_{ij}(\theta) + \mathcal{O}(r^0)

The Mode I Stress Intensity Factor (KIK_I) scales the entire local crack-tip stress field and is defined macroscopically as:

KI=YσπaK_I = Y \sigma \sqrt{\pi a}

where:

  • σ\sigma is the nominal gross far-field tensile stress acting normal to the crack plane.
  • aa is the characteristic crack dimension (a=depth of an edge cracka = \text{depth of an edge crack}; 2a=total length of an internal crack2a = \text{total length of an internal crack}, so a=half-crack lengtha = \text{half-crack length}).
  • YY (also denoted β\beta) is the dimensionless geometry correction factor, which accounts for crack geometry, component thickness, finite width (a/Wa/W), and boundary constraints.
    • Central through-crack of length 2a2a in an infinitely wide plate: Y=1.0Y = 1.0.
    • Single edge crack of depth aa in a semi-infinite plate: Y=1.12Y = 1.12 (the 12%12\% increase reflects the loss of mechanical constraint at the free edge).
  • Units: MPam\text{MPa}\sqrt{\text{m}} in SI; ksiin\text{ksi}\sqrt{\text{in}} in USCS. Conversion factor: 1 ksiin=1.09884 MPam1\text{ ksi}\sqrt{\text{in}} = 1.09884\text{ MPa}\sqrt{\text{m}}

Plane Strain Fracture Toughness (KIcK_{Ic}) & Fracture Criterion

When the applied stress intensity factor reaches a critical threshold, crack growth becomes unstable and propagates through the component at high speed (a sizable fraction of the material's sound speed, on the order of 1 km/s1\text{ km/s} in steel). This critical threshold under maximum triaxial constraint is the Plane Strain Fracture Toughness (KIcK_{Ic})—a fundamental material property measured experimentally in accordance with ASTM E399.

The fracture mechanics design criterion is:

KI≥KIc  ⟹  Unstable Brittle FractureK_I \ge K_{Ic} \implies \text{Unstable Brittle Fracture}

The factor of safety against brittle fracture is:

n=KIcKI=KIcYσπan = \frac{K_{Ic}}{K_I} = \frac{K_{Ic}}{Y \sigma \sqrt{\pi a}}

ASTM E399 Valid Plane Strain Thickness Requirement

Fracture toughness is strongly dependent on specimen thickness (BB):

  • Thin Plates (Plane Stress): At a free surface, out-of-plane normal stress is zero (σz=0\sigma_z = 0). The material is free to contract laterally (ϵz≠0\epsilon_z \ne 0), allowing extensive plastic yielding in a large yield zone ahead of the crack tip. This blunts the crack, dissipates substantial plastic work, and produces high apparent toughness (KcK_c).
  • Thick Plates (Plane Strain): Deep inside a thick plate, surrounding elastic material prevents lateral contraction (ϵz=0\epsilon_z = 0). This constraint generates a large triaxial tensile stress (σz≈ν(σx+σy)\sigma_z \approx \nu(\sigma_x + \sigma_y)). Triaxial tension severely suppresses plastic slip, producing a microscopic plastic zone and an invariant, lower-bound material toughness plateau: KIcK_{Ic}.

To ensure that experimental measurements or design calculations are valid lower-bound plane strain values, textbooks use the classic ASTM E399 size criterion:

B≥2.5(KIcSy)2anda≥2.5(KIcSy)2B \ge 2.5\left(\frac{K_{Ic}}{S_y}\right)^2 \quad \text{and} \quad a \ge 2.5\left(\frac{K_{Ic}}{S_y}\right)^2

where SyS_y is the 0.2%0.2\% offset yield strength of the material. Since the 2012 revision, ASTM E399 makes the uncracked ligament the mandatory size check, W−a≥2.5(KIc/σYS)2W - a \ge 2.5(K_{Ic}/\sigma_{YS})^2, together with standard specimen proportions (such as B≈W/2B \approx W/2) and a Pmax⁡/PQ≤1.10P_{\max}/P_Q \le 1.10 limit. These rules qualify a laboratory test specimen; in design, KIcK_{Ic} is used as a conservative lower-bound property for thick, highly constrained parts.

Irwin Plastic Zone Size

Because actual engineering metals yield at stresses exceeding SyS_y, a small zone of plastic deformation develops at the crack tip. Irwin estimated the radius of this plastic zone (ryr_y) as:

ry=12π(KISy)2(Plane Stress)r_y = \frac{1}{2\pi}\left(\frac{K_I}{S_y}\right)^2 \quad (\text{Plane Stress})

ry=16π(KISy)2(Plane Strain)r_y = \frac{1}{6\pi}\left(\frac{K_I}{S_y}\right)^2 \quad (\text{Plane Strain})

Under plane strain, triaxial constraint reduces the plastic zone volume by a factor of 3, promoting sudden cleavage.

Non-Destructive Evaluation (NDE) Crack Detection Thresholds

In modern damage-tolerant aerospace and pressure vessel design, structures are assumed to contain initial flaws (a0a_0) equal to the largest flaw that could reliably escape detection during quality inspection (typically established at a 90%90\% Probability of Detection with a 95%95\% confidence limit, 90/95 POD):

Inspection MethodologyAcronymTypical Surface Flaw Limit (adetecta_{\text{detect}})Subsurface Flaw Limit (adetecta_{\text{detect}})
Visual TestingVT5.0–10.0 mm5.0\text{--}10.0\text{ mm} (0.20–0.40 in0.20\text{--}0.40\text{ in})Not Detectable
Radiographic TestingRT1.5–3.0 mm1.5\text{--}3.0\text{ mm} (0.06–0.12 in0.06\text{--}0.12\text{ in})1.5–3.0 mm1.5\text{--}3.0\text{ mm} (0.06–0.12 in0.06\text{--}0.12\text{ in})
Ultrasonic TestingUT1.0–2.0 mm1.0\text{--}2.0\text{ mm} (0.04–0.08 in0.04\text{--}0.08\text{ in})1.0–2.0 mm1.0\text{--}2.0\text{ mm} (0.04–0.08 in0.04\text{--}0.08\text{ in})
Magnetic Particle TestingMT0.5–1.5 mm0.5\text{--}1.5\text{ mm} (0.02–0.06 in0.02\text{--}0.06\text{ in})Shallow only (<2 mm< 2\text{ mm})
Liquid Penetrant TestingPT0.5–1.5 mm0.5\text{--}1.5\text{ mm} (0.02–0.06 in0.02\text{--}0.06\text{ in})Not Detectable
Eddy Current TestingET0.2–0.8 mm0.2\text{--}0.8\text{ mm} (0.008–0.03 in0.008\text{--}0.03\text{ in})Near-surface only

Leak-Before-Break (LBB) Design Criterion

In the design of high-pressure piping, gas storage cylinders, and reactor pressure vessels (ASME Section XI), Leak-Before-Break (LBB) is the primary safety criterion against catastrophic explosion.

At a design operating stress σ\sigma, the critical crack size for unstable fast fracture is:

KI=KIc  ⟹  Yσπac=KIc  ⟹  ac=1π(KIcYσ)2K_I = K_{Ic} \implies Y \sigma \sqrt{\pi a_c} = K_{Ic} \implies a_c = \frac{1}{\pi}\left(\frac{K_{Ic}}{Y \sigma}\right)^2

Let tt represent the wall thickness of the pressure vessel:

  1. Leak-Before-Break Satisfied (ac>ta_c > t): If the critical crack length exceeds the wall thickness, a growing crack will breach the exterior surface (a=ta = t) while remaining stable. The pressurized fluid escapes through the through-wall slit, producing a detectable drop in system pressure or visible fluid venting. Operators can execute an orderly emergency shutdown before catastrophic fracture occurs.
  2. Leak-Before-Break Violated (ac<ta_c < t): If the critical crack depth is less than the wall thickness, the crack reaches unstable fast-fracture velocity while still partially embedded inside the wall. The vessel ruptures catastrophically into multiple flying fragments without prior leakage warning.

5. Worked Numerical Engineering Scenario: Linear Elastic Fracture Mechanics of a Heavy Pressure Wall

Problem Statement

A thick-walled tension plate forming part of a high-pressure hydraulic accumulator is fabricated from high-strength quenched and tempered alloy steel (AISI 4340). Quality control records document the following mechanical properties:

  • Yield strength: Sy=1200 MPaS_y = 1200\text{ MPa} (174 ksi174\text{ ksi})
  • Ultimate tensile strength: Sut=1400 MPaS_{ut} = 1400\text{ MPa} (203 ksi203\text{ ksi})
  • Plane strain fracture toughness: KIc=65.0 MPamK_{Ic} = 65.0\text{ MPa}\sqrt{\text{m}} (59.2 ksiin59.2\text{ ksi}\sqrt{\text{in}})

The plate has a rectangular cross-section with width W=200 mmW = 200\text{ mm} (0.200 m0.200\text{ m}) and thickness B=t=25.0 mmB = t = 25.0\text{ mm} (0.025 m0.025\text{ m}). It is subjected to a uniform nominal gross tensile stress of σ=350 MPa\sigma = 350\text{ MPa} (50.8 ksi50.8\text{ ksi}).

During a scheduled ultrasonic maintenance inspection (UT), an edge crack of depth a=6.0 mma = 6.0\text{ mm} (0.0060 m0.0060\text{ m}) is detected normal to the tensile stress axis. For this single-edge notched geometry, the geometry correction factor is governed by:

Y=1.12−0.231(aW)+10.55(aW)2−21.72(aW)3+30.39(aW)4Y = 1.12 - 0.231\left(\frac{a}{W}\right) + 10.55\left(\frac{a}{W}\right)^2 - 21.72\left(\frac{a}{W}\right)^3 + 30.39\left(\frac{a}{W}\right)^4

Perform the following engineering evaluations:

  1. Verify whether the component satisfies the ASTM E399 minimum thickness requirement for valid plane strain conditions.
  2. Calculate the operating Mode I stress intensity factor (KIK_I) and the resulting factor of safety against catastrophic brittle fracture.
  3. Determine the critical crack depth (aca_c) that would cause sudden brittle fracture under the current operating stress.
  4. Determine whether this component satisfies the Leak-Before-Break (LBB) safety criterion.

Step-by-Step Solution

Step 1: ASTM E399 Valid Plane Strain Thickness Verification

The minimum thickness required for valid plane strain fracture toughness behavior is:

Bmin⁡=2.5(KIcSy)2B_{\min} = 2.5\left(\frac{K_{Ic}}{S_y}\right)^2

Substitute the material properties:

Bmin⁡=2.5(65.0 MPam1200 MPa)2=2.5(0.054167m)2=2.5(0.002934 m)=0.007335 m=7.34 mmB_{\min} = 2.5\left(\frac{65.0\text{ MPa}\sqrt{\text{m}}}{1200\text{ MPa}}\right)^2 = 2.5(0.054167\sqrt{\text{m}})^2 = 2.5(0.002934\text{ m}) = 0.007335\text{ m} = 7.34\text{ mm}

Comparing with the actual plate dimensions:

  • Actual thickness: B=25.0 mmB = 25.0\text{ mm}, more than three times 7.34 mm7.34\text{ mm}, so the plate is thick enough for high crack-tip constraint.
  • Actual crack depth: a=6.0 mma = 6.0\text{ mm}, slightly less than 7.34 mm7.34\text{ mm}, so a specimen this shallow would not meet the classic size rule. A shallow crack loses some constraint, which raises the apparent toughness, so using the lower-bound KIcK_{Ic} is conservative.

The size criteria qualify a laboratory test; for this component the conclusion is that KIc=65.0 MPamK_{Ic} = 65.0\text{ MPa}\sqrt{\text{m}} is an appropriate (conservative) design toughness.

Step 2: Operating Stress Intensity Factor (KIK_I) & Factor of Safety

First, determine the crack-to-width ratio:

aW=6.0 mm200.0 mm=0.030\frac{a}{W} = \frac{6.0\text{ mm}}{200.0\text{ mm}} = 0.030

Calculate the geometry correction factor YY:

Y=1.12−0.231(0.030)+10.55(0.030)2−21.72(0.030)3+30.39(0.030)4Y = 1.12 - 0.231(0.030) + 10.55(0.030)^2 - 21.72(0.030)^3 + 30.39(0.030)^4

Y=1.12−0.00693+0.009495−0.000586+0.000025=1.1220≈1.122Y = 1.12 - 0.00693 + 0.009495 - 0.000586 + 0.000025 = 1.1220 \approx 1.122

Now calculate the applied Mode I stress intensity factor:

KI=Yσπa=(1.122)(350 MPa)π(0.0060 m)K_I = Y \sigma \sqrt{\pi a} = (1.122)(350\text{ MPa})\sqrt{\pi (0.0060\text{ m})}

π(0.0060)=0.0188496=0.137294m\sqrt{\pi (0.0060)} = \sqrt{0.0188496} = 0.137294\sqrt{\text{m}}

KI=(1.122)(350)(0.137294)=53.915 MPam≈53.92 MPamK_I = (1.122)(350)(0.137294) = 53.915\text{ MPa}\sqrt{\text{m}} \approx 53.92\text{ MPa}\sqrt{\text{m}}

The factor of safety against static brittle fracture is:

n=KIcKI=65.0 MPam53.92 MPam=1.205≈1.21n = \frac{K_{Ic}}{K_I} = \frac{65.0\text{ MPa}\sqrt{\text{m}}}{53.92\text{ MPa}\sqrt{\text{m}}} = 1.205 \approx 1.21

Note on Yielding: The nominal stress is σ=350 MPa\sigma = 350\text{ MPa}, giving a static tensile yield factor of safety ny=Sy/σ=1200/350=3.43n_y = S_y / \sigma = 1200 / 350 = 3.43. However, the presence of the 6 mm6\text{ mm} crack reduces the true safety margin to a precarious 1.211.21!

Step 3: Critical Crack Depth (aca_c) Under Operating Stress

At the point of catastrophic fracture, KI=KIc=65.0 MPamK_I = K_{Ic} = 65.0\text{ MPa}\sqrt{\text{m}}. Using Y≈1.12Y \approx 1.12 for shallow edge cracks:

KIc=Yσπac  ⟹  ac=1π(KIcYσ)2K_{Ic} = Y \sigma \sqrt{\pi a_c} \implies a_c = \frac{1}{\pi}\left(\frac{K_{Ic}}{Y \sigma}\right)^2

ac=1π(65.0 MPam(1.12)(350 MPa))2=1π(65.0392.0)2=1π(0.165816)2=0.027495π=0.008752 m=8.75 mma_c = \frac{1}{\pi}\left(\frac{65.0\text{ MPa}\sqrt{\text{m}}}{(1.12)(350\text{ MPa})}\right)^2 = \frac{1}{\pi}\left(\frac{65.0}{392.0}\right)^2 = \frac{1}{\pi}(0.165816)^2 = \frac{0.027495}{\pi} = 0.008752\text{ m} = 8.75\text{ mm}

Refining with a/W=8.75/200=0.0438a/W = 8.75/200 = 0.0438 and all five terms of the geometry polynomial gives Y=1.128Y = 1.128:

ac=1π(65.0(1.128)(350))2=1π(65.0394.8)2=1π(0.16464)2=0.00863 m=8.63 mma_c = \frac{1}{\pi}\left(\frac{65.0}{(1.128)(350)}\right)^2 = \frac{1}{\pi}\left(\frac{65.0}{394.8}\right)^2 = \frac{1}{\pi}(0.16464)^2 = 0.00863\text{ m} = 8.63\text{ mm}

(A second iteration changes nothing at this precision.) Catastrophic fracture will occur if the crack grows by an additional 8.63−6.00=2.63 mm8.63 - 6.00 = 2.63\text{ mm}.

Step 4: Leak-Before-Break (LBB) Evaluation

Treat the 25.0 mm25.0\text{ mm} plate thickness as the pressure wall and apply the simplified textbook LBB check, comparing the critical crack size with the wall thickness tt (a full LBB assessment uses surface-crack and through-wall-crack solutions for the real vessel geometry):

ac=8.63 mm<t=25.0 mma_c = 8.63\text{ mm} < t = 25.0\text{ mm}

  • Engineering Finding: The critical crack size (8.63 mm8.63\text{ mm}) is far smaller than the 25 mm25\text{ mm} wall thickness (ac/t=0.345a_c / t = 0.345, roughly one-third of the wall).
  • Safety Conclusion: The component fails the Leak-Before-Break criterion. A growing crack would reach the fast-fracture size long before it could breach the wall and leak, so the wall would rupture without a warning leak.
  • Corrective Engineering Action: For the crack to reach the full 25 mm25\text{ mm} wall before going unstable (Y≈1.22Y \approx 1.22 at a/W=0.125a/W = 0.125), the material would need KIc≥(1.22)(350)π(0.025)≈120 MPamK_{Ic} \ge (1.22)(350)\sqrt{\pi(0.025)} \approx 120\text{ MPa}\sqrt{\text{m}} at the present stress (for example, a tougher grade or a lower-strength temper of the steel). Alternatively, keeping KIc=65 MPamK_{Ic} = 65\text{ MPa}\sqrt{\text{m}}, the operating stress would have to drop to about 65/(1.22π(0.025))≈190 MPa65/(1.22\sqrt{\pi(0.025)}) \approx 190\text{ MPa}.
Test Your Knowledge

A machine base fabricated from brittle ASTM Class 30 gray cast iron (Sut=214 MPaS_{ut} = 214\text{ MPa}, Suc=750 MPaS_{uc} = 750\text{ MPa}) is subjected to a state of plane stress with in-plane principal stresses σA=+150 MPa\sigma_A = +150\text{ MPa} and σB=−75 MPa\sigma_B = -75\text{ MPa}. According to the Modified Mohr (MM) theory, what is the factor of safety against static failure?

A

1.18

B

1.43

C

2.85

D

1.00

Test Your Knowledge

A forged alloy steel structural tension tie has a yield strength of Sy=900 MPaS_y = 900\text{ MPa} and a plane strain fracture toughness of KIc=45.0 MPamK_{Ic} = 45.0\text{ MPa}\sqrt{\text{m}}. Using the classic ASTM E399 plane-strain size criterion B≥2.5(KIc/Sy)2B \ge 2.5(K_{Ic}/S_y)^2, what minimum thickness is needed for plane-strain (lower-bound) fracture behavior?

A

1.25 mm

B

3.15 mm

C

12.5 mm

D

6.25 mm

Test Your Knowledge

In the structural design of pressurized pipelines and pressure vessels, what condition defines compliance with the Leak-Before-Break (LBB) design criterion?

A

The vessel wall thickness must exceed the critical crack length (t>act > a_c) so that surface flaws cannot penetrate to the exterior.

B

The operating pressure must be low enough that the nominal hoop stress remains strictly below one-third of the yield strength.

C

The critical crack length for fast fracture must exceed the wall thickness (ac>ta_c > t), so a crack leaks before it ruptures.

D

The fracture toughness must be high enough that the crack tip plastic zone radius exceeds the wall thickness.

Sections you finish are checked off in the contents.