应力强度因子 模拟器 返回
断裂力学模拟器

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

I型应力强度因子 K_I=Yσ√(πa) 实时计算。通过滑块调整形状因子、名义应力、裂纹长度、断裂韧性,在裂纹示意图和K_I-a曲线图中可视化安全系数、许用裂纹长度、断裂传播应力。

材料·形状预设

计算结果
应力强度因子 K_I
安全系数 SF
许用裂纹长度 a_crit
断裂传播应力 σ_crit
裂纹模型示意图(应力集中可视化)
K_I vs 裂纹长度 a — 断裂韧性线与曲线交点为 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 从 5mm 变为 10mm,你会看到 K_I 从约 28 跳到 40 MPa·√m。
🙋
那么当「断裂韧性 K_IC」被超过时,材料就会断裂?当 SF 低于 1 时就危险了?
🎓
完全正确。当 K_I ≥ K_IC 时,会发生不稳定破裂(一下子就断了)。SF = K_IC/K_I 是安全系数。在实际工程中,SF = 1 太危险了,通常要求 SF = 2~3。试着用工具把应力从 200 变到 300 MPa,你会看到 SF 从 1.78 降到 1.19。同一个裂纹,应力越高危险性越大。
🙋
「许用裂纹长度 a_crit」是什么意思?是说在检查中发现的裂纹长度可以到多少还安全?
🎓
完全就是这个意思。a_crit = (K_IC/(Yσ))²/π 表示「在当前应力下,裂纹长度达到多少就会破裂」。默认条件下 200MPa·Y=1.12 时,大约是 15.9mm 的破断界线。在实务中,结合检查周期,通常规定「下次检查前裂纹增长也不超过 a_crit/2」作为管理标准。图上你可以看到红色水平线(K_IC)与蓝色曲线(K_I 曲线)的交点,那就是 a_crit。
🙋
形状因子 Y 在边缘裂纹时是 1.12,中央裂纹时是 1.0,为什么边缘的更大呢?
🎓
因为自由表面(边界)存在时,裂纹尖端的应力只能向一侧释放,应力场就更强。这叫「自由表面修正」,半无限平板边缘裂纹的经典解就是 Y≈1.1215。压制品的角或焊接缝端部实际上相当于边缘裂纹,设计上采用 Y=1.12 就能得到更安全的评估。

常见问题

理论上 SF≥1 表示不会破裂,但在实际工程中存在荷载波动、温度变化、腐蚀、材料偏差、裂纹检测精度等不确定性,所以设计安全系数通常要求 SF≥2~3。压力容器规范(ASME、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 积分方法。
可以查 Tada/Paris/Irwin 的《应力分析裂纹手册》、BS 7910、NASGRO 等数据库,里面有各种试样形状和裂纹位置的 Y(a/W) 解。不清楚时,保守做法是取 Y=1.5 左右,或者用有限元方法直接计算 K_I。

实际应用

压力容器·管道完整性评估:定期检查发现焊接缺陷、SCC 或腐蚀坑时,计算运行压力下的 K_I,与规范 K_IC 比较,判断下次检查前的安全性。这是 ASME BPVC Section XI 断裂力学评估的核心。

航空器结构损伤容限设计:对机翼或机身蒙皮的「设想初始缺陷」计算 K_I,结合 Paris 准则估计到达许用裂纹长度所需的飞行周数。检查间隔设为该周数的 1/2。

核电站反应堆压力容器低温脆性评估:考虑照射后 K_IC 下降,计算 PTS(加压热冲击)工况下 K_I 是否超过低温脆性迁移后的 K_IC。这是 PTS 分析的关键。

桥梁·储罐维护管理:焊接部位发现疲劳裂纹时,比较当前 K_I 与 K_IC,当 SF<2 时列为补修对象,优先级高。

常见误解和注意事项

最常见的误解是「σ 是裂纹尖端的局部应力」。这是错的。σ 应该是远离裂纹的名义应力(平均应力或无缺陷时的应力),不是应力集中系数 K_t 乘以的局部应力。应力集中(缺口)和应力强度扩大(裂纹)是不同的概念。

其次是「形状因子 Y 的选择错误」。本工具直接设置 Y,但实际结构中 Y 往往随 a/W(裂纹长与部材宽比)变化。裂纹扩展时 Y 会增大,K_I 增长就超过单纯的 √a 规律,会急速加速。必须查表确认有限宽度补正。

最后是「K_IC 不是单一常数」。它随温度、板厚、应变速率、环境(氢、海水)大幅变化。比如低于脆性转变温度时 K_IC 急剧下降,同一部件在室温安全但冬季破裂的事故曾发生过。本工具数值代表「室温·大气·足厚平面应变」条件,实际运用必须根据环境选用材料数据。

使用指南

  1. 调整「名义应力 σ」滑块(0~500 MPa 范围),输入部材承受的拉伸应力
  2. 改变「裂纹长度 a」滑块(0.1~50 mm 范围),设置既存缺陷的尺寸
  3. 操作「几何因子 Y」滑块(0.5~3.0 范围),反映裂纹位置(表面裂纹、贯通裂纹)导致的应力集中
  4. 调整「断裂韧性 K_IC」滑块(20~100 MPa√m),输入材料的耐裂纹性能
  5. 查看实时输出的「应力强度因子 K_I」「安全系数 SF」「许用裂纹长度 a_crit」「断裂传播应力 σ_crit」

具体计算示例

以铝合金 7075-T73 航空器翼面为例验证。输入名义应力 σ=250 MPa、既存裂纹长度 a=5 mm、几何因子 Y=1.12(表面裂纹)、断裂韧性 K_IC=35 MPa√m,则 K_I = σ × Y × √(π×a) = 250 × 1.12 × √(π×5) = 43.8 MPa√m,安全系数 SF = K_IC / K_I = 35 / 43.8 = 0.80。断裂传播应力为 σ_crit = K_IC / (Y×√(π×a_crit)),许用裂纹长度 a_crit = (K_IC/(Y×σ))² / π = 2.94 mm。此时 SF<1.0,需要更换部件。

实务注意事项