手机版

ANSYS有限元仿真分析 — 焊接残余应力与变形的热-力顺序耦合模拟

发布时间：2026-07-17 来源：科研学术网

小中大

字号：

焊完之后的变形，谁算谁头疼

焊接残余应力和变形是结构制造中最让工程师头疼的问题之一。一道看似简单的焊缝，由于局部高温加热和随后的快速冷却，在结构中”储存”了复杂的残余应力场——峰值通常达到材料屈服强度水平。这些残余应力与服役载荷叠加，可能导致疲劳寿命的显著降低；焊接变形则直接影响装配精度。

ANSYS有限元仿真分析中的焊接模拟属于顺序热-力耦合分析：首先进行瞬态热分析——模拟焊接热源的移动和热传导过程，获得温度场的历史；然后将温度历史作为热载荷输入到结构分析中，计算应力和变形。

但这条看似清晰的路线在实践中充满了数值陷阱和物理复杂性。我们的项目是一个船舶结构中典型的T型接头——面板厚20毫米、腹板厚15毫米，材料为AH36船用钢，使用CO₂气体保护焊进行多层多道焊接。

热源建模：Goldak双椭球热源的参数怎么定

焊接仿真的起点是热源模型——焊接电弧的热输入如何分布在焊缝区域。Goldak双椭球热源模型是应用最广泛的选择，它用两个半椭球（前半椭球代表电弧前缘的集中加热区，后半椭球代表尾部的扩散热区）来描述焊接热源的三维分布。

双椭球模型有6个几何参数（前/后椭球的三个半轴长度a_f/a_r, b, c），加上热效率η（电弧功率有多大比例实际输入到工件中）。这些参数的选择直接影响温度场的准确性。

我们的参数标定策略是结合焊缝宏观金相来反推。实际的焊缝熔合线（fusion boundary）在金相照片中清晰可见——熔合线内部的金属经历了熔化，外部保持固态。在仿真中，通过调整双椭球参数使计算出的熔池形状（由超过液相线温度1520°C的节点定义）与金相照片中的熔合线吻合。经过几轮迭代，最终确定的参数组合使熔池宽度和高度的仿真值与实测值偏差在10%以内。

一个常被忽视的参数是热源的移动速度。在手工焊中，焊接速度不是恒定的——焊工在起弧、收弧和经过板边时的速度变化可达±20%。在仿真中需要将这些变化纳入考虑——最简单的方法是将焊接速度定义为时间的函数，在不同位置使用不同的移动速度。

相变效应：SSPT不等于”加热再冷却”

钢材在焊接热循环中经历复杂的固态相变（SSPT）：加热时铁素体+珠光体→奥氏体转变（约750-850°C），冷却时奥氏体可能转变为铁素体、珠光体、贝氏体或马氏体——取决于冷却速率。不同的相变产品具有不同的比容（密度），这就产生了相变体积应变——除了热膨胀/收缩之外的一个额外变形来源。

对于AH36船用钢（相当于DH36级），在焊接的典型冷却速率下（10-100°C/s），奥氏体主要转变为贝氏体+少量铁素体——相变体积应变约+0.3%（贝氏体的比容比奥氏体大约0.3%）。这个看似微小的体积应变，在焊缝及热影响区的局部约束条件下，可以产生高达100-200 MPa的相变应力——与热应力同一数量级。

如果在仿真中忽略相变效应，焊接残余应力的预测误差可能高达30%-50%。在ANSYS中，相变效应可以通过定义随温度变化的”热应变+相变应变”曲线来近似处理——在相变温度区间内人为引入一个额外的应变增量。更精确的方法是耦合冶金模型（如Avrami方程描述扩散型相变、K-M方程描述马氏体相变），但这在大型结构焊接仿真中计算成本过高。