单晶钒中的动态应变时效：材料表征与基于位错的单晶塑性框架

《METALLURGICAL AND MATERIALS TRANSACTIONS A-PHYSICAL METALLURGY AND MATERIALS SCIENCE》：Dynamic Strain Aging in Single Crystal Vanadium: Material Characterization and Dislocation-Based Crystal Plasticity Framework

　　

研究人员针对单晶钒（V）的力学行为开展研究，并引入一种改进的基于位错的晶体塑性框架以模拟其在不同加载条件下的响应。研究人员在宽温度与应变率范围内对单晶钒样品进行实验，揭示了低应变率与较高温度下渐进的热硬化现象，该行为指示了动态应变时效（DSA），此前该现象已在多晶钒中被观测到。受DSA影响，常规单晶塑性模型——对其他体心立方（BCC）金属（如钽或钼）有效——无法准确捕捉DSA区间内钒的力学响应。为克服该限制，研究人员改进了饱和位错密度公式以纳入DSA效应，同时保持对其他BCC晶体的预测能力。该方法相对简洁且稳健，对材料在DSA区内外的强度行为均给出满意预测，成功预测了已报道的多晶钒流变应力温度依赖性，并预测了率跳跃实验中与负应变率敏感性相关的应力调整。

研究背景与立项依据

钒（V）作为强金属键结合的过渡金属，因高熔点与高温性能优异，被视为聚变堆、航天器及舰船等极端环境候选结构材料。然而相较于钽、钼等同为体心立方（BCC）结构的过渡金属，钒的强度行为研究仍不充分：已有剪切强度数据多来自多晶钒准静态至10⁴ s^-1的单轴实验，单晶钒数据稀疏且限于准静态；低温以变形孪晶为主，近室温以上转为滑移主导，且对间隙杂质敏感、微纳尺度存在尺寸效应。多晶钒在500–900 K、低应变率（≤10^-3 s^-1）下已确认存在动态应变时效（DSA），表现为热硬化与负应变率敏感性，源于溶质原子（O、N等）扩散钉扎运动位错导致硬化模量上升、总位错密度增加；但单晶钒此前未明确报道DSA，且常规BCC单晶塑性模型（不含DSA机制）在DSA温区会错误预测“高温饱和应力低于低温”。因此，研究人员旨在（i）通过单轴压缩实验表征单晶钒宽温域、宽应变率剪切强度，（ii）在既有位错基晶体塑性模型上嵌入简洁DSA项，实现DSA区内外统一预测。

主要关键技术方法

研究人员购自Princeton Scientific的>99.99%纯单晶钒，取向[100]与[111]，Laue背反射校验取向，加工为?5×5 mm圆柱。准静态压缩用MTS 880伺服液压机（193–673 K，10^-3或10^-1 s^-1，LVDT与引伸计测应变）；高应变率用分离式霍普金森压杆（SHPB，293–673 K，1.8–2.5×10³ s^-1）。模型端基于Nguyen等已发表的位错基BCC单晶塑性框架（乘法分解变形梯度、Green-Lagrange弹性、Mie-Grüneisen状态方程、位错密度分移动/不动ρ_m/ρ_im、Orowan关系、热激活等待时间动力学），将饱和位错密度ρ^sat由原MTS型公式扩展为“基准项+双Weibull型DSA项”；参数用pyBASS构建应力-应变主成分代理模型、Impala做贝叶斯校准（10链、5万迭代、1万燃烧），以单晶实验应力曲线拟合ρ₀^DSA与形状参数ω，γ_DSA=10^-2 s^-1与T_DSA=700 K借多晶钒先验固定。

研究结果

2 材料表征与测试

研究人员对[100]与[111]单晶钒做193–673 K准静态及SHPB压缩，发现总体符合BCC趋势：193→473 K热软化、10^-3→10³ s^-1应变率硬化、[111]强于[100]；但在473→673 K、10^-3 s^-1出现明确热硬化（DSA信号），且[100]取向在应变<10%时673 K应力即超473 K，[111]取向则迟至>25%应变才超越，显示DSA具取向依赖性。与文献微柱压缩（~3 μm）比，毫米级样品初始屈服高（低初始位错密度致黏性应力大），但室温流动应力与Han等3420 nm微柱吻合；193 K下流动应力仅为微柱一半，说明低温尺寸效应临界直径D_c(193 K)?3 μm。因供应商未提供批特异杂质浓度，DSA中间隙溶质贡献未能定量分离，但热硬化跨两取向一致，构成单晶钒DSA指示性证据。

3 考虑DSA的动力学位错基晶体塑性模型

3.1 理论模型构建

基准模型饱和位错密度ρ_l^sat,α=ρ₀^sat(γ^α/γ₀)^k_BT/A随T升高而降，与DSA矛盾。研究人员改写为ρ_l^sat,α=ρ₀^sat(γ^α/γ₀)^k_BT/A+ρ₀^DSAexp(-γ^α/γ_DSA)(T/T_DSA)^ω-1exp(-(T/T_DSA)^ω)。新增项为应变率Weibull（γ^α?γ_DSA时消失）与温度Weibull（ω>1时峰值在T_DSA(1-1/ω)^1/ω）之积，直接抬升饱和位错密度而非另设时效时间变量，从而避免求解溶质浓度演化刚性问题，也不改变总位错密度物理含义。γ_DSA=10^-2 s^-1（单晶与多晶DSA消失边界一致），T_DSA=700 K取多晶峰中位，ω与ρ₀^DSA由单晶数据校准。

3.2 校准模型响应

模型在LANL的FLAG（任意拉格朗日-欧拉多物理码）实现。校准后整体复现实验：非DSA区随T升高流变应力降；DSA区（673 K vs 473 K，10^-3 s^-1）热硬化被捕获，[111]预测优于[100]（后者略高估），SHPB区正应变率敏感与热软化正确，但未捕捉[100]/[111]的DSA onset应变差异——因仅两取向、为避免过拟合未给非等价滑移系分配独立DSA参数，模型给出滑移系平均DSA响应。早期屈服尖峰（低初位错密度特征）未被拟合，研究人员将其归为样品预处理状态变量而非材料常数，故只校大应变段。

讨论与结论翻译

研究人员得出结论：单晶钒在473–673 K、10^-3 s^-1呈现与多晶一致的DSA指示性热硬化，且具取向依赖表现；通过将DSA效应嵌入饱和位错密度（双Weibull型温度-应变率函数，不引入额外内变量），改进后的位错基BCC单晶塑性框架在DSA区内可复现热硬化与负应变率敏感，区外保持对Ta、Mo类BCC金属适用形式，并成功外推预测多晶钒流变应力温度曲线及率跳跃应力调整。该框架计算开销低、便于大尺度多晶模拟，为钒基合金在聚变堆等高温结构的力学行为预测提供了可校准、物理一致的建模基础；后续若获更多取向与宽溶质浓度单晶数据，可进一步分辨滑移系层面DSA异性或耦合溶质浓度项。