🧮 直接使用计算器:跳过推导?进入 BS 7910 断裂韧性估算计算器 ,加载标准算例一键验证。


1. 为什么需要从 CVN 估算 K_mat?

在结构完整性评估(FFS/ECA)中,断裂韧性 $K_{mat}$(Fracture Toughness)是评判裂纹是否会发生脆断的核心参数。然而在大量工程实际场景中:

  • 构件已在役多年,无法补测断裂韧性试样
  • 规范仅要求 Charpy 冲击试验,历史数据仅有 CVN 值
  • 全尺寸断裂韧性试验成本极高,工程上不可行

此时,间接转换成为唯一可行的路径。BS 7910:2019 Annex J 正是为此目的提供了两条经过充分验证的转换方法,专用于铁素体钢(Ferritic steels)。

⚠️ 本方法仅适用于铁素体钢(屈服强度 $\sigma_{ys} \leq 690$ MPa),不适用于奥氏体不锈钢(无韧脆转变行为)。


2. 两条转换路径概览

路径 规范条款 适用条件 特点
A:下平台公式 BS 7910:2019, Annex J.2.1 3 J < $C_v$ ≤ 27 J,断口结晶度 ≥ 80% 保守,适用于极低温脆性区
B:Master Curve BS 7910:2019, Annex J.2.2 韧脆转变区 统计物理模型,更精确

3. 路径 A:下平台公式 (Annex J.2.1)

3.1 适用前提

下平台公式要求材料处于"接近下平台"状态,表现为:

  • Charpy 冲击功 $C_v$ 在 3 J 到 27 J 之间
  • 断口中的解理(脆断)面积比例(结晶度,Crystallinity)≥ 80%

当 $C_v > 27$ J 时,材料已进入韧脆转变区,J.2.1 公式将高估韧性,此时必须改用 Master Curve 方法。

3.2 计算公式

$$K_{mat} = \left[(12\sqrt{C_v} - 20) \cdot \left(\frac{25}{B}\right)^{0.25}\right] + 20 \quad [\text{MPa}\sqrt{\text{m}}]$$

(来源:BS 7910:2019, Annex J.2.1, Eq J.1)

各项物理含义

  • 常数 20:铁素体钢理论最低断裂韧性 $K_{min}$ = 20 MPa√m(渐近下限,来自 ASTM E1921)
  • $12\sqrt{C_v}$ 项:大量铁素体钢试验数据统计拟合的下包络线经验系数
  • $(25/B)^{0.25}$:最弱链理论(Weakest-link Theory)尺寸修正——更厚的截面包含更多微裂纹萌生位点,统计上 $K_{mat}$ 偏低;指数 0.25 = 1/β,β = 4 为 Weibull 形状参数

3.3 验证算例 2

输入:$C_v = 20$ J,$B = 50$ mm

$$K_{mat} = [(12 \times \sqrt{20} - 20) \times (25/50)^{0.25}] + 20 = [33.66 \times 0.8409] + 20 = \mathbf{48.3 \text{ MPa}\sqrt{\text{m}}}$$

4. 路径 B:Master Curve 方法 (Annex J.2.2)

Master Curve 方法由 Kim Wallin(1984)基于大量铁素体钢断裂韧性试验数据提出,核心思想是:铁素体钢在韧脆转变区的断裂韧性分布遵循三参数 Weibull 分布,且分布形状具有材料无关性(仅位置参数 $T_0$ 因材料而异)。

4.1 第一步:确定参考转变温度 $T_0$

方式 1(推荐):若已有 ASTM E1921 实测 Master Curve 数据,直接使用实测 $T_0$,此时无需统计裕度。

方式 2:由 $T_{27J}$ 估算(本计算器默认路径):

$$T_0 = T_{27J} - 18 \quad [°\text{C}]$$

(来源:BS 7910:2019, Annex J.2.2)

$T_{27J}$ 是 CVN 转变曲线上 $C_v = 27$ J 对应的温度,经验偏移 −18°C 得到 $T_0$。

4.2 第二步:T_K 置信裕度(⚠️ 关键易错点)

当 $T_0$ 由 $T_{27J}$ 估算时,BS 7910 规定必须附加 $T_K = 25°\text{C}$ 置信裕度,以覆盖 CVN → $T_0$ 转换过程的统计散布(90% 置信度):

$$T_{eff} = T - T_0 - T_K = T - T_0 - 25°\text{C}$$

这是最常见的计算错误来源!若忽略 $T_K$,将高估 $K_{mat}$。在标准算例 1 中,忽略 $T_K$ 时 $K_{mat}$ 将从 90.8 MPa√m 虚高到约 128 MPa√m。

4.3 第三步:Master Curve 计算

Weibull 尺度参数 $K_0$(63.2% 失效分位数,对应 1T 标准试样):

$$K_0 = 31 + 77 \cdot \exp[0.019 \cdot T_{eff}] \quad [\text{MPa}\sqrt{\text{m}}]$$

(来源:BS 7910:2019, Annex L.9.5.3, Eq L.13;31 = $K_{min}$ + 11 = 20 + 11)

指定失效概率下的 $K_{mat}$(1T 基准):

$$K_{mat,1T} = 20 + \left[\ln\left(\frac{1}{1-P_f}\right)\right]^{0.25} \cdot (K_0 - 20)$$

(来源:BS 7910:2019, Annex J.2.2, Eq J.5;ASTM E1921 三参数 Weibull 反函数)

当 $P_f = 0.05$:概率项 = $[ln(1/0.95)]^{0.25} = (0.05129)^{0.25} = 0.4759$

厚度修正(最弱链尺寸效应):

$$K_{mat} = 20 + (K_{mat,1T} - 20) \cdot \left(\frac{25}{B}\right)^{0.25}$$

(来源:BS 7910:2019, Annex J.2.2;ASTM E1921)

4.4 验证算例 1(标准算例)

输入:$T_{27J} = -50°\text{C}$,$T = 0°\text{C}$,$B = 60$ mm,$P_f = 0.05$

步骤 计算 结果
B1: T_0 = T_27J - 18 = -50 - 18 -68°C
T_eff (含 T_K=25°C) 0 - (-68) - 25 43°C
K_0 = 31 + 77·exp(0.019×43) = 31 + 174.3 205.3 MPa√m
K_mat_1T = 20 + 0.4759×(205.3-20) = 20 + 88.2 108.2 MPa√m
K_mat = 20 + 88.2×(25/60)^0.25 = 20 + 70.9 90.9 MPa√m ✅

规范参考值:90.8 MPa√m(BS 7910:2019, Annex J.2.5 Note 2)


5. API 579 对比路径

API 579-1:2021, 9F.4.3.2, Eq.(9F.74) 提供了一个更精细的多参数 $T_0$ 估算公式:

$$T_0 = T_{28J} - 79 + \frac{\sigma_{ys}}{9} - \frac{C_{V-US}}{59} + \Delta T_0 \quad [°\text{C}]$$

其中 $\Delta T_0 = 18°\text{C}$(固定保守裕度),$T_{28J} \approx T_{27J}$(27 J ≈ 20 ft-lbs)。

与 BS 7910 的 $T_0 = T_{27J} - 18°\text{C}$ 相比,API 579 额外考虑了:

  • 屈服强度项 $+\sigma_{ys}/9$:强度越高→脆性倾向越大→T_0 向高温偏移
  • 上平台能量项 $-C_{V-US}/59$:延性越好→T_0 向低温偏移

注意:API 579 的 $T_0$ 视为"直接测定值",代入 Master Curve 时 $T_K = 0$($\Delta T_0 = 18°\text{C}$ 已起到同等作用)。

算例 3 对比结果

参数 BS 7910 API 579
T_0 -18°C -24.1°C
K_mat 44.7 MPa√m 61.0 MPa√m

API 579 因额外考虑屈服强度和上平台能量,T_0 更低,K_mat 更高(更乐观)。


6. 工程使用建议

  1. 优先顺序:实测 $T_0$(ASTM E1921)> API 579 多参数转换 > BS 7910 单参数 $T_{27J}$ 估算
  2. 必须检查 $T_K = 25°C$ 是否已代入(CVN 路径专用)
  3. 下平台边界:$C_v > 27$ J 时禁止使用 J.2.1,强制改用 Master Curve
  4. 厚度效应:当构件壁厚 $B > 25$ mm 时,尺寸修正使 $K_{mat}$ 显著下降,不可忽略

7. 快速上手计算器

  1. 进入 断裂力学计算器 → 选择"BS 7910 断裂韧性估算"
  2. 点击"算例 1 (Master Curve, BS 7910 标准算例)“按钮自动填入参数
  3. 点击”计算“验证结果:$K_{mat} = 90.85$ MPa√m ✅
  4. 右侧面板显示 Master Curve 三条概率曲线($P_f$ = 5%/50%/95%)及当前评定点位置

📖 参考引用:

  • BS 7910:2019+A1:2020, Guide to methods for assessing the acceptability of flaws in metallic structures
  • API 579-1/ASME FFS-1:2021, Part 9 Annex 9F — Fracture Toughness
  • ASTM E1921, Standard Test Method for Determination of Reference Temperature T_0
  • Wallin, K. (1984). The scatter in K_Ic results. Engineering Fracture Mechanics, 19(6), 1085-1093.