1、 高功率激光熔覆 TiCNiCrBSiC 复合涂层的温度分布和 TiC 的生长动力学的数值模拟摘要建立了 Ti6Al4V 合金高功率激光熔覆 TiCNiCrBSiC 复合涂层的三维模型。从模型中可以得到熔池和基底热影响区(HAZ)的温度分布，不同节点的温度曲线，三维形状和熔化区 TiC 的尺寸。基于熔池模拟温度曲线的动力学理论分析，解释了 TiC 颗粒的成核，生长速度，和尺寸，清晰的呈现了熔覆涂层在激光熔覆过程中的相变和微观组织演化。用 TiC 和 NiCrBSiC 的粉末混合物为熔覆材料，矩形Ti6Al4V 合金为基体，以优化的工艺参数制备了稀释率低，冶金结合良好的高质量 TiCNiCrBS
2、iC 复合涂层，通过理论结果和涂层微观组织的比较分析，验证所建立模型的可靠性。结果表明，理论结果与实验例证的一致性良好。引言钛基合金如今被广泛的用作结构材料，如航天器，化学和海洋相关工业，这是由于他们具有极好的强度-重量比和出色的抗腐蚀能力。然而，较弱的滑动摩擦抗力限制了他们更宽的应用。经过激光熔覆陶瓷强化的金属基体母体覆层常常铸造在钛合金表面以提高他们的耐磨性，该覆层使具有非常高的强度和延展性金属和具有独特硬度和耐磨性的陶瓷相结合。为了生产出具有冶金结合的高质量，低的稀释率和最小量的空洞和裂纹的覆层，熔覆过程应该被精确控制。然而，参数（材料特征，激光能量，扫描速度，被取代层的厚度）之间复杂的
3、关系，局部化的激光能量特征以及冷却速度使得很难监测这样的过程。能够准确模拟这个过程的方法对于制造高质量的覆层具有重要价值。其中一个被广泛采用的方法是有限元模型。在近几年，几个数学模型被提出来模拟激光熔覆过程。Sowdari和 Majumdar 研发出了一个有限元模型来分析金属的激光加热和融化过程。激光能量对于液-固交界面和熔池的影响得到了说明。Kumar和 Roy 提出了一个个模拟扫描速度、Marangoni 量对熔池和微观结构的影响的模型。扫描速度高时，熔池体积减少，在基体附近的冷却速度在增加，微观结构也更细。Roberts 等人使用出生和死亡元技术模拟了多个覆层的温度场，发现被加热的区域经
4、理快速的热循环。Tan 等人建立一个模型来估算该荣覆层的几何结构。通过控制扫描速度和粉末传递速率可以有效限制荣覆层的高度。Zhu 等人发现熔池的最高温度（1919）是低于热动力学温度（1927），这是为了能够获得羟磷灰石生物陶瓷的复合覆层。在这些文献中，这些参数的作用，例如热传导率，扫描速度，激光能量，和激光束的聚焦等这些加工特征都已经得到说明。在本文中，为了得到高功率激光熔覆工艺的物理现象，我们提出了用三维模型来预测涂层内特定区域的温度分布和三维形状与大小。除此此外，为了阐明在激光熔覆过程中的涂层相变和微观组织的演化过程，我们将基于温度曲线熔池的和理论的动力学来分析成核，生长速度，和 TiC
5、 颗粒的大小。最后将通过比较模拟结果与实验结果来验证模型的可靠性。建立数学模型我们使用商业有限元分析软件（ANSYS ver10.0），来计算在高功率激光熔覆过程中，随时间变化的温度分布。采用 SOLID70 单元来实体进行网格划分，这种划分得到三维热传导元素的八个节点和其在每一个节点的自由度，而这种元素适用于稳态或瞬态热分析。有一种生死单元技术，这种技术是利用激活和停用元素对时间的依赖性来模拟熔覆工艺。由于考虑到相变潜热，在热传导熔覆过程中的偏微分方程 7 ，该方程是其一般的形式表示，例如其中的 是密度，c 热容量，是预置层的导热系数，T 是温度分布函数，Q 是相变潜热，t 是时间。样品的底
6、面被设为绝热和分配绝热边界传导。样品的表面部分只有底面被设为与环境的对流边界。样品表面的边界条件通过该公式得出，其中 是对流换热系数， 是环境温度， 是试样hoTn表面法线方向。为了获得唯一解，该热传导偏微分方程应由牛顿-拉夫逊法求解。图 1 显示了这个过程的示意图。如图所示，当激光束照射到预置层的顶面，激光功率的一部分被反射，剩余部分被吸收。吸收的功率使预置层和基底的薄薄一层融化。在激光熔覆过程中，预置层通过冶金结合形成并包覆在基体表面。图一.激光熔覆过程示意图激光熔覆过程改变了激光照射的表面的边界形状并将预置层变为多孔的、薄弱的材料。除了考虑到热传导，对流（在较低的温度下的不熔化区）过程中
7、的热量流失，辐射（在熔化区接近高温度时）也应被考虑到，以此来得到在激光熔覆过程中正确的温度分布描述。这些因素都大大增加模型过程的难度。我们为了简化问题，采用了以下的假设： 整个样品的初始温度均假定为 20 摄氏度。在这个过程中，使用高斯公式来计算一个移动的激光束的强度分布：其中 是总的入射激光功率， 是预置层的吸收系数， 是激光Pov束， 扫描速度是激光束直径。r 为了使问题易于处理，这被认为是一个单一的复合轨道。有固定的参照系，因此这个过程可以看作是一个平稳过程。 将预置层假定为是均质且各向同性的。混合粉末的特性计算如下【10】：P( )代表的混合物 的热性能。 和 分别是组分的热性能和kf
8、kfkPf分数。 我们引入了吸收系数 来说明预置层的反射和散射所导致的激光能量损失 11-14 。它取决于表面条件 15 的性质。我们总是将 假定为一个常数，它的范围从 0.2 到 0.3 16 。在这个模型中， 为 0.25。 液态金属在熔池中的对流效应带来的影响是被我们忽视的。因为（a）热源的高度集中和移动（b）与熔池的样品相比是非常小的，（c）液体的对流速度热传导速度相比也非常慢，因此影响较小 17 。温度依赖于热容量，在表 1 和表 2 分别给出了覆层材料和基底的热导率和密度 。实验程序使用 20mm15mm5mm 大小的 Ti6Al4V 长方体样品来当作基体。将体积比为 1:3 的
9、TiC 和 NiCrBSiC 粉末混合物作为复合材料。其中NiCrBSiC 粉末的组成为 17.0%的 cr，3.5%的 B，4.0%的 si，1.0%的C，小于 12%的 Fe（质量百分比），剩余的是 Ni。混合粉末预先用有机粘合剂（聚乙烯醇和水）粘在基体表面，形成一层 0.6 毫米厚的涂层。激光熔覆是在一个 5 千瓦的二氧化碳激光系统中进行。工艺参数：激光功率 2500 W，扫描速度 4 毫米/秒，激光束的直径 4.5毫米。将样品切割，镶嵌，抛光，用 HF/HNO =1:15 腐蚀剂腐蚀，最3后用 Quanta 200 型扫描电子显微镜观察涂层的微观结构。4.结果与讨论4.1 温度分布的计
10、算结果如图 2 所示，在尺寸为 20mm15mm5mm 的基体上进行计算。为了保证计算效率和模拟精度，沿中心线的区域为致密网 19 剩余的为粗网。三维有限元素实体占有 42000 个单元和 47773 个节点。图 2.有限元网格模型在图 3 中有 7 个不同的温度曲线节点（2762，2785，2831，2877，2969，3107 和 3245）。这些节点位于沿激光扫描方向的涂层表面。所选择的节点之间的距离和节点数差异具有线性关系。当激光束从起始节点 2762 开始扫描时，峰值温度增加。然而，在节点 2969，3107 和 3245 上，峰值温度和温度曲线都几乎是稳定的。结果表明，在激光束开始
11、照射预置层后，热导率和温度分布在短距离内达到准静态。在高温期时，拥有更快的加热和冷却速率，表明了导热系数在较高的温度范围内更有效。图 4 显示了一个复合轨道的纵向截面的三维准稳态温度分布云图（在 0.2 秒熔覆时间内）。在激光照射区域存在较大的温度梯度，样品上的温度等温线呈现半椭圆形。激光束的前端相比激光束的末端，等温线更密集，温度梯度较大。由于最大模拟温度高于金属材料的熔点，因此熔池是在激光熔覆过程中形成的。图 3.在不同的节点温度曲线4.2.三维形状和特定区域的大小图 5 显示了基体的复合轨道纵断面（熔覆时间 2.0 秒内）的温度分布，TiC 熔化区的三维形状，熔池和热影响区。如图 6 所
12、示，为了深入了解激光熔覆过程中的物理现象，我们也对特定路径上的温度分布进行了研究。路径的位置如下：路径 1，2 和 3 分别定位在纵向截面并通过 TIC 熔化区的最深处，基体的熔池和热影响区。如图 5a 所示，路径 4 则位于沿激光扫描方向的涂层表面。根据温度分布，可以将三个特定区域按特征分为：TiC 熔化区，基体熔池和热影响区。在图 5a 黄色部分表示，在 TiC 熔点等温边界（3140 摄氏度）作为 TiC 熔化区 20 。这个区域在温度高于 3140摄氏度时位于涂层的上部。如图 5b，6a 和 d 所示，TiC 熔化区出现深度为 0.22 毫米且在激光扫描方向下长度为 1.4 毫米的凹陷
13、，.这表明涂层上层的 TiC 颗粒会在激光照射下融化。图 5a 绿色的部分表示，在 NiCrBSi 熔点等温边界（1668 摄氏度）作为熔池。熔池显示为深度为 0.76 毫米且在激光扫描方向长度为 2.87 毫米的的半椭圆形的形状，如图 5c，6b 和 d 所表示。熔池的深度略大于预置层的厚度（0.6 毫米）。这意味着，融化的是预置层和一层薄的基板（0.16mm），且基体和涂层之间的冶金结合可能是在激光熔覆过程中形成的。如图 5a 绿色的部分所示，在熔池中有一个区域的温度范围是 16683140之间。这表明，在激光照射期间，底层的 TiC颗粒部分溶解在了熔池。同上，我们用钛的同素异晶转变温度等温边界（882）来表征基体的热影响区，如在图 5a 蓝色部分所示。由图 5d，6b 和 d 可知，基体的焊接热影响区的形状也是深度为1.35 毫米且激光扫描方向长度为 4.83 毫米长的半椭圆形的形状。Fig.4 4. Fig.5
