在骨科植入物领域，可生物降解镁合金被誉为“革命性的一代材料”。它们能在人体内完成使命后自行吸收，避免了二次取出手术的痛苦与风险。然而，过快且不可控的腐蚀速率，一直是横亘在其临床广泛应用前的一道天堑。

为了调控降解速率，等离子体电解氧化（PEO）涂层技术应运而生。PEO涂层能在镁合金表面形成一层陶瓷质保护层，成为抵抗生理腐蚀的第一道防线。然而，一个更隐蔽、更具破坏性的“隐形杀手”正在潜伏——应力腐蚀开裂。

当植入物在体内承受日常活动带来的持续力学载荷，同时又处于复杂的生理腐蚀环境中时，裂纹会在远远低于材料屈服强度的应力下萌生、扩展，最终导致植入物的过早、无预警失效。这种现象并非单一的电化学腐蚀或力学断裂，而是两者耦合作用的结果，其背后更是涉及阳极溶解、氢致开裂等多种复杂机制的博弈。

图1. 所提出模型中多物理场相互作用的示意图。

一、挑战：应力腐蚀的多物理场耦合困境

要理解应力腐蚀开裂的复杂性，我们首先要看这场“战争”发生的战场。PEO涂层并非完美无瑕，其表面存在微孔和微裂纹等固有缺陷。

• 腐蚀前沿：电解质通过这些微缺陷渗透，接触镁合金基体。基体中的金属间化合物与镁基体之间形成微电偶，驱动局部腐蚀，形成腐蚀坑。
• 力学前沿：人体活动产生的载荷导致涂层与基体界面产生应力集中，诱发裂纹萌生。
• 氢脆危机：腐蚀过程中产生的氢原子扩散进入镁基体，降低晶界结合力，使材料变脆，加剧裂纹扩展。

这三个过程——腐蚀、断裂、氢扩散——相互交织、互为因果，形成了一个极其复杂的多物理场耦合问题。传统的实验手段难以实时、原位地观测这一动态过程，而现有的计算模型存在以下关键瓶颈：

1. 氢脆效应的缺失：先前Zhang等人[31]提出的断裂-腐蚀耦合相场模型虽已取得突破，但忽略了氢原子的扩散及其对断裂韧性的劣化作用，无法完整刻画应力腐蚀的全貌。
2. 微观结构与宏观性能的脱节：Gazenbiller等人[33]的宏观损伤模型无法捕捉涂层孔隙、界面粗糙度等微观特征对局部腐蚀与断裂的诱发作用。
3. 唯象模型的机理局限：连续损伤力学（CDM）模型虽能描述损伤演化，但缺乏对电化学-力学耦合机理的物理解释；而近场动力学（PD）模型虽能模拟裂纹萌生，但其计算成本高昂且边界条件处理复杂。

图2. 缺口板的计算结果。(a) 在加载时间 
 s 和 
 s 时，考虑氢脆与忽略氢脆两种情况下预测的腐蚀坑演变；(b) 在加载时间 
 s 和 
 s 时，考虑氢脆情况下预测的裂纹场、von Mises应力场和归一化氢浓度场的云图。

二、破局：一个热力学一致的相场框架

该研究的核心创新在于，提出了一个热力学一致的多物理场相场模型，首次将氢扩散效应、PEO涂层微观结构与应力腐蚀断裂统一在一个理论框架下。模型引入了五个关键的全局变量来刻画系统的演化：

1. 位移场：描述材料的变形。
2. 断裂相场
   ：用一个在0到1之间连续变化的变量来表征裂纹，0代表完好，1代表完全断裂。这避免了传统断裂力学中需要追踪复杂裂纹尖端的难题。
3. 腐蚀相场
   ：类似地，用一个变量表征材料从未腐蚀（1）到完全腐蚀（0）的连续转变。
4. 镁浓度场
   ：追踪镁离子的扩散，刻画腐蚀的化学进程。
5. 氢浓度场
   ：追踪氢原子的扩散与富集，定量描述氢脆效应。

通过严谨的热力学推导，模型清晰地刻画了各物理场之间的耦合关系：

• 变形驱动断裂：弹塑性变形能驱动断裂相场的演化。
• 裂纹加速腐蚀：裂纹的出现为离子和氢的扩散提供了“高速公路”，显著加快了腐蚀和氢扩散速率。
• 腐蚀与氢脆协同弱化：腐蚀产物和富集的氢原子会共同降低材料的断裂韧性，使得裂纹更易扩展。
• 断裂促进反应：裂纹尖端新鲜暴露的金属表面，会加剧局部电化学反应。

理论创新性的核心体现在于氢脆效应的引入方式。模型采用了氢致减聚（HEDE）机理，通过将临界断裂能
表达为氢浓度的函数：

其中。这一方程首次在相场框架下实现了氢原子对材料断裂韧性劣化的定量刻画，使得模型能够捕捉氢致裂纹扩展的动力学特征。

图3. (a) 不同应力水平下静态加载1天（240 MPa加载工况时为2183 s）的腐蚀坑形貌与裂纹场分布；(b) 加载240 MPa下失效试样的PEO涂层表面带有定向微裂纹的扫描电镜图像。

三、验证：从现象到机理的量化洞察

模型的强大不仅在于理论构建，更在于其与实验结果的惊人吻合以及对复杂机制的深入洞察。

1. 氢脆效应的定量揭示
研究首先对一个含初始缺陷的缺口板进行模拟。结果显示，考虑氢脆效应时，相同加载时间下的裂纹深度是不考虑时的两倍（图4）。环境中的氢浓度从0.1 wtppm提高到0.3 wtppm，裂纹萌生时间提前约300秒，最终失效时间加速约600秒。这首次在计算层面揭示了氢原子对应力腐蚀开裂的“催化”作用，并量化了氢损伤系数
对裂纹扩展速率的显著影响——当
从0.5增至0.9时，相同加载时间下的裂纹深度从64 μm增至166 μm。

2. 与经典模型的完美对标
研究将模拟结果与经典的“膜破裂/溶解/再钝化”模型进行对比。结果显示，模拟得到的应力腐蚀裂纹扩展速率与应力强度因子
之间遵循完美的幂律关系：

这一结果与经典的实验观察高度一致，不仅验证了模型的正确性，更表明该相场模型能够从底层机理上重现宏观的断裂力学行为。这一对比的意义在于：首次证明了基于热力学原理的相场模型能够自洽地“涌现”出唯象模型中需要人为假定的幂律关系，体现了模型在机理层面的完备性。

3. 真实PEO涂层的实验验证
研究进一步将模型应用于真实的医用镁合金WE43MEO及其PEO涂层。通过构建包含随机孔隙和初始微裂纹的微观结构模型，模拟结果与体外实验数据完美匹配（图11-12）。

更关键的是，研究揭示了不同载荷水平下的失效机制（图14-15）：

• 低载荷（120 MPa）：涂层虽然完全破裂，但基体损伤有限，腐蚀坑深度约30 μm。
• 中载荷（180 MPa）：涂层破裂加剧，腐蚀产物破坏明显，加速了坑蚀。
• 高载荷（240 MPa）：仅2183秒后，裂纹就深入基体达89 μm，预示材料即将完全失效。

值得注意的是，在高载荷下，涂层-基体界面发生大面积脱粘，消除了涂层对基体的残余保护作用。这一现象首次从计算角度证实了：在力学主导的失效模式下，PEO涂层的保护效果几乎可以忽略，基体本身的抗应力腐蚀性能成为决定性因素。

图4. 后封孔处理且无初始缺陷的PEO涂层样品的代表性计算模型与边界条件。坐标系原点位于模型的左下角。

四、创新性再审视：本研究对领域的三大贡献

1. 首次将氢脆效应纳入PEO涂层镁合金的应力腐蚀相场模型
与Zhang等人[31]的前期工作相比，本研究通过引入氢浓度场及其与断裂韧性的耦合关系，实现了对氢致裂纹扩展的定量刻画。这一突破使得模型能够捕捉氢原子在裂纹尖端的富集效应，以及其对裂纹扩展动力学的显著加速作用。

2. 建立了微观结构与宏观性能的定量关联
通过对涂层孔隙率、界面粗糙度等微观参数的敏感性分析，研究揭示了这些微结构特征与应力腐蚀速率之间的定量关系。例如，涂层孔隙率从0.05增至0.20，平均腐蚀深度增加78%；界面粗糙度从2 μm增至5 μm，腐蚀深度增加58%。这些结果首次为PEO涂层的微观结构优化提供了量化依据。

3. 实现了从理论到实验的完整闭环验证
研究不仅通过缺口板模型验证了氢脆效应的理论预测，还通过CT试样与经典膜破裂模型的对比验证了裂纹扩展动力学的普适性，最终通过与WE43MEO合金的体外实验数据对比，完成了从理论到实验的完整闭环。这种多层次的验证策略，确保了模型的可靠性与普适性。

图5. 加载速率为 时，后封孔涂层样品在加载时间 

、
和 
下的裂纹相场与腐蚀相场云图。

五、展望：后处理工艺的“数字孪生”

该模型的真正价值在于其预测能力。研究最后对后封孔处理后的PEO涂层进行了数值分析（图17-24）。通过模拟不同加载速率、不同环境氢浓度、不同涂层孔隙率和界面粗糙度下的应力腐蚀行为，得出了一系列具有工程指导意义的结论：

• 加载速率决定失效模式：更慢的加载速率（1×10⁻⁷ s⁻¹）意味着更长的腐蚀暴露时间，导致腐蚀损伤更严重，在更小的整体应变下即可引发失效。
• 环境氢浓度是“关键变量”：骨腔环境（0.16 wtppm）下的点蚀深度比皮下组织环境（0.06      wtppm）高出38.7%，为植入部位的选择和风险评估提供了量化依据。
• 微观结构决定宏观性能：涂层孔隙率从0.05增加到0.20，平均腐蚀深度从19.4 μm增加到34.5 μm。而更平滑的涂层-基体界面能有效缓解应力集中，提升涂层的保护效果