高级检索

留言板

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

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

选区激光熔化316L扫描特征参数对表面形貌影响

朱春霞 邱渤杰

朱春霞, 邱渤杰. 选区激光熔化316L扫描特征参数对表面形貌影响[J]. 焊接学报, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
引用本文: 朱春霞, 邱渤杰. 选区激光熔化316L扫描特征参数对表面形貌影响[J]. 焊接学报, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
ZHU Chunxia, QIU Bojie. Effect of scanning characteristic parameters on surface morphology of selective laser melting 316L[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
Citation: ZHU Chunxia, QIU Bojie. Effect of scanning characteristic parameters on surface morphology of selective laser melting 316L[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002

选区激光熔化316L扫描特征参数对表面形貌影响

doi: 10.12073/j.hjxb.20220426002
详细信息
    作者简介:

    朱春霞,教授,博士研究生导师;主要从事3D打印及其相关技术研究;已发表论文40余篇;Email: 11905269@qq.com

  • 中图分类号: TG 401; TH 162; TG 141; TG665

Effect of scanning characteristic parameters on surface morphology of selective laser melting 316L

  • 摘要: 为研究不同扫描特征参数组合对选区激光熔化(selective laser melting, SLM)表面形貌的影响,以316L不锈钢粉末为例,进行介观尺度的单层双道数值模拟研究.基于离散元法建立单层的粉床数值模型.使用流体体积法对粉床受热部分粉末的熔化过程中的熔化、流动和凝固过程进行计算.考虑激光功率、扫描速度和扫描间距3个扫描特征参数,设计并进行正交试验,从熔道形貌特征和熔道宽度2个方面研究所选扫描特征参数对成形件表面的熔道形貌影响.依据数值模拟中的参数进行实际打印及形貌观察试验,验证数值模拟的有效性.结果表明,在313 ~ 500 J/m的线能量密度和50 ~ 90 μm的扫描间距范围内,可以得到平整连续局部缺陷少的熔道形貌,且该区间内的参数组合依次线性对应;对熔道形貌的完整性影响由大到小依次为扫描速度、扫描间距和激光功率.
  • 图  1  316L热物参数

    Figure  1.  Theromphysicl parameters of 316L. (a) liquid viscosity; (b) specific heat; (c) thermal conductivity; (d) density

    图  2  数值模拟熔道形貌

    Figure  2.  Numerical simulation of melt tracks morphology. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2

    图  3  理想熔道形貌示意

    Figure  3.  Ideal melt track morphology

    图  4  线能量密度(LED)与扫描间距对熔道形貌特征影响

    Figure  4.  Influence of linear energy density (LED) and scanning spacing on the morphology characteristics of the melt tracks

    图  5  不同扫描间距下的线能量密度与熔宽

    Figure  5.  Line energy density and melting width under different scanning distances

    图  6  316L粉末形状轮廓图像

    Figure  6.  316L powder shape outline image

    图  7  316L粉末粒径分布

    Figure  7.  316L powder particle size distributio

    图  8  不同参数组合下的成形件表面形貌

    Figure  8.  Forming surface topography under the different parameters combination conditions. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2

    图  9  熔道宽度及误差

    Figure  9.  Width and error of the melt track

    表  1  316L材料参数

    Table  1.   316L material parameters

    固相密度ρ/(kg·m−3) 固相温度T/K 液相温度T/K 沸点温度T/K 熔化潜热Hf /(J·kg−1) 蒸发潜热Hv /(J·kg−1) 表面张力r /(N·m−1) 吸收率A
    7 850 1 658 1 723 3 090 2.7 × 105 7.45 × 106 1.6 0.25
    下载: 导出CSV

    表  2  正交试验设计

    Table  2.   Experimental design using orthogonal table

    参数变量A功率P/WB速度v/(m·s−1)C间距S/μm
    水平12000.650
    水平22500.970
    水平33001.290
    下载: 导出CSV

    表  3  试验参数及编组

    Table  3.   Experimental parameters and grouping

    组号功率 P/W速度 v /(m·s−1)间距 S /μm编号
    12000.650A1B1C1
    22000.970A1B2C2
    32001.290A1B3C3
    42500.670A2B1C2
    52500.990A2B2C3
    62501.250A2B3C1
    73000.690A3B1C3
    83000.950A3B2C1
    93001.270A3B3C2
    下载: 导出CSV
  • [1] Niendorf T, Leuders S, Riemer A, et al. Highly anisotro pic steel processed by selective laser melting[J]. Metallurgical & Materials Transactions B, 2013, 44(4): 794 − 796.
    [2] 高倩. 激光选区熔化316L不锈钢微观组织与机械性能研究[D]. 西安科技大学, 2020.

    Gao Qian. Microstructure and mechanical properties of selective laser melting 316L stainless steel[D]. Xi’an University of Science and Technology, 2020.
    [3] 宋亚男. SLM增材制造表面激光抛光机理与工艺研究[D]. 山东大学, 2021.

    Song Yanan. Research on laser polishing mechanism and process of SLM additive manufacturing surface[D]. Shandong University, 2021.
    [4] Liu Binqi, Fang Gang, Lei Liping. An analytical model for rapid predicting molten pool geometry of selective laser melting (SLM)[J]. Applied Mathematical Modelling, 2021, 92: 505 − 524. doi:  10.1016/j.apm.2020.11.027
    [5] Körner C, Attar E, Heinl P. Mesoscopic simulation of selective beam melting processes[J]. Journal of Materials Processing Technology, 2011, 211(6): 978 − 987. doi:  10.1016/j.jmatprotec.2010.12.016
    [6] Wu Y C, Hwang W S, San C H, et al. Parametric study of surface morphology for selective laser melting on Ti6Al4V powder bed with numerical and experimental methods[J]. International Journal of Material Forming, 2018, 11(6): 807 − 813. doi:  10.1007/s12289-017-1391-2
    [7] Wu Y C, San C H, Chang C H, et al. Numerical modeling of melt-pool behavior in selective laser melting with random powder distribution and experimental validation[J]. Journal of Materials Processing Technology, 2018, 254: 72 − 78. doi:  10.1016/j.jmatprotec.2017.11.032
    [8] 向羽, 张树哲, 李俊峰, 等. Ti6Al4V的激光选区熔化单道成形数值模拟与实验验证[J]. 浙江大学学报(工学版)(自然科学版), 2019, 53(11): 2102 − 2109.

    Xiang Yu, Zhang Shuzhe, Li Junfeng, et al. Numerical simulation and experimental verification for selective laser single track melting forming of Ti6Al4V[J]. Journal of Zhejiang University (Engineering Science)(Natural Science Edition), 2019, 53(11): 2102 − 2109.
    [9] 单雨薇. 选择性激光熔化中表层粉末熔池热力学特性3维数值模拟研究[D]. 上海交通大学, 2020.

    Shan Yuwei. 3D numerical modeling of thermodynamic characteristics of surface powder melt pool in selective laser melting[D]. Shanghai Jiao Tong University, 2020.
    [10] 叶唯娟. 激光选区熔化过程中熔池温度场及形貌演化的数值模拟[D]. 西安理工大学, 2019.

    Ye Weijuan. Numeriacl simulation on temperature field and morphology evolution of molten pool during selective laser melting processing[D]. Xi’an University of Technology, 2019.
    [11] Hen Xuehui, Mu Weihao, Xu Xin, et al. Numerical analysis of double track formation for selective laser melting of 316L stainless steel[J]. Applied Physics A, 2021, 127(8): 1 − 13.
    [12] Xiang Y, Zhang S, Wei Z, et al. Forming and defect analysis for single track scanning in selective laser melting of Ti6Al4V[J]. Applied Physics A, 2018, 124(10): 1 − 12.
    [13] Cao L, Yuan X. Study on the numerical simulation of the SLM molten pool dynamic behavior of a nickel-based superalloy on the workpiece scale[J]. Materials, 2019, 12(14): 2272. doi:  10.3390/ma12142272
    [14] 胡国明. 颗粒系统的离散元素法分析仿真: 离散元素法的工业应用与EDEM 软件简介[M]. 武汉: 武汉理工大学出版社出版, 2010.

    Hu Guoming. Analysis and simulation of granular system by discrete element method using EDEM[M]. Wuhan: Wuhan University of Technology Press, 2010.
    [15] Dai Donghua, Gu Dongdou. Influence of thermodynamics within molten pool on migration and distribution state of reinforcement during selective laser melting of AlN/AlSi10Mg composites[J]. International Journal of Machine Tools and Manufacture, 2016, 100: 14 − 24. doi:  10.1016/j.ijmachtools.2015.10.004
    [16] 梁祖磊, 孙中刚, 张少驰, 等. 数值模拟在激光选区熔化中的应用及研究现状[J]. 航空制造技术, 2018, 61(22): 87 − 91. doi:  10.16080/j.issn1671-833x.2018.22.087

    Liang Zulei, Sun Zhonggang, Zhang Shaochi, et al. Application and research status of numerical simulation in laser selective melting[J]. Aeronautical Manufacturing Technology(Natural Science Edition), 2018, 61(22): 87 − 91. doi:  10.16080/j.issn1671-833x.2018.22.087
    [17] 邱云基, 申梦清, 牛晓峰, 等. 选区激光熔化316L粉末过程数值模拟研究[J]. 特种铸造及有色合金, 2021, 41(9): 1092 − 1096. doi:  10.15980/j.tzzz.2021.09.007

    Qiu Yunji, Shen Mengqing, Niu Xiaofeng, et al. Numerical simulation of the melting process of selective laser melting 316L powders[J]. Special Casting & Nonferous Alloys(Natural Science Edition), 2021, 41(9): 1092 − 1096. doi:  10.15980/j.tzzz.2021.09.007
    [18] 曹志伟. 激光选区熔化单熔道成型数值模拟及缺陷分析[D]. 天津理工大学, 2019.

    Cao Zhiwei. Numerical simulation of single track formation and defects analysis in selective laser melting[D]. Tianjin University of Technology, 2019.
    [19] Wang Wenhao, Liu Xinyu, Pan Lu. Effect of linear energy density on bubble-defect of 316L stainless steel by selective laser melting[J]. Key Engineering Materials, 2019, 4834: 32 − 36.
    [20] Shi Xiaojie, Liu Xin, Ren Shuai, et al. Selective laser melting fabricated tungsten with thin-walled structure: role of linear energy density on temperature evolution and manufacturing quality[J]. International Journal of Material Forming, 2022, 15(1): 1 − 11. doi:  10.1007/s12289-021-01647-3
    [21] Gusarov A V, Smurov I. Modeling the interaction of laser radiation with powder bed at selective laser melting[J]. Physics Procedia, 2010, 5: 381 − 394. doi:  10.1016/j.phpro.2010.08.065
    [22] Olakanmi E O, Cochrane R F, Dalgarno K W. A review on selective laser sintering/melting (SLS/SLM) of aluminium alloy powders: Processing, microstructure, and properties[J]. Progress in Materials Science, 2015, 74: 401 − 477. doi:  10.1016/j.pmatsci.2015.03.002
    [23] Yuan W H, Chen H, Cheng T, et al. Effects of laser scanning speeds on different states of the molten pool during selective laser melting: Simulation and experiment[J]. Materials & Design, 2020, 189: 108542.
    [24] Yao N N, Peng X H. The preparation method of metal powder for 3D printing[J]. Sichuan Nonferrous Metals, 2013, 4: 48 − 51.
    [25] Seyda V, Herzog D, Emmelmann C. Relationship between powder characteristics and part properties in laser beam melting of Ti–6Al–4V, and implications on quality[J]. Journal of Laser Applications, 2017, 29(2): 022311. doi:  10.2351/1.4983240
  • [1] 张楠, 张海武, 王淼辉, 杜兵, 张平, 张志豪.  微米级选区激光熔化316L不锈钢拉伸变形中Σ3n特殊晶界的分布 . 焊接学报, 2023, 44(1): 33-39. doi: 10.12073/j.hjxb.20220106003
    [2] 李俐群, 何平, 宫建锋.  铝合金激光-MIG复合焊熔滴对匙孔作用的模拟 . 焊接学报, 2022, 43(8): 1-7. doi: 10.12073/j.hjxb.20220304002
    [3] 徐榕蔚, 张振杰, 刘清原, 张光辉, 龙芋宏.  选区激光熔化制备多孔结构的成形偏差及力学性能与压缩失效分析 . 焊接学报, 2022, 43(10): 49-56. doi: 10.12073/j.hjxb.20211005001
    [4] 巴培培, 董志宏, 张炜, 彭晓.  选区激光熔化成形12CrNi2合金钢的显微组织和力学性能 . 焊接学报, 2021, 42(8): 8-17. doi: 10.12073/j.hjxb.20210323003
    [5] 张宇, 姜云, 胡晓安.  选区激光熔化成形Inconel 625合金的激光焊接头组织及高温蠕变性能 . 焊接学报, 2020, 41(5): 78-84. doi: 10.12073/j.hjxb.20191211001
    [6] 唐琪, 陈静青, 陈鹏, 陈勇, 赵颖.  基于有限元的激光增材过程熔化热积累模拟 . 焊接学报, 2019, 40(7): 100-104. doi: 10.12073/j.hjxb.2019400189
    [7] 刘剑1,樊丁1,陈秀娟2,黄健康1,余淑荣2.  异厚度铝钢电弧辅助激光对接熔钎焊温度场和应力应变场数值模拟 . 焊接学报, 2018, 39(2): -. doi: 10.12073/j.hjxb.2018390036
    [8] 尹燕1,刘鹏宇1,路超2,肖梦智1,3,张瑞华2,3选区激光熔化成形316L不锈钢微观组织及拉伸性能分析 . 焊接学报, 2018, 39(8): 77-81. doi: 10.12073/j.hjxb.2018390205
    [9] 樊丁, 黄自成, 黄健康, 郝珍妮, 王新鑫, 黄勇.  氧元素分布模式与AA-TIG焊熔池形貌的数值模拟 . 焊接学报, 2016, 37(2): 38-42.
    [10] 陈树海, 张茗馨, 李猛, 黄继华.  钢/铝异种金属光纤激光熔钎焊数值模拟 . 焊接学报, 2015, 36(10): 21-24.
    [11] 李美艳, 韩彬, 蔡春波, 王勇, 宋立新.  激光熔覆镍基合金温度场和应力场数值模拟 . 焊接学报, 2015, 36(5): 25-28,32.
    [12] 周广涛, 郭广磊, 方洪渊.  激光诱导焊接温度场的数值模拟 . 焊接学报, 2014, 35(7): 22-26.
    [13] 周一俊, 邓德安, 冯可, 毕涛.  低碳钢薄板单道堆焊焊接变形的数值模拟 . 焊接学报, 2013, (12): 101-104.
    [14] 邢淑清, 郝飞, 闫波, 麻永林.  DP590钢两道次焊接温度场数值模拟与试验验证 . 焊接学报, 2012, (12): 81-84.
    [15] 李慧琴, 刘轶暄, 韩强, 麻永林.  DP590双相钢焊接的数值模拟与试验测试 . 焊接学报, 2011, (10): 9-12.
    [16] 伊鹏, 刘衍聪, 石永军, 姜浩.  铸铁表面激光熔凝行为及温度场数值模拟 . 焊接学报, 2011, (8): 81-84.
    [17] 董文超, 陆善平, 李殿中, 李依依.  焊接电弧与活性组元对TIG焊熔池形貌影响的数值模拟 . 焊接学报, 2009, (11): 49-52,56.
    [18] 齐海波, 杨明辉, 齐芳娟.  扫描路径对电子束选区熔化TC4成形件性能影响的数值模拟 . 焊接学报, 2009, (8): 5-8.
    [19] 张广军, 赵琳琳, 冷雪松.  双钨极氩弧焊焊缝成形的数值模拟 . 焊接学报, 2008, (8): 29-31.
    [20] 徐庆鸿, 郭伟, 田锡, 唐李智.  激光熔覆三维温度场数值模型的建立与验证 . 焊接学报, 1997, (2): 58-62.
  • 加载中
图(9) / 表 (3)
计量
  • 文章访问数:  76
  • HTML全文浏览量:  17
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-04-26
  • 网络出版日期:  2023-03-20
  • 刊出日期:  2023-03-25

选区激光熔化316L扫描特征参数对表面形貌影响

doi: 10.12073/j.hjxb.20220426002
    作者简介:

    朱春霞,教授,博士研究生导师;主要从事3D打印及其相关技术研究;已发表论文40余篇;Email: 11905269@qq.com

  • 中图分类号: TG 401; TH 162; TG 141; TG665

摘要: 为研究不同扫描特征参数组合对选区激光熔化(selective laser melting, SLM)表面形貌的影响,以316L不锈钢粉末为例,进行介观尺度的单层双道数值模拟研究.基于离散元法建立单层的粉床数值模型.使用流体体积法对粉床受热部分粉末的熔化过程中的熔化、流动和凝固过程进行计算.考虑激光功率、扫描速度和扫描间距3个扫描特征参数,设计并进行正交试验,从熔道形貌特征和熔道宽度2个方面研究所选扫描特征参数对成形件表面的熔道形貌影响.依据数值模拟中的参数进行实际打印及形貌观察试验,验证数值模拟的有效性.结果表明,在313 ~ 500 J/m的线能量密度和50 ~ 90 μm的扫描间距范围内,可以得到平整连续局部缺陷少的熔道形貌,且该区间内的参数组合依次线性对应;对熔道形貌的完整性影响由大到小依次为扫描速度、扫描间距和激光功率.

English Abstract

朱春霞, 邱渤杰. 选区激光熔化316L扫描特征参数对表面形貌影响[J]. 焊接学报, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
引用本文: 朱春霞, 邱渤杰. 选区激光熔化316L扫描特征参数对表面形貌影响[J]. 焊接学报, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
ZHU Chunxia, QIU Bojie. Effect of scanning characteristic parameters on surface morphology of selective laser melting 316L[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
Citation: ZHU Chunxia, QIU Bojie. Effect of scanning characteristic parameters on surface morphology of selective laser melting 316L[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2023, 44(3): 114-121, 128. doi: 10.12073/j.hjxb.20220426002
    • 选区激光熔化(selective laser melting, SLM)是激光增材技术的一种,以金属或陶瓷粉末为原料,在激光的熔化作用下快速成形出具有高致密度的零件[1-2].316L是常用的铁基材料因为出色的力学性能与耐腐蚀性等.然而,该工艺下的零件表面常极为粗糙,影响应用[3].鉴于SLM成形过程的迅速性和微小性,不少研究者开始建立激光热源与金属粉末的相互作用数值计算模型—介观尺度模型,研究SLM中工艺参数对熔道形貌特征的影响来合理调节扫描参数,以对成形件表面形貌进行优化[4].

      在介观尺度层面,基于离散元法的激光增材工艺介观数值模拟中,最早由德国Carolin Körner等人[5]于2011年基于二维玻尔兹曼模型实现,可以对粉末熔化再凝固的过程进行预测.然而,这仅是二维模型,无法对成形件的表面的三维形貌进行模拟及实际试验验证.后续的研究者主要焦距于三维化的基于离散元与流体体积法的介观数值模型的合理性.吴宇车等人[6-7]以TC4和H13为例建立SLM工艺下基于离散元法和流体体积法的三维介观数值模型,对在不同功率和速度下的熔道形貌进行研究,证实了数值模拟中加入蒸发反冲力模型的必要性;向羽等人[8]在基于离散元和流体体积法的介观数值模拟中,加入蒸汽反冲力模型的情况下进行TC4的SLM单熔道成形计算,有效得出了成形理想熔道形貌的工艺参数组合;单雨薇[9]对激光粉床增材的介观模拟中的常被采用的多种热源模型进行比较,发现表面热源模型模拟下的熔池形状相比体积热源模型更加具有合理性.叶唯娟[10]从上述被逐渐完善的介观数值模型出发,对两种不同路径的单层多道扫描进行研究,发现有效熔宽的扫描间距和Z形的扫描方式下能得到更为稳定的温度场. Hen等人[11]也在此类介观模型上对SLM下的单层多熔道下的成形件熔道形貌特征进行了研究,但未涉及到不同扫描间距对熔道形貌的影响性.

      可见,当前基于离散元和流体体积法的SLM介观尺度数值模拟,大都局限于与单熔道成形效果有关的扫描工艺参数研究,即使有单层多道的数值相关研究,对不同扫描间距下成形的熔道形貌的关注也较少.SLM的本质是层层累积叠加成形,每一层皆由复数条单熔道扫描搭接而成,不同的扫描间距对于成形件的表面形貌必定存在较大影响,客观上单熔道的数值模拟研究存在一定缺陷.文中也将基于离散元法和流体体积法,从单层双道的角度出发,对316L粉末的成形过程进行瞬态过程的数值模拟及验证,从表面的熔道形貌和熔道宽度两个方面,研究SLM工艺下激光功率、扫描速度和扫描间距3个扫描参数对成形件表面形貌特征的影响.

    • 数值模型的建立要符合基本的工艺过程,SLM的基本工艺过程可分为粉末铺设和激光扫描2个阶段.文中将直接建立粉末铺设完成后的粉床模型和介观层面下的瞬态的熔道成形模型.为能够合理阐明SLM工艺下的诸如热传导、蒸发、凝固、表面张力等效应,在模型建立过程中需要一定简化假设,具体如下:①熔化流动的金属为不可压缩牛顿流体;②所有不锈钢粉末在单个上均为形貌均匀的球体;③不考虑蒸发带来的质量损失且气液界面始终稳定[12-13].

    • 使用离散元模拟工具EDEM,建立单层堆积密度0.4,计算域为1 020 μm × 408 μm × 61.2 μm的长方体的球状单层粉末局部堆积模型,该计算域的选择满足熔道设计长宽、粉末层的最高粉末颗粒高度以及后处理过程中对于成形过程物理现象的观测.球状粉末D50为30 μm且在10 ~ 60 μm范围内呈正态分布.使用只考虑颗粒间的法向与切向作用力Hertz-Mindlin(no slip)模型作为颗粒生成过程中的接触模型[14].

    • SLM工艺的热源输入是激光束,材料因激光所带的自由光子与材料的原子碰撞被加热.在数值模拟中[15-16],激光通常被定义为移动的高斯热源,自由光子的热量较低,易被反射难以穿透到材料表面以下较深的位置,因此在基于激光束的加热模式中,采用平面热源模型更具有合理性.

      $$ Q\left( r \right) = \frac{{2AP}}{{{\text{π}} {R^2}}}\exp \left( {\frac{{ - 2{r^2}}}{{{R^2}}}} \right) $$ (1)

      式中:A为粉末的热量的吸收率;P为激光功率;R表示半径为50 μm的光斑其数据源于设备说明;r为熔化区域任一点到激光中心距离.

      在SLM成形过程中,热量交互行为主要有热传导,热对流、热辐射和蒸发散热.基本传热方程为

      $$ \rho c\frac{{\partial T}}{{\partial Z}} = Q\left( r \right) - {q_{\rm{c}}} - k\nabla T - {q_{{\rm{e}}{\rm{v}}}} $$ (2)

      式中:$ \rho $为材料密度;$ c $为材料比热;$ T $为初始温度;$ k $为材料热导率;${q}_{{\rm{c}}}$为气体的对流换热;$ k\nabla T $为表面热辐射.

      在熔池的自由界面与上方保护气体间存在的热量交换边界条件为

      $$ {q_{{\rm{in}}}} = q - {h_{\rm{c}}}\left( {T - {T_0}} \right) - {\varepsilon _{\rm{r}}}{\sigma _{\rm{s}}}\left( {{T^4} - {T_0}^4} \right) - {q_{{\rm{ev}}}} $$ (3)

      式中:$ {q_{{\rm{in}}}} $为熔池自由界面上的热流变化;$ {h_{\rm{c}}} $为自然对流传热系数;$ {\varepsilon _{\rm{r}}} $为发射率;$ {\sigma _{\rm{s}}} $为斯特凡−玻尔兹曼常数;$ T $为熔池上表面温度;$ {T_0} $为环境温度.

      金属粉末熔化发生蒸发散热控制方程为

      $$ {q_{{\rm{ev}}}} = 0.82\frac{{\Delta H^* {P_0}}}{{\sqrt {2{\text{π}} MRT} }}{e^{\left( {\Delta * \tfrac{{T - {T_{{\rm{lv}}}}}}{{RT{T_{{\rm{lv}}}}}}} \right)}} $$ (4)

      式中:$ M $为材料摩尔质量;$ R $为理想气体常数;$ {P}_{0} $为环境压力;$ {T}_{\mathrm{l}\mathrm{v}} $为材料沸点;$ \Delta {H}^{\mathrm{*}} $为损失金属蒸气有效焓.

      为简化研究模型和提高数值计算可靠性,文中未考虑热力耦合问题.由于热量分布不均涉及到的热变形涉及高度非线性的热弹塑性问题,同时残余应力的试验测试精度也难以得到准确验证.

    • 粉末材料在热源范围内被加热至熔点以上形成熔池,熔池上部和下部的密度以及温度不同使得内部存在自然对流,其流动遵循N-S方程,在对流上主要考虑重力,表面张力以及蒸汽反冲力对熔池稳定性的影响[17].

      表面张力在熔池上主要体现为Marangoni效应,即

      $$ \left\{ \begin{gathered} {f_{\rm{s}}} = \sigma \kappa n\left| {\nabla \alpha } \right|\frac{{2\rho }}{{{\rho _{316{\rm{L}}}} + {\rho _{\text{α}} }}} \\ {f_{\rm{m}}} = \frac{{d\sigma }}{{dT}}\left[ {\nabla T - \left( {n \cdot \nabla T} \right)n} \right]\left| {\nabla \left( {1 - \alpha } \right)} \right|\frac{{2\rho }}{{{\rho _{316{\rm{L}}}} + {\rho _{\text{α}} }}} \\ \end{gathered} \right. $$ (5)

      式中:${f}_{{\rm{s}}}$为熔池表面张力;$ \sigma $为表面张力系数;n为熔池表面单位法向量;$ \kappa $表示曲率;$ \left|\nabla \alpha \right| $ 将界面力转化为体积表面力;${{f}}_{{\rm{m}}}$为熔池中Marangoni力.

      当熔池最高温度超过金属沸点时,金属在蒸发溢出同时对熔池表面产生一个反冲压力,亦称蒸汽反冲力,表达式为

      $$ {p_{\rm{r}}} = 0.54{P_0}\exp \left( {{L_{{\rm{lv}}}}\frac{{T - {T_{{\rm{lv}}}}}}{{RT{T_{{\rm{lv}}}}}}} \right) $$ (6)

      在FLOW-3D软件中,对传热与流动相关的控制方程进行求解;使用流体体积法对熔池流动及凝固行为实现跟踪.

    • 在SLM中,材料的部分热物参数难以通过试验手段获取,文中将材料的热物属性看作仅是与温度有关的变量,采用JMATPRO软件进行计算得到材料的密度、比热、粘度和热导率的热物理参数,如图1所示.

      图  1  316L热物参数

      Figure 1.  Theromphysicl parameters of 316L. (a) liquid viscosity; (b) specific heat; (c) thermal conductivity; (d) density

      在基于有限元的数值模拟中,需要考虑同种材料的粉末情况与实体块情况下的导热差异等物性参数.文中的数值模拟研究模型是基于离散元法的建模,不再需考虑粉体与块体材料的物理差异,表1为材料的参数设置[18].

      表 1  316L材料参数

      Table 1.  316L material parameters

      固相密度ρ/(kg·m−3) 固相温度T/K 液相温度T/K 沸点温度T/K 熔化潜热Hf /(J·kg−1) 蒸发潜热Hv /(J·kg−1) 表面张力r /(N·m−1) 吸收率A
      7 850 1 658 1 723 3 090 2.7 × 105 7.45 × 106 1.6 0.25

      试验涉及参数自变量多,为提高数值计算效率,设置三水平三因素正交试验,见表2,不同水平按照经验范围进行选定.

      表 2  正交试验设计

      Table 2.  Experimental design using orthogonal table

      参数变量A功率P/WB速度v/(m·s−1)C间距S/μm
      水平12000.650
      水平22500.970
      水平33001.290

      选区激光熔化的本质是若干条预定扫描轨迹下形成多条熔道的搭接集合,为简化数值模型的复杂度和运算时间,采用同扫描层相邻的双熔道数值模拟的特征效果来分析实际成形件上表面可能出现的熔道形貌特征.为方便后续分析,对正交试验不同组合进行编号,见表3.

      表 3  试验参数及编组

      Table 3.  Experimental parameters and grouping

      组号功率 P/W速度 v /(m·s−1)间距 S /μm编号
      12000.650A1B1C1
      22000.970A1B2C2
      32001.290A1B3C3
      42500.670A2B1C2
      52500.990A2B2C3
      62501.250A2B3C1
      73000.690A3B1C3
      83000.950A3B2C1
      93001.270A3B3C2
    • 图2为数值模拟的正交试验设计下的9组单层双道形貌图.图2a2d2g2h这4组参数组合下的熔道成形效果更为接近理想状态的熔道形貌,如图3所示.这4组参数下的双道形貌,搭接部分平滑,未熔合缺陷特征少且熔道平直,是良好的的表面形貌.

      图  2  数值模拟熔道形貌

      Figure 2.  Numerical simulation of melt tracks morphology. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2

      图  3  理想熔道形貌示意

      Figure 3.  Ideal melt track morphology

      图2e虽然熔道整体平直,但有明显岛屿球化特征;图2f图2i参数组合,激光功率更大,虽然粉床受热源照射部分熔化效果好,但扫描速度较快,也存在明显未熔合缺陷和岛屿球化效应;图2b图2c下的参数组合成形的熔道形貌存在明显缩颈、熔道小角度扭曲和未熔合特征,相邻熔道的搭接不明显,激光功率低,扫描速度较快,粉末熔化时较小的熔池尺寸导致湿润角增大,在表面张力作用下极易发生岛屿球化形成不连续的熔道,不生成规律性纹路,最终形成缺陷特征多的熔道形貌.

    • 研究中所涉及的扫描参数包括扫描速度、激光功率和扫描间距.涉及参数较多且参数之间对于熔道形貌特征的影响存在一定耦合.

      为更好评估各扫描参数与熔道良好形貌特征之间的关系,引入综合参数激光线性能量密度(line energy density,LED),即激光功率与扫描速度的比值,评估激光的能量输入[19-20].

      图4所示,成形效果良好状态表面的熔道形貌特征需要较高的线能量密度和扫描间距的参数组合,这4组参数组合从左往右依次是1,8,4,7组,对应的线能量密度依次为333,333,416,500 J/m,对应扫描间距为0.05,0.05,0.07,0.09 mm.其中1组和8组在图中数据重合.在非良好熔道形貌区域中,线能量密度较低,在线能量密度较低的情况下,由于Plateau-rayleigh不稳定性,熔体趋向于从圆柱状向球状变化以减少其表面能,导致出现球化、不连续的熔道等不良的表面形貌特征的风险增加.虽然也有研究表示,当线能量密度较高时,熔池面积变大且表面会受到较大的反冲作用,凝固后的熔道表面粗糙[21].文中的数值结果中并无相关体现,这可能是受材料、堆积密度、参数设置、统计方法等因素影响下造成的结果误差.

      图  4  线能量密度(LED)与扫描间距对熔道形貌特征影响

      Figure 4.  Influence of linear energy density (LED) and scanning spacing on the morphology characteristics of the melt tracks

      图4的良好熔道形貌区域可知,扫描间距与线能量密度在所设置的数值模拟参数范围内满足某种线性关系,其表达式为

      $$ {x_1} = 0.239\;5x - 29.719 $$ (7)

      式中:$ {x_1} $为扫描间距,单位 μm;$ x $为线能量密度,单位J/m.当线能量密度和扫描间距满足该式时,可以得到较好的熔道形貌.

      对以上4组参数组合进行各参数的变化分析:7组和8组,功率一致,扫描速度上升50%,扫描间距下降44%;对比1组和4组,在速度一致情况下,第4组功率提升25%,间距提升40%;对比1组和8组,在间距一致情况下,第8组功率提升33%,速度提升50%.说明在影响熔道形貌完整度上,扫描速度、扫描间距和扫描功率是依次递减的. 同时,高的扫描速度下液相时间缩短,有利于抑制球化现象,提高致密度[22-23].

    • 图5所示,在所有9组计算结果中,熔道的宽度基本随着扫描间距和线能量密度的增加而增加,且相比扫描间距,线能量密度对熔宽的影响更大.在50 μm间距下,有2组数据点高度重合,说明在不同功率与速度的相同比值下的线能量密度对于熔宽影响微乎.

      图  5  不同扫描间距下的线能量密度与熔宽

      Figure 5.  Line energy density and melting width under different scanning distances

      具有良好熔道形貌的4组参数组合均在3组数据曲线的最末端,可显著看出线能量密度和扫描间距对熔宽的影响也具备线性的特征.

    • 为保证对数值模拟结果的可靠性,由表3参数组号进行实际试验,通过打印样件,观察成形件表面熔道形貌与非直接的测量熔道宽度的方式与数值结果进行对比,以验证数值模拟结果的有效性.

    • 使用易博三维公司研发的IGAM-2金属打印机和商业化316L粉末分别对9组不同的扫描参数组合进行单个试样打印,样件为长宽高为7 mm × 7 mm × 5 mm的长方体块.

      局限于设备因素,无法单独在每个试件表面进行指定熔道数量的扫描及作断面切割观测.双道的本质是多道的局部体现且不涉及成形后合金组织形貌研究,因此在结果观察时将直接观察成形件最表层的形貌特征.使用超景深三维显微系统VHX-1000E进行对最表层熔道形貌图像和熔道宽度数据进行采集,图6为超景深系统下的316L粉末图像.

      图  6  316L粉末形状轮廓图像

      Figure 6.  316L powder shape outline image

      在基于粉床的金属增材工艺中,粉末须满足粒径细小、粒度分布窄、球形度高等要求.从图6可以看出,粉末基本为球形.经激光粒度测试,该批316L不锈钢金属粉末的粒径分布为D10 = 24.7 μm,D50 = 35.5 μm,D90 = 50.7 μm,其中平面平均径D[3,2]与体积平均径D[4,3]分别为34.1,36.7 μm,说明该批粉末形状较为规则,粒径较为集中.图7为粉末的粒径分布图像,基本符合正态分布,这有利于粉层颗粒分布的均一性和粉末熔化与铺设过程中良好的流动性[24-25].

      图  7  316L粉末粒径分布

      Figure 7.  316L powder particle size distributio

    • 图8为试验后不同成形件在超景深显微镜下的表面形貌图像,其中各个分图的数字和字母编号对应表3中的正交设计表.

      图  8  不同参数组合下的成形件表面形貌

      Figure 8.  Forming surface topography under the different parameters combination conditions. (a) group No. 1 No. A1B1C1; (b) group No. 2 No. A1B2C2; (c) group No. 3 No. A1B3C3; (d) group No. 4 No. A2B1C2; (e) group No. 5 No. A2B2C3; (f) group No. 6 No. A2B3C1; (g) group No. 7 No. A3B1C3; (h) group No. 8 No. A3B2C1; (i) group No. 9 No. A3B3C2

      当激光功率为200 W时,在不同扫描速度和扫描间距下,成形件表面的熔道形貌如图8a ~ 8c所示.图8a与数值结果的熔道与实际结果具有较高一致性,但也存在少量球化特征,试件表面熔道并不完全平整,这可能与基板平台与刮刀的定位精度和平整度有关.图8b中熔道存在明显缩颈及熔道小角度扭曲状态,且同层熔道之间的搭接线不再明显,所对应的数值模拟结果中并未出现熔道不连续状态,可能是受限于模拟的熔道长度过短而没有展现.在图8c中,熔道特征难以辨认,大量岛屿球化和不连续缺陷.图8a ~ 8c的表面熔道形貌特征结果与数值模拟基本一致.

      当激光功率为250 W时,对应不同扫描速度和扫描间距的成形件表面形貌如图8d ~ 8f所示.相比图8a8d的表面更为平坦,在更高的激光功率下,粉末的熔化效果更好.图8e中由于扫描间距过大,相邻熔道搭接线不明显,较快的扫描速度使得熔化效果不足,熔道不连续且有弯曲与岛屿球化现象.图8f相比图8e有更为平直连续的熔道,这是更高的扫描速度下较小的扫描间距可以有效弥补前一条熔道的因扫描过程中能量密度不足而造成的未熔合特征.

      图8e图8f所对应的数值结果与试验验证结果略有出入,在数值模拟结果中,图8e中的熔道形貌相比于图8f更为平直连续且不存在明显缩颈;但在实际打印结果中,图8f中的熔道形貌略比图8e中平直连续,这说明数值模型在有在细节上有待完善的地方:一是模拟的熔道长度仅有几百微米长的区域,在范围尺度上可能缺少一定广泛性;二是数值模型的成形底面是假设平整的基板平面,实际成形时,仅有第一层是在较为平整的基板面上进行,其余层均是在上一层具有相对于基板平面存在高低起伏的面上进行粉末层铺设、扫描、熔化凝固过程.

      当激光功率为300 W时,对应不同扫描速度和扫描间距的成形件表面形貌如图8g ~ 8i所示. 图8g图8h中,在低扫描速度下(0.6 m/s和0.9 m/s),熔道平直完整,球化特征少.即使在0.6 m/s的低扫描速度下搭配较大的扫描间距,由于较高的能量密度输入,熔池面积大、粘度低,有利于熔池凝固前的铺展成形,熔道成形效果优于低扫描间距(50 μm),高扫描速度(0.9 m/s)参数组合.图8i中,成形件表面明显熔道不连续,扭曲现象严重且伴随大量岛屿球化特征,与图8c的特征基本一致,但图8i下的熔道比图8c中更为连续,这说明在高的扫描速度下,大幅提升激光功率和适当减小扫描间距对成形件表面形貌完整度的提升不太明显.

    • 在实际成形的熔宽方面,在共聚焦系统下采集到的较为完整的单道和双道的熔宽数据和数值结果所测量结果对比如图9所示.

      图  9  熔道宽度及误差

      Figure 9.  Width and error of the melt track

      不同线能量密度和扫描间距下,无论是数值模拟的计算结果还是试验结果,熔道宽度变化趋势基本一致,平均相对误差分别为单道15.91%,双道20.98%,由于多种客观原因比如刮刀的不平整、基板平面校准的误差、光路系统定位与聚焦平面的准确性等,试验与仿真之间存在误差是不可避免的,其中最大的误差体现为在9组参数组合中存在2组在表面的熔道形貌成形效果不佳的参数组合的缺陷特征的体现上.其中,数值结果下的表面熔道特征与实际结果误差较小的1,4,7,8组双道成形及搭接效果良好,熔道平整连续,搭接区域过渡平整,与数值模拟部分的结论吻合.

      综上,文中的数值模拟可以为SLM工艺下不同激光功率、扫描间距和扫描速度参数组合下的成形件的表面的熔道形貌特征提供有效参考.

    • (1)对于成形效果良好表面的熔道形貌,在线能量密度313 ~ 500 J/m与扫描间距50 ~ 90 μm之间存在明显的线性关系,虽然这种线性关系可能是参数梯度设置不够小的结果.

      (2)扫描速度、扫描间距和激光功率对成形理想表面形貌的影响由大到小依次为扫描速度、扫描间距、扫描功率.

      (3)熔道宽度增大与减小主要受线能量密度影响.

参考文献 (25)

目录

    /

    返回文章
    返回