引言
传统强度设计以”应力不超过许用值”为准则，但大量工程事故——压力容器爆破、桥梁坍塌、飞机结构失事——发生在名义应力远低于屈服强度的工况下，根源在于结构中预先存在的裂纹类缺陷。断裂力学正是研究含裂纹体力学行为的学科：Griffith 于 1920 年从能量平衡角度解释了玻璃的实际强度为何远低于理论强度，Irwin 在 1950 年代将其发展为以应力强度因子为核心的定量工程理论。本文推导裂纹尖端渐近场、应力强度因子与能量释放率的概念及其相互关系，讨论断裂判据的适用条件，并以含边缘裂纹的高强钢构件为例完成安全评估。
裂纹尖端渐近场与应力强度因子
按裂纹面位移模式，裂纹分为 I 型（张开型）、II 型（滑开型）与 III 型（撕开型），其中 I 型最危险也最常用。对弹性裂纹体做渐近分析，I 型裂纹尖端附近的应力场具有普遍形式：
$$\sigma_{ij} = \frac{K_{\mathrm{I}}}{\sqrt{2\pi r}},f_{ij}(\theta) + O(r^{0})$$
其中 $r$、$\theta$ 为以裂尖为原点的极坐标，$f_{ij}(\theta)$ 是已知的角度分布函数（例如裂尖延长线上 $\theta = 0$ 处 $\sigma_{yy} = K_{\mathrm{I}}/\sqrt{2\pi r}$）。应力随 $r^{-1/2}$ 趋于无穷，说明”裂尖应力”本身不能作为强度指标——真正刻画场强弱的是奇异项的幅值 $K_{\mathrm{I}}$，即应力强度因子，量纲为 $\text{MPa}\cdot\sqrt{\text{m}}$。其通式为：
$$K_{\mathrm{I}} = Y,\sigma\sqrt{\pi a}$$
$\sigma$ 为远场应力，$a$ 为裂纹尺寸，$Y$ 为几何修正系数：无限大板中心裂纹 $Y = 1$，半无限板边缘裂纹 $Y \approx 1.12$，有限尺寸与孔边、表面裂纹等情形的 $Y$ 可查手册或用权函数法、有限元法计算。
能量释放率与 Griffith–Irwin 关系
Griffith 的能量观点认为，裂纹扩展需要消耗表面能，驱动力来自弹性应变能的释放。定义能量释放率 $G$ 为裂纹扩展单位面积时系统势能的下降率（量纲 J/m²）。对无限大板中心裂纹（总长 $2a$），能量平衡给出临界断裂应力：
$$\sigma_f = \sqrt{\frac{2E\gamma_s}{\pi a}}$$
其中 $\gamma_s$ 为表面能密度。该式正确预言了”强度与裂纹尺寸平方根成反比”，但对金属材料严重低估断裂阻力，因为金属裂尖伴随塑性变形，实际耗散能 $\gamma_{\mathrm{eff}}$ 比表面能大几个量级。Irwin 的关键贡献是证明线弹性条件下 $G$ 与 $K$ 严格等价：
$$G = frac{K_{mathrm{I}}^2}{E’}, qquad E’ = begin{cases} E, & text{平面应力} \ E/(1-nu^2), & text{平面应变} end{cases}$$
于是能量准则 $G \geq G_c$ 与应力强度因子准则 $K_{\mathrm{I}} \geq K_{\mathrm{Ic}}$ 完全等价。$K_{\mathrm{Ic}}$ 称平面应变断裂韧度，是材料常数，须按标准试验在足够厚度（满足平面应变与小范围屈服）的试样上测定。工程常用材料的典型断裂韧度范围如下表。
| 材料 | $K_{\mathrm{Ic}}$ 典型范围（MPa·√m） | 屈服强度量级（MPa） |
|---|---|---|
| 高强度钢（如 4340） | 50–120 | 1200–1900 |
| 铝合金 7075-T6 | 24–29 | 500 左右 |
| 钛合金 Ti-6Al-4V | 55–115 | 800–1100 |
| 铸铁 | 6–20 | — |
| 有机玻璃 PMMA | 0.7–1.6 | 50–75 |
| 普通混凝土 | 0.2–1.4 | —（抗压 20–40） |
断裂判据与适用条件
线弹性断裂力学（LEFM）的断裂判据为 $K_{\mathrm{I}} \geq K_{\mathrm{Ic}}$，其成立有两个前提。其一为小范围屈服：裂尖塑性区必须远小于裂纹尺寸与韧带尺寸。Irwin 估计塑性区尺寸为：
$$r_p = frac{1}{alphapi}left(frac{K_{mathrm{I}}}{sigma_y}right)^2, qquad alpha = begin{cases} 1, & text{平面应力} \ 3, & text{平面应变} end{cases}$$
其二为平面应变约束：试样厚度 $B \geq 2.5(K_{\mathrm{Ic}}/\sigma_y)^2$，否则测得的韧度偏高、偏于不保守。当塑性区与裂纹尺寸相当时 LEFM 失效，须转入弹塑性断裂力学，采用 $J$ 积分或裂纹尖端张开位移（CTOD）作为特征参数。此外，实际结构中的裂纹多在循环载荷下由疲劳萌生扩展，常用 Paris 公式 $\mathrm{d}a/\mathrm{d}N = C(\Delta K)^m$ 描述亚临界扩展，把”剩余寿命”问题同样归结为应力强度因子的计算。
工程案例分析：含边缘裂纹构件的安全评估
某高强度结构钢构件承受远场工作应力 $\sigma = 350$ MPa，材料参数为 $E = 210$ GPa、$\nu = 0.30$、$K_{\mathrm{Ic}} = 77$ MPa·√m、屈服强度 $\sigma_y = 1200$ MPa。表面检查发现深度 $a_0 = 1$ mm 的边缘裂纹，按半无限板边缘裂纹取 $Y = 1.12$。临界裂纹尺寸为：
$$a_c = \frac{1}{\pi}\left(\frac{K_{\mathrm{Ic}}}{Y\sigma}\right)^2 = \frac{1}{\pi}\left(\frac{77}{1.12\times350}\right)^2 \approx 12.3 \text{mm}$$
以下程序计算当前裂纹的 $K_{\mathrm{I}}$、能量释放率与塑性区尺寸，并扫描裂纹长大过程中各量的演化：
import numpy as np
E, nu = 210e3, 0.30 # 弹性模量 MPa，泊松比
K_Ic = 77.0 # 平面应变断裂韧度 MPa·√m
sig_y = 1200.0 # 屈服强度 MPa
sigma = 350.0 # 远场工作应力 MPa
Y = 1.12 # 半无限板边缘裂纹几何修正系数
K = lambda a: Y * sigma * np.sqrt(np.pi * a) # a 单位 m
a_c = (K_Ic / (Y * sigma))**2 / np.pi
print(f"临界裂纹尺寸 a_c = {a_c*1e3:.2f} mmn")
Ep = E / (1 - nu**2) # 平面应变等效模量
print(" a (mm) K_I (MPa√m) G (kJ/m^2) r_p (mm)")
for a_mm in [1, 2, 4, 6, 8, 10, 12]:
a = a_mm * 1e-3
Ki = K(a)
G = Ki**2 / Ep * 1e3 # MPa·m → kJ/m²
rp = (Ki / sig_y)**2 / np.pi * 1e3 # 平面应力塑性区，m→mm
print(f"{a_mm:7.1f} {Ki:12.2f} {G:10.3f} {rp:9.4f}")
K0 = K(1e-3)
print(f"n初始裂纹 a0=1 mm：K_I = {K0:.2f} MPa√m，断裂安全系数 n = {K_Ic/K0:.2f}")计算结果表明：初始 $a_0 = 1$ mm 裂纹的 $K_{\mathrm{I}} \approx 22.0$ MPa·√m，断裂安全系数约 3.5，对应能量释放率 $G \approx 2.1$ kJ/m²（临界值 $G_c = K_{\mathrm{Ic}}^2/E' \approx 25.7$ kJ/m²），塑性区半径仅约 0.11 mm，远小于裂纹尺寸与构件尺寸，LEFM 适用；但当裂纹扩展到 12 mm 量级时 $K_{\mathrm{I}}$ 迅速逼近 $K_{\mathrm{Ic}}$，说明安全评估对裂纹尺寸高度敏感。工程上应结合无损检测能力与 Paris 公式估算从检出尺寸到临界尺寸的剩余寿命，据此制定检修周期，这正是”损伤容限设计”的基本逻辑。
结论与展望
断裂力学把含裂纹结构的强度问题归结为单一参数的计算与测定：应力强度因子 $K$ 描述驱动力，断裂韧度 $K_{\mathrm{Ic}}$ 描述材料阻力，能量释放率 $G$ 则架起与 Griffith 能量观点的桥梁。小范围屈服与平面应变约束是 LEFM 的适用边界，超出边界须借助 $J$ 积分与 CTOD 等弹塑性参数。当前的研究前沿包括非均质材料与焊接残余应力场中的裂纹分析、多轴与变幅载荷下的疲劳裂纹扩展、以及与相场断裂模拟等数值方法的结合，这些方向使断裂力学持续成为结构完整性评估的核心理论。
发表回复