What is the stress intensity factor K_I?

K_I characterizes the magnitude of the stress field near a mode-I crack tip. It is calculated as K_I = σ Y √(πa), where σ is the nominal stress, Y is the geometry factor, and a is the crack half-length. Fracture occurs when K_I reaches the material's fracture toughness K_Ic. Units are MPa√m.

Why does the geometry factor Y matter?

Different specimen geometries — center-crack tension (CCT), single-edge notch tension (SENT), compact tension (CT), or three-point bend (SENB) — produce different local stress fields at the crack tip. Y is a dimensionless correction factor that depends on a/W and accounts for these geometric effects. Accurate Y values are essential for correct K_I calculation.

How is the plastic zone radius calculated?

Using Irwin's first-order correction, the plastic zone radius is r_p = (1/2π)(K_I/σ_y)² for plane stress or r_p = (1/6π)(K_I/σ_y)² for plane strain, where σ_y is the yield stress. LEFM is valid when r_p ≪ a and r_p ≪ (W−a).

What is the critical crack size a_c?

a_c is the crack length at which unstable fracture occurs at the applied stress level: K_I(a_c) = K_Ic. Solving for a gives a_c = (1/π)(K_Ic / (σY))². Comparing a_c to the NDT detection threshold provides a margin of safety for structural integrity assessments.

› 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) 100

Crack length a (mm) 10

Specimen width W (mm) 100

Material Properties

Fracture toughness K_Ic (MPa√m) 50

Yield stress σ_y (MPa) 350

Key Equations

$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'$

—

K_I (MPa√m)

—

Safety Factor K_Ic/K_I

—

Plastic zone r_p (mm)

—

Critical a_c (mm)

—

Geometry factor Y

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).

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?"

Related Engineering Fields

The concept of stress intensity factor K is fundamental not only to the brittle fracture you've learned here but also active at the core of various engineering fields. The first that comes to mind is "fatigue analysis". The propagation rate of a crack under cyclic loading strongly correlates with the range of the stress intensity factor ΔK (=K_max - K_min). Using this relationship (Paris' law), you can predict the number of cycles for a crack to grow from a subcritical length to the critical length a_c, and design the machine's life accordingly.

Another important field is "stress corrosion cracking". When stress and a corrosive environment coexist, fracture can propagate at a stress intensity factor much lower than the material's inherent K_Ic. This threshold is called the "Stress Corrosion Cracking Threshold Stress Intensity Factor, K_Iscc". For example, K_Iscc can be a concern for aluminum alloys or high-strength steels in specific environments. You can assess the risk of stress corrosion cracking by calculating K_I with this tool and comparing it to the K_Iscc for that material and environment.

Looking further, it's also applied to "interface fracture" and "composite material fracture". To describe the delamination of adhesive interfaces or coating films, or crack propagation in CFRP (Carbon Fiber Reinforced Plastic), you need to consider not only Mode I (tensile opening) but also Mode II (in-plane shear) and Mode III (out-of-plane shear) stress intensity factors. The foundation of this tool is the first step towards learning these more complex fracture modes.

For Further Learning

If you're interested in the background of this tool's formulas or want to take the next step, I recommend studying in the following order. First, for the mathematical foundation, try following the derivation process of the stress field equation at the crack tip. The formula used in the tool, $$K_I = \sigma Y \sqrt{\pi a}$$, is actually derived from a solution in elasticity theory using complex function theory, known as "Westergaard's stress function". Understanding this derivation will make it clear why it's proportional to √a. Textbook-wise, the first few chapters of "Fundamentals of Fracture Mechanics" (for example, the one edited by the Japan Society of Mechanical Engineers) are quite thorough.

Next, move on to "Elastic-Plastic Fracture Mechanics (EPFM)". Here, new parameters like the J-integral and CTOD (Crack Tip Opening Displacement) take center stage. These are essential concepts for evaluating fracture in "tough" materials (e.g., pipeline steels) where a large plastic zone forms at the crack tip. Cases where the plastic zone radius r_p becomes large in this tool are precisely where EPFM comes into play.

Finally, for the practical stage, the quickest path is to experience "fracture analysis using the Finite Element Method (FEA)". Modern CAE software (like Abaqus, ANSYS, FrontSTR) comes equipped with functions to model cracks and directly calculate the J-integral or stress intensity factors. Once you can use the tool for simple shapes and verify complex real structures with FEA, you'll be a competent practitioner of fracture mechanics.