摇动电弧窄间隙GMAW焊温度场数值分析模型

胥国祥,潘海潮,王加友

doi:10.12073/j.hjxb.20151210001

摘要: 从宏观传热学出发，依据横断面几何形状特征，综合考虑电弧摇动、接头几何特征及焊缝表面形状对电弧热流分布特征的影响，建立了摇动电弧窄间隙GMAW焊瞬态温度场数值分析模型. 利用ANSYS软件对摇动电弧瞬态温度场及热循环曲线进行了模拟计算，分析了其分布特征. 结果表明，模型能够合理、准确反映摇动电弧的移动轨迹及分布特性，且焊缝横断面的计算结果与试验结果吻合较好，从而证明了模型的准确性. 在焊接电流为300 A，摇动频率为2 Hz时，熔池形态变化有限.

关键词:

0. 序言

焊接是造船的最重要工艺之一，但焊接烟尘、弧光和金属飞溅等存在，使得焊接环境非常恶劣. 焊接自动化已经成为船舶先进制造的必由之路^[1-2]. 机器人高质量、高效率焊接的前提是焊接工艺参数的准确、快速选取，而这需要对焊接件数字模型中的每道焊缝进行特征识别及参数提取，并通过与焊接工艺数据库中参数进行智能匹配实现最优工艺参数获取，进而提高焊接效率和质量^[3]. 因此，开展焊缝特征识别以及焊缝特征参数提取是实现机器人全自动焊接的核心和难点. 但目前已有研究主要是针对焊缝图像或数字点云模型进行处理为主. Zhu等人^[4]基于线结构光的三维视觉测量系统，提出了一种焊缝特征点提取算法，实现了对接接头的焊缝宽度和余高的提取. 岳建锋等人^[5]对采集到的焊缝图像进行特征分析，获取焊缝特征点的图像坐标，进而提取出V形坡口的焊缝特征参数. Wen等人^[6]采用结构光三角测量方法，实现了对窄间隙平板对接接头的坡口尺寸提取，包括坡口宽度、中心线位置和坡口夹角等信息. Wang等人^[7]提出了一种基于激光焊接的焊缝高精度测量，通过灰度值统计和平面拟合的方法，完成对焊缝中心线、焊缝宽度和焊缝法向量的提取. 通过传感器搭建的视觉系统，虽能实现坡口特征参数的实时提取，但会受制于传感器扫描范围的限制，且特征参数提取耗时较长. 此外，张泽南等人^[8]通过SolidWorks API接口进行焊接件标注模块的二次开发，针对平板对接和T形接头，利用遍历特征树技术访问到切割清单，从而提取出焊缝厚度和长度信息. 从已有研究可以看出，目前焊缝特征参数提取所针对的接头或坡口比较单一，应用范围较窄. 而所提取的焊缝特征参数并不完整，所获取特征信息有限，不具备对焊接工艺参数进行全自动化选择条件. 为此文中提出了基于空间位置和轮廓线距离的船舶焊缝特征参数提取，首先基于接头空间位置完成接头形式和坡口类型的焊缝特征识别；然后根据已建的3种数学模型进行特征参数提取，进而获取坡口处焊缝间隙、夹角和板厚等信息，数据测试结果验证了文中算法的有效性.

1. 基于STL焊缝特征识别

1.1 基本概念

焊缝特征^[9-10]主要包括焊接接头形式和焊接坡口类型. 焊接接头形式通常可分为平板对接、角接、T形和搭接4类接头. 焊接坡口类型与焊接接头形式直接相关，对于平板对接接头，常见坡口类型^[11-12]主要有I形坡口、V形坡口和Y形坡口；对于角接接头和T形接头，主要有不开坡口、单边坡口和双边坡口；对于搭接接头则常为不开坡口.

焊缝特征识别主要基于三维模型的STL文件进行焊缝接头形式和坡口类型识别，而两接头的待判定面位置关系直接影响接头形式. 如图1所示，由于待判定面相对于坡口端面的面积较大，利用海伦公式^[13]先求出每个轮廓面面积并取最大的前四组，即可得到待判定面；然后根据每个待判定面处法向量(虚线)平行或垂直关系进行4类接头的空间位置关系判别. 若法向量出现平行次数为6次且方向相反次数为4次，则两接头处于平行状态；若垂直次数为4次且方向相反次数为2次，则处于垂直状态.

图 1 待判定面示意图

Figure 1. Schematic diagram of the face to be determined. (a) schematic diagram of flat joint; (b) schematic diagram of corner joint; (c) schematic diagram of T-joint

下载: 全尺寸图片幻灯片

1.2 识别原理与算法

文中针对4类焊接接头形式，共10种焊接坡口进行焊缝特征识别. 根据接头轮廓线总条数，通常可分为24条、27条和30条3大类.

当轮廓线条为24条时，主要有平板对接接头的I形坡口、V形坡口；角接接头的不开坡口、单边坡口；T形和搭接接头的不开坡口六种类型^[14-15]，图2a即属于24条轮廓线类型中平板对接接头所对应的V形坡口，此类型所涉及焊缝特征识别，主要是先取出前8组长度最短的轮廓线. 当轮廓线条为27条时，只有角接接头双边坡口和T形接头单边坡口两种^[16]，图2b即属于27条轮廓线类型中角接接头所对应的双边坡口，焊缝特征识别主要是针对前10组最短轮廓线. 当轮廓线条为30条时，同样只有平板对接接头Y形坡口和T形接头双边坡口两种，图2c即属于30条轮廓线类型中T形接头所对应的双边坡口，焊缝特征识别主要是针对前12组最短轮廓线. 然后针对上述所取的最短轮廓线组，计算两两轮廓线的端点间所有距离，去除轮廓线本身长度后，将最小距离与次短长度的轮廓线对比，若大于次短轮廓线，则称轮廓线距离为大间距；反之称为小间距. 最后根据接头空间位置和最小轮廓线距离进行焊缝特征识别，若接头平行，当距离为小间距时，则为平板对接接头，反之为搭接接头；若接头垂直，当距离为小间距时，则为角接接头，反之为T形接头.

图 2 轮廓线条选取示意图

Figure 2. Schematic diagram of contour line selection. (a) first 8 groups of contour lines; (b) first 10 groups of contour lines; (c) first 12 groups of contour lines

下载: 全尺寸图片幻灯片

当轮廓线条为24条时，搭接和T形接头只有不开坡口一种类型. 平板对接和角接接头的坡口则通过两条最小距离轮廓线的夹角α进行判别. 若接头形式为平板对接，当α = 0时，则为I形坡口，反之为V形坡口；若接头形式为角接，当α = 90°时，则为不开坡口；反之为单边坡口. 当轮廓线条为27条时，若轮廓线距离为小间距，则为角接接头的双边坡口；反之为T形接头的单边坡口. 当轮廓线条为30条时，若接头处于平行状态且小间距，则为平板对接接头的Y形坡口；若接头处于垂直状态且大间距，则为T形接头的双边坡口.

2. 焊缝特征参数提取

基于上述识别出的接头形式和坡口类型，文中针对平板对接、角接和T形接头建立3种数学模型，进而提取出每类接头中的3种坡口焊缝特征参数. 根据所分的3类轮廓线总条数进行不同形式的处理，如图2所示，若轮廓线条为24条，取前8组长度最短的轮廓线组；若轮廓线条为27条，取前10组最短轮廓线组；若轮廓线条为30条，则取前12组最短轮廓线组. 将上述取出的最短轮廓线组存储到数组空间P，计算两两轮廓线的端点间所有距离并去除轮廓线本身长度. 若轮廓线条为24条，则将求距结果(附带距离两端点坐标值)存储到数组空间Q；若轮廓线条为27或30条，则去除两条轮廓线平行时，轮廓线的4个端点及4个端点所在轮廓线其余端点的求距结果，将剩余求距结果(附带距离两端点坐标值)存储到数组空间Q.

2.1 平板对接接头三种坡口提取

为获取焊缝坡口处特征尺寸，先对焊接件数字模型中的焊缝进行特征点及线提取. 图3为平板对接接头尺寸标注示意图，θ₁为向量A₁B₁和向量B₁C₁的夹角；θ₂为向量A₂B₂和向量B₂C₂的夹角. 平板对接接头共有三种坡口，若为I形坡口，则点B₁，C₁两点重合，点B₂，C₂重合，即l_C1B1 = 0，l_C2B2 = 0，只需确定A₁，A₂，C₁，C₂ 4点坐标；若为V形坡口，则点B₁，A₁重合，点B₂，A₂重合，即l_A1B1 = 0，l_A2B2 = 0，只需确定A₁，A₂，C₁，C₂四点坐标；若为Y形坡口，则需确定A₁，A₂，B₁，B₂，C₁，C₂六点坐标.

图 3 平接接头尺寸示意图

Figure 3. Schematic diagram of butt welding joint size

下载: 全尺寸图片幻灯片

2.1.1 特征点及线提取

根据轮廓线总条数确定平板对接接头的3种坡口类型. 若为24条，即平板对接接头为I形或V形坡口，则取数组空间Q最小距离所对应的两端点，即为A₁，A₂. 遍历数组空间P，点A₁所在轮廓线的另一端点即为C₁，点A₂所在轮廓线的另一端点即为C₂.

轮廓线总条数若为30条，即平板对接接头Y形坡口，则遍历数组空间P所有轮廓线的两端点，将重复点两两求距，最小距离所对应的两点，即为B₁，B₂. 分别将组成交点B₁，B₂的轮廓线其余端点两两求距，最小距离所对应的两点即为A₁，A₂，剩余两端点则为C₁，C₂.

2.1.2 特征尺寸提取

坡口夹角θ根据所提取坡口处特征点及线，计算出夹角θ₁，θ₂，由式(1)可求得坡口夹角θ，即

$$ \theta = {\theta _1} + {\theta _2} $$

(1)

式中：θ₁为左坡角；θ₂为右坡角.若轮廓线总条数为24条，则向量B₁C₁等同于向量A₁C₁，向量A₂B₂等同于向量A₂C₂，此时θ₁和θ₂均为0.

坡口底部垂直间隙l₁若为平板对接接头I形或V形坡口，可由式(2)求得^[17]，即

$$ {l_1} = \frac{{{S_{{A_1}{A_2}{C_2}{C_1}}}}}{{{l_{{A_1}{C_1}}}\cos {\theta _1}}} = \frac{{1/2(\left| {\overrightarrow {{A_1}{C_2}} } \right| \cdot \left| {\overrightarrow {{A_2}{C_1}} } \right| \cdot \sin \alpha )}}{{{l_{{A_1}{C_1}}}\cos {\theta _1}}} $$

(2)

式中：$\alpha $为向量A₁C₂和向量A₂C₁的锐夹角；S_{A₁A₂C₂C₁}为四边形A₁A₂C₂C₁的面积；l_A₁C₁、l_A₁C₂和l_A₂C₁分别为线段A₁C₁、A₁C₂和A₂C₁的长度.

坡口底部垂直间隙l₁若为Y形坡口，可由公式(3)求得，即

$$ {l_1} = \frac{{{S_{{A_1}{A_2}{B_2}{B_1}}}}}{{{l_{{A_1}{B_1}}}}} = \frac{{1/2(\left| {\overrightarrow {{A_1}{B_2}} } \right| \cdot \left| {\overrightarrow {{A_2}{B_1}} } \right| \cdot \sin \alpha )}}{{{l_{{A_1}{B_1}}}}} $$

(3)

式中：$\alpha $为向量A₁B₂和向量A₂B₁的锐夹角；S_{A₁A₂B₂B₁}为四边形A₁A₂B₂B₁的面积；l_A₁B₁、l_A₁B₂和l_A₂B₁分别为线段A₁B₁、A₁B₂和A₂B₁的长度.

坡口顶部垂直间隙l₂可由式(4)求得，即

$$ {l_2} = {l_1} + {l_{{C_1}{B_1}}}\sin {\theta _1} + {l_{{C_2}{B_2}}}\sin {\theta _2} $$

(4)

式中：l_C₁B₁、l_C₂B₂分别为线段C₁B₁、C₂B₂的长度.

板厚h₁、h₂可由下式(5)，式(6)求得，即

$$ {h_1} = {l_{{A_1}{B_1}}} + {l_{{B_1}{C_1}}}\cos {\theta _1} $$

(5)

$$ {h_{\text{2}}} = {l_{{A_{\text{2}}}{B_{\text{2}}}}} + {l_{{B_{\text{2}}}{C_{\text{2}}}}}\cos {\theta _{\text{2}}} $$

(6)

2.2 角接接头三种坡口提取

为获取角接接头焊缝坡口处特征尺寸，先提取焊接件数字模型中坡口处特征点及线. 图4为角接接头尺寸标注示意图. θ₁为向量AD与过点A水平向右向量的夹角；θ₂为向量BC与过点B水平向右向量的夹角；θ₃为向量BA与向量BC的夹角；θ₄为向量BA与过点B水平向右向量的夹角；θ₅为向量BC与向量ED的夹角；θ₆为向量AD与向量BC的夹角. 角接接头共有三种坡口，若为不开坡口，则点A、D重合，点B、C水平重合，即l₃ = 0，l_ED= l_EA，θ₁ = 0，θ₂ = 0，θ₆ = 0，θ₃ = θ₄，只需确定A，B，C，E四点坐标；若为单边坡口，则点A，D重合，即l₃ = 0，l_ED= l_EA，θ₁ = 0,只需确定A，B，C，E四点坐标；若为双边坡口，则需确定A，B，C，D，E五点坐标.

图 4 角接接头尺寸示意图

Figure 4. Schematic diagram of fillet welded joint size

下载: 全尺寸图片幻灯片

2.2.1 特征点及线提取

步骤1：根据轮廓线总条数确定角接接头的三种坡口类型. 若为24条，即角接接头不开坡口或单边坡口，转步骤2；若为27条，即角接接头双边坡口，转步骤6.

步骤2：将数组空间Q最小距离的两端点遍历数组空间P，输出重复点所在的轮廓线，得到的两条轮廓线即AE和BC，计算出夹角θ₅^[18]，即

$$ {\theta _5} = \arccos \frac{{\overrightarrow {{\boldsymbol{CB}}} \cdot \overrightarrow {{\boldsymbol{DE}}} }}{{\left| {\overrightarrow {{\boldsymbol{CB}}} } \right| \cdot \left| {\overrightarrow {{\boldsymbol{DE}}} } \right|}} $$

(7)

步骤3：根据夹角θ₅确定角接接头的坡口类型. 若θ₅ = 90°，则为不开坡口，转步骤4；否则为单边坡口，并以数组空间Q最小距离的两端点为起点，所在轮廓线为长度且方向背离各自起点的两向量夹角记为β，转步骤5.

步骤4：根据两个夹角提取出特征点. 如图5所示存在两种情况，数组空间Q最小距离的两端点之一为左端点，或为右端点. 此时，取数组空间Q最小距离的两端点(对应图5a中为A、B，图5b中为A、C)，组成该两点的相交线段AE、AB(AC)和BC所构成的两夹角中，大于90度的为夹角α，小于90度的为夹角θ₃. 根据构成夹角α的两线段AE、AB(AC)可获取点A、E；根据构成夹角θ₃的两线段AB(AC)、BC可获取轮廓线BC.

图 5 不开坡口的两种情形

Figure 5. Two situations of not opening the weld groove. (a) Minimum distance from the left end; (b) minimum distance from the right end

下载: 全尺寸图片幻灯片

步骤5：根据夹角β确定单边坡口的两种情形. 若β > 90°，如图6a所示，数组空间Q最小距离的两端点即为A、B. 遍历数组空间P，点A所在轮廓线的另一端点即为E；点B所在轮廓线的另一端点即为C. 若β < 90°，如图6b所示，数组空间Q最小距离的两端点即为A、C. 同理，可获得点E、B.

图 6 单边坡口的两种情形

Figure 6. Two situations for single-side weld groove. (a) Minimum distance from the left end; (b) minimum distance from the right end

下载: 全尺寸图片幻灯片

步骤6：取数组空间Q最小距离的两端点即为A、B. 遍历数组空间P所有轮廓线的两端点，若重复点所在轮廓线的另一端点为A，则该重复点即为D，并输出该点在数组空间P中所在轮廓线的其余两端点，其中不等于点D和点A的另一端点即为E. 同理，遍历数组空间P，点B所在轮廓线的另一端点即为C.

2.2.2 特征尺寸提取

坡口夹角θ可由式(8)求得,即

$$ \theta = {\theta _6} = {\theta _1} + {\theta _2} = \arccos \frac{{\overrightarrow {{\boldsymbol{AD}}} \cdot \overrightarrow {{\boldsymbol{BC}}} }}{{\left| {\overrightarrow {{\boldsymbol{AD}}} } \right| \cdot \left| {\overrightarrow {{\boldsymbol{BC}}} } \right|}} $$

(8)

坡口底部垂直间隙l₁可由式(9)求得,即

$$ \left\{ \begin{gathered} {\theta _2} = {90^ \circ } - {\theta _5} \\ {\theta _3} = \arccos \frac{{\overrightarrow {{\boldsymbol{BA}}} \cdot \overrightarrow {{\boldsymbol{BC}}} }}{{\left| {\overrightarrow {{\boldsymbol{BA}}} } \right| \cdot \left| {\overrightarrow {{\boldsymbol{BC}}} } \right|}} \\ {\theta _4} = {\theta _3} - {\theta _2} \\ {l_1} = {l_{AB}}\sin {\theta _4} \\ \end{gathered} \right. $$

(9)

式中：l_AB为线段AB的长度.

坡口顶部垂直间隙l₂可由式(10)求得,

$$ \left\{\begin{array}{l} \theta_2=90^{\circ}-\theta_5 \\ \theta_1=\theta_6-\theta_2 \\ l_4=l_{B C} \sin \theta_2 \\ l_3=l_{A D} \sin \theta_1 \\ l_2=l_1+l_3+l_4 \end{array}\right. $$

(10)

式中：l_BC，l_AD分别为线段BC，AD的长度.

板厚h₁、h₂可式(11)式(12)求得,即

$$ {h_1} = {l_{ED}} + {l_3} $$

(11)

$$ {h_2} = {l_{BC}}\cos {\theta _2} $$

(12)

式中： l_ED为线段ED的长度；若轮廓线总条数为24条，则l_ED= l_EA.

2.3 T形接头三种坡口提取

图7为T形接头尺寸标注示意图. θ₁是向量BA与向量CB的夹角；θ₂是向量CD与向量BC的夹角；θ₃为向量EB与过点E水平向右向量的夹角. T形接头共有三种坡口，若为不开坡口，则A与B两点重合，C与D两点重合，即l_AB = 0，l_CD= 0，l_AD=l_BC，θ₁ = 0，θ₂ = 0，只需确定B，C，E，F四点坐标；若为单边坡口，则A与B两点重合(C与D两点重合)，即l_AB= 0，l_AC = l_BC(l_BD= l_BC)，θ₁ = 0 (θ₂ = 0)，只需确定B，C，D(A)，E，F五点坐标；若为双边坡口，则需确定A，B，C，D，E，F六点坐标.

图 7 T形接头尺寸示意图

Figure 7. Schematic diagram of T type welded joint size

下载: 全尺寸图片幻灯片

2.3.1 特征点及线提取

为获取焊缝坡口处特征尺寸，先对焊接件数字模型中的焊缝进行特征点及线的提取. 首先根据轮廓线总条数确定T形接头的三种坡口类型.

若轮廓线总条数为24条，即T形接头不开坡口，取数组空间Q最小距离的两端点即为点B、E. 遍历数组空间P，点B所在轮廓线的另一端点即为点C；点E所在轮廓线的另一端点即为点F.

若轮廓线总条数为27条，即T形接头单边坡口，进行以下处理.

步骤1：遍历数组空间P所有轮廓线的两端点，会有两个重复点，若两重复点与数组空间Q最小距离的两端点不相同，则将两端点遍历数组空间P，获取轮廓线BC(CD)和轮廓线EF，并判断两条轮廓线的位置关系，转步骤2.

步骤2：若两条轮廓线垂直，则数组空间Q最小距离的两端点之一为左端点，如图8a所示. 此时两端点即为点B、E，点B所在轮廓线的另一点即为点C. 若两轮廓线不垂直，或如图8b所示的两个重复点与数组空间Q最小距离的两端点有相同，则将两重复点与数组空间Q最小距离的两端点中，不等于重复点及其所在轮廓线另一端点的点进行两两求距，最小距离的两端点中不等于两重复点的即为点E，另一点即为点C.

图 8 单边坡口的两种情形

Figure 8. Two situations for single-side weld groove. (a) Minimum distance from the left end; (b) minimum distance from the right end

下载: 全尺寸图片幻灯片

步骤3：遍历数组空间P，点E所在轮廓线的另一端点即为F，轮廓线含有点C且垂直于EF的另一端点为点B，不垂直EF的另一端点即为点D.

若轮廓线总条数为30条，即T形接头双边坡口，进行以下处理：

步骤A：遍历数组空间P所有轮廓线的两端点，会有四个重复点，进行两两求距，去除轮廓线本身长度后，最小距离的两端点，即为B、C. 如图9所示，有两种情况分别为，数组空间Q最小距离的两端点之一为左端点，或为右端点. 其中轮廓线两端点含有点B且不等于点C的即为点A，同理，轮廓线两端点含有点C且不等于点B的即为点D.

图 9 双边坡口的两种情形

Figure 9. Two situations of bilateral weld groove. (a) Minimum distance from the left end; (b) minimum distance from the right end

下载: 全尺寸图片幻灯片

步骤B：将点B，C分别与数组空间Q最小距离的两端点中，不等于点A，B，C，D的端点进行求距，最小距离的两点中不等于点B，C的即为点E，另一点则为点B.

步骤C：将组成交点B的三条轮廓线和组成交点E的三轮廓线，构成六条方向背离交点B(E)的向量，其中与向量EB的最小夹角，即为θ₃.

2.3.2 特征尺寸提取

左坡角θ₁和右坡角θ₂可由式(13)求得，即

$$ {\theta _{\text{1}}}{\text{ = }}\arccos \frac{{\overrightarrow {{\boldsymbol{BA}}} \cdot \overrightarrow {{\boldsymbol{CB}}} }}{{\left| {\overrightarrow {{\boldsymbol{BA}}} } \right| \cdot \left| {\overrightarrow {{\boldsymbol{CB}}} } \right|}}{\theta _2}{\text{ = }}\arccos \frac{{\overrightarrow {{\boldsymbol{CD}}} \cdot \overrightarrow {{\boldsymbol{BC}}} }}{{\left| {\overrightarrow {{\boldsymbol{CD}}} } \right| \cdot \left| {\overrightarrow {{\boldsymbol{BC}}} } \right|}} $$

(13)

坡口底部垂直间隙l₁可由式(14)求得，即

$$ {l_1} = {l_{EB}}\sin {\theta _3} $$

(14)

式中：l_EB为线段EB的长度.

坡口顶部垂直间隙l₂，l₃如式(15)式(16)所示，即

$$ {l_2} = {l_1} + {l_{AB}}\sin {\theta _1} $$

(15)

$$ {l_3} = {l_1} + {l_{CD}}\sin {\theta _2} $$

(16)

式中：l_AB，l_CD分别为线段AB，CD的长度.

板厚h₁，h₂可由下式求得，即

$$ {h_1} = {l_{BC}} + {l_{AB}}\cos {\theta _1} + {l_{CD}}\cos {\theta _2} $$

(17)

$$ {h_2} = \sqrt {{{\left( {{x_E} - {x_F}} \right)}^2} + {{\left( {{y_E} - {y_F}} \right)}^2} + {{\left( {{z_E} - {z_F}} \right)}^2}} $$

(18)

式中：l_BC为线段BC的长度；(x_E, y_E, z_E)为点E坐标；(x_F, y_F, z_F) 为点F坐标.

3. 试验测试

3.1 焊缝特征识别测试

为了验证文中船舶焊缝特征参数提取算法的有效性和稳定性，在CPU主频4.6 GHz、内存16 GB的PC平台上，使用Visual Studio 2019编程语言针对图10两舱段上的13处焊缝进行了特征识别测试，并与基于焊缝特征的焊接工件库建模(特征库设计法)^[19]和基于SolidWorks的遍历特征树技术(遍历特征法)^[6]进行了比较. 3种方法的对比测试结果如表1所示.

图 10 焊缝特征识别及参数提取测试

Figure 10. Tests on weld feature recognition and parameter extraction. (a) cabin section; (b) cargo hold assembly section

下载: 全尺寸图片幻灯片

表 1 焊缝特征识别对比结果

Table 1. Comparison results of weld feature recognition

焊缝编号及时间指标	接头测试结果			坡口测试结果
焊缝编号及时间指标	特征库设计法	遍历特征法	文中算法	特征库设计法	遍历特征法	文中算法
1	平接	平接	平接	I形	无法识别	I形
2	平接	平接	平接	V形	无法识别	V形
3	T形	T形	T形	不开	无法识别	不开
4	角接	无法识别	角接	不开	无法识别	不开
5	角接	无法识别	角接	双边	无法识别	双边
6	T形	T形	T形	双边	无法识别	双边
7	平接	平接	平接	I形	无法识别	I形
8	平接	平接	平接	Y形	无法识别	Y形
9	角接	无法识别	角接	单边	无法识别	单边
10	角接	无法识别	角接	双边	无法识别	双边
11	T形	T形	T形	不开	无法识别	不开
12	T形	T形	T形	单边	无法识别	单边
13	搭接	无法识别	搭接	不开	无法识别	不开
平均运行时间t/s	1.118	—	0.895	1.291	—	1.127

下载: 导出CSV

| 显示表格

由表可以看出，遍历特征法仅能识别出平板对接和T形接头两种形式，且无法识别坡口类型，识别效果较差，同时受制于建模平台，应用范围有限. 特征库设计法虽然与文中算法一样，都能准确识别出13处焊缝特征，但在接头和坡口识别的平均运行时间上均高于文中算法，执行效率低，此外，特征库设计法需要先对船舱段上的所有焊缝，进行特征造型设计并建立数据库，使得开发工作量大大增加，且当数据库规模较大在匹配焊缝特征时可能出现耗时较长的情况. 而文中基于STL的焊缝特征识别算法的坡口识别准确率为100%，耗时相比特征库设计法提升了16.06%，平均运行时间更短，执行效率更高，且只需获取焊缝模型的STL文件，前期处理量较小，识别的接头形式和坡口类型多，准确率高.

3.2 焊缝特征参数提取测试

为了进一步验证文中算法的有效性，针对图10两舱段上的13处焊缝进行了焊缝特征参数提取测试，并与基于斜率分析的轮廓线特征点(斜率分析法)^[20]和遍历特征法提取进行了比较. 3种方法的对比测试结果如表2，表3所示.

表 2 焊缝特征参数提取对比结果

Table 2. Comparison results of the extraction of weld feature parameters

焊缝编号	坡口类型	斜率分析法						遍历特征法
焊缝编号	坡口类型	板厚 h₁ /mm	板厚 h₂ /mm	底部间隙 l₁/mm	顶部间隙 l₂/mm	坡角θ₁/ (°)	坡角θ₂/ (°)	板厚 h₁ /mm	板厚 h₂ /mm	焊缝长度 d/mm
1	平接I形	30.15	30.15	5.47	5.47	0	0	30.15	30.15	1350
2	平接V形	29.45	29.45	5.89	26.59	17.50	17.50	29.45	29.45	1401
3	T形不开	21.34	—	6.18	6.18	0	0	21.34	30.51	1545
4	角接不开	—	30	4.34	4.34	0	0	—	—	—
5	角接双边	—	20.57	7.51	26.39	21.00	21.00	—	—	—
6	T形双边	—	—	—	—	—	—	24.35	40.15	1755
7	平接I形	21.5	21.5	4.15	4.15	0	0	21.5	21.5	875
8	平接Y形	27.38	27.38	4.96	26.71	20.00	20.00	27.38	27.38	1003
9	角接单边	—	21.5	5.62	17.13	25.00	0	—	—	—
10	角接双边	—	22.35	4.76	23.18	22.50	22.50	—	—	—
11	T形不开	50.98	—	5.37	25.7	0	0	50.98	24.5	2572
12	T形单边	26.23	—	22.35	4.18	18.22	0	26.23	22.35	718
13	搭接不开	30	21.5	4.5	4.5	0	0	—	—	—

下载: 导出CSV

| 显示表格

由表可以看出，遍历特征法只能针对平板对接和T形接头进行参数提取，所能提取的对象有限，实用性较低，同时，只能提取坡口处两焊接件的板厚和焊缝长度，而与焊接工艺参数密切相关的焊缝间隙、坡口夹角等信息无法获取，参数提取效果差，不利于焊接工艺参数的快速、准确选择. 斜率法相比遍历特征法，提取的接头形式和坡口类型较多，获取的焊缝特征参数信息较多，但对角接、T形接头提取时，无法获取部分焊接件的板厚信息，且无法对T形接头双边坡口的焊缝特征参数进行提取，提取效果较差. 而文中算法不仅能针对平板对接、角接、T形和搭接接头4类接头及10种坡口进行提取，而且获取的焊缝特征参数最多，包括两焊接件板厚、底部和顶部坡口间隙和坡口夹角等信息，为焊缝工艺参数的快速、自适应选择提供了充分的参考依据. 这主要得益于基于STL的焊缝特征识别，能准确识别出焊缝模型STL文件中的接头形式和坡口类型，且执行效率高，在焊缝特征识别的基础上，对3类接头9种坡口进行共性分析，进而建立3种焊缝特征参数数学模型完成坡口信息的快速精确提取，验证了文中算法的有效性.

表 3 文中算法焊缝特征参数提取结果

Table 3. Results of weld characteristic parameters calculated in this paper

焊缝编号	坡口类型	板厚h₁ /mm	板厚h₂ /mm	底部间隙l₁/mm	顶部间隙l₂/mm	坡角θ₁/(°)	坡角θ₂/(°)
1	平接I形	30.15	30.15	5.47	5.47	0	0
2	平接V形	29.45	29.45	5.89	26.59	17.50	17.50
3	T形不开	21.34	30.5	6.18	6.18	0	0
4	角接不开	18.39	30	4.34	4.34	0	0
5	角接双边	19.45	20.57	7.51	26.39	21.00	21.00
6	T形双边	21.34	40.58	8.76	23.81	70.25	69.55
7	平接I形	21.5	21.5	4.15	4.15	0	0
8	平接Y形	27.38	27.38	4.96	26.71	20.00	20.00
9	角接单边	29	21.5	5.62	17.13	25.00	0
10	角接双边	50.98	22.35	4.76	23.18	22.50	22.50
11	T形不开	50.98	24.5	5.37	25.7	0	0
12	T形单边	26.23	22.35	4.18	18.22	0	61.75
13	搭接不开	30	21.5	4.5	4.5	0	0

下载: 导出CSV

| 显示表格

4. 结论

为了解决目前焊缝特征识别单一、焊缝特征参数信息提取有限的问题，文中提出了基于空间位置和轮廓线距离的船舶焊缝特征参数提取算法. 根据测试结果可以得出以下结论：

(1) 基于轮廓线总条数和最小轮廓线距离的焊缝特征识别算法，能准确识别4类接头形式和10种坡口类型，识别种类多、耗时少、准确率高，大大提高了算法的实用性、广泛性.

(2) 基于模型轮廓线总条数，能快速区分出三类接头中的每种坡口类型；根据最小轮廓线距离的两端点，能快速提取出数学模型所需的特征点、线及关键夹角，有利于焊缝特征参数的提取.

(3) 基于三类接头形式所建立的三种焊缝特征参数数学模型，能适应目前常见的坡口类型，如I形、V形和Y形、不开、单边和双边等坡口，所能提取的对象较广，获取的特征参数信息量多，为焊接工艺参数的选择提供了充分依据.

测试结果表明文中算法具有焊缝特征识别种类多且精确、特征参数提取信息量多且高效的优势，为船舶焊接机器人自适应、快速准确选择焊接工艺参数，提高焊接质量和效率提供了保障.

[1]

王加友，杨峰，韩伟. 摇动电弧窄间隙熔化极气体保护焊接方法及焊炬:中国, 200810236274.5[P], 2008-11-19.[2] 任志鹏. 电弧摆动式NG-GMAW焊枪及其工艺性能[D]. 上海：上海交通大学，2013.[3] 朱杰. 摇动电弧窄间隙焊接系统的单机控制研究[D]. 镇江：江苏科技大学, 2009.[4] Wang J Y, Zhu J, Zhang C,et al. Development of swing arc narrow gap vertical welding process[J]. ISIJ International, 2015, 55(5): 1076-1082.[5] Wang J Y, Zhu J, Fu P,et al. A swing arc system for narrow gap GMA welding[J]. ISIJ International, 2012, 52(1):110-114.[6] 姬书得，方洪渊，刘雪松，等. 基于串状热源的手工摆动焊接压力场的数值模拟[J]. 焊接学报，2005, 26(5):46-48.Ji Shude, Fang Hongyuan, Liu Xuesong,et al. Numerical simulation of stress field of manual swing welding on basis of cluster heat source[J]. Transactions of the China Welding Institution, 2005, 26(5): 46-48.[7] 张华军，张广军，蔡春波，等. 摆动焊接动态过程温度场数值模拟[J]. 焊接学报，2008, 29(2): 65-68.Zhang Huajun, Zhang Guangjun, Cai Chunbo,et al. Numerical simulation of temperature field of dynamic welding process with weaving[J]. Transactions of the China Welding Institution, 2008, 29(2): 65-68.[8] 兰虎，张华军，陈阿静，等. 窄间隙GMAW立焊动态过程模拟及热物理特性[J]. 焊接学报，2015，36(7): 77-82.Lan Hu, Zhang Huajun, Chen Ajing,et al. Numerical simulation of dynamic process and thermal physical properties of narrow gap MAG vertical welding[J]. Transactions of the China Welding Institution, 2015, 36(7): 77-82.[9] 方臣富，陈志伟，胥国祥，等. 缆式焊丝CO2气体保护焊工艺研究[J]. 金属学报，2012,48(11): 1299-1305.Fang Chenfu, Chen Zhiwei, Xu Guoxiang,et al. Study on the process of CTWW CO2gas shielded welding[J]. Acta Metallurgica Sinca, 2012, 48(11): 1299-1305.

期刊类型引用(2)

1.	张敏，王新宝，王浩军，马克，朱子越，张志强. 激光熔覆TC4/Inconel 625/316L不锈钢梯度材料组织与性能. 焊接学报. 2023(07): 16-23+129-130 . 本站查看
2.	王永东，宫书林，汤明日，宋闽. 激光熔覆工艺对高熵合金组织与性能影响. 焊接学报. 2023(08): 116-122+136 . 本站查看