氢脆之所以难以预测，在于它并不是一种单一的失效方式，而是一个随氢浓度变化的连续谱。在不少高强钢中，低氢水平下的断口仍能看到韧窝，裂纹依靠微孔形核、长大和聚合推进，但起裂阻力和扩展阻力已经被削弱。氢浓度继续升高，沿晶断裂才逐渐占主导。同一种材料因此可能跨过延性、混合和脆性三种状态。仅仅描述界面脱黏的模型解释不了前一段，只追踪微孔演化的模型又难以跨过后面的模式转变。

准确预测氢脆，是氢气输送与储存系统完整性评价的前提。然而氢脆研究最关键的短板，恰恰是缺少预测模型。原因在于HELP、HESIV与HEDE这些基础机制本身与断裂之间并没有直接的定量联系。本文介绍一篇发表于Corrosion Science的工作，作者以Complete Gurson model为基础构建H-CGM+预测框架，让氢促进的空洞形核与氢致脱黏在同一计算中竞争，再把氢扩散、应力驱动和应变相关陷阱一并结合。

▍论文信息

▍研究背景

氢进入金属后会改变位错活动、空位形成、界面结合和裂纹尖端局部应力场。HELP（氢增强局部塑性）可以让变形集中在更窄的区域，HESIV（氢增强应变诱导空位）会增加空位及空位团簇，HEDE（氢增强脱黏）则降低界面或晶格的分离阻力。三种机制经常同时出现，但在连续介质模型中往往各自起不同作用。尤其在低氢区，试样呈延性断裂而韧性下降，简单采用临界应力随氢降低的脆性判据并不足够。

传统Gurson类模型擅长描述空洞长大和聚合，却常把初始空洞或形核参数当作固定输入，氢对形核阶段的作用因而被弱化。作者据此提出hydrogen-enhanced void nucleation（HEVN，氢增强空洞形核），HELP会造成局部滑移和位错塞积，HESIV会增加空位，两者共同提高第二相颗粒、夹杂或局部弱区附近的空洞形核概率。约束效应决定了这一框架还要回答另一个问题。空气中，深裂纹或弯曲载荷通常带来更高的裂尖三轴应力，测得韧性也更低。含氢后可能出现三种关系：原有排序保持（H1），不同约束下的韧性趋于接近（H2），或者低约束试样的韧性反而更低（H3）。已有AISI 4340数据分别呈现过H1和H2，近期少量结果提示H3也可能存在。后者若被证实，会直接影响用高约束试样数据评估薄壁管线的保守性。本文用机制模型检查HELP、HESIV与HEDE的现有组合能否自行产生这种反转。

图1 氢脆随氢浓度演化的全谱示意与典型断口。

▍研究思路

H-CGM+的骨架来自Complete Gurson model。均匀塑性阶段采用GTN多孔塑性描述空洞长大，局部塑性失稳后用Thomason条件判断空洞聚合。作者把氢对延性损伤的影响放在形核项中，空洞形核速率乘以1+kₙCH，其中 CH为局部总氢浓度，kₙ表示微观组织对氢促进形核的敏感性。kₙ越大，同一氢水平下空洞出现得越早、数量越多，延性断裂阻力随之降低。这种写法避免把HELP简化为材料整体软化，也保留了后续用显微统计或原位实验标定kₙ的可能。

脆性通道采用应力控制的脱黏判据，局部最大主应力达到氢降低后的临界结合强度时，材料点失效。氢覆盖度用Langmuir型关系计算，结合强度退化参考氢/铁体系的DFT结果，参数γ控制材料对HEDE 的敏感程度。这样一来，低氢时HEVN推动微孔断裂，高氢时脱黏判据可能先被触发；最终模式由两条通道在裂尖的竞争决定，而不是预先指定。

氢输运部分以化学势为主变量，包含浓度梯度、静水应力驱动和可逆陷阱。化学势形式不必直接对静水应力求空间梯度，也便于统一处理通量与吸氢边界。作者在Abaqus/Explicit中通过VUMAT和VUMATHT实现耦合，选用一组通用高强铁素体BC管线钢参数，对裂纹深度比a/W=0.1–0.5 的SENT试样开展参数研究。初始晶格氢为0–2wppm，断裂韧性以裂纹扩展0.1 mm 时的CTOD表征。

图2 H-CGM+的失效路径与HEVN概念。

▍主要结果

模型先在材料点层面比较两种断裂通道。无氢时，空洞体积分数持续增长并在局部失稳后完成聚合；加入氢后，形核项被放大，延性损伤更早积累。当裂尖主应力先达到氢退化后的结合强度，计算会转入脱黏型脆性失效。这个竞争关系解释了一个常被忽略的事实，氢致韧性下降并不等同于断口立即变脆。低氢区完全可能仍由微孔聚合控制，只是形核提前、可用塑性耗散减少。

图3 延性空洞演化与脆性脱黏判据的竞争，以及SENT试样与裂尖网格。

kₙ决定低氢区的韧性下降方式。在a/W=0.3 的计算中，kₙ=0 时氢不加速空洞形核，CTOD随氢浓度呈现两个较清楚的平台，随后由HEDE判据触发陡降。kₙ增大后，CTOD从低氢区就开始连续下降，延性向脆性的过渡也被拉宽。换句话说，HEVN不只是降低某个临界值，它改变了整条韧性-氢浓度曲线的形状。该结果与低氢仍有韧窝，但韧性已损失的实验观察相容，不过kₙ的绝对数值仍需针对具体夹杂、第二相和位错结构标定。γ控制另一端的脆性阈值。γ越低，无氢结合强度的基准越低，相同氢覆盖度下的临界阈值也越低，延性向脆性转变因而会在更低氢浓度发生；令γ→∞ 等于关闭脆性判据，此时模型只给出HEVN导致的渐进式延性韧性衰减。kₙ与γ因而分工明确，前者控制氢怎样加快微孔形核，后者控制脱黏何时压过延性通道。

图4 两个关键参数对CTOD-氢浓度关系的影响。

Kumnick-Johnson、Dietzel和Sofronis三种关系预测的裂尖总氢分布差别明显，但本组工况下得到的CTOD曲线几乎重合。作者没有据此宣称陷阱不重要，而是指出初始氢较高、所选结合能较低，可逆陷阱很快趋于饱和，导致不同模型在断裂响应上暂时收敛。降低初始氢、提高陷阱结合能、延长暴露或改变加载速率后，差异仍可能被放大。裂尖氢分布对陷阱模型敏感，并不表示某一特定韧性指标在所有边界条件下都会同样敏感。

当前损伤方程使用晶格氢与陷阱氢之和 CH。这样处理便于把可扩散氢和可逆陷阱中的氢共同纳入局部损伤，但活性氢究竟主要来自晶格、弱可逆陷阱还是特定界面，文献中仍有争论；强陷阱甚至可能降低可参与脆化的氢。论文把这一点列为模型不确定性，因此不能把总氢越高、损伤越大当作脱离陷阱类型和停留时间的普适定律。

图5 三种陷阱关系给出的裂尖氢分布及对应CTOD曲线。

氢浓度扫描清楚展示了模式转换。以a/W=0.3 为例，0–0.25 wppm 时CTOD仅小幅下降，裂纹仍按延性微孔机制扩展；约 0.75 wppm 时，裂尖先出现脆性失效，周围区域仍以延性损伤为主，形成混合模式；升至 1.5–2.0 wppm 后，脆性通道控制更长的裂纹路径，CTOD明显陡降。

约束效应的计算结果总体支持H1。无氢时，a/W越大，裂尖约束越高，起裂CTOD越低；加入 1 wppm 氢后，所有几何的韧性都下降，原有排序仍保持。更值得注意的是，模型给出的延性向脆性转变所需初始晶格氢随约束显著降低，从a/W=0.1 时约 1.0 wppm，降到a/W=0.5 时约 0.1 wppm。中等氢水平下，高约束会更早触发脱黏；当氢足够高、所有试样都进入脆性区后，约束差异对绝对韧性的影响反而减弱，接近 H2。模型没有得到H3，也就是没有出现低约束试样比高约束试样更低韧性的反转。这个负结果说明，在当HELP-HESIV-HEDE表达下，几何约束仍然通过三轴应力主导断裂，单靠应力辅助扩散、HEVN和氢致脱黏不足以制造反转。如果后续实验确认H3稳定存在，就需要寻找低约束条件下特有或更强的局部机制，例如高度局域的氢诱导形核带、更复杂的塑性-输运反馈、微观组织各向异性，或试样尺度与边界充氢的耦合。

图6 氢浓度与裂尖约束共同决定断裂模式和韧性。

▍文章亮点

1.本文引入HEVN框架，将原本定性的物理机制转化为可计算的本构关系。HELP 与 HESIV 长期停留在定性层面，难以直接进入断裂预测。作者把二者的共同后果归结为孔洞形核的加速，用一个形核放大系数加以量化，既给氢压制韧性断裂提供了物理依据，也纠正了此前模型对氢效应的经验式处理。

2. 让延性损伤和脆性脱黏在同一裂尖场中相互竞争。模型不预先规定裂纹必须按哪一种模式扩展，而由空洞聚合与氢退化后的临界结合强度谁先达到决定。因此，延性、混合和脆性断裂可以随氢浓度连续切换，两个参数kₙ和γ的物理职责也相对清楚。

3. 系统扫描约束与氢浓度，而不是事后解释一条曲线。五种a/W、0–2 wppm 初始晶格氢以及两类损伤参数构成了清楚的参数空间。结果同时报告韧性、裂纹路径和延性向脆性转变阈值，揭示高约束会显著下移该阈值。

4. 没有把解释不了的H3填进结论，将负结果讲清楚并转化为问题。模型只能给出H1，并在高氢区趋向H2，无法产生低约束更低韧性的H3。作者明确承认这一边界，并据此提出更彻底的空洞局域化路线。对机制模型而言，能排除一组现有解释同样有价值，它为后续实验应该测什么、理论还缺什么划出了范围。

▍小结

这项工作把氢脆预测从各自为政的韧性模型与脆性模型，推进到一个允许两者竞争的统一框架。其中最具方法学意义的一步，是用氢增强孔洞形核把HELP与HESIV落到可计算的形核过程上，使描述性的机理第一次有了通向断裂韧性预测的路径。计算给出的约束效应对工程评定有直接参考价值，约束越高，不仅韧性越低，进入脆性失效所需的氢含量也越低。与此同时，模型未能复现实验上出现的反转趋势，这一坦率的结论把问题清楚地留给了领域，也说明要建立可靠的氢环境断裂韧性评定方法，还需要针对氢控孔洞形核的定向实验与更自洽的理论描述。

▍原文链接

https://doi.org/10.1016/j.corsci.2026.114152