LEFM Calculator Back
Fracture Mechanics Calculator

Linear Fracture Mechanics (LEFM) Calculator

Input crack length, specimen geometry, and material properties to compute $K_I = \sigma Y\sqrt{\pi a}$ in real time. Visualize safety factor, plastic zone, and critical crack size.

Specimen Geometry
Loading & Crack
Applied stress σ (MPa)
MPa
Crack length a (mm)
mm
Specimen width W (mm)
mm
Material Properties
Fracture toughness K_Ic (MPa√m)
MPa√m
Yield stress σ_y (MPa)
MPa
Results
K_I (MPa√m)
Safety Factor K_Ic/K_I
Plastic zone r_p (mm)
Critical a_c (mm)
Geometry factor Y
Crack
Stress Intensity K vs Crack Length
Theory & Key Formulas
$K_I = \sigma\,Y\sqrt{\pi a}$
Plastic zone: $r_p = \dfrac{1}{2\pi}\!\left(\dfrac{K_I}{\sigma_y}\right)^{\!2}$
Critical crack: $a_c = \dfrac{1}{\pi}\!\left(\dfrac{K_{Ic}}{\sigma Y}\right)^{\!2}$
Energy release: $G = K_I^2/E'$

What is Linear Elastic Fracture Mechanics (LEFM)?

🙋
What exactly is the "stress intensity factor" this simulator calculates? It sounds like it measures how dangerous a crack is.
🎓
Exactly right! It's the key metric in fracture mechanics. Basically, K_I (pronounced "K-one") quantifies the intensity of the stress field right at the tip of a crack. In practice, it combines the applied stress, the crack size, and the shape of the component into one number. Try moving the "Applied Stress σ" and "Crack Length a" sliders above—you'll see K_I change instantly, showing how a bigger load or a longer crack makes the situation more critical.
🙋
Wait, really? So fracture happens when this K_I number hits a certain limit? What's the "Fracture Toughness K_Ic" parameter for then?
🎓
That's the crucial material property! Think of K_Ic as the material's inherent resistance to crack propagation. For instance, a brittle material like glass has a very low K_Ic, while a tough steel has a high one. Fracture occurs when the driving force (K_I from your load and crack) equals or exceeds the material's resistance (K_Ic). In the simulator, when the calculated K_I bar turns red and exceeds the K_Ic you've set, it predicts failure.
🙋
That makes sense. But the simulator also shows a "Plastic Zone Size". What's that, and why does the "Yield Stress σ_y" matter if we're talking about cracks?
🎓
Great question! Even in "Linear Elastic" analysis, the stress at the crack tip is theoretically infinite. In reality, the material yields a little bit, creating a small plastic zone. The size of that zone, r_p, tells us if LEFM is valid—if it's too large, the simple linear model breaks down. A common case is a very ductile material (low σ_y); you'll see a huge plastic zone. Adjust the σ_y control and watch r_p change dramatically. It's a check on the theory's own limits.

Physical Model & Key Equations

The fundamental equation of LEFM for a through-thickness crack. The stress intensity factor K_I scales with the applied stress and the square root of the crack length. The geometry factor Y (often around 1 for a simple edge crack) accounts for the shape of the specimen and crack.

$$K_I = \sigma\,Y\sqrt{\pi a}$$

Where:
K_I = Stress Intensity Factor (MPa√m)
σ = Applied Remote Stress (MPa)
Y = Dimensionless Geometry Correction Factor
a = Crack Length (m)

These derived equations provide critical safety insights. The plastic zone size estimates the region where material yields, and the critical crack length is the maximum flaw size allowable before fast fracture under a given stress.

$$r_p = \dfrac{1}{2\pi}\!\left(\dfrac{K_I}{\sigma_y}\right)^{\!2}\quad\quad\quad a_c = \dfrac{1}{\pi}\!\left(\dfrac{K_{Ic}}{\sigma Y}\right)^{\!2}$$

Where:
r_p = Radius of the estimated plastic zone at the crack tip.
σ_y = Material Yield Stress (MPa).
a_c = Critical Crack Length for fast fracture.
K_Ic = Material Fracture Toughness (MPa√m).

Frequently Asked Questions

Y depends on the specimen shape and the position/orientation of the crack. In this tool, typical shapes (such as center crack, single-edge crack, compact specimen, etc.) can be selected. If an appropriate Y is not chosen, the calculation accuracy of K_I will vary significantly, so be sure to select the one that matches the actual specimen shape.
When rp exceeds the specimen thickness, the plane strain condition is no longer valid, and the applicability limit of linear elastic fracture mechanics is exceeded. In such cases, elastic-plastic fracture mechanics (such as the J-integral) should be used. The K_I-based evaluation in this tool should be treated as a reference value, and caution is required for actual design.
If ac < a, linear elastic fracture mechanics indicates a high possibility of unstable fracture. However, actual fracture is also influenced by loading rate and environmental effects (such as corrosion). Therefore, it is recommended to consider a safety factor and, if necessary, conduct non-destructive testing or more detailed evaluation.
The nominal stress σ should be input as the stress acting on the cross-section assuming no crack exists (load/cross-sectional area). Since the stress concentration at the crack tip is accounted for by the shape correction factor Y, σ itself does not need to consider the presence or absence of a crack.

Real-World Applications

Aircraft Structural Inspection: LEFM is used to set inspection intervals for critical components. Engineers calculate the critical crack length (a_c) for parts like wing spars or landing gear. Based on predicted crack growth rates, they determine how often these parts must be inspected with ultrasound or X-rays to find cracks before they reach a dangerous size.

Pipeline and Pressure Vessel Design: For gas pipelines or nuclear reactor pressure vessels, engineers must ensure that even if a manufacturing flaw exists, it will not cause catastrophic failure under operating pressure. They use LEFM to define acceptable flaw sizes during quality control and to set safe operating pressures relative to the material's K_Ic.

Failure Analysis & Forensics: After a structural failure like a bridge collapse or a broken crankshaft, investigators use LEFM principles. By examining the fracture surface and estimating the stress at failure, they can work backwards to determine the initial crack size that caused the failure, helping to identify the root cause.

Material Selection for Critical Components: When designing a component where cracks are inevitable (e.g., turbine blades subjected to fatigue), engineers compare candidate materials based on their fracture toughness K_Ic. A higher K_Ic means the component can tolerate larger cracks before failing, leading to safer and more durable designs.

Common Misunderstandings and Points to Note

There are several key points you should be especially mindful of when starting to use this tool. First and foremost, the calculation results are only estimates based on an idealized model. The formulas used in the tool are fundamentally based on a simple model of a crack in an infinite plate. Actual components often have corners, holes, or multiple cracks, in which case the correction factor Y becomes more complex. For example, a surface crack on a 100mm diameter shaft and a center crack in a wide plate will have completely different K_I values even with the same σ and a. The golden rule is to first get a feel using the tool, then verify with more detailed standards (like JSME or ASTM codes) or FEA (Finite Element Analysis) in actual design.

Next, the interpretation of the "plastic zone radius". If the calculated r_p becomes larger than, say, 1/10 of the crack length a, that's a yellow flag. The premise of Linear Elastic Fracture Mechanics (LEFM)—that the plastic zone at the crack tip is sufficiently small—is starting to break down. For example, calculating K_I=100 MPa√m for a mild steel with σ_y=300MPa gives an r_p of about 1.4mm. This is fine for a crack with a=5mm, but it's a significant size for a crack with a=2mm. In such cases, you need the concepts of Elastic-Plastic Fracture Mechanics (EPFM).

Finally, the point that "the fracture toughness K_Ic varies depending on conditions". The K_Ic values in material tables are for room temperature and static loading. Actual structures may be at low temperatures, subject to cyclic loading, or in corrosive environments. For instance, the same steel can have a significantly lower K_Ic at -40°C compared to room temperature (embrittlement). For dynamic loading, you'll separately need the da/dN-ΔK curve (fatigue crack growth characteristics) which represents crack propagation rate. Even after getting an answer from the tool, always ask yourself, "What are the actual service conditions?"

How to Use

  1. Enter applied stress (MPa) in valSigNum or drag sliderSig; typical range 50–500 MPa for steel structures.
  2. Input crack length (mm) in valANum or sliderA; common inspection limits are 0.5–5 mm surface cracks.
  3. Set plate width (mm) in valWNum or sliderW; standard sheet widths are 50–500 mm depending on component.
  4. Enter material fracture toughness K_Ic (MPa√m) in valKicNum or sliderKic; structural steel typically 50–100 MPa√m, aluminum 25–45 MPa√m.
  5. Read K_I stress intensity factor, safety factor (K_Ic/K_I ratio), plastic zone radius r_p, critical crack length a_c, and geometry factor Y in real time.

Worked Example

ASTM A36 steel plate under tension: applied stress 200 MPa, surface crack 2.0 mm, plate width 150 mm, K_Ic = 65 MPa√m. Using Newman-Raju correction for edge crack, Y ≈ 1.26. Calculated K_I = 200 × 1.26 × √(π × 0.002) ≈ 9.98 MPa√m. Safety factor = 65 / 9.98 ≈ 6.5 (acceptable). Plastic zone r_p = (1/6π)(K_I/σ_y)² ≈ 0.024 mm (elastic-plastic regime valid). Critical crack length a_c = (K_Ic / (Y × σ))² / π ≈ 5.2 mm before fracture.

Practical Notes

  1. LEFM validity requires K_I < 0.6 K_Ic and plastic zone < 10% of specimen thickness; validate material thickness ≥ 12.5(K_Ic/σ_y)².
  2. Geometry factor Y depends on crack location (center vs. edge) and aspect ratio a/W; edge cracks require Y-correction multiplier 1.12–1.4.
  3. Fatigue cracks in aluminum 7075-T73 (K_Ic ≈ 28 MPa√m) demand lower stress thresholds than ferritic steel; periodic inspections at 50% critical crack length.
  4. Environmental corrosion reduces K_Ic by 15–25% in seawater; apply knockdown factor 0.75–0.85 for offshore structures.