多相高强钢拉伸应变分配的晶体塑性有限元模拟与准原位实验研究

《Coatings》：Crystal Plasticity Finite Element Simulation and Quasi-In-Situ Experimental Study of Tensile Strain Partitioning in Multiphase High-Strength Steel

　　

研究人员采用准原位拉伸表征和基于电子背散射衍射（EBSD）的晶体塑性有限元（CPFE）方法，研究了在260 °C下等温淬火24 h的多相高强钢。实验观察揭示局部塑性变形具有显著异质性：冯·米塞斯（von Mises）应变优先集中在贝氏体铁素体束、相界以及残余奥氏体/马氏体-奥氏体（M/A）区域附近，而块状残余奥氏体在早期变形阶段有助于应变协调。为量化潜在的应力-应变分配，研究人员从EBSD数据重建了准二维代表性体积单元（RVE），并通过用户材料子程序（UMAT）在ABAQUS中实现。模型包含了24 h等温淬火试样的真实晶粒形貌、相分布和晶体取向信息。研究人员校准了具有体心立方（BCC）基体、面心立方（FCC）残余奥氏体和转变马氏体分支的率相关晶体塑性本构框架，以匹配宏观拉伸曲线。模拟的拉伸响应在宏观失稳前与实验曲线吻合良好，预测的局部场与准原位应变图一致。结果表明，局部塑性应变首先在M/A相关区域和相界相邻区域累积，而高米塞斯（Mises）应力随滑移活动和应力诱发马氏体相变动态迁移。残余奥氏体转变提高了局部承载能力，改变了相间载荷传递，并延迟了应变局域化带的直接连通。本工作阐明了多相高强钢中残余奥氏体稳定性、相变诱导塑性（TRIP）辅助载荷再分配和微观结构应变分配之间的耦合关系，为微观结构指导的强度-延性优化提供了细观尺度基础。

**研究背景与问题**
多相高强钢因其通过贝氏体铁素体、残余奥氏体、马氏体及马氏体-奥氏体（M/A）组元的协同变形实现高强度与良好延展性及韧性的结合，在承载构件领域备受关注。残余奥氏体是相变诱导塑性（TRIP）辅助和无碳化物贝氏体钢中的关键组元，其机械稳定性受尺寸、形貌、碳富集、周围约束及晶体取向控制，决定TRIP效应是逐步发生还是过早发生。然而，宏观拉伸曲线无法揭示局部塑性应变起始位置、不同相之间的载荷分担及损伤倾向于在特定相界区域萌生的机制。准原位变形实验虽能观察到高度异质的局部应变，但难以获得全场应力状态与相特异性塑性的时变再分配。晶体塑性有限元模拟（CPFEM）为连接晶体学滑移、相拓扑与局部力学响应提供了有效细观框架，但针对具有真实EBSD相形貌的贝氏体高强钢，仍需将局部场预测与准原位应变测量进行定量验证。为此，本研究旨在结合准原位表征与EBSD基CPFEM，阐明260 °C×24 h等温淬火高强钢中应变局域化起始及TRIP辅助载荷再分配的演化。研究论文发表在《Coatings》。

**主要技术方法**
研究人员使用0.44 V高强钢（化学成分见表1），经900 °C奥氏体化30 min后转入260 °C盐浴等温淬火24 h，空冷至室温。关键技术包括：①准原位拉伸表征，在配备拉伸台的扫描电镜（SEM）下以1 mm/min速率进行，采集不同工程应变下的SEM图像、相图、局部取向差（KAM）图及应变场信息；②EBSD基RVE构建，直接从EBSD数据中提取真实晶粒形貌、相分布和晶体取向，构建准二维有限元模型（约10.496 μm×10.166 μm×0.0939 μm，含168个晶粒/相区、35936个C3D6楔形单元和36411个节点），右面耦合参考点施加位移载荷至1.784 μm（对应约0.17工程应变）；③基于有限变形运动学的率相关晶体塑性本构框架（在ABAQUS/UMAT中实现），BCC基体激活{110}<111>滑移系，FCC残余奥氏体激活{111}<110>滑移系，转变马氏体分支M1、M2使用{110}<111>滑移系，并引入现象学应力辅助相变分支，相变驱动力D_i = |τ^ir_i| ? λ(T ? θ_T)/θ_T，参数λ=546，T=293 K，θ_T=246 K，F_c=26，V_max=4.0×10^?5，V_e=0.17；④参数标定采用物理约束的试错法，调整τ₀、h₀、τ_s以匹配宏观屈服、早期硬化斜率及高应变流动应力。

**研究结果**
**3.1 准原位变形与应变分配**：准原位观察显示，随工程应变从0.01增至0.11，局部应变优先集中在贝氏体铁素体束交汇处和相界相邻区域，块状残余奥氏体区在早期应变较低。EBSD相图统计表明残余奥氏体面积分数由未变形时的14.7%降至6%应变时的2.7%和12%应变时的1.2%，证实残余奥氏体逐步消耗，部分残余奥氏体在12%应变时仍保持稳定。

**3.2 EBSD基RVE与模型验证**：重建的RVE保留了实验微观结构的主要特征，包括板条状和块状相形貌、晶体取向差异及残余奥氏体空间分布。CPFE预测的拉伸曲线与实验曲线在0.2~16%应变范围内对比，均方根误差（RMSE）28.21 MPa，平均绝对百分比误差（MAPE）1.51%，R²=0.9906。在最大载荷及高均匀变形区域模拟曲线略低于实验，归因于τ_s、h₀的保守标定及转变诱发硬化处理简化。

**3.3 局部米塞斯应力演化**：米塞斯应力场高度异质，高应力区首先出现在BCC基体部分和相界附近，随应变增加向残余奥氏体相关区域扩展并形成互联路径。高应力位置动态迁移，当残余奥氏体转变为马氏体时承载能力提高，高应力区转移至邻近未充分变形区域，维持持续加工硬化。

**3.4 等效塑性应变与TRIP辅助载荷再分配**：PEEQ（等效塑性应变）沿路径呈现多个快速放大的峰值，对应相界相邻及几何约束区域。M/A相关区域在早期塑性应变高于BCC基体。马氏体体积分数（SDV2376）随应变显著升高，相变后局部承载能力增强，使M/A区平均PEEQ增长速率变化，而整体平均PEEQ稳步增长。模拟的PEEQ峰值沿局域化路径集中，高应力区沿晶界和相界迁移，与准原位von Mises应变图一致。

**3.5 细观变形机制**：24 h等温淬火高强钢的拉伸变形分为三阶段：早期，M/A相关及相界相邻区域因稳定性低和约束强优先启动局部塑性；中期，残余奥氏体转变为马氏体增强局部承载，应力从初始变形区向相邻晶粒再分配；后期，高应力和高PEEQ区逐渐互联，为损伤萌生和裂纹扩展提供潜在路径。残余奥氏体并非单纯软相，马氏体也非理想刚性相，局部响应受晶体滑移、相变诱发硬化和相间约束的动态竞争控制。

**讨论与结论**
讨论部分指出，模拟与实验在局部场特征上的一致性表明CPFE模型不仅拟合宏观拉伸曲线，还捕捉到局部应变集中和TRIP诱导载荷再分配的关键细观机制。残余奥氏体稳定性、形貌及周围约束共同决定其转变行为，部分保留奥氏体在12%应变时仍稳定支持了稳定性的非均匀性。结论如下：(1) 建立了260 °C×24 h等温淬火高强钢的EBSD基准二维RVE模型，保留真实晶粒形貌、相分布和晶体取向，含168个晶粒/相区、35936个C3D6单元和36411个节点。(2) CPFE模型再现了实验拉伸曲线的主要特征，在0.2–16%工程应变范围内RMSE为28.21 MPa，MAPE为1.51%，R²=0.9906，表明标定的相依赖参数适用于细观应力-应变分配分析。(3) 单轴拉伸过程中局部米塞斯应力和PEEQ高度异质，高应力区优先在晶界、相界及取向差显著区域形成，PEEQ峰值沿与相界约束和转变诱发不相容性相关的局域化路径发展。(4) 残余奥氏体在拉伸过程中逐步转变为马氏体，EBSD相分数统计显示残余奥氏体面积分数从未变形的14.7%降至6%应变时的2.7%和12%应变时的1.2%，证实其消耗，而12%应变时的残留表明部分残余奥氏体保持机械稳定。(5) 预测的局部场演化与准原位应变观测一致，证明当前EBSD基CPFEM框架能同时捕捉多相高强钢的宏观拉伸行为和细观应变分配机制。