先筛选，再精算：剪切带预测的两步法

DOI：10.1016/j.ijplas.2024.104084推荐理由氢不总让镁更脆。Ji等在特定原子氢环境下测得纯Mg延伸率增加，同时强度下降。论文把表面入氢、晶内孪生和晶界开裂联系起来，说明判断氢的作用，必须看它在哪里、怎样参与变形。作者研究思路作者结合拉伸/EBSD、第一性原理（DFT）、原子模拟和晶体塑性（CPFEM）。贡献是用晶面能量变化约束含氢参数，并借助裂尖与多晶模型检验机制，再连接宏观响应。CPFEM依据EBSD建立587晶粒、500×500μm二维模型；空气参数由曲线标定，含氢参数结合DFT缩放。它并未显式求解氢扩散与裂纹演化。核心公式 式（7）：相对位移引起的能量变化除以界面面积，得到广义层错能；路径峰值用于比较变形势垒。 式（13—15）：采用已有晶体塑性框架，分别描述母相滑移、孪生转变与孪晶内再滑移。 式（16—18）：g为临界分解剪应力，孪生只取正向驱动；0.80是模型上限，不是实测孪晶含量。核心结果图 图 2(a)：10⁻⁶s⁻¹、pH14、−1.6V相对SCE，各条件3次。延伸率19.46%→22.06%，增加2.60个百分点，抗拉强度降低。 图 3(b)(c)(e)：依次为空气、含氢断后组织，白色标尺均为100μm；孪晶面积分数1.06%→2.84%，支持孪生参与增塑的解释。 图 6：柱面I入氢能垒0.156eV/H，低于所算基面及锥面I；不能仅凭能垒判断整体耐蚀性。 图 9(a)(b)：基面势垒降18.5%，稳定层错能为−7.13mJ/m²；柱面I势垒反增4.5%，作用有晶面选择性。 图 10(b)(c)：近0K裂尖模型中，含氢后裂纹偏折并发射位错。原σxy色条单位为eV/ų。 图 11(a)(b)(d)(e)：2H沿路径与6H环绕裂尖的设置同时改变数量、位置，后者层错发射反而延后。 图 12(a)(f)：近0K将氢保留在晶内，含氢模型至50%模拟应变未见孔洞；这不是实测延伸率。 图 13(a)(c)：10%应变下，含氢模型孪生更多、应力较低；曲线吻合仍需局部实验进一步检验。 图 14(b)(c)：全局孪生体积分数低于1.8%；增长速率差异主要集中在前4%应变，勿与局部场混同。 图 16(c)：Mg-H仍取约200ppm。300K下氢向晶界富集，24%模拟应变出现晶界孔洞，与近0K情景相反。作者核心结论晶内氢可降低部分变形势垒、促进孪生；晶界富集又可能诱发损伤。氢效应取决于位置、晶面和变形条件，延伸率提高不代表强度与耐蚀性同步改善。对我们的启示可借鉴能量约束本构、局部组织检验机制的思路。原子模拟约10⁸s⁻¹，与慢拉伸相差巨大；应用时还需可测氢分布与晶界证据，避免把特定情景直接当成工程加氢方案。来源：Ji等，International Journal of Plasticity 181（2024），104084。公式及结果图为原文节选、重排，版权归原权利人。来源：我的博士日记