j-integral

name	j-integral
description	J-integral fracture mechanics — elastic-plastic fracture, ASTM E1820 J-R curve, η-factor method, validity requirements, J_IC to K_IC conversion, crack tip opening displacement (CTOD), ductile tearing, constraint effects (T-stress, Q-factor), piping flaw assessment, API 579/BS 7910.
metadata	{"priority":7,"promptSignals":{"phrases":["J-integral","J-R curve","elastic plastic fracture","J_IC","CTOD fracture","ductile fracture mechanics"],"minScore":3}}

J-Integral Fracture Mechanics — Complete Skill

J-Integral Definition

Rice's Path-Independent Integral

J-integral (Rice, 1968): J = ∮_Γ (W dy - T_i × ∂u_i/∂x ds) [N/m or J/m²; Γ = any contour surrounding crack tip] W = strain energy density = ∫₀^ε σ_ij dε_ij [N/m²] T_i = traction vector; u_i = displacement vector; ds = arc length element

Physical interpretation: J = energy release rate for crack extension per unit crack front length J = G (energy release rate) for linear elastic case For linear elastic: G = J = K²/E' where E' = E/(1-ν²) for plane strain; E' = E for plane stress

Equivalent energy definition (useful for experimental measurement): J = -(1/B) × dU/da [B = specimen thickness; U = total strain energy; a = crack length] Measured from area under load-displacement curve as crack extends

Elastic-Plastic Fracture (EPFM)

J-Dominance and HRR Field

HRR singularity (Hutchinson-Rice-Rosengren): σ_ij = σ_0 × (E × J / (α × σ_0² × I_n × r))^(1/(n+1)) × σ̃_ij(θ, n) [MPa; r = distance from tip; α, n = Ramberg-Osgood parameters; I_n = integration constant]

Ramberg-Osgood hardening law: ε/ε_0 = σ/σ_0 + α × (σ/σ_0)^n [ε_0 = σ_0/E; σ_0 = yield stress; n = hardening exponent; typically n = 5–20 for steels]

J-dominance zone radius: r_J = (1/2π) × J / (σ_0 × ε_0) [region where HRR singular field dominates] Valid EPFM: plastic zone size r_p >> r_J (J-dominance maintained)

Plastic Zone Size

Irwin plastic zone (plane stress): r_y = (1/2π) × (K/σ_y)² [m; r_y = plane stress plastic zone radius]

Plane strain plastic zone: r_y,3D = (1/6π) × (K/σ_y)² [3× smaller than plane stress; triaxial constraint suppresses plasticity]

Large-scale plasticity: when r_y > 0.1 × (W-a) or 0.1 × a → LEFM invalid → use EPFM with J

Experimental J-R Curve (ASTM E1820)

Specimen Types

Compact tension (CT): W/B = 2; a₀/W = 0.45–0.70; width W = 25–50 mm typical Single edge notch bend (SENB or 3PB): W/B = 2; span S = 4W Disk-shaped compact (DCT): for small material volumes

Fatigue pre-cracking requirement: ΔK at final 2.5% of precrack: ΔK ≤ 0.6 × K_max_valid [prevent large reversed plastic zone] Pre-crack length: 0.45 ≤ a₀/W ≤ 0.70 (per ASTM E1820)

J Calculation Methods

η-factor method (single specimen — compliance method): J = J_el + J_pl [elastic + plastic component]

Elastic component: J_el = K²/E' = (F × f(a/W))² / (E × B × W) × (1 - ν²) [K from stress intensity factor formula]

Plastic component (CT specimen): J_pl = η × U_pl / (B_N × b) [η = dimensionless factor; U_pl = plastic area under P-δ curve; B_N = net thickness; b = W - a = ligament] For CT: η = 2 + 0.522 × b/W For SENB: η = 2

**Incremental J (normalized): During crack growth, J updated for each Δa: J_i = (J_(i-1) + (η/b_(i-1)) × ΔU_pl / B_N) × (1 - γ × Δa_i / b_(i-1)) [crack growth correction; γ = 1 + 0.76 b/W for CT]

Compliance rotation: C = δ / F [elastic compliance used to infer crack length a at each unload] a = f(C, E, B, W) from ASTM E1820 compliance equation

J-R Curve Determination

Blunting line: J = 2 × σ_flow × Δa [σ_flow = (σ_y + σ_UTS)/2; linear portion before crack growth initiation; "blunting" of crack tip]

J_Q (provisional J at initiation): Intersection of J-R curve with blunting line at Δa = 0.2 mm exclusion line offset Exclusion lines: 0.15 mm offset (lower) and 1.5 mm offset (upper) Data within these lines used for power-law fit: J = C₁ × Δa^C₂

Validity checks for J_IC (per ASTM E1820): B ≥ 10 × J_Q / σ_flow [thickness sufficiency] b₀ ≥ 10 × J_Q / σ_flow [ligament sufficiency] a₀/W: 0.45–0.70 Δa_max ≤ 0.25 × b₀ [maximum crack extension for valid data] If valid: J_IC = J_Q

J_IC to K_IC equivalent: K_IC(J) = √(J_IC × E / (1 - ν²)) [MPa√m; valid for plane strain conditions] Compare to LEFM K_IC from ASTM E399 — should agree within 10–15%

Crack Tip Opening Displacement (CTOD)

CTOD (δ) definition: δ = 2 × u_y at r = r* from crack tip (usually measured at crack mouth or defined as displacement at original tip)

CTOD-J relationship (general): δ = J / (m × σ_y) [m = dimensionless factor; m = 1 for plane stress; m = 2 for plane strain hardening materials] For Prandtl fan: m = 1 + π/2 for non-hardening; m ≈ 2 for typical engineering steels

British Standard (BS 7910) CTOD: K² = E × σ_y × δ [simplified; relates K to CTOD for linear elastic case]

CTOD critical (δ_c): measured at fracture instability; equivalent to K_IC for LEFM CTOD at ductile initiation (δ_i): lower bound; used in BS 7910 assessments

BS 7448 — CTOD Testing

Similar to ASTM E1820: CT or SENB specimens; fatigue pre-crack; measure V-clip gage CTOD formula from BS 7448: δ = K²(1-ν²)/(2σ_ys × E) + 0.4(W-a₀)V_p / (0.4W + 0.6a₀ + z) [V_p = plastic component of clip gage; z = knife-edge height above surface]

Constraint Effects

T-Stress and Two-Parameter Fracture Mechanics

T-stress (second term in Williams expansion): σ_x = K/(√(2πr)) cos(θ/2)[1 - sin(θ/2)sin(3θ/2)] + T [T parallel to crack plane] T > 0: high constraint (plane strain-like); T < 0: low constraint (plane stress-like)

Biaxiality ratio: B = T × √(πa) / K [dimensionless; B = -1 for SENB; B = 0 for remote tension; B varies with a/W]

Q-factor (O'Dowd-Shih): σ_θθ(r,0) = σ_HRR(r,0) + Q × σ_0 [Q = constraint parameter; Q < 0 = lower constraint → higher apparent toughness] High constraint (CT specimen): Q ≈ 0 (conservative) Low constraint (large components with shallow crack): Q < -0.5 → J at fracture can be higher

Practical implication: Laboratory specimens (high constraint) give lower bound toughness Large structural members (often lower constraint) may have higher effective toughness SENB with a/W = 0.1: Q ≈ -0.5; SENB with a/W = 0.5: Q ≈ 0

Engineering Applications

Piping Flaw Assessment (API 579-1/ASME FFS-1 and BS 7910)

Failure Assessment Diagram (FAD): Two axes: K_r = K/K_IC (fracture axis); L_r = σ/σ_y (plasticity axis) Failure locus (Level 2): K_r = [f(L_r)]^(-1) = (1 - 0.14L_r²) × (0.3 + 0.7e^(-0.65L_r^6)) [for L_r ≤ L_r^max] Assessment point (K_r^a, L_r^a) inside locus → safe; outside → unsafe

K_r^a: calculated K_I from flaw in structure / K_IC L_r^a: reference stress σ_ref / σ_ys (normalized applied load) Margin factor: distance from assessment point to failure locus; target factor of safety ≥ 1.5

For pressure vessels and pipelines: RSTRENG (API 579 Level 1): remaining strength factor for corrosion; J-based Level 3 for complex flaws BS 7910 Level 3: full J-R curve tearing analysis; most rigorous; used for critical defect disposition

Ductile Tearing Resistance

Tearing modulus (T_mat): T_mat = (E/σ_0²) × (dJ/da) [dimensionless; measures slope of J-R curve; high T_mat = stable tearing] Typical values: 80–250 for structural steels; < 50 indicates risk of unstable tearing

Instability condition: dJ_applied/da > dJ_material/da → unstable tearing (crack runs) For fixed displacement loading: more stable (decreasing applied J as crack extends) For dead weight loading: less stable (applied J can increase with a)

Standards

Standard	Scope
ASTM E1820	J-integral and CTOD testing (unified standard)
ASTM E399	Linear elastic plane strain K_IC
BS 7448	CTOD testing
BS 7910	Guide to methods for assessing flaws in metallic structures
API 579-1/ASME FFS-1	Fitness-for-service (contains FAD procedures)
ISO 12135	Unified method for determining quasi-static fracture toughness
ASTM E1290	CTOD testing (now incorporated into E1820)

Output

Provide: test standard (ASTM E1820 or BS 7448), specimen type (CT/SENB) and dimensions (W [mm], B [mm], a₀ [mm], a₀/W), η-factor value for chosen specimen, J_el [kJ/m²] and J_pl [kJ/m²] at initiation, J_Q [kJ/m²] and validity check result (B/b₀ ≥ 10J_Q/σ_flow), J_IC [kJ/m²] (if valid), equivalent K_IC(J) [MPa√m], CTOD δ_i [mm], J-R curve fit parameters (C₁, C₂) and tearing modulus T_mat, constraint assessment (Q-factor, biaxiality B), FAD assessment point (K_r^a, L_r^a) vs. Level 2 locus (pass/fail), safety margin factor, and applicable standard (ASTM E1820, API 579, BS 7910).