Abaqus焊接仿真：热力耦合分析的建模策略与收敛技巧

Abaqus焊接仿真是典型的移动热源瞬态热力耦合问题，技术门槛在”耦合”两个字上——热场和应力场的前后衔接、时间步长在两个物理场之间的平衡、高温区材料属性的剧烈非线性，任何一个环节掉链子，结果要么死在半路，要么出来一组物理上说不通的应力分布。

移动热源的建模是开局第一关。双椭球热源模型（Goldak）在Abaqus中通过DFLUX子程序实现，参数包括前后半椭球的长度比例、热源宽度、熔池深度以及热输入效率。这五个参数决定了焊缝区域从熔化到凝固的整个热历史——温度曲线的峰值、冷却速率、以及热影响区（HAZ）的宽度都由这五个参数间接决定。经验是：熔池宽度以实际焊缝横截面的金相观测值为参照，前后椭球长度比取1.5~2.0之间，热输入效率对电弧焊取0.7~0.85。参数确定之后必须在无应力条件下先跑一遍纯热传导分析，确认热循环的峰值温度和冷却速率在与实测数据可比的范围内，再联入应力分析。

如果实际焊接实验有热电偶测温数据，可以用这些数据来反演校准Goldak模型参数——固定熔池几何参数用金相结果，然后把热输入效率作为唯一的自由参数做调节，直到数值预测的热循环曲线与热电偶实测数据在峰值温度和T₈/₅冷却时间上对齐。这种单参数校准比盲目试多参数组合有效得多。

分析步的时间离散需要分阶段差异化处理。焊接加热阶段的时间步长取0.1~0.2秒——这个量级是为了捕捉移动热源的局部梯度。冷却阶段可以把时间步长逐渐拉开，因为温度梯度随时间递减，不需要维持同样的瞬态精度。总的模拟物理时间必须覆盖到焊件冷却至室温附近——残余应力在冷却至200℃以下之后基本定型，再往后变化的量级不大。

材料属性随温度变化的非线性需要全套输入。从室温到熔点以上两百度的杨氏模量、屈服强度、热膨胀系数、热导率、比热容，缺一不可。部分高温段的参数——尤其是进入塑性流变区后的应力-应变数据——文献里往往不完整，需要通过JMatPro等材料性能计算软件做补充或外推。缺失了高温塑性参数，焊后残余拉应力的预测值会严重偏低。

收敛困难最常出现在焊缝附近的单元进入塑性软化阶段时。标准Newton-Raphson在这一区域的切线刚度矩阵接近奇异，迭代失败率很高。有效的应对措施包括：打开自动稳定（stabilization）并控制耗散能占比在5%以内；在建模时对焊缝附近的网格做加密但单元质量必须保证，长宽比不建议超过3:1；如果仍然不收敛，回到材料定义里检查高温段的强化系数是否为零或接近零，必要时在应力-应变曲线的尾部给一个微量的非零斜率防止矩阵零奇异。
焊后残余应力的验证是焊接仿真的一条底线。数值计算给出的焊缝纵向残余拉应力通常在屈服强度附近——这和理论预期一致，因为焊缝区在冷却收缩过程中被周围的冷金属约束，必然进入屈服。但如果数值结果显著低于屈服强度（比如低了30%以上），那大概率是高温塑性参数不全或者单元在塑性软化区数值失稳导致的。一个简单的验证方法：在仿真结果中沿焊缝中心线取一条路径，提取纵向应力分布，和文献中同类材料、同类焊接工艺的残余应力分布做定性对比。峰值位置、应力符号转换点、HAZ宽度——这三个特征对上了，才能说仿真结果是”物理上合理的”。

https://www.keyanxueshu.com/category/fea/