Microstructure transformation and variant selection of SLM and SLM + HT formed TA17 alloy
-
摘要: 为了系统认识选区激光熔化(selective laser melting, SLM)成形TA17合金微观组织特征,对SLM和SLM + 热处理(heat treatment, HT)两种状态下TA17合金的微观组织进行表征和分析. 结果表明,在极快的冷却条件下,SLM成形的TA17钛合金形成了以细小针状α相为主的微观组织,其织构类型主要为{0 0 0 1} < 1 0 −1 0 >、{1 0 −1 0} < 1 1 −2 0 >及{1 1 −2 0} < 1 0 −1 0 >3类板织构,且变体具有明显选择性特征,倾向于析出60°/[1 1 −2 0]和60.83°/[−1.377 −1 2.377 0.359]两类变体. 在950 ℃/2 h + 空冷后,TA17合金形成双态α相组织,组织形貌为典型的等轴状 + 板条状,织构散漫程度增加,强度显著降低,变体选择发生变化,60°/[1 1 −2 0]和60.83°/[−1.377 −1 2.377 0.359]两类变体减少,转变为63.26°/[−10 5 5 −3]类变体.Abstract: In order to systematically investigate the microstructure transformation mechanism of selective laser melting(SLM) formed TA17 alloy, the microstructures of TA17 alloy were characterized and analyzed under SLM and SLM + heat treatment(HT). The results show that under extremely fast cooling conditions, the microstructure of TA17 is dominated by fine needle-shape α phase, with plate textures {0 0 0 1} < 1 0 −1 0 >, {1 0 −1 0} < 1 1 −2 0 > and {1 1 −2 0} < 1 0 −1 0 > as main types; and two types of variants, namely the 60°/[1 1 −2 0] and 60.83°/[−1.377 −1 2.377 0.359], are inclined to be generated. While after a treatment of 950 °C/2 h + air cooling, a dual state α structure is formed on the TA17 alloy with typical equiaxed grains and a plate-like morphology. With the increase in the degree of texture dispersion, the strength decreased significantly. There was also a change in variant selection, as the proportion of 60°/[1 1 −2 0] and 60.83°/[−1.377 −1 2.377 0.359] variants decreased significantly, 63.26°/[−10 5 5 −3] become the major variant.
-
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.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形接头.
当轮廓线条为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为平板对接接头尺寸标注示意图,θ1为向量A1B1和向量B1C1的夹角;θ2为向量A2B2和向量B2C2的夹角. 平板对接接头共有三种坡口,若为I形坡口,则点B1,C1两点重合,点B2,C2重合,即lC1B1 = 0,lC2B2 = 0,只需确定A1,A2,C1,C2 4点坐标;若为V形坡口,则点B1,A1重合,点B2,A2重合,即lA1B1 = 0,lA2B2 = 0,只需确定A1,A2,C1,C2四点坐标;若为Y形坡口,则需确定A1,A2,B1,B2,C1,C2六点坐标.
2.1.1 特征点及线提取
根据轮廓线总条数确定平板对接接头的3种坡口类型. 若为24条,即平板对接接头为I形或V形坡口,则取数组空间Q最小距离所对应的两端点,即为A1,A2. 遍历数组空间P,点A1所在轮廓线的另一端点即为C1,点A2所在轮廓线的另一端点即为C2.
轮廓线总条数若为30条,即平板对接接头Y形坡口,则遍历数组空间P所有轮廓线的两端点,将重复点两两求距,最小距离所对应的两点,即为B1,B2. 分别将组成交点B1,B2的轮廓线其余端点两两求距,最小距离所对应的两点即为A1,A2,剩余两端点则为C1,C2.
2.1.2 特征尺寸提取
坡口夹角θ根据所提取坡口处特征点及线,计算出夹角θ1,θ2,由式(1)可求得坡口夹角θ,即
$$ \theta = {\theta _1} + {\theta _2} $$ (1) 式中:θ1为左坡角;θ2为右坡角.若轮廓线总条数为24条,则向量B1C1等同于向量A1C1,向量A2B2等同于向量A2C2,此时θ1和θ2均为0.
坡口底部垂直间隙l1若为平板对接接头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 $ 为向量A1C2和向量A2C1的锐夹角;SA1A2C2C1为四边形A1A2C2C1的面积;lA1C1、lA1C2和lA2C1分别为线段A1C1、A1C2和A2C1的长度.坡口底部垂直间隙l1若为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 $ 为向量A1B2和向量A2B1的锐夹角;SA1A2B2B1为四边形A1A2B2B1的面积;lA1B1、lA1B2和lA2B1分别为线段A1B1、A1B2和A2B1的长度.坡口顶部垂直间隙l2可由式(4)求得,即
$$ {l_2} = {l_1} + {l_{{C_1}{B_1}}}\sin {\theta _1} + {l_{{C_2}{B_2}}}\sin {\theta _2} $$ (4) 式中:lC1B1、lC2B2分别为线段C1B1、C2B2的长度.
板厚h1、h2可由下式(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为角接接头尺寸标注示意图. θ1为向量AD与过点A水平向右向量的夹角;θ2为向量BC与过点B水平向右向量的夹角;θ3为向量BA与向量BC的夹角;θ4为向量BA与过点B水平向右向量的夹角;θ5为向量BC与向量ED的夹角;θ6为向量AD与向量BC的夹角. 角接接头共有三种坡口,若为不开坡口,则点A、D重合,点B、C水平重合,即l3 = 0,lED= lEA,θ1 = 0,θ2 = 0,θ6 = 0,θ3 = θ4,只需确定A,B,C,E四点坐标;若为单边坡口,则点A,D重合,即l3 = 0,lED = lEA,θ1 = 0,只需确定A,B,C,E四点坐标;若为双边坡口,则需确定A,B,C,D,E五点坐标.
2.2.1 特征点及线提取
步骤1:根据轮廓线总条数确定角接接头的三种坡口类型. 若为24条,即角接接头不开坡口或单边坡口,转步骤2;若为27条,即角接接头双边坡口,转步骤6.
步骤2:将数组空间Q最小距离的两端点遍历数组空间P,输出重复点所在的轮廓线,得到的两条轮廓线即AE和BC,计算出夹角θ5[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:根据夹角θ5确定角接接头的坡口类型. 若θ5 = 90°,则为不开坡口,转步骤4;否则为单边坡口,并以数组空间Q最小距离的两端点为起点,所在轮廓线为长度且方向背离各自起点的两向量夹角记为β,转步骤5.
步骤4:根据两个夹角提取出特征点. 如图5所示存在两种情况,数组空间Q最小距离的两端点之一为左端点,或为右端点. 此时,取数组空间Q最小距离的两端点(对应图5a中为A、B,图5b中为A、C),组成该两点的相交线段AE、AB(AC)和BC所构成的两夹角中,大于90度的为夹角α,小于90度的为夹角θ3. 根据构成夹角α的两线段AE、AB(AC)可获取点A、E;根据构成夹角θ3的两线段AB(AC)、BC可获取轮廓线BC.
步骤5:根据夹角β确定单边坡口的两种情形. 若β > 90°,如图6a所示,数组空间Q最小距离的两端点即为A、B. 遍历数组空间P,点A所在轮廓线的另一端点即为E;点B所在轮廓线的另一端点即为C. 若β < 90°,如图6b所示,数组空间Q最小距离的两端点即为A、C. 同理,可获得点E、B.
步骤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) 坡口底部垂直间隙l1可由式(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) 式中:lAB为线段AB的长度.
坡口顶部垂直间隙l2可由式(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) 式中:lBC,lAD分别为线段BC,AD的长度.
板厚h1、h2可式(11)式(12)求得,即
$$ {h_1} = {l_{ED}} + {l_3} $$ (11) $$ {h_2} = {l_{BC}}\cos {\theta _2} $$ (12) 式中: lED为线段ED的长度;若轮廓线总条数为24条,则lED = lEA.
2.3 T形接头三种坡口提取
图7为T形接头尺寸标注示意图. θ1是向量BA与向量CB的夹角;θ2是向量CD与向量BC的夹角;θ3为向量EB与过点E水平向右向量的夹角. T形接头共有三种坡口,若为不开坡口,则A与B两点重合,C与D两点重合,即lAB = 0,lCD = 0,lAD =lBC,θ1 = 0,θ2 = 0,只需确定B,C,E,F四点坐标;若为单边坡口,则A与B两点重合(C与D两点重合),即lAB = 0,lAC = lBC(lBD = lBC),θ1 = 0 (θ2 = 0),只需确定B,C,D(A),E,F五点坐标;若为双边坡口,则需确定A,B,C,D,E,F六点坐标.
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.
步骤3:遍历数组空间P,点E所在轮廓线的另一端点即为F,轮廓线含有点C且垂直于EF的另一端点为点B,不垂直EF的另一端点即为点D.
若轮廓线总条数为30条,即T形接头双边坡口,进行以下处理:
步骤A:遍历数组空间P所有轮廓线的两端点,会有四个重复点,进行两两求距,去除轮廓线本身长度后,最小距离的两端点,即为B、C. 如图9所示,有两种情况分别为,数组空间Q最小距离的两端点之一为左端点,或为右端点. 其中轮廓线两端点含有点B且不等于点C的即为点A,同理,轮廓线两端点含有点C且不等于点B的即为点D.
步骤B:将点B,C分别与数组空间Q最小距离的两端点中,不等于点A,B,C,D的端点进行求距,最小距离的两点中不等于点B,C的即为点E,另一点则为点B.
步骤C:将组成交点B的三条轮廓线和组成交点E的三轮廓线,构成六条方向背离交点B(E)的向量,其中与向量EB的最小夹角,即为θ3.
2.3.2 特征尺寸提取
左坡角θ1和右坡角θ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) 坡口底部垂直间隙l1可由式(14)求得,即
$$ {l_1} = {l_{EB}}\sin {\theta _3} $$ (14) 式中:lEB为线段EB的长度.
坡口顶部垂直间隙l2,l3如式(15)式(16)所示,即
$$ {l_2} = {l_1} + {l_{AB}}\sin {\theta _1} $$ (15) $$ {l_3} = {l_1} + {l_{CD}}\sin {\theta _2} $$ (16) 式中:lAB,lCD分别为线段AB,CD的长度.
板厚h1,h2可由下式求得,即
$$ {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) 式中:lBC为线段BC的长度;(xE, yE, zE)为点E坐标;(xF, yF, zF) 为点F坐标.
3. 试验测试
3.1 焊缝特征识别测试
为了验证文中船舶焊缝特征参数提取算法的有效性和稳定性,在CPU主频4.6 GHz、内存16 GB的PC平台上,使用Visual Studio 2019编程语言针对图10两舱段上的13处焊缝进行了特征识别测试,并与基于焊缝特征的焊接工件库建模(特征库设计法) [19]和基于SolidWorks的遍历特征树技术(遍历特征法)[6]进行了比较. 3种方法的对比测试结果如表1所示.
表 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 由表可以看出,遍历特征法仅能识别出平板对接和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焊缝
编号坡口类型 斜率分析法 遍历特征法 板厚
h1 /mm板厚
h2 /mm底部间隙
l1/mm顶部间隙
l2/mm坡角θ1/
(°)坡角θ2/
(°)板厚
h1 /mm板厚
h2 /mm焊缝长度
d/mm1 平接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 — — — 由表可以看出,遍历特征法只能针对平板对接和T形接头进行参数提取,所能提取的对象有限,实用性较低,同时,只能提取坡口处两焊接件的板厚和焊缝长度,而与焊接工艺参数密切相关的焊缝间隙、坡口夹角等信息无法获取,参数提取效果差,不利于焊接工艺参数的快速、准确选择. 斜率法相比遍历特征法,提取的接头形式和坡口类型较多,获取的焊缝特征参数信息较多,但对角接、T形接头提取时,无法获取部分焊接件的板厚信息,且无法对T形接头双边坡口的焊缝特征参数进行提取,提取效果较差. 而文中算法不仅能针对平板对接、角接、T形和搭接接头4类接头及10种坡口进行提取,而且获取的焊缝特征参数最多,包括两焊接件板厚、底部和顶部坡口间隙和坡口夹角等信息,为焊缝工艺参数的快速、自适应选择提供了充分的参考依据. 这主要得益于基于STL的焊缝特征识别,能准确识别出焊缝模型STL文件中的接头形式和坡口类型,且执行效率高,在焊缝特征识别的基础上,对3类接头9种坡口进行共性分析,进而建立3种焊缝特征参数数学模型完成坡口信息的快速精确提取,验证了文中算法的有效性.
表 3 文中算法焊缝特征参数提取结果Table 3. Results of weld characteristic parameters calculated in this paper焊缝编号 坡口类型 板厚h1 /mm 板厚h2 /mm 底部间隙l1/mm 顶部间隙l2/mm 坡角θ1/(°) 坡角θ2/(°) 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 4. 结论
为了解决目前焊缝特征识别单一、焊缝特征参数信息提取有限的问题,文中提出了基于空间位置和轮廓线距离的船舶焊缝特征参数提取算法. 根据测试结果可以得出以下结论:
(1) 基于轮廓线总条数和最小轮廓线距离的焊缝特征识别算法,能准确识别4类接头形式和10种坡口类型,识别种类多、耗时少、准确率高,大大提高了算法的实用性、广泛性.
(2) 基于模型轮廓线总条数,能快速区分出三类接头中的每种坡口类型;根据最小轮廓线距离的两端点,能快速提取出数学模型所需的特征点、线及关键夹角,有利于焊缝特征参数的提取.
(3) 基于三类接头形式所建立的三种焊缝特征参数数学模型,能适应目前常见的坡口类型,如I形、V形和Y形、不开、单边和双边等坡口,所能提取的对象较广,获取的特征参数信息量多,为焊接工艺参数的选择提供了充分依据.
测试结果表明文中算法具有焊缝特征识别种类多且精确、特征参数提取信息量多且高效的优势,为船舶焊接机器人自适应、快速准确选择焊接工艺参数,提高焊接质量和效率提供了保障.
-
表 1 SLM成形用TA17合金粉末成分(质量分数,%)
Table 1 Compositions of TA17 powder for SLM
Al V C O N H Zr Fe 其它 Ti 3.8 ~ 5.0 1.4 ~ 2.5 ≤0.04 ≤0.15 ≤0.04 ≤0.006 ≤0.30 ≤0.25 ≤0.30 余量 表 2 SLM制造工艺参数及粉末粒度范围
Table 2 SLM process and particle size range of TA17 powder
激光功率P/W 打印速率v/(mm·s−1) 层间距δ/mm 粉末粒度d/μm 240 1200 0.14 16 ~ 51 表 3 与β相具有Burgers取向关系的12个α相变体[17,19]
Table 3 12 α variants with Burgers orientation relationship to β phase
序号 取向关系 与A的旋转角θ/(°) 与A的旋转轴 A (0 0 0 1)α//(1 −1 0)β, [1 1 −2 0]α//[1 1 1]β — — B (0 0 0 1)α//(1 0 −1)β, [1 1 −2 0]α//[1 1 1]β 60 [1 1 −2 0] C (0 0 0 1)α//(0 1 −1)β, [1 1 −2 0]α//[1 1 1]β 60 [1 1 −2 0] D (0 0 0 1)α//(1 1 0)β, [1 1 −2 0]α//[1 1 1]β 90 [1 −2.38 1.38 0] E (0 0 0 1)α//(1 0 1)β, [1 1 −2 0]α//[1 1 1]β 63.26 [−10 5 5 −3] F (0 0 0 1)α//(0 1 −1)β, [1 1 −2 0]α//[1 1 1]β 60.83 [−1.377 −1 2.377 0.359] G (0 0 0 1)α//(1 1 0)β, [1 1 −2 0]α//[1 1 1]β 90 [1 −2.38 1.38 0] H (0 0 0 1)α//(1 0 −1)β, [1 1 −2 0]α//[1 1 1]β 60.83 [−1.377 −1 2.377 0.359] I (0 0 0 1)α//(0 1 1)β, [1 1 −2 0]α//[1 1 1]β 63.26 [−10 5 5 −3] J (0 0 0 1)α//(1 −1 0)β, [1 1 −2 0]α//[1 1 1]β 10.53 [0 0 0 1] K (0 0 0 1)α//(1 0 1)β, [1 1 −2 0]α//[1 1 1]β 60.83 [−1.377 −1 2.377 0.359] L (0 0 0 1)α//(0 1 1)β, [1 1 −2 0]α//[1 1 1]β 60.83 [−1.377 −1 2.377 0.359] -
[1] 张喜燕, 赵永庆, 白晨光. 钛合金及应用[M]. 北京: 化工工业出版社, 2005. Zhang Xiyan, Zhao Yongqing, Bai Chenguang. Titanium alloys and applications[M]. Beijing: Chemical Industry Press, 2005.
[2] Short A B. Gas tungsten arc welding of α + β titanium alloys: a review[J]. Metal Science Journal, 2009, 25(3): 309 − 324.
[3] 金和喜, 魏克湘, 李建明, 等. 航空用钛合金研究进展[J]. 中国有色金属学报, 2015, 25(2): 280 − 292. Jin Hexi, Wei Kexiang, Li Jianming, et al. Research development of titanium alloy in aerospace industry[J]. The Chinses Journal of Nonferrous Metals, 2015, 25(2): 280 − 292.
[4] 孟圣昊, 司昌健, 任逸群, 等. 中厚板TC4钛合金真空环境激光焊接特性[J]. 焊接学报, 2021, 42(8): 40 − 48. Meng Shenghao, Si Changjian, Ren Yiqun et al. Study on laser welding characteristics of thick wall TC4 titanium alloy in vacuum environment[J]. Transactions of the China Welding Institution, 2021, 42(8): 40 − 48.
[5] 习文顺, 任鑫焱, 张金元, 等. 高速列车TC4钛合金焊接构架强度及寿命评估[J]. 焊接学报, 2022, 43(5): 29 − 35. Xi Wenshun, Ren Xinyan, Zhang Jinyuan, et al. Strength and life assessment of TC4 titanium alloy welded frame for high-speed railway vehicles[J]. Transactions of the China Welding Institution, 2022, 43(5): 29 − 35.
[6] 张义文. 增材制造用钛合金粉末和高温合金粉末的生产[J]. 粉末冶金工业, 2016, 26(1): 76. Zhang Yiwen. Production of titanium alloy powder and high-temperature alloy powder for additive manufacturing[J]. Powder Metallurgy Industry, 2016, 26(1): 76.
[7] 魏明炜, 陈岁元, 郭快快, 等. EIGA法制备激光3D打印用TA15钛合金粉末[J]. 材料导报, 2017, 31(12): 64 − 67,78. doi: 10.11896/j.issn.1005-023X.2017.012.014 Wei Mingwei, Chen Suiyuan, Guo Kuaikuai, et al. Preparation of TA15 titanium alloy powder by EIGA for laser 3D printing[J]. Materials Reports, 2017, 31(12): 64 − 67,78. doi: 10.11896/j.issn.1005-023X.2017.012.014
[8] 蔡笑宇, 董博伦, 王俊哲, 等. 热处理对GTA增材制造TiAl合金组织与性能的调控[J]. 焊接学报, 2022, 43(3): 7 − 12. Cai Xiaoyu, Dong Bolun, Wang Junzhe, et al. Control of the microstructure and mechanical properties of GTA-based wire arc additive manufactured TiAl alloys using post heat treatment[J]. Transactions of the China Welding Institution, 2022, 43(3): 7 − 12.
[9] 谢波. EIGA雾化法制备激光3D打印用TC4合金粉末工艺研究[J]. 钢铁钒钛, 2019, 40(3): 7 − 12. Xie Bo. Preparation of TC4 alloy powders used for laser 3D printing via EIGA method[J]. Iron Steel Vanadium Titanium, 2019, 40(3): 7 − 12.
[10] 郭快快, 刘常升, 陈岁元, 等. 功率对EIGA制备3D打印用TC4合金粉末特性的影响[J]. 材料科学与工艺, 2017, 25(1): 16 − 22. Guo Kuaikuai, Liu Changsheng, Chen Suiyuan, et al. Effect of EIGA power parameter on the characteristics of TC4 alloy powder for 3D printing[J]. Materials Science & Technology, 2017, 25(1): 16 − 22.
[11] Dutta B, Froes F H. The additive manufacturing (AM) titanium alloys[J]. Metal Powder Report, 2017, 72(2): 96 − 106. doi: 10.1016/j.mprp.2016.12.062
[12] 付超, 叶义海, 陶涛, 等. EIGA法制备TA17合金粉末的激光增材适应性研究[J]. 材料科学与工艺, 2021, 29(5): 91 − 96. doi: 10.11951/j.issn.1005-0299.20200324 Fu Chao, Ye Yihai, Tao Tao, et al. Study on laser additive adaptability of TA17 titanium alloy powder prepared by EIGA method[J]. Materials Science & Technology, 2021, 29(5): 91 − 96. doi: 10.11951/j.issn.1005-0299.20200324
[13] 伍建文, 芦丽莉, 陶涛, 等. 退火热处理对TA17钛合金激光选区熔化成形件的力学性能影响研究[J]. 电焊机, 2021, 51(1): 111 − 114. Wu Jianwen, Lu Lili, Tao Tao, et al. Study on the effect of annealing heat treatment on the mechanical properties of selective laser melting of TA17 titanium alloy[J]. Electric Welding Machine, 2021, 51(1): 111 − 114.
[14] Zhuo L, Jing L, Zhu Y, et al. Variant selection in laser melting deposited α + β titanium alloy[J]. Journal of Alloys & Compounds, 2016, 661: 126 − 135.
[15] 魏子淦, 杨平. TA10钛合金相变过程中的变体选择规律[C]// 第十六届中国体视学与图像分析学术会议——交叉、融合、创新. 海口, 中国, 2019: 058018. Wei Zigan, Yang Ping. Variant selection in the phase transformation process of TA10 titanium alloy[C]//The 16th Conference on Stereology and Image Analysis-Intersection, Fusion, and Innovation. Haikou, China, 2019: 058018.
[16] Liu Q, Qiu C. Variant selection of α precipitation in a beta titanium alloy during selective laser melting and its influence on mechanical properties[J]. Materials Science & Engineering: A, 2020, 784: 139336. doi: 10.1016/j.msea.2020.139336
[17] Lei L, Zhao Q Y, Cong W C, et al. Variant selection, coarsening behavior of α phase and associated tensile properties in an α + β titanium alloy[J]. Journal of Materials Science & Technology, 2021, 99: 101 − 113.
[18] 左玉婷, 王书明, 李聪, 等. 激光增材制造钛合金织构跨尺度分析[J]. 稀有金属材料与工程, 2021, 50(4): 1365 − 1370. Zuo Yuting, Wang Shuming, Li Cong, et al. Multi-scale texture analysis of titanium alloy made by laser additive manufacturing[J]. Rare Metal Materials and Engineering, 2021, 50(4): 1365 − 1370.
[19] 郑国明, 李磊, 毛小南, 等. 钛合金BCC-HCP相变的变体选择及其对晶体取向的影响[J]. 材料导报, 2019, 33(17): 2910 − 2917. doi: 10.11896/cldb.18070038 Zheng Guoming, Li Lei, Mao Xiaonan, et al. Variant selection during titanium alloy BCC-HCP phase transformation and its effect on crystal orientation[J]. Materials Reports, 2019, 33(17): 2910 − 2917. doi: 10.11896/cldb.18070038
[20] Ahmed T, Rack H J. Phase transformations during cooling in α + β titanium alloys[J]. Materials Science & Engineering: A, 1998, 243(1): 206 − 211.
[21] 谢英杰, 付文杰, 王蕊宁, 等. 热处理对 TA15钛合金中厚板材组织及力学性能的影响[J]. 钛工业进展, 2013(6): 26 − 29. Xie Yingjie, Fu Wenjie, Wang Ruining, et al. Effect of heat treatment on microstructure and mechanical properties of TA15 plates[J]. Titanium Industry Progress, 2013(6): 26 − 29.
-
期刊类型引用(2)
1. 张敏,王新宝,王浩军,马克,朱子越,张志强. 激光熔覆TC4/Inconel 625/316L不锈钢梯度材料组织与性能. 焊接学报. 2023(07): 16-23+129-130 . 本站查看
2. 王永东,宫书林,汤明日,宋闽. 激光熔覆工艺对高熵合金组织与性能影响. 焊接学报. 2023(08): 116-122+136 . 本站查看
其他类型引用(1)