1、第3 7卷第6期2 0 2 2年1 2月安 徽 工 程 大 学 学 报J o u r n a l o fA n h u iP o l y t e c h n i cU n i v e r s i t yV o l.3 7N o.6D e c.2 0 2 2文章编号:1 6 7 2-2 4 7 7(2 0 2 2)0 6-0 0 2 5-0 7收稿日期:2 0 2 2-0 4-2 6 基金项目:河南省重点研发与推广专项科技攻关基金资助项目(2 1 2 1 0 2 2 1 0 3 7 0);安徽省高等学校自然科学研究重点基金资助项目(K J 2 0 2 1 A 1 3 2 4);安徽省教育厅教学研
2、究基金资助项目(2 0 2 0 J Y XM 0 6 4 3)作者简介:郭文举(1 9 8 1-),男,河南鹿邑人,讲师,硕士。通信作者:郭 楠(1 9 8 4-),女,河南扶沟人,高级工程师,博士。热输入与薄板材T I G焊接面内外变形关系研究郭文举1,2,郭 楠3*,马喜强3,杨 芳3(1.芜湖市经济和信息化局 投资与技术创新科,安徽 芜湖 2 4 1 0 6 0;2.芜湖职业技术学院 智能制造学院,安徽 芜湖,2 4 1 0 0 0;3.河南科技大学 机电工程学院,河南 洛阳 4 7 1 0 4 9)摘要:热输入对薄板材T I G焊接面外变形的影响较大,本文采用基于数字图像相关的三维非接
3、触检测试验与热弹塑性有限元预测相结合的方法对不同热输入下的薄板材焊接变形进行研究,研究在不同焊接热输入下平板堆焊的面内、面外变形规律。结果表明,热输入需要达到一定的门槛值才能产生焊缝收缩力;随着焊接热输入的增大,焊缝处的面外变形也随之增大;热输入与薄板材焊接最大横向塑性应变值呈线性关系。关 键 词:热输入;焊接;面外变形;焊缝中图分类号:T G 4 0 4 文献标志码:A焊接热输入是焊接工艺的一项重要技术指标,对焊接接头1-2、残余变形3-4等的影响非常显著。为了获得较小的焊接残余变形,工程中通常取较小的热输入量,虽然减小了变形量,但同时影响了生产效率,而过大的热输入会产生焊接缺陷,影响焊接接
4、头强度。在研究焊接失稳变形时,将焊接薄板材焊缝收缩变形的产生设想成在焊缝位置存在一个收缩力,这一收缩力作用在原始无应力的薄板材上,使薄板材产生压缩变形。根据弹性理论,可以近似地由焊缝及其附近区域的焊缝收缩力决定。因此,基于弹性理论的固有应变法认为5,焊接失稳变形是由纵向固有收缩力产生的,纵向固有收缩是由薄板材周边材料的限制引起的,阻止了焊缝的自由收缩,因此,纵向固有变形在焊缝区域出现拉力,也就是焊缝收缩力。焊接循环完成后,在远离焊缝处的纵向压应力在横截面上的积分值就等于固有收缩力。因此,国内外众多学者研究焊缝固有收缩力与热输入之间的关系。郭楠等6采用基于热弹塑性有限元法对薄板材T I G堆焊屈
5、曲变形进行预测并获得较高精度。周佳惠7提出无量纲热输入参数,该参数考虑不同材料和焊接条件的热输入,建立了焊接工艺参数与无量纲热输入参数之间的联系。邵晴等3研究焊接热输入对高速动车组转向架侧梁焊接变形的影响,基于热机耦合仿真分析方法和具有交互作用的正交试验设计方法建立代理模型,获得了最优的焊接工艺参数组合。曹宇等8提出基于热弹塑性有限元法与固有应变法对薄板焊接过程进行数值仿真分析,以船用薄板为研究对象验证了该方法预报焊接变形的合理性与准确性。Wh i t e等9提出了焊缝固有收缩力与净热输入之间的关系。S a t o h等1 0提出了焊缝固有收缩力和净热输入之间关系的试验经验公式。T e r a
6、 s a k i等1 1研究了焊缝固有收缩力和净热输入的关系与薄板材的厚度有关系。L u o等1 2采用热弹塑性有限元法分析研究了薄板材焊缝固有收缩力与净热输入的关系,认为热输入必须超过一定的门槛值才会出现焊缝固有收缩力。本文通过改变焊接电流来研究不同热输入下焊接及冷却过程中薄板材的位移和应变分布,提出热输入与纵向收缩力的数学关系模型,阐明热输入对薄板材焊接面内以及面外变形的影响。1 数字图像相关技术数字图像相关法是对散斑图像和变形后的散斑图像进行相关性计算以获取待匹配C点位移坐标的变化,如图1所示。在参考图像中的选定区域划分网格,取一个以待匹配点C为中心的矩形子网格作为参考子图像。在变形图像
7、中,根据已定义的相关系数以及一定的搜索方法与参考子图像的像素灰度进行相关计算,以确定与参考子图像相关性最大的目标子图像,并计算出变形图像中心点C的坐标,那么变形图像中的点C就是参考图像中点C的对应点。图1 数字图像相关原理图数字图像相关法的关键是通过求取参考子图像与目标子图像的相关系数极值来完成图像匹配,进而得到相应的位移计算。假设参考子图像上任意一点P的坐标为(x,y),变形后的对应点P的图像坐标为(x,y),如果物体在只发生刚体运动的情况下,两点之间的映射关系为零阶映射函数:x=x+u,y=y+v,(1)式中,u为x方向的位移分量/像素;v为y方向的位移分量/像素。当物体受力变形时,除了会
8、产生简单的平移和转动,还可能会产生伸缩、扭曲等变形,因此,对于物体表面测量点坐标变化的计算,除了要考虑刚体位移本身,还需要考虑伸缩、扭转等变形后的坐标变化,也就是需要映射函数的引入导数项,一阶映射函数如下:x=x+u+uxx+uyy,y=x+v+vxx+vyy,(2)式中,u为子图像中心点C变形后在x方向的位移;v为子图像中心点C变形后在y方向的位移;ux、uy、vx、vy为参考子区的位移梯度。同样的,在变形更为复杂的情况下,为了提高测量精度,还可以采用更为复杂的二阶映射函数:x=x+u+uxx+uyy+12u2x2 x2+12u2y2 y2+uxuyx y,y=y+v+vxx+vyy+12v
9、2x2 x2+12v2y2 y2+vxvyx y,(3)本文综合考虑焊接变形计算精度与速度,选择一阶映射函数进行相关性计算。2 热输入与焊缝收缩力关系模型 图2 薄板材焊接方案示意图对于不同热输入的焊接变形检测试验,采用基于数字图像相关技术的三维检测装置进行薄板材变形测量,试验方案如图2所示。在试验准备阶段,需要测量变形的薄板材一侧制备高温黑白散斑,根据薄板材幅面对装置进行标定,将T I G焊枪固定在焊接机器人上,通过焊接机器人来控制焊接路径。T I G焊接机器人和测量装置分置于薄板材两侧,试验开始时焊接和测量同时作业,焊接完成薄板材冷却一定时间后,测量装置停止对薄板材变形图像采集。将所有采集
10、的变形图像导入装置后处理软件中进行位移场、应变场以及关键点的动态变形62安 徽 工 程 大 学 学 报第3 7卷计算。检测装置相机型号A c a 1 6 0 0-2 0 g m、分辨率16 2 6p i x 12 3 6p i x、像素大小4.4m4.4m,标定板视场4 0 0mm3 0 0mm,设置子区大小为3 0p i x3 0p i x,步长为1 5p i x1 5p i x,焊接冷却时间为6m i n。不同热输入焊接参数如表1所示。表1 热输入参数薄板材材料焊件尺寸mmmmmm热效率焊接电流A焊接电压U/(V)热输入Q n e t/(J/mm)焊接速度v/(mm/s)氩气流量q/(L/
11、m i n)Q 2 3 53 0 02 0 030.89 01 3.61 9 651 0Q 2 3 53 0 02 0 030.89 01 52 1 651 0Q 2 3 53 0 02 0 030.81 1 01 52 6 451 0Q 2 3 53 0 02 0 030.81 1 01 62 8 1.651 0Q 2 3 53 0 02 0 030.81 2 01 52 8 851 0Q 2 3 53 0 02 0 030.81 3 01 53 1 251 0Q 2 3 53 0 02 0 030.81 3 01 6.53 4 3.251 0Q 2 3 53 0 02 0 030.81 5
12、 01 7.54 2 051 0采用变形检测装置获得数据计算热输入与焊缝收缩力之间的关系,在焊缝中心线上,受边界效应影响,薄板材5 02 0 0mm处冷却变形较均匀,而且焊缝收缩力与压缩塑性应变和塑性变形区有关,与焊缝长度无关,因此可以在此范围内计算焊接收缩变形量。在薄板材散斑参考图像上全场划分网格,焊前经过一系列网格点L1(x1,y1,z1)、L2(x2,y2,z2)、L3(x3,y3,z3),Ln(xn,yn,zn),则每一小段的长度为LiLi+1=(xi+1-xi)2+(yi+1-yi)2+(zi+1-zi)2,(4)焊前焊缝中心线上的面内长度为L焊前=n-1i=1(xi+1-xi)2+
13、(yi+1-yi)2+(zi+1-zi)2。(5)焊接冷却后,焊缝中心线上经过的一系列网格点坐标变成L1(x1,y1,z1)、L2(x2,y2,z2)、L3(x3,y3,z3),Ln(xn,yn,zn),焊后焊缝横截线中心的面内长度为L焊后=n-1i=1(xi+1-xi)2+(yi+1-yi)2+(zi+1-zi)2,(6)则焊缝中心线上焊接前后的纵向面内收缩计算公式为L=n-1i=1(xi+1-xi)2+(yi+1-yi)2+(zi+1-zi)2-n-1i=1(xi+1-xi)2+(yi+1-yi)2+(zi+1-zi)2,(7)式中,LiLi+1为第i点和第i+1点的距离(mm);(xi,
14、yi,zi)为Li点的坐标;(xi+1,yi+1,zi+1)为Li+1点的坐标;L焊前为截线焊接前的面内长度(mm);L焊后为截线焊接后的面内长度(mm)。同理,焊缝横截线上的面内收缩也可以用上述相同的方法进行计算得到:T=n-1j=1(xj+1-xj)2+(yj+1-yj)2+(zj+1-zj)2-n-1j=1(xj+1-xj)2+(yj+1-yj)2+(zj+1-zj)2。(8)在薄板材的焊缝收缩变形区,假设薄板材存在收缩力,在收缩力Ff的作用下,构件在焊缝处产生收缩,在不考虑边界效应的情况下,焊缝收缩力为=FfLE A,(9)式中,A为薄板材的截面积;L为薄板材长度;E为薄板材材料的弹性
15、模量。根据上述分析,计算不同热输入与焊缝收缩量之间的关系如图3所示。经过线性拟合后得到热输入与焊缝收缩量之间的数学模型关系如下:=0.0 0 14Qn e t-0.2 5 26,(1 0)72第6期郭文举,等:热输入与薄板材T I G焊接面内外变形关系研究 图3 热输入与焊缝收缩量的关系结合式(9),可得焊缝收缩力与热输入的关系如下:Ff=0.8 6 1(Qn e t-1 8 0),(1 1)式中,4 2 0 J/mm Qn e t1 9 6 J/mm。对于薄板材堆焊在实际工程中的应用,可以根据两者之间的相互关系,从而控制薄板材焊接电流、焊接电压以及焊接速度,使得薄板材经过焊接之后不会出现薄板
16、材的失稳变形,从而提高焊接结构件的使用寿命。3 热输入对薄板材面外变形的影响采用基于热弹塑性有限元法对不同热输入下的薄板材焊接全场面外变形进行研究,之前的研究1 0已全面验证了薄板材焊接失稳变形有限元预测模型的精度,基于高斯热源模型、非线性瞬态热传导边界条件、材料高温性能参数等的热-力耦合热弹塑性预测模型具有较高的精度。焊接参数分别为薄板材尺寸3 0 0mm2 0 0mm3mm,焊接电压1 7.5V,焊接速度5mm/s,焊接电流取7 0A、9 0A、1 1 0A、1 3 0A、1 5 0A,对应的热输入分别为1 9 6 J/mm、2 5 2 J/mm、3 0 8J/mm、3 6 4J/mm、4 2 0 J/mm。薄板材焊接完成后的全场变形云图如图4所示。从图4中可以看出,薄板材在热输入范围内均呈反马鞍状,随着热输入的增大,最大负向面外变形和最大正向面外变形也随之增大。其中,最大负向面外位移增大速度较慢,仅从-0.8 7 5mm增加到-2.6 4 9mm,而最大正向面外位移增大速度较快,从1.3 3 6mm增大到8.3 4 1mm,这是因为薄板材放置于支架上,薄板材两侧除受支撑力外其余无
