3. 缺陷深度与壁厚剩余强度计算

缺陷深度与剩余强度计算的任务边界

本节针对已经完成尺寸量化的单个腐蚀缺陷(长度 L、深度 d、壁厚 t、外径 D、屈服强度 SMYS),完成两件事:一是按规范选择计算参数并算出失效压力 Pf,二是把 Pf 与最大允许操作压力 MAOP 比较,给出强度校核结论。核心难点不在公式本身,而在于参数怎么取(流变应力、缺陷深度是否含制造公差、是否考虑内压波动)以及结果怎么判读(安全系数落在哪个区间、是否需要降压或修复)。

方法选择与参数取值

工程上常用两条路线,选择依据是缺陷形态与数据质量:

  • ASME B31G(含改进版 RSTRENG):适合中低强度钢、缺陷形状较规则的场合。原始 B31G 用 2/3·SMYS 作流变应力、用抛物线近似金属损失面积,结果偏保守;改进版改用 0.8·SMYS 流变应力与更精确的面积因子,保守度下降。
  • DNV RP-F101(Part B,单缺陷):适合高强度钢、需要更精细评估的场合,采用分项安全系数法,把内压、缺陷深度、材料强度分别乘以分项系数,结果更贴近实际承载能力。

本节示例采用 DNV RP-F101 单缺陷方法,因为它把“参数选择”显式化为分项系数,便于讲清判读逻辑。示例参数如下(均为教学假设值,非实测):

参数符号取值说明
外径D323.9 mm12 英寸管
壁厚t12.7 mm名义壁厚
缺陷深度d3.8 mmILI 量化值,已扣除检测误差上限
缺陷长度L120 mm轴向投影长度
屈服强度SMYS450 MPaX65 钢
抗拉强度UTS535 MPa用于流变应力
运行压力Pop9.5 MPaMAOP

计算过程

DNV RP-F101 单缺陷失效压力公式(Part B,仅内压载荷):

Pf = (2·t·f_u / (D - t)) · [1 - (d/t)] / [1 - (d/t)/Q]
Q  = sqrt(1 + 0.31·(L/sqrt(D·t))²)
f_u = UTS(分项系数法中用抗拉强度,不用屈服强度)

代入示例数据:

D·t = 323.9 × 12.7 = 4113.5 mm²,sqrt = 64.14 mm
L/sqrt(D·t) = 120 / 64.14 = 1.871
Q = sqrt(1 + 0.31 × 1.871²) = sqrt(1 + 1.085) = 1.443
d/t = 3.8 / 12.7 = 0.2992
2·t·f_u/(D-t) = 2 × 12.7 × 535 / (323.9 - 12.7) = 13589 / 311.2 = 43.67 MPa
Pf = 43.67 × (1 - 0.2992) / (1 - 0.2992/1.443)
   = 43.67 × 0.7008 / 0.7927
   = 38.61 MPa

得到无分项系数的失效压力约 38.6 MPa。接下来按 DNV 分项系数法做设计校核,取内压分项系数 γp=1.05、缺陷深度分项系数 γd=1.20、材料强度分项系数 γu=1.15(示例取值,实际应按规范表格与目标安全等级选取):

P_safe = Pf / (γ_p · γ_d · γ_u)
       = 38.61 / (1.05 × 1.20 × 1.15)
       = 38.61 / 1.449
       = 26.6 MPa

把 P_safe 与 MAOP = 9.5 MPa 比较:安全裕度 = 26.6 / 9.5 ≈ 2.8,远大于 1,说明该缺陷在当前尺寸下不构成强度威胁。若把缺陷深度改为 7.6 mm(d/t=0.6),重算可得 Pf≈22.4 MPa、P_safe≈15.5 MPa,裕度降到 1.6,仍合格但已进入需要加密监测的区间;若 d/t 超过 0.8,多数规范直接判定为不合格,不再走公式。

参数选择三个易错点:① 流变应力用 UTS 还是 SMYS 取决于所选方法,B31G 原始版用 2/3·SMYS,DNV 分项系数法用 UTS,混用会系统性偏差;② 缺陷深度必须取“检测值 + 检测误差上限”,否则低估风险;③ 长度 L 应取轴向投影长度,环向缺陷不适用本公式。
单缺陷剩余强度计算与判读链路 输入参数 D = 323.9 mm t = 12.7 mm d = 3.8 mm(含误差) L = 120 mm UTS = 535 MPa MAOP = 9.5 MPa d/t = 0.299 L/√(Dt) = 1.871 DNV RP-F101 计算 Q = √(1+0.31·(L/√Dt)²) = 1.443 Pf = 2t·UTS/(D−t) ×(1−d/t)/(1−d/t/Q) Pf ≈ 38.6 MPa 分项系数:γp=1.05 γd=1.20 γu=1.15 设计安全压力 P_safe = Pf / (γp·γd·γu) = 38.6 / 1.449 ≈ 26.6 MPa 裕度 = 26.6/9.5 ≈ 2.8 判读 裕度 > 2 → 合格 1 < 裕度 < 2 → 监测 裕度 < 1 → 降压/修复 缺陷深度 d/t 对安全压力与裕度的影响(L=120 mm 固定) 缺陷深度比 d/t 压力 (MPa) 0.2 0.4 0.6 0.8 1.0 0 20 40 MAOP 9.5 Pf 失效压力 P_safe 设计安全压力 当前 d/t=0.30 d/t > 0.8 区域 规范判定不合格

结果判读与工程动作

把计算结果落到决策上,需要区分三个层次:

  1. 强度是否合格:P_safe ≥ MAOP 即合格。本例 26.6 MPa 远高于 9.5 MPa,结论为合格。
  2. 裕度落在哪个区间:裕度 > 2 属正常;1~2 之间说明缺陷已接近临界,应缩短复检周期、纳入重点监测清单;< 1 则必须降压运行或安排修复。
  3. 参数敏感性:从上方曲线可见,d/t 从 0.3 增至 0.6 时 P_safe 从约 26.6 MPa 降到约 15.5 MPa,降幅接近 42%。这意味着深度量化的误差会被放大,因此 ILI 深度检测误差上限必须计入 d,而不能直接用读数。

若同一管段存在多个缺陷,本节只做单缺陷评估;相邻缺陷是否需要组合成复合缺陷、组合后等效尺寸如何取,属于交互作用判定的范畴,不在本节展开。

快速复核口诀:先算 d/t 与 L/√(Dt) 两个无量纲量,再套 Q 因子,最后乘分项系数。若手算结果与软件差 10% 以上,优先检查流变应力取值和 d 是否含检测误差,而不是怀疑公式。