3. 应力强度因子K
在上一节中，我们通过裂纹尖端附近的应力场分析看到：无论远场载荷如何，裂尖附近的应力分布都呈现出 \(1/\sqrt{r}\) 的奇异性。这个奇异场的“强度”由一个参数唯一决定，它就是应力强度因子 K。K 是线弹性断裂力学（LEFM）的核心参量，它把远场载荷、裂纹几何和构件尺寸统一压缩成一个数，用来判断裂纹是否会失稳扩展。
一、K 的定义与物理意义
以 I 型（张开型）裂纹为例，裂尖附近应力场可写为：
σ_ij(r,θ) = K_I / √(2πr) · f_ij(θ) + 高阶项其中 r 是到裂尖的距离，θ 是以裂纹延长线为基准的极角，f_ij(θ) 是只与角度有关的无量纲函数。当 r→0 时，高阶项可忽略，应力场完全由 K_I 控制。因此 K 的物理意义是：它表征裂纹尖端弹性应力场的奇异强度，而不是某一点的应力值。K 越大，裂尖附近同样距离处的应力越高，裂纹越容易扩展。
K 的单位由定义式可推出：σ 的单位是 MPa，√r 的单位是 √mm，因此 K 的量纲为 MPa·√m 或 MPa·√mm。工程中常用 MPa·√m，注意换算：1 MPa·√m ≈ 31.62 MPa·√mm。K 不是应力，也不是能量，而是一个“场强度参量”。
二、K 的三种基本类型
按裂纹面相对位移方式，K 分为三类：
- K_I：张开型，裂纹面垂直于拉应力方向分离，最常见、最危险；
- K_II：滑开型，裂纹面在面内相对滑移，由剪应力驱动；
- K_III：撕开型，裂纹面在面外相对错开，由反平面剪应力驱动。
工程中绝大多数断裂问题以 I 型为主，因此后续讨论默认 K = K_I。当多种载荷同时存在时，可用等效应力强度因子或混合型准则处理，但本节聚焦 I 型。
三、常见裂纹构型下的 K 表达式
K 的一般形式可写成：
K = Y · σ · √(πa)其中 σ 为远场名义应力，a 为裂纹尺寸（边缘裂纹取裂纹长度，中心裂纹取半长），Y 为几何修正因子，也称形状因子，它集中体现了裂纹位置、构件边界、加载方式的影响。Y 通常通过手册、权函数法或有限元标定获得。下面给出几种典型构型：
| 构型 | K 表达式 | 说明 |
|---|---|---|
| 无限大板中心穿透裂纹 | K = σ√(πa) | Y = 1，a 为半裂纹长度 |
| 无限大板单边穿透裂纹 | K = 1.12σ√(πa) | Y ≈ 1.12，自由表面修正 |
| 无限大板双边穿透裂纹 | K = 1.12σ√(πa) | a 为单侧裂纹长度 |
| 圆孔边单侧裂纹 | K = F·σ√(πa) | F 随 a/R 变化，孔边应力集中放大 |
| 三点弯曲 SE(B) 试样 | K = f(a/W)·P·S/(B·W^{3/2}) | f 为标准化函数，见 ASTM E399 |
| 紧凑拉伸 CT 试样 | K = f(a/W)·P/(B·√W) | f 由标准表格给出 |
以中心裂纹为例，若 σ = 200 MPa，a = 5 mm = 0.005 m，则 K = 200×√(π×0.005) ≈ 200×0.1253 ≈ 25.1 MPa·√m。若材料断裂韧度 K_IC = 30 MPa·√m，则 K < K_IC，裂纹暂不扩展；若 a 增大到 8 mm，K ≈ 31.7 MPa·√m > K_IC，裂纹将失稳扩展。这说明 K 把“载荷—裂纹尺寸—材料抗力”放在同一尺度上比较。
四、K 控制断裂的适用条件
K 判据 K = K_IC 并非无条件成立，它建立在线弹性断裂力学框架下，必须满足以下条件：
- 材料宏观线弹性：裂纹尖端以外区域应力—应变关系近似线性，卸载可恢复。
- 小范围屈服：裂尖塑性区尺寸 r_p 远小于裂纹尺寸 a 和韧带宽度 W−a。对 I 型裂纹，平面应变下 r_p ≈ (1/6π)(K/σ_y)²。若 r_p 与 a 同量级，K 场被塑性区严重扰动，需改用弹塑性断裂参量 J 积分或 CTOD。
- 裂纹尺寸足够大：a 需大于材料晶粒尺寸和塑性区尺寸，否则裂纹尖端无法形成稳定的 K 主导区。
- 平面应变状态：厚板中心区域处于平面应变，K_IC 才是真正的材料常数；薄板平面应力下测得的 K_C 随厚度变化，不是材料常数。
- 载荷单调或准静态：K 判据适用于单调加载下的起裂与失稳，疲劳裂纹扩展需用 Paris 公式等描述。
五、工程判断示例
某压力容器筒体材料 K_IC = 55 MPa·√m，屈服强度 σ_y = 450 MPa，工作应力 σ = 180 MPa，检测发现一条长 2a = 12 mm 的中心穿透裂纹。先算 K：a = 6 mm = 0.006 m，K = 180×√(π×0.006) ≈ 180×0.1373 ≈ 24.7 MPa·√m。再算塑性区：r_p ≈ (1/6π)(24.7/450)² ≈ 0.00016 m = 0.16 mm，远小于 a = 6 mm，满足小范围屈服。K = 24.7 < K_IC = 55，裂纹不会失稳扩展。但若裂纹扩展到 2a = 40 mm，则 K ≈ 180×√(π×0.02) ≈ 45.1 MPa·√m，仍小于 55；当 2a ≈ 60 mm 时，K ≈ 55.2 MPa·√m，达到临界，裂纹将失稳。由此可反推临界裂纹尺寸 a_c = (1/π)(K_IC/(Yσ))² ≈ (1/π)(55/180)² ≈ 0.0297 m ≈ 29.7 mm，即半长 29.7 mm、全长约 59.4 mm。这个 a_c 就是安全评估和检测周期制定的关键依据。
掌握 K 的定义、量纲、常见构型表达式和适用条件，是后续学习 K 的有限元计算、J 积分以及裂纹扩展仿真的基础。在实际仿真中，K 通常通过裂尖奇异单元、围线积分或位移外推法提取，而本节给出的解析表达式正是验证仿真结果是否合理的第一把标尺。