断裂评定的第一步，是算出裂纹尖端的应力强度因子 $K_I$。本文讲清 $K_I$ 到底是什么、BS 7910 附录 M（Annex M）怎么把它算出来——从通用母公式、各修正因子，到最常用的半椭圆表面裂纹解与焊趾修正。它对应失效评定图（FAD）的纵轴 $K_r$，与 Annex P（参考应力） 的横轴 $L_r$ 一起完成 Clause 7 断裂评定 。
引子：裂纹尖端的应力，是"无穷大"
含缺陷结构做断裂评定，第一步要回答：这条裂纹驱动力有多大？
传统强度设计用「工作应力 < 许用应力」。可一旦有了裂纹，按线弹性理论，裂纹尖端的应力会趋于无穷大——再也没法拿一个有限的应力值去和许用应力比较。断裂力学换了一把新尺子，就是应力强度因子 $K_I$。
一、为什么用 $K_I$，而不用普通应力
裂纹尖端附近的应力场有一个通用形式（线弹性断裂力学）：
$$\sigma_{ij} = \frac{K_I}{\sqrt{2\pi r}}\, f_{ij}(\theta) + \cdots$$其中 $r$ 是到裂尖的距离、$\theta$ 是方位角。当 $r \to 0$ 时应力按 $1/\sqrt{r}$ 发散——这就是裂尖应力奇异性。
关键在于：尽管尖端应力无穷大，但整个奇异应力场的"强弱"只由一个系数决定，这个系数就是 $K_I$。所以：
- 普通应力在裂尖无限大、没法直接用；$K_I$ 是有限、可算、可测的量。
- 只要知道 $K_I$，裂尖附近的应力/应变/位移分布就完全确定（形状一样，只差幅度）。
- 不同几何、不同裂纹尺寸、不同载荷，都能统一用 $K_I$ 来度量"这条裂纹有多危险"。
$K_I$ 的量纲是 $\mathrm{MPa}\sqrt{\mathrm m}$。注意别把它和无量纲的应力集中系数 $k_t$ 混为一谈——后者描述无裂纹几何突变处的应力放大，是两码事。
二、$K_I$ 与断裂韧性：脆断判据
$K_I$ 是裂纹的驱动力；材料抵抗裂纹扩展的能力叫断裂韧性 $K_{mat}$（由试验测定，随温度/约束/厚度变化）。二者一比，就得到脆断判据：
$$K_I \ge K_{mat} \quad\Rightarrow\quad \text{裂纹失稳扩展（脆性断裂）}$$在失效评定图（FAD）里，纵轴正是断裂比 $K_r = K_I/K_{mat}$，衡量"离脆断有多近"。所以 Annex M 算出的 $K_I$，就是喂给 FAD 纵轴的原料。
三、Annex M 的通用计算框架
Annex M 把所有几何的 $K_I$ 解都套进同一个母公式（BS 7910:2019, §M.1, Eq. M.1）：
$$K_I = (Y\sigma)\sqrt{\pi a}$$三个要素：$\sqrt{\pi a}$ 是裂纹特征尺寸 $a$ 的开方项；$\sigma$ 是驱动应力；$Y$ 是几何修正因子的总称，把「理想无限大平板」修正到「真实有限构件」。Annex M 的全部工作，就是给不同几何写出具体的 $Y\sigma$。
对一次应力，$Y\sigma$ 展开为（BS 7910:2019, §M.1, Eq. M.4）：
$$(Y\sigma)_p = M f_w\left[k_{tm}M_{km}M_m\,P_m + k_{tb}M_{kb}M_b\big(P_b+(k_m-1)P_m\big)\right]$$看着吓人，其实就是膜应力项 + 弯曲应力项各自被一串修正因子放大。二次应力则简单得多（BS 7910:2019, §M.1, Eq. M.5）：$(Y\sigma)_s = M_m Q_m + M_b Q_b$。
各修正因子的物理含义：
| 符号 | 名称 | 物理含义 | 何时 = 1 |
|---|---|---|---|
| $M_m$ / $M_b$ | 膜 / 弯曲应力修正系数 | 裂纹形状与位置对 $K_I$ 的几何放大 | 由各几何解给出 |
| $M$ | 鼓胀效应系数 | 曲壳内压下裂纹处壳壁外鼓的附加放大 | 平板 / 无鼓胀 |
| $f_w$ | 有限宽修正系数 | 裂纹占截面比例大时净截面减小、应力集中 | 裂纹相对很小 |
| $M_{km}$ / $M_{kb}$ | 焊趾应力放大系数 | 裂纹位于焊趾局部应力集中区的额外放大 | 不在焊趾 |
| $k_{tm}$ / $k_{tb}$ | 结构应力集中系数 | 接管 / 节点等粗大结构突变的应力集中 | 无结构突变 |
| $k_m$ | 错边因子 | 错边 / 角变形使膜应力派生附加弯曲 $(k_m-1)P_m$ | 无错边 |
这里的 $P_m$、$P_b$ 是一次应力（外载荷产生：内压、拉力、弯矩），$Q_m$、$Q_b$ 是二次应力（自平衡：焊接残余应力、热应力）。二者在 FAD 里角色不同：
图1：一次应力 vs 二次应力。一次应力（Primary，Pm/Pb）由外载荷维持、不自平衡，既进 K_r 又进 L_r；二次应力（Secondary，Qm/Qb，如焊接残余应力）自平衡、随塑性流动松弛，只进 K_r（且要乘塑性交互修正），不进 L_r。
一次应力 $K_I^P$ 与二次应力 $K_I^S$ 叠加进纵轴：$K_r=(K_I^P+K_I^S)/K_{mat}+\rho$。当 $K_I^S<0$（二次应力为压、有利）时，把 $K_I^S$ 与 $\rho$ 一并置零（保守）。
四、最常用的解：半椭圆表面裂纹
工程上绝大多数表面缺陷（焊趾裂纹、疲劳裂纹、腐蚀坑）都被表征为半椭圆表面裂纹，用 Newman–Raju 解（BS 7910:2019, §M.4.1）。它也是 Annex M 的主力解。
图2：三种基本缺陷的尺寸定义。左上为表面裂纹（深度 a、全长 2c、壁厚 B）；右上为埋藏裂纹（高 2a、长 2c、到近表面距离 p）；下为穿透裂纹（全长 2a）。半椭圆表面裂纹用深度 a 与半长 c 描述，长深比 a/c 刻画裂纹的圆扁。
膜载荷下的膜应力放大因子（BS 7910:2019, §M.4.1, Eq. M.10）：
$$M_m = \left[M_1 + M_2\left(\tfrac{a}{B}\right)^2 + M_3\left(\tfrac{a}{B}\right)^4\right]\frac{g\,f_\theta}{\Phi}$$其中方括号内是裂纹深度 $a/B$ 的多项式（系数 $M_1,M_2,M_3$ 只依赖长深比 $a/c$，且按 $a/2c\le 0.5$ 与 $>0.5$ 分两套）；$\Phi$ 是缺陷形状系数，$\Phi=\big[1+1.464(a/c)^{1.65}\big]^{0.5}$；$g$、$f_\theta$ 描述沿椭圆前缘的分布。弯曲载荷不另起炉灶，直接令 $M_b=H\cdot M_m$（BS 7910:2019, §M.4.1, Eq. M.12），$H$ 把弯曲应力沿厚度的线性衰减折进去。
求值点：最深点 A 与表面点 C
$K_I$ 沿半椭圆前缘是变化的，工程上评定两个特殊点：
- 最深点 A（$\theta=\pi/2$）：裂纹最深处，决定「会不会穿透壁厚」。
- 表面点 C（$\theta=0$）：裂纹与自由表面相交端，决定「会不会沿表面变长」。
表面裂纹两点都要算（BS 7910:2019, §7.2.12），取较危险者进 FAD——究竟哪点先临界，取决于长深比与载荷类型（拉伸 / 弯曲）。
有限宽修正对表面裂纹取 $f_w=\big\{\sec\big[(\pi c/W)(a/B)^{0.5}\big]\big\}^{0.5}$：裂纹相对板宽越小 $f_w\to 1$；越长越深则 $f_w$ 越大。适用上限约 $2c/W\le 0.8$。
五、焊接接头：焊趾应力放大系数 $M_k$
裂纹若位于焊趾（焊缝与母材相交处），那里有强烈的局部应力集中，$K_I$ 要额外乘一个焊趾应力放大系数 $M_k$（BS 7910:2019, §M.11）。$M_k$ 从表面的高值随深度快速衰减，到约 30% 壁厚处回到 1：
图3：焊趾放大因子 Mk 随裂纹深度的衰减。裂纹很浅（靠近焊趾表面应力集中区）时 Mk 明显大于 1；随裂纹加深穿出焊趾局部场，Mk 迅速回落到 1（母材场）。膜、弯曲各有一条 Mkm/Mkb 曲线。
只有裂纹在焊趾时才用 $M_k$；不在焊趾则取 1。同理：$f_w$ 仅当裂纹占截面比例较大时启用，$k_m$ 仅当有错边时，鼓胀因子 $M$ 仅对曲壳内压轴向裂纹——随意乱加修正会把结果搞坏。
六、不同几何的 $K_I$ 解（速览）
Annex M 按几何族给解，选解就是「对号入座」（BS 7910:2019, Annex M）：
| 几何 | 覆盖裂纹类型 | 关键特点 |
|---|---|---|
| 平板 | 穿透 / 边缘 / 表面半椭圆 / 埋藏 / 角 | $M_m/M_b$ 多项式 + $f_w$ |
| 圆筒·轴向 | 穿透 / 内外表面 / 埋藏 | 内压下鼓胀因子 $M$ |
| 圆筒·环向 | 穿透 / 内表面 | 全局弯矩折算等效膜应力 |
| 球壳 | 赤道穿透 | 鼓胀 + 影响系数 |
| 圆棒 / 螺栓 | 表面 / 环向 | 圆截面 + 螺纹几何 |
| 焊接接头 | 焊趾 / 焊根 | 焊趾放大因子 $M_k$ |
选解的四个判断：曲壳还是平板（有无鼓胀）？表面还是穿透（要不要评两点）？有限长还是扩展（三维半椭圆 vs 二维）？在不在焊趾（要不要 $M_k$）？
⚠️ 两条纪律：① 不能外推——每个解都有 $a/B$、$a/c$、$2c/W$ 适用范围，超范围会严重失真，须改用权函数多项式解或有限元；② 鼓胀不能算两遍——「平板解另乘 $M$」与「已含鼓胀的圆筒专用解」二选一，绝不混用。
七、计算步骤
从「几何 + 载荷 + 裂纹尺寸」到「算出 $K_I$」（BS 7910:2019, §M 与 §7.2）：
- 表征缺陷：把实测缺陷规则化为标准形（表面 / 埋藏 / 穿透 / 角），量出 $a$、$c$、$B$、$W$。（必做）
- 分解应力：把应力拆成膜 $P_m$ + 弯曲 $P_b$（非线性廓线则拟合为多项式）；区分一次 $P$ 与二次 $Q$。（必做）
- 选几何解：按缺陷几何在 Annex M 中对号入座，确认适用范围。（必做）
- 算修正因子：按需算 $M_m/M_b$、$f_w$（占截面 >10% 才用）、$M_k$（焊趾）、$M$（曲壳内压）、$k_m$（错边）。（按需）
- 组装 $K_I$：按 Eq. M.4 逐项乘累得 $K_I^P$，按 Eq. M.5 得 $K_I^S$。（必做）
- 在求值点取值：表面裂纹分别算最深点 A、表面点 C，取危险者。（必做）
- 进 FAD：$K_r=(K_I^P+K_I^S)/K_{mat}+\rho$ 作纵坐标，配 Annex P 的 $L_r$ 完成评定。（必做）
对沿厚度非线性的应力（如焊接残余应力廓线），用 5 阶多项式 + 权函数解（BS 7910:2019, §M.4.2 / §M.4.4）：$K_I=\sqrt{\pi a}\sum_{i=0}^{5}\sigma_i(a/B)^i f_i$，各 $f_i$ 由 $a/B$ 与 $2c/a$ 查表。
八、常见误区
- $K_I$ 与应力集中系数 $k_t$ 混淆：$K_I$ 有量纲（$\mathrm{MPa}\sqrt{\mathrm m}$）、描述裂尖奇异场；$k_t$ 无量纲、描述无裂纹几何突变——两码事。
- 只算一个求值点：表面裂纹必须最深点 A、表面点 C 都算，取危险者。
- 二次应力算进 $L_r$：残余 / 热应力只进 $K_r$，绝不进 $L_r$。
- 修正因子乱加：$f_w$ 要占截面比例够大才用、$M_k$ 只在焊趾、$M$ 只对曲壳内压、$k_m$ 只在错边；不满足条件时取 1。
- 超范围外推：解有适用域，越界宁可换权函数 / 有限元，不硬套。
用 MechCalc 在线计算 $K_I$
理解了原理，动手最快。用在线的 BS 7910 Annex M 计算器：选缺陷几何、填裂纹尺寸与应力，它按 Annex M 组装各修正因子、在最深点/表面点算出 $K_I$，并给出白盒推导步骤；支持多项式应力（残余应力廓线）与焊趾 $M_k$。
🧮 在线计算器：BS 7910 Annex M — 应力强度因子计算器 — 选几何、填 a/c/B/W 与 Pm/Pb，即可算出最深点 A 与表面点 C 的 K_I。
本文图示为 MechCalc 原创示意图，仅作教学说明，不替代 BS 7910 规范原文。工程评定请以 BS 7910:2019+A1:2020 现行版为准。