上海理工大学学报  2020, Vol. 42 Issue (5): 441-447   PDF    
超快激光切割金属材料中快速相变的数值模拟
张琦, 李凌     
上海理工大学 能源与动力工程学院,上海 200093
摘要: 采用二维双温度模型对多脉冲激光照射金薄板的相变传热情况进行研究,通过等效比热容方法确定固液界面的位置并研究了激光参数对传热过程的影响。结果表明,当激光垂直照射金薄板时,表面整体温度不断上升,而在脉冲间隔时间内略有下降,但温度的峰值相比激光作用的中心有延迟。随着激光的移动,激光照射点处的温度会出现一个峰值,之后会回落。激光作用时间内热影响区的横向变化比纵向大,移动光源向内部的热传递占主要作用,脉冲间隔时间内,热影响区的纵向变化比横向大,金薄板内部导热占主要作用。随着激光的照射,熔化状态横向不断增加的同时,熔化的深度也在不断增加。提高入射脉冲激光能量,会导致熔化的时间提前且熔化深度增加。
关键词: 多脉冲激光     二维双温度模型     固液界面     等效比热容     热影响    
Numerical simulation on the rapid phase transformation of metal under the action of moving light source
ZHANG Qi, LI Ling     
School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
Abstract: The phase change heat transfer of multi-pulse laser irradiated gold thin plates was studied by using a two-dimensional dual temperature model. The position of solid-liquid interface was determined by the equivalent specific heat capacity method and the influence of laser parameters on the heat transfer process was studied. The results show that when the laser illuminates the gold thin plate vertically, the overall temperature of the surface increases continuously, but decreases slightly during the pulse interval, and the temperature peak is slightly delayed compared to that at the center of the laser action. As the laser moves, the temperature at the laser spot will show a peak and then fall back. The lateral change of the heat-affected zone during laser action is larger than the longitudinal one, where the heat transfer from the moving light source to the inside plays a major role. During the pulse interval, the longitudinal change of the heat-affected zone is larger than the lateral one, where the internal heat conduction of the gold thin plate plays a major role. As the laser is irradiating, the melting state extends continuously in the lateral direction, while the depth of melting is also constantly increasing. Increasing the incident pulse laser energy causes the melting time to advance and the melting depth to increase.
Key words: multi-pulse laser     two-dimensional two-temperature model     solid-liquid interface     equivalent specific heat capacity     thermal influence    

激光的微加工技术是利用激光束与物质相互作用的特性对材料(包括金属与非金属)进行切割、焊接、表面处理及打孔等的一门加工技术。超短脉冲激光切割材料时显现出的“冷加工[1]”特性,使其成为近年来激光微加工领域的研究热点。利用超短脉冲激光峰值功率密度极高、热效应极小的特点[2],可以有效地避免热效应所导致的各种影响。在激光与材料相互作用过程中,会出现熔化以及凝固[3]等快速相变传热现象,与固体激光器和CO2激光器的“热”激光切割相比,超短激光的加工精度更高。目前对于超短激光和金属材料的研究很多[4-6],但大多数是单脉冲且光源不动,Zhou等[7]研究了超短激光与纳米材料的微尺度的传热机理,金方圆等[8]采用有限元的方法研究了铝膜内电子和晶格温度的分布变化。他们的研究结果表明,超短脉冲激光与金属相互作用过程中激光的能量通过光子沉积在金属材料的表面,而金属材料表面中的电子首先通过吸收光子的能量,并通过电子间的碰撞使电子系统在大约几百fs的时间内得到加热,而此时金属晶格温度并没有太大变化[9]。随后,能量通过电子与离子的碰撞传递给晶格,使晶格系统得到加热。实际脉冲激光切割过程中,激光是以一定的速度移动的,激光能量的变化情况、脉冲激光的移动速度以及脉冲之间的间隔时间都对加工过程中材料的相变传热有很大影响,并进而影响成型件的质量。本文将双温度模型[10]与等效比热容方法[11]相结合,建立二维模型对多脉冲激光照射下的金薄板相变传热过程进行模拟研究,并探讨激光参数对该过程的影响,研究有助于推广激光微加工技术在工业中的应用。

1 物理模型与数学方法 1.1 物理模型

图1为物理模型,一束超短脉冲激光沿上表面从左侧对金薄板进行照射,并以速度v向右移动。为了简化计算,建立二维模型,如图1所示,金的长度为L,高度为H,初始温度为Ti,其他表面均按绝热条件处理。


图 1 物理模型 激光照射表面 Fig. 1 Physical model. Laser irradiated surface
1.2 数学方法

当超短脉冲激光照射到金薄板表面时,电子和晶格之间处于不平衡传热,采用双温度模型(TTM)对上述非平衡过程进行描述,直角坐标系下的非线性方程为

$\begin{split} {{{C}}_{\rm{e}}}\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{t}}}} =& \frac{\partial }{{\partial {{x}}}}\left( {{{{k}}_{\rm{e}}}\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{x}}}}} \right) + \frac{\partial }{{\partial {{y}}}}\left( {{{{k}}_{\rm{e}}}\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{y}}}}} \right) -\\ & {{G}}\left( {{{{T}}_{\rm{e}}} - {{{T}}_{\rm{l}}}} \right) + {{S}} \end{split}$ (1)
$\begin{split} {{{C}}_{\rm{l}}}\frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{t}}}} =& \frac{\partial }{{\partial {{x}}}}\left( {{{{k}}_{\rm{l}}}\frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{x}}}}} \right) + \frac{\partial }{{\partial {{y}}}}\left( {{{{k}}_{\rm{l}}}\frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{y}}}}} \right) +\\ & {{G}}\left( {{{{T}}_{\rm{e}}} - {{{T}}_{\rm{l}}}} \right) + {\rm{\rho }}{{{L}}_{\rm{f}}}\frac{{{\rm{d}}f}}{{{\rm{d}}t}} \end{split}$ (2)

式中:C为体积热容;T为温度;t为时间;k为导热系数;G为电子−声子耦合因子;S为热源,由于激光被用作边界条件,在本研究中S=0;下标e和l代表电子和晶格;ρ为密度;Lf为熔化潜热;f为液体分数。

文献[12]可以获得金的物性参数。假设液体分数f随温度线性变化。

$ {{f}} = \left\{ {\begin{array}{*{20}{l}} 1,&{{{{T}}_{\rm{l}}} > \left( {{{{T}}_{\rm{m}}} + {\rm{\delta }}{{{T}}_{\rm{l}}}} \right)}\\ {\dfrac{{{\rm{}}{{{T}}_{\rm{l}}} - {{{T}}_{\rm{m}}} + {\rm{\delta }}{{{T}}_{\rm{l}}}}}{{2\,{\rm{\delta }}\,{{{T}}_{\rm{l}}}}}},&{{{T}}_{\rm{m}}} - \delta {{{T}}_{\rm{l}}} < {{{T}}_{\rm{l}}} < {{{T}}_{\rm{m}}} + \delta {{{T}}_{\rm{l}}}\\ 0,&{{{{T}}_{\rm{l}}}<{{{T}}_{\rm{m}}} - \delta {{{T}}_{\rm{l}}}}\!\!\!\!\!\! \end{array}} \right.\!\!\!\!\!\! $ (3)

Tm是熔点温度的情况下,δ是一个小数,以防止除以零。

图2是超快脉冲激光在时间及空间上的分布示意图,其函数表达式为


图 2 激光能量强度在时间和空间上的分布示意图 Fig. 2 Schematic diagram of laser intensity distribution in time and space
$\begin{split} {{q}} = {{A}}\frac{{2{{E}}}} {{{{\text{π} }}{{{R}}^2}{{{t}}_{\rm{p}}}}}&{\rm{exp}}\left( { - 2\frac{{{{{r}}^2}}}{{{{{R}}^2}}}} \right) \left\{ \exp \left[ - 2.77\frac{\left( {{t}} - 2{{{t}}_{\rm{p}}} \right)^2}{{{t}}_{\rm{p}}^2} \right] \right. + \cdots\\ & \left. \exp \left[ - 2.77\frac{\left( {{t}} - {{n}}{{t}}_{\rm{p}} \right)^2}{{t_{\rm{p}}^2}} \right] \right\} \\[-15pt] \end{split}$ (4)
$ {{r}} = \sqrt {{{\left( {{{x}} - {{{x}}_0}} \right)}^2} + {{\left( {{{y}} - {{{y}}_0}} \right)}^2}} x = vt $

式中:x0y0为激光中心点;A为金对激光的吸收率;E为单脉冲能量;tp为脉冲宽度;R为激光源的半径;n的值与脉冲宽度和脉冲间隔以及脉冲数有关。

t=0开始计算,初始条件:

$ {{{T}}_{\rm{e}}}\left( {{{x}},\ {{y}},\ 0} \right) = {{{T}}_{\rm{l}}}\left( {{{x}},\ {{y}},\ 0} \right) = 300\;{\rm{K}} $ (5)
$ \frac{{\partial {T_{\rm e}}}}{{\partial t}} = \frac{{\partial {T_{\rm l}}}}{{\partial t}} = 0 $ (6)

金薄板上表面收到激光的辐照,忽略表面热辐射热损失,左、右、下设置为绝热条件:

$ {{y}} = {{H}},\qquad {{q}} = - {{k}}\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{y}}}} $ (7)
$ {{x}} = 0,\quad{{x}} = {{L}},\quad\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{x}}}} = \frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{x}}}} = \frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{y}}}} = \frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{y}}}} = 0 $ (8)
$ {{y}} = 0,\qquad{\rm{}}\frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{x}}}} = \frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{x}}}} = \frac{{\partial {{{T}}_{\rm{e}}}}}{{\partial {{y}}}} = \frac{{\partial {{{T}}_{\rm{l}}}}}{{\partial {{y}}}} = 0 $ (9)
2 结果与讨论 2.1 模型验证

为了对模型进行验证,首先对单脉冲激光照射金膜的传热过程利用软件comsol5.3进行模拟。金膜的初始温度为300 K,tp=600 fs,模拟从−2tp开始,计算得到了金膜右侧刚发生熔化时的能量密度,将其损失阈值模拟结果与实验数据[12]进行了对比,结果如图3所示。从图3中可以看出,模型与实验数据吻合良好。


图 3 损伤阈值的模型和实验数据比较 Fig. 3 Model and experimental data comparison of the damage threshold
2.2 激光移动作用下金薄板的熔化过程

E=0.3 μJ,tp=10 ps,R=10 μm,脉宽间隔dt=40 ps的移动激光照射金薄板的过程进行了模拟研究。金薄板L=25 μm,H=0.5 μm,采用自适应网格技术生成网格,初始网格数为100×100,初始温度为300 K。图4给出了不同时刻其温度场的变化情况。从图4中可以看出,在激光移动过程中,表面温度T整体在不断地上升,而在2个脉冲间隔内即t=40~80 ps以及t=120~160 ps时,表面的最高温度有所降低,这是因为热量不断地向低温区传递。由图4中也可以发现,在脉冲间隔内,热影响区纵向的变化要比横向的变化大,而在激光照射期间,热影响区横向的变化要比纵向的变化大。这是因为脉冲间隔内没有能量的进入,金薄板内部导热占主要作用,而激光照射时则移动光源向内部的热传递占主要作用。


图 4 不同时刻金薄板晶格的温度分布 Fig. 4 Lattice temperature at different time

金薄板的相变熔化情况如图5所示。其中,红色部分是发生熔化区域。计算发现,在t=22 ps时开始发生熔化,此时熔化最深为0.007 μm。而在2个脉冲间隔内,即t=40~80 ps以及t=120~160 ps时,熔化的部分开始缓慢地再凝固,熔化深度也由0.075 μm到0.07 μm以及0.105 μm到0.102 μm,这是由于热量向低温区域传递,温度有所降低。在激光作用时间内,随着激光的移动,熔化状态横向不断增加的同时,熔化的深度也在不断增加,在脉冲激光结束时,t=200 ps,此时熔化深度为0.11 μm。


图 5 不同时刻金薄板熔化情况 Fig. 5 Melting at different time

图6是金薄板表面点(5,0.5),(12.5,0.5),(20,0.5)处电子和晶格温度随时间的变化曲线。从图6中可以看出,随着激光的移动,不同点的温度都会出现3个峰值,之后会回落。这是因为激光包括3个脉冲,同一个点会受到3个脉冲的叠加影响,激光与金属作用过程中,激光能量在金属表面沉积,电子通过吸收能量从而转换为激发态,又由于电子温度高于附近的晶格温度,因此,电子的温度先是快速升高而后又迅速降低。而晶格接受电子传递的能量后温度会不断上升,之后由于向低温区传热温度逐渐降低。


图 6 表面不同位置的温度 Fig. 6 Surface temperature at different position
2.3 激光能量对照射过程中金薄板快速相变的影响

研究激光能量对照射过程中快速相变的影响,图7为当激光脉宽tp=10 ps,激光半径R=10 μm,脉宽间隔为40 ps时,金薄板在不同能量E的激光照射下界面的温度及切割深度。从图7(a)可以看出,随着入射激光能量的增大,界面的温度峰值越来越高。当3个脉冲作用结束后,温度的最高值均在右侧,但温度的峰值相比激光作用的中心略有延迟,这是由于激光在照射过程中上一个脉冲对下一时刻激光中心位置热影响的累积,使得界面温度的峰值位置在激光中心点之后;在实际工程中,切面最高点与最低点之间的高差垂直距离越小,说明切割效果越好。从图7(b)可以看出,随着入射激光能量的增大,切割的深度也越来越大,但相比较而言,E=0.3 μJ时边缘处切割相对整齐,切割的效果相对较好。


图 7 不同入射能量的温度 Fig. 7 Temperature with different incident energy

不同能量E照射下金薄板相变熔化情况如图8所示。随着入射激光能量的增大,熔化的深度会增加且靶材的熔化体积越来越大,熔化开始的时间也从25 ps提前到19 ps。这是因为,当其他因素都不变的情况下,提高激光能量,相同时间内表面接受能量多,使得温度上升速度变快,温度的峰值也提高。当激光作用结束后,激光能量越大,靶材熔化体积到达最大值的时间越缓慢。在整个过程中靶材的熔化体积会有3个“台阶”,这是因为3个脉冲激光是间歇作用。在E=0.2,0.3 μJ时,脉冲间隔时间内其熔化体积会略有减小,这是因为在这个时间内没有能量注入,无法使得熔化区域的温度持续高于熔化温度,而E=0.4 μJ时,由于能量过大,在最后一个脉冲结束后内部导热进行的时间比较长,但随着时间的推移,会发生再凝固现象。


图 8 不同入射能量的熔化体积 Fig. 8 Melting volume with different incident energy
3 结 论

当移动激光垂直照射金薄板时,在激光作用时间内随着激光的移动热影响区横向的变化要比纵向的变化大,熔化状态横向不断增加的同时熔化的深度也在不断增加。在脉冲间隔时间内,表面的最高温度都有所降低,随着热量不断地向低温区传递,热影响区纵向的变化要比横向的变化大,熔化的部分会有缓慢的再凝固过程。随着入射激光能量的增加,会导致熔化的时间提前且熔化的深度增加。

参考文献
[1]
SUGIOKA K, CHENG Y. Ultrafast laser processing: from micro-to nanoscale[M]. Singapore: Pan Stanford Publishing, 2013.
[2]
邵俊峰, 郭劲, 王挺峰. 飞秒双脉冲激光照射金属薄膜的热行为[J]. 强激光与粒子束, 2014, 26(9): 92-97.
[3]
王志坚, 王宗园, 宋鸿武, 等. TC4钛合金激光熔覆熔池凝固传热研究[J]. 机械设计与制造, 2018(10): 85-88. DOI:10.3969/j.issn.1001-3997.2018.10.023
[4]
QIU T Q, TIEN C L. Heat transfer mechanisms during short-pulse laser heating of metals[J]. Journal of Heat Transfer, 1993, 115(4): 835-841. DOI:10.1115/1.2911377
[5]
SHIOMI M, YOSHIDOME A, ABE F, et al. Finite element analysis of melting and solidifying processes in laser rapid prototyping of metallic powders[J]. International Journal of Machine Tools and Manufacture, 1999, 39(2): 237-252. DOI:10.1016/S0890-6955(98)00036-4
[6]
AN W, ZHU T. Finite element method with diffusion synthetic acceleration for short pulsed laser in turbid media[J]. Chinese Journal of Computational Physics, 2012, 29(1): 121-126.
[7]
ZHOU L, LI L, ZHANG Y W, et al. Source distribution based on Mie’s scattering theory for heat conduction in nanoparticle subject to ultrashort laser irradiation[J]. Numerical Heat Transfer, Part A: Applications, 2014, 66(6): 605-621. DOI:10.1080/10407782.2014.894394
[8]
金方圆, 陈波, 鄂书林, 等. 皮秒激光烧蚀Al膜的理论分析与实验验证[J]. 发光学报, 2014, 35(6): 754-760.
[9]
RETHFELD B, KAISER A, VICANEK M, et al. Femtosecond laser-induced heating of electron gas in aluminium[J]. Applied Physics A, 1999, 69(7): S109-S112. DOI:10.1007/s003390051363
[10]
ANISIMOV S I, KAPELIOVICH B L, PERELMAN T L. Electron emission from metal surfaces exposed to ultra-short laser pulses[J]. Soviet Physics, JETP, 1974, 39(2): 375-377.
[11]
周业涛, 关振群, 顾元宪. 求解相变传热问题的等效热容法[J]. 化工学报, 2004, 55(9): 1428-1433. DOI:10.3321/j.issn:0438-1157.2004.09.008
[12]
STUART B C, FEIT M D, HERMAN S, et al. Optical ablation by high-power short-pulse lasers[J]. Journal of the Optical Society of America B, 1996, 13(2): 459-468. DOI:10.1364/JOSAB.13.000459