焊接残余应力与变形的热-力耦合模拟

316L不锈钢板（200×150×6mm³）的TIG焊接，单道焊，电流120A、电压12V、焊接速度2mm/s。焊后在距焊缝中心30mm处用盲孔法（ASTM E837）测了表面残余应力：纵向约+280 MPa（拉伸），横向约+95 MPa。板的最大挠曲变形约2.7mm（中心拱起）。

用ABAQUS跑了一轮顺序热-力耦合分析：先用Goldak双椭球热源模型跑热传导（非线性，考虑对流+辐射边界），再把温度历史作为载荷导入应力分析。热分析和应力分析解耦（顺序求解），不考虑力学变形对温度场的影响（位移-温度的单向耦合）。

第一轮结果：纵向残余应力+355 MPa（偏高27%），横向+140 MPa（偏高47%），变形3.8mm（偏高41%）。三个量都在高估，幅度还不一样——变形偏差最大。接下来逐环节排查偏差的来源。

测试一：热源模型参数——Goldak参数对残余应力的敏感性比预期低

Goldak双椭球模型的四个几何参数（前椭球半轴a_f、后椭球半轴a_r、宽度b、深度c）决定了熔池的热输入分布。标准推荐值（a_f=3mm, a_r=8mm, b=4mm, c=4mm）在这个板厚（6mm）上偏保守——c=4mm意味着热源深度只覆盖了板厚的2/3，底部1/3靠热传导。

把c调到6mm（全板厚覆盖），熔池底部温度从约1450°C涨到约1580°C，但残余应力的变化幅度只有3-5%。纵向从+355变成+343 MPa，横向从+140变成+135 MPa。热源参数的调整对温度场有明显效果，但对最终的残余应力场影响极其有限——原因在于残余应力是由冷却阶段的温度梯度和塑性应变累积决定的，不是由熔池的峰值温度决定的。熔池大了或小了，冷却阶段的梯度分布差别不大，应力场也不会跟着大动。

测试二：相变潜热——有必要，但不决定结果

316L在凝固温度区间（约1370-1400°C）和固态相变温度区间（奥氏体→铁素体不稳定区，P92耐热钢更明显）有潜热效应。忽略相变潜热会让熔池附近的冷却速率偏高约15%——因为潜热释放会减缓熔池的冷却速度，让高温停留时间延长。

加上相变潜热（凝固潜热约270 kJ/kg，从JMatPro的热力学数据库提取），残余应力的变化幅度约6-8%。纵向从+355→+332 MPa。效果比热源参数大，但仍然不是主要误差源。相变潜热对焊缝金属在高温段的应力松弛有帮助（塑性应变累积时间更长），但这个效应在316L这种单相奥氏体钢上不如在P92这种有固态相变的钢上明显。

测试三：高温力学性能——屈服强度的温度曲线是个黑箱

316L的屈服强度-温度曲线在高温段（>800°C）很难测——高温拉伸试验的氧化、蠕变和动态回复会干扰测量，不同文献的数据分散性可达±20%。原本用的是一组文献数据（800°C时σ_y≈120 MPa），换成另一组更保守的数据（800°C时σ_y≈85 MPa），纵向残余应力从+355降到+310 MPa——下降了13%，是之前几项调整中幅度最大的一次。

物理逻辑也很清晰：残余应力来自热膨胀/收缩被约束后产生的塑性变形，而塑性变形的门槛就是温度相关的屈服强度。高温段屈服越低，冷却时越早进入屈服，塑性应变累积越多，最终的残余应力越低——因为塑性卸载把弹性应力释放了。所以高温屈服强度每降低10%，残余应力大约降低5-8%。而这个参数的真实值在不同文献之间差了很多，是一个被严重低估的误差源。

测试四：拘束条件——这才是大头

前面的几项调整加起来让残余应力从+355回到+310 MPa，离实验值+280还有30 MPa的差距。变形从3.8mm降到3.2mm，离实验值2.7mm还有0.5mm。这30MPa和0.5mm藏在哪？

出在拘束条件上。初始模型用了简化的完全刚性约束（焊板底面固定），但在实际焊接工装中，压板对焊板的约束是可滑动的——焊板在热膨胀阶段可以沿面内方向滑移，只在法线方向上被约束。把约束改为接触-摩擦边界（摩擦系数0.3，允许面内滑移）之后，残余应力曲线整体下移了约25 MPa，变形从3.2mm降到了2.9mm。

拘束条件的建模精度对残余应力的影响超过了热源参数和相变潜热之和。焊接模拟里最大的误差源往往不是材料模型或者热源模型，而是边界条件和你以为的边界条件之间的差距。

排序和结论

热源参数的敏感性低到这个程度，多少有点反直觉——毕竟焊接模拟的第一步就是热源。但如果目标是残余应力而不是温度场，热源参数的精细标定回报率远不如花时间搞清楚实际的拘束条件。这个结论不只适用于316L TIG焊，对激光焊、电子束焊这类高能束焊接同样成立——热源能量越集中，残余应力受拘束条件的支配性越强。