各向异性是关键：304L不锈钢中氢如何“定向削弱”原子键合力？

1 引言

氢脆是制约结构合金在储氢、海工与焊接接头等场景中长期安全服役的关键因素之一。当前的机理研究主要聚焦于塑性变形阶段，代表性框架包括氢增强局部塑性(HELP)和氢增强脱粘(HEDE)模型。然而，氢作为间隙原子对材料弹性阶段——尤其是单晶弹性常数(Cij)——的定量影响至今仍不清晰。由于晶格层次的弹性行为直接决定了后续位错形核、应力场分布和微观结构稳定性，因此精确量化氢致单晶弹性常数变化是揭示氢脆早期驱动力的前提。现有方法在直接实验测量方面存在显著局限：超声波测试无法与原位拉伸同步，纳米压痕仅反映极表面行为，而第一性原理(DFT)模拟则受限于0 K与无热近似，且忽略了多晶集合体的晶间弹性耦合。

针对上述空白，釜山大学文彬若、成白石、崔东贤等人在Materials期刊（2026年，第19卷）发表了题为“Hydrogen-Induced Anisotropy in Single-Crystal Elastic Constants of 304L Stainless Steel via In Situ Neutron Diffraction and Kröner Modeling"的论文。本研究以304L奥氏体不锈钢为对象，结合原位中子衍射与Kröner自洽模型，从实验上反演得到充氢前后的单晶弹性常数，系统定量评估氢对体积效应与模量效应的影响，并将其与HELP机制的微观起源相关联。

2 材料与研究方法

2.1 材料与微观组织表征

实验采用20 mm厚的304L奥氏体不锈钢热轧板，经1060 ℃/1 h固溶处理后水淬。化学成分(质量分数，%)为：C 0.02、Cr 18.1、Ni 8.0、Mn 1.6、Mo 0.1、Si 0.4，余Fe。固溶态试样屈服强度(0.2%偏移)为341 MPa，充氢与未充氢试样之间屈服强度差异不显著。试样经自动抛光(胶态二氧化硅悬浊液，6 N，70 r/min)后，使用配备EBSD探测器的场发射SEM进行微观组织观测，EBSD扫描参数为15 kV、17 mm工作距离、1.0 μm步长。

图1 拉伸试样尺寸规格(单位：mm，R6表示6 mm圆弧半径)

2.2 气相预充氢与热脱附分析

拉伸试样按ASTM E8标准加工，采用800号碳化硅砂纸精磨表面。充氢工艺为气相氢环境，温度280 ℃、压力16 MPa、保温168 h。为排除热历史差异，未充氢对照试样进行同等条件(280 ℃/168 h)的等温热暴露。根据扩散系数D=7.2×10⁻⁸ cm²/s，1.5 mm厚试样从两侧同时扩散，所需最大扩散距离仅0.75 mm，可保证氢在试样截面内均匀分布。充氢后采用气相色谱TDS定量氢含量，以200 ℃/h恒速率升温，通过不对称高斯(SplitGaussian)对TDS谱进行分峰拟合，识别不同氢陷阱位点。

2.3 原位中子衍射与中断拉伸

原位中子衍射实验在韩国原子能研究院(KAERI)残余应力衍射仪上完成，集成最大载荷20 kN的拉伸台。入射狭缝5 mm×10 mm，探测器狭缝5 mm宽，采用弯曲Si(220)单色器获得波长λ=1.46 Å的单色中子。散射矢量Q平行于加载方向LD，从而实现轴向应变的同时测量。所有测量在固定试样方位(恒定ψ角)下完成，拉伸过程中不进行试样倾斜。

拉伸应变速率恒定为1×10⁻⁴ s⁻¹。为保证测量严格处于弹性区，每次衍射采集前中断加载，采集应力水平均低于300 MPa，每条衍射峰计数时间为135 s。数据采集后恢复加载，采用非接触式激光引伸计实时监测试样标距长度。

2.4 晶格参数与晶格应变的反演

衍射峰位与峰形采用Fityk软件VoigtA算法拟合，扣除背景后在GSAS-II中以8项Chebyshev-1多项式拟合背景。304L试样的晶格参数通过Cohen法线性化Bragg定律后三维最小二乘回归求得;晶格应变ε由实测晶面间距d与无应力参考间距d₀之差归一化获得。与反射的峰位根据Warren理论按层错概率加权修正，以避免堆垛层错引入的表观晶格膨胀干扰。

2.5 Kröner自洽模型与Cij反演

Voigt与Reuss近似分别给出多晶弹性模量的上界(等应变)与下界(等应力)，但两者均无法同时满足晶界处的应力平衡与应变连续。Kröner自洽模型将每颗晶粒视为嵌入均匀等效基体中的椭球夹杂，通过取向因子Γ=(h²k²+k²l²+l²h²)/(h²+k²+l²)²与单晶弹性参数μ=(C11−C12)/2、η=C44、K=(C11+2C12)/3建立三次方程，逐晶面族反算衍射剪切模量G、杨氏模量E与泊松比ν。最终通过加权最小二乘(χ²最小)迭代调整C11、C12与C44，使模型预测的1/EK与实验1/EEX达成匹配。宏观各向同性下Γ=1/5即可恢复多晶剪切模量GKM，进而由标准关系得到EKM与νKM。

3 氢陷阱分布、微观组织与晶格膨胀

图2 预充氢试样的TDS谱及其分峰结果(对应晶界与晶内两类陷阱位点)

TDS分析得到试样总氢含量为14.9±0.1 mass ppm(两次重复测量)，分峰后识别出晶界陷阱与晶内陷阱两类位点，对应氢含量分别为9.4±0.5 mass ppm和5.5±0.5 mass ppm。低温峰(~400 ℃)对应晶界捕集氢，高温峰(~550 ℃)对应晶内间隙氢，与既有304L奥氏体不锈钢气相充氢研究的TDS谱一致，印证了氢在不同微结构位置的差异化分布。

图3 (a，b)304L试样初始组织的EBSD图，分别为反极图与相图;(c，d)非充氢与充氢试样中断拉伸过程中的原位中子衍射谱

EBSD统计表明，304L试样平均晶粒尺寸为41±19 μm(面积加权平均)，组织以γ奥氏体为主(99.9±0.1%)，仅含可忽略的α'马氏体相(0.1±0.1%)。由于充氢温度(280 ℃)远低于固溶温度(1060 ℃)，充氢前后晶粒尺寸与相比例基本保持不变。原位中子衍射谱(c、d)在弹性阶段(≤300 MPa)未出现应变诱发α'马氏体对应的新峰，说明充氢与拉伸加载均未引起显著的宏观相变，排除了相变对弹性常数测量的干扰。

通过Cohen法三维最小二乘回归，得到未充氢与充氢试样的晶格参数分别为3.558±0.008 Å与3.583±0.006 Å。间隙氢的引入使FCC晶格产生约0.7%(0.025 Å)的膨胀，该结果与图2中5.5±0.5 mass ppm的晶内间隙氢含量一致，表明晶格膨胀是间隙氢嵌入的直接结果。该膨胀为后续衍射弹性常数(DEC)的反演提供了关键基准。

4 单晶弹性常数与各向异性分析

4.1 衍射弹性常数的晶面族依赖性

图4 非充氢与充氢304L试样不同晶面族晶格应变随真应力的演化

图4展示了弹性区内不同晶面族晶格应变随真应力的演化，所有线性拟合的R²均高于0.97，表明应力-应变响应高度线性。各晶面族拟合斜率由小到大依次为：≈ < < < ，与FCC晶体弹性各向异性的取向依赖相符。充氢前后，、、的差异较小，但与面族的DEC出现明显下降，降幅分别约为3.2%与6.7%。这一结果表明氢对弹性常数的影响具有显著的取向依赖性，提示充氢进一步放大了晶体的弹性各向异性。

4.2 单晶弹性常数Cij的反演结果

将Kröner模型应用于未充氢与充氢两组试样的实测DEC数据，通过加权最小二乘反演得到单晶弹性常数。未充氢试样：C11=216 GPa、C12=126 GPa、C44=173 GPa;充氢试样(0.1 at.% H)：C11=243 GPa、C12=161 GPa、C44=171 GPa。充氢后，C11升高约13%，C12升高约28%，而C44仅下降约1%(几乎不变)。充氢与非充氢拟合的约化χ²分别为1.42与3.39，均处于统计可接受范围，表明反演结果稳健。

Zener各向异性比A=2C44/(C11−C12)由未充氢的3.8升高至充氢后的4.2，显著增大的弹性各向异性将通过三种协同机制放大氢脆敏感性：其一，加剧相邻晶粒间的弹性失配，促使晶界局部应力集中并萌生沿晶裂纹;其二，放大不同晶面族刚度差异，使后续塑性变形集中于主滑移面，形成强烈的应变局部化;其三，Zener比升高会按γ∝A⁻⁰·³⁷关系降低层错能，从而促进应变诱发马氏体的形成，进一步损害材料塑性。

4.3 与既有DFT/IET结果的比较

将本研究结果与已有DFT计算及IET宏观模量实验对比：Teus等针对FCC Fe的DFT结果显示C44随氢含量增加呈下降趋势;Shi等对BCC Fe的DFT研究表明氢使C11、C12、C44及Zener比单调降低;Hachet等对FCC Ni的DFT计算亦显示氢使各弹性常数轻微下降;Scott通过IET对316L和XM-19的宏观测量则与本工作一致，氢使杨氏模量轻微升高。综合对比可发现，DFT模拟大多预测氢使弹性常数降低，而本实验及Scott的宏观测量均显示氢使杨氏模量轻微升高、C11/C12上升。这一系统性差异主要源于DFT在0 K与无热近似下忽略了温度振动及氢原子的动态重排，且DFT通常针对孤立无约束单晶，无法反映多晶集合体的晶间约束。本研究基于Kröner模型反演得到的Cij内禀地考虑了多晶粒间弹性相互作用，更贴近真实工程材料。

5 多晶弹性模量与HELP机制的晶格起源

5.1 宏观模量与泊松比的差异化响应

基于Kröner模型求得的多晶弹性模量显示：充氢前后杨氏模量EKM由251 GPa升至252 GPa，几乎不变;剪切模量GKM由102 GPa轻微下降至99 GPa;而泊松比νKM则由0.23显著升高至0.28(升高约22%)。这一“宏观不变、泊松比显著升高"的表象本质上是体积模量K(由156 GPa升至188 GPa，升高20.5%)与几乎不变的剪切模量之间竞争的结果，K/G比值的显著增大反映材料抵抗体积变化的能力相对强于抵抗形状变化的能力。

宏观弹性模量的近似恒定解释了为何常规SSRT应力-应变曲线在断裂前几乎观察不到软化，但泊松比的显著变化提示多晶集合体内部的应力-应变分配在塑性流动启动前已经悄然改变。这种体积效应主导的弹性响应是后续应变局部化的先兆，也是HELP机制得以发挥作用的弹性基础。

5.2 剪切模量与Peierls-Nabarro应力

由单晶Cij计算得到<110>主滑移系的剪切模量G：非充氢为60 GPa，充氢后降至55 GPa，降幅约8.3%。基于此并结合Peierls-Nabarro应力公式τP−N=2G(1−ν)exp(−4πζ/b)估算，得到τP−N从未充氢的0.198 GPa降至充氢后的0.123 GPa，降幅高达约38%。尽管该经典公式仅针对位错并可能高估绝对值，但因其与实际面层错介导部分位错在同一滑移面上运动，G下降必使部分位错的P-N应力同向减小，从而实验证实HELP所必需的位错易滑移特征在弹性阶段就已埋下伏笔。

G的下降同时降低了层错能，抑制交滑移并促进平面滑移主导，塑性应变易在面局部化;P-N应力的下降则进一步降低位错在面上滑移的临界分切应力，使位错在应变集中区域更易启动并堆积。这些效应与Birnbaum与Sofronis所提出的HELP理论体积效应(volumetric effect)和模量效应(modulus effect)在晶格尺度上得到实验定量验证。

5.3 HELP机制的晶格尺度验证

结合上述结果与多种氢脆机理的对比，本研究支持以下论断：氢致C11上升从晶格层面否定了晶内HEDE(否则纵向刚度应下降);晶间HEDE与AIDE仅在表面或晶界局部起作用，与中子衍射测得的体弹性响应无关;而G与P-N应力的显著下降则直接佐证HESIV等依赖位错运动的机理。本研究是从单晶弹性常数层面实验证实HELP机制的体积效应与模量效应，为通过合金设计抑制HELP、改善奥氏体钢氢脆抗力提供了量化依据。

6 结论

本研究结合原位中子衍射与Kröner自洽模型，实验定量获得了304L奥氏体不锈钢充氢前后的单晶弹性常数及其对HELP机制的晶格起源意义，主要结论如下：

(1)气相充氢(280 ℃/16 MPa/168 h)后，304L试样的FCC晶格常数由3.558 Å膨胀至3.583 Å，膨胀量约0.7%，总氢含量14.9 mass ppm，其中晶界陷阱约9.4 mass ppm、晶内间隙约5.5 mass ppm。

(2)反演得到的单晶弹性常数表明：充氢使C11升高约13%、C12升高约28%，C44基本不变，Zener各向异性比由3.8升至4.2，体模量K升高约20.5%，而多晶杨氏模量与剪切模量几乎不变，泊松比由0.23显著升至0.28，呈现“体积刚化、剪切不变"的差异化响应。

(3)基于Cij进一步求得<110>滑移系剪切模量G由60 GPa降至55 GPa(−8.3%)，Peierls-Nabarro应力由0.198 GPa降至0.123 GPa(−38%)，从晶格尺度实验证实HELP机制中Birnbaum-Sofronis所提出的体积效应与模量效应。

(4)上述弹性阶段的变化早于宏观应力-应变曲线的任何软化迹象，可视为氢脆塑性局部化的早期前兆，为通过合金成分设计抑制HELP、改善奥氏体不锈钢氢脆抗力提供了关键定量依据。