应力强度因子

断裂评定的第一步，是算出裂纹尖端的应力强度因子 $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 纵轴的原料。 ...

在 [[bs7910-a533b-residual-stress-fad|A533B 焊接平板四道题]]的 FAD 评定里，焊态件那一项残余应力强度因子 $K_I^S \approx 46\ \mathrm{MPa\cdot m^{0.5}}$ 一直是直接输入的——它由 FITNET 原文把实测残余应力廓线积分得到，我们只是把这个现成数喂进纵坐标 $K_r$。 一个自然的追问：这个 46 到底是怎么从一条残余应力曲线变出来的？mechCalc 自己能不能算？ 能。这正是 BS 7910 Annex M.4.2（平板有限表面多项式裂纹） 解干的活：给定一条沿壁厚分布的应力多项式 $\sigma(x)$，把它对裂纹前缘做权函数积分，得到该裂纹的应力强度因子。本文就用 mechCalc 的 Annex M 计算器单独跑这一步，把焊态残余廓线积分成最深点 SIF，再和 FITNET 的 46 对一对——这是一次纯粹的应力强度因子单点计算，与 FAD 评定无关。 原理：多项式应力廓线 → SIF 被评定对象是焊缝区的半椭圆表面裂纹（深度 $a$、表面半长 $c$）。二次应力（焊接残余）沿壁厚不是常数，而是一条曲线，BS 7910 用一个不超过五次的多项式来描述它（坐标 $x$ 从开裂表面起算、以壁厚 $B$ 归一）： $$ \sigma(x) = \sum_{n=0}^{5} \sigma_n \left(\frac{x}{B}\right)^n $$最深点 $K_I$ 按 Annex M.4.2.2 的 Eq. M.13 装配——每一阶应力分量配一个 Fett 几何函数 $f_i^d$（查 Table M.1，按 $a/B$ 与 $a/2c$ 双线性插值），求和后乘 $\sqrt{\pi a}$： ...

在压力容器、石油管道、海洋钢结构以及核电站承压构件的结构完整性评定中，应力强度因子 (Stress Intensity Factor, SIF, 记为 $K_I$) 是线弹性断裂力学 (Linear Elastic Fracture Mechanics, LEFM) 的核心物理量。它精确表征了裂纹尖端奇异应力场的强弱程度。 英国焊接研究所 (The Welding Institute, TWI) 颁布的缺陷评定标准 BS 7910:2019+A1:2020 在其附录 M (Annex M) 中，针对工程中常见的八大几何分组，整理出了一套极其经典、高精度且自洽的应力强度因子解析手算求解体系。本博文剖析其物理本质、规范计算链与工程应用，并介绍完全基于纯 Python 开发的交互式在线计算器。 1. 物理本质：通用应力强度因子叠加合成模型 BS 7910 附录 M 采用基于叠加原理的通用应力强度因子合成模型 (BS 7910:2019, Annex M, Eq. M.1)： $$ K_I = Y\sigma \sqrt{\pi a} $$更具体地，为了同时涵盖一次载荷 (Primary Stress) 和二次载荷 (Secondary Stress，如残余应力) 以及宏观和微观的局部应力集中效应，规范将该公式细化展开为 (BS 7910:2019, Annex M, Clause M.2.1) 形式： $$ K_I = \left[ M \cdot f_w \left( k_{tm} \cdot M_{km} \cdot M_m \cdot P_m + k_{tb} \cdot M_{kb} \cdot M_b \cdot (P_b + (k_m - 1) P_m) \right) + (M_m Q_m + M_b Q_b) \right] \sqrt{\frac{\pi a}{Q}} $$核心物理乘子链解析： $P_m, P_b$ 与 $Q_m, Q_b$：一次薄膜/弯曲应力与二次薄膜/弯曲应力 (Clause M.2.1)。 $M_m, M_b$：无界板几何修正因子 (Membrane & Bending Correction Factors)。 $f_w$：有限截面宽度修正系数 (Finite Width Correction Factor)，如穿透裂纹的割线宽度修正 (Eq. M.8)。 $Q$：缺陷形状参数 (Crack Shape Parameter, $Q = \Phi^2$)，用以将直线开裂修正为半椭圆或全椭圆奇异前缘 (Eq. M.11)。 $M$：薄壳拉伸大跨度物理鼓胀因子 (Bulging Factor, 亦称 Folias 因子)，用于补偿圆筒管道和弯曲壳体的壳张力释放 (Clause M.6)。 $k_m, k_{tm}, k_{tb}$：制造错边 (Misalignment) 与宏观几何应力集中修正链。 $M_{km}, M_{kb}$：焊缝趾部 (Weld Toe) 微观物理缺口应力集中因子，用于捕捉焊脚几何的微观应力奇异性 (Clause M.4.1.2)。 2. 经典 Newman-Raju 平板半椭圆解 (M.4) 对于工程中最具有代表性的平板表面半椭圆裂纹 (Surface Semi-elliptical Crack)，BS 7910 Annex M.4 完全采纳了 Newman-Raju (1986) 多项式拟合解。该解是断裂力学发展史上的里程碑，它首次实现了对缺陷前缘最深点 (Deepest Point，记为 A 点) 与表面自由端点 (Surface Point，记为 C 点) 双前缘特征的并行精确评估： ...

含缺陷结构的断裂评定要同时看两件事——离脆断多近、离塑性失稳多近。本文用一张失效评定图（FAD）把断裂力学的评定原理、计算链路与 BS 7910 的评定步骤串成一条清晰的主线。读完你应能看懂一个评定点是怎么算出来的、又是怎么判定合格与否的。 引子：一条裂纹被发现之后 在役的压力容器、管道或焊接结构，检验时常会发现裂纹一类的缺陷。此刻工程师要回答的，不是「这条裂纹能不能用」——裂纹本身不会被「使用」——而是：含着这条缺陷的结构，还能不能在当前载荷下继续安全服役？ 这正是合于使用评价（Fitness-for-Service, FFS）要解决的问题。BS 7910:2019 的第 7 章（Clause 7）给出了其中最核心的一环——断裂评定：判断一个已知尺寸的缺陷，在给定载荷与材料下是否可接受。 传统的强度设计（「工作应力 < 许用应力」）在这里失灵了。因为按线弹性理论，裂纹尖端的应力趋于无穷大，根本没法直接和许用应力比较。断裂力学换了一套语言来回答这个问题。 一、两种失效的角力：为什么需要一张"评定图" 含缺陷结构可能沿两条完全不同的路径失效（BS 7910:2019, §7.1.1）： 脆性断裂（脆断）：裂纹尖端的驱动力超过材料的断裂韧性，裂纹突然、快速扩展，几乎无预兆。低温、高强度、厚截面材料尤其危险。 塑性失稳（塑性垮塌）：含缺陷的剩余承载截面整体进入屈服、失去承载能力而垮塌。高韧性、延展性好的材料倾向于此。 真实结构往往介于两者之间：材料越脆越偏向前者，越韧越偏向后者，而在韧脆转变区两种风险并存。单一判据覆盖不了这种竞争——只看韧性会漏掉垮塌，只看强度会漏掉脆断。 BS 7910 的解法是用一张二维图同时刻画这两种失效，这张图就是失效评定图（Failure Assessment Diagram, FAD）。 二、一张图看懂断裂评定：失效评定图（FAD） 先看整张图的全貌，后面各节都是在把这张图讲透： 图1：BS 7910 断裂评定全景。左侧是计算链——三个已验证的构件 K_I（Annex M）、σ_ref（Annex P）、K_mat（Annex J）装配出横坐标 L_r 与纵坐标 K_r，得到一个评定点 (L_r, K_r)；右侧是失效评定图 FAD——评定点落在评定曲线 FAL 之内（绿色安全域 SAFE）即缺陷可接受，落在曲线之上或之外（红色 UNSAFE）即不可接受。 这张图的两根轴各自量一种失效，评定点越靠近原点越安全： 纵轴 $K_r$（断裂比）：衡量离脆断有多近。 横轴 $L_r$（载荷比）：衡量离塑性失稳有多近。 失效评定曲线 FAL（蓝线）：一条 $K_r = f(L_r)$ 的曲线，把两种失效的竞争合成为一条边界。评定点在曲线内 → 缺陷可接受；在曲线上或曲线外 → 不可接受。 三、断裂力学评定原理 纵轴 $K_r$：离脆断有多近 $$K_r = \frac{K_I}{K_{mat}}$$ $K_I$ 是应力强度因子，度量裂纹尖端应力奇异场的强度——它是裂纹的驱动力。 $K_{mat}$ 是材料的断裂韧性，度量材料抵抗裂纹扩展的能力——它是抗力。 因此 $K_r$ 是「驱动力 / 抗力」。$K_r = 1$ 意味着驱动力恰好等于抗力，即临界脆断。（BS 7910:2019, §7.1.1） 横轴 $L_r$：离塑性失稳有多近 $$L_r = \frac{\sigma_{ref}}{\sigma_Y}$$ $\sigma_{ref}$ 是参考应力，代表含缺陷截面所承受的「等效应力」水平。 $\sigma_Y$ 是材料的屈服强度。 $L_r = 1$ 表示缺陷截面整体达到屈服、开始塑性失稳。注意这里是含缺陷的净截面先失效，而非整根结构一起屈服——这正是「净截面屈服」与「塑性失稳」的区别所在。（BS 7910:2019, §7.1.1） 失效评定曲线（FAL）：把两种失效合成一条线 FAL 不是两根轴的简单拼接，而是一条随 $L_r$ 升高而下压的曲线： ...