引言
表面渐进淬火作为局部热处理核心技术,通过移动热源实现材料表层硬化与芯部韧性保留,广泛应用于齿轮、轴类等高承载零部件制造。然而,该工艺涉及强非稳态温度场、固态相变、热-相变-力学耦合等多重物理机制,传统经验公式难以定量预测淬硬层深度、相组成分布及残余应力状态。懿朵科技基于Abaqus平台,构建热-相变-力学耦合模型,为工艺优化提供精准仿真解决方案。
一、Abaqus多物理场耦合仿真核心技术
(1) 瞬态热传导与移动热源建模
建立三维热传导模型,集成材料内部热传导与表面对流换热效应
通过FORTRAN用户子程序实现高斯分布热源的连续移动,精准模拟热源轨迹与温度场动态演化
图1 表面渐进淬火热源作用系统示意及温度分布
支持复杂瞬态问题的稳定求解,避免传统方法简化误差
(2)相变动力学模型与用户子程序嵌入
在温度场计算基础上,研究进一步引入了固态相变动力学模型,用于描述钢材在加热与冷却阶段的组织演化行为。模型基于连续加热转变(CHT)与连续冷却转变(CCT)曲线,对以下相变过程进行了描述:
铁素体/珠光体→奥氏体(加热阶段)
奥氏体→铁素体、珠光体、贝氏体(扩散型相变)
奥氏体→马氏体(无扩散型相变)
相变模型通过Abaqus用户子程序嵌入求解流程,在单元尺度上实时更新各相体积分数,使得相变行为与温度演化保持严格一致。
图2 C60钢连续冷却转变(CCT)示意图
该方法避免了将相变结果作为后处理的简化做法,显著提升了仿真结果的物理一致性。
(3)热-相变-力学全耦合分析
相变过程不仅改变材料组织形态,同时伴随着体积变化,从而引入显著的相变应变。研究在Abaqus中进一步建立了结构分析模型,将以下应变分量统一纳入力学求解:
热应变
相变应变
弹塑性应变
在材料本构关系中,弹性模量、屈服应力等参数均设置为温度与相组成的函数,从而真实反映材料性能随热处理过程的动态演化。
这一建模策略充分发挥了Abaqus在非线性材料建模与多物理场耦合分析方面的工程能力,使得残余应力预测具备直接工程参考价值。
二、数值算例与温度场仿真结果
研究以一块尺寸为120×25×10mm的钢板为研究对象,利用几何对称性有效降低了计算规模。仿真结果显示,Abaqus能够清晰给出不同截面上的温度分布特征,并准确刻画热源移动过程中形成的高温区及热影响区范围。这些结果为后续相变区域及淬硬层形态分析提供了直接依据。
图 3 YZ 截面上的温度场分布
图 4 YX 截面上的温度场分布
在温度历史基础上,模型进一步预测了不同截面上的相组成分布情况。仿真结果表明,马氏体主要集中于靠近表面的区域,其分布形态与热源移动轨迹高度相关。通过对相变动力学过程的分析,可以定量评估淬硬层深度及不同相的体积分数变化规律。该结果表明,基于Abaqus的耦合建模方法能够有效预测热处理后的组织状态。
图 5 截面相组成分布及相变动力学过程
三、工程意义与Abaqus应用价值
研究结果表明:
(1)Abaqus可作为表面渐进淬火过程仿真的统一平台
(2)热–相变–力学耦合模型能够准确预测淬硬层形态
(3)用户子程序为复杂制造过程建模提供了高度灵活性
(4)数值仿真可显著降低试验成本,提高工艺设计效率
从工程应用角度看,该研究充分展示了Abaqus在热处理、材料相变及制造过程仿真领域的工程级能力,为复杂热加工工艺的数字化设计提供了可靠工具。
四、懿朵科技全流程仿真服务
(1)工艺建模服务
热源参数优化(功率、移动速度、轨迹设计)
材料相变动力学模型定制(CHT/CCT曲线嵌入)
多物理场耦合边界条件设置
(2)工程验证服务
温度场、相组成、残余应力实验对标
淬硬层深度无损检测(X射线衍射/硬度测试)
疲劳寿命加速试验验证
(3)数字化设计解决方案
工艺参数反演优化(基于遗传算法)
批量仿真自动化(Python脚本二次开发)
工艺知识库构建(参数-性能映射)
【参考文献】Tomasz Domański, Alžbeta Sapietová, Milan Sága,Application of Abaqus software for the modeling of surface progressive hardening,Procedia Engineering, 177 (2017) 64–69,DOI: 10.1016/j.proeng.2017.02.184.