X90M管线钢多丝直缝埋弧焊接模拟与分析

X90M管线钢多丝直缝埋弧焊接模拟与分析 Simulation of Multi-wire Longitudinal Submerged Arc Welding of X90M Pipeline Steel 1 first-author 于　磊（1983—），男，硕士，高级工程师，主要从事压力管道、焊接以及无损方面的检验检测及研究等工作。E-mail：467511091@qq.com。 于　磊（1983—），男，硕士，高级工程师，主要从事压力管道、焊接以及无损方面的检验检测及研究等工作。E-mail：467511091@qq.com。 467511091@qq.com 1 1 1 1 2 安徽省特种设备检测院,安徽 合肥 230051 安徽省特种设备检测院,安徽 合肥 230051 Anhui Special Equipment Inspection Institute, Hefei 230051, China Anhui Special Equipment Inspection Institute, Hefei 230051, China 合肥紫金钢管股份有限公司,安徽 合肥 230051 合肥紫金钢管股份有限公司,安徽 合肥 230051 Hefei Zijin Steel Pipe Co., Ltd., Hefei 230051, China Hefei Zijin Steel Pipe Co., Ltd., Hefei 230051, China 以大型通用有限元软件ABAQUS作为模拟分析工具，基于施加在DFLUX子程序中的双椭球热源模型，采用间接耦合方式和“生死单元”算法，模拟X90M钢板内、外三丝直缝埋弧焊接过程，获得了相应的焊接温度场及残余应力场。结果表明：内、外多丝直缝埋弧焊接在焊接工艺参数一定的条件下，每丝的起始焊接时间（或焊丝间距）决定是否形成共同熔池；多丝共熔池焊接因多热源温度场叠加作用剧烈，熔池峰值温度高于多丝单熔池焊接时的，较高的焊接温度场使得焊接过程中能够较为充分地进行冶金反应，且有效降低焊接接头及其附近的冷却速度，使得多丝共熔池焊接形成的焊接接头的焊后残余应力也相对较小。在兼顾质量和生产效率的前提下，X90M多丝直缝埋弧焊接生产实践推荐采用多丝共熔池焊接工艺。 By using large general finite element software ABAQUS as a simulation analysis tool in this paper,based on the double ellipsoidal heat source applied in DFLUX subroutine, the process of internal and external three wire longitudinal submerged arc welding of X90M steel plate is simulated through indirect coupling method and “birth-death element” algorithm,and the welding temperature field and residual stress field are obtained.The results show that when the welding process parameters of internal and external multi wire longitudinal submerged arc welding are certain,the initial welding time (or wire spacing) of each wire determine whether to form a common weld pool; the peak temperature of multi wire eutectic pool welding is higher than that of single wire pool welding due to the severe superposition of multi heat source temperature field，the higher welding temperature field makes the metallurgical reaction in the welding process more fully, and can effectively reduce the cooling rate of the welded joint and its vicinity, resulting in the welded joint formed by multi wire eutectic pool welding has smaller residual stress.On the premise of considering both quality and production efficiency, the multi wire eutectic pool welding process is recommended in the production practice of X90M multi wire longitudinal submerged arc welding. X90M管线钢 多丝直缝埋弧焊 双椭球热源模型 残余应力场 温度场 X90M submerged arc welding residual stress straight seam temperature field 前言 1 随着石油、天然气等能源工业的迅猛发展，世界各国都大规模地进行长距离输送管道的铺设。为了增大输送压力，提高输送效率，高钢级、大口径、大壁厚、高韧性的焊管成为输送管道的首选管材，直缝埋弧焊钢管因质量好、可靠性高 ［ 1 1 - 2 2 1-2 ］ 得到了广泛使用。 随着钢管壁厚的越来越大，传统的单丝焊接已不能满足生产需要。为兼顾生产效率和焊接质量，多丝埋弧焊接工艺开始广泛应用于实际生产。多丝直缝埋弧焊根据每根焊丝的焊接工艺参数和起始焊接时间（或焊丝间距）不同，可形成多丝共熔池焊接和多丝单熔池焊接。多丝共熔池焊接相对于多丝单熔池焊接，其熔池长度长、焊接质量高，如三丝单行纵向排列共熔池焊接时，熔池长度通常在几十毫米甚至上百毫米，且熔池存在时间长，使得焊接冶金反应能够充分进行，在熔池金属凝固前熔池中的气体和熔渣等杂质有充足的时间上浮逸出，焊缝中气孔、夹渣等焊接缺陷很少，焊接一次合格率较高。但多丝埋弧焊高度集中的瞬时热输入，在焊后将产生相当大的残余应力，从而影响焊接结构的可焊性和构件的脆性断裂强度。祁帅 ［ 3 3 ］ 等基于ANSYS软件建立了X80管线钢（Nb-Mo系）三丝埋弧焊的三维热力耦合数值模型，指出温度的分布不均匀是导致残余应力出现的主要原因，增大焊接速度或焊丝间距均会引起焊后纵向残余应力的升高。刘笑笑 ［ 4 4 ］ 等采用ANYSY有限元软件模拟X100管线钢的四丝埋弧焊的温度场和应力场，研究焊接电流、焊接速度、焊丝间距对焊接温度场和残余应分布的影响，结果表明，增加每根焊丝的焊接电流，焊件残余应力水平略有上升，增加焊接速度和焊丝间距时，残余应力峰值水平下降。李丽 ［ 5 5 ］ 等对80 mm E36钢厚板分别进行双丝、三丝、四丝埋弧焊试验，发现多丝埋弧焊接头冲击性能较差的部位在焊缝，随着焊接丝数的增加，焊缝组织中先共析铁素体含量增加，晶粒尺寸增大，冲击性能降低。 本文以X90M钢板内外三丝单行纵向排列直缝埋弧焊接为研究对象，采用ABAQUS有限元软件，研究在焊接工艺参数一定的情况下，内、外每根焊丝起始焊接时间（或焊丝间距）不同时的多丝共熔池、多丝单熔池直缝埋弧焊接温度场和残余应力场，得出起始焊接时间（或焊丝间距）对焊接温度场和焊后残余应力场的影响，旨在为X90M多丝直缝埋弧焊接生产实践提供指导与参考。 试验材料及焊接工艺参数 1 试验用母材为X90M钢板，尺寸为600 mm×300 mm×16 mm，选用MK800GX-Ⅲ焊丝，直径4 mm，母材及焊丝主要化学成分、力学性能如 表1 表1 、 表2 表2 所示。坡口为带钝边的X型坡口，外焊坡口角度为88°，内焊坡口角度为86°，坡口形式如 图1 图1 所示。焊接时先进行内坡口焊接，然后进行外坡口焊接，内、外坡口均采用三丝单行纵向排列直缝埋弧焊接一次成形。内、外坡口焊接的具体焊接工艺参数如 表3 表3 所示。 材质 C Si Mn S P Cr Ni Mo Ti X90M 0.05 0.26 1.83 0.001 0.01 0.12 0.21 0.11 0.004 MK800GX-Ⅲ 0.06 0.27 1.79 0.004 0.01 0.14 0.22 0.41 0.36 材质 屈服强度 R eL /MPa 抗拉强度 R m /MPa 延伸率 A /% X90M 648 757 28 MK800GX-Ⅲ 705 805 23 焊接控制方程及热源模型 1 温度场控制方程 2 热传导方程式是一个重要的偏微分方程，描述一个区域内的温度如何随时间变化。在区域Ω中的热传导控制方程为： <math display="block" id="M1"><mfrac><mrow><mo>∂</mo><mi>T</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>λ</mi></mrow><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mi>ρ</mi></mrow></mfrac><msup><mrow><mo>∇</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub><mi>ρ</mi></mrow></mfrac><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>v</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034348&type= 38.26933289 8.63599968 式中　 <math id="M2"><mi>T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034341&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034340&type= 2.11666679 2.62466669 为温度； <math id="M3"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034357&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034355&type= 1.10066664 2.62466669 为时间； <math id="M4"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034361&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034359&type= 1.60866666 2.62466669 为热导率； <math id="M5"><msub><mrow><mi>C</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034365&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034364&type= 3.64066648 2.62466669 为比热容； <math id="M6"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034347&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034346&type= 1.77800000 3.38666677 为密度； <math id="M7"><msup><mrow><mo>∇</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034373&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034379&type= 3.72533321 2.62466669 为拉普拉斯算子； <math id="M8"><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>v</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034376&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034382&type= 6.34999990 7.61999989 为内热源强度。 应力场控制方程 2 在区域Ω中的力平衡方程为： <math display="block" id="M9"><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034394&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034386&type= 8.29733276 2.62466669 式中　 σ ij 为已包括热应力项的应力分量。 在区域Ω中的热应变平衡方程为： <math display="block" id="M10"><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow><mrow><mi>t</mi></mrow></msup><mo>=</mo><mi>α</mi><mo stretchy="false">(</mo><mi>T</mi><mo>-</mo><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034395&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034389&type= 23.28333473 3.38666677 式中　 <math id="M11"><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow><mrow><mi>t</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034412&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034408&type= 3.89466691 2.62466669 为热应变张量； <math id="M12"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034402&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034398&type= 1.94733346 2.62466669 为线膨胀系数； <math id="M13"><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034424&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034404&type= 3.21733332 2.62466669 为环境参考温度； <math id="M14"><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034439&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034422&type= 3.04800010 2.62466669 为 <math id="M15"><msub><mrow><mi>δ</mi></mrow><mrow/></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034433&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034440&type= 2.62466669 2.62466669 算子。 应力应变之间的本构关系为： <math display="block" id="M16"><mi mathvariant="normal">d</mi><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msub><mrow><mi>D</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><mi mathvariant="normal">d</mi><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub><mo>-</mo><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>p</mi></mrow></msup><mo>-</mo><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>c</mi></mrow></msup><mo>-</mo><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>t</mi></mrow></msup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034454&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034436&type= 49.10666656 3.38666677 式中　 <math id="M17"><msub><mrow><mi>D</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034459&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034456&type= 3.97933316 2.62466669 为热应变张量； <math id="M18"><mi mathvariant="normal">d</mi><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034446&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034462&type= 4.91066647 2.62466669 、 <math id="M19"><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>p</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034451&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034466&type= 6.43466663 2.62466669 、 <math id="M20"><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>c</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034473&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034484&type= 6.34999990 2.62466669 、 <math id="M21"><mi mathvariant="normal">d</mi><msup><mrow><msub><mrow><mi>ε</mi></mrow><mrow><mi>k</mi><mi>l</mi></mrow></msub></mrow><mrow><mi>t</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034476&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034474&type= 6.01133299 2.62466669 分别为总应变、塑性应变、蠕变应变、热应变。 双椭球热源模型 2 多丝直缝埋弧焊接因其熔深较大，故采用双椭球热源模型来近似描述熔池形貌和尺寸，表达式为 <math display="block" id="M22"><msub><mrow><mi>q</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mroot><mrow><mn mathvariant="normal">3</mn></mrow><mrow/></mroot><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>Q</mi></mrow><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>b</mi><mi>c</mi><mi mathvariant="normal">π</mi><mroot><mrow><mi mathvariant="normal">π</mi></mrow><mrow/></mroot></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034493&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034490&type= 70.01933289 11.26066685 <math display="block" id="M23"><msub><mrow><mi>q</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mroot><mrow><mn mathvariant="normal">3</mn></mrow><mrow/></mroot><msub><mrow><mi>f</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>Q</mi></mrow><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>b</mi><mi>c</mi><mi mathvariant="normal">π</mi><mroot><mrow><mi mathvariant="normal">π</mi></mrow><mrow/></mroot></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">(</mo><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msup><mrow><mi>c</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034496&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27034494&type= 70.01933289 11.26066685 式中　 a 1 、 a 2 、 b 、 c 分别为双椭球体积热源三维形状分布参数； f 1 、 f 2 分别为总的焊接热输入功率在熔池前后两部分的能量分布系数，且有 f 1 + f 2 =2。 设置内焊三丝形成的双椭球热源形状参数比例关系为 a 1 ∶ a 2 ∶ b ∶ c =1∶3∶（1.22~1.68）∶（1.45~1.06），外焊三丝形成的双椭球热源形状参数比例关系为 a 1 ∶ a 2 ∶ b ∶ c =1∶3∶（1.38~1.68）∶（1.54~1.22）。 建立有限元模型 1 材料热物性能参数 2 假设母材与焊材的材料性能参数相同， X90M对应的具体材料性能参数如 图2 图2 所示。设泊松比μ为0.3，且不随温度变化而变化。 有限元网格模型 2 建立与试验用钢板试板尺寸完全相同的三维有限元模型。为兼顾模拟求解精度和运算效率，模型网格的划分策略为：对焊缝及近缝区域的网格划分较细，远离焊缝及近缝区域的网格则划分较粗，中间采用1∶3、1∶2网格过渡方式。试验用钢板试板及有限元网格模型如 图3 图3 所示，该有限元网格模型共有单元数64 000个，共有节点数71 563个。 假定及边界条件 2 假设材料服从热弹塑性和双线性本构模型以及Von Mises屈服准则 ［ 6 6 -7］ 。模拟采用间接耦合，总体思路为：焊接过程中必须考虑温度场对应力场的影响，而忽略应力场对温度场的微小影响，在实际模拟计算过程中先进行温度场的计算，然后将计算结果作为初始条件导入相应应力场的分析计算过程。 为较好地模拟焊缝金属的填充过程，运用“生死单元”算法 ［8］ 。实施步骤为：首先将内、外焊道处焊缝金属总体作为一个大的单元集合全部被“杀死”，然后将内焊道处焊缝、外焊道处焊缝分别作为独立单元依次“激活”。 在焊接阶段，双椭球热源模型的加载是通过FORTRAN语言编写DFLUX子程序实现，首先基于内、外焊三丝每根焊丝的焊接工艺参数，利用双椭球热源模型方程计算模型上各节点的热流密度值，并施加于所选节点上，随着电弧中心位置的移动，重复在各节点上施加相应的热流密度，从而实现热源的移动。 假设焊件初始温度为20 ℃，将对流换热、辐射换热效果叠加考虑并将总的换热系数定为20 W/（m 2 ·℃），在焊接应力场的求解过程中设置相应的力学边界条件使模型不产生刚体位移以达到分析收敛 ［9］ 。 计算结果及分析 1 焊丝间距对峰值温度场的影响 2 温度场的数值模拟是残余应力场数值模拟的前提，温度场对残余应力和残余变形影响极大 ［10］ 。三丝直缝埋弧焊每根焊丝对应独立弧焊电源，故在内、外焊接过程中三丝中的每丝均形成独立熔池，而独立熔池会在时空维度相互影响，在满足一定条件的情况下形成共熔池状态，能否形成共熔池状态取决于每丝的焊接工艺参数和起始焊接时间（或焊丝间距）。 图4 图4 为在 表3 表3 焊接工艺参数范围内（焊丝间距除外），每根焊丝起始焊接时间间隔1 s（内焊焊丝间距均为21.7 mm，外焊焊丝间距均为20 mm）、3 s（内焊焊丝间距均为65.1 mm，外焊焊丝间距均为60 mm）、5 s（内焊焊丝间距均为108.5 mm，外焊焊丝间距为100 mm）、8 s（内焊焊丝间距为173.6 mm，外焊焊丝间距均为160 mm）的焊接峰值温度场云图，该峰值温度场是内外焊接连续过程中的峰值温度，内焊完成后间隔1 min进行了外焊。 1,2;3,4; 2 焊接位置 电源极性 焊接电流 I /A 电弧电压 U /V 焊接速度 v /（m·min -1 ） 焊丝倾角/（°） 焊丝间距/mm 3 内坡口 前丝 DC+ 580~600 32~34 1.3 90 3 前中丝、中后丝：18~28 中丝 AC 530~550 35~37 1.3 85 后丝 AC 490~510 38~40 1.3 80 3 外坡口 前丝 DC+ 630~650 32~34 1.2 90 3 前中丝、中后丝：18~28 中丝 AC 580~600 35~37 1.2 85 后丝 AC 530~550 38~40 1.2 80 由 图4 图4 可知，在焊接工艺参数一定的条件下，每丝起始焊接时间间隔1 s焊接过程中形成了共同熔池，为多丝共熔池，熔池峰值温度为2 625 ℃；焊丝起始焊接时间分别间隔3 s、5 s、8 s时，每根焊丝形成的熔池依旧独立，为多丝单熔池，熔池峰值温度分别为2 080 ℃、1 864 ℃、1 687 ℃。多丝共熔池因是每根焊丝独立熔池互相叠加形成，故峰值温度高于多丝单熔池；无论是否形成共熔池状态，焊接过程中熔池的峰值温度随每根焊丝起始焊接时间间隔（焊丝间距）增大而降低，如 图5 图5 所示。 选取外焊道焊缝表面中心同一节点（节点编号：19588）处，该节点经历的焊接热循环如 图6 图6 所示。从 图6 图6 可以看出：（1）焊接热源的输入是一个极不均匀的过程，先受热位置处的温度急剧升高，然后随着热源的离开温度逐渐降低，这种局部加热—熔化—冷却的剧烈变化过程决定了焊后复杂的残余应力场；（2）焊丝起始焊接时间分别间隔1s、3 s、5 s、8 s时，对应的热循环峰值温度逐渐降低，最高为2 373 ℃、最低为1 485 ℃。 焊接残余应力 2 电弧焊接过程中非均匀、极不平衡的热量输入是导致焊后存在残余应力的根本原因，焊后残余应力分布场基本由焊接温度场决定，内、外焊每丝起始焊接时间间隔1s、5s焊接的Von Mises等效残余应力分布如 图7 图7 所示。 1,2; 由 图7 图7 可知，内、外焊时每丝起焊间隔1 s、5 s时焊接残余应力分布特征为焊缝及其近缝区全部为拉应力，从横向上看远离该区域残余应力逐渐降低；每丝起焊间隔1 s的多丝共熔池焊接的残余应力较间隔5 s时的多丝单熔池焊接更小，这是因为前者较高的熔池温度和较长的电弧总体停留时间使得其冶金反应充分，这样在冷却过程中能够有效降低焊缝及近缝区的临界冷却速度，利于形成淬硬倾向较小的组织，所以焊后残余应力相对较小。较小的焊接残余应力可以降低直缝埋弧焊接缺陷产生的概率。 图8 图8 为内、外焊每丝起焊间隔1 s形成多丝共熔池焊接时Von Mises等效残余应力沿试板起焊端面厚度方向（自内焊中心至外焊中心）的分布，可以看出，从内表面焊缝中心到外表面焊缝中心残余应力逐渐降低。 图9 图9 为内、外焊每丝起焊间隔1 s形成多丝共熔池焊接的侧向弯曲后的试样，原始弯曲试样为矩形试样，尺寸为350 mm（长）×38.1 mm（宽）×16 mm（厚），弯轴直径为90 mm，可以看出，侧向弯曲180°后未发现任何裂缝及其他缺陷，塑性良好。 结论与展望 1 （1）X90M管线钢内、外多丝直缝埋弧焊接在每丝焊接工艺参数一定的条件下，是否形成多丝共熔池焊接取决于每丝起焊时间（或间距）；多丝共熔池因每丝熔池叠加形成的峰值温度较多丝单熔池高；焊接过程中熔池峰值温度随每丝起焊时间（或间距）增大而降低。 （2）X90M管线钢内、外多丝共熔池直缝埋弧焊接相对于多丝单熔池焊接其峰值应力较低，主要是因为多丝共熔池焊接时较高的焊接温度场在冷却过程中其冷却速率较低，焊后接头性能较好，在X90M等高钢级钢板、钢管生产实践中推荐采用多丝共熔池直缝埋弧焊接工艺。 （3）本文研究对象为X90M钢板试板的焊接，钢板焊接的温度场及应力场不完全等同于钢管的温度场及应力场，对X90M管线钢直缝埋弧焊的温度场及应力场研究仅作参考；另本文未模拟在焊丝起始焊接时间（或间距）一定情况下，改变焊接工艺参数后给温度及应力场分布带来的影响，对此问题在今后需进行进一步研究。 参考文献 宋艾玲 ， 梁光川 ， 王文耀 . 世界油气管道现状与发展趋势 ［J］. 油气储运 ， 2006 （ 10 ）： 3，8 - 13，69 . 世界油气管道现状与发展趋势 Song ailing ， Liang Guangchuan ， Wang Wenyao . Current situation and development trend of oil and gas pipelines in the world ［J］. Oil and gas storage and transportation ， 2006 （ 10 ）： 3， 8 - 13， 69 . Current situation and development trend of oil and gas pipelines in the world 晁利宁 ， 毕宗岳 ， 余晗 ， 等 . 直缝埋弧焊管焊接残余应力的数值模拟 ［J］. 焊管 ， 2015 ， 38 （ 8 ）： 24 - 28，33 . 直缝埋弧焊管焊接残余应力的数值模拟 Chao lining ， Bi Zongyue ， Yu Han ， et al . Numerical simulation of welding residual stress of LSAW pipe ［J］. Welded pipe ， 2015 ， 38 （ 8 ）： 24 - 28， 33 . Numerical simulation of welding residual stress of LSAW pipe 祁帅 ， 严春妍 ， 史志丹 ， 等 . X80管线钢三丝埋弧焊纵向残余应力数值模拟 ［J］. 电焊机 ， 2015 ， 45 （ 04 ）： 153 - 157，169 . X80管线钢三丝埋弧焊纵向残余应力数值模拟 Qi shuai ， Yan Chunyan ， Shi Zhidan ， et al . Numerical simulation of longitudinal residual stress of X80 Pipeline Steel by three wire submerged arc welding ［J］ Electric welding machine ， 2015 ， 45 （ 04 ）： 153 - 157，169 . Numerical simulation of longitudinal residual stress of X80 Pipeline Steel by three wire submerged arc welding 刘笑笑 ， 颜志勇 ， 倪川皓 . X100管线钢四丝埋弧焊温度应力场的数值分析 ［J］. 矿冶工程 ， 2019 ， 39 （ 03 ）： 123 - 127 . X100管线钢四丝埋弧焊温度应力场的数值分析 Liu Xiaoxiao ， Yan Zhiyong ， Ni Chuanhao . Numerical analysis of temperature stress field of X100 Pipeline Steel by four wire submerged arc welding ［J］ Mining and Metallurgy Engineering ， 2019 ， 39 （ 03 ）： 123 - 127 . Numerical analysis of temperature stress field of X100 Pipeline Steel by four wire submerged arc welding 李丽 ， 陈亚魁 ， 吴君明 ， 等 . 海洋平台用特厚板大线能量焊接接头力学性能研究 ［J］. 电焊机 ， 2018 ， 48 （ 11 ）： 27 - 30，35 . 海洋平台用特厚板大线能量焊接接头力学性能研究 LI Li ， CHEN Yakui ， WU Junming ， et al . Study on mechanical properties of large line energy welded joints of extra thick plates for offshore platforms ［J］ Electric welding machine ， 2018 ， 48 （ 11 ）： 27 - 30，35 . Study on mechanical properties of large line energy welded joints of extra thick plates for offshore platforms 张玉凤 ， 霍立兴 ， 董俊慧 . 低碳钢管管道焊接残余应 低碳钢管管道焊接残余应 10.7512/j.issn.1001-2303.2022.03.06 TG457 A 1001-2303（2022）03-0046-07 于磊，李志宏，姜璐，等.X90M管线钢多丝直缝埋弧焊接模拟与分析［J］.电焊机，2022，52（3）：46-52. 于磊，李志宏，姜璐，等.X90M管线钢多丝直缝埋弧焊接模拟与分析［J］.电焊机，2022，52（3）：46-52. YU Lei, LI Zhihong, JIANG Lu, et al.Simulation of Multi-wire Longitudinal Submerged Arc Welding of X90M Pipeline Steel[J].Electric Welding Machine, 2022, 52(3): 46-52. YU Lei, LI Zhihong, JIANG Lu, et al.Simulation of Multi-wire Longitudinal Submerged Arc Welding of X90M Pipeline Steel[J].Electric Welding Machine, 2022, 52(3): 46-52. 2021-07-19 2021-10-26 2022-03-30 2022-04-26 电焊机 3 52