一种铸锭中宏观偏析缩孔疏松的预测模拟方法与流程

本发明涉及一种关于偏析的预测模拟方法，具体涉及一种铸锭中宏观偏析缩孔疏松的预测模拟方法，属于黑色金属材料制备领域。

背景技术：

凝固过程中成分偏析是合金固有特性。因此，凝固后铸锭普遍存在成分不均匀现象，这种现象叫做宏观偏析；宏观偏析不仅和材料特性还与凝固过程中液体流动、等轴晶沉积、凝固收缩流动等原因有关；此外凝固过程中还存在收缩现象，这是由材料在固、液状态下的密度差异引起的。铸锭凝固过程是一个复杂的过程，如何高效准确地预测铸锭的成分、缩孔缩松分布一直是各国铸造工作者努力的方向。

目前国内外用于铸锭偏析预测的模拟方法中，存在以下主要问题：很少同时考虑液相、等轴晶、柱状晶和空气的相互作用；没有考虑凝固收缩时的流动对最终偏析的影响；几乎没有跟踪缩松形成。

技术实现要素：

本发明的目的在于提供一种钢锭中的偏析预测、缩孔缩松预测；解决实际生产中凝固后成分、缩孔、疏松难以预测问题；模拟结果为实际生产提供有效的参考价值。本发明的技术方案如下：

一种钢锭中偏析的预测模拟方法，包括如下步骤：

(1)获取材料的热物性参数；

(2)计算求解所编写的凝固时所需的质量、动量、能量、溶质传输；

(3)设置计算边界条件，初始化计算条件，开始循环，耦合求解液相、等轴晶中固相和柱状晶相的质量、动量、能量、溶质的传输方程，同时求解空气相的动量及能量方程；

(4)跟踪空气相和液相的界面，根据凝固收缩量，求解出空气相分布域；

(5)更新材料物性参数，同时跟踪缩松判据，判断收敛条件，如果不满足进入下一个循环，如果满足进入下一个时间步计算；

(6)采用钢锭模实际生产钢锭；

(7)将模拟结果与实际生产结果进行对照，验证模拟结果的可靠性，然后用于后续铸锭的预测。

本发明的有益效果：

1、本发明采用四相模型，通过耦合求解凝固过程中金属液相、等轴晶相、柱状晶相以及气相之间的质量、能量、动量及溶质传输方程，预测铸锭凝固过程中的偏析、缩孔及疏松分布；

2、本发明使用温度梯度跟踪法确定柱状晶生长方向和生长区域，解决了柱状晶向等轴晶转变的预测难题；

3、本发明采用枝晶状而非球状等轴晶，解决了枝晶结构对预测等轴晶的生长及运动的难题；

4、本发明适用于预测各种材料、吨位铸锭的宏观成分分布，对铸锭浇铸工艺优化具有参考作用。

附图说明

图1是枝晶状等轴晶示意图；

图2是55吨钢模铸锭示意图；图中，1钢模；2铸锭；3空气；

图3是500s液相分布图；

图4是5000s液相分布图；

图5是10000s液相分布图；

图6是20000s液相分布图；

图7是500s缩孔形状图；

图8是5000s缩孔形状图；

图9是10000s缩孔形状图；

图10是20000s缩孔形状图；

图11是最终偏析值分布图；

图12是缩松分布图；

图13是中心线碳分布与实验值比较图。

具体实施方式

下面结合附图对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

一种钢锭中偏析的预测模拟方法，包括如下步骤：

(1)获取材料的热物性参数；

(2)计算求解所编写的凝固时所需的质量、动量、能量、溶质传输

(3)设置铸锭凝固的初始条件，包括初始温度，初始溶质含量，各相初始含量、速度；设置铸锭凝固的边界条件，包括铸锭与铸模的换热系数；

(4)开始迭代计算，计算每个单元的温度梯度，确定柱状晶、等轴晶的生长区域；

(5)通过fluent中的相耦合算法，耦合求解各相的质量、动量、能量、溶质的传输方程；

(6)跟踪固液界面预测凝固缩孔形成；

(7)跟踪缩松预测判据计算缩松量、更新热物性参数；

(8)判断是否收敛，如果满足条件进入下一个时间步长，如果不满足条件，从第(4)步开始重复迭代计算；

(9)计算结束，分析计算结果。

所述的钢锭中偏析预测方法的具体步骤如下：

1)等轴晶、柱状晶生长条件判定

只要相应的过冷度存在，等轴晶被允许在所有的单元中形核与长大；而柱状晶只从铸模壁开始生长，此后跟踪每个单元的温度梯度来判断柱状晶的生长方向，同时该单元的温度梯度满足一个阀值并且柱状晶前端沉积少于一定量的等轴晶时，柱状晶才能沿着热流反方向生长；

2)等轴晶、柱状晶生长速度及质量转换量计算

本发明所述的枝晶状等轴晶结构如图1所示，包络线所包围部分为一个枝晶，包含两部分：包络线内部固相fs，包络线内部液相等轴晶的生长包含两个部分：等轴晶的生长速度vRs，包络线的生长速度vtip。等轴晶枝生长速度计算如下：

等轴晶包络线的生长速度计算如下：

a＝0.4567；b＝1.195；k为溶质再分配系数；m相图中液相线斜率；为Γ为吉布斯-汤姆逊系数；

等轴晶包络线内部固相凝固速率计算如下：

Mls＝vRs·Ss·ρs·fl，

其中其中：fsi为等轴晶包络线内部固相率fsi＝fs/fenv；λ2为二次枝晶间距；

等轴晶包络线的凝固速率计算如下：

其中：fc柱状晶体积分数；

其中：Φe等轴晶球形度，根据合金体系确定；

柱状晶生长速率计算如下：

ρc为柱状晶密度

柱状晶凝固速率为：

Mlc＝vRc·Sc·ρc·fl

其中，并且有λ1为一次枝晶间距；

3)分别显性求解等轴晶(env)、等轴晶中固相(s)、柱状晶(c)、空气相(a)、金属液相(l)的质量、动量、能量及必要的溶质传输方程；使用质量守恒方程计算质量，使用动量守恒方程计算动量，使用能量守恒方程计算能量，使用成分守恒方程计算溶质传输。

质量守恒方程：

其中：Mls液相向等轴晶相的质量传输率，Mlc液相向柱状轴晶相的质量传输率；

动量守恒方程：

其中：P为压强；为应力-应变张量；

其中t为l，s，c三相的下角标；Ata＝6ftfa/dta；为l，s，c三相的速度；dta为t相与空气相的最小特征直径；μta为t相与空气相的最小粘度；

能量守恒方程：

其中：kl,ks,kc分别是液相、等轴晶中固相、柱状晶相的热导率；Tl,Ts,Tc分别是液相、等轴晶中固相、柱状晶相温度；

成分守恒方程：

其中：

3)更新热物性参数、缩松判据及判断收敛

粘度求解更新：

其中：μl为液相粘度，

单元平均溶质跟踪计算：

缩松判据跟踪：

其中SPC为缩松判据；Cλ为材料属性常数；G为单元的温度梯度；为单元的冷却速度；ΔPcr为压力下降阀值，大气中浇铸时为101325帕斯卡；μl为液相粘度；β凝固收缩率；ΔTf液固相线温度差。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。