-
等离子弧切割作为一种热切割方式,是以电离击穿的离子气为工作介质,将熔融金属吹除进而形成割缝的一种切割技术,激发的电弧收缩质量好,穿透能力强 [1-2]. 同等条件下,相对于气割与电弧切割,等离子弧切割不仅效率高,而且在切割质量和切割精度上也优于前两种切割方式. 此外切割材料的种类可以覆盖绝大部分的金属[3-4]. 等离子切割电弧与TIG(Tungsten inert gas)电弧本质相同,TIG电弧被认为是等离子态,两者本质上是阴极与阳极之间的一种剧烈的气体电离行为. 根据等离子体温度的分类,TIG 焊接电弧属于低温等离子体,而等离子切割电弧属于高温等离子体,因此等离子弧被认为是焊接电弧的升级[5-6]. Matus等人[7]使用等离子弧切割技术在0.3 MPa高压环境下完成了油管和套管的切割,提高了作业效率,降低了人力成本. 陈林柯等人[8]采用高频高压引弧技术,实现了高压环境(0.4 MPa)下等离子弧的稳定切割,并通过试验得出了随环境压力的升高,电弧长度逐渐变短的结论. 谷孝满等人[9]通过对高压环境(0.1 ~ 0.5 MPa)下脉冲MIG焊接电弧的研究,得出弧长与压力的变化量不成正比的结论.
电弧电离行为的理论研究主要是应用数学建模的方法进行数值模拟,数学模型主要有磁流体动力学模型、流体模型、粒子模型等. 胡辉等人[10]依据磁流体动力学理论,建立了等离子体的数学模型,对脉冲放电的等离子弧温度场以及电流密度分布进行了数值计算. Baimbetov等人[11]引出等离子体的粘性以及辐射能量等概念,并提出热导率、电导率等物性参数的相关计算方法,由此建立了等离子体流体模型. Carman等人[12]通过等离子体的粒子化构思,构建了碰撞动力学模型,得到了氩气在电离过程中电子在鞘层区的分布状态.
文中研究了非高频引弧技术下的等离子弧,建立并优化了电弧磁流体模型,通过COMSOL Multiphysics软件实现了数值模拟;并搭建了实验平台,实现了环境压力为0.1 ~ 0.7 MPa的稳定起弧,最终通过高压梯度试验验证了电弧形态及能量变化的规律.
-
磁流体动力学(magnetohydrodynamics,MHD)模型本质是一种求解基本流体方程以及耦合的电磁场方程的数学模型. 切割电弧等离子体不同于一般流体,需要考虑其中存在的磁场和电场作用. MHD模型符合电弧流体的基本参数要求,通过构建合理的电弧几何模型,结合电弧的控制方程,进而达到对电弧流体研究的目的. MHD方程组是由一组非线性偏微分方程组成,其中包含基本流体方程、Maxwell电磁方程以及广义欧姆定律[13-14].
-
等离子切割电弧是一种压缩电弧,具有电离程度高、挺度好、放电持久、电弧形态稳定的优点. 为便于研究高压环境等离子切割电弧的特性,在不影响仿真目标的前提下,假定切割电弧流体域呈轴对称分布,进而达到电弧模型简化的目的. 等离子弧割炬模型如图1所示.
为便于进一步分析高压环境下等离子切割电弧的电离行为,对等离子切割电弧的模型做如下基本假设[15-17].
假设1:等离子切割电弧是由逆变电源产生的直流电弧,电弧是稳定的二维轴对称结构,忽略电弧流动的粘性,其电场为有势场;
假设2:认为等离子切割电弧是连续的流体,其流动与传热特性可用N-S方程来描述. 等离子体电弧的热力学参数以及电导率是仅为温度的函数,不考虑等离子体辐射的自吸收效应,认为等离子电弧是光学薄的;
假设3:等离子切割电弧处于化学平衡状态,满足准电中性条件,电弧整体对外不显电性;
假设4:等离子切割电弧的阳极处于不活跃状态,不考虑金属蒸汽对电弧的影响;
假设5:等离子切割电弧为理想气体,不考虑重力作用.
-
等离子切割电弧不仅具有一般流体的流动特性,而且因其受到电磁场的作用,其流动特性更加复杂. MHD方程组的合理简化与求解,有利于问题的简化,电弧数学模型的构建依赖于MHD方程组的求解[15-17].
(1)状态方程,认为等离子切割电弧符合理想气体状态方程,每种组分的分压
$ {p}_{j} $ ,数密度$ {n}_{j} $ ,温度$ {T}_{j} $ 满足式(1),即$$ p=\sum {p}_{j}=\sum ({n}_{j}{\mathrm{k}}_{\mathrm{B}}{T}_{j})=n{\mathrm{k}}_{\mathrm{B}}T $$ (1) 式中:
$ {\mathrm{k}}_{\mathrm{B}} $ 是玻尔兹曼常数($ {\mathrm{k}}_{\mathrm{B}}=1.381\times 1{0}^{-23}J/K $ );$ p $ 为气体总压;$ n $ 是总粒子数密度;$ T $ 是气体总温度.(2)质量守恒方程可表示为
$$ \frac{\partial \rho }{\partial t} + \nabla \cdot \left(\rho \nu \right)=0 $$ (2) 式中:
$ \rho $ 是气体混合物的质量密度;$ \nu $ 为气体混合物速度.(3)动量守恒方程可表示为
$$ \frac{\partial \left({\rho }_{j}{\nu }_{j}\right)}{\partial t} + \nabla \cdot \left({\rho }_{j}\nu {\nu }_{j}\right) =\nabla \cdot {P}_{j} + {\rho }_{j}{F}_{j}-{M}_{j} $$ (3) 式中:
$ {P}_{j} $ 是$ j $ 组分的应力张量;$ {M}_{j} $ 为$ j $ 组分动量损失率;$ {F}_{j} $ 是作用于$ j $ 粒子的体积力;$ {\nu }_{j} $ 为速度分量;$ {\rho }_{j} $ 是$ j $ 组分的质量密度 .(4)能量守恒方程可表示为
$$ \begin{split} & \rho \frac{D}{{D}_{t}}\left(\varepsilon + {v}^{2}/2\right)=-\nabla \cdot q + \nabla \cdot \left(p\cdot v\right) + v\cdot {\sum }_{j}{\rho }_{j}{F}_{j} +\\&\qquad {\sum }_{j}{\rho }_{j}{F}_{j}\cdot {U}_{j} + {S}_{h} \end{split}$$ (4) 式中:
$\varepsilon$ 是单位质量体积的内能,单位$ J/kg $ ;其中$ {S}_{h} $ 为体热源;$ {U}_{j} $ 为扩散速度;$ q $ 为热流量.(5)广义欧姆定律可表示为
$$ \overrightarrow{j}=\sigma (E + V\times B) $$ (5) 式中:
$ \overrightarrow{j} $ 为电流密度;$ \sigma $ 为电导率;$ E $ 为电场强度;$ V $ 为速度矢量;$ B $ 为磁场强度. -
电弧模型的求解需要不同的物理边界条件和初始值. 边界条件是电弧流体方程具有解的基本条件,描述的是待求量和其相关导数在流体边界上的取值情况,此外控制方程在计算域中的各个边界条件也均有不同. 边界条件如下表1所示.
表 1 求解域边界条件
Table 1. boundary conditions of solution domain
边界 电场条件 磁场条件 传热条件T1/K 流体条件 AB 绝缘 绝缘 300 壁 BC 绝缘 绝缘 300 壁 CD 绝缘 绝缘 热通量 出口 DE 接地 绝缘 热通量 — EF 接地 绝缘 热通量 — FG 接地 绝缘 热通量 — GH 绝缘 绝缘 热通量 出口 GI — — 热源 壁 ID — — 热源 壁 HJ 对称 对称 对称 对称 JK — — 热源 壁 KM — — 热源 壁 JL 对称 对称 对称 — LM 电流密度 绝缘 300 — AM 绝缘 绝缘 300 入口 表1中电流密度可根据电流值与电极截面积计算得出;热通量选用对流热通量,自定义传热系数和外部温度值;热源选择自定义广义源;流体场中边界壁类型均为无滑移;入口为速度边界,出口为压力边界,两者可根据压力差与环境压力等具体条件计算得出.
-
等离子切割电弧作为特殊流体,其电弧由压缩空气电离产生,其物性参数不同于空气,文中采用空气等离子体的热物性参数作为电弧的基本参量. 仿真过程中高质量的网格剖分有利于提高流体动力模型的收敛性和计算效率,因此选择用户自定义网格,其中空气域选择极细化网格剖分;电极域和工件域选择细化网格剖分;边界JK,KM,DI和GI选择超细化网格剖分;空气域全部边界选择角细化网格剖分. 网格剖分如图2所示.
-
采用基于环境压力的控制变量法进行仿真. 对比0.1 ~ 0.7 MPa环境压力下的电弧形态有利于研究电弧电离行为. 由温度场分布图3和图4可以看出,0.1 ~ 0.7 MPa的电弧在整体温度分布上呈现逐渐收缩趋势. 电弧整体呈射流状,经喷嘴射出电弧呈发散状,几何平面形态呈梯形. 此外,电弧中心温度较高,径向与轴向梯度较大.
图 3 0.1 ~ 0.4 MPa环境压力下电弧温度场
Figure 3. Arc temperature field under 0.1 ~ 0.4 MPa ambient pressure. (a) 0.1 MPa; (b) 0.2 MPa; (c) 0.3 MPa; (d) 0.4 MPa
图 4 0.5 ~ 0.7 MPa环境压力下电弧温度场
Figure 4. Arc temperature field under 0.5 ~ 0.7 MPa ambient pressure. (a) 0.5 MPa; (b) 0.6 MPa; (c) 0.7 MPa
电弧的速度场分布如图5和图6所示,切割电弧中心处速度最大,径向梯度较大. 随着环境压力的升高,速度逐渐收敛,整体速度下降;速度收敛,进而导致了电弧的动量降低,吹除能力下降,因此导致工件上的挂渣增多. 最终可知环境压力越大,对其熔融金属流动的阻碍作用力越强.
图 5 0.1 ~ 0.4 MPa环境压力下电弧速度场
Figure 5. Arc velocity field under 0.1 ~ 0.4 MPa ambient pressure. (a) 0.1 MPa; (b) 0.2 MPa; (c) 0.3 MPa; (d) 0.4 MPa
图 6 0.5 ~ 0.7 MPa环境压力下电弧速度场
Figure 6. Arc velocity field under 0.5 ~ 0.7 MPa ambient pressure. (a) 0.5 MPa; (b) 0.6 MPa; (c) 0.7 MPa
等离子切割电弧主要依靠其高温熔化金属,进而吹除形成切口,最终达到切割工件的目的. 因此,分析不同环境压力下切割试件的温度场分布规律,有助于研究工件切口宽度等切割质量指标的变化情况,从而进一步研究电弧电离行为.
工件与电弧的接触面上温度较高,温度由高到低,由近及远整体呈现扩散式分布. 由图7和图8可知,随着环境压力的提高,工件温度最大值超过5 000 K,远超过了Q235碳钢的熔点. 此外结合图3和图4可知,0.1 ~ 0.7 MPa的温度范围逐渐增大,这就导致了在切割试验过程中,上部金属熔化较多,切口增大;而0.1 ~ 0.7 MPa环境压力下的工件底部温度范围逐渐减小,但整体差异性较小. 由此可见,在0.1 ~ 0.7 MPa环境压力下,底部熔融金属逐渐减小. 进而导致在切割试验过程中,底部金属熔化较少,底部切口较小. 因此,工件的整体缺口一般呈现凸状. 随着环境压力的增大,倾斜角增大,平面度增大,切割质量下降,电弧整体的切割能力下降. 印证了随着环境压力的增大,工件切割质量变差的结果. 同时也指出了高压环境下工艺参数的重要性.
-
高压环境下的等离子切割试验系统主要由空气等离子切割电源、三轴运动平台、焊接相机系统、立式高压舱与气路系统等组成. 为实现对高压环境下的引弧切割进行控制,同时保障实验操作的安全性,采用管线穿舱的连接形式达到舱外控制的目的. 以高压焊接试验舱为基础,采用非高频引弧技术,搭建试验平台如图9所示.
在不同的环境压力下,设置相同的切割工艺参数:切割速度1.2 m/min、切割电流120 A、送气压力1.5 MPa、喷嘴高度7 mm. 调整环境压力,通过焊接相机拍摄得到高压环境下的电弧图像,如图10所示.
在不同环境压力下电弧整体形态差异较小,但随着环境压力的升高,弧柱有收缩的趋势,亮度有所增加. 从火花飞溅程度可知,电弧稳定性下降. 如图10所示,随着环境压力的升高,电弧挺度下降,切割试件表面处的电弧形态呈收缩状态,电弧整体有收缩趋势但并不明显,弧柱扭曲程度增大,电弧稳定性变差. 随着环境压力的提高,受环境压力的挤压作用,能量相对集中,此时的弧柱等离子体密度较大,亮度较高. 根据光电效应理论可知,当大量的电子同时向低能级跃迁时,会激发大量光子. 因此距离弧柱中心越近,电子的密度越大,同时电子跃迁的程度越高,导致电弧的电离程度就越大,能量越大,电弧的亮度越高.
通过分析切割质量指标的变化规律,可进一步验证上述规律的正确性. 采用基于环境压力的控制变量法进行试验,相关评价指标的数值如表2所示,试验工件割面如图11所示.
表 2 切割质量参数
Table 2. Cutting quality parameters
环境压力P/MPa 平面度F/mm 切口宽度W/mm 后拖量N/mm 0.1 0.7 3.5 1.8 0.2 0.8 3.5 3.4 0.3 1.1 3.6 4.5 0.4 1.4 3.7 3.9 0.5 1.6 3.8 5.5 0.6 1.7 3.8 5.2 0.7 2.0 3.9 4.8 由图12可知切割质量与电弧行为变化的关系:随着环境压力的增大,电弧所受阻力增大,气体密度越大,电弧与周围气体之间作用力越大,导致两者间的热交换速率增加,进而能量输出率增加,最终导致切口处板材熔化量增多,出现切口宽度增加的情况,但整体变化趋势较平缓;随着环境压力的升高,电弧受气压的向上阻力增大,试件反面电弧受压收缩,电弧挺度变差,导致平面度增大;随着环境压力的升高,电弧受气流影响较大,后拖量处于波动状态,故呈增加趋势.
-
(1)通过COMSOL Multiphysics构建电弧的几何模型,利用平衡放电等离子模型,实现了包含电场、磁场、传热和流体的多物理场耦合仿真. 最终通过稳态求解完成了等离子切割电弧的数值模拟. 通过分析不同环境压力下电弧的温度场、速度场以及工件的温度场可知:随着环境压力的升高,电弧在温度和速度分布上都呈现收缩的态势;随着环境压力的升高,工件的温度分布表明,工件融化的温度范围具有增大的趋势.
(2)采用非高频引弧技术,基于相同工艺参数,实现0.1 ~ 0.7 MPa的引弧试验. 通过高压梯度试验可知:随着环境压力的提高,切口宽度、平面度和后拖量均有增加,间接验证了文中关于等离子切割电弧电离行为的分析.
Study on behavior and process of plasma arc in high pressure environment
-
摘要: 等离子弧切割技术因其高效稳定的工艺优势被广泛应用于工业领域. 文中以空气等离子弧为研究对象,通过COMSOL Multiphysics软件建立了喷嘴结构的二维轴对称有限元数学模型,并对电弧磁流体模型进行了优化. 基于磁流体动力学和电弧等离子体理论,选用等离子平衡放电多物理场接口,并确立了空气等离子电弧模型的控制方程和边界条件,实现对电弧模型的编译求解. 仿真结果表明,在引弧电流一致的条件下,随着环境压力的增加,电弧在温度分布和速度分布上均呈现收缩的态势. 基于3 MPa 高压焊接试验舱,搭建了高压环境等离子弧切割实验系统,通过对气路和非高频引弧电路的优化设计,实现了环境压力为0.1 ~ 0.7 MPa的稳定起弧. 并基于此开展了高压梯度下的等离子弧切割实验,并结合切割质量指标研究了环境压力对等离子弧电离行为的影响.Abstract: Plasma arc cutting technology is widely used in industrial fields because of its efficient and stable process advantages. Taking the air plasma arc as the research object, this paper establishes a two-dimensional axisymmetric finite element mathematical model of the nozzle structure through COMSOL multiphysics software, and optimizes the arc magnetic fluid model. Based on magnetohydrodynamics and arc plasma theory, the multi physical field interface of plasma equilibrium discharge is selected, and the control equations and boundary conditions of air plasma arc model are established to realize the compilation and solution of the arc model. The simulation results show that under the condition of consistent arc starting current, the arc shrinks in temperature distribution and velocity distribution with the increase of ambient pressure. Based on the 3Mpa high-pressure welding test chamber, an experimental system of plasma arc cutting in high-pressure environment is built. Through the optimization design of gas circuit and non high-frequency arc starting circuit, the stable arc starting with environmental pressure of 0.1 ~ 0.7 MPa is realized. Based on this, plasma arc cutting experiments under high pressure gradient were carried out, and the influence of environmental pressure on plasma arc ionization behavior was studied combined with cutting quality index.
-
表 1 求解域边界条件
Table 1. boundary conditions of solution domain
边界 电场条件 磁场条件 传热条件T1/K 流体条件 AB 绝缘 绝缘 300 壁 BC 绝缘 绝缘 300 壁 CD 绝缘 绝缘 热通量 出口 DE 接地 绝缘 热通量 — EF 接地 绝缘 热通量 — FG 接地 绝缘 热通量 — GH 绝缘 绝缘 热通量 出口 GI — — 热源 壁 ID — — 热源 壁 HJ 对称 对称 对称 对称 JK — — 热源 壁 KM — — 热源 壁 JL 对称 对称 对称 — LM 电流密度 绝缘 300 — AM 绝缘 绝缘 300 入口 表 2 切割质量参数
Table 2. Cutting quality parameters
环境压力P/MPa 平面度F/mm 切口宽度W/mm 后拖量N/mm 0.1 0.7 3.5 1.8 0.2 0.8 3.5 3.4 0.3 1.1 3.6 4.5 0.4 1.4 3.7 3.9 0.5 1.6 3.8 5.5 0.6 1.7 3.8 5.2 0.7 2.0 3.9 4.8 -
[1] 胡义勇, 吴波, 顾先凯, 等. 浅谈等离子弧切割在管道开孔中应用[J]. 中国石油和化工标准与质量, 2021, 41(15): 118 − 119. Hu Yiyong, Wu Bo, Gu xiankai, et al. Discussion on the application of plasma arc cutting in pipe opening[J]. China Petroleum and chemical industry standards and quality, 2021, 41(15): 118 − 119. [2] 胡义勇, 李夏喜, 李松, 等. 带压燃气环境下等离子弧切割技术可行性研究[J]. 煤气与热力, 2019, 39(5): 37 − 39. Hu Yiyong, Li Xiaxi, Li Song, et al. Feasibility study on plasma arc cutting technology in pressurized gas environment[J]. Gas and heat, 2019, 39(5): 37 − 39. [3] 陈月峰, 严得忠, 程定富, 等. 国内数控等离子弧切割技术的发展与应用[J]. 焊接技术, 2017, 46(8): 1 − 4. Chen Yuefeng, Yan Dezhong, Cheng Dingfu, et al. Development and application of CNC plasma arc cutting technology in China[J]. Welding technology, 2017, 46(8): 1 − 4. [4] 陈传杰, 孔维宾, 方忠庆, 等. 基于电磁模型的大气压表面波等离子体数值模拟[J]. 电子科技大学学报, 2022, 51(2): 305 − 313. Chen Chuanjie, Kong Weibin, Fang Zhongqing, et al. Numerical simulation of atmospheric pressure surface wave plasma based on electromagnetic model[J]. Journal of University of Electronic Science and technology, 2022, 51(2): 305 − 313. [5] 张达, 叶凯, 唐政刚, 等. 等离子体冶金的现状与发展[J]. 中国有色金属学报, 2021, 31(7): 1907 − 1921. Zhang Da, Ye Kai, Tang Zhenggang, et al. Current situation and development of plasma metallurgy[J]. Chinese Journal of nonferrous metals, 2021, 31(7): 1907 − 1921. [6] 刘吉宇, 张帆, 陈阳, 等. 低温等离子体辅助加工综述[J]. 航空学报, 2021, 42(10): 74 − 86. Liu Jiyu, Zhang Fan, Chen Yang, et al. Overview of low temperature plasma assisted machining[J]. Acta Aeronautica Sinica, 2021, 42(10): 74 − 86. [7] Matus Gajdos, Tomas Kristofic, Slavomir Jankovic, et al. Use of plasma-based tool for plug and abandonment[C]. Aberdeen, Scotland, UK: SPE Offshore Europe Conference and Exhibition, 2015. [8] 陈林柯. 环境压力对等离子切割电弧的影响研究[D]. 北京石油化工学院, 2019. Chen Linke. Study on the influence of ambient pressure on plasma cutting arc [D]. Beijing, Beijing Institute of petrochemical technology, 2019. [9] 谷孝满, 黄松涛, 焦向东, 等. 高气压环境下脉冲MIG焊频率优化改善电弧稳定性的研究[J]. 电焊机, 2015, 45(3): 69 − 74. Gu Xiaoman, Huang Songtao, Jiao Xiangdong, et al. Research on improving arc stability by optimizing the frequency of pulsed MIG welding in high pressure environment[J]. Electric welding machine, 2015, 45(3): 69 − 74. [10] 胡辉, 杨旗, 包斌, 等. 基于ANSYS的空气电弧放电等离子体温度数值模拟[J]. 电工电能新技术, 2009, 28(2): 20 − 23. Hu Hui, Yang Qi, Bao bin, et al. Numerical simulation of plasma temperature in air arc discharge based on ANSYS[J]. New technology of electrical energy, 2009, 28(2): 20 − 23. [11] Baimbetov F B, Giniyatova S G. Interparticle interaction and transport processes in dense semiclassical plasmas[J]. Physics of plasmas, 2005, 12(4): 1 − 4. [12] Herrebout D, Bogaerts A, Yan M, et al. Modeling of a capacitively coupled radio-frequency methane plasma: Comparison between a one-dimensional and two-dimensional fluid model[J]. Journal of Applied Physics, 2002, 92(5): 2290 − 2295. doi: 10.1063/1.1500789 [13] Annemie Bogaerts, Luís L Alves. Special issue on numerical modelling of low-temperature plasmas for various applications — part II: Research papers on numerical modelling for various plasma applications[J]. Plasma Processes and Polymers, 2017, 14(4-5). [14] 王路明, 高辉, 周灿丰, 等. 高压环境等离子切割电弧数值模拟[J]. 焊接, 2021(5): 15 − 19. Wang Luming, Gao Hui, Zhou Canfeng, et al. Numerical simulation of plasma cutting arc in high pressure environment[J]. Welding & Joining, 2021(5): 15 − 19. [15] Rong M, Li M, Wu Y, et al. 3-D MHD modeling ofinternal fault arc in a closed container[J]. IEEE Transactions on Power Delivery, 2017, 32(3): 1220 − 1227. doi: 10.1109/TPWRD.2014.2375891 [16] 向凌峰. 基于磁流体力学模型的稳态真空电弧建模与仿真[D]. 华中科技大学, 2019. Xiang Lingfeng. Modeling and simulation of steady-state vacuum arc based on MHD model [D]. Wuhan, Huazhong University of science and technology, 2019. [17] 李国明. 真空断路器开断过程中电弧等离子体的特性研究[D]. 华北电力大学, 2019. Li Guoming. Study on the characteristics of arc plasma during the breaking process of vacuum circuit breaker [D]. Beijing. North China Electric Power University, 2019. [18] 李天庆, 陈璐, 张宇, 等. 气流再压缩等离子弧焊接电弧行为[J]. 焊接学报, 2020, 41(5): 50 − 55. Li Tianqing, Chen Lu, Zhang Yu, et al. Arc behavior of gas flow recompression plasma arc welding[J]. Transactions of The China Welding Institution, 2020, 41(5): 50 − 55. [19] 李正农, 李炜, 张学文. 应用等离子体激励抑制风力机叶片噪声的数值仿真研究[J]. 太阳能学报, 2021, 42(12): 220 − 229. Li Zhengnong, Li Wei, Zhang Xuewen. Numerical simulation study on suppression of wind turbine blade noise by plasma excitation[J]. Journal of solar energy, 2021, 42(12): 220 − 229. [20] 李岩, 王领, 张冀翔, 等. 基于电弧-熔池耦合的小孔型等离子弧焊接简化模型[J]. 机械工程学报, 2021, 57(18): 144 − 152. Li Yan, Wang Ling, Zhang Jixiang, et al. A simplified model of small hole plasma arc welding based on arc pool coupling[J]. Journal of mechanical engineering, 2021, 57(18): 144 − 152. -