旁路等离子-MIG 复合电弧及耦合熔池作用机理与数值分析

王子然，苗玉刚，王林，马旭伯，魏超，张本顺

(1.哈尔滨工程大学,水下机器人技术国家级重点实验室,哈尔滨,150001；2.中国机械总院集团哈尔滨焊接研究所有限公司,哈尔滨,150028；3.江苏自动化研究所,连云港,222006)

0 序言

目前在工业制造领域，包括船舶建造、航空航天、轨道交通等，对高质量和高效率的新型焊接技术需求越来越迫切[1].实际生产中，等离子弧焊(plasma arc welding，PAW)和熔化极惰性气体保护焊(metal inert gas welding，MIG)由于前者的高能量密度[2]和后者的高柔性度、高效率[3]等优点在工业领域广泛应用.旁路等离子体-MIG 同轴复合电弧焊是一种新型工艺，既结合了这两种方法的优点，同时又可以通过调节旁路电流的大小精确控制电弧能量，从而能够适应不同加工场景及工作需求.

近年来，等离子-MIG 复合电弧焊技术在国内外科研与工程领域屡见报道.Bai 等人[4]对其熔滴过渡特性进行了研究，证明了亚射流过渡和射滴过渡是这种工艺最稳定的熔滴过渡方式；Yang 等人[5-6]则对其工艺特性进行了试验研究，发现外侧的等离子电流有助于增加焊接熔深,另外还对等离子体-MIG 复合电弧的作用机制进行了研究，发现随着外部等离子电流的增加，电弧的稳定性得到加强，从而有助于进一步提高焊接质量；而在基于旁路电流的等离子体-MIG 复合工艺方面，Miao 等人[7]研究了其在加热和冷却过程中应力和变形的演变过程和作用机理，发现这种方法有利于缓解焊后残余应力和变形；需要指出的是旁路电流的概念是由美国肯塔基大学Zhang 教授[8]首先提出，他们使用钨极作为分流电极分走MIG 电弧的电流，以减少母材的热输入；Huang 等人[9]提出将这一思想与等离子热丝工艺方法相结合进行镍基合金增材制造，发现熔覆层的高度随着旁路电流的增加而增加.前期的研究报道主要集中在等离子-MIG 同轴复合技术在改善焊接增材过程中材料成形、电弧及熔滴特性等方面的研究，但是基于旁路电流的等离子-MIG 复合电弧和熔池的物理作用机制研究却未见报道.

文中借助商用流体计算软件Fluent，采用数值模拟[10]和物理试验相结合方法，针对旁路电流对等离子体-MIG 复合电弧和耦合熔池影响机理进行深入研究，相关研究结果不仅丰富了电弧焊接理论及方法，同时为焊接过程能量精确控制提供更多可行性.

1 数值模型

1.1 原理与计算域

图1 是旁路等离子体-MIG 复合工艺的原理和2D 计算模型.旁路等离子-MIG 复合工艺的设计原理是通过与电源负极相连的等离子焊枪铜嘴和与电源正极相连的焊丝间形成分流弧，分走焊丝与母材间主路电弧的电流，如图1a 所示，为了防止铜嘴过热烧损，铜嘴采用循环水冷却；图1b 为旁路等离子-MIG 复合焊接电弧与熔池的作用原理；图1c 为根据原理设计的数值模型包括焊丝区、铜嘴区、电弧区和母材区.模型的尺寸为24 mm × 12.8 mm(KJ × JH).焊丝的半径为0.6 mm，铜嘴极的尺寸设置为1 mm，焊丝端部距离母材上表面的距离为5 mm，母材的厚度为6 mm，铜嘴出口半径1.5 mm，以上参数均是根据物理试验情况设定.全计算域采用四边形网格离散，由于焊丝下方、铜嘴下方以及气液交界面的区域传热传质行为较为剧烈，因此采用较小的网格尺寸0.1 mm × 0.1 mm，为了提高计算效率，其它部分的网格尺寸相对粗化.

图1 旁路等离子-MIG 复合工艺原理及计算域Fig.1 Principle and calculation domain of plasma-MIG hybrid technology.(a) principle; (b) hybrid arc and molten pool;(c) calculation domain

1.2 基本假设

为了合理地计算、分析和证明旁路等离子-MIG 复合电弧和耦合熔池的复杂流体动力学行为，建立了如下重要假设.

(1) 设置的模型是轴向对称的.

(2) 流体是不可压缩的牛顿流体，流体流动为层流.

(3) 在计算过程中，旁路等离子体-MIG 复合电弧直接作用于熔池的上表面，这里不考虑熔池变形.

(4) 气体为纯氩气，其物理特性包括密度、比热、热导率等随温度变化，高温等离子体属于局部热动力平衡(LTE).

(5) 不考虑金属液滴与电弧和熔池的相互作用.

1.3 控制方程

计算涉及的控制方程主要包括连续方程、动量守恒、能量守恒方程.二维坐标系下,连续方程为[11]

式中：ρ为密度；t为时间；u和v分别为径向r和轴向a的速度.

径向和轴向的动量守恒方程为[11]

式中：p为压力；μ为动态粘度；Sr和Sa为径向和轴向源项；Bθ为磁场强度，可由麦克斯韦方程得到；ja和jr为轴向和径向电流密度.

能量守恒方程为[11]

式中：H为焓；k为导热系数；T为温度；SE为能量源项.

1.4 界面追踪

焓-孔技术[11]是用来处理固相向液相转化过程中的热量变化，可以表示为

式中：h为金属的焓；Cp为比热；fl为液相分数；L为潜热；Ts和Tl分别为固相和液相温度.

1.5 源项

源项包括动量和能量源项，动量源项包括热浮力、马兰戈尼力，分别表示为[11]

式中：β为热膨胀系数；A为糊状区常数；τ为马兰戈尼力；γ为表面张力系数；n为局部表面的切线方向.电弧区的能量源项为

式中：σe为电导率；Hlos为热损失.另外界面处的能量源项包括阴极中的电子激活、能量传递以及其表面的热损失被添加，即

式中：jc为阴极表面的电流密度；ϕw为阴极的功函数；qc为传导能量；σ为Stefan-Boltzmann 常数；ε为辐射率.

1.6 边界条件

阳极的边界(AB，BC)、等离子铜嘴的边界(DE，EF，FG)以及两个气体速度入口(CD，GH)的温度被设定为1 000 K，KA 为轴向对称轴，HI 为压力-出口边界，在基底区域的边界(IJ 和JK)有两个热交换过程，即对流和辐射，可以表示为[12]

式中：hconv为对流传热系数；Tamb为环境温度.

1.7 材料模型

试验使用的母材是典型船用5083 铝合金，其热物理性能[13]如表1 所示.气相是纯氩气且不考虑空气的影响[14]，其热物理性能见文献[15].

表1 5 083 铝合金热物理性能Table 1 5 083 aluminum alloy property parameters

2 计算结果与分析

2.1 复合电弧与耦合熔池演变过程

计算过程采用的工艺参数如表2 所示.计算结果如图2 所示，其中图2a 为施加旁路电流时的复合电弧和耦合熔池的演变过程，图2b 为不施加旁路电流时的电弧和耦合熔池的演变过程.可以看出t=0.2 s 时，这两者的电弧和熔池都处于长大状态，其中图2a 中的母材仍处于尚未熔化状态；t=0.5 s 时，二者的电弧都已经趋于稳定，但熔池的尺寸随着热通量的增加还在继续增长;随着时间的增加，t=0.7 s 时，熔池内部的液态金属流动行为已经相对剧烈；t=1.2 s 时，图2b 的母材已经熔透，而图2a 的熔池形态还处于稳定的状态.

表2 计算过程工艺参数Table 2 Process parameters

图2 旁路等离子-MIG 复合电弧与熔池演变过程Fig.2 Plasma-MIG hybrid arc and molten pool evolution.(a) with bypass-current; (b) without bypass-current

从计算结果来看，施加旁路电流时，母材的热输入降低效果是明显的，此外还可以看出施加旁路电流后，电弧的高温等离子体形态产生了明显的变化，这种变化主要体现在电弧下方出现了明显的收缩趋势，而电弧加热面积的减少也是母材热输入降低的主要原因之一.由于进入电弧区的气相带电粒子部分向阴极扩散，导致图2a 的电弧自身的热效应也有所降低，这可以从图3 的电弧最高温度曲线得出.相比图2b 所示结果，图2a 的电弧最高温度从初始时刻就较低，这个最大温差接近1 000 K，且随着时间变化一直稳定存在，如图3 复合电弧温度曲线所示.差异的电弧行为最终导致了熔池内部液态金属的对流换热行为也产生了差异.

图3 电弧最高温度与熔池最大流速随时间变化趋势Fig.3 Trend of maximum arc temperature and maximum melt pool flow velocity over time

整体来看，熔池内部液态金属的流速均呈现先增加后稳定的变化趋势，且施加旁路电流时熔池内部液态金属的流速更低，如图3 中最大流速曲线所示.当施加旁路分流时，从初始时刻到t=0.2 s，母材几乎不熔化，而不施加旁路分流时，母材已经熔化，且液相金属的最大流速达到0.42 m/s；t=0.5 s，这两者的液态金属流速基本达到最大值，以后时刻增长幅度都很小；但是随着热输入增加，母材的活化金属量继续增加，熔池的面积继续增大，不施加旁路电流时的母材金属活化速度要远远高于施加旁路电流时，这可以通过图2 中母材熔池面积的增加趋势看出，最终t=1.2 s，不施加旁路电流时，母材出现熔透现象，而施加旁路电流时，熔池仍处于稳定状态.

2.2 电弧—熔池界面热通量分析

图4 为电弧与熔池的交界面处的有效热通量对比情况.由于电弧的本质是高温等离子体，因而在焊接加工过程中电弧本身要向周围低温环境散失能量，而真正由电弧传递给母材的能量是焊接加工过程中的有效能量，即除去热损失而真正形成熔池的能量，此处通过有效热通量的方式对有无旁路电流施加时有效能量区别进行对比分析.从计算结果来看，当无旁路电流施加时，电弧与熔池界面处的有效热通量整体更高，但是当施加旁路电流时，界面中心的能量峰值更高，这是因为旁路电流施加后，电弧的下部分产生压缩态，因而在中心部分的能量密度更高，但是由于部分带电粒子没有向母材运动，趋向母材运动的能量粒子水平更低，因而随着径向距离的增加，能量的分布快速下降，即有效加热面积减少，整体的能量水平也就更低.

图4 复合电弧与耦合熔池界面处的有效热通量对比Fig.4 Comparison of the effective heat flux at the interface between the hybrid arc and the molten pool

热量由电弧通过气液交界面向母材传递，经过一段时间的能量演变和积累逐渐形成熔池.图5 为t=1.2 s 时熔池内部液态金属径向和轴向的流速以及它们对熔池尺寸的影响对比.在轴向上，液态金属在温度梯度的作用下由表面向下运动，在一定范围内形成回流，从计算结果可以看出，当旁路电流施加时，轴向的液态金属最大流动速度明显降低，最终也导致了熔深尺寸的下降；在径向上，液态金属的回流驱动力还包括马兰戈尼力，从计算结果可以看出径向上的液态金属回流运动最大速度整体水平较高；当施加旁路电流时，径向上的液态金属流动速度相也有所降低，最终也导致熔池的熔宽尺寸有所下降.整体来看，旁路电流的施加有减小熔池尺寸的作用，在熔深方向的减小程度更明显，这是因为旁路电流间接导致电弧中心附近能量传递水平下降剧烈，例如在距离中心1.5 mm 左右，有效热通量已经下降了约60 × 106W/m2，因此母材积累的热量水平也就下降，液化的母材体积更少，以上计算结果和分析说明这种新型工艺在电弧能量精确调控和减少母材热输入方面是有优势的.

图5 熔池内部液态金属流速对熔池尺寸的影响Fig.5 Influence of the liquid metal flow velocity on sizes of the molten pool

2.3 电弧与熔池电磁效应分析

图6 为t=1.2 s 时整个计算模型内，包括电弧区和母材区的电流密度矢量分布.从计算结果来看，当施加旁路电流时，在焊丝端和分流极端都有较高的电流密度分布，且矢量方向相反，这与试验设置是吻合的.

图6 计算域内电流密度矢量分布情况对比Fig.6 Comparison of the current density vector distribution in the calculation domain.(a) with bypass-current;(b) without bypass-current

在电弧区，无论是否施加旁路电流，最大电流密度区均出现在焊丝极下方的电弧中心区，且方向均是由焊丝极指向母材极，但是最高电流密度的数值有所不同.当施加旁路电流时，电弧区的最大电流密度为1.32 × 108A/m2，熔池区的最大电流密度为1.3 × 107A/m2；不施加旁路电流时，电弧区的最大电流密度为1.71 × 108A/m2，熔池区的最大电流密度为2.25 × 107A/m2.旁路电流施加前后电弧区最大电流密度相差接近1.3 倍，而母材区相差约1.7 倍，这是由于部分带电粒子在旁路电流施加前后分别向分流电极和向母材极运动的结果，而电流密度的不同又会对磁场强度及分布产生影响.

图7 为t=1.2 s 时，计算域内电磁矢量的分布和数值对比情况.从计算结果来看，无论是否施加旁路电流，最大磁力的位置和方向并没有改变，但是其在电弧区的数值因为旁路电流的施加从7.1 ×106N/m3下降到5.1 × 106N/m3；而其在熔池区的数值从8.4 × 106N/m3下降到5.6 × 106N/m3，下降的趋势基本接近，这说明旁路电流的施加不会改变最大磁场强度的方向，但会降低电弧和熔池内部的电流密度的数值，进而降低电磁力的数值，即旁路电流有弱化磁场作用.

图7 计算域内电磁矢量分布情况对比Fig.7 Comparison of electromagnetic vector distribution in the calculation domain.(a) with bypasscurrent;(b) without bypass-current

2.4 试验验证

图8 为施加旁路分流时的等离子-MIG 复合电弧和熔池在t=1.2 s 时的计算与试验结果对比，其中类似的双侧分流机制下的高温等离子体图像已经由兰州理工大学石玗教授团队[15]成功拍摄(图8 左上所示).从对比结果看出电弧下方的等离子体具有更高的温度，而两侧的分流弧的等离子体温度较低，图像中其亮度也较弱.熔池的试验结果如图8 左下所示.总体来看，模拟结果与试验结果在一定程度上能够匹配，但是也存在一些误差：其一是电弧根部的发散；其二是熔宽尺寸误差约0.5 mm，造成这些误差的原因是由于没有考虑熔滴的影响，其次数据处理过程也会产生一些误差.

图8 模拟与试验结果对比[15]Fig.8 Comparison between simulation and experiment results

3 结论

(1) 其它试验条件相同，施加旁路电流后，电弧有压缩趋势，此外复合电弧和母材交界面处的有效热通量整体减少.

(2) 施加旁路电流后，熔池内部液态金属在轴向和径向上的回流速度有所降低，最终导致熔池面积减小.

(3) 旁路电流的施加使电弧和熔池的电流密度下降约23%和42%，因此导致等离子-MIG 复合电弧和耦合熔池的最大电磁力数值减小约28%和33%，但其方向没有改变.