-
$\begin{array}{l} \mathit{\rho c}\frac{{\partial T}}{{\partial t}} - \frac{\partial }{{\partial x}}\left( {{\mathit{\kappa }_x}\frac{{\partial T}}{{\partial x}}} \right) - \frac{\partial }{{\partial \mathit{y}}}\left( {{\kappa _\mathit{y}}\frac{{\partial T}}{{\partial y}}} \right) - \\ \;\;\;\;\;\;\;\;\;\;\;\frac{\partial }{{\partial \mathit{z}}}\left( {{\mathit{\kappa }_z}\frac{{\partial T}}{{\partial \mathit{z}}}} \right) - \mathit{\Phi = }{\rm{0}} \end{array} $
(1) 式中, ρ为材料的密度;c为材料的比热容;κx, κy, κz分别为x, y, z方向上随温度变化的热导系数;Φ为内部热源的热流率; T为瞬时温度值。
-
利用ANSYS软件建立有限元模型。熔覆层采用solid70单元,基体采用solid90单元。该单元有8个节点,每个节点只有一个温度自由度。该单元可用于3维的稳态或瞬态热分析问题,并可以补偿由于恒定速度场质量输运带来的热流损失[14]。图 1为solid70 3维热实体单元的示意图。
为了研究多道熔覆过程温度场变化,建立如图 2所示有限元模型。下层较大面积的网格单元为板状基体,基体尺寸为110mm×60mm×8mm。为提高运算速度,简化模型,忽略零件的圆角等细节。为减少计算量,要尽量减少网格单元数量。在基材下层远离熔覆层的区域采用尺寸较大的网格,单元尺寸为2mm。最上部网格代表熔覆沉积层,熔覆层尺寸30mm×40mm。单道熔覆宽度4mm,搭接率取50%,熔覆14道,熔覆1层。熔覆层网格较细,单元尺寸1mm。熔覆层与基材结合区域采用与熔覆层尺寸相似的单元过渡。因此基材单元尺寸由上至下逐渐变大。
-
目前大部分激光熔覆轨迹中的扫描顺序都按照“单向逐次顺序”、“往复逐次顺序”等简单模式如图 3所示。每一道激光扫描位置的改变将会影响整个基体上的温度场。每产生一道合金层,整体基材上的温度场都会随之变化。所以可以通过规划扫描轨迹,再结合不同材料的凝固成形规律。使每一次熔覆都能在相对最合适(优)的区域进行以获得最大、最好的控形控性效果。本文中称这种方法为评估选择法。
评估选择法需要在每一道激光扫描之前,对基材温度场进行评估,并按照凝固理论来择优选择扫描位置。根据上面的结论,数值模拟方法可以很好地预测激光熔覆过程温度场分布趋势及温度随时间变化情况。所以可以先应用数值模拟方法来得到基体的瞬态温度场,再根据不同材料快速凝固过程中组织成形规律来设计相应的评估方法,通过评估选择激光熔覆位置。从而产生最好的组织形貌及最少的缺陷,以达到优化轨迹的目的。
下面对评估选择方法做详细说明。由于温度是影响凝固过程的最重要因素,所以首先将温度值做为评估的特征参量。如图 4a所示,在正方形基体上进行单向多道熔覆,一共需要熔覆n道。每一道的位置已确定,需要确定的是熔覆扫描顺序。
首先可以选择最边缘位置进行第1次熔覆,见图 4b中阴影位置。应用上面数值模拟方法,在第1次熔覆结束后,求解整体温度场瞬态解。以剩余未熔覆的区域为研究对象(2, 3, 4, …,n道), 取出每道区域中的有限元模型单元的温度瞬态值。比如,第2道区域,包含j个单元,单元编号为2-1, 2-2, …, 2-j,如图 4c所示。此刻每个单元的温度值瞬态解为t2-1, t2-1, …, t2-j。第2道区域中所有单元温度值瞬态解的平均值为:
$\overline {{t_2}} = \left( {{\mathit{t}_{2 - 1}} + {t_{2 - 2}} + {\mathit{t}_{2 - 3}} + \cdots + {t_{2 - \mathit{j}}}} \right)/j $
(2) 式中,j为2道区域单元总数。同理第n道区域的单元温度值瞬态解的平均值为:
$\overline {{t_n}} = \left( {{\mathit{t}_{n - 1}} + {t_{n - 2}} + {\mathit{t}_{n - 3}} + \cdots + {t_{n - \mathit{j}}}} \right)/j $
(3) 式中,j为n道区域单元总数。根据以上方法得到了各个待熔覆区域的温度特征参量$\overline {{t_\mathit{n}}} $。
同理,温度在基材区域上分布的均匀程度也可以作为评估特征参量。设任意第n道区域的单元温度值瞬态解的标准差值为:
$\mathit{\alpha = }\sqrt {\frac{{\sum\limits_{i = 1}^j {\left( {\overline {{t_n}} - {t_{n - i}}} \right)} }}{{j - 1}}} $
(4) 式中,标准差α反映了温度的分布均匀性,所以将α也定为评估特征参量。
根据凝固原理,较大的温度梯度可以使凝固组织形成细小的柱状晶和枝晶。同时为了缓解组织开裂的趋势,应尽量使基材温度分布均匀。如果以此为评估准则,那么评估参量温度瞬态解均值$\overline {{t_\mathit{n}}} $和标准差α必须为最小。根据这个准则,图 4中所示的例子要将最后一道区域作为第2次熔覆区域, 如图 4d所示。需要注意的是, 根据不同的控形控性要求和不同材料的凝固规律所选用的评估参量及评估准则并不相同, 必须以工艺试验结果作为依据。
-
激光熔覆再制造系统[15-16]包含激光制造系统、送粉系统、运动执行机构及工作台、软件控制系统以及其它辅助装置,本实验中采用的是美国IPG公司生产的4kW光纤激光器和德国KUKA机器人KR30HA,能实现柔性加工,系统布局如图 5a所示。在所选位置采用K型热电偶对温度变化进行测量。
-
为了研究几种方法对温度场的影响,将单向逐次、往复逐次、评估选择几种方法所形成的温度场关键特征参量进行了对比。根据熔覆过程的温度场瞬态解来研究每一次熔覆前目标区域的基材温度。图 6、图 7表示的是第2次~第14次熔覆前,3种方法所面对的熔覆目标区域(道)的基材平均温度和温度标准差。比如图 6和图 7中,对于单向逐次法,第3次熔覆的目标就是第3道。熔覆前第3道的温度均值为573℃,温度标准差为150℃左右。图 6中上面的箭头引出的是单向逐次法熔覆第3次前的温度场云图。可以看出,第3道(云图中白线)紧靠上一次完成的第2道,所以基础温度较高,温度值分布也不均匀,温度标准差大。而评估选择法第3次的熔覆目标区域是第2道(云图中白线),远离上一次熔覆的第14道(图 6下面的箭头所引出云图)。熔覆前第2道的温度均值为190℃,温度标准差为16℃左右。都远低于单向逐次法。从图 6中可以看出,在前7次,评估选择法每一次所选区域的基材(础)平均温度远低于单向逐次法和往复逐次法。之后,几种方法的基础温度开始逐步接近。第14次熔覆,评估法的基础平均温度大于单向逐次法和往复逐道法。这是由于评估法的原理决定的。
图 7中的基础温度标准差表示的是熔覆区域基础温度分布的均衡状况。从图上可以看出,评估选择法在每一次熔覆前目标区域的温度标准差都是最低的,特别是在熔覆前期比其它两种方法均匀很多。说明评估选择法每一次熔覆都是在一个温度场分布非常均匀的区域展开。
-
为了研究不同激光扫描顺序对激光熔覆沉积层组织的影响,应用上面的单向逐次、双向逐次和评估选择法对图 5b所示的45#钢试验板进行多道激光熔覆。基材成分见表 1。采用铁基粉末Fe50(HRC50~HRC58),粉末成分见表 2。激光功率1600W, 光斑直径4mm,送粉速率30g/min,载气速率600L/h,激光扫描速率5mm/s,搭接率50%。熔覆单层,14道。
Table 1. Composition(mass fraction w) of 45# steel
Cr S Si C P Mn Ni Cu ≤0.003 ≤0.001 0.002~0.004 0.004~0.005 ≤0.001 0.005~0.008 ≤0.003 ≤0.003 Table 2. Composition(mass fraction w) of Fe50-based powder
Cr B Si C Ni Mo Fe 0.150~0.180 0.030~0.040 0.035~0.045 0.006~0.010 0.280~0.300 0.040~0.060 balance 图 8所示为试验熔覆层的宏观形貌,单向逐次法形成的样件形貌最平整。而评估选择法由于最后两次熔覆都在中部区域(图中黑圈处),所以中部凸起比较高。双向逐次法外部形貌最差,大部分区域形成缺陷,这是由于试验采用了侧向送粉装置。因为侧向送粉装置对激光扫描方向有严格要求,而双向逐次法中必然有一个方向与要求相反,所以形成缺陷。随着逐步搭接,这个区域的缺陷趋势放大,最终形成如图所示的形貌。
将单向逐次法和评估选择法得到的试验样件沿中部截面切开, 如图 8所示,并在相应位置取样进行显微硬度测试。图 8a中的S1, S2, S3处为单向逐次法试样显微硬度测试位置;图 8c中的X1, X2, X3处为评估选择法试样显微硬度测试位置;每个位置测试3次求均值。测试结果见表 3。
Table 3. Micro-hardness of one way successive method and evaluation selection method
number S1 S2 S3 X1 X2 X3 microhardness/HV 625.38 620.58 623.34 680.09 673.58 683.01 从表中数据可以看出, 评估选择法试样的显微硬度高于单向逐次法。这是因为对于铁基材料熔池的冷却速度对于相变强化起到重要作用。较快的冷却速度使熔覆层的硬度提升。将S3和X3处切下的试样做电镜扫描,电镜照片如图 9所示。图中显示了熔覆层与基体结合处、熔覆层中部、熔覆层顶部的组织。从组织形态看单向逐次法形成了柱状枝晶,侧向分枝较为发达。生长方向基本垂直于结合界面,呈现典型的外延生长特征。评估选择法择形成了较多的胞状晶。这是由于评估选择法的熔池周围温度梯度更大,影响到熔池内的温度梯度,控制凝固过程更加稳定。从尺度上看, 枝晶的间距较大。
-
为了对比3种方法对基材变形的影响,在试样熔覆面的背后画上网格,测量网格节点的变形情况,试样网格测量点分布如图 10a所示。由于平板试样使用前经过精磨,可以近似假定激光熔覆前整个试样的各点变形量均为零。以基材未熔覆面矩形角点的平均值为参考零点,即其变形为0.000mm。这样,其它测量点所测数值与参考零点之差即可认为是激光熔覆所引起的弯曲变形。
如图 10b所示,将激光熔覆后的试样放在专用的夹具上,样件的侧面可以作为水平基准与工作台接触。在三坐标测量机(其最小分辨率为1μm)上,测量激光熔覆后试样背面网格点的坐标值。每块样件共计测量77个点。为对比3种方法所熔覆样件变形的大小,将所测数据中在厚度方向最大的10个值列出,再将每列数据求均值, 如表 4所示。从数据中可以看出评估选择法的基材在厚度方向上变形最小,往复逐次法最大,单向逐次法在两者之间。三者之间的差距见图 11。
Table 4. Deformation data of thickness direction
measurement point number one way successive/mm evaluation selection/mm reciprocating successive/mm 1 0.840 0.593 1.328 2 0.866 0.596 1.363 3 0.953 0.619 1.385 4 0.955 0.628 1.398 5 0.956 0.638 1.401 6 0.999 0.659 1.410 7 1.034 0.672 1.432 8 1.036 0.675 1.553 9 1.039 0.681 1.554 10 1.044 0.697 1.607 average 0.9722 0.6458 1.4431(max) 为了直观地观察每种方法所产生的变形大小,把所测数据拟合成曲面。图 12是将3种情况放在一起对比。从图上可以看出, 3种方法形成的试验件的两端都沿基材中轴线发生了弯曲; 板状基材向着激光束的方向形成一个弯曲角;评估选择变形量小于往复逐次及单向逐次的变形量。这是由于逐次法,特别是往复逐次法能量输入的位置集中,在短时间内形成热积累效应。
基于温度场评估的激光熔覆顺序决策方法研究
Research of decision method of laser cladding sequence selection based on temperature field evaluation
-
摘要: 为了在激光熔覆再制造过程中得到更优的激光熔覆顺序选择决策方法,采用有限元法对平面基材多道激光熔覆传热学模型的温度场瞬态解进行了理论分析,并利用基于热电偶的测温系统验证了整个数值模拟过程的可靠性。提出了一种评估选择法,即利用数值模拟来分析和评估基体瞬态温度场,根据评估准则选择熔覆过程激光扫描顺序的轨迹优化方法。结果表明,取得单向逐次和评估选择样件的实验硬度数据分别为625.38HV,620.58HV,623.34HV,680.09HV,673.58HV和683.01HV,变形均值为0.9722mm和0.6458mm;评估选择法有最均匀的温度场,熔池周围有最大的温度梯度,其能产生较大的熔覆层硬度及较小的组织尺度,同时测量变形较小。该方法为激光熔覆的顺序选择提供了重要的参考价值。Abstract: In order to obtain a better decision method of laser cladding sequence selection in laser cladding remanufacturing process, finite element method was used to analyze transient solution of temperature field of multichannel laser cladding heat transfer model of flat substrate. The reliability of whole numerical simulation process was verified by temperature measurement system based on thermocouple. An evaluation selection method was proposed. Numerical simulation was used to analyze and evaluate the transient temperature field of the matrix, and select laser scanning sequence during the cladding process according to evaluation criteria. It was the trajectory optimization method. The results show that experimental hardness data are 625.38HV, 620.58HV, 623.34HV, 680.09HV, 673.58HV and 683.01HV of one way successive method and evaluation selection method respectively. The average deformation values are 0.9722mm and 0.6458mm respectively. Evaluation selection method has the most uniform temperature field and the maximum temperature gradient around the pool, which can produce larger cladding layer hardness, smaller microstructure scale and smaller measurement deformation. This method provides an important reference value for sequential selection of laser cladding.
-
Table 1. Composition(mass fraction w) of 45# steel
Cr S Si C P Mn Ni Cu ≤0.003 ≤0.001 0.002~0.004 0.004~0.005 ≤0.001 0.005~0.008 ≤0.003 ≤0.003 Table 2. Composition(mass fraction w) of Fe50-based powder
Cr B Si C Ni Mo Fe 0.150~0.180 0.030~0.040 0.035~0.045 0.006~0.010 0.280~0.300 0.040~0.060 balance Table 3. Micro-hardness of one way successive method and evaluation selection method
number S1 S2 S3 X1 X2 X3 microhardness/HV 625.38 620.58 623.34 680.09 673.58 683.01 Table 4. Deformation data of thickness direction
measurement point number one way successive/mm evaluation selection/mm reciprocating successive/mm 1 0.840 0.593 1.328 2 0.866 0.596 1.363 3 0.953 0.619 1.385 4 0.955 0.628 1.398 5 0.956 0.638 1.401 6 0.999 0.659 1.410 7 1.034 0.672 1.432 8 1.036 0.675 1.553 9 1.039 0.681 1.554 10 1.044 0.697 1.607 average 0.9722 0.6458 1.4431(max) -
[1] XU B Sh. Development status and prospect of green remanufacturing engineering[J]. Engineering Sciences, 2011, 13(1):4-10(in Chinese). [2] ZHANG A F, LI D Ch, LIANG Sh D, et al. Development of laser additive manufacturing of high-performance metal parts[J]. Aeronautical Manufacturing Technology, 2016, 517(22):16-22(in Chinese). [3] LUO K L, GUO Sh Q, HE Y, et al. Repairing TC4 titanium stator blade of aero-engine fan casing by laser cladding[J]. China Surface Engineering, 2015, 28(6):141-146(in Chinese). [4] HAN Y Y, LU J J, LI J F, et al. Lathe spindle remanufacturing based on laser cladding technology[J]. China Surface Engineering, 2015, 28(6):147-153(in Chinese). [5] SON S, KIM S, LEE K H. Path planning of multi-patched freeform surfaces for laser scanning[J]. The International Journal of Advanced Manufacturing Technology, 2003, 22(5):424-435. [6] DENG Sh Sh, YANG Y Q, LI Y, et al. Planning of area-partition scanning path and its effect on resi-dual stress of SLM molding parts[J]. Chinese Journal of Lasers, 2016, 43(12):1202003(in Chinese). doi: 10.3788/CJL [7] BIAN H Y, FAN Q Ch, LI Y, et al. Scan path generating method based on temperature subarea of laser deposition shaping[J]. Journal of Mechanical Engineering, 2015, 51(24):57-62(in Chinese). doi: 10.3901/JME.2015.24.057 [8] XU H Y, LI T, LI H B, et al. Study on quality prediction and path selection of 316L laser cladding[J].Laser Technology, 2018, 42(1):53-59(in Chinese). [9] LI Sh Y, TIAN X G, LI Ch B. Study on wear-resistance and scanning path of laser alloying on end cap working face[J]. Chinese Journal of Lasers, 2013, 40(2):0203004(in Chinese). doi: 10.3788/CJL [10] LONG R Sh, LIU W J, BIAN H Y, et al. Effects of scanning methods on thermal stress during laser metal deposition shaping[J]. Journal of Mechanical Engineering, 2007, 43(11):74-81(in Chinese). doi: 10.3901/JME.2007.11.074 [11] CHEN L, XIE P L. Theory and experimental research on controlling crack in double-scanning laser cladding process[J]. Transactions of the China Welding Institution, 2011, 32(2):65-68(in Chinese). [12] QUAN X M, DING L, WEI X. Analysis of temperature field of laser cladding Ni-based alloy[J]. Laser Technology, 2013, 37(4):547-550(in Chinese). [13] ZHAO H Y, SHU F Y, ZHANG H T, et al. Numerical simulation on temperature field of laser cladding based on birth-death element method[J]. Transactions of the China Welding Institution, 2010, 31(5):81-84(in Chinese). [14] LIU H, YU G, HE X L, et al. Three-dimensional numerical simulation of transient temperature field and coating geometry in powder feeding laser cladding[J]. Chinese Journal of Lasers, 2013, 40(12):1203007(in Chinese). doi: 10.3788/CJL [15] WANG X L, SUN W L, ZHANG J J, et al. Study on damage boundary extraction and shape reduction of laser cladding parts[J]. Laser Technology, 2017, 41(5):675-679(in Chinese). [16] HUANG Y, SUN W L, CHEN Y. Research on laser remanufacturing method of complex shaft parts based on NURBS interpolation[J]. Infrared and Laser Engineering, 2017, 46(5):45-51(in Chinese).