2. 鞍山正发表面技术工程股份有限公司 鞍山 114044
2. Anshan Zhengfa Surface Technology Engineering Company Limited by Shares, Anshan 114044, China
激光熔覆是用不同的送料方式在熔覆基体表面放置所需的涂层材料,经激光束辐照使涂层和基体同时熔化,并快速凝固形成冶金结合的熔覆层,改善基体表面的耐磨、耐蚀、耐热、抗氧化性及电气特性。具有稀释率小、组织致密、结合性好、效率高、速度快、绿色环保等特点,在航空航天、汽车、海洋、石油、化工等领域应用广泛[1-2]。
20世纪60到70年代,随着大功率激光器的发展,激光束开始进入工业表面处理领域。2016年,Gao W建立了三维热场模型计算熔池温度、冷却速率和固液界面凝固速率,但未考虑Marangoni流对熔池的影响[3]。2017年,甘政涛利用COMSOL建立了包括温度场、流场的熔覆模型,但未对残余应力进行分析,且建模中未考虑基体和粉材物性参数的温变影响,忽略了粉末对激光能量的衰减作用[4]。2018年,YU TB采用Taguchi-grey关联法对激光熔覆工艺参数进行了多响应优化,并验证了Taguchi-grey关联法的可行性[5]。
激光熔覆是动态物理冶金过程,涉及激光、粉末、工件三者之间的交互作用,存在复杂的多场耦合变化,多场耦合相互作用规律如图1。目前对激光熔覆模拟主要集中在对单一温度场的模拟,对热应力场和残余应力场的瞬变规律研究较少。且建模中均忽略了粉末流与激光束间的相互作用,同时没有考虑基体和熔覆粉的物性参数随温度变化。而激光熔覆过程是热-弹-塑-流多场耦合演变过程,各物理场间相互耦合影响,建模中只考虑单一场或部分场必然会造成误差,同时基体和熔覆粉的物性参数温变影响较大,激光束与粉末相互作用明显,忽略这些影响将造成计算误差。因此,文中建立了碟片激光器激光熔覆过程热-弹-塑-流多场耦合模型,综合考虑了粉末流与激光束之间的相互作用,熔池表面张力、浮力对熔池液态金属流动的影响,温度的自熔性及熔覆带形状的瞬时变化,以相图计算法(Calculation of phase diagram, 简称CALPHAD)得出基体和粉材的温变物性参数,计算得出了熔覆温度场、速度场、凝固行为、应力场分布与演变规律。考虑可移动熔池的对流和扩散,通过温度梯度G和凝固生长速率S预测其形态和凝固组织,运用Axioskop 2电子显微镜和Zeiss-∑IGMA HD扫描电镜对熔覆层微观组织进行观察,验证了分析方法的可行性。
1 试验材料与设备基体为45钢,退火处理,成分如表1,熔覆粉为Fe60,成分如表2。试件如图2,熔覆原理如图3。试验选取TruDisk4002激光器,配以六自由度机器人实现激光熔覆。碟片激光器优势在于将棒状晶改为碟片晶,传统棒状晶通过晶体棒径向热传导实现侧面冷却,棒内温度呈抛物线分布,导致热透镜效应影响光束质量,且这种影响随着抽运功率增大而增大。碟片晶体薄且径厚比大,可及时有效冷却,一维热传导使晶体内温度分布非常均匀,解决了热透镜问题,改善了光束质量、转换效率及功率稳定性,比传统激光器具有优势。试验激光器波长为1030 nm,激光功率80~4000 W,输出光束0.2 mm时,光束质量为8 mm·mard。其转换功率可达30%。
将试件沿垂直光束扫描方向剖开,打磨抛光处理,采用4%的硝酸酒精溶液腐蚀,用Axioskop 2电子显微镜、Zeiss-∑IGMA HD场发射电子显微镜观察熔覆层显微组织形貌。
2 激光熔覆工艺过程理论建模与求解激光熔覆过程理论建模基于以下假设:①熔池内金属流为层流、不可压缩牛顿流体;②激光束能量在光斑内呈高斯分布且功率恒定;③材料为各向同性;④粉末流浓度服从高斯分布,掉落到熔池中粉末立即熔化;⑤熔覆材料服从Von Mises屈服准则。
2.1 激光熔覆过程总体控制方程激光熔覆过程总体控制方程为:
式中,ρ为密度,t为时间,u为熔池内金属流动速度,μ为流体动力粘度,p为压强,T为温度,Cp为比热容,k为导热系数,ξ是位移矢量,F为约束力,S是总位移量。方程(1)为连续质量方程。方程(2)为Navier-Stokes动量方程,K0是一个由多孔形态确定的常数,B是一个避免分母为零的极小参数。方程(3)是能量方程,H是金属熔化潜焓,即为∆H=Lfl,其中fl为液体质量分数,可计算为:
式中,下标s、l分别表示固相和液相。方程(4)是应力方程。
2.2 模型求解初始条件与边界条件为 2.2.1 激光热源碟片激光器的高斯热源与普通激光器相比,有效热影响区面积较大且温度均匀,通过改变x与y的取值范围,可调整有效热影响区的大小,以针对碟片激光器高斯热源模型做出适应性改变。采用改进后的移动高斯热源,热通量为:
式中,P为激光功率,a为激光能量的吸收率,R是激光光斑半径,V为扫描速度,x、y分别是激光光斑中心的瞬时坐标值。hc是热传递系数,σb是Stefan–Boltzmann常量,ε是发射率,T0是初始环境温度。
2.2.2 液体/气体动量方程的边界条件式中,等号右侧两项分别代表毛细力和热毛细力,σ为表面张力,n*是曲面法向参数,κ是曲面曲率。
以同轴送粉为例,考虑液/气界面流动,边界移动速度为:
式中,u是液/气界面流速;Vp代表因熔覆粉材的增加导致液/气界面的移动速度。
式中,mf是质量流量,ηm是粉末流率,ρm是粉末密度,而Rp是粉末流半径,z是在z方向的单位向量。
2.2.3 熔覆粉末与激光束的相互作用激光熔覆中,当高能激光束穿过粉末流时,其部分能量被粉末颗粒吸收、反射和散射,剩余能量穿过粉末流到达熔池表面。同时粉末颗粒因吸收激光能量,在落入熔池前温度已经升高[6-8]。根据Lambert-Beer光透射定理和米氏理论,粉末流对激光束的衰减与激光功率、波长、粉末颗粒特性、粉末流浓度及光束穿越的粉末流长度有关,可表达为:
式中,
粉末颗粒在落入熔池前,因吸收激光束能量的温升由热平衡式得出:
式中,Radd为粉末颗粒半径,ap为粉末颗粒对激光能量的吸收率,s为喷嘴到熔池表面的距离,vf为粉末颗粒的速度,ρm为密度,cpm为比热容,∆T'为粉末颗粒的温升,可计算为:
屈服准则如式(13)所示。式中,σ1、σ2、σ3分别为不同方向主应力,σs为屈服极限。
模型中设置基体为45钢,熔覆粉末为Fe60。基于COMSOLMultiphysics平台搭建熔覆过程热-弹-塑-流多场耦合有限元模型。熔覆工件左右对称,为节省计算时间,建立1/2模型,尺寸20 mm×10 mm×6 mm,采用自由四面体划分网格,打开网格自动加密功能,网格包含315 989个域单元、11 325个边界元和368个边单元,如图4。选用热应力、层流、两相流水平集模块瞬时计算熔覆过程中温度场、流场和热应力场的多场耦合变化规律,采用基于任意拉格朗日-欧拉法(ALE)的动网格描述熔池表面动态形状变化过程。模型参数选取如表3。工作初始温度为环境温度,熔池初速为0 m/s。熔覆过程中,基体对激光束能量吸收率α与激光波长、基体粗糙度有关,且随温度变化,由Hagen-Rubens公式得出吸收率。
Parameters | Values |
Laser power, P / W | 900 |
Beam radius, R / mm | 3 |
Scanning speed, V / (mm·s−1) | 1.8 |
Mass flow rate, mf / (g·min−1) | 6.2 |
Mass flow radius, Rp / mm | 2.5 |
Convective heat transfer coefficient, hc / (W·m−2·K −1K) | 100 |
Laser wavelength, λ0 / nm | 1030 |
Environmental temperature, T0 / K | 298.15 |
式中,re为电阻率,λ0为激光波长。熔池物性参数为:
Wbath是熔池热物性参数,Wmatrix和Wpowder分别是基体和熔覆粉的热物性参数。β是混合分数,取决于二者在熔覆层横截面积占比。用CALPHAD法计算得基体、粉材的温变物性参数,如图5,用插值函数导入COMSOL中。设置求解和边界条件,利用SIMPLE求解器求解。
3.2 模拟计算与结果分析 3.2.1 温度场计算结果熔覆温度对熔覆层的相变凝固行为、晶体生长、形核率及微观组织具有重要影响[9-11]。计算得出从0时刻到700 ms温度场变化规律,选取3个时刻的计算结果分别如图6,0 ms时刻熔覆从x=0 mm处开始,100 ms时熔覆0.18 mm,500 ms时熔覆0.9 mm,700 ms时熔覆1.26 mm。激光熔覆过程信息采集轨迹如图7,绘制熔覆温度瞬时变化曲线,分别如图8。图6(a)~(c)计算结果表明:熔覆初始阶段,因激光辐照产生急热效应,加热区温度迅速升高,到700 ms时刻熔池温度达最高值2600 K,最高温度出现在移动光斑中心偏后位置。当最高温度超过了固相线后,基体发生固液相转变,开始形成熔池,最终形成了近似2 mm×1.5 mm×1 mm的椭球体熔池。不同时刻熔池前方温度等高线分布较为密集,温度变化较为剧烈,后方温度等高线分布较为稀疏,温度变化较为平缓。图6(d)~(i)是将熔池等距剖分观察到的温度分布,结果表明:随着粉末颗粒不断添加到熔池内,导致熔覆层逐渐变厚,但每一时刻熔池最高温度始终处于熔覆层的顶端表面,不同剖面上温度围绕光斑呈近似椭圆形分布,温度值由内向外逐渐递减。图8是沿着图7不同轨迹采集的温度曲线,计算结果表明:不同时刻,沿着熔池深度1号轨迹方向,靠近光斑温度呈递增趋势,增长变化梯度以熔池热影响区边界呈现明显不同阶段变化,接近影响区前温度变化梯度逐渐增大,进入热影响区后,温度变化梯度逐渐减小,最终趋于一致。沿着熔池y向2号轨迹,靠近光斑温度呈递增趋势,增长变化梯度随着时间逐渐减小,最终趋于一致。沿着扫描方向3号轨迹,不同时刻的温度先增大后减小,增长梯度要高于减小梯度,光斑偏后位置温度最高,最高温度随熔覆过程逐渐增大。
3.2.2 速度流场计算结果Marangoni效应是熔池中存在的对流现象。用Marangoni数表征流动特性,取决于几何形状[12-15]。熔池内Marangoni力会影响流体流动和温度分布,如图9,进而改变熔池的生长,导致材料内部产生应力,引起变形裂纹。
计算得出从不同时刻流场变化规律,选取3个时刻的计算结果分别如图10所示,100 ms时熔覆0.18 mm,500 ms时熔覆0.9 mm,700 ms时熔覆1.26 mm。图10计算结果表明:由于表面张力的温度系数为负,熔池金属流从激光束中心向熔池的边缘移动,熔池中心流速小,熔池与基体交界处流速大,最大流速出现在熔池前方。t=700 ms时,熔池内金属流速最大值达到0.23 m/s,熔池被强烈的Marangoni流所控制。绘制熔覆流速瞬时变化曲线,如图11所示。图11(a)表明:沿着1轨迹线,远离熔池区流速为0,越接近光斑,流速呈现先增大后减小的变化,最大流速逐渐变大。图11(b)为沿着2轨迹线熔池内部流速曲线,直到300 ms以前,熔池内流速变化趋势相近,流速最大值不断升高,300 ms时刻最大流速达流的增强,t=700 ms时,熔池内金属流速最大值达到0.19 m/s,呈现微小变动。300 ms以后流速呈现先增大后减小,然后继续增大,最大流速仍随着时间逐渐变大,700 ms时刻最大流速达到0.243 m/s。图11(c)表明:沿着扫描方向,熔池流速呈现两边高、中间低的“双峰”分布,在0.0116 m位置出现最高峰值,流速为0.284 m/s,最高峰位置滞后熔池前端0.0004 m。熔池外仍不存在流速,且不同位置流速最大值随着时间逐渐增大。
3.2.3 应力场计算结果图12为激光熔覆过程中不同时刻的应力状态。激光熔覆中,因高能激光束辐照,基体因局部受热熔化而发生热膨胀,而局部体积热膨胀受限于周围低温区域的约束,产生弹性热应力;同时,急热导致材料屈服极限随温度的升高而下降,当局部热应力值超过材料屈服极限时,熔覆区产生局部弹塑性热压缩变形;熔覆后冷却阶段,熔覆区相对周围区域将缩短、变窄形成塑性热压缩变形[16-18]。当温度降至初始温度后会形成残存于零件中的残余应力。如图12所示,100 ms时熔覆0.18 mm,700 ms时熔覆1.26 mm,1000 ms时熔覆1.8 mm。图12(a)~(f)表明:靠近扫描前端的熔池边界应力较大,熔覆热应力呈先快速增大后趋于稳定的趋势;700 ms以后,等效热应力在548 MPa左右趋于稳定;熔池边界热应力明显高于熔池内部,熔池底部热应力最大。图12(g)~(i)表明:沿1轨迹线,熔池应力逐渐呈现峰值,熔池顶部应力达550 MPa,明显高于底部;沿着2轨迹线,应力呈现先缓慢增大而后递减趋势,最终接近0 MPa,熔池中没有应力。熔池边缘处应力达最大值,为460 MPa;沿3轨迹线,应力呈现明显“驼峰”分布,熔池中应力为0 MPa,熔池边界应力偏高,远离边界应力逐渐递减。
3.3 枝晶生长试验多场耦合作用对熔池的凝固行为有显著影响。图13描述了冷却速率、温度梯度和凝固速度对凝固微观组织结构和尺寸的影响,表明冷却速率G×S对凝固组织尺寸大小的影响,较高的冷却速率导致组织更加精细[19]。G/S决定了微观组织的形态,即从平面晶、胞晶、柱状枝晶、等轴晶的形态变化。影响凝固微观组织的两个关键因素是凝固前的温度梯度G和固/液界面推移速率S。
∆Td是单位距离上温度的变化值,∆d是单位距离,θ为固/液界面推移速率(即凝固生长速率)方向与激光扫描速度方向的夹角。接近熔池与基体相交接的底部,S与V近似垂直,越靠近熔池顶端,θ角越小,S越大。可断定熔池的凝固生长速率从熔池底端到熔池表面呈增大趋势。G由模拟计算得出,如图14。根据模拟得出的G与计算得出S可导出冷却速率G×S和形状控制因子G/S,进而预测熔覆过程中熔池凝固组织的形态和尺寸。
图15为熔覆层试件的试验剖面与A、B、C位置的微观组织,该试件是利用碟片激光器在45钢基体熔覆Fe60所形成的熔覆层,分别运用Axioskop 2电子显微镜、Zeiss-∑IGMA HD型场发射扫描电子显微镜观察得出的微观组织。模拟计算与试验结果对比表明:A位置处于熔池底部,具有较高的G/R和较低的G×R,处于平面结晶区。C位置处于熔池顶部,具有较低G/R和较高G×R,为等轴晶与柱状枝晶共存区。B区处于枝晶区,为柱状的枝晶。熔覆层试验剖面B位置处观察到柱状和粗枝状的枝晶。C位置处为等轴晶和较细的枝晶组织,试验结果与仿真吻合良好,验证了文中所建模型的有效性。
4 结 论(1)采用移动的高斯热源,综合考虑了传导、对流、粉光匹配影响、相变、凝固、应力的相互作用,基于随温度变化的热物理参数,建立了三维热-弹-塑-流多场耦合数值模型,使模拟结果更接近实际。
(2)采用碟片激光器在45钢基体熔覆Fe60粉,加热区温度迅速升高,形成近似2 mm×1.5 mm×1 mm的椭球体熔池,最高温度出现在移动光斑中心偏后位置。熔覆初期,熔池金属流速最大值为0.13 m/s,熔池内热传导起主要作用;700 ms最大金属流速为0.3 m/s,熔池热对流起主要作用。700 ms以后,等效热应力稳定在548 MPa左右。熔池边界热应力明显高于熔池内部应力,熔池底部热应力最大。
(3)对相同工艺和材料下的熔覆过程进行材料学试验分析,模拟结果与试验数据吻合良好,验证了文中方法的有效性。
[1] | DHARANI S. Pradip majumdar finite element analysis of laser irradiated metal heating and melting processes[J]. Optics & Laser Technology, 2010, 42: 855-865. |
点击浏览原文 | |
[2] | LEI Y W, SUN R L, TANG Y, et al. Numerical simulation of temperature distribution and TiC growth kinetics for high power laser clad TiC/NiCrBSiC composite coatings[J]. Optics & Laser Technology, 2012, 44: 1141-1147. |
点击浏览原文 | |
[3] | GAO W, ZHAO S, WANG Y. Numerical simulation of thermal field and Fe based coating doped Ti[J]. International Journal of Heat and Mass Transfer, 2016, 92: 83-90. |
点击浏览原文 | |
[4] | GAN Z T, YU G, HE X L, et al. Numerical simulation of thermal behavior and multicomponent mass transfer in direct laser deposition of Co-base alloy on steel[J]. International Journal of Heat and Mass Transfer, 2017, 104: 28-38. |
点击浏览原文 | |
[5] | YU T B, YANG L, ZHAO Y, et al. Experimental research and multi-response multi-parameter optimization of laser cladding Fe313[J]. Optics and Laser Technology, 2018, 108: 321-332. |
[6] | HUANG Y L, YANG Y Q, WEI G Q, et al. Boundary coupled dual-equation numerical simulation on mass transfer in the process of laser cladding[J]. Chinese Optics Letters, 2008, 6: 356-360. |
[7] | HOFMANA J T, LANGEC D F, PATHIRAJ B, et al. FEM modeling and experimental verification for dilution control in laser cladding[J]. Journal of Materials Processing Technology, 2011, 211: 187-196. |
点击浏览原文 | |
[8] | ALI E, MASOUD A, AMIR K. Correlation between temperature distribution and in situ formed microstructure of Fe-TiC deposited on carbon steel using laser cladding[J]. Applied Surface Science, 2012, 258: 9025-9031. |
点击浏览原文 | |
[9] | BEDENEO D V, KOVALEV O B, SMURVA I, et al. Numerical simulation of transport phenomena, formation the bead and thermal behavior in application to industrial DMD technology[J]. International Journal of Heat and Mass Transfer, 2016, 95: 902-912. |
点击浏览原文 | |
[10] | HAO M Z, SUN Y W. A FEM model for simulating temperature field in coaxial laser cladding of TI6AL4V alloy using an inverse modeling approach[J]. International Journal of Heat and Mass Transfer, 2013, 64: 352-360. |
点击浏览原文 | |
[11] | NIE P, OJO O A, LI Z G. Modeling analysis of laser cladding of a nickel-based superalloy[J]. Surface & Coatings Technology, 2014, 258: 1048-1059. |
点击浏览原文 | |
[12] | WU J Z, WEI H Y, YUAN F B, et al. Effect of beam profile on heat and mass transfer in filler powder laser welding[J]. Journal of Materials Processing Technology, 2018, 258: 47-57. |
[13] | LEE Y S, NORDIN M, BABU S S, et al. Influence of fluid convection on weld pool formation in laser cladding[J]. Welding Research, 2014, 93: 292s-300s. |
点击浏览原文 | |
[14] |
周佳平. 激光沉积制造应力演化及其控制[D]. 沈阳: 沈阳航空航天大学, 2016: 15-20.
ZHOU J P. Research on stress evolution mechanism and control of laser deposition manufacturing[D]. Journal of Shenyang Aerospace University, 2016:15-20(in Chinese). |
[15] |
王福雨, 刘伟军, 赵宇辉. 复杂薄壁零件激光快速成型过程热力耦合场的数值模拟[J]. 机械工程学报, 2013, 49(5): 192-198.
WANG F Y, LIU W J, ZHAO Y H. Thermo-mechanical coupling field simulation of complex thin-wall part laser rapid prototype process[J]. Journal of mechanical engineering, 2013, 49(5): 192-198 (in Chinese). |
[16] |
龙日升, 刘伟军, 卞宏友. 扫描方式对激光金属沉积成形过程热应力的影响[J]. 机械工程学报, 2007, 43(11): 74-81.
LONG R S, LIU W J, BIAN H Y. Effects of scanning methods on thermal stress during laser metal deposition shaping[J]. Chinese journal of mechanical engineering, 2007, 43(11): 74-81 (in Chinese). |
点击浏览原文 | |
[17] | CHEN C R, LIAN G F, JIANG J B, et al. Simplification and experimental investigation of geometrical surface smoothness model for multi-track laser cladding processes[J]. Journal of Manufacturing Processes, 2018, 36: 621-628. |
点击浏览原文 | |
[18] | LIN J B, CHEN C J, ZHANG M, et al. The effects of heat treatment on microstructure and mechanical properties of Tribaloy 400 coatings deposited by laser cladding[J]. Journal of Materials Engineering and Performance, 2018, 27(12): 6339-6348. |
点击浏览原文 | |
[19] |
王丽. 激光熔覆YSZ@Ni涂层的温度场和应力场有限元研究[D]. 南昌: 南昌航空大学, 2017: 31-32.
WANG L. Finite element study of temperature field and stress field of laser cladding YSZ@Ni coating[D]. Nanchang: Journal of Nanchang Hangkong University, 2017: 31-32 (in Chinese). |
点击浏览原文 |