Thermal-fluid behavior and distribution of ceramic particles during laser melting deposition of TiC/TC4 composite materials
-
摘要:
颗粒增强钛基复合材料在航空航天等领域具有广阔的应用前景，采用仿真与试验相结合的方法，对激光熔化沉积TiC/TC4复合材料过程中的热-流行为和沉积层中陶瓷颗粒的分布状态进行了研究. 结果表明，工艺参数对熔池的热-流行为及沉积层形貌具有显著的影响，熔池表面的流体流动呈波纹状，熔融金属从熔池中心向四周流动；熔池底部等温面上熔融金属从四周流向熔池底部，并出现漩涡现象. 在激光熔化沉积过程中，TiC颗粒穿透熔池表面的Marangoni对流区，在熔池中与流体交互作用，最终在沉积层中出现沉积层中上方部位团聚、沉积层中下方部位团聚、沉积层中均匀分布等分布状态.
Abstract:Particulate reinforced titanium matrix composites have broad application prospects in aerospace and other fields. The thermal-fluid behavior and the distribution of the ceramic particles in the deposited layer during the laser melting deposition of TiC/TC4 composites were studied by combining simulation and experiment. The results show that the process parameters have a significant effect on the thermal-fluid behavior of the molten pool and the morphology of the deposited layer. With the increase of laser power, the molten pool size and the maximum fluid flow rate increase. As the scanning speed increases, the molten pool size and the maximum fluid flow rate decrease. The fluid flow on the surface of the molten pool is corrugated. The molten metal flows from the center of the molten pool to the periphery. The molten metal on the isothermal surface of the bottom of the molten pool flows from the periphery to the bottom, and the vortex flow phenomenon occurs. In the process of laser melting deposition, TiC particles penetrate the Marangoni convection zone on the surface of the molten pool and interact with the fluid in the molten pool. Finally, they may appear in the distribution of the agglomeration in middle-upper part, the agglomeration in middle-lower part and uniform distribution in the deposited layer.
-
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神经网络能够精准得反映多层多道焊填充层焊接工艺参数与焊后残余应力之间的非线性映射关系，相对误差小，具有较强的预测能力和泛化能力，该方法对奥氏体不锈钢中厚板多层多道焊填充层焊接工艺参数具有指导作用，并对焊后残余应力和变形控制提供了一种新的思路和方法，在焊接过程中发挥了重要作用.
-
![]()
图 1 TiC/TC4复合材料激光熔化沉积模型
Figure 1. Laser melting deposition model of TiC/TC4 composite materials. (a) geometric model structure; (b) geometric model meshing; (c) partial meshing
![]()
图 2 激光熔化沉积过程熔池形成及长大过程
Figure 2. Formation and growth of molten pool in laser melting deposition process. (a) t = 40 ms; (b) t = 80 ms; (c) t = 120 ms; (d) t = 160 ms; (e) t = 200 ms; (f) t = 240 ms; (g) t = 280 ms; (h) t = 320 ms; (i) t = 360 ms
![]()
图 3 不同激光功率下激光熔化沉积过程中熔池形貌及表面流场
Figure 3. Morphology and surface flow field of molten pool during laser melting deposition at different laser powers. (a) P = 1 200 W, vs = 10 mm/s; (b) P = 1 400 W, vs = 10 mm/s; (c) P = 1 600 W, vs = 10 mm/s; (d) P = 1 400 W, vs = 8 mm/s; (e) P = 1 400 W, vs = 12 mm/s
![]()
图 4 激光熔化沉积试验设备示意图
Figure 4. Schematic diagram of laser melting deposition test equipment.
![]()
图 5 不同工艺参数下熔覆层横截面形貌
Figure 5. Cross-section morphology of cladding layer under different process parameters. (a) P = 1 400 W, vs = 8 mm/s; (b) P = 1 400 W, vs = 10 mm/s; (c) P = 1 400 W, vs = 12 mm/s; (d) P = 1 200 W, vs = 10 mm/s; (e) P = 1 600 W, vs = 10 mm/s
![]()
图 6 不同工艺参数下沉积层横截面尺寸
Figure 6. Cross-section dimensions of the deposited layer under different process parameters. (a) schematic diagram of cross-section dimension measurement; (b) variation of cross-section dimension with laser power; (c) variation of cross-section dimension with scanning speed
![]()
图 7 不同工艺参数下沉积层中TiC颗粒的分布状态
Figure 7. Distribution of TiC particles in the deposited layer under different process parameters. (a) P = 1 400 W, vs = 8 mm/s; (b) P = 1 400 W, vs = 10 mm/s; (c) P = 1 400 W, vs = 12 mm/s; (d) P = 1 200 W, vs = 10 mm/s; (e) P = 1 600 W, vs = 10 mm/s
![]()
图 8 熔池中粒子迁移机制
Figure 8. Mechanism of particle migration in molten pool. (a) particles enter the gas/liquid interface; (b) particle migration mechanism near the back end of the molten pool; (c) particle migration mechanism near the front end of the molten pool; (d) migration mechanism of cross-section of the molten pool
表 1 TiC/TC4复合材料激光熔化沉积试验工艺参数
Table 1 Design table of technological parameters for laser melting deposition test of TiC/TC4 composite materials
编号 激光功率
P/W扫描速度
vs /(mm·s−1)TiC 含量
w(%)送粉率
vf /(g·min−1)1 1400 8 20 9.6 2 1400 10 20 9.6 3 1400 12 20 9.6 4 1200 10 20 9.6 5 1600 10 20 9.6 
下载: 导出CSV
-
[1] Gao Zhuanni, Bu Hengchang, Feng Yu, et al. Strengthening mechanism of Y2O3 nanoparticles on microstructure and mechanical properties of the laser additive manufacturing joint for large thickness TC4 titanium alloy[J]. Journal of Manufacturing Processes, 2021, 71: 37 − 55. doi: 10.1016/j.jmapro.2021.09.011
[2] 冯进宇, 肖华强, 肖易, 等. TC4表面激光熔覆Ti-Al-(C, N)复合涂层的组织与力学性能[J]. 中国激光, 2022, 49(2): 159 − 167. Feng Jinyu, Xiao Huaqiang, Xiao Yi, et al. Microstructure and mechanical properties of laser cladded Ti-Al-(C, N) composite coating on TC4 surface[J]. Chinese Journal of Lasers, 2022, 49(2): 159 − 167.
[3] Chen Changrong, Liu Chang, Wang Qianting, et al. Effects of heat input on layer heterogeneity of selective laser melting Ti-6Al-4V components[J]. China Welding, 2023, 32(3): 51 − 66.
[4] 顾冬冬, 张红梅, 陈洪宇, 等. 航空航天高性能金属材料构件激光增材制造[J]. 中国激光, 2020, 47(5): 32 − 55. Gu Dongdong, Zhang Hongmei, Chen Hongyu, et al. Laser additive manufacturing of high-performance metallic aerospace components[J]. Chinese Journal of Lasers, 2020, 47(5): 32 − 55.
[5] 周亦人, 沈自才, 齐振一, 等. 中国航天科技发展对高性能材料的需求[J]. 材料工程, 2021, 49(11): 41 − 50. doi: 10.11868/j.issn.1001-4381.2021.000129 Zhou Yiren, Shen Zicai, Qi Zhenyi, et al. Demand for high performance materials in development of China's aerospace science and technology[J]. Journal of Materials Engineering, 2021, 49(11): 41 − 50. doi: 10.11868/j.issn.1001-4381.2021.000129
[6] 吴诚福, 李新意, 陈洪胜, 等. 激光增材制造WCp钛基复合材料界面连接机理及力学性能[J]. 焊接学报, 2023, 44(3): 44 − 53. doi: 10.12073/j.hjxb.20220425003 Wu Chengfu, Li Xinyi, Chen Hongsheng, et al. Interface connection mechanism and mechanical properties of WCp titanium matrix composites fabricated by laser additive[J]. Transactions of the China Welding Institution, 2023, 44(3): 44 − 53. doi: 10.12073/j.hjxb.20220425003
[7] Yang Chen, Zhao Zhanyong, Bai Peikang, et al. Nano-SiC whisker-reinforced Ti6Al4V matrix composites manufactured by selective laser melting: fine equiaxed grain formation mechanism and mechanical properties[J]. Journal of Materials Processing Technology, 2023, 317: 117981. doi: 10.1016/j.jmatprotec.2023.117981
[8] Xiao Yunmian, Yang Yongqiang, Wang Di, et al. In-situ synthesis of high strength and toughness TiN/Ti6Al4V sandwich composites by laser powder bed fusion under a nitrogen-containing atmosphere[J]. Composites Part B:Engineering, 2023, 253: 110534. doi: 10.1016/j.compositesb.2023.110534
[9] 钦兰云, 门继华, 赵朔, 等. TiB2含量对选区激光熔化TiB/Ti6Al4V复合材料组织及力学性能的影响[J]. 中国激光, 2021, 48(6): 49 − 58. Qin Lanyun, Men Jihua, Zhao Shuo, et al. Effect of TiB2 content on microstructure and mechanical properties of TiB/Ti6Al4V composites formed by selective laser melting[J]. Chinese Journal of Lasers, 2021, 48(6): 49 − 58.
[10] Huo Pengcheng, Zhao Zhanyong, Du Wenbo, et al. Deformation strengthening mechanism of in situ TiC/TC4 alloy nanocomposites produced by selective laser melting[J]. Composites Part B:Engineering, 2021, 225: 109305. doi: 10.1016/j.compositesb.2021.109305
[11] 郭顺, 王朋坤, 顾介仁, 等. 电弧熔炼Ti6Al4V/B4C复合材料微观组织与力学性能[J]. 焊接学报, 2022, 43(9): 62 − 68,117. doi: 10.12073/j.hjxb.20220403002 Guo Shun, Wang Pengkun, Gu Jieren, et al. Microstructure and mechanical properties of Ti6Al4V/B4C composite prepared by arc melting[J]. Transactions of the China Welding Institution, 2022, 43(9): 62 − 68,117. doi: 10.12073/j.hjxb.20220403002
[12] 王建东. 激光熔化沉积TiC/Ti6Al4V复合材料的组织性能调控[D]. 哈尔滨: 哈尔滨工业大学, 2018. Wang Jiandong. Control of microstructure and properties of TiC/Ti6Al4V composites manufactured by laser melting deposition[D]. Harbin: Harbin Institute of Technology, 2018.
[13] Ma Guangyi, Yu Chao, Tang Bokai, et al. High-mass-proportion TiCp/Ti6Al4V titanium matrix composites prepared by directed energy deposition[J]. Additive Manufacturing, 2020, 35: 101323. doi: 10.1016/j.addma.2020.101323
[14] Hu Zhengyang, Cheng Xingwang, Li Shenglin, et al. Investigation on the microstructure, room and high temperature mechanical behaviors and strengthening mechanisms of the (TiB + TiC)/TC4 composites[J]. Journal of Alloys and Compounds, 2017, 726: 240 − 253. doi: 10.1016/j.jallcom.2017.08.017
[15] Wang Tao, Liu Xingyu, Chen Siyu, et al. Study on microstructure and tribological properties of nano/micron TiC/TC4 composites fabricated by laser melting deposition[J]. Journal of Manufacturing Processes, 2022, 82: 296 − 305. doi: 10.1016/j.jmapro.2022.07.068
[16] 王少伟. 选区激光熔化成形TiC/Ti6Al4V复合材料及热处理过程组织性能演化[D]. 太原: 中北大学, 2022. Wang Shaowei. TiC/Ti6Al4V composites by selective laser melting and microstructure evolution following different heat treatment[D]. Taiyuan: North University of China, 2022.
[17] Zhang Y M, Lim C W J, Tang C, et al. Numerical investigation on heat transfer of melt pool and clad generation in directed energy deposition of stainless steel[J]. International Journal of Thermal Sciences, 2021, 165: 106954. doi: 10.1016/j.ijthermalsci.2021.106954
[18] Song Boxue, Yu Tianbiao, Jiang Xingyu, et al. The relationship between convection mechanism and solidification structure of the iron-based molten pool in metal laser direct deposition[J]. International Journal of Mechanical Sciences, 2020, 165: 105207.
[19] Gan Zhengtao, Yu Gang, He Xiuli, et al. Numerical simulation of thermal behavior and multicomponent mass transfer in direct laser deposition of Co-base alloy on steel[J]. International Journal of Heat and Mass Transfer, 2017, 104: 28 − 38. doi: 10.1016/j.ijheatmasstransfer.2016.08.049
[20] Gan Zhengtao, Liu Hao, Li Shaoxia, et al. Modeling of thermal behavior and mass transport in multi-layer laser additive manufacturing of Ni-based alloy on cast iron[J]. International Journal of Heat and Mass Transfer, 2017, 111: 709 − 722. doi: 10.1016/j.ijheatmasstransfer.2017.04.055
[21] Song Boxue, Yu Tianbiao, Jiang Xingyu, et al. Effect of laser power on molten pool evolution and convection[J]. Numerical Heat Transfer, Part A:Applications, 2020, 78(2): 48 − 59. doi: 10.1080/10407782.2020.1777795
[22] 刘伟嵬, 王怡文, 王灏, 等. 不同工作模式下激光熔化沉积过程熔池演变数值研究[J]. 大连理工大学学报, 2021, 61(4): 344 − 354. doi: 10.7511/dllgxb202104003 Liu Weiwei, Wang Yiwen, Wang Hao, et al. Numerical study of evolution of molten pool in laser melting deposition process under different working modes[J]. Journal of Dalian university of Technology, 2021, 61(4): 344 − 354. doi: 10.7511/dllgxb202104003
[23] Wei Shaopeng, Wang Gang, Shin Y C, et al. Comprehensive modeling of transport phenomena in laser hot-wire deposition process[J]. International Journal of Heat and Mass Transfer, 2018, 125: 1356 − 1368. doi: 10.1016/j.ijheatmasstransfer.2018.04.164
[24] Sun Zhe, Guo Wei, Li Lin. Numerical modelling of heat transfer, mass transport and microstructure formation in a high deposition rate laser directed energy deposition process[J]. Additive Manufacturing, 2020, 33: 101175. doi: 10.1016/j.addma.2020.101175
[25] Du Jun, Wu Yunxiao, Wang Daqing, et al. Droplet impact-induced molten pool dynamics and particle migration patterns during TIG-assisted droplet deposition manufacturing of SiC particle-reinforced Al-Si alloy[J]. Case Studies in Thermal Engineering, 2022, 33: 101962. doi: 10.1016/j.csite.2022.101962
[26] 胡勇, 王梁, 娄复兴, 等. 稳态磁场对激光熔注球形WC颗粒分布的影响机理研究[J]. 机械工程学报, 2021, 57(1): 240 − 248. doi: 10.3901/JME.2021.01.240 Hu Yong, Wang Liang, Lou Fuxing, et al. Mechanism study of steady magnetic field effect on spherical WC particle distribution during laser melt injection[J]. Journal of Mechanical Engineering, 2021, 57(1): 240 − 248. doi: 10.3901/JME.2021.01.240
[27] Vreeling J A, Ocelík V, Pei Y T, et al. Laser melt injection in aluminum alloys: On the role of the oxide skin[J]. Acta Materialia, 2000, 48(17): 4225 − 4233. doi: 10.1016/S1359-6454(00)00278-0
[28] Xu Boan, Jiang Ping, Wang Yilin, et al. Multi-physics simulation of wobbling laser melting injection of aluminum alloy with SiC particles: SiC particles gradient distribution in fusion zone[J]. International Journal of Heat and Mass Transfer, 2022, 182: 121960. doi: 10.1016/j.ijheatmasstransfer.2021.121960
[29] Wang Liang, Yao Jianhua, Hu Yong, et al. Influence of electric-magnetic compound field on the WC particles distribution in laser melt injection[J]. Surface and Coatings Technology, 2017, 315: 32 − 43. doi: 10.1016/j.surfcoat.2017.01.116
计量
- 文章访问数: 128
- HTML全文浏览量: 34
- PDF下载量: 36
