留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

熔滴填充位置对激光焊接熔池动态行为的影响

彭进 王星星 许红巧 杨嘉佳 张芙蓉

彭进, 王星星, 许红巧, 杨嘉佳, 张芙蓉. 熔滴填充位置对激光焊接熔池动态行为的影响[J]. 红外与激光工程, 2020, 49(6): 20200025. doi: 10.3788/IRLA20200025
引用本文: 彭进, 王星星, 许红巧, 杨嘉佳, 张芙蓉. 熔滴填充位置对激光焊接熔池动态行为的影响[J]. 红外与激光工程, 2020, 49(6): 20200025. doi: 10.3788/IRLA20200025
Peng Jin, Wang Xingxing, Xu Hongqiao, Yang Jiajia, Zhang Furong. Effect of droplet filling position on dynamic behavior of molten pool in laser welding[J]. Infrared and Laser Engineering, 2020, 49(6): 20200025. doi: 10.3788/IRLA20200025
Citation: Peng Jin, Wang Xingxing, Xu Hongqiao, Yang Jiajia, Zhang Furong. Effect of droplet filling position on dynamic behavior of molten pool in laser welding[J]. Infrared and Laser Engineering, 2020, 49(6): 20200025. doi: 10.3788/IRLA20200025

熔滴填充位置对激光焊接熔池动态行为的影响

doi: 10.3788/IRLA20200025
基金项目: 河南省自然科学基金(162300410191);河南省高等学校重点科研项目(17A430021);华北水利水电大学高层次人才科研启动基金(201705006);华北水利水电大学大学生创新创业训练计划(2019XB005)
详细信息
    作者简介:

    彭进(1984-),男,讲师,博士,主要从事激光焊接方面的研究。Email:pengjin1985120@163.com

  • 中图分类号: TG456.7

Effect of droplet filling position on dynamic behavior of molten pool in laser welding

  • 摘要: 借助熔滴作用下的三维瞬态激光焊接热-流耦合有限元模型,对不同的熔滴填充位置下熔滴进入熔池过程的匙孔三维形貌、熔池金属流动特性进行研究。数值模拟计算结果表明,熔滴填充位置对激光焊接过程中匙孔三维形貌及熔池液态金属的流动行为的影响较大。当熔滴填充位置由0.5 mm增大到1.8 mm时,对匙孔三维形貌变化的影响减弱,熔池内部挤压匙孔前壁和后壁驱使匙孔闭合的流动趋势减弱而维持匙孔壁张开的流动趋势增强,匙孔底部液态金属流动速度的波动幅度减弱。
  • 图  1  模型计算流程图

    Figure  1.  Flow chart of model calculation

    图  2  熔滴填充位置图

    Figure  2.  Location of liquid droplet filling

    图  3  熔滴填充位置为0.5 mm时的匙孔形态

    Figure  3.  Keyhole morphology of liquid droplet filling position at 0.5 mm

    图  4  熔滴填充位置为1.8 mm时的匙孔形态

    Figure  4.  Keyhole morphology of liquid droplet filling position at 1.8 mm

    图  5  熔滴填充位置为0.5 mm和1.8 mm的匙孔深度变化

    Figure  5.  Keyhole depth change for liquid droplet filling positions at 0.5 mm and 1.8 mm

    图  6  熔滴填充位置为0.5 mm的熔池流场

    Figure  6.  Molten pool flow field for liquid droplet at 0.5 mm filling position

    图  7  熔滴填充位置为1.8 mm的熔池流场

    Figure  7.  Molten pool flow field for liquid droplet at 1.8 mm filling position

    图  8  熔滴填充位置为0.5 mm和1.8 mm的熔池内部H点流动速度变化

    Figure  8.  Flow velocity change at the H point inside the molten pool with liquid droplet filling position of 0.5 mm and 1.8 mm

    图  9  熔滴对匙孔稳定性影响的原理图

    Figure  9.  Schematic diagram for the keyhole stability influence by liquid droplet

  • [1] 杨志斌, 李浩哲, 王斌, 等. 铝合金T型接头双侧激光焊接气孔缺陷影响规律研究[J]. 红外与激光工程, 2018, 47(9): 0906004. doi:  10.3788/IRLA201847.0906004

    Yang Zhibin, Li Haozhe, Wang Bin, et al. Laser process for high strength aluminum alloys [J]. Infrared and Laser Engineering, 2018, 47(9): 0906004. (in Chinese) doi:  10.3788/IRLA201847.0906004
    [2] 石俊凯, 王国名, 纪荣祎, 等. 结构紧凑的双波长连续波掺铒光纤激光器[J]. 中国光学, 2019, 12(4): 811−819.

    Shi Junkai, Wang Guoming, Ji Rongwei, et al. Compact dual-wavelength continuous-wave Er-doped fiber laser [J]. Chinese Opics, 2019, 12(4): 811−819. (in Chinese)
    [3] 孟佳, 张伟, 赵开祺 等. 国产化掺铥光纤激光振荡器性能研究[J]. 中国光学, 2019, 12(5): 1109−1116. doi:  10.3788/co.20191205.1109

    Meng Jia, Zhang Wei, Zhao Kaiqi, et al. Investigation on the performance of a homemade thulium-doped fiber laser oscillator [J]. Chinese Opics, 2019, 12(5): 1109−1116. (in Chinese) doi:  10.3788/co.20191205.1109
    [4] Wu Y, Cai Y, Wang H, et al. Investigation on microstructure and properties of dissimilar joint between SA553 and SUS304 made by laser welding with filler wire [J]. Materials & Design, 2015, 87: 567−578.
    [5] 王韵致, 谢芳, 陈龙辉, 等. 用于高精度测量位移等参量的光纤多波长激光器[J]. 光学 精密工程, 2019, 27(9): 2036−2042. doi:  10.3788/OPE.20192709.2036

    Wang Yunzhi, Xie Fang, Chen Longhui, et al. Research on optical fiber multi-wavele-ngth laser for measuring displaceme-nt precisely [J]. Optics and Precision Engineering, 2019, 27(9): 2036−2042. (in Chinese) doi:  10.3788/OPE.20192709.2036
    [6] 李彦洲, 石岩. 铝合金表面激光沉积AlCrFeCoNiCu高熵合金组织与工艺优化[J]. 光学 精密工程, 2019, 27(4): 795−806. doi:  10.3788/OPE.20192704.0795

    Li Yanzhou, Shi Yan. Microstructure and process optimization of AlCrFeCoNiCu high-entropy alloy by laser depositi-on [J]. Optics and Precision Engineering, 2019, 27(4): 795−806. (in Chinese) doi:  10.3788/OPE.20192704.0795
    [7] 余阳春. 激光填丝焊的焊丝熔入行为及工艺研究[D]. 武汉: 华中科技大学, 2010: 60-83.

    Yu Yangchun. Study on the technology and filler wire melting dynamics during the laser welding with filler wire[D]. Wuhan: Huazhong University of Science and Technology, 2010: 60-83. (in Chinese)
    [8] 彭进, 李俐群, 林尚扬 等. 激光液态填充焊的熔滴熔化与过渡稳定性[J]. 焊接学报, 2016, 37(7): 9−12.

    Peng Jin, Li Liqun, Lin Shangyang, et al. Stability of filler melting and transfer with pre-melting liquid filler [J]. Transactions of the China Welding Institution, 2016, 37(7): 9−12. (in Chinese)
    [9] Sohail M, Han S W, Na S J, et al. Characteristics of weld pool behavior in laser welding with various power inputs [J]. Welding in the World, 2014, 58(3): 269−277. doi:  10.1007/s40194-014-0112-4
    [10] Wang H, Shi Y W, Gong S L. Numerical simulation of laser keyhole welding processes based on control volume methods [J]. Journal of Physics D: Applied Physics, 2006, 39(21): 4722−4730. doi:  10.1088/0022-3727/39/21/032
    [11] 彭进, 张文洁, 王星星, 等. 焊丝熔化填充方式对激光焊接熔池影响的数值模拟[J]. 红外与激光工程, 2018, 47(3): 0306005. doi:  10.3788/IRLA201847.0306005

    Peng Jin, Zhang Wenjie, Wang Xingxing, et al. Numerical simulation of effect of filler wire melting and filling mode on molten pool in laser welding [J]. Infrared and Laser Engineering, 2018, 47(3): 0306005. (in Chinese) doi:  10.3788/IRLA201847.0306005
    [12] 杨志斌. 铝合金机身壁板结构双侧激光焊接特征及熔池行为研究[D]. 哈尔滨: 哈尔滨工业大学, 2013: 74-90.

    Yang Zhibin. Study on double-sided laser beam welding characteristics and weld pool behavior for aluminium alloy fuselage panel structure[D]. Harbin: Harbin Institute of Technology, 2003, 17(6): 431-437. (in Chinese)
    [13] Matsunawa A, Kim J D, Seto N, et al. Dynamics of keyhole and molten pool in laser welding [J]. Journal of Laser Applications, 1998, 10(6): 247−254. doi:  10.2351/1.521858
  • [1] 张月来, 龙维峰, 姚杞, 林吉, 刘泽辉, 刘博文.  铝合金车体焊接同步激光清洗技术研究 (特邀) . 红外与激光工程, 2023, 52(2): 20220534-1-20220534-8. doi: 10.3788/IRLA20220534
    [2] 杨正伟, 陈金树, 田干, 张炜, 陈家威, 肖俊玲.  2195-T8铝锂合金激光填丝双面焊熔池行为与接头力学性能研究 . 红外与激光工程, 2023, 52(2): 20220350-1-20220350-9. doi: 10.3788/IRLA20220350
    [3] 彭进, 许红巧, 王永彪, 王星星, 张永振, 龙伟民, 张定宇.  激光光斑尺寸对激光焊接熔池、匙孔行为的影响 . 红外与激光工程, 2023, 52(7): 20220130-1-20220130-11. doi: 10.3788/IRLA20220130
    [4] 李耀明, 陈淑琴, 张煌.  基于激光谐波调制的深孔内壁三维面型分析系统 . 红外与激光工程, 2022, 51(3): 20210862-1-20210862-6. doi: 10.3788/IRLA20210862
    [5] 巴一, 韩善果, 杨永强, 郑世达, 师文庆.  激光摆动路径对钢/铝点焊接头组织与性能的影响 . 红外与激光工程, 2022, 51(8): 20210799-1-20210799-10. doi: 10.3788/IRLA20210799
    [6] Xu Baoxiang, Xiong Zhi, Huang Jixun, Yu Haicheng.  Research on method of improving magnetic field adaptability of high-precision IFOG . 红外与激光工程, 2021, 50(4): 20200239-1-20200239-8. doi: 10.3788/IRLA20200239
    [7] 钦兰云, 徐丽丽, 杨光, 刘奇, 王维.  激光沉积成形熔池尺寸分析与预测 . 红外与激光工程, 2018, 47(11): 1106009-1106009(7). doi: 10.3788/IRLA201847.1106009
    [8] 王聪, 石世宏, 方琴琴, 石拓, 夏志新.  封闭空腔回转薄壁件的激光熔覆成形 . 红外与激光工程, 2018, 47(1): 106006-0106006(8).
    [9] 李超, 申力, 施斐斐.  同轴空心激光检测圆柱内孔圆度误差的方法研究 . 红外与激光工程, 2018, 47(S1): 1-4. doi: 10.3788/IRLA201847.S106001
    [10] 杨志斌, 李浩哲, 王斌, 陈禹昕.  铝合金T型接头双侧激光焊接气孔缺陷影响规律研究 . 红外与激光工程, 2018, 47(9): 906004-0906004(8). doi: 10.3788/IRLA201847.0906004
    [11] 彭进, 张文洁, 王星星, 郭国全, 张芙蓉.  焊丝熔化填充方式对激光焊接熔池影响的数值模拟 . 红外与激光工程, 2018, 47(3): 306005-0306005(7). doi: 10.3788/IRLA201847.0306005
    [12] 杜道中, 刘婷婷, 廖文和, 张长东, 张凯.  激光选区熔化熔池光强监测系统设计 . 红外与激光工程, 2017, 46(12): 1206002-1206002(6). doi: 10.3788/IRLA201746.1206002
    [13] 邓志强, 石世宏, 周斌, 贾帅, 王聪.  弧形悬垂实体结构激光熔覆成形 . 红外与激光工程, 2017, 46(10): 1006004-1006004(6). doi: 10.3788/IRLA201762.1006004
    [14] 李东康, 刘玉娟, 田杏霞, 王涛, 付长宝.  多色相干光学烧孔的动态生成与调控 . 红外与激光工程, 2017, 46(11): 1106007-1106007(4). doi: 10.3788/IRLA201746.1106007
    [15] 杨光, 赵恩迪, 钦兰云, 李长富, 王维.  电磁搅拌对激光熔池熔体流速及其凝固组织影响研究 . 红外与激光工程, 2017, 46(9): 906006-0906006(8). doi: 10.3788/IRLA201746.0906006
    [16] 黄磊, 陈希章, 陈兴, 马宏波, 林涛.  镀锌高强钢激光焊接气孔形成机理 . 红外与激光工程, 2017, 46(S1): 1-6. doi: 10.3788/IRLA201746.S106001
    [17] 李东康, 刘玉娟.  激光场作用下原子系统中光学烧孔和慢光效应 . 红外与激光工程, 2016, 45(12): 1206005-1206005(4). doi: 10.3788/IRLA201645.1206005
    [18] Sekou Singare, 陈盛贵, 钟欢欢.  激光透射焊接聚碳酸酯的有限元分析 . 红外与激光工程, 2016, 45(2): 206005-0206005(6). doi: 10.3788/IRLA201645.0206005
    [19] 赵娟宁, 董晓娜, 达争尚, 李红光.  单脉冲高能激光系统动态焦斑位置测量 . 红外与激光工程, 2014, 43(6): 1925-1929.
    [20] 闫世兴, 董世运, 徐滨士, 王玉江, 任维彬, 方金祥.  激光熔覆过程中熔池对流运动对熔覆层气孔和元素分布的影响 . 红外与激光工程, 2014, 43(9): 2832-2839.
  • 加载中
图(9)
计量
  • 文章访问数:  349
  • HTML全文浏览量:  125
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-06
  • 修回日期:  2020-03-31
  • 刊出日期:  2020-07-01

熔滴填充位置对激光焊接熔池动态行为的影响

doi: 10.3788/IRLA20200025
    作者简介:

    彭进(1984-),男,讲师,博士,主要从事激光焊接方面的研究。Email:pengjin1985120@163.com

基金项目:  河南省自然科学基金(162300410191);河南省高等学校重点科研项目(17A430021);华北水利水电大学高层次人才科研启动基金(201705006);华北水利水电大学大学生创新创业训练计划(2019XB005)
  • 中图分类号: TG456.7

摘要: 借助熔滴作用下的三维瞬态激光焊接热-流耦合有限元模型,对不同的熔滴填充位置下熔滴进入熔池过程的匙孔三维形貌、熔池金属流动特性进行研究。数值模拟计算结果表明,熔滴填充位置对激光焊接过程中匙孔三维形貌及熔池液态金属的流动行为的影响较大。当熔滴填充位置由0.5 mm增大到1.8 mm时,对匙孔三维形貌变化的影响减弱,熔池内部挤压匙孔前壁和后壁驱使匙孔闭合的流动趋势减弱而维持匙孔壁张开的流动趋势增强,匙孔底部液态金属流动速度的波动幅度减弱。

English Abstract

    • 与传统的熔焊方法相比,激光焊接具有很高的能量密度[1-4]。而且激光焊接具有焊接效率高,热输入量小等优点[5-6]。由于焊丝的加入可以降低激光焊接对接板对接间隙的苛刻要求。通过选用合适的焊丝可以改善焊缝组织,进而提高焊缝的力学性能。

      然而激光填丝焊由于需要依靠激光的能量熔化焊丝,焊丝通常距离匙孔的距离很近,容易影响匙孔的稳定性。为了提高激光填丝焊的稳定性,国内的学者Yu[7]研究了焊接工艺参数对焊缝形貌的影响,研究发现焊接过程中熔池、匙孔的稳定性对焊缝成形有直接的影响。哈尔滨工业大学的彭进等[8]提出了激光液态填充焊新方法。通过研究发现激光液态填充焊焊丝前端的熔化金属(熔滴)与匙孔边缘的距离大于激光填丝焊焊丝前端的熔化金属(熔滴)与匙孔边缘的距离,熔化后可以沿熔池边缘流入,进而避免了焊丝对匙孔的直接冲击,匙孔未出现闭合的情况。

      为获得激光焊接过程的熔池流动及匙孔形态,很多学者采用数值模拟的方法对激光自熔焊过程进行相应的研究。国外的学者Sohail[9]等人采用数值模拟的方法研究了激光自熔焊熔池动态行为,研究发现提高焊接速度后熔池的尺寸减小,匙孔闭合的趋势减弱。国内的学者Wang[10]等采用数值模拟的方法研究了激光深熔焊过程的熔池行为,当激光的能量使焊接工件实现单面焊双面成形时熔池的上方和下方都会产生液态金属的流动旋涡。

      熔滴的加入会使激光深熔焊过程的熔池、匙孔动态行为更加复杂,然而国内外关于填充金属对熔池动态行为影响的研究相对较少。

      基于以上分析,结合作者前期建立的熔滴送入过程的三维瞬态激光焊接热-流耦合模型,研究了熔滴不同送进位置下液态金属进入熔池对匙孔形貌、熔池流场的影响。

    • 激光焊接过程的物理量要满足以下三大控制方程组。

      连续性控制方程

      $$ \frac{{\partial \rho }}{{\partial t}} + \frac{{\partial \left( {\rho u} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho v} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho w} \right)}}{{\partial x}} = 0 $$ (1)

      Navier-Stokes控制方程

      x方向:

      $$\begin{split} & \frac{{\partial \left( {\rho u} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho uu} \right)}}{{\partial x}} - \frac{{\partial \left( {\rho u{u_0}} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho uv} \right)}}{{\partial y}} + \frac{{\partial \left( {\rho uw} \right)}}{{\partial z}} = - \frac{{\partial P}}{{\partial x}}+\\ & \frac{\partial }{{\partial x}}\left( {\mu \frac{{\partial u}}{{\partial x}}} \right) + \frac{\partial }{{\partial y}}\left( {\mu \frac{{\partial u}}{{\partial y}}} \right) + \frac{\partial }{{\partial z}}\left( {\mu \frac{{\partial u}}{{\partial z}}} \right) + {S_u} \end{split} $$ (2)

      y方向:

      $$ \begin{split} & \frac{{\partial \left( {\rho v} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho vu} \right)}}{{\partial x}} - \frac{{\partial \left( {\rho v{u_0}} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho vv} \right)}}{{\partial y}} + \frac{{\partial \left( {\rho vw} \right)}}{{\partial z}} = - \frac{{\partial P}}{{\partial y}}+\\ & \frac{\partial }{{\partial x}}\left( {\mu \frac{{\partial v}}{{\partial x}}} \right) + \frac{\partial }{{\partial y}}\left( {\mu \frac{{\partial v}}{{\partial y}}} \right) + \frac{\partial }{{\partial z}}\left( {\mu \frac{{\partial v}}{{\partial z}}} \right) + {S_v} \end{split} $$ (3)

      z方向:

      $$\begin{split} & \frac{{\partial \left( {\rho w} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho wu} \right)}}{{\partial x}} - \frac{{\partial \left( {\rho w{u_0}} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho wv} \right)}}{{\partial y}} + \frac{{\partial \left( {\rho ww} \right)}}{{\partial z}} = - \frac{{\partial P}}{{\partial z}}+\\ & \frac{\partial }{{\partial x}}\left( {\mu \frac{{\partial w}}{{\partial x}}} \right) + \frac{\partial }{{\partial y}}\left( {\mu \frac{{\partial w}}{{\partial y}}} \right) + \frac{\partial }{{\partial z}}\left( {\mu \frac{{\partial w}}{{\partial z}}} \right) + {S_w} \end{split} $$ (4)

      能量守恒控制方程

      $$\begin{split} & \frac{{\partial \left( {\rho H} \right)}}{{\partial t}} + \frac{{\partial \left( {\rho uH} \right)}}{{\partial x}} - \frac{{\partial \left( {\rho {u_0}H} \right)}}{{\partial x}} + \frac{{\partial \left( {\rho vH} \right)}}{{\partial y}} + \frac{{\partial \left( {\rho wH} \right)}}{{\partial z}}=\\ & \frac{\partial }{{\partial x}}\left( {k\frac{{\partial T}}{{\partial x}}} \right) + \frac{\partial }{{\partial y}}\left( {k\frac{{\partial T}}{{\partial y}}} \right) + \frac{\partial }{{\partial z}}\left( {k\frac{{\partial T}}{{\partial z}}} \right) + {S_E} \end{split} $$ (5)

      式中:SuSvSw分别为动量方程沿xyz方向的源项;SE是能量源项;字母$u$$v$$w$分别为文中所建立的数学模型的xyz 方向的速度矢量;$\rho $是焊接工件的密度(文中所用材质为铝合金);P为压力;$\mu $为焊接材质熔化为液态后的动力粘度;H是混合焓。

    • 在焊接初始阶段,焊接工件在激光热源的作用下不断熔化,当激光功率密度达到一定值时,焊接匙孔形成,由于匙孔壁面为气-液界面,也就是说在匙孔边界存在液相和气相的不断动态变化(即能量的不断变化),文中采用VOF方法(Volume of Fluid)来处理匙孔壁面的气-液转化问题。

      VOF控制方程为:

      $$ \frac{{\partial {{F}}}}{{\partial {{t}}}}{{ + u}}\frac{{\partial {{F}}}}{{\partial {{x}}}}{{ + v}}\frac{{\partial {{F}}}}{{\partial {{y}}}}{\rm{ + w}}\frac{{\partial {{F}}}}{{\partial {{z}}}}{\rm{ = }}0 $$ (6)

      式中:uvw分别是液态金属沿xyz方向的流动速度。

      激光焊接过程中,匙孔壁面的焊接驱动力对匙孔稳定性起至关重要的作用。

      匙孔表面受液态金属蒸发的反冲压力Pr的作用,对于所建的模型可以将反冲压力以源项的方式引入到Navier-Stokes方程,反冲压力的方程如公式(7)所示:

      $$ {{{P}} _r} = A{B_0}T_s^{ - \frac{1}{2}}\exp ( - \frac{U}{{{T_s}}}) $$ (7)

      式中:符号B0为焊接工件的蒸发常数;U为焊接工件单位质量的能量;Ts为焊接工件的局部温度;A为常数。

      表面张力${p_\sigma }$的数学表达式如下:

      $$ {p_\sigma } = k\sigma $$ (8)

      式中:$\sigma $为表面张力系数;$k$是自由表面曲率,其数学表达式为:

      $$ k = - \left[ {\nabla .\left( {\frac{{\mathop n\limits^ \to }}{{\left| {\mathop n\limits^ \to } \right|}}} \right)} \right] = \frac{1}{{\left| {\mathop n\limits^ \to } \right|}}\left[ {\left( {\frac{{\mathop n\limits^ \to }}{{\left| {\mathop n\limits^ \to } \right|}}.\nabla } \right)\left| {\mathop n\limits^ \to } \right| - \left( {\nabla \mathop n\limits^ \to } \right)} \right] $$ (9)

      液态熔池的热浮力为:

      $$ \delta (T) = {\delta _0} + {A_\delta }(T - {T_m}) $$ (10)

      式中:${\delta _0}$的物理含义为纯金属熔点的表面张力系数;Tm的物理含义为模型中焊接母材的熔点。

      模型中熔滴模型的表达式以及热源模型、边界条件见作者发表的参考文献[11]。

    • 采用流体动力学软件FLUENT15.0对模型进行数值模拟计算。模型采用C/C++计算机语言对模型软件进行二次开发,将热源模型、焊接驱动力等引入计算模型中。由于文中主要分析焊丝端部熔化成熔滴后,熔滴在匙孔前壁不同位置填充对匙孔及熔池的影响,为了减少计算机计算过程的工作负担,不考虑焊丝端部由固态熔化成液态的过程,假设焊丝端部已经熔化并且以熔滴的形式流入熔池,因此通过采用FLUENT软件在建立的模型采用功能键“patch”熔滴。

      文中进行模拟采用的实际激光类型为光纤激光器。

      图  1  模型计算流程图

      Figure 1.  Flow chart of model calculation

      计算机计算过程中不断进行迭代求解,最终判定计算收敛并终止迭代的判断依据为:

      $$ \frac{{{\max}\left( {\left| {\phi _i^{n + 1} - \phi _i^n} \right|} \right)}}{{\phi _i^{n + 1}}} \leqslant {R_\phi } $$ (11)

      式中:$\phi _i^{n + 1}$为模型计算迭代的实时迭代值;$\phi _i^n$n+1迭代步的上一步迭代值。

      所建立的模型中考虑了温度对材料热物性参数的影响。焊接母材采用的是6056铝合金,其比热容、热导率以及动力粘度参数随温度变化的数学表达式如下[12]

      T≥573 K时,比热容的数学表达式为:

      $$ {C_{\rm{p}}}\left( \rm{J \cdot\rm {k{g^{ - 1}}}} \right) = - 0.001 \times {T^2} + 1.160\;9 \times T + 267.71 $$ (12)

      当573 KT≤913 K时,比热容的数学表达式为:

      $$ {C_{\rm{p}}}\left( \rm{J \cdot k{g^{ - 1}}} \right) = 0.000\;9 \times {T^2} - 0.390\;1 \times T + 514.45 $$ (13)

      T>913 K时,比热容的数学表达式为:

      $$ {C_{\rm{p}}}\left( \rm{J \cdot k{g^{ - 1}}} \right) = - 0.000\;1 \times {T^2} + 0.583\;2 \times T + 435.14 $$ (14)

      T>860 K时,热导率的数学表达式为:

      $$ k\left(\rm {W \cdot {m^{ - 1}} \cdot {K^{ - 1}}} \right) = 0.000\;1 \times {T^2} - 0.069\;7 \times T + 95.334 $$ (15)

      当860 K<T≤917 K时,热导率的数学表达式为:

      $$ k\left(\rm {W \cdot {m^{ - 1}} \cdot {K^{ - 1}}} \right) =-0.004\;8 \times {T^2}+9.281\;2 \times T-4\;275.6 $$ (16)

      当917 K<T≤2 740 K时,热导率的数学表达式为:

      $$ k\left( \rm{W \cdot {m^{ - 1}} \cdot {K^{ - 1}}} \right) =- 0.000\;01 \times {T^2}+0.058\;2 \times T + 148.74 $$ (17)

      T>2 740 K时,热导率的数学表达式为:

      $$ k\left(\rm {W \cdot {m^{ - 1}} \cdot {K^{ - 1}}} \right) =-0.000\;01 \times {T^2}+0.085 \times T + 71.111 $$ (18)

      当897 KT≤937 K时,动力粘度的数学表达式为:

      $$ \mu \left(\rm {kg \cdot {m^{ - 1}} \cdot {s^{ - 1}}} \right) = 1 \times {10^{ - 7}} \times {T^2}-0.000\;2 \times T + 0.120\;2 $$ (19)

      当937 K<T≤2 650 K时,动力粘度的数学表达式为:

      $$ \mu \left( \rm{kg \cdot {m^{ - 1}} \cdot {s^{ - 1}}} \right) = 2 \times {10^{ - 11}} \times {T^2}-5 \times {10^{ - 7}} \times T + 0.003\;8 $$ (20)

      当2 650 K<T≤2 720 K时,动力粘度的数学表达式为:

      $$ \mu \left(\rm {kg \cdot {m^{ - 1}} \cdot {s^{ - 1}}} \right) = - 6 \times {10^{ - 8}} \times {T^2}+0.000\;3 \times T - 0.415\;1 $$ (21)

      模型计算的流程图如图1所示。

      图  2  熔滴填充位置图

      Figure 2.  Location of liquid droplet filling

      模型计算的工艺参数:激光功率4 500 W,熔滴流入熔池的初始速度为3 m/min,焊接速度为3.0 m/min,熔滴送进位置距离匙孔前壁L分别为0.5 mm和1.8 mm。焊丝材料为6056铝合金,焊丝直径为1.2 mm。熔滴填充位置图如图2所示。

    • 文中首先研究了熔滴填充位置分别为0.5 mm和1.8 mm的匙孔三维形貌波动特性,分别如图3图4所示。由图3可知,熔滴填充位置为0.5 mm时,由于液态金属距离匙孔很近,当熔滴不断填送进熔池时,匙孔三维形貌变化较大,而且容易出现匙孔底部闭合的现象;而熔滴填充位置为1.8 mm时,由于熔滴与匙孔的距离较远,当熔滴流入熔池时,对匙孔的冲击较小,匙孔的三维形貌也会出现一定程度的波动,但未出现匙孔闭合的现象。

      图  3  熔滴填充位置为0.5 mm时的匙孔形态

      Figure 3.  Keyhole morphology of liquid droplet filling position at 0.5 mm

      图  4  熔滴填充位置为1.8 mm时的匙孔形态

      Figure 4.  Keyhole morphology of liquid droplet filling position at 1.8 mm

      图  5  熔滴填充位置为0.5 mm和1.8 mm的匙孔深度变化

      Figure 5.  Keyhole depth change for liquid droplet filling positions at 0.5 mm and 1.8 mm

      进一步对不同焊接时刻的熔滴填充位置分别为0.5 mm和1.8 mm的匙孔深度的变化进行了对比分析,如图5所示。由图5可知,熔滴填充位置为0.5 mm时,匙孔深度变化的幅度比熔滴填充位置为1.8 mm时匙孔深度变化的幅度要大。熔滴填充位置为0.5 mm时,匙孔深度变化的最大值为2.9 mm,最小值为1.1 mm;熔滴填充位置为1.8 mm时,匙孔深度变化的最大值为2.9 mm,最小值为1.7 mm。

      图  6  熔滴填充位置为0.5 mm的熔池流场

      Figure 6.  Molten pool flow field for liquid droplet at 0.5 mm filling position

    • 激光焊接过程中熔池内部液态金属的流动行为会影响到匙孔形态的稳定性,匙孔闭合时会产生焊接气泡,进而会影响焊缝质量。文中进一步对比分析熔滴填充位置分别为0.5 mm和1.8 mm的熔池纵截面的流场,分别如图6图7所示。由图6可知,当熔滴距离匙孔较近时,熔池内部产生了挤压匙孔前壁和后壁的流动趋势(图6(a)),在熔池金属的挤压匙孔壁的作用下匙孔底部出现了闭合的现象(图6(b)(c)),而且液态金属对匙孔底部的挤压趋势十分明显;由图7可知,当熔滴距离匙孔较远时,液态金属填充进入熔池,虽然也会产生液态金属流动挤压匙孔壁的现象,但这种挤压趋势很弱,维持匙孔壁张开的流动趋势较为明显。由图6图7可知,在Marangoni流力的作用下,在熔池上部出现液态金属流向熔池后方的逆时针流动现象。

      图  7  熔滴填充位置为1.8 mm的熔池流场

      Figure 7.  Molten pool flow field for liquid droplet at 1.8 mm filling position

      参考文献[13]表明激光深熔焊过程中匙孔形态不断出现波动,当匙孔下部波动剧烈时,匙孔底部容易产生焊接气泡。匙孔壁面液态金属的流动速度波动情况可以间接反映匙孔的稳定性。激光焊接过程中,填充金属的存在一定程度上增大了小孔前壁的流体静压力作用,对匙孔后壁也有一定程度的影响,这是由于当匙孔前壁内部有凸起时会反射一部分激光能量到匙孔后壁进而引起匙孔形态的波动,而匙孔的波动较大时会出现匙孔中下部断开的情况,进而产生气泡,因此文中对匙孔后壁下部的H点(如图8所示)的不同焊接时刻的流动速度进行观察,发现液态金属不同填充位置下的H点的流动速度也会出现一定程度的波动。由图8可知,熔滴填充位置为0.5 mm时熔池内部H点的流动速度波动较大。熔滴填充位置为0.5 mm时熔池内部H点的流动速度最大值和最小值分别为2 989 mm/s,118 mm/s。熔滴填充位置为0.5 mm时,熔池内部H点的流动速度平均值为967 mm/s;熔滴填充位置为1.8 mm时,熔池内部H点的流动速度波动相对较小。熔滴填充位置为1.8 mm时,熔池内部H点的流动速度最大值和最小值分别为1 953 mm/s,243 mm/s。熔滴填充位置为1.8 mm时,熔池内部H点的流动速度平均值为848 mm/s。

      图  8  熔滴填充位置为0.5 mm和1.8 mm的熔池内部H点流动速度变化

      Figure 8.  Flow velocity change at the H point inside the molten pool with liquid droplet filling position of 0.5 mm and 1.8 mm

    • 对于激光深熔焊,匙孔壁面处于实时波动状态,维持匙孔张开的作用力为反冲压力(Pr)和匙孔内部金属元素蒸发的蒸气压力(бPg);而驱使匙孔闭合的作用力为表面张力(Pб)以及熔池内部的液态金属流体静压力(Ph)。当匙孔处于稳定状态时,有如下公式:

      $$ {P_{{r}}} + \delta {P_g} = {P_\sigma } + {P_h} $$ (22)

      激光深熔焊过程中,由于熔滴的加入,很容易造成对匙孔前壁的冲击作用,而当匙孔前壁波动较大时,会影响激光束照射到匙孔内部的激光能量密度分布的均匀性,进而会影响激光深熔焊匙孔内部金属元素的蒸发量,在匙孔内部激光束照射到的区域,金属气化的蒸气反冲作用力增大,而激光束无法照射到的区域,金属气化的蒸气反冲作用力减小,进一步会造成匙孔形貌出现剧烈的波动,熔滴填充进入熔池对匙孔影响的原理图,如图9所示。

      图  9  熔滴对匙孔稳定性影响的原理图

      Figure 9.  Schematic diagram for the keyhole stability influence by liquid droplet

    • (1) 熔滴填充位置为0.5 mm时,由于熔滴距离匙孔很近,当液态金属不断填送进熔池时,匙孔三维形貌变化较大,而且容易出现匙孔底部闭合的现象;而熔滴填充位置为1.8 mm时,由于熔滴与匙孔的距离较远,当液态金属流入熔池时,对匙孔的冲击较小,匙孔的三维形貌也会出现一定程度的波动,但未出现匙孔闭合的现象。

      (2) 当熔滴距离匙孔较近时(L=0.5 mm),熔池内部产生了挤压匙孔前壁和后壁的流动趋势,在熔池金属的挤压匙孔壁的作用下匙孔底部出现了闭合的现象,这样容易产生焊接气泡,进而造成焊缝工艺气孔率增高;当熔滴距离匙孔较远时(L=1.8 mm),液态金属填充进入熔池时产生的液态金属流动挤压匙孔壁的挤压趋势很弱,维持匙孔壁张开的流动趋势较为明显,降低了因匙孔闭合形成的焊缝工艺气孔数量。

      (3) 熔滴填充位置为0.5 mm时熔池内部H点的流动速度波动较大,其流动速度最大值和最小值分别为2 989 mm/s,118 mm/s,熔池内部H点的流动速度平均值为967 mm/s;熔滴填充位置为1.8 mm时,熔池内部H点的流动速度波动相对较小,其流动速度最大值和最小值分别为1 953 mm/s,243 mm/s,熔池内部H点的流动速度平均值为848 mm/s。

参考文献 (13)

目录

    /

    返回文章
    返回