1.一种激光熔覆过程中的元素分布预测模拟方法,其特征在于:包括以下步骤:(1)采用单元赋值方法离散计算域,设置空气域和基体域两个计算域,初始化计算条件;
(2)循环迭代开始,根据单元坐标和单元记录的空气体积分数标定计算域,设置激光扫描域、送粉域、热源域,分别赋值2、3和4,激光束域照在基体域上方,并沿着扫描方向移动,随着熔池的形成,凝固,再形成,基体域呈现梯度增长,能量吸收域和粉末沉积域的位置也随之发生变化,其中热源域中的热源视为均匀热源,其公式为:其中:q为能量密度,Q为激光功率,r为激光光斑半径;
(3)送粉域内的质量源根据单位时间送粉质量以静态液相形式记录在计算单元中,实时更新计算域;
(4)根据单元温度值和材料热物性参数求解凝固及熔化速率,当单元温度低于凝固点时,液相以柱状晶的形式凝固,柱状晶平均尺寸计算公式为:其中:dc为柱状晶直径,Rc为柱状晶半径,λ1为一次枝晶间距,fc-柱状晶体积分数;
柱状晶的凝固速率计算公式为:
Mlc=vRc·Sc·ρc·fl (3)其中: Mlc为柱状
晶凝固速率,vRc为柱状晶生长速率,Sc为柱状晶侧面积,ρc为柱状晶密度,fl为液相体积分数,t为单位时间,Dl为液相扩散系数, 为固液界面处液相溶质浓度, 为固液界面处固相溶质密度,cl为液相溶质浓度,Rf为柱状晶远场半径,dc为柱状晶直径;
(5)耦合求解液相和柱状晶相的质量、动量、能量、溶质的传输方程,其中质量守恒方程为:其中:ρl为液相密度,fc—柱状晶体积分数, 为液相流动速度;
动量守恒方程计算公式为:
其中:p为压强, 为液相应力-应变张量, 为重力加速度, 为凝固熔化导致的动量变化, 为固液间阻力导致的动量变化, 为固相流动速率, K为渗透率;
能量守恒方程计算公式为:
熔化:
凝固:
其中:hf为潜热,kl,kc分别是液相、固相的热导率,Tl、Tc分别是液相、固相得温度,hl、hc分别是液相、固相的焓, 分别是由相变导致的液相、固相的能量变化, 是由液相和固相之间的热传递导致的能量变化,H*为导热系数;
成分守恒方程计算公式为:
其中: 为相变导致的溶质浓度变化,cc为固相溶质浓度,k为溶质再分配系数;
(6)判断收敛条件,如果不满足进入下一个循环,如果满足进入下一个时间步计算;
(7)设计激光熔覆实验,采用金相制样制备熔覆层纵截面试样;
(8)将模拟结果与实际样品结果进行对照,验证模拟结果的可靠性,然后用于后续激光熔覆过程中元素分布的预测,通过能谱分析测量出试样的成分分布,与模拟结果进行对比,模拟结果的元素分布使用单元平均溶质表征,其计算公式如下:其中:cmix为单元平均溶质浓度。
2.如权利要求1所述的一种激光熔覆过程中的元素分布预测模拟方法,其特征在于:步骤(1)所述的设置空气域和基体域两个计算域,分别赋值1和5,空气域中的气相体积分数设为99.998%,基体域中固相体积分数设为99.998%,设置计算边界条件,基体域的边界类型设置为wall,边界温度为Tw=300K,热传导系数H=10W·m-2·K-1,空气域的边界类型设置为pressure outlet。
3.如权利要求1所述的一种激光熔覆过程中的元素分布预测模拟方法,其特征在于:步骤(7)设计激光熔覆实验,采用金相制样制备熔覆层纵截面试样,激光熔覆实验采用45钢为基体,316L为金属粉末,激光功率为1500W,扫描速度为10mm/s,送粉速率为0.6r/min;金相制样采用线切割,磨,抛,腐蚀的技术制备出熔覆层纵截面试样。