高级检索

留言板

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

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

表面张力对激光焊接熔池及匙孔的影响

彭进 许红巧 杨晓红 王星星 李帅 李俐群 张健

彭进, 许红巧, 杨晓红, 王星星, 李帅, 李俐群, 张健. 表面张力对激光焊接熔池及匙孔的影响[J]. 焊接学报. doi: 10.12073/j.hjxb.20220502001
引用本文: 彭进, 许红巧, 杨晓红, 王星星, 李帅, 李俐群, 张健. 表面张力对激光焊接熔池及匙孔的影响[J]. 焊接学报. doi: 10.12073/j.hjxb.20220502001
PENG Jin, XU Hongqiao, YANG Xiaohong, WANG Xingxing, LI Shuai, LI Liqun, ZHANG Jian. Effect of surface tension on molten pool and keyhole during laser welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION. doi: 10.12073/j.hjxb.20220502001
Citation: PENG Jin, XU Hongqiao, YANG Xiaohong, WANG Xingxing, LI Shuai, LI Liqun, ZHANG Jian. Effect of surface tension on molten pool and keyhole during laser welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION. doi: 10.12073/j.hjxb.20220502001

表面张力对激光焊接熔池及匙孔的影响

doi: 10.12073/j.hjxb.20220502001
基金项目: 国家自然科学基金面上项目(52071165);河南省优秀青年基金项目(202300410268);河南省博士后科研项目(202003076);河南省重点研发与推广专项(科技攻关)项目(222102220082);河南省高校科技创新人才(22HASTIT026);河南省高层次人才特殊支持“中原英才计划-中原青年拔尖人才”(ZYBJ202126)
详细信息
    作者简介:

    彭进,博士,讲师;主要从事激光、电弧焊接方面的科研和教学工作;Email: pengjin1985120@163.com

    通讯作者: 杨晓红,教授. Email: 20050626@jhc.edu.cn.
  • 中图分类号: TG 456.7

Effect of surface tension on molten pool and keyhole during laser welding

  • 摘要: 基于FLUENT19.0软件,建立了激光焊接热-流耦合模型,对比分析了不同表面张力温度系数(为负值)对熔池流场的影响. 结果表明,随着表面张力温度系数的减小,熔池后方顺时针漩涡的流动趋势逐渐减弱,甚至消失,而且焊接飞溅的数量增多. 纵截面熔池长度逐渐增加,纵截面熔池流体最大流动速度逐渐增大,熔池横截面的面积逐渐减小.当表面张力温度系数为−2.5 × 10−4 N/(m·K)时,熔池长度平均值为3.28 mm、熔池流体最大流动速度的平均值为2.89 m/s、熔池横截面面积的平均值为4.52 mm2;当表面张力温度系数为−3.5 × 10−4 N/(m·K)时,熔池长度平均值为3.73 mm、熔池流体最大流动速度的平均值为3.53 m/s、熔池横截面面积的平均值为4.03 mm2;当表面张力温度系数为−4.9 × 10−4 N/(m·K)时,熔池长度平均值为4.14 mm、熔池流体最大流动速度的平均值为4.09 m/s、熔池横截面面积的平均值为3.28 mm2.
  • 图  1  熔池纵截面和横截面截取示意图

    Figure  1.  Schematic diagram of longitudinal and cross sections of molten pool

    图  2  激光自熔焊纵截面熔池流场(−2.5 × 10−4 N/(m·K))

    Figure  2.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −2.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  3  激光自熔焊纵截面熔池流场(−3.5 × 10−4 N/(m·K))

    Figure  3.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −3.5 × 10−4N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  4  激光自熔焊纵截面熔池流场(−4.9 × 10−4 N/(m·K))

    Figure  4.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −4.9 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  5  不同时间步的熔池长度尺寸

    Figure  5.  Length of molten pool at different time steps

    图  6  不同时间步的熔池流体最大流动速度

    Figure  6.  Maximum flow velocity of melt in molten pool at different time steps

    图  7  激光自熔焊纵截面熔池流场(−2.5 × 10−4 N/(m·K))

    Figure  7.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −2.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  8  激光自熔焊纵截面熔池流场(−3.5 × 10−4 N/(m·K))

    Figure  8.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −3.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  9  激光自熔焊纵截面熔池流场(−4.9 × 10−4 N/(m·K))

    Figure  9.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −4.9 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

    图  10  不同时间步的熔池面积尺寸

    Figure  10.  Area of molten pool at different time steps

    图  11  试验熔池图像与数值模拟熔池图像

    Figure  11.  Weld pool images obtained by experiment and numerical simulation

    表  1  6056铝合金热物理性能参数

    Table  1.   Thermophysical property parameters of 6056 aluminum alloy

    固相密度
    ${\rho _{\rm{s}}}$/(kg·m−3
    液相密度
    ${\rho _{\rm{l}}}$/(kg·m−3
    固相线温度
    ${T_{\rm{s}}}$/K
    液相线温度
    ${T_{\rm{L}}}$/K
    气相线温度
    ${T_{\rm{v}}}$/K
    熔化潜热
    ${L_{\rm{m}}}$/(J·kg−1
    蒸发潜热
    ${L_{\rm{v}}}$/(J·kg−1
    热膨胀系数
    ${\beta _{\rm{k}}}$/$ {{\rm{K}}^{ - 1}} $
    传热系数
    ${h_{\rm{c}}}$/(${\rm{W} } \cdot { {\rm{K} }^{ - 1} } \cdot { {\rm{m} }^{ - 2} }$
    发射率
    $ \varepsilon $
    环境温度
    ${T_{{\rm{ref}}} }$/K
    2 7002 5908609172 7403.87 × 1051.08 × 1071.92 × 105150.08300
    下载: 导出CSV
  • [1] Cao Huijun, Li Shiqin, Zhang Yinggan, et al. Mg2Sn-induced whisker growth on the surfaces of Mg/Sn/Mg ultrasonic-assisted soldering joints[J]. China Welding, 2022, 31(1): 47 − 59.
    [2] Lu Xiaohong, Sun Xudong, Sun Shixuan, et al. Numerical simulation of thermal field of FSW 2219 aluminum alloy thick plate[J]. China Welding, 2021, 30(4): 1 − 8.
    [3] Jiang Chao, Long Weimin, Feng Jian, et al. Thermal fatigue behavior of copper/stainless steel explosive welding joint[J]. China Welding, 2021, 30(4): 25 − 29.
    [4] 王传洋, 姜沐晖, 龙庆, 等. 激光工艺参数对PC/Cu/PC焊接性能及残余应力影响[J]. 焊接学报, 2021, 42(1): 24 − 29. doi:  10.12073/j.hjxb.20201019002

    Wang Chuanyang, Jiang Muhui, Long Qing, et al. Influence of laser process parameters on PC/Cu/PC welding performance and residual stress[J]. Transactions of the China Welding Institution, 2021, 42(1): 24 − 29. doi:  10.12073/j.hjxb.20201019002
    [5] 王磊, 许雪宗, 王克鸿, 等. 中厚板7A52铝合金光纤激光焊接接头组织与性能[J]. 焊接学报, 2021, 41(10): 28 − 31.

    Wang Lei, Xu Xuezong, Wang Kehong, et al. Microstructures and mechanical properties of fiber laser beam welded 7A52 alloy joint[J]. Transactions of the China Welding Institution, 2021, 41(10): 28 − 31.
    [6] Ye Xiaohu, Chen Xi. Three-dimensional modelling of heat transfer and fluid flow in laser full-penetration welding[J]. Journal of Physics D:Applied Physics, 2002, 35(10): 1049 − 1056.
    [7] Ribic B, Tsukamoto S, Rai R, et al. Role of surface-active elements during keyhole-mode laser welding[J]. Journal of Physics D:Applied Physics, 2011, 44: 485203. doi:  10.1088/0022-3727/44/48/485203
    [8] 樊丁, 郝珍妮, 黄勇, 等. 氧含量对TIG焊瞬态熔池行为影响的数值分析[J]. 兰州理工大学学报, 2013, 39(3): 18 − 21. doi:  10.3969/j.issn.1673-5196.2013.03.005

    Fan Ding, Hao Zhenni, Huang Yong, et al. Numerical analysis of effect of oxygen content on transient welding pool behavior on TIG welding[J]. Journal of Lanzhou University of Technology, 2013, 39(3): 18 − 21. doi:  10.3969/j.issn.1673-5196.2013.03.005
    [9] 董文超, 陆善平, 李殿中, 等. 微量活性组元氧对焊接熔池Marangoni 对流和熔池形貌影响的数值模拟[J]. 金属学报, 2008, 44(2): 249 − 256. doi:  10.3321/j.issn:0412-1961.2008.02.023

    Dong Wenchao, Lu Shanping, Li Dianzhong, et al. Numerical simulation of effects of the minor active-element oxygen on the marangoni convection and the weld shape[J]. Acta Metallurgica Sinica, 2008, 44(2): 249 − 256. doi:  10.3321/j.issn:0412-1961.2008.02.023
    [10] 庞盛永. 激光深熔焊接瞬态小孔和运动熔池行为及相关机理研究[D]. 武汉: 华中科技大学, 2011.

    Pang Shengyong. A study on the transient keyhole and moving weld pool behaviors and mechanisms of deep penetration laser welding[D]. Wuhan: Huazhong University of Science and Technology, 2011.
    [11] 李俐群, 郝雨, 彭进. 表面张力对激光深熔焊熔池流动的影响[J]. 焊接学报, 2019, 40(2): 13 − 19.

    Li Liqun, Hao Yu, Peng Jin. Effect of surface tension on flow in laser deep penetration welding molten pool[J]. Transactions of the China Welding Institution, 2019, 40(2): 13 − 19.
    [12] Yang Zhibin, Zhao Xin, Tao Wang, et al. Effects of keyhole status on melt flow and flow-induced porosity formation during double-sided laser welding of AA6056/AA6156 aluminium alloy T-joint[J]. Optics and Laser Technology, 2019, 109: 39 − 48. doi:  10.1016/j.optlastec.2018.07.065
    [13] Yan Zhihong, Shi Bowei, Sun Luping, et al. Surface defect detection of aluminum alloy welds with 3D depth image and 2D gray image[J]. The International Journal of Advanced Manufacturing Technology, 2020, 110: 741 − 752. doi:  10.1007/s00170-020-05882-x
    [14] Voller V R, Prakash C. A fixed grid numerical modeling methodology for convection-diffusion mushy region phase-change problems[J]. International Journal of Heat and Mass Transfer, 1987, 30(8): 1709 − 1719. doi:  10.1016/0017-9310(87)90317-6
    [15] Peng Jin, Xu Hongqiao, Yang Xiaohong, et al. Numerical simulation of molten pool dynamics in laser deep penetration welding of aluminum alloys[J]. Crystals, 2022, 12(6): 873. doi:  10.3390/cryst12060873
    [16] 张高磊, 孔华, 邹江林, 等. 高功率光纤激光深熔焊接飞溅特性以及离焦量对飞溅的影响[J]. 中国激光, 2021, 48(22): 2202008. doi:  10.3788/CJL202148.2202008

    Zhang Gaolei, Kong Hua, Zou Jianglin, et al. Spatter char-acteristics of high-power fibre laser deep penetration weldi-ng and effect of defocus on spatter[J]. Chinese Journal of Lasers, 2021, 48(22): 2202008. doi:  10.3788/CJL202148.2202008
  • [1] 康悦, 赵艳秋, 李悦, 陈丹, 綦娜, 占小红.  驱动力对2219铝合金DLBSW焊接熔池行为的影响 . 焊接学报, doi: 10.12073/j.hjxb.20210901001
    [2] 顾玉芬, 边春红, 李春凯, 石玗.  不同组元活性剂对不锈钢脉冲TIG焊熔池表面张力的影响 . 焊接学报, doi: 10.12073/j.hjxb.20190911001
    [3] 韩晓辉, 马寅, 马国龙, 杨海锋, 徐良.  双光束激光焊匙孔动态特征分析 . 焊接学报, doi: 10.12073/j.hzxb.20190811002
    [4] 黄瑞生, 邹吉鹏, 宫建锋, 杨义成, 梁晓梅.  激光扫描焊接熔池及等离子体动态行为 . 焊接学报, doi: 10.12073/j.hjxb.20191016004
    [5] 李俐群, 郝雨, 彭进.  表面张力对激光深熔焊熔池流动的影响 . 焊接学报, doi: 10.12073/j.hjxb.2019400034
    [6] 王宏宇,丁瑞,黄爱国,阚鹏.  背反射增效激光焊接熔池匙孔相变及流场分析 . 焊接学报, doi: 10.12073/j.hjxb.2018390137
    [7] 刘政军,王晓慧,苏允海.  磁场和活性焊联合作用对焊接熔池形态的影响 . 焊接学报, doi: 10.12073/j.hjxb.20150410002
    [8] 黄勇, 李慧, 王新鑫, 姚宇航.  不同驱动力对熔池表面变形行为影响的数值模拟 . 焊接学报,
    [9] 王厚勤, 张秉刚, 王廷, 冯吉才.  电子束定点焊接304不锈钢熔池流动行为数值模拟 . 焊接学报,
    [10] 刘成财, 刘铖丹, 刘巍巍, 何景山.  2219铝合金电子束焊接匙孔演变过程的数值模拟 . 焊接学报,
    [11] 占小红, 米高阳, 陶汪, 陈洁, 魏艳红, 陈彦宾.  薄板铝合金激光深熔焊熔池流动数值模拟 . 焊接学报,
    [12] 胡仙, 刘南生, 魏昇, 魏义庆.  钨极氩弧焊熔池表面高度测量系统纵向标定 . 焊接学报,
    [13] 高向东, 张勇, 游德勇, KatayamaSeiji.  大功率光纤激光焊熔池形态及焊接稳定性分析 . 焊接学报,
    [14] 汪任凭, 雷永平, 史耀武.  激光深熔焊中匙孔形成过程的动态模拟 . 焊接学报,
    [15] 王志江, 张广军, 高洪明, 吴林.  焊接现场脉冲GTAW熔池图像动态检测 . 焊接学报,
    [16] 王宏, 史耀武, 巩水利.  激光深熔焊熔池流动的换热特性 . 焊接学报,
    [17] 芦凤桂, 姚舜, 楼松年, 张毓新.  熔池表面变形对电弧行为特征的影响 . 焊接学报,
    [18] 周琦, 刘方军.  电子束深穿透焊接过程匙孔动力学研究进展 . 焊接学报,
    [19] 雷永平, 史耀武, 村川英一.  表面张力与合金汽化对激光加热表面温度的影响 . 焊接学报,
    [20] 冯雷, 陈树君, 殷树言.  高速焊接时焊缝咬边的形成机理 . 焊接学报,
  • 加载中
图(11) / 表 (1)
计量
  • 文章访问数:  147
  • HTML全文浏览量:  29
  • PDF下载量:  26
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-02
  • 网络出版日期:  2022-09-13

表面张力对激光焊接熔池及匙孔的影响

doi: 10.12073/j.hjxb.20220502001
    基金项目:  国家自然科学基金面上项目(52071165);河南省优秀青年基金项目(202300410268);河南省博士后科研项目(202003076);河南省重点研发与推广专项(科技攻关)项目(222102220082);河南省高校科技创新人才(22HASTIT026);河南省高层次人才特殊支持“中原英才计划-中原青年拔尖人才”(ZYBJ202126)
    作者简介:

    彭进,博士,讲师;主要从事激光、电弧焊接方面的科研和教学工作;Email: pengjin1985120@163.com

    通讯作者: 杨晓红,教授. Email: 20050626@jhc.edu.cn.
  • 中图分类号: TG 456.7

摘要: 基于FLUENT19.0软件,建立了激光焊接热-流耦合模型,对比分析了不同表面张力温度系数(为负值)对熔池流场的影响. 结果表明,随着表面张力温度系数的减小,熔池后方顺时针漩涡的流动趋势逐渐减弱,甚至消失,而且焊接飞溅的数量增多. 纵截面熔池长度逐渐增加,纵截面熔池流体最大流动速度逐渐增大,熔池横截面的面积逐渐减小.当表面张力温度系数为−2.5 × 10−4 N/(m·K)时,熔池长度平均值为3.28 mm、熔池流体最大流动速度的平均值为2.89 m/s、熔池横截面面积的平均值为4.52 mm2;当表面张力温度系数为−3.5 × 10−4 N/(m·K)时,熔池长度平均值为3.73 mm、熔池流体最大流动速度的平均值为3.53 m/s、熔池横截面面积的平均值为4.03 mm2;当表面张力温度系数为−4.9 × 10−4 N/(m·K)时,熔池长度平均值为4.14 mm、熔池流体最大流动速度的平均值为4.09 m/s、熔池横截面面积的平均值为3.28 mm2.

English Abstract

彭进, 许红巧, 杨晓红, 王星星, 李帅, 李俐群, 张健. 表面张力对激光焊接熔池及匙孔的影响[J]. 焊接学报. doi: 10.12073/j.hjxb.20220502001
引用本文: 彭进, 许红巧, 杨晓红, 王星星, 李帅, 李俐群, 张健. 表面张力对激光焊接熔池及匙孔的影响[J]. 焊接学报. doi: 10.12073/j.hjxb.20220502001
PENG Jin, XU Hongqiao, YANG Xiaohong, WANG Xingxing, LI Shuai, LI Liqun, ZHANG Jian. Effect of surface tension on molten pool and keyhole during laser welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION. doi: 10.12073/j.hjxb.20220502001
Citation: PENG Jin, XU Hongqiao, YANG Xiaohong, WANG Xingxing, LI Shuai, LI Liqun, ZHANG Jian. Effect of surface tension on molten pool and keyhole during laser welding[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION. doi: 10.12073/j.hjxb.20220502001
    • 焊接技术在国内工业领域具有较好应用和发展前景[1-3]. 激光焊接在激光加工制造中发挥着巨大的作用. 与钨极氩弧焊、熔化极气体保护焊等相比,激光焊接具有能量密度高、热影响区小、焊接速度高等优点[4-5].

      在激光焊接过程中,表面张力对熔池流动行为有一定程度的影响. 熔池中不同位置的表面张力的数值不同,这是因为表面张力与温度有直接关系,温度越高表面张力的值越低[6]. 在焊接试验中,添加氧等活性元素可增大熔深[7]. 一些学者基于FLUENT软件研究了表面张力对TIG焊接过程熔池的影响[8-9],发现氧含量对Marangoni流力的影响较大,进而影响焊缝的熔深和熔宽. 庞盛永[10]对表面张力分别为0.2, 0.6和1.0 N/m时的匙孔、熔池行为进行了研究,发现表面张力较小时匙孔的最大深度和最小深度均最大,匙孔表面受到的表面张力越小,匙孔越不容易闭合. 李俐群等人[11]研究了表面张力作用下的熔池流动行为,发现表面张力温度系数为负值时熔池表面以匙孔为中心由内向外流动,表面张力温度系数为正值时,熔池表面为由外向内流动.

      铝合金具有较高的热导率、线膨胀系数以及较低的熔点. 激光焊接铝合金时易产生气孔、咬边等焊接缺陷[12-13].

      基于以上分析,文中通过FLUENT19.0软件建立激光焊接热-流耦合模型,研究不同表面张力温度系数(为负值)对激光焊接熔池行为的影响规律.

    • 通过FLUENT19.0软件对激光自熔焊熔池行为进行数值模拟. 激光焊接模型的控制方程为[14]

      连续性控制方程

      $$ \frac{\partial \rho}{\partial {t}} + \frac{\partial(\rho u)}{\partial {x}} + \frac{\partial(\rho v)}{\partial {y}} + \frac{\partial(\rho w)}{\partial {\textit{z}}}=0 $$ (1)

      动量守恒方程,x方向为

      $$ \frac{{\partial \left( {\rho {{u}}} \right)}}{{\partial t}} + \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 {\textit{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 {\textit{z}}}}\left( {\mu \frac{{\partial u}}{{\partial {\textit{z}}}}} \right) + {S_{\text{u}}} $$ (2)

      y方向为

      $$ \frac{{\partial \left( {\rho {{v}}} \right)}}{{\partial t}} + \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 {\textit{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 {\textit{z}}}}\left( {\mu \frac{{\partial v}}{{\partial {\textit{z}}}}} \right) + {S_{\rm{v}}} $$ (3)

      z方向为

      $$ \frac{{\partial \left( {\rho {{w}}} \right)}}{{\partial t}} + \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 {\textit{z}}}} = $$
      $$ - \frac{{\partial P}}{{\partial {\textit{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 {\textit{z}}}}\left( {\mu \frac{{\partial w}}{{\partial {\textit{z}}}}} \right) + {S_{\text{w}}} $$ (4)

      能量守恒方程为

      $$ \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 {\textit{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 {\textit{z}}}}\left( {k\frac{{\partial T}}{{\partial {\textit{z}}}}} \right) + {S_{\rm{E}}} $$ (5)

      式中:u, v, w分别为沿x, y, z方向的速度;P为压力;$ {u_0} $为焊接速度;Su, Sv, Sw分别为动量方程沿x, y, z方向的源项;H0为混合焓;SE为能量方程的源项;$ \rho $为母材的密度;$ T $为温度;t为焊接时间;$ k $为母材导热系数;$ \mu $为液态流体的动力粘度.

      模型计算的工艺参数为:激光功率2 800 W,焊接速度为3.0 m/min. 激光束热源模型公式为

      $$ {q_{{\rm{laser}}}} = \frac{{9{\alpha _{{\rm{abs}}}}Q}}{{{\text{π}} R_0^2{H_0}\left( {1 - {e^3}} \right)}}\exp \left[ {\frac{{ - 9\left( {{x^2} + {y^2}} \right)}}{{R_0^2\log \left( {{H_0}/{\textit{z}}} \right)}}} \right] $$ (6)

      式中:Q为激光束的能量;αabs为相应的系数;H0为激光热源高度;R0为激光有效半径. 数学模型中添加的反冲压力、表面张力、浮力见文献[15].

      为提高计算的精确度,软件中设置的时间步长为1 × 10−5. 焊接工件的材料为6056铝合金,其热物理性能参数如表1所示.

      表 1  6056铝合金热物理性能参数

      Table 1.  Thermophysical property parameters of 6056 aluminum alloy

      固相密度
      ${\rho _{\rm{s}}}$/(kg·m−3
      液相密度
      ${\rho _{\rm{l}}}$/(kg·m−3
      固相线温度
      ${T_{\rm{s}}}$/K
      液相线温度
      ${T_{\rm{L}}}$/K
      气相线温度
      ${T_{\rm{v}}}$/K
      熔化潜热
      ${L_{\rm{m}}}$/(J·kg−1
      蒸发潜热
      ${L_{\rm{v}}}$/(J·kg−1
      热膨胀系数
      ${\beta _{\rm{k}}}$/$ {{\rm{K}}^{ - 1}} $
      传热系数
      ${h_{\rm{c}}}$/(${\rm{W} } \cdot { {\rm{K} }^{ - 1} } \cdot { {\rm{m} }^{ - 2} }$)
      发射率
      $ \varepsilon $
      环境温度
      ${T_{{\rm{ref}}} }$/K
      2 7002 5908609172 7403.87 × 1051.08 × 1071.92 × 105150.08300
    • 为分析表面张力温度系数为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)对熔池流场的影响,分别截取熔池纵截面和熔池横截面,沿匙孔中心截取获得熔池横截面. 熔池纵截面和熔池横截面示意图如图1所示.

      图  1  熔池纵截面和横截面截取示意图

      Figure 1.  Schematic diagram of longitudinal and cross sections of molten pool

    • 图2 ~ 图4分别为表面张力温度系数为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)的不同时间步熔池纵截面流场. 可以发现,随着表面张力温度系数的减小,熔池后方顺时针流动漩涡的趋势逐渐减弱,甚至消失,这主要是因为随着表面张力温度系数的减小,Marangoni作用力增强造成的. 而且随着表面张力温度系数的减小,焊接飞溅的数量增多. 分析认为,激光入射到焊接工件,激光作用区域的材料熔化、蒸发形成匙孔,在蒸发反冲压力、表面张力、浮力等力的作用下,熔池表面产生凸起. 凸起的液态金属脱离熔池形成飞溅需要满足下式[16],即

      图  2  激光自熔焊纵截面熔池流场(−2.5 × 10−4 N/(m·K))

      Figure 2.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −2.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      图  3  激光自熔焊纵截面熔池流场(−3.5 × 10−4 N/(m·K))

      Figure 3.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −3.5 × 10−4N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      图  4  激光自熔焊纵截面熔池流场(−4.9 × 10−4 N/(m·K))

      Figure 4.  Flow field of longitudinal section molten pool in autogenous laser welding: surface tension temperature coefficient was −4.9 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      $$ \rho_{{\rm{m}}} v_{{\rm{m}}}^{2}>\frac{2 \delta}{r} $$ (7)

      式中:δ为表面张力;r为熔池表面凸起金属的半径;ρm为熔池金属密度;vm为熔池表面凸起金属的流动速度.

      当焊接过程处于稳定状态时,匙孔内壁的反冲压力pr与熔池表面凸起金属流动速度为

      $$ p_{{\rm{r}}}=\frac{\rho_{{\rm{m}}} v_{{\rm{m}}}^{2}}{2} $$ (8)

      由公式(7)和公式(8)可得到

      $$ p_{{\rm{r}}}>\frac{\delta}{r} $$ (9)

      随着表面张力系数的减小,熔池相同位置的表面张力减小,因此与表面张力温度系数为−2.5 × 10−4, −3.5 × 10−4 N/(m·K)相比,表面张力温度系数为−4.9 × 10−4 N/(m·K)的熔池更容易产生焊接飞溅(图4).

      图2 ~ 图4的熔池表面均出现了由匙孔开口向熔池边缘流动的趋势,这主要是Marangoni作用力引起的. 当表面张力温度系数为−2.5 × 10−4 N/(m·K)时,匙孔底部出现了闭合,产生了焊接气泡(图2c),而当表面张力温度系数为−3.5 × 10−4, −4.9 × 10−4 N/(m·K)时,匙孔底部未产生焊接气泡.

      对表面张力温度系数分别为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)的不同时间步熔池长度进行统计,如图5所示(图中虚线为熔池长度的平均值). 由图可知,随着表面张力温度系数的减小,熔池长度逐渐增加. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池长度平均值为3.28 mm;当表面张力温度系数为−3.5 × 10−4 N/(m·K),熔池长度平均值为3.73 mm;当表面张力温度系数为−4.9 × 10−4 N/(m·K),熔池长度平均值为4.14 mm.

      图  5  不同时间步的熔池长度尺寸

      Figure 5.  Length of molten pool at different time steps

      对表面张力温度系数分别为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)的不同时间步熔池流体最大流动速度进行统计,如图6所示(图中虚线为熔池流体最大流动速度的平均值). 由图可知,随着表面张力温度系数的减小,熔池流体最大流动速度逐渐增大. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池流体最大流动速度的平均值为2.89 m/s;当表面张力温度系数为−3.5 × 10−4N/(m·K),熔池流体最大流动速度的平均值为3.53 m/s;当表面张力温度系数为−4.9 × 10−4N/(m·K),熔池流体最大流动速度的平均值为4.09 m/s.

      图  6  不同时间步的熔池流体最大流动速度

      Figure 6.  Maximum flow velocity of melt in molten pool at different time steps

    • 图7 ~ 图9分别为表面张力温度系数为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)的不同时间步熔池横截面流场(左侧为熔池流场方向,右侧为熔池流线). 可以发现,随着表面张力温度系数的减小,由熔池边缘指向匙孔的流动趋势变强,这主要是由于表面张力温度系数的变小,Marangoni作用力增强造成的.

      图  7  激光自熔焊纵截面熔池流场(−2.5 × 10−4 N/(m·K))

      Figure 7.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −2.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      图  8  激光自熔焊纵截面熔池流场(−3.5 × 10−4 N/(m·K))

      Figure 8.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −3.5 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      图  9  激光自熔焊纵截面熔池流场(−4.9 × 10−4 N/(m·K))

      Figure 9.  The longitudinal section flow field of the molten pool in autogenous laser welding: surface tension temperature coefficient was −4.9 × 10−4 N/(m·K). (a) step 2100; (b) step 2300; (c) step 2600; (d) step 2900

      对表面张力温度系数分别为−2.5 × 10−4, −3.5 × 10−4, −4.9 × 10−4 N/(m·K)的不同时间步熔池横截面面积进行统计,如图10所示(图中虚线为熔池面积的平均值). 由图可知,随着表面张力温度系数的减小,熔池横截面的面积逐渐减小. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池面积的平均值为4.52 mm2;当表面张力温度系数为−3.5 × 10−4 N/(m·K),熔池面积的平均值为4.03 mm2;当表面张力温度系数为−4.9 × 10−4 N/(m·K),熔池面积的平均值为3.28 mm2.

      图  10  不同时间步的熔池面积尺寸

      Figure 10.  Area of molten pool at different time steps

    • 为验证数值模拟模型的正确性,采用型号为$ {\text{camrecord 5000}} \times {\text{2}} $的高速摄像机对激光焊接过程进行实时监测,并采用辅助光源对熔池区域进行照亮. 试验参数与数值模拟参数相同(表面张力温度系数为−3.5 × 10−4 N/(m·K)),图11为高速摄像机拍摄的熔池图像与数值模拟熔池图像的对比验证. 通过对比发现,高速摄像机拍摄的熔池图像与数值模拟熔池图像的轮廓基本一致.

      图  11  试验熔池图像与数值模拟熔池图像

      Figure 11.  Weld pool images obtained by experiment and numerical simulation

    • (1) 随着表面张力温度系数的减小,熔池后方顺时针流动漩涡的趋势逐渐减弱,甚至消失. 而且随着表面张力温度系数的减小,焊接飞溅的数量增多.

      (2) 随着表面张力温度系数的变小,纵截面熔池长度逐渐增加. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池长度平均值为3.28 mm;当表面张力温度系数为−3.5 × 10−4 N/(m·K),熔池长度平均值为3.73 mm;当表面张力温度系数为−4.9 × 10−4 N/(m·K),熔池长度平均值为4.14 mm.

      (3) 随着表面张力温度系数的减小,纵截面熔池流体最大流动速度逐渐增大. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池流体最大流动速度的平均值为2.89 m/s;当表面张力温度系数为−3.5 × 10−4 N/(m·K),熔池流体最大流动速度的平均值为3.53 m/s;当表面张力温度系数为−4.9 × 10−4 N/(m·K),熔池流体最大流动速度的平均值为4.09 m/s.

      (4) 随着表面张力温度系数的减小,熔池横截面的面积逐渐减小. 当表面张力温度系数为−2.5 × 10−4 N/(m·K),熔池面积的平均值为4.52 mm2;当表面张力温度系数为−3.5 × 10−4N/(m·K),熔池面积的平均值为4.03 mm2;当表面张力温度系数为−4.9 × 10−4 N/(m·K),熔池面积的平均值为3.28 mm2.

参考文献 (16)

目录

    /

    返回文章
    返回