Study on dynamic behavior of plasma in deep penetration laser welding of TC4 titanium alloy based on triple-probe detection
-
摘要:
在激光深熔焊中，激光诱导等离子体的产生具有很强的随机性和时变性.利用新型无源三探针电信号采集装置，采集高时间分辨率的等离子体电信号，获得大样本的等离子体温度等热力学信息，采用无监督学习方法中的K-均值聚类对复杂的等离子体动态行为聚类分析.结果表明，在文中试验条件下，等离子体的喷发频率主要集中在
1240 -2200 Hz范围，等离子体喷发的温度变化率主要分布在5000 -30000 K/ms的范围内；等离子体的摆动频率主要在1050 -1410 Hz范围，等离子体摆动的温度变化率主要分布在200-2000 K/ms的范围内，喷发和摆动的温度变化率差异明显.Abstract:In deep penetration laser welding, the generation of the laser-induced plasma exhibits significant randomness and temporal variability. A novel passive triple-probe electrical signal acquisition device was utilized to collect the plasma electrical signals with high temporal resolution, thereby obtaining the large sample of thermodynamic information, including plasma temperature. K-means clustering, an unsupervised learning method, was used to conduct the clustering analysis of the complex dynamic behaviors of the plasma. The results indicated that under the experimental conditions of present work, the frequency of the plasma eruption mainly fell within the range of
1240 -2200 Hz, and the temperature change rate was primarily distributed between5000 -30000 K/ms. Additionally, the sway frequency of the plasma was mainly observed in the range of1050 -1410 Hz, and the temperature change rate concentrated between 200-2000 K/ms. Consequently, there was a significant difference in the temperature change rates between plasma eruption and sway. -
0. 序言
残余应力和变形是焊接过程中常见的问题，会对焊接工件的质量和性能产生重要影响，因此，预测焊接残余应力和变形已经成为焊接领域中的重要研究内容. 焊接工艺参数，焊后残余应力及变形之间是一种非常复杂的非线性映射关系，在厚板焊接中对焊接残余应力与变形影响较大的填充层焊接工艺参数主要有焊接电流、电弧电压、焊接速度以及层间温度[1-3]. 焊接残余应力与变形的预测是高度非线性问题，很难用简单的数学模型来描述，需要通过各种方法和经验对其中的权重和因子进行计算.
人工神经网络可实现复杂非线性系统的建模和预测，具有强大的联想能力、容错能力以及自适应学习能力，近年来在焊接质量预测中得到了广泛关注[4-7]，其中BP神经网络是人工神经网络一种常用的深度学习模型，可以在没有大量数据的情况下进行预测，BP神经网络具有较强的自适应学习能力和非线性拟合能力，已有诸多学者将其应用于非线性领域的预测中，例如Shanavas等人[8] 通过建立反向传播神经网络模型，来预测铝合金搅拌摩擦焊焊缝的抗拉强度和伸长率，模型预测误差小于5%；马思群等人[9] 采用BP神经网络的方法对铝合金车体的焊接残余应力进行预测，预测误差在8%以内.
传统的BP神经网络在训练过程中容易出现梯度消失或梯度爆炸等问题，导致模型训练困难或预测效果不佳，遗传算法是一种优化算法，可用于解决复杂问题时寻找最优解，近年来遗传算法在神经网络优化方面得到了广泛应用，研究表明遗传算法可以在一定程度上解决BP神经网络训练困难的问题，从而提高模型的预测精度[10-11]，例如张天鹏等人[12] 基于GA-BP神经网络来预测FSW接头的性能，结果表明，实际值与预测值之间的误差在5%左右，具有较高的预测精度；关子奇等人[13] 采用GA遗传算法来优化BP神经网络模型来预测焊接熔池照度，试验结果表明，优化后的GA-BP神经网络模型对焊接熔池照度的预测更为精准；张永波等人[14] 基于GA-BP神经网络对带钢焊接质量进行预测，结果表明，与传统的BP神经网络模型相比，通过遗传算法对BP神经网络的权值和阈值进行优化，改良后的模型无论是在精度上还是在收敛速度上都有了大幅度的提升.
文中探讨基于遗传算法优化的BP神经网络在316L中厚板TIG焊的焊接残余应力和变形预测中的应用. 首先，采用正交试验的方法获取焊接工艺参数，残余应力和变形之间的典型样本数据；然后，介绍传统BP神经网络的原理和在焊接残余应力及变形预测中的应用；其次，介绍遗传算法的基本原理和应用优势，并探讨其在BP神经网络优化中的作用；最后，结合具体的焊接试验对两种神经网络模型的预测结果分别进行验证和误差计算，并对基于遗传算法优化的BP神经网络，在焊接残余应力和变形预测中的预测精度以及该方法的收敛速度等优缺点进行讨论.
1. 试验设计
1.1 试验方法
试验使用的焊机为松下YC-500WX5，焊丝牌号为TGS-316L，保护气体为全氩气，焊接试板材料为316L不锈钢，尺寸为300 mm × 300 mm × 20 mm，采用U形7.5°坡口，钝边为1.5 mm，装配间隙为2.5 mm，施焊20层，每层为多道焊，每层焊厚不超过2 mm. 打底层采用的焊接电流为60 A，电弧电压为12 V，焊接速度为50 mm/min，盖面层采用的焊接电流为180 A，电弧电压为17 V，焊接速度为80 mm/min.
为了测量不同填充层焊接工艺参数对焊接变形的影响，在考虑焊接工装夹的情况下，最终采用反变形的方法来测量计算每块试验样板最终的角变形，如图1所示，并通过点焊的方式将每块试板的反变形角度统一固定为20°，此外，为了焊接时能够避免起弧和落弧留在试板上，造成对接板件边缘焊接缺陷，保证所有试板宽度范围内的焊缝都能够达到相同的强度，在每对焊接试板的焊缝两端各放置一块引弧板，通过点焊连接，见图1(b).
![]()
图 1 焊接试板反变形示意图Figure 1. Welding indication. (a) Sketch of striking arc, ending arc and spot welding position; (b) Sketch of the anti-deformation of the welding test plate对焊缝填充层焊接工艺参数，采用正交试验设计L16(44)因素水平表，并以填充层采用的焊接电流、电弧电压、焊接速度及层间温度为因素，经过综合分析试验材料的板厚、材质以及焊接顺序等确定焊接工艺参数的水平值，具体参数见表1.
表 1 因素水平表L16(44)Table 1. Factor level table for L16(44)水平
编号焊接电流
I/A电弧电压
U/V焊接速度
v/(mm·min−1)层间温度
T/℃1 A1(140) B1(14) C1(50) D1(80) 2 A2(165) B2(16) C2(90) D2(120) 3 A3(190) B3(18) C3(120) D3(160) 4 A4(215) B4(20) C4(150) D4(200) 1.2 残余应力及变形测量
试板上的残余应力并不是处处相同，测量时采用焊缝方向和垂直于焊缝方向上均匀取点的方法进行测量，图2为沿焊缝横向和纵向的残余应力测试点分布示意图，其中测试点有13个，沿垂直焊缝方向的残余应力测试点分别在距焊缝中心20，50，80，110，140 mm处；而沿着焊缝方向的测试点，即A，B，C，D，E，F，G，H，I，每个点之间间距为30 mm. 由于盲孔法测残余应力的局限性，文中所测残余应力均为试板近表面残余应力，将测试结果进行比较，得出焊缝近表面横纵向最大残余应力.
![]()
图 2 试板的残余应力测试取点示意图(mm)Figure 2. Sketch of the point selection for residual stress testing of the test plate如图3所示，试验采用盲孔法测量316L试板焊缝近表面的残余应力，使用角度尺等工具测量试板焊接角变形，并将测得的数据作为样本数据集应用于焊接残余应力和变形预测模型的训练.
![]()
图 3 盲孔法测量残余应力示意图Figure 3. Schematic diagram of the blind hole method for measuring residual stress1.3 正交试验
按照L16(44)正交试验的方法测得的试验结果见表2，在打底层和盖面层的焊接参数相同的条件下，填充层不同的焊接工艺参数对焊接残余应力和角变形的影响程度是显著不同. 横向最大残余应力第6组最大，为444.5 MPa；第18组最低，为123.9 MPa，相差320.6 MPa；纵向最大残余应力第10组最大，为494.3 MPa；第17组最小，为136.7 MPa，相差357.6 MPa；角变形第14组最大，为9.84°；第一组最小，为6.18°，相差3.66°.
表 2 正交试验结果Table 2. Results of the orthogonal experiment序号 焊接电流
I/A电弧电压
U/V焊接速度
v/(mm·min−1)层间温度
T/℃横向最大残余应力
σx-max/MPa纵向最大残余应力
σy-max/MPa角变形
θ/(°)1 140 14 50 80 124.4 386.7 6.18 2 140 16 90 120 287.6 275.8 7.56 3 140 18 120 160 237.3 241.9 6.93 4 140 20 150 200 413.3 187.6 7.42 5 165 14 90 200 397.6 284.3 7.12 6 165 16 50 160 444.5 365.8 7.88 7 165 18 120 120 315.3 397.2 9.33 8 165 20 150 80 274.3 426.3 8.52 9 190 14 120 120 356.1 468.4 8.78 10 190 16 150 80 333.1 494.3 8.81 11 190 18 50 200 297.6 412.6 9.24 12 190 20 90 160 183.7 471.8 8.38 13 215 14 150 160 244.2 356.1 9.56 14 215 16 120 200 174.3 337.1 9.84 15 215 18 90 120 211.3 387.5 9.36 16 215 20 50 80 193.5 434.8 9.21 17 163 17 120 88 163.1 136.7 8.41 18 175 18 110 165 123.9 256.8 8.65 19 185 15 100 95 263.4 313.7 8.93 20 210 19 135 145 309.8 351.8 9.43 2. 神经网络预测
2.1 BP神经网络
BP神经网络是一种按照误差逆向传播算法训练的多层前馈神经网络，是应用最广泛的神经网络模型之一[15]，具有任意复杂的模式分类能力和优良的多维函数映射能力. 在网络结构上，BP神经网络由输入层、隐含层、输出层组成；在计算方法上，BP神经网络以网络误差平方为目标函数，利用梯度下降法来计算目标函数的最小值，因此，BP神经网络在解决非线性拟合的问题上具有巨大的优势[16] .
BP神经网络的主要特点是信号前向传播，误差反向反馈传播，设定期望输出，网络根据预测输出与期望输出的误差对比结果进行反向传播，调整各神经元的权值和阈值，整个过程反复进行，直到网络样本输出误差小于期望误差精度或训练达到最大迭代次数为止[17] .
图4为神经网络结构示意图，x1，x2，x3，x4是神经网络模型的输入, y1，y2，y3是网络的输出，整个神经网络模型可以看成一个非线性函数，网络的输入x1，x2，x3，x4是自变量，网络的输出y1，y2，y3为因变量，自变量和因变量构成了函数映射关系. 假设
$ f $ 为激活函数的映射关系，将该神经网络的隐藏层神经元激活函数设定为Sigmoid函数, 将输出层神经元激活函数设定为Purelin线性函数.Sigmoid函数公式为
$$ f\left(x\right)=\frac{1}{1 + {e}^{-x}} $$ (1) Purelin线性函数公式为
$$ f\left(x\right)=x $$ (2) 故此时前向传播公式为
$$ {u}_{j}=f\left({\sum }_{i=1}^{n}{v}_{ij}{x }_{i} + {\theta }_{j}^{u}\right)，j={\mathrm{1,2},\cdots ,}m $$ (3) $$ y=f\left({\sum }_{j=1}^{m}{w}_{j}{u}_{j} + {\theta }^{y}\right) $$ (4) 为了推演方便，先拿第k个对象作为目标，假设
$ {y}^{\left(k\right)} $ 为经过神经网络预测出的值，$ {y}^{{\left(k\right)}^{'}} $ 为真实的输出值, 为了得到真实值与预测值之间误差尽可能小的目标，将目标函数设定为$$ J=\sum _{k}^{\infty }{\left({y}^{\left(k\right)}-{y}^{{\left(k\right)}^{'}}\right)}^{2} $$ (5) 式中假设
$ {v}_{ij} $ 为第i个输入变量与第j个隐藏神经元的权重，$ {\mathrm{\theta }}_{j}^{u} $ 为隐藏层u第j个神经元的阈值，或者称为偏置项，$ {w}_{j} $ 为第j个神经元与y连接的权重，$ {\mathrm{\theta }}^{y} $ 为y的偏置，$ {u}_{j} $ 为隐藏层第j个神经元，$ {v}_{ij} $ ，$ {\theta }_{j}^{u} $ ，$ {w}_{j} $ ，$ {\theta }^{y} $ 为需要确定的网络参数. 可以采用梯度下降法对网络参数和目标函数$ J $ 进行权重更新和优化，其中$ \mu $ 为学习率.$$ {\nu }_{ij}^{'}={v}_{ij}-\mu \frac{dJ}{d{v}_{ij}} $$ (6) $$ {w}_{j}^{'}={w}_{j}-\mu \frac{dJ}{d{w}_{j}} $$ (7) $$ {\theta }_{j}^{u'}={\theta }_{j}^{u}-\mu \frac{dJ}{d{\theta }_{j}^{u}} $$ (8) $$ {\theta }^{y'}={\theta }^{y}-\mu \frac{dJ}{d{\theta }^{y}} $$ (9) 式中：
$ {v}_{ij}\mathrm{'} $ ，$ {\theta }_{j}^{u'} $ ，$ {{w}_{j}}' $ ，$ {\theta }^{y'} $ 为更新以后的网络参数，通过不断迭代更新网络参数，当样本训练输出误差小于期望误差或者达到最大迭代次数时，得到的网络参数就是训练之后的最佳神经网络模型.目标函数对各神经网络中参数的微分由反向梯度下降公式求得
$$ \frac{\partial {J}^{\left(k\right)}}{\partial {w}_{j}}=2\left({y}^{\left(k\right)}-{y}^{\left(k\right)'}\right){u}_{j} $$ (10) $$ \frac{\partial {J}^{\left(k\right)}}{\partial {\theta }^{y}}=2\left({y}^{\left(k\right)}-{y}^{\left(k\right)'}\right) $$ (11) $$ \frac{\partial {J}^{\left(k\right)}}{\partial {\nu }_{ij}}=2\left({y}^{\left(k\right)}-{y}^{\left(k\right)'}\right){w}_{j}{u}_{j}\left(1-{u}_{\dot{j}}\right){x}_{i}^{\left(k\right)} $$ (12) $$ \frac{\partial {J}^{\left(k\right)}}{\partial {\theta }_{j}^{u}}=2\left({y}^{\left(k\right)}-{y}^{\left(k\right)'}\right){w}_{j}{u}_{j}\left(1-{u}_{\dot{j}}\right) $$ (13) 在BP神经网络训练学习过程中，网络的初始权值对学习效果和收敛速度有着至关重要的影响，合适的初始值可以使网络模型较快收敛到最优解，如果初始权值和阈值选取不当，会导致网络的收敛速度慢并且容易陷入局部最小值，由于BP神经网络的初始权值是随机赋予的，因此非常容易出现收敛速度慢且易陷入局部最小值的问题[18] .
采用遗传算法来优化BP神经网络的初始权值和阈值，以提高网络输出的准确性和速度，从而加速BP神经网络的收敛速度，解决局部极值问题，再通过BP神经网络的反向传播算法来寻找整个网络参数的最优解，以实现更好的预测效果.
2.2 GA-BP神经网络
遗传算法是一种模拟自然界遗传机制，以达尔文的生物进化论为理论基础，通过数学利用仿真运算的过程搜索最优解算法. 利用遗传算法，首先需要创建初始种群，选择种群规模为
$ M $ ，取网络所有权值和阈值组成$ N $ 维向量，随机生成一个$ M\cdot N $ 的矩阵作为初始种群，其中N=inputnumber·hiddennumber + hiddennumber + hiddennumber·outputnumber + outputnumber，适应度是遗传算法中衡量个体性能的主要指标，根据适应度的大小对个体进行筛选，因此，适应度是遗传算法的动力. 由于遗传算法的复杂度主要取决于适应度函数的复杂度，因此适应度函数的设计应尽可能简单，以减少算法的计算时间.遗传算法首先需要对初始值进行编码，由于BP神经网络采用了Sigmoid函数作为激活函数，其权值和阈值介于−1 ~ 1，因此采用实数编码会比二进制编码更合适. 误差与适应度为相反的关系，当适应度越好时，预测的结果越准确，误差值越小，因此采用均方误差的倒数作为适应度函数.
适应度函数的公式为
$$ F\left(x\right)=\frac{1}{J\left(x\right)} $$ (14) 遗传算法主要包括选择操作、交叉操作和变异操作，在进行选择操作时，个体被选中的概率与适应度相关，当个体的适应度越大，被选中的可能性就越大，因此最常用到的方法就是轮盘赌法，每个个体被选择的概率为
$$ {p}_{i}=\frac{{F}_{i}}{\displaystyle\sum _{i=1}^{N}{F}_{i}} $$ (15) 式中：
$ {F}_{i} $ 为个体$ i $ 的适应度函数值；$ N $ 为种群个体数量.在进行交叉操作时，将第
$ k $ 条染色体$ {a}_{k} $ 和第$ l $ 条染色体$ {a}_{l} $ 同时在$ j $ 位置进行交换，交换方式如下$$ {a}_{kj}={a}_{kj}\left(1-b\right) + {a}_{lj}b $$ (16) $$ {a}_{lj}={a}_{lj}\left(1-b\right) + {a}_{kj}b $$ (17) 式中：b是[0,1]的随机数.
在进行变异操作时，通过选择第
$ i $ 个个体的第$ j $ 个基因$ {a}_{ij} $ 进行变异，具体操作方法如下$$ {a}_{ij}=\left\{\begin{array}{c}{a}_{ij} + {(a}_{ij}-{a}_{\max})\cdot f(g),r > 0.5\\ {a}_{ij} + {(a}_{\min}-{a}_{ij})\cdot f(g),r\leqslant 0.5\end{array}\right. $$ (18) $$ f\left(g\right)=r{\cdot \left(1-\frac{g}{{G}_{\max}}\right)}^{2} $$ (19) 式中：
$ {a}_{\max} $ 与$ {a}_{\min} $ 分别是基因$ {a}_{ij} $ 的上限和下限；$ r $ 是0 ~ 1之间的随机数；$ g $ 是当前迭代次数；$ {G}_{\max} $ 是最大进化次数.在传统的BP神经网络之前，使用遗传算法进行编码、选择、交叉和变异操作，以优化网络的初始权值和阈值，这种改进避免了传统生物神经网络随机初始化带来的问题，例如网络陷入局部最小值. 遗传算法可以确定最佳的初始权值和阈值，从而解决了传统BP神经网络的局限性，图5为GA-BP神经网络的算法流程图.
3. 仿真试验验证
3.1 数据预处理
经过正交试验L16(44)共获得16组残余应力和角变形的实测数据，再分别采用不同的焊接工艺参数随机进行4组试验，一共获得20组试验数据. 试验取前16组数据作为训练样本，用于GA-BP神经网络的学习和训练，取后4组随机的试验数据作为验证样本，用于验证神经网络的预测精度和泛化能力.
考虑到输入层各个参数之间差异比较大，故对输入层的样本数据采用Mapminmax函数进行归一化处理，将数值变换到0 ~ 1，目的在于消除因量纲不同对预测结果造成的影响，具体消除方法如下
$$ {x}^{*}=\frac{\left(x-{x}_{\min}\right)({y}_{\max}-{y}_{\min})}{{x}_{\max}-{x}_{\min}} + {y}_{\min} $$ (20) 式中：
$ {y}_{\max} $ 和$ {y}_{\min} $ 为设定的上下限，分别为1和0；$ x $ 是待归一化的样本数据；$ {x}_{\max} $ 和$ {x}_{\min} $ 分别为样本数据的最大值和最小值；$ {x}^{*} $ 是样本数据归一化处理后得到的结果.3.2 仿真试验
为了验证改进后的遗传算法对焊接残余应力和变形有良好的预测效果，因此需要进行仿真试验，数据集为上述预处理后的数据，在Matlab2021b中对预测模型进行训练，结果表明，利用遗传算法优化后的BP神经网络模型，对焊接残余应力和变形的预测在预测精度和泛化能力方面比传统的BP神经网络效果要好.
文中的仿真试验，通过焊接电流、电弧电压、焊接速度以及层间温度来预测焊接之后的残余应力和变形，因此输入层节点为4个，输出层节点为3个，关于隐含层的节点数根据经验公式手动设定
$$ m_1=\sqrt{a + b} + k $$ (21) 式中：
$ m_1 $ 为隐含层节点数；$ a $ 和$ b $ 分别为输入层和输出层的节点数；$ k $ 为1 ~ 10之间的常数.3.3 参数设定
采用神经网络对模型训练时，隐藏层设定为两层，在Matlab中的基本参数设置如下.
(1) net.trainParam.epochs=100，迭代次数为100次.
(2) net.trainParam.lr=0.1，学习率为0.1.
(3) net.trainParam.goal=0.001，训练目标误差为0.001.
设置遗传算法的迭代次数为50，种群规模为20，发生交叉的初始概率为0.6，变异概率为0.4.
3.4 仿真结果及分析
使用Matlab2021b软件分别对BP神经网络和GA-BP神经网络进行训练，预测结果与实际值比较结果如图6所示.
![]()
图 6 BP和GA-BP神经网络训练和预测结果Figure 6. Training and prediction results of BP and GA-BP neural networks. (a) maximum lateral residual stress of the training sample； (b) maximum longitudinal residual stress of training sample； (c) training sample deformation angle； (d) test the maximum lateral residual stress of the sample； (e) longitudinal maximum residual stress of the sample is tested； (f) test sample deformation angle通过对比优化前后神经网络模型的预测输出以及期望输出分析神经网络的预测效果，从图6中可以看出，BP神经网络模型的预测较不稳定，在个别样本点的预测误差较大，在少数样本点以及测试点的预测精度较高，因此BP神经网络模型的数据拟合能力和预测性都不够稳定. 对比图6中这6张图不难看出，无论是训练样本还是测试样本，GA-BP神经网络的预测值比BP神经网络更接近于实际测量值，拟合程度更高，预测效果更好.
BP和GA-BP神经网络适应度曲线如图7所示，GA-BP神经网络的收敛速度和训练效率都明显高于BP神经网络，通过对比两种模型的适应度曲线，可以看出BP网络的曲线波动更大，GA-BP网络的曲线更加平滑，这表明优化后的网络模型收敛更快，训练更加稳定，优化后网络的最佳适应度值更低，模型的预测结果与实际结果之间的差异更小，预测结果更好.
![]()
图 7 适应度变化曲线Figure 7. Curves of fitness changes. (a) BP neural network; (b) GA-BP neural network神经网络预测平均相对误差见表3，BP神经网络对横、纵向残余应力和角变形的训练样本预测误差分别为12.15%，12.52%，6.8%，GA-BP神经网络的训练样本预测误差分别为2.51%，2.01%，1.28%，由此可见优化后的网络模型预测效果显著提高. BP神经网络的测试集的预测误差分别为14.56%，14.35%，8.53%，GA-BP神经网络的预测误差分别为2.38%，2.84%，1.32%，优化后的网络模型无论是在训练集还是测试集上的表现都优于传统的BP神经网络. GA-BP神经网络预测模型在预测结果上有了一定的改善，在预测精度上明显优于传统的BP神经网络模型，平均相对误差在3%以内.
表 3 神经网络预测平均相对误差Table 3. Average relative error of neural network prediction输出结果 BP GA-BP 训练样本 测试样本 训练样本 测试样本 横向最大残余应力 12.15% 14.56% 2.51% 2.38% 纵向最大残余应力 12.52% 14.35% 2.01% 2.84% 角变形 6.8% 8.53% 1.28% 1.32% 试验结果表明，根据正交试验和预测值不难发现相对于电弧电压，焊接电流对残余应力和角变形的影响更大，随着电流电压的增大，残余应力呈现上升趋势，同时变形也在增大.
通过对比表中的平均相对误差表值，可以得出GA-BP模型比BP网络模型的误差小，预测的精度更高，通过采用随机试验的方式验证模型的预测精度和误差，可以发现模型的泛化能力很好，并且优化后模型的泛化能力比传统BP神经网络模型表现更加优异，通过对比两组模型的训练迭代次数和适应度曲线图，可以发现GA-BP的收敛速度要明显好于BP，因此GA-BP的收敛性更好，总体效果更好.
4. 结论
(1) BP神经网络的预测相对误差在15%以内，而GA-BP神经网络的预测相对误差能控制在3%以内，因此GA-BP神经网络比传统的BP神经网络预测效果更好，拟合能力更突出，并且模型的收敛速度更快，训练过程更稳定.
(2) GA-BP神经网络模型对横、纵向残余应力以及角变形的真实值随机测试预测误差分别为2.38%，2.84%，1.32%，结果表明，GA-BP神经网络能够精准得反映多层多道焊填充层焊接工艺参数与焊后残余应力之间的非线性映射关系，相对误差小，具有较强的预测能力和泛化能力，该方法对奥氏体不锈钢中厚板多层多道焊填充层焊接工艺参数具有指导作用，并对焊后残余应力和变形控制提供了一种新的思路和方法，在焊接过程中发挥了重要作用.
-
![]()
图 2 无源电探针检测示意图
Figure 2. Schematic diagram of passive electrical probe detection. (a) simultaneous acquisition location of the triple-probe; (b) probe structure diagram; (c) probe detection distance
![]()
图 3 三路探针同步采集到的原始信号示意图
Figure 3. Schematic diagram of the original signals captured simultaneously by the triple-probe
![]()
图 4 等离子体电信号波形提取示意图
Figure 4. Schematic diagram for the extraction of plasma electrical signal waveforms. (a) plasma electrical signal waveform; (b) schematic diagram of feature point extraction for different behaviors of plasma
![]()
图 5 等离子体不同动态行为与探针的相对位置变化示意图
Figure 5. Schematic diagram between dynamic behaviors of plasma and changes in the relative position of the probe. (a) eruption; (b) sway; (c) deflective eruption or sway
![]()
图 6 等离子体偏置喷发和摆动行为示意图
Figure 6. Schematic diagram of plasma deflective eruption and sway. (a) deflective eruption; (b) sway
![]()
图 7 不同焊接速度下等离子体喷发温度变化率的概率分布
Figure 7. Probability distribution of the temperature change rate of plasma eruption
![]()
图 8 不同焊接速度下等离子体摆动温度变化率的概率分布
Figure 8. Probability distribution of the temperature change rate of plasma sway
![]()
图 9 不同等离子体动态行为的温度变化率概率分布
Figure 9. Probability distribution of the temperature change rates for different dynamic behaviors of plasma
![]()
图 10 等离子体不确定行为温度变化率的聚类形状
Figure 10. Cluster shapes of temperature change rates for uncertain plasma behavior
![]()
图 11 K-均值聚类分析后的不确定行为特征分布
Figure 11. Distribution of uncertain behavior after K-means clustering analysis
![]()
图 12 不同焊接速度下等离子体动态行为的变化趋势
Figure 12. Trends in plasma dynamic behavior at different welding speeds
![]()
图 13 不同焊接速度下等离子体喷发的时频域分布
Figure 13. Frequency of plasma eruption at different welding speeds. (a) time-domain variation diagram of eruption frequency at different welding speeds; (b) probability distribution diagram of eruption frequency
![]()
图 14 不同焊接速度下等离子体摆动的时频域分布
Figure 14. Frequency of plasma sway at different welding speeds. (a) time-domain variation diagram of sway frequency at different welding speeds; (b) probability distribution diagram of sway frequency
表 1 不同焊接速度下等离子体各动态行为的平均频率
Table 1 The average frequency of dynamic behaviors of plasma at different welding speeds
焊接速度v/(mm·s−1) 喷发频率f1/Hz 摆动频率f2/Hz 15 1310 951 20 1510 1026 25 1712 1131 
下载: 导出CSV
-
[1] Beck M, Berger P, Hugel H. The effect of plasma formation on beam focusing in deep penetration welding with CO2 lasers[J]. Journal of Physics D: Applied Physics, 1995, 28: 2430 − 2442. doi: 10.1088/0022-3727/28/12/007
[2] Zhao C, Parab N D, Li X, et al. Critical instability at moving keyhole tip generates porosity in laser melting[J]. Science, 2020, 370(6520): 1080 − 1086. doi: 10.1126/science.abd1587
[3] Seely J. Quantum theory of inverse bremsstrahlung absorption and pair production[J]. Laser Interaction and Related Plasma Phenomena, 1974: 835 − 847.
[4] Meng W, Yin X, Fang J, et al. Dynamic features of plasma plume and molten pool in laser lap welding based on image monitoring and processing techniques[J]. Optics and Laser Technology, 2019, 109: 168 − 177. doi: 10.1016/j.optlastec.2018.07.073
[5] Yusof M, Ishak M, Ghazali M. Weld depth estimation during pulse mode laser welding process by the analysis of the acquired sound using feature extraction analysis and artificial neural network-ScienceDirect[J]. Journal of Manufacturing Processes, 2020, 63: 163 − 178
[6] Li J, Zhang Y. Prediction of penetration based on plasma plume and spectrum characteristics in laser welding[J]. Journal of Manufacturing Processes, 2022, 75: 593 − 604. doi: 10.1016/j.jmapro.2022.01.032
[7] Tenner F, Brock C, Klämpfl F, et al. Analysis of the correlation between plasma plume and keyhole behavior in laser metal welding for the modeling of the keyhole geometry[J]. Optics and Lasers in Engineering, 2015, 64: 32 − 41. doi: 10.1016/j.optlaseng.2014.07.009
[8] Pang S, Chen X, Shao X, et al. Dynamics of vapor plume in transient keyhole during laser welding of stainless steel: Local evaporation, plume swing and gas entrapment into porosity[J]. Optics and Lasers in Engineering, 2016, 82: 28 − 40. doi: 10.1016/j.optlaseng.2016.01.019
[9] Zhao S, Yang L, Liu T, et al. Analysis of plasma oscillations by electrical detection in Nd: YAG laser welding[J]. Journal of Materials Processing Technology, 2017, 249: 479 − 489. doi: 10.1016/j.jmatprotec.2017.06.032
[10] 贾晨鹏, 黄一鸣, 杨立军, 等. TC4钛合金激光深熔焊等离子体羽流热力学行为研究[J]. 中国激光, 2023, 50(20): 2002105. Jia Chenpeng, Huang Yiming, Yang Lijun, et al. Plasma plume thermodynamic behavior during laser deep-penetration welding of TC4 titanium alloy[J]. Chinese Journal of Lasers, 2023, 50(20): 2002105.
[11] 赵圣斌, 杨立军, 刘桐, 等. 不同焊接模式下的A304不锈钢YAG激光焊等离子体的电信号特征[J]. 中国激光, 2016, 43(12): 1202005. Zhao Shengbin, Yang Lijun, Liu Tong, et al. Electrical signal characteristics of plasma in YAG laser welding of A304 stainless steels under different modes[J]. Chinese Journal of Lasers, 2016, 43(12): 1202005.
[12] Kim C H, Ahn D C. Coaxial monitoring of keyhole during Yb: YAG laser welding[J]. Optics and Laser Technology, 2012, 44(6): 1874 − 1880. doi: 10.1016/j.optlastec.2012.02.025
[13] Qiu W, Yang L, Zhao S, et al. A study on plasma plume fluctuation characteristic during A304 stainless steel laser welding[J]. Journal of Manufacturing Processes, 2018, 33: 1 − 9.
[14] 段爱琴, 陈俐, 王亚军, 等. CO2激光焊接不锈钢光致等离子体动态特性分析[J]. 焊接学报, 2005, 26(11): 17 − 20. Duan Aiqin, Chen Li, Wang Yajun, et al. Analysis of dynamic characteristics of photoinduced plasma of CO2 laser-welded stainless steel[J]. Transactions of the China welding institution, 2005, 26(11): 17 − 20.
[15] 高向东, 吕威兴, 游德勇, 等. 大功率盘形激光焊接过程等离子体图像特征分析[J]. 焊接学报, 2011, 32(12): 5 − 8. Gao Xiangdong, Lü Weixing, You Deyong, et al. Analysis of plasma image characteristics in high-power disc laser welding process[J]. Transactions of the China welding institution, 2011, 32(12): 5 − 8.
[16] Iranitalab A, Khattak A. Comparison of four statistical and machine learning methods for crash severity prediction[J]. Accident Analysis & Prevention, 2017, 108: 27 − 36.
[17] Dash B, Mishra D, Rath A, et al. A hybridized K-means clustering approach for high dimensional dataset[J]. International Journal of Engineering Science and Technology, 2010, 2: 59 − 66.
计量
- 文章访问数: 83
- HTML全文浏览量: 4
- PDF下载量: 28
