高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

Zr65Al7.5Ni10Cu17.5非晶合金激光熔凝的热效应模拟

陈星 葛亚琼

引用本文:
Citation:

Zr65Al7.5Ni10Cu17.5非晶合金激光熔凝的热效应模拟

    作者简介: 陈星(1994-), 男, 硕士研究生, 主要从事金属材料激光加工的研究.
    通讯作者: 葛亚琼, geyaqiong@tyust.edu.cn
  • 基金项目:

    国家自然科学基金资助项目 51405324

    太原科技大学博士科研启动基金资助项目 20152010

    山西省自然科学基金面上青年基金资助项目 201701D221068

  • 中图分类号: TG665

Simulation of thermal effect of Zr65Al7.5Ni10Cu17.5 amorphous alloy by laser melting

    Corresponding author: GE Yaqiong, geyaqiong@tyust.edu.cn
  • CLC number: TG665

  • 摘要: 为了研究激光增材制备非晶合金过程中熔池和热影响区的成形机制, 利用有限元软件ANSYS, 对激光增材制造技术的基础过程-激光快速熔凝Zr65Al7.5Ni10Cu17.5非晶合金的热效应进行了数值模拟分析。结果表明, 激光单点熔凝时, 熔池的平均冷却速率为6.3×104K/s, 热影响区的平均冷却速率为1.4×104K/s, 远高于Zr65Al7.5Ni10Cu17.5非晶合金的临界冷却速率1.5K/s, 说明激光单点熔凝的热变化满足非晶合金的生长条件; 激光单道熔凝过程中, 熔池的平均冷却速率仍比较高, 为2.11×102K/s, 但热影响区的平均冷却速率较低, 为74K/s, 且热影响区会产生弛豫累积, 可能造成一定程度的晶化。此研究为激光增材制备非晶合金材料提供热效应的理论基础。
  • Figure 1.  Finite element model

    Figure 2.  Simulation of temperature distribution of Zr65 amorphous alloy by laser single point melting at different time

    a—after 3ms irradiation b—after 4ms irradiation

    Figure 3.  Time-temperature curves at A, B and C shown in Fig. 2a

    Figure 4.  Temperature-time curves of side and bottom of molten pool with the same peak temperature of 777K

    Figure 5.  Temperature profile of longitudinal section after 2s laser irradiation

    Figure 6.  Cross section temperature cloud map after 2s laser irradiation

    Figure 7.  Time-temperature curves at E, D, F and G shown in Fig. 6

    Table 1.  Material property parameters of amorphous alloys

    tempera-ture/K dencity/(kg·m-3) specific heat/(J·kg-1·K-1) coefficient of heat conduction/(W·m-1·K-1) coefficoent of convection/(W·mm-2·K-1)
    293.15 7141 286 8.1 62.3
    373.15 7127 295 9.2 62.3
    573.15 7088 317 12.1 62.3
    873.15 7024 369 16.6 62.3
    1173.15 6954 210 21.4 62.3
    1473.15 6473 342 28.7 62.3
    1773.15 6302 383 32.5 62.3
    下载: 导出CSV
  • [1]

    PAULY S, SCHRICKER C, SCUDINO S, et al. Processing a glass-forming Zr-based alloy by selective laser melting[J]. Materials & Design, 2017, 135:133-141. 
    [2]

    YE X, BAE H, SHIN Y C, et al. In situ synthesis and characterization of Zr-based amorphous composite by laser direct deposition[J]. Metallurgical and Materials Transactions, 2015, A46(9):4316-4325. 
    [3]

    PAULY S, LBER L, PETTERS R, et al. Processing metallic glasses by selective laser melting[J]. Materials Today, 2013, 16(1/2):37-41.
    [4]

    JUNG H Y, CHOI S J, PRASHANTH K G, et al. Fabrication of Fe-based bulk metallic glass by selective laser melting: A parameter stu-dy[J]. Materials & Design, 2015, 86:703-708.
    [5]

    LI X P, ROBERTS M P, O'KEEFFE S, et al. Selective laser melting of Zr-based bulk metallic glasses: Processing, microstructure and me-chanical properties[J]. Materials & Design, 2016, 112:217-226.
    [6]

    ZHANG Y, LIN X, WANG L, et al. Microstructural analysis of Zr55Cu30Al10Ni5 bulk metallic glasses by laser surface remelting and laser solid forming[J]. Intermetallics, 2015, 66:22-30. doi: 10.1016/j.intermet.2015.06.007
    [7]

    SCHMIDT M, MERKLEIN M, BOURELL D, et al. Laser based additive manufacturing in industry and academia[J]. CIRP Annals—Manufacturing Technology, 2017, 66(2):561-583. doi: 10.1016/j.cirp.2017.05.011
    [8]

    LU Y, ZHANG H, LI H, et al. Crystallization prediction on laser three-dimensional printing of Zr-based bulk metallic glass[J]. Journal of Non-Crystalline Solids, 2017, 461:12-17. doi: 10.1016/j.jnoncrysol.2017.01.038
    [9]

    GU D D, MA Ch L, XIA M J, et al. A multiscale understanding of the thermodynamic and kinetic mechanisms of laser additive manufacturing[J]. Engineering, 2017, 3(5):675-684. doi: 10.1016/J.ENG.2017.05.011
    [10]

    LIU Sh B, LIU J Ch, QI L J, et al. Numerical simulation of laser remelting process based on smoothed particles hydrodynamics method[J]. Transactions of Materials and Heat Treatment, 2017, 38(1):184-190(in Chinese). 
    [11]

    GAN Y. Study on microstructure and property of Zr-Al-Ni-Cu amorphous material prepared by laser additive manufacturing[D]. Taiyuan: Taiyuan University of Technology, 2015: 23-29(in Chinese).
    [12]

    CAI Ch B, LI M Y, HAN B, et al. Numerial simulation iron-based cladding coating with wide-band laser at different preheating tempe-ratures[J]. Applied Laser, 2017, 37(1):66-71(in Chinese).
    [13]

    TSENG W C, AOH J N. Simulation study on laser cladding on preplaced powder layer with a tailored laser heat source[J]. Optics & Laser Technology, 2013, 48(11):141-152. 
    [14]

    HU Q, LIN X, YANG G L, et al. Crystallization behavior of Zr55Al10Ni5Cu30 amorphous alloys with different morphologies and thermal history conditions[J]. Acta Metallurgica Sinica, 2012, 48(12):1467-1473(in Chinese). doi: 10.3724/SP.J.1037.2012.00433
    [15]

    YANG G L, LIN X, HU Q, et al. Effect of specimen temperature on crystallization during laser remeltitng Zr55Cu30Al10Ni5 bulk metallic glass[J]. Acta Metallurgica Sinica, 2013, 49(8):925-931(in Ch-inese). doi: 10.3724/SP.J.1037.2013.00166
    [16]

    LI G, JIA M D, LIU L, et al. Microstructure and properties of laser cladding Ni40Zr30Ta28 amorphous alloy coating[J]. Laser Technology, 2010, 34(3):306-308(in Chinese).
    [17]

    HE Y Y, LIU Y J, CHEN M, et al. Study on the law of laser scanning rate and microstructure and properties of cladding layer[J]. Laser Technology, 2019, 43(2):201-204(in Chinese).
    [18]

    XUA H D, LUA Y Zh, et al. Laser 3-D printing of Zr-based bulk metallic glass[J]. Journal of Manufacturing Processes, 2019, 39:102-105. doi: 10.1016/j.jmapro.2019.02.020
    [19]

    DI O, NING L, LIN L. Structural heterogeneity in 3-D printed Zr-based bulk metallic glass by selective laser melting[J]. Journal of Alloys & Compounds, 2018, 740:603-609. 
    [20]

    DI O, NING L, WEI X, et al. 3-D printing of crack-free high strength Zr-based bulk metallic glass composite by selective laser melting[J]. Intermetallics, 2017, 90:128-134. doi: 10.1016/j.intermet.2017.07.010
  • [1] 李明海柳爱国宋耀祖 . 激光放大介质温度场和热应力场的数值模拟. 激光技术, 2002, 26(2): 86-89.
    [2] 高立张永康 . 镍-钴合金杆料激光冲击强化的实验和数值模拟. 激光技术, 2006, 30(5): 507-510.
    [3] 王琪琪任乃飞任旭东 . GH4037镍基高温合金激光打孔相变过程数值模拟. 激光技术, 2018, 42(6): 764-768. doi: 10.7510/jgjs.issn.1001-3806.2018.06.007
    [4] 花银群蔡峥嵘陈瑞芳姜辉吉光 . TC4钛合金激光搭接冲击强化的实验和数值模拟. 激光技术, 2010, 34(5): 632-635,639. doi: 10.3969/j.issn.1001-3806.2010.O5.015
    [5] 禹烨牛燕雄王秀生刘杰姜楠 . 强激光稳态热晕效应的数值模拟研究. 激光技术, 2007, 31(2): 182-184.
    [6] 陶向阳季小玲吕百达 . 光束控制系统中热效应的数值模拟. 激光技术, 2003, 27(6): 514-516,596.
    [7] 赵臻龙芋宏黄宇星张光辉周辽蔺泽焦辉黄平 . 水束导引高功率激光的热效应分析. 激光技术, 2023, 47(5): 672-677. doi: 10.7510/jgjs.issn.1001-3806.2023.05.015
    [8] 陈子伦侯静姜宗福 . 高功率掺镱双包层光纤激光器热效应理论研究. 激光技术, 2007, 31(5): 544-547,550.
    [9] 熊征曾晓雁 . 化学反应热效应对激光熔凝区几何特征的影响. 激光技术, 2007, 31(5): 462-464,468.
    [10] 高俊超朱长虹李正佳 . 固态激光介质的热效应与光泵浦极限研究. 激光技术, 2004, 28(3): 271-274.
    [11] 李隆史彭李东亮甘安生白晋涛 . 高功率CO2激光器输出窗热效应的研究. 激光技术, 2004, 28(5): 510-513.
    [12] 陈子伦马厉克姜宗福曹涧秋 . LD端面抽运激光介质热效应的有限元分析. 激光技术, 2005, 29(5): 543-545.
    [13] 冯祝万云芳 . LD端面抽运Nd:GGG激光器热效应研究. 激光技术, 2014, 38(3): 360-363. doi: 10.7510/jgjs.issn.1001-3806.2014.03.016
    [14] 冯杰范宗学单常亮魏欣芮吴琴杨永佳周自刚 . 基于热效应飞秒激光诱导LiNbO3表面结构的研究. 激光技术, 2015, 39(6): 869-872. doi: 10.7510/jgjs.issn.1001-3806.2015.06.029
    [15] 李隆董武威史彭甘安生许启明 . 高功率Yb:YAG微片激光器热效应研究. 激光技术, 2010, 34(1): 8-12. doi: 10.3969/j.issn.1001-3806.2010.01.003
    [16] 黄云火黄磊张海涛柳强闫平巩马理 . 水柱导引抽运光纤激光器端面热效应分析. 激光技术, 2009, 33(3): 225-227.
    [17] 徐仰彬凌亚文 . 复合YVO4-Nd:YVO4激光晶体的超高斯热效应研究. 激光技术, 2009, 33(1): 80-82.
    [18] 崔莉军陈星明杨美霞代俊刘文兵黄守彬 . 高能激光系统中不同材料合束镜的热效应分析. 激光技术, 2024, 48(2): 159-165. doi: 10.7510/jgjs.issn.1001-3806.2024.02.003
    [19] 胡席远胡伦骥熊建钢刘建华骆红陈祖涛 . 激光焊接中聚焦透镜热效应研究. 激光技术, 1998, 22(3): 188-190.
    [20] 刘顺洪吉巧杰扬晶 . 钢管激光弯曲成形的数值模拟. 激光技术, 2006, 30(4): 355-359.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  5902
  • HTML全文浏览量:  4580
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-03-27
  • 录用日期:  2019-05-28
  • 刊出日期:  2020-03-25

Zr65Al7.5Ni10Cu17.5非晶合金激光熔凝的热效应模拟

    通讯作者: 葛亚琼, geyaqiong@tyust.edu.cn
    作者简介: 陈星(1994-), 男, 硕士研究生, 主要从事金属材料激光加工的研究
  • 太原科技大学 材料科学与工程学院, 太原 030024
基金项目:  国家自然科学基金资助项目 51405324太原科技大学博士科研启动基金资助项目 20152010山西省自然科学基金面上青年基金资助项目 201701D221068

摘要: 为了研究激光增材制备非晶合金过程中熔池和热影响区的成形机制, 利用有限元软件ANSYS, 对激光增材制造技术的基础过程-激光快速熔凝Zr65Al7.5Ni10Cu17.5非晶合金的热效应进行了数值模拟分析。结果表明, 激光单点熔凝时, 熔池的平均冷却速率为6.3×104K/s, 热影响区的平均冷却速率为1.4×104K/s, 远高于Zr65Al7.5Ni10Cu17.5非晶合金的临界冷却速率1.5K/s, 说明激光单点熔凝的热变化满足非晶合金的生长条件; 激光单道熔凝过程中, 熔池的平均冷却速率仍比较高, 为2.11×102K/s, 但热影响区的平均冷却速率较低, 为74K/s, 且热影响区会产生弛豫累积, 可能造成一定程度的晶化。此研究为激光增材制备非晶合金材料提供热效应的理论基础。

English Abstract

    • 块体非晶合金有着优异的力学性能、良好的物理、化学、电化学性能及良好的加工性能等,可应用于航空航天、精密机械制造、生物医学等诸多领域,但其传统制备方法——如铜模铸造法、喷射吸铸法等有着明显的尺寸缺陷,限制了非晶合金应用的发展。因此,寻找能够突破尺寸限制的制备方法,是研究块体非晶合金的一个重要课题[1]

      激光增材制造技术有着快速升温-冷却的特点,满足了非晶态组织产生的外部条件,这为探索块体非晶合金的制备方法提供了新思路。YE等人[2]制备了具有高非晶相体积分数且无裂纹的Zr65Al10Ni10Cu15块体非晶合金复合材料。PAULY等人[3]通过激光选区熔化(selective laser melting, SLM)技术制备了具有复杂几何结构的Fe74Mo4P10C7.5B2.5Si2块体非晶合金。随后,JUNG等人[4]研究了SLM工艺参量对微观组织演变和材料热、磁性能的影响。LI等人[5]利用优化后的线扩散函数(line spread function, LSF)技术成功制备出了尺寸达35mm×5mm的致密Zr52.5Ti5Cu17.9Ni14.6Al10块体非晶合金齿轮零件。ZHANG等人[6]研究了不同粒度粉末在LSF中的晶化规律。晶化率是目前激光增材制备块体非晶合金的主要研究内容之一,在块体非晶的制备过程中,复杂的热历史和热影响区结构弛豫的累积是导致晶化率上升的主要原因[7]

      激光熔凝作为激光增材制造的基础过程,有着与其相似的热循环过程,因此研究激光熔凝Zr65Al7.5Ni10Cu17.5(以下简称Zr65)非晶合金,对探究激光增材制备过程Zr65非晶合金中的热效应、晶化规律以及组织演变规律有着重要的意义。而激光熔凝Zr65非晶合金的过程中,高能量激光的高升温和冷却速率,使得熔池存在时间很短且尺寸很小,无法在实验中直接观测[8]。本文中选用数值模拟的方法对激光增材制备的基础过程-激光快速熔凝Zr65非晶合金的温度场进行模拟[9-10],分析激光单点熔凝和单道熔凝非晶合金过程中的温度场变化,计算熔池和热影响区的冷却速率,探讨激光熔凝非晶合金过程中不同部位的晶化问题。

    • 本文中设计了激光单点熔凝Zr65非晶合金和激光单道熔凝Zr65非晶合金两个模拟实验。分析激光单点和单道熔凝非晶合金过程中熔池和热影响区的热历史和组织变化。重点模拟了在单道熔凝过程中,后置位的熔池成形对先置位熔池和热影响区的热效应,研究再升温引起的非晶合金的晶化效应。

      为了方便模拟计算,在一定程度上对模型提出以下假设:(1)假设材料为各向同性;(2)忽略熔池流体的流动作用;(3)忽略材料的汽化作用;(4)忽略材料相变潜热。

      模型的建立如下。激光功率P=3200W,光斑直径D=3mm,对流换热加载在除激光辐照面的其余表面上。对于单点熔凝实验,模拟总时长为1s,热源加载3ms,然后自然空冷。对于单道熔凝实验,模拟总时长为6s,激光扫描速率v=180mm/min,激光行走18mm。

    • 图 1是用DesignModeler软件建立的有限元模型。几何尺寸为50mm×20mm×5mm,越靠近热源位置的网格越细。模拟计算采用的Zr65非晶合金材料性能如表 1所示[11]

      Figure 1.  Finite element model

      Table 1.  Material property parameters of amorphous alloys

      tempera-ture/K dencity/(kg·m-3) specific heat/(J·kg-1·K-1) coefficient of heat conduction/(W·m-1·K-1) coefficoent of convection/(W·mm-2·K-1)
      293.15 7141 286 8.1 62.3
      373.15 7127 295 9.2 62.3
      573.15 7088 317 12.1 62.3
      873.15 7024 369 16.6 62.3
      1173.15 6954 210 21.4 62.3
      1473.15 6473 342 28.7 62.3
      1773.15 6302 383 32.5 62.3
    • 采用高斯热源模型,材料表面激光辐照的功率密度为[12]:

      $ F(r)=\frac{P}{\pi w^{2}} \exp \left(-\frac{r^{2}}{w^{2}}\right) $

      (1)

      式中,P为入射激光功率,w为高斯形式的光束半径,r为其余点离加热斑点中心的距离。

      材料表面受激光辐照处的温度分布为[13]

      $ \rho c \frac{\partial T}{\partial t}=\frac{\partial}{\partial x}\left(\kappa \frac{\partial T}{\partial x}\right)+\frac{\partial}{\partial y}\left(\kappa \frac{\partial T}{\partial y}\right)+\frac{\partial}{\partial z}\left(\kappa \frac{\partial T}{\partial z}\right)+Q $

      (2)

      式中,κ为传导率,ρ为材料密度,c为材料比热容,T为温度,t为时间,Q为加热速率。

    • 图 2是激光单点熔凝Zr65合金分别在3ms(见图 2a)和4ms(见图 2b)的温度场。可以看出,受激光光斑中心辐照区的温度最高,距离光斑中心越远,温度越低;在激光辐照3ms时,材料表面中心最高温度达到2210.05K;4ms时,由于在3ms时激光已停止辐照,材料表面温度开始迅速下降。以升温速率为10K/s的Zr65非晶合金的组织转变温度作参考[14],将峰值温度高于Tm的区域定义为熔池,峰值温度在TmTg之间的区域定义为热影响区,峰值温度低于Tg的区域定义为热稳定区,激光单点熔凝实验中,熔池宽为3.17mm,熔深为0.24mm。

      Figure 2.  Simulation of temperature distribution of Zr65 amorphous alloy by laser single point melting at different time

      选取了图 2aA, BC 3点绘制激光单点熔凝过程的温度-时间曲线(见图 3),A点和B点位于熔池,C点位于热影响区。图中可知,越靠近激光光斑中心的位置,升降温速率越高。计算得A点从TmTg的平均冷却速率为6.5×104K/s,B点的平均冷却速率为6.3×104K/s,C点的平均冷却速率为1.4×104K/s。Zr65非晶合金的临界冷却速率为1.5K/s,铜模铸造时合金的中心平均冷却速率为2×103K/s[15],故A, B, C3点的冷却速率远大于铜模铸造冷却速率和Zr65非晶合金的临界冷却速率,所以在激光单点熔凝非晶合金时,热影响区可能会避免晶化。

      Figure 3.  Time-temperature curves at A, B and C shown in Fig. 2a

      由于熔池侧面和下方的散热条件不同,导致熔池侧边和熔池底部的冷却速率产生一定的差异。图 4为熔池下方和熔池侧方的热影响区的温度-时间曲线。

      Figure 4.  Temperature-time curves of side and bottom of molten pool with the same peak temperature of 777K

      选取峰值温度均为777K的两点进行比较。可以看出,熔池侧边的平均冷却速率为1.2×105K/s,而在熔池正下方的冷却速率为4.3×104K/s。熔池侧边的冷却速率明显高过熔池下方的冷却速率,所以熔池侧边的晶化程度相对熔池下方应该要小一些。

    • 图 5为激光辐照2s后的纵截面温度云图。中心温度达到2256.65K,周围温度依次降低,温度曲线为椭圆状。图 6显示了激光辐照2s后的横截面温度云图。激光单道熔凝中熔宽5.34mm,熔深1.98mm。

      Figure 5.  Temperature profile of longitudinal section after 2s laser irradiation

      Figure 6.  Cross section temperature cloud map after 2s laser irradiation

      图 7图 6E, D, FG 4个点的温度变化曲线图。E点、D点和F点位于熔池,G点位于热影响区。E点的平均冷却速率为2.75×102K/s,D点的平均冷却速率为2.65×102K/s,F点的平均冷却速率为2.11×102K/s,G点的平均冷却速率为74K/s。可以看出,无论是熔池还是热影响区,其冷却速率远低于单点熔凝的冷却速率,略高于Zr65非晶合金的临界冷却速率。当温度降到Tm以下时,熔融态的非晶合金转为过冷液相,过冷液相会在出现显著晶化之前产生初生相或利于晶化的结构弛豫[16-18],由于冷却速率较低,非晶合金在过冷液相转变为非晶态的时间较长,产生了较多的初生相和结构弛豫,导致更多晶化相的产生。在峰值温度高于Tg低于Tm的热影响区,升温过程中产生的结构弛豫无法被破坏[19-20],在冷却的过程中,与新的弛豫累加起来,导致更严重的晶化倾向。因此热影响区更易产生晶化。

      Figure 7.  Time-temperature curves at E, D, F and G shown in Fig. 6

    • (1) 激光单点熔凝模拟实验中,熔池的冷却速率为6.3×104K/s,热影响区为1.4×104K/s,都远高于铜模铸造法的冷却速率2×103K/s,也都远高于Zr65非晶合金的临界冷却速率1.5K/s,热影响区在过冷液相区停留时间短,难以生成稳定长大的晶核和结构的弛豫。因此, 理论上激光单点熔凝可以制备出无晶化的高纯度非晶合金。

      (2) 激光单点熔凝非晶合金过程中,熔池下方的冷却速率4.3×104K/s,小于熔池侧边的冷却速率1.2×105K/s,因此熔池侧边相比下方更难发生晶化。

      (3) 激光单道熔凝非晶合金过程中,熔池冷却速率为2.11×102K/s,热影响区冷却速率为74K/s,都小于铜模铸造法冷却速率。热影响区在过冷液相存在时间较长,可能会形成晶核并生长,复杂的热历史又会导致结构弛豫的累积,从而发生晶化。

参考文献 (20)

目录

    /

    返回文章
    返回