管道环焊缝焊接及焊后热处理过程的有限元模拟.docxVIP

管道环焊缝焊接及焊后热处理过程的有限元模拟

刘永，王苹，方洪渊

哈尔滨工业大学先进焊接与连接国家重点实验室（哈尔滨，150001）

摘要：使用大型通用有限元软件ABAQUS计算了管道环焊缝的焊接过程，获得焊接过

程的温度场和应力应变场，并对焊接残余应力的分布状态进行了分析；之后进行了管道

焊后热处理（PWHT）工艺过程的有限元计算，获得了PWHT后的残余应力场，并分析

了残余应力消除机理。结果表明：管道环焊缝周向焊接残余应力整体为较高的拉应力，

峰值达到甚至超过屈服强度；而轴向应力在内外壁均为拉应力，在壁厚中部为压应力，

即轴向应力在壁厚方向呈现“自平衡”的趋势。经PWHT后，各方向残余应力均得到大幅

度的降低；PWHT过程的塑性应变和高温蠕变是残余应力消除的力学机制，后者起到主

要的作用；因而在制定PWHT工艺时，加热温度应高于蠕变激活温度，以有效消除构件

的焊接残余应力。

关键词：ABAQUS；管道环焊缝；焊接残余应力；焊后热处理

1热-弹塑性理论模型[1]

1.1焊接热过程的计算

在焊接过程中，热分析的控制方程为：

??????

c(x,y,z,t)q(x,y,z,t)Q(x,y,z,t)

T

?t

（1）

式中：?为材料密度；c为材料的比热容；T为瞬时温度；q为热流矢量；Q为热源产热

率；x，y，z为参考坐标系中的坐标，t为时间，?为拉普拉斯算子。文中使用的是轴对

称模型，选用等密度热源模型。

在焊接和冷却过程中，工件由于对流和辐射而存在热量散失。其中，对流换热在低

温区域起主要作用，辐射换热在高温区域起主要作用，例如在800℃时约占总放热量的

80%。对流换热遵循牛顿定律：

q??h(T?T)（2）

cfsur0

式中，hf为对流表面换热系数；Tsur为固体表面温度；T0为环境温度。辐射换热遵循斯

蒂芬-玻尔兹曼定律：

q???C0(T?T0)（3）

44rsur

C???Wmm?K，适用于绝对黑体；?为黑度系数，??1。式中：1424

05.6710/()1.2焊接力学过程

热弹塑性分析中，总的应变增量??total符合下列公式：

??total???e???p???th???tr???cr（4）

式中：等号左侧为总应变增量；等号右边依次为弹性应变增量、塑性应变增量、热应变

增量、相变应变增量和蠕变应变增量。对于一般材料而言，焊接过程中固态和蠕变对焊

-1-

接应力应变的演化作用较小，可将其忽略，因而焊接过程的的热弹塑性增量公式为：

???????????（5）

totalepth

而焊后热处理计算中需要将高温蠕变机制考虑进去，因此PWHT计算过程中总的应

变增量??total为：

??total???e???p???th???cr（6）

热应变增量??th????T；蠕变过程使用Norton公式，即幂律蠕变模型：

??A?t（7）

nm

2焊接过程的有限元模拟

2.1管道有限元模型及材料参数

管道的内径??1.5m，壁厚为22mm；坡口角度44°，钝边1mm，装配间隙3mm，

余高2mm，共6层11道，焊接方法为MIG焊，建立的轴对称管道环焊缝的几何模型如

图2-1a所示。选用四边形缩减积分单元CAX4RT，焊缝局部有限元网格如图2-1b所示。

（a）焊缝几何模型（b）有限元网格

图2-1管道焊接有限元模型

管道材料为普通碳素钢，并假设焊缝材料性能与母材相同，材料硬化模型选择理想

弹塑性，并考虑材料非线性。随温度变化的材料参数如图2-2所示。热学边界条件设置

对流与辐射，初始温度和环境温度均设置为20℃。

（a）热-物理材料参数（b）热-力学材料参数

图2-2随温度变化的非线性材料参数[2]

2.2焊接温度场及热循环

图2-3为焊接过程中第5道和第6道焊缝截面的温度场云图，可见等密度热源的加

-2-

热温度为材料熔点（1500℃）。图2-4为管道坡口内外壁边缘位置的热循环曲线。

（a）第5道焊接

（b）第6道焊接

图2-3焊道截面的温度场云图

1600

1400

1200

1000

800

600

400

200

0

020004000600080001000012000

时间t/s

（a）内壁坡口边缘

-3-

1600

1400

1200

1000

800

600

400

200

0

020004000600080001000012000

时间t/s