Effect of joint property differences on residual stresses in electron beam welded joints of Zr-Sn-Nb alloys
-
摘要:
为了提升Zr-Sn-Nb合金电子束焊残余应力计算精度,对Zr-Sn-Nb合金及焊接接头热物性参数与力学性能进行测量,研究接头性能差异对电子束焊接头残余应力影响规律,建立了电子束焊复合热源模型,采用热弹塑性有限元方法对4.45mm厚Zr-Sn-Nb合金电子束焊应力进行数值模拟,并进行焊缝形貌验证. 结果表明,考虑焊缝与母材属性差异时,试板上表面焊缝及热影响区区域的横向残余应力峰值高于二者属性相同时的模拟结果,峰值分别为319 MPa和296 MPa,进一步考虑热影响区属性与母材及焊缝差异对接头残余应力影响不大,峰值纵向应力为318 MPa.
-
关键词:
- Zr-Sn-Nb合金 /
- 接头性能差异 /
- 电子束焊 /
- 残余应力
Abstract:In order to improve residual stress numerical simulation accuracy of Zr-Sn-Nb alloy electron beam, thermophysical parameters and mechanical properties of Zr-Sn-Nb alloys and welded joints were measured. and the influence of the difference in joint properties on the residual stress of electron beam welded joints was investigated. The hybrid heat source model of electron beam welding is established, and the thermoelastic-plastic finite element method is used to numerically simulate the stress of electron beam welded 4.45 mm thick Zr-Sn-Nb alloy. The simulation model was verify by weld morphology. The results show that the peak transverse residual stresses in the weld and heat-affected zone region on the upper surface of the test plate are higher than those simulated when the properties of the two are the same, with peaks of 319 MPa and 296 MPa, respectively. The peak longitudinal stress was 318 MPa when the differences in the properties of heat-affected zone and the base material and weld are considered, which has little effect on the residual stresses in the joint.
-
0. 序言
增材制造技术是一种互联网 + 传统工艺相结合发展起来的,是20世纪80年代提出的,是一种基于离散-堆积原理的新型制造技术[1],随着机器人和计算机的不断发展,更容易实现工业智能制造,已经成为工业4.0的核心技术之一[2],相比于减材制造和等材制造具有生产效率高、加工成本低、材料利用率高、用途广泛等优点[3]. 在航空航天、在医疗器械、汽车制造和消费电子等行业应用较多,用于高精度、复杂模具成批制造[4],而受热源本身的局限性,难以实现大型复杂构件的一次精密成形,因此,增材制造技术在大型3D部件整体制造面临不少挑战[5],因此在重大装备制造业领域应用较少.
电弧熔丝增材制造作为最有可能实现大型复杂金属制件批量工业制造的一种增材制造技术[6],具有材料利用率高、无需密闭真空腔和设备制造成本低等特点. 相比于 MIG 电弧,TIG 电弧以其焊接过程稳定、飞溅小、空气污染小的特点深受研究学者们的青睐,但是效率低制约了其在大尺寸构件制造中的应用[7-9].
为了提高单电弧制造效率,围绕TIG 电弧衍生出来多种新型TIG复合热源制造方法,包括双钨极TIG[10]、三TIG电弧[11]、空心钨极TIG[12]、GPCA-TIG焊[13]、A-TIG[14]、激光-双TIG电弧复合[15]等焊接方法,研究表明利用复合电弧见得相互作用,降低了电弧压力,提高了电弧的稳定性,解决了在大熔覆电流条件下因电弧压力过大导致咬边、凹坑等问题[16],但现有技术仍难以实现大型复杂构件大电流、高焊速、高质量增材制造. 多热源增材制造过程是一个复杂多物理场耦合过程,涉及到多热源耦合能量吸收与转化、热量/质量传输、热源与溶体交互作用等复杂的冶金过程,其中多电弧耦合物理特性及高效稳定燃烧机制是影响热源耦合能量吸收与转化的关键因素,而关于这方面的研究未见报道. 特别是,对于大型复杂金属制件大多采用单电弧增材制造方法,金属熔敷效率低、生产制造周期长,极大限制了电弧增材制造技术在大型构件的应用推广,而多电弧复合焊接方法为其提供了思路.
为此,本文针对大尺寸复杂构件高效优质制造提出了一种四钨极TIG电弧增材制造方法,借助高速摄像采集了四钨极TIG电弧形态图像,通过提取从电弧引燃放电、稳定燃烧过程中电弧等离子体面积及波动等特征参数,掌握关键工艺参数对焊道成形的影响规律,分析四钨极TIG电弧耦合物理特性及稳定燃烧机制,为获得适合于大尺寸复杂构件增材制造的熔宽大、熔深浅、尺寸规则的焊道成形提供理论依据.
1. 试验方法
试验以30mm厚的304不锈钢作为基板,尺寸为300mm × 200mm × 30mm,金属丝材选用哈尔滨焊接研究所自研的HS13/5焊丝,焊丝直径为1.2mm,基板和焊丝的化学成分如表1所示. 设计了一套四钨极TIG电弧双丝增材制造系统,主要由焊接电源、送丝装置、四钨极焊接电源控制系统、三维行走试验台及自主研发的四钨极焊枪组成.
表 1 304不锈钢与HS13/5焊丝成分(质量分数,%)Table 1. Chemical composition of 304 stainless steel and HS13/5 welding wire材料 C Si Mn S P Cr Ni Mo 304不锈钢 0.040 0.43 1.17 0.0014 0.0280 18.05 8.08 0.054 HS13/5焊丝 0.016 0.46 0.54 0.0083 0.0019 12.30 4.51 0.480 焊接电源采用四台EWM Tetrix 352直流焊接电源;送丝机为两台EWM TS drive送丝机,采用前送丝的送丝方式,两根金属丝材采用同向两侧位置排布,金属丝材间夹角为20°,金属丝材与基板表面夹角为15°,其中增材制造设备示意图如图1所示. 四钨极热源呈空间对称分布,钨极沿z轴方向的倾角为15°,四钨极热源空间排布图如图2所示. 每个钨极的参数相同,相互间位置参数一致.
通过观察四钨极TIG电弧引燃阶段、稳定燃烧阶段电弧等离子体变化过程,采用单一变量法,定量表征不同钨极间距、沉积电流、电弧弧长下四钨极TIG电弧的变化规律,提取四钨极TIG电弧等离子体面积特征参数如图3所示,其中S:电弧等离子体面积(单位:mm2、pixel);h:电弧高度(单位:mm、pixel);l:电弧弧长(单位:mm、pixel).
2. 试验结果与分析
2.1 四TIG电弧引燃放电过程
图4为四钨极间距为2mm时多TIG电弧引燃放电过程及稳定燃烧过程电弧等离子体变化图像,可以看出,每个电弧燃烧过程都包括引燃放电、稳定燃烧两个阶段,且稳定燃烧过程随着脉冲电流峰值、基值交替形成峰值稳定燃烧、基值稳定燃烧两种电弧形态. 第一个电弧起弧以基板为触发点以高频起弧方式引燃,并沿着钨极中心延长线方向稳定燃烧,后面的第二、三、四个电弧起弧并没有以基板为触发点起弧,均是沿着已燃烧电弧的边缘起弧,起弧均伴随着电弧等离子体的扩大和电弧公共导电区域的形成两个阶段,另外,随着不同电弧的陆续燃烧,多电弧弧柱区偏离了钨极中心延长线均有向着四钨极TIG电弧中心相互吸引的趋势,且随着电弧的增多,偏离角度也逐渐增大,具体表现为多电弧向中心汇聚
为了清晰描述多电弧耦合作用机制,选取了第一电弧及第二电弧起弧过程典型图像并通过Image J图像处理软件对其进行二值化处理,如图5所示. 图a1为第一电弧稳定时的电弧形态;b1为第二电弧起弧瞬间电弧等离子体的扩大;c1为第一电弧与第二电弧公共导电通道形成前;d1为第一电弧与第二电弧公共导电通道形成,两电弧处于稳定状态,a2、b2、c2、d2则表示二值化处理后的图片. 此外,对4个阶段电弧等离子体面积(黑色部分)进行数据统计,其结果如图6所示,当电弧状态由a至b时,第一电弧稳定燃烧、第二电弧开始引燃的瞬间,首先发生的是电弧等离子体面积的扩大;随后由b至c时,第二电弧弧柱区形成,电弧等离子体面积瞬间扩张;最后由c至d,两个电弧间形成了稳定的公共导电区域,提高了弧柱区的能量密度,四钨极TIG电弧的自磁收缩效应增强,电弧等离子体面积减小[17],因此,由b至d阶段,电弧等离子体面积先增加后减小.
图7所示为四个电弧的引燃时间,由图可知,随引燃时序,四钨极TIG电弧的引燃时间逐渐缩短. 这是因为电弧1在引燃前没有导电通道提供“便利”,只能在基板处起弧,而当后续电弧起弧时,按照最小电压原理,电弧总是在最有利于电弧引燃的位置起弧,前一个已经存在的电弧为其提供了公共导电通道,使其跳过基板在已存在电弧处直接引燃,且能量密度的增加也为电弧引燃提供了“便利”,缩短了后续电弧的引燃时间.
2.2 钨极间距对四TIG电弧形态的影响
图8为不同钨极间距下四钨极TIG电弧形态的正面拍摄图像,由图可知,当钨极间距为2mm、4mm时,四钨极TIG电弧呈现出类曲边三角形的电弧等离子体形态,而不是传统单钨极电弧等离子的钟罩形,当钨极间距增加到为2mm、4mm时;四钨极TIG电弧在钨极上方仍存在电弧等离子体,并且可以明显观察到随钨极间距的增加,电弧沿钨极中心线方向的偏转角度逐渐增加. 图9为不同钨极间距下四钨极TIG电弧面积及面积方差,由图可知,随钨极间距增加,四钨极TIG电弧面积及波动均增加,当钨极间距≤6mm时,电弧面积及波动变化较小,当钨极间距增加至8mm、10mm时,电弧面积及波动骤增,说明电弧稳定性下降;结合电弧形态图像可知,虽然电弧面积整体增加,但电弧间公用导电通道的亮度显著降低,说明此处能量密度减小.
2.3 沉积电流对四TIG电弧形态的影响
图10为不同沉积电流下四钨极TIG电弧形态,由图可知,在同一钨极间距下,沉积电流对四钨极TIG电弧形态影响不大,同一钨极间距下沉积电流增加时,四钨极TIG电弧高亮区面积略有增加,电弧间的相互吸引增强,公共导电区域增加,多电弧更容易耦合成一个大电弧,且多电弧沿四钨极几何中心向上汇聚的等离子高亮区亦增加,但与钨极间距相比,沉积电流对多电弧耦合作用效果相对较小.
图11为不同沉积电流下四钨极TIG电弧面积及面积方差,由图可知,随沉积电流增加,四钨极TIG电弧面积及波动略有增加,与钨极间距作用相比较小,且当沉积电流峰值达到180A时,电弧面积方差显著增加,说明电弧稳定性下降. 分析认为,当沉积电流增加时,四钨极TIG电弧间耦合作用力增加,特别是四钨极TIG电弧沿几何中心向上汇聚的等离子聚集区域扩大,电子碰撞几率增加引起电弧稳定性下降,这可通过文中2.6节四钨极TIG电弧稳定机制中电子迁移现象解释.
2.4 电弧弧长对四TIG电弧形态的影响
图12为不同电弧弧长的四钨极TIG电弧图像,可以观察到随电弧弧长增加,四钨极TIG电弧等离子体面积明显增加,当电弧亮度下降,代表电弧能量密度降低,当钨极间距为2mm时,随电弧弧长变化,四钨极TIG电弧形态由类“陀螺”形转变为类曲边三角形,而当钨极间距为10mm时,随电弧弧长增加,四钨极TIG电弧形态变化不明显. 图13为不同电弧弧长下四钨极TIG电弧面积及面积方差,由图可知,随电弧弧长增加,四钨极TIG电弧面积及波动显著增加,特别是,当电弧弧长达到7mm时,电弧面积方差显著增加,说明电弧稳定性下降.
2.5 四钨极TIG电弧热源有效利用率分析
为了更好的表征四钨极TIG电弧有效利用率,引入熔化能的概念,并以此作为四钨极TIG电弧有效利用率的表征方式,其中熔化能公式[17]如式1所示.
熔化能公式:
$$ {\mathrm{E}}_{\mathrm{m}} = \left(\mathrm{\nu }\mathrm{S}\right)\mathrm{\rho }\left({\int }_{{\mathrm{T}}_{0}}^{\mathrm{T}\mathrm{m}}\mathrm{C}d\mathrm{T} + \mathrm{Y}\right) $$ (1) 式中 Em——熔化能 J
V——扫描速度 mm/min
S——横截面熔化面积 mm2
ρ——材料密度 g/mm3
T——温度 ℃
T0——室温 ℃
Tm——熔点 ℃
C——比热容 J/(g·℃)
Y——熔化潜热 J/g
由该公式可知,熔化能的大小与横截面熔化面积密切相关,即当横截面面积较大时,即可认为此时四钨极TIG电弧有效利用率较高. 因此实验条件相同条件下以自熔的方式对不同钨极间距下四钨极TIG电弧有效利用率进行了焊接试验,以此确定最佳的耦合距离,试验及测量结果见表2和图14所示.
表 2 不同钨极间距下四钨极TIG电弧形态及接头横截面Table 2. Cross-sectional area of four tungsten electrode heat source fusion joint for different tungsten electrode spacing)钨极间距(d/mm) 2 4 6 8 10 电弧形态 接头横截面 由表2可知,四钨极间距为2mm 时,焊缝熔宽与接头横截面面积均最大,而间距为10mm时两者均最小,且随着钨极间距的增大,熔宽没有明显增加,且熔化面积逐渐减小,这是因为焊道尺寸不仅与作用在基体的热源能量面积有关,还与能量密度有关,即热源能量有效利用率,当钨极间距较小时,四钨极TIG电弧能量面积虽然较小,但能量密度较大,此时能量有效利用率较高,而当钨极间距较大时,虽然四钨极TIG电弧能量面积较大,但能量密度较小,而能量有效利用率较低. 因此,钨极间距为2mm时,焊道熔化面积最大;钨极间距为10mm时,熔化面积最小,这是由不同钨极间距时四钨极TIG电弧热源有效利用率的不同造成的. 利用熔化能公式计算得出:钨极间距为2mm的熔化能约为钨极间距为10mm时熔化能的9.2倍.
2.6 四钨极TIG电弧稳定机制
从电子迁移角度解释了上述关键参数对四钨极TIG电弧形态的影响结果,从而揭示四钨极TIG电弧稳定燃烧机制. 图15为四钨极TIG电弧电子迁移情况示意图,如图所示,单钨极时,电子运动方向基本是沿着钨极中心线延长线运动的,没有发生偏转,而当四个钨极共同放电时,每个钨极的电子运动方向发生偏转并向四钨极几何中心汇聚,钨极间距增大,偏转现象更加明显. 分析认为,钨极间距为2-6mm时,四钨极TIG电弧在安培力的作用下形成公共导电通道[18],显著提高四钨极TIG电弧能量密度,增强自磁收缩效应,大部分电子沿钨极延长线方向迁移,最终到达基板位置,只有少部分电子在电弧之间安培力的作用下向上运动,最终形成高度较小的弧峰,此时四钨极TIG电弧能量利用率和稳定性均下降;当钨极间距为8-10mm时,由于距离较远不能形成公共导电通道,但仍然在安培力的作用下发生偏转,因此多电弧向钨极几何中心处汇聚,相互撞击,并最终向钨极上方迁移形成弧峰,此时四钨极TIG电弧能量利用率和稳定性均下降,此外,沉积电流与电弧弧长都会影响电子密度与电子运动方向,因此,两者对电弧形态的影响也可利用该理论进行解释.
3. 结论
(1)针通过对多电弧耦合关键参数对四钨极TIG电弧形态及稳定性的研究表明:钨极间距、电弧弧长可以显著的影响四钨极TIG电弧形态,沉积电流对四钨极TIG电弧形态影响相对较小;其中,钨极间距对四钨极TIG电弧稳定性的影响最为明显,当钨极间距为2mm时,四钨极TIG电弧稳定性最好.
(2)对不同钨极间距下四钨极TIG电弧热源有效利用率进行研究,结果表明,当钨极间距为2mm时,四钨极TIG电弧热源有效利用效率最高,相较于钨极间距为10mm时提高近9.2倍.
(3) 四钨极TIG电弧在自磁收缩和安培力的作用下相互吸引,当钨极间距较小时(2-6mm),四钨极TIG电弧弧柱区相互重叠形成公共导电通道,此时四钨极TIG电弧稳定性最好,热源利用效最高;当钨极间距较大时(8-10mm),四钨极TIG电弧不能形成稳定的导电通道,大部分电子在安培力的作用下相互碰撞,最终向钨极上方运动,造成四钨极TIG电弧稳定性下降、热源有效利用效下降,并且对焊枪的烧损降低了其使用寿命.
-
表 1 Zr-Sn-Nb合金电子束焊接工艺参数
Table 1 Zr-Sn-Nb alloy electron beam welding parameters
加速电压
U/kV聚焦电流
If /mA焊接束流
Ib /mA焊接速度
vs /(mm·s−1)工作距离
d/mm60 706* + 8 27 480 190.5 *706 mA为试板上表面聚焦电流. 表 2 Zr-Sn-Nb合金焊接残余应力数值模拟与测试对比
Table 2 Comparison between numerical simulation and test of welding residual stress of Zr-Sn-Nb alloy
验证位置 应力方向 应力测量值δ1/MPa 应力模拟值δ2/MPa 第1组 第2组 第3组 焊缝 σx −109.3 −60.8 −112.8 −82.6 σy 114.8 103 151.6 143.3 热影响区 σx 67.5 46 52.2 87.3 σy 265.8 268 273.8 259.7 母材 σx 46.7 31 22.3 34.8 σy −47.3 −48 −50.4 −47.5 -
[1] 贾豫婕, 林希衡, 邹小伟, 等. 锆合金的研发历史、现状及发展趋势[J]. 中国材料进展, 2022, 41(5): 354 − 370. Jia Yujie, Lin Xiheng, Zou Xiaowei, et al. Research & development history, status and prospect of zirconium alloys[J]. Materials China, 2022, 41(5): 354 − 370.
[2] 王旭峰, 李中奎, 周军, 等. 锆合金在核工业中的应用及研究进展[J]. 热加工工艺, 2012, 41(2): 71 − 74. Wang Xufeng, Li Zhongkui, Zhou Jun, et al. Application and research progress of zirconium alloy in nuclear industry[J]. Hot Working Technology, 2012, 41(2): 71 − 74.
[3] Zhao Wanqian, Song Pengcheng, Peng Xiaoming, et al. Fracture mode analysis of Zr-Sn-Nb alloy under simulated LOCA condition in advanced nuclear reactors[J]. Engineering Failure Analysis, 2022, 137: 106392. doi: 10.1016/j.engfailanal.2022.106392
[4] 邢硕, 蒲曾坪, 张坤, 等. 新型锆合金包壳蠕变性能评价方法研究[J]. 核动力工程, 2023, 44(4): 234 − 239. Xing Shuo, Pu Zengping, Zhang Kun, et al. Study on evaluation method for creep performance of new zirconium alloy cladding[J]. Nuclear Power Engineering, 2023, 44(4): 234 − 239.
[5] 钟建伟, 安军靖, 丁怀博, 等. Zr-Sn-Nb-Fe-Cr与Zr-Nb-Fe锆合金电阻点焊工艺及显微组织[J]. 焊接学报, 2021, 42(8): 82 − 90. Zhong Jianwei, An Junjing, Ding Huaibo, et al. Welding processes and microstructures of weld bead of Zr-Sn-Nb-Fe-Cr and Zr-Nb-Fe zirconium alloy[J]. Transactions of the China Welding Institution, 2021, 42(8): 82 − 90.
[6] Slobodyan M S. Arc welding of zirconium and its alloys: a review[J]. Progress in Nuclear Energy, 2021, 133: 103630. doi: 10.1016/j.pnucene.2021.103630
[7] Slobodyan M S. Resistance, electron and laser-beam welding of zirconium alloys for nuclear applications: A review[J]. Nuclear Engineering and Technology, 2021, 53(4): 1049 − 1078. doi: 10.1016/j.net.2020.10.005
[8] 王博, 包陈, 魏连峰, 等. 氢化物对锆合金薄板焊缝断裂行为的影响[J]. 机械工程学报, 2021, 57(20): 133 − 140. doi: 10.3901/JME.2021.20.133 Wang Bo, Bao Chen, Wei Lianfeng, et al. Effect of hydride on fracture behavior of zirconium alloy platy welds[J]. Journal of Mechanical Engineering, 2021, 57(20): 133 − 140. doi: 10.3901/JME.2021.20.133
[9] 杨建国, 谢浩, 闫德俊, 等. 随焊干冰激冷冷源尺寸对焊接残余应力影响的有限元分析[J]. 焊接学报, 2017, 38(2): 14 − 18. Yang Jianguo, Xie Hao, Yan Dejun, et al. FEM analysis about effect of cooling source size during welding with dry ice on welding residual stress[J]. Transactions of the China Welding Institution, 2017, 38(2): 14 − 18.
[10] 杨帆, 陈芙蓉. A-UIT处理对7075铝合金焊接应力影响的数值模拟[J]. 焊接学报, 2021, 42(12): 91 − 96. Yang Fan, Chen Furong. Numerical simulation of effect of A-UIT treatment on welding stress of 7075 aluminum alloy[J]. Transactions of the China Welding Institution, 2021, 42(12): 91 − 96.
[11] 许谦, 胡广旭, 董志波, 等. 基于相变与收缩耦合的1Cr12焊缝冷却动态力学行为[J]. 中国机械工程, 2021, 32(3): 341 − 347. Xu Qian, Hu Guangxu, Dong Zhibo, et al. Dynamic mechanics behaviors of 1Cr12 weld cooling based on coupling of phase transformation and contraction[J]. China Mechanical Engineering, 2021, 32(3): 341 − 347.
[12] 张开元. 固态相变对超高强钢焊接-热处理过程应力变形的影响[D]. 合肥: 中国科学技术大学, 2022. Zhang Kaiyuan. Effect of solid-state phase transformation on stress deformation of ultra-high strength steel during weld-heat treatment [D]. Hefei: University of Science and Technology of China, 2022.
[13] Banik S D, Kumar S, Singh P K, et al. Distortion and residual stresses in thick plate weld joint of austenitic stainless steel: Experiments and analysis[J]. Journal of Materials Processing Technology, 2021, 289: 116944. doi: 10.1016/j.jmatprotec.2020.116944
[14] 吴婧, 陈静, 蒋亦岚, 等. 冷变形N18合金再结晶过程中的织构演变[J]. 金属热处理, 2021, 46(5): 38 − 46. Wu Jing, Chen Jing, Jiang Yilan, et al. Texture evolution of cold-rolled N18 alloy during recrystallization[J]. Heat Treatment of Metals, 2021, 46(5): 38 − 46.
[15] 张建军, 李中奎, 周军, 等. 新锆合金高温变形行为[C]//中国有色金属学会合金加工学术委员会. 中国有色金属学会第十二届材料科学与合金加工学术年会论文集. 北京: 中国有色金属学会, 2007: 374 − 378. Zhang Jianjun, Li Zhongkui, JiangYilan, et al. High temperature deformation behavior of new zirconium alloy [C]// Alloy Processing Academic Committee of China Nonferrous Metals Society. Proceedings of the 12th Annual Conference of Materials Science and Alloy Processing of the Chinese Non-Ferrous Metals Society. Beijing: Nonferrous Metals in China, 2007: 347 − 378.