什么是应力强度因子 K_I？

应力强度因子 K_I 是线弹性断裂力学中用来表征裂纹尖端应力场强度的参数，定义为 K_I = Y σ √(πa)。Y 是形状（几何）因子，σ 是远离裂纹处的名义应力，a 是裂纹长度。当 K_I 达到材料的断裂韧度 K_IC 时即发生失稳断裂。

形状因子 Y 如何确定？

Y 取决于试件几何、载荷形式与裂纹位置。对无限大板中央贯穿裂纹 Y = 1.0；对半无限大板边缘裂纹 Y ≈ 1.12（自由表面修正）。对于有限宽板与曲面结构，应查阅 SIF 手册（Tada/Paris/Irwin）的 Y(a/W) 表或使用有限元法直接计算。

如何使用临界裂纹长度 a_crit 与断裂应力 σ_crit？

a_crit = (1/π)(K_IC/(Y σ))^2 是当前应力下能引起断裂的裂纹长度，用作检验管理限值。σ_crit = K_IC/(Y √(πa)) 是当前裂纹长度下导致断裂的应力。两者结合并取裕度（如 a_crit/2）作为运行管理基准。

K_I 与 Paris 公式有何关系？

K_I 用于静力（一次性）断裂评估，超过 K_IC 即发生断裂。Paris 公式 da/dN = C (ΔK)^m 用于循环载荷下裂纹扩展，输入量为应力强度因子幅 ΔK = K_max − K_min。二者在同一 K 框架下，分别对应静力断裂评估与疲劳寿命评估。

应力强度因子模拟器 — 线弹性断裂力学

几何与材料预设

名义应力 σ 200 MPa

裂纹长度 a 5.0 mm

形状因子 Y 1.12

断裂韧度 K_IC 50 MPa·√m

计算结果

—

应力强度因子 K_I

—

安全系数 SF

—

临界裂纹长度 a_crit

—

断裂应力 σ_crit

裂纹示意图（应力集中可视化）

K_I 与裂纹长度 a 的关系 — 与 K_IC 交点即 a_crit

理论与主要公式

$$K_I = Y\,\sigma\,\sqrt{\pi a}$$

Y 为形状因子（中央贯穿裂纹取 1.0，边缘裂纹取 1.12），σ 为远端名义应力 (MPa)，a 为裂纹长度 (m)。

$$\mathrm{SF} = \frac{K_{IC}}{K_I},\quad a_{crit} = \frac{1}{\pi}\left(\frac{K_{IC}}{Y\,\sigma}\right)^2,\quad \sigma_{crit} = \frac{K_{IC}}{Y\sqrt{\pi a}}$$

当 K_I ≥ K_IC 时发生失稳断裂。SF<1 表示即时断裂；a_crit 是当前应力下的临界裂纹长度；σ_crit 是当前裂纹长度下导致断裂的应力。

应力强度因子模拟器是什么

🙋

应力强度因子是什么意思？K_I 和普通的应力 σ 到底有什么区别？听说只要有裂纹一切都会改变。

🎓

很棒的问题。在尖锐裂纹尖端，线弹性理论给出的应力会发散到无穷大，所以单靠 σ 就没法描述。我们用 K_I = Y σ √(πa) 来衡量这个奇异场的强度。a 是裂纹长度，Y 是几何因子。把工具里的 a 从 5 mm 滑到 10 mm 看看：K_I 会从大约 28 跳到 40 MPa·√m。

🙋

那只要 K_I 超过断裂韧度 K_IC 就会断裂吗？安全系数 SF 小于 1 就意味着已经出事？

🎓

正是如此。K_I ≥ K_IC 时发生失稳断裂，结构会瞬间裂开。SF = K_IC/K_I 就是安全系数。工程上不会让 SF 接近 1，一般要求 2～3。把 σ 从 200 滑到 300 MPa 试试，SF 会从 1.78 掉到 1.19。裂纹一样，但风险完全不同。

🙋

那临界裂纹长度 a_crit 呢？是不是表示「在役条件下还能容忍的最长裂纹」？

🎓

对。a_crit = (K_IC/(Y σ))^2 / π 是当前应力下能引起断裂的裂纹长度。在 σ = 200 MPa、Y = 1.12 时大约 15.9 mm。工程上会把检验管理限值定为 a_crit/2 或更小，并按裂纹扩展速度设定检验周期。图中红色虚线就是 K_IC，与蓝色 K_I 曲线的交点正好就是 a_crit。

🙋

为什么边缘裂纹的 Y 是 1.12，而中央裂纹的 Y 才 1.0？

🎓

边缘裂纹紧贴自由表面，应力没法对称地绕过去，裂纹尖端的奇异性会更强，这就是「自由表面修正」。半无限大板边缘裂纹的经典解给出 Y ≈ 1.1215。焊趾、机加工棱角和冲压角通常都按边缘裂纹处理，工程上用 Y = 1.12 是一种现实的保守取值。

常见问题

理论上 SF ≥ 1 即可避免断裂，但实际部件会面临载荷波动、温度变化、腐蚀、材料数据分散和裂纹检测精度不确定性。设计中一般要求 SF = 2～3。ASME BPVC 与 JIS B 8265 等压力容器标准多采用 SF = 3。

按 ASTM E399 或 JIS Z 2284，使用紧凑拉伸（CT）或三点弯曲试件在平面应变条件下测定 K_IC。试件较薄时塑性区扩展，得到的 K_C 偏大；当试件足够厚时收敛得到的就是 K_IC（平面应变断裂韧度）。薄板设计应取更保守的数据。

本工具基于线弹性断裂力学（LEFM），假定裂纹尖端塑性区相对 a 足够小。当塑性区超过 a/50 时，应使用 Irwin 修正（有效裂纹长度 a + r_y），或直接采用弹塑性断裂力学（J 积分、CTOD）方法。

可查阅 Tada/Paris/Irwin 的《Stress Analysis of Cracks Handbook》或 BS 7910、NASGRO 等数据库，其中列有各种试件几何和裂纹位置下的 Y(a/W) 解。无数据时可保守取 Y ≈ 1.5，或直接用有限元法计算 K_I。

实际应用

压力容器与管道完整性评估：定期检验中发现的焊接缺陷、SCC、腐蚀坑等，根据工作压力计算 K_I 并与规范 K_IC 比较，判断到下次检验前是否安全。LEFM 是 ASME BPVC Section XI 的基础评估框架。

飞机结构的损伤容限设计：从「假设初始缺陷尺寸」算 K_I，结合 Paris 公式估算到达 a_crit 所需的飞行循环数，再把检验间隔设为该寿命的 1/2。

核电压力热冲击（PTS）评估：反应堆压力容器（RPV）经过中子辐照后 K_IC 下降，PTS 分析需验证冷水注入瞬态过程中的 K_I 不会超过低温下的 K_IC。

桥梁与储罐运维管理：对焊接节点发现的疲劳裂纹，比较 K_I 与 K_IC，对 SF<2 的部位优先采取止裂孔、贴板补强等修复措施。

常见误区与注意事项

最常见的误解是「σ 是远场名义应力，不是裂纹尖端的局部应力」。不要再乘以应力集中系数 K_t。应力集中适用于钝缺口，应力强度适用于尖锐裂纹，两者概念不同，重复计入会严重高估 K_I。

其次形状因子 Y 通常并非常数。本工具用滑块直接设定 Y，但实际结构中 Y 是 a/W（裂纹长度/部件宽度）的函数。裂纹增长后 Y 也增长，K_I 比单纯的 √a 律加速更快，必须考虑有限宽修正。

最后K_IC 不是单一的材料常数。它随温度、板厚、应变速率和环境（氢、海水）显著变化。在韧脆转变温度以下 K_IC 会急剧降低，二战时期 Liberty 货船和密尔沃基 Hoan 桥的脆性断裂事故就是惨痛教训。务必使用与实际工作条件相符的 K_IC 数据。