厚板多层多道对接焊残余应力轮廓法测量
及热-弹-塑性有限元分析
张庆亚1，卓子超1，周宏2，刘建成3，王江超1
（1.华中科技大学船舶与海洋工程学院，武汉430074；2.江苏科技大学船舶与海洋工程学院，
江苏镇江212100；3.招商局重工（江苏）有限公司，江苏南通226100）
摘要：船体外板结构多采用多层多道焊工艺建造，焊接产生的残余应力复杂，易导致焊接结构断裂和疲劳失效。本文采用轮廓法与基于并行计算技术的热-弹-塑性有限元研究Q235厚板多层多道对接接头内部残余应力分布及其变化过程，预测结果与测量结果吻合较好。预测结果表明，纵向残余应力在焊缝区为拉应力，沿接头宽度方向逐渐减小最终转变为压应力，正面焊缝中部区域拉应力值明显降低；横向残余应力在焊缝区上表面及背面填充处为拉应力，沿接头宽度方向逐渐降低，在正面焊缝中部区域为压应力；Z 向残余应力主要存在于焊缝区，正面焊缝以压应力为主，背面焊缝以拉应力为主。厚板多层多道焊接过程中，残余应力变化是由于后一道焊缝对已焊的焊缝起到热处理作用导致的，焊缝接头残余应力分布由最后一道焊缝决定。
关键词：多层多道对接焊；焊接残余应力；轮廓法；热-弹-塑性有限元
中图分类号：U661.4文献标识码：A doi:10.3969/j.issn.1007-7294.2021.05.011Welding residual stress in multi-pass butt joint of thick plate through contour method and thermal-elastic-plastic finite element analysis
ZHANG Qing-ya 1,ZHUO Zi-chao 1,ZHOU Hong 2,LIU Jian-cheng 3,WANG Jiang-chao 1(1.School of Naval Architecture and Ocean Engineering,Huazhong University of Science and Technology,Wuhan 430074,China;2.School of Naval Architecture and Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang 212100,China;3.China Merchants Group Heavy Industry (Jiangsu)Co.Ltd.,Nantong 226100,China)
Abstract ：The hull plate structure is often fabricated by multi-pass welding process.However,the welding-induced residual stress will lead to abrupt fracture and fatigue-related failure of welded structure.In this pa⁃per,the internal residual stress distribution in a multi-pass butt joint of a Q235thick plate was investigated by using the contour method and Thermal Elastic Plastic Finite Element Analysis (TEP FEA).The measured results were approximated with the predicted results.The longitudinal residual stress in the weld vicinity was tensile stress,it decread gradually towards the ba metal and finally turned into compressive stress.The stress value in the middle ction of the front weld decread obviously.The transver residual stress was the tensile stress near the surfac
e and then gradually decread along the weld width.The compressive stress 第25卷第5期船舶力学Vol.25No.52021年5月
Journal of Ship Mechanics May 2021
文章编号：1007-7294（2021）05-0627-10收稿日期：2020-11-28
基金项目：工业与信息化部高技术船舶科研专项（159****9518）
作者简介：张庆亚（1990-），男，博士研究生，E-mail:******************* ；
王江超（1983-），男，博士，副教授，通讯作者，E-mail:************** 。
628船舶力学第25卷第5期was in the middle ction of the front weld.The Z-direction residual stress mainly existed in the weld zone.
The compressive stress was dominant in the front weld and the tensile stress was dominant in the back weld. During the multi-pass butt welding process of the thick plate,the latter weld had heat treatment effect on the former weld,leading to the change of residual stress.The residual stress distribution of the weld joint was de⁃termined by the last weld.
Key words:multi-pass butt welding;welding residual stress;contour method;
thermal elastic plastic finite element analysis(TEP FEA)
0引言
船舶建造过程中，船体外板结构往往采用厚板多层多道焊接而成，多次焊接热循环导致焊接接头内部残余应力分布复杂，此外焊接残余应力易导致焊接结构断裂，对焊接结构的完整性和使用性能危害极大[1-3]。因此，准确获取厚板多层多道焊接头内部焊接残余应力分布是焊接构件寿命评估的基础。
虽然厚板多层多道焊接结构残余应力的测定困难，但一些有效的方法已被应用，如小孔法、衍射法等，其中小孔法测量精度有限，衍射法测量极其昂贵[4]，且只能测量焊缝近表面处残余应力。Prime[5]基于叠加原理提出轮廓法测试构件内部应力，诸多研究表明该方法能完整得到焊缝接头某一截面上的内部残余应力分布云图[6-10]。由于厚板焊接时每一道焊缝均会产生残余应力，采用试验方法难以合理地评估接头内部的残余应力变化过程。而对于已知的焊接参数以及接头形式，热-弹-塑性有限元法可以预测焊接瞬态温度、残余应力和变形，成为焊接接头及小型焊接结构的常用焊接数值模拟方法[11-13]。周宏等[14]采用基于生死单元技术的热-弹-塑性有限元法模拟70mm厚EH47高强钢多层多道焊残余应力，计算结果与小孔法实验测量吻合较好。孙加民等[15]对板厚为30mm的Q390高强钢多层多道接头的焊接残余应力进行了数值模拟，采用小孔法实测了接头的焊接残余应力，证明焊接
残余应力有较高的计算精度。黄俊等[16]采用SYSWELD软件预测了超高强钢厚板多层多道对接焊残余应力，小孔法验证了数值模拟结果的合理性。此外，由于厚板多层多道焊焊道密集，传统的热-弹-塑性分析消耗大量的计算时间，具有一定的局限性。近年来，迭代子结构法和并行计算技术大大提高了厚板多层多道焊数值模拟的计算效率[17-20]。
厚板多层多道焊残余应力数值模拟及测量已取得显著成果，然而对于厚板接头内部的残余应力分布及其变化过程关注较少。本文针对船用Q235钢多层多道焊对接接头，采用轮廓法测量与基于并行计算技术的热-弹-塑性有限元研究接头内部残余应力分布，并分析焊接过程中残余应力的变化过程，为评定船用钢多层多道对接焊接头的断裂性能及结构寿命提供一定的理论依据。
1数值计算方法
本文采用热-弹-塑性有限元法研究焊接过程中的传热以及应力问题，采用并行计算技术提高热-弹-塑性有限元计算效率。
1.1热-弹-塑性有限元分析
热-弹-塑性有限元计算分析主要包括热分析和力学分析两个过程，其中热分析结果对力学分析结果具有决定性的作用，反过来力学分析结果对热分析结果的影响很小，可忽略不计[11-13,19]。本文采用非
耦合的热-力分析过程，即只考虑焊接瞬态温度对力学分析过程的影响。热分析过程为使用热传导理论和材料的热物理性能参数求解非线性传热方程，进而得到焊接瞬态温度场；力学分析过程是以
热分析得到的温度场为载荷，同样考虑到材料的高温力学性能参数，从而得到焊缝接头的应变、应力和变形。
1.2并行计算技术
由于多层多道焊接头有限元模型节点及单元数量庞大，而传统的热-弹-塑性有限元分析多采用串行模式，计算耗时过长，难以满足复杂的多层多道焊求解。本文基于Dell Power Edge T420服务器，Ubuntu14.04版本操作系统以及Intel Compiler 编译器，实现热-弹-塑性有限元计算程序的并行化。在焊接温度场以及力学响应分析的计算程序中，对于大型矩阵求解的子循环，通过调用OpenMP 代码，编译器自动将程序并行化处理；而对于数据读入和处理过程和分支判断则采用串行处理提高代码执行效果[20]。在进行厚板多层多道热分析以及力学分析过程时，实时调用多核及多线程同时进行计算，可大大提高计算效率。
2焊接及残余应力测量
焊接实验包括焊前准备、焊接参数记录以及焊缝检测；焊后，采用轮廓法测量焊缝接头的残余应力分布。
2.1焊接实验
焊接实验所用母材为30mm 厚船用Q235B 低碳钢板，焊缝长度为300mm ，焊接坡口如图1（a ）所示；焊接方法为焊条电弧焊，填充金属为直径4mm 的J507焊条。焊前，烘干焊条并打磨坡口；焊接时的焊接工艺参数为打底焊电流150~170A ，填充焊电流170~190A ，盖面焊电流160~180A ，焊接电压为26~27V ，焊接速度为180~230mm/min 。此外，焊接时先焊接正面坡口，碳弧气刨清根后再焊接背面坡口。
焊缝成形如图1（b ）所示，可以看出焊缝成形良好，无咬边、表面气孔、未熔合等缺陷。采用线切割垂直切割焊缝试样，依次用600目以及1000目金相砂纸打磨焊缝区域，用4%硝酸酒精溶液腐蚀打磨区域获得焊缝宏观金相。
2.2残余应力测量根据Bueckner 叠加原理[5]，对于待测构件任意平面的应力，将其完整切为两半，因切割面处应力释放导致轮廓发生变形；若施加的外力将变形后的切割面恢复到原始的平面状态，那么所得的应力即等效为该平面的原始残余应力。轮廓法测量顺序为切割、轮廓变形测量、数据处理以及应力重构分析。具体流程如下：
（1）切割：试样切割质量决定着最终的应力分布云图。为提高切割面光洁度，本文通过去除表面第5期张庆亚等：厚板多层多道对接焊残余应力轮廓法测量…
629（a ）坡口
（b ）焊缝成形
图1厚板对接接头坡口及焊缝成形Fig.1Groove feature and weld appearance of butt welded joint
余高实现贴面切割加工；然后将待测接头对称装夹到慢走丝线切割机床中，水平校准后通入去离子水充满切割空间；设定切割工艺参数。为避免切割时应力释放引起的试样移动，切割前需采用夹具固定试样，焊缝试样约束如图2（a ）所示。本文采用的切割设备为日本Sodick AQ400LS 慢走丝线切割机床，实验中切割速度为0.2mm/min 。
（2）轮廓变形测量：切割完成后，准确获取切割面变形值是轮廓法流程中重要的环节。切割面变形值可通过非接触式的光学仪器测量或者非接触式的三坐标测量仪测量获得，其中三坐标测量仪可得到均匀的测量点变形数据，成为测量切割面轮廓变形数据的常用测量设备。本文中，切割面轮廓面外变形所采用的测量设备为Hexagon 三坐标测量仪，见图2（b ），设备测量精度可达1μm 。测量过程中采用线扫描点接触测量，以保证每个测量点的精度，测量点间距为1mm 。
（3）数据处理：因试样表面并非理想的光洁度，测量得到的异常值易造成计算过程中应力局部突变，因此需对测量数据进行处理。数据处理流程包括数据对齐、求和以及拟合。数据对齐是为了保证轮廓面数据点起始点、测量点间距和数据分布一致；然后通过两个表面的变形数据平均以消除轮廓切应力的影响，获得纵向残余应力释放引起的变形量。采用Sgolay 方法去除噪音数据，实现测量数据光滑处理。该方法在数据光滑滤波的同时也对数据进行拟合，处理后的轮廓数据如图2（c ）所示。
（4）应力重构分析：应力重构分析是轮廓法流程中的最后一步。本文采用ABAQUS 建立1/2焊缝有限
元模型如图2（d ）所示，所用单元类型为C3D8R ；为了反映切割面轮廓变形，将变形放大185倍。为了保证有限元求解的精度和速度，靠近轮廓面的网格较密，而远离轮廓面的网格较稀疏；切割面单元尺寸与测量点间距相等，为1mm 。图中的红色箭头表示约束条件，约束刚体移动。采用线性静力学分析，只考虑材料的弹性行为，弹性模量取210GPa ，泊松比取0.3，通过施加反向的轮廓变形进行逆有限元分析进而求解得到切割面焊接残余应力。
图3所示为焊缝内部纵向残余应力重构分析结果，纵向残余应力在焊缝区域为拉应力，正面焊缝中部拉应力明显降低，最大拉应力值为269.9MPa 。残余应力由焊缝区拉应力沿接头宽度方向逐渐减小并转变为压应力。此外，云图左侧灰色区域压应力突变是由于边缘效应导致的[8]。630船舶力学第25卷第5
期（a ）慢走丝线切割设备及试样约束（b ）三坐标测量仪
（c ）轮廓变形
（d ）应力重构有限元模型
图2轮廓法实验流程Fig.2Experimental procedure of contour method
3热-弹-塑性有限元分析
采用基于并行计算技术的热-弹-塑性有限元法预测厚板多层多道对接接头温度场以及内部残余应力分布，并分析厚板多层多道焊接过程中的残余应力变化过程。
3.1有限元模型
图4所示为热-弹-塑性有限元模型及焊道布置，该模型包括25172个节点，23130个单元。此外，在进行焊接力学热-弹-塑性有限元计算时，接头有限元模型网格尺寸对焊接温度场、塑性应变、焊接变形以及残余应力的计算精度起着重要的作用。精细网格虽可以得到较好的预测结果，但同时将消耗大量的计算机资源和计算时间。因此，为了平衡计算精度和计算机资源消耗，通常焊接区（焊缝及热影响区）采用细网格，远离焊缝的区域采用粗网格。在热-弹-塑性有限元分析过程中，材料的热物理性能参数对计算结果的影响很大。本文采用的Q235相关热物理性能参数通过JMatPro 专业软件获得，如图5所示。
第5期张庆亚等：厚板多层多道对接焊残余应力轮廓法测量…
631
图3纵向残余应力分布云图
Fig.3Contour of longitudinal residual stress distribution 热传导(10J/(mm ·s))
比热(Ｊ/(kg ·℃））
换热系数(10W/(mm ·K))
密度(10g/mm )
材料传热性能参数温度／℃泊松比
杨氏模量(10Pa)
屈服强度(10Pa)
线膨胀系数(10/℃)温度/℃
材料力学性能参数图4多层多道对接接头热-弹-塑性有限元模型
Fig.4TEP FEA finite element model of multi-pass butt joint
图5Q235材料热物理性能参数Fig.5Material properties of
Q235（a ）材料传热性能（b ）材料力学性能
