首页 / 专利库 / 金工 / 冶金 / 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法

一种涉及冶金领域凝固过程中强制对流对AlN枝晶受计算方法

阅读:805发布:2022-12-28

专利汇可以提供一种涉及冶金领域凝固过程中强制对流对AlN枝晶受计算方法专利检索,专利查询,专利分析的服务。并且本 发明 提供一种涉及 冶金 领域 钢 液 凝固 过程中强制 对流 对AlN枝晶受 力 计算方法,涉及到冶金领域。该方法首先收集钢种成份以及凝固条件,然后计算凝固过程中 传热 与传质,凝固过程中界面胞的生长及凝固过程中AlN形核与生长,进而建立Fe-C-Al-N四元 合金 凝固中AlN析出模型。根据浇铸 温度 、溶质成分、冷却速率等 连铸 工艺条件,通过建立的AlN析出模型对其析出规律进行预测,并利用数据分析和 可视化 处理 软件 数据图像化显示AlN的析出 位置 、大小、形状和尺寸,以及定量化AlN析出数量。本发明提供的钢液凝固过程中AlN夹杂物析出情况的预测方法,为优化凝固技术、控制钢中AlN析出物尺寸和提高 铸坯 质量 提供了理论指导。,下面是一种涉及冶金领域凝固过程中强制对流对AlN枝晶受计算方法专利的具体信息内容。

1.一种涉及冶金领域凝固过程中强制对流对AlN枝晶受计算方法,其特征在于:
包括以下步骤:
步骤1收集钢种成分、凝固过程中热力学、动力学参数,以及流场边界条件;
步骤2基于金属凝固理论,采用元胞自动机模型计算界面生长方向、液相溶质场与温度场分布以及固相区溶质场分布;同时耦合格子Boltzmann方法中经典模型D2Q9模型去计算流场速度分布;
步骤2.1为使元胞自动机模型具有一定物理意义,引入金属凝固理论,采用尖锐界面模型去计算界面生长方向、曲率过冷
步骤2.2确定流体流动边界,采用D2Q9模型去计算计算流场的速度分布;
步骤2.3结合流场速度分布,采用元胞自动机模型去计算液相温度场;
步骤2.4结合流场速度分布,采用元胞自动机模型去计算液相区和固相区的溶质分布;
步骤3对于复杂边界流动,采用Mei修正F-H格式,沿x方向速度,在右侧添加自由边界条件,结合格子Boltzmann方法建立Fe-C-Al-N四元合金流场下枝晶受力生长模型。利用数据和分析和可视化处理软件图像显示AlN枝晶的形状、尺寸和受力情况。
2.根据权利要求1所述的一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于:所述步骤2.1的具体方法为:
首先假设固液界面处于热力学平衡状态,根据界面处溶质守恒定律,界面生长方向由如下公式计算:
其中 表示界面生长方向,曲率过冷K由如下公式计算:
其中, 和 分别表示固相率在x轴和y轴上的一阶偏导; 为固相率在先在x轴
上求偏导后对y轴求二阶偏导; 和 分别为固相率在x轴和y轴上的二阶偏导;
根据步骤2.2所描述,探测元胞状态,确定流体流动边界及其边界条件。在二维数值模拟计算中,通常采用D2Q9模型计算流场分布,同时忽略外力作用,流场的速度分布采用如下公式计算:
其中τ为无量纲单步松弛时间,fi(x+ciΔt,t+Δt)、fi(x,t)为液相流体粒子分布函数,fieq(x,t)为液相流体粒子平衡分布函数;其中ci为液相流体粒子在格子i方向的迁移速度,Δt为时间步长;由于采用D2Q9 模型,其中ci计算公式和为权重系数wi如下公式计算:
其中,c与cs分别为格子速度、格子声速;宏观密度ρ、宏观速度u以及流体动力学粘度v由如下公式计算:
其中fi为入口节点分布函数,i(取值为0,1,2,3,4,5,6,7,8)为步节点序数,入口密度ρin、出口密度ρout、流体粒子在x轴方向速度ux和流体粒子在y轴方向速度uy计算公式由如下式计算:
f1=f3,f2=f4,f8=f6           (19)
3.根据权利要求1所述的一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,步骤2.3结合流场速度分布,采用元胞自动机模型去计算液相温度场;具体方法由如下方法计算:
首先液相温度场的计算是传热过程中的计算,故液相区温度场采用如下公式计算:
其中,t为时间,ρ为基体密度,λ为导热系数,T为元胞温度,cp为基体比热容,fs为元胞固相率,L为潜热;qw,w、qw,e、qw,n、qw,s分别为左边界、右边界、上边界和下边界四个方向的热流密度;为了简化计算,假设基体在x轴和y轴方向上的导热系数λ相等;
固液界面处过冷度由成分过冷、曲率过冷、热过冷度和动力学过冷度组成,由于动力学过冷度在凝固过程中相对于其他过冷度较小,故忽略不计;因此,过冷度由如下公式计算:
ΔT=ΔTc+ΔTr+(Tbulk-TL)           (22)
其中,ΔTr为曲率过冷、Tbulk为凝固前沿的液相温度、TL为液相线温度、ΔTc为成分过冷;
液相线温度、曲率过冷和成分过冷分别由如下公式求解:
TL=1535-65[%C]-2.7[%Al]-90[%N]        (23)
其中,Γ为Gibbs-Thomson系数;mL,i为元素i的液相线斜率;cL,i表示元素i的液相浓度;
为元素i的固液界面前沿浓度;
由于在强制对流AlN枝晶受力计算方法中传热与传质密不可分,故根据步骤2.4结合流场速度分布,采用元胞自动机模型去计算液相区和固相区的溶质分布;液相中考虑了溶质之间的相互作用对枝晶的生长,由于固相中的溶质传输比液相中小几个量级,故忽略固相中的溶质相互作用,固相和液相的溶质传输分别采用如下公式进行计算:
其中,CL,i、CS,i分别表示液相和固相中元素i的浓度;CS,i为元素i在固相基体中的传输系数;i=1,2,3...n-1,n为钢液中元素的总个数,第n个元素表示溶剂; 表示液相中Darken系数矩阵,为简化计算假设基体中x轴和y轴方向 取值相同,利用下式求解:
其中,R表示气体常数,ak表示元素k的活度,xk表示元素k的摩尔分数,xj表示元素j的摩尔分数,δki表示Kronecter delta函数,当k=i时,δki取1,否则δki取0;Mk表示元素k在体系中的迁移率,其根据爱因斯坦公式求解,如下公式所示:
其中, 表示跟踪元素k的扩散系数。假设固液界面处于热力学平衡状态,界面处满足溶质分配定律:
根据界面处溶质守恒定律,在热力学平衡状态下界面胞的生长速度按下式求解:
其中,右边的 为溶质传输项,其只考虑液相溶质间相互作
用对枝晶生长的影响; 表示界面胞凝固前沿生长速度;通过求解出的界面胞凝固前沿生长速度,单位时间步长内界面胞的生长通过固相率增加来计算,如下公式所示:
其中, 和 分别为上一时刻和这一时刻AlN的元胞固相率Δfs,AlN的固相率的增加;Δt为单位时间步长; 表示沿 方向穿过元胞中心的单位长度;Δl=1μm为网格单元长度;θ表示枝晶生长方向和x轴方向的夹角。
4.根据权利要求3所述的一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于:
所述步骤2.4的具体方法为:对强制对流对AlN枝晶受力计算模拟计算以下假设:钢中析出的AlN不含有其他相;仅考虑液相中的AlN析出;忽略生长过程中的界面能增加;忽略AlN析出时的热量变化;基于以上假设,首先计算AlN在钢中形核的热力学和动力学条件;当满足形核条件时在元胞中形核并对周围溶质场进行更新计算,由于AlN界面处保持热力学平衡状态,溶质不断向界面处扩散,根据热力学平衡状态对强制对流对AlN枝晶受力计算模拟,固液界面法相与生长方向的函数 由如下公式计算:
其中,θ和 分别为界面生长法向和择优生长方向与x轴正方向的夹角。ε表示各项异性参数,其中由择优生长方向与x轴正方向的夹角 如下公式计算:
固液两相区生成AlN的反应公式和活度与活度系数计算公式由如下所示:
[Al]+[N]=(AlN)               (39)
a[M]=[%M]f[M]             (40)
其中,a[M]表示元素M的Henry活度;M表示元素Al或N,[M]表示元素M溶解在钢中;[%M]表示元素M基于质量1%标准的浓度;f[Al]和f[N]分别表示元素Al和N基于质量1%标准的活度系数;eij表示元素j对元素i的相互作用系数;当计算区域内某一元胞满足AlN析出条件时,标记该元胞并对该元胞进行3×3的网格细化分,同时根据动态化学平衡对其生长进行计算,如下公式所示:
MAlN表示AlN的相对分子质量;Δx表示反应量。
5.根据权利要求4所述的一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于采用Mei修正F-H格式,沿x方向速度,在右侧添加自由边界条件,结合格子Boltzmann方法建立Fe-C-Al-N四元合金流场下枝晶受力生长模型。
所以根据Mei修正F-H格式,首先对固相边界点xb,假设其有一个虚拟态平衡分布函数如下:
其中uf是xf处的流体速度,ubf为待定虚拟速度,为了求解ubf,构建一插值因子α,α与xw位置有关,ubf由如下公式求得:
其中,β是为了计算xf反弹分布构造的插值因子,uw是xw处界面移动速度;根据上述公式可求解虚拟态平衡分布函数;故Mei修正F-H格式下溶质场分布以及固相区溶质分布由如下公式求得:
所述根据动态化学平衡对元胞生长进行计算时,为减少由时间步长偏大带来的计算误差,对元胞进行空间上的3×3细化分时同时对其时间上进行细化分,在计算强制对流对AlN枝晶受力计算生长时在一个时间步长内的生长进行多次循环计算,以减少时间步长带来的计算误差;当基体元胞中AlN体积增加且接触到邻近液相胞时,邻近液相胞将被细划分为3×3的F-H格式下边界AlN析出元胞,继续AlN生长。利用数据和分析和可视化处理软件图像显示AlN枝晶的形状、尺寸和受力情况。

说明书全文

一种涉及冶金领域凝固过程中强制对流对AlN枝晶受

计算方法

技术领域

背景技术

[0002] 钢铁工业是国民经济重要支柱型基础产业。在钢液凝固的过程中,随着温度的降低铸坯表层开始形核生长形成柱状晶,随柱状晶的生长固液界面前沿溶质富集严重。不同尺寸的AlN夹杂物对钢材性能起到不同的作用。为此,钢液凝固过程中夹杂物析出预测对于控制铸坯裂纹,提高铸坯质量具有重要意义。连铸过程中,铸坯内的钢液会不断流动,对微观枝晶的生长和形貌产生一定影响。同时,钢液流动也会因为枝晶的存在而变得更加复杂。传统的流场数值模拟主要是通过Navier-Stokes(N-S)方程对压力场进行迭代求解。但该方法求解繁琐,计算量大,数值稳定性差,并且处理枝晶曲边边界时需要借助壁面函数,增大求解误差。因此,准确描述流场分布对枝晶求解极其重要。本文介绍一种近年兴起的格子Boltzmann方法来对枝晶周围的流场进行求解的方法。格子Boltzmann方法是源自于Boltzmann方程。 Boltzmann方程主要是用来描述系统非平衡态变化规律的方程。

发明内容

[0003] 本发明要解决的问题是针对现有凝固技术上的不足,提供一种钢液凝固过程中强制对流对AlN枝晶受力计算方法。预测在不同工艺的情况下强制对流对AlN枝晶受力情况,优化了凝固技术以及枝晶在强制对流下受力生长的理论指导。
[0004] 为解决上述技术问题,本发明所采取的技术方案是:一种钢液凝固过程中强制对流对AlN枝晶受力计算方法,包括以下步骤:
[0005] 步骤1收集钢种成分、凝固过程中热力学、动力学参数,以及流场边界条件;
[0006] 步骤2基于金属凝固理论,采用元胞自动机模型计算界面生长方向、液相溶质场与温度场分布以及固相区溶质场分布;同时耦合格子Boltzmann方法中经典模型D2Q9模型去计算流场速度分布;
[0007] 步骤2.1为使元胞自动机模型具有一定物理意义,引入金属凝固理论,采用尖锐界面模型去计算界面生长方向、曲率过冷
[0008] 步骤2.2确定流体流动边界,采用D2Q9模型去计算计算流场的速度分布;
[0009] 步骤2.3结合流场速度分布,采用元胞自动机模型去计算液相温度场;
[0010] 步骤2.4结合流场速度分布,采用元胞自动机模型去计算液相区和固相区的溶质分布;
[0011] 步骤3对于复杂边界流动,采用Mei修正F-H格式,沿x方向速度,在右侧添加自由边界条件,结合格子Boltzmann方法建立Fe-C-Al-N四元合金流场下枝晶受力生长模型。利用数据和分析和可视化处理软件图像显示AlN枝晶的形状、尺寸和受力情况。
[0012] 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于:所述步骤2.1的具体方法为:
[0013] 首先假设固液界面处于热力学平衡状态,根据界面处溶质守恒定律,界面生长方向 由如下公式计算:
[0014]
[0015] 其中 表示界面生长方向,曲率过冷K由如下公式计算:
[0016]
[0017] 其中, 和 分别表示固相率在x轴和y轴上的一阶偏导; 为固相率在先在x轴上求偏导后对y轴求二阶偏导; 和 分别为固相率在x轴和y轴上的二阶偏导;
[0018] 根据步骤2.2所描述,探测元胞状态,确定流体流动边界及其边界条件。在二维数值模拟计算中,通常采用D2Q9模型计算流场分布,同时忽略外力作用,流场的速度分布采用如下公式计算:
[0019]
[0020]
[0021] 其中τ为无量纲单步松弛时间,fi(x+ciΔt,t+Δt)、fi(x,t)为液相流体粒子分布函数, 为液相流体粒子平衡分布函数;其中ci为液相流体粒子在格子i方向的迁移速度,Δt为时间步长;由于采用D2Q9 模型,其中ci计算公式和为权重系数wi如下公式计算:
[0022]
[0023]
[0024]
[0025] 其中,c与cs分别为格子速度、格子声速;宏观密度ρ、宏观速度u以及流体动力学粘度v由如下公式计算:
[0026]
[0027]
[0028]
[0029] 其中fi为入口节点分布函数,i(取值为0,1,2,3,4,5,6,7,8)为步节点序数,入口密度ρin、出口密度ρout、流体粒子在x轴方向速度ux和流体粒子在y轴方向速度uy计算公式由如下式计算:
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038] f1=f3,f2=f4,f8=f6  (19)
[0039]
[0040] 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,步骤2.3结合流场速度分布,采用元胞自动机模型去计算液相温度场;具体方法由如下方法计算:
[0041] 首先液相温度场的计算是传热过程中的计算,故液相区温度场采用如下公式计算:
[0042]
[0043] 其中,t为时间,ρ为基体密度,λ为导热系数,T为元胞温度, cp为基体比热容,fs为元胞固相率,L为潜热;qw,w、qw,e、qw,n、qw,s分别为左边界、右边界、上边界和下边界四个方向的热流密度;为了简化计算,假设基体在x轴和y轴方向上的导热系数λ相等;
[0044] 固液界面处过冷度由成分过冷、曲率过冷、热过冷度和动力学过冷度组成,由于动力学过冷度在凝固过程中相对于其他过冷度较小,故忽略不计;因此,过冷度由如下公式计算:
[0045] ΔT=ΔTc+ΔTr+(Tbulk-TL)  (22)
[0046] 其中,ΔTr为曲率过冷、Tbulk为凝固前沿的液相温度、TL为液相线温度、ΔTc为成分过冷;液相线温度、曲率过冷和成分过冷分别由如下公式求解:
[0047] TL=1535-65[%C]-2.7[%Al]-90[%N]  (23)
[0048]
[0049]
[0050] 其中,Γ为Gibbs-Thomson系数;mL,i为元素i的液相线斜率; cL,i表示元素i的液相浓度; 为元素i的固液界面前沿浓度;
[0051] 由于在强制对流AlN枝晶受力计算方法中传热与传质密不可分,故根据步骤2.4结合流场速度分布,采用元胞自动机模型去计算液相区和固相区的溶质分布;液相中考虑了溶质之间的相互作用对枝晶的生长,由于固相中的溶质传输比液相中小几个量级,故忽略固相中的溶质相互作用,固相和液相的溶质传输分别采用如下公式进行计算:
[0052]
[0053]
[0054] 其中,CL,i、CS,i分别表示液相和固相中元素i的浓度;CS,i为元素 i在固相基体中的传输系数;i=1,2,3...n-1,n为钢液中元素的总个数,第n个元素表示溶剂; 表示液相中Darken系数矩阵,为简化计算假设基体中x轴和y轴方向 取值相同,利用下式求解:
[0055]
[0056]
[0057]
[0058] 其中,R表示气体常数,ak表示元素k的活度,xk表示元素k的摩尔分数,xj表示元素j的摩尔分数,δki表示Kronecter delta函数,当k=i时,δki取1,否则δki取0;Mk表示元素k在体系中的迁移率,其根据爱因斯坦公式求解,如下公式所示:
[0059]
[0060] 其中, 表示跟踪元素k的扩散系数。假设固液界面处于热力学平衡状态,界面处满足溶质分配定律:
[0061]
[0062] 根据界面处溶质守恒定律,在热力学平衡状态下界面胞的生长速度按下式求解:
[0063]
[0064] 其中,右边的 为溶质传输项,其只考虑液相溶质间相互作用对枝晶生长的影响; 表示界面胞凝固前沿生长速度;通过求解出的界面胞凝固前沿生长速度,单位时间步长内界面胞的生长通过固相率增加来计算,如下公式所示:
[0065]
[0066]
[0067]
[0068] 其中, 和 分别为上一时刻和这一时刻AlN的元胞固相率Δfs,AlN的固相率的增加;Δt为单位时间步长; 表示沿 方向穿过元胞中心的单位长度;Δl=1μm为网格单元长度;θ表示枝晶生长方向和x轴方向的夹角。
[0069] 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于:
[0070] 所述步骤2.4的具体方法为:对强制对流对AlN枝晶受力计算模拟计算以下假设:钢中析出的AlN不含有其他相;仅考虑液相中的 AlN析出;忽略生长过程中的界面能增加;
忽略AlN析出时的热量变化;基于以上假设,首先计算AlN在钢中形核的热力学和动力学条件;当满足形核条件时在元胞中形核并对周围溶质场进行更新计算,由于AlN界面处保持热力学平衡状态,溶质不断向界面处扩散,根据热力学平衡状态对强制对流对AlN枝晶受力计算模拟,固液界面法相与生长方向的函数 由如下公式计算:
[0071]
[0072] 其中,θ和 分别为界面生长法向和择优生长方向与x轴正方向的夹角。ε表示各项异性参数,其中由择优生长方向与x轴正方向的夹角 如下公式计算:
[0073]
[0074] 固液两相区生成AlN的反应公式和活度与活度系数计算公式由如下所示:
[0075] [Al]+[N]=(AlN)  (39)
[0076] a[M]=[%M]f[M]  (40)
[0077]
[0078] 其中,a[M]表示元素M的Henry活度;M表示元素Al或N,[M] 表示元素M溶解在钢中;[%M]表示元素M基于质量1%标准的浓度; f[Al]和f[N]分别表示元素Al和N基于质量1%标准的活度系数; 表示元素j对元素i的相互作用系数;当计算区域内某一元胞满足AlN析出条件时,标记该元胞并对该元胞进行3×3的网格细化分,同时根据动态化学平衡对其生长进行计算,如下公式所示:
[0079]
[0080] MAlN表示AlN的相对分子质量;Δx表示反应量;
[0081] 一种涉及冶金领域钢液凝固过程中强制对流对AlN枝晶受力计算方法,其特征在于采用Mei修正F-H格式,沿x方向速度,在右侧添加自由边界条件,结合格子Boltzmann方法建立Fe-C-Al-N四元合金流场下枝晶受力生长模型。
[0082] 所以根据Mei修正F-H格式,首先对固相边界点xb,假设其有一个虚拟态平衡分布函数如下:
[0083]
[0084] 其中uf是xf处的流体速度,ubf为待定虚拟速度,为了求解ubf,构建一插值因子α,α与xw位置有关,ubf由如下公式求得:
[0085]
[0086]
[0087]
[0088] 其中,β是为了计算xf反弹分布构造的插值因子,uw是xw处界面移动速度;根据上述公式可求解虚拟态平衡分布函数;故Mei修正 F-H格式下溶质场分布以及固相区溶质分布由如下公式求得:
[0089]
[0090]
[0091] 所述根据动态化学平衡对元胞生长进行计算时,为减少由时间步长偏大带来的计算误差,对元胞进行空间上的3×3细化分时同时对其时间上进行细化分,在计算强制对流对AlN枝晶受力计算生长时在一个时间步长内的生长进行多次循环计算,以减少时间步长带来的计算误差;当基体元胞中AlN体积增加且接触到邻近液相胞时,邻近液相胞将被细划分为3×3的F-H格式下边界AlN析出元胞,继续AlN 生长。利用数据和分析和可视化处理软件图像显示AlN枝晶的形状、尺寸和受力情况。
[0092] 采用上述技术方案采用上述技术方案所产生的有益效果在于:本发明提供的一种钢液凝固过程中强制对流对AlN枝晶受力计算方法。在计算强制对流下枝晶生长,通过对时间步长的细化提高了计算精度。显示AlN枝晶形状、尺寸和受力情况。为发明优化凝固技术,对枝晶在强制对流下受力分析进行预测,提高铸坯质量提供理论指导。附图说明
[0093] 图1为本发明实施例提供的一种钢液凝固过程中强制对流对AlN 枝晶受力计算方法的流程图
[0094] 图2为本发明实施例提供的强制对流对AlN枝晶受力计算程序建立的流程图;
[0095] 图3为本发明实施例提供的格子Boltzmann方法常用D2Q9模型示意图的示意图;
[0096] 图4为本发明实施例提供的Mei修正F-H格式示意图;
[0097] 图5为本发明实施例提供的流场中枝晶形貌演变和Al溶质分布模拟图;

具体实施方式

[0098] 下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
[0099] 本实例以普通低钢为为例,采用本发明的钢液凝固过程中强制对流对AlN枝晶受力计算方法对该钢液凝固过程中的强制对流对 AlN枝晶受力计算。
[0100] 一种钢液凝固过程中强制对流对AlN枝晶受力计算方法,如图1 所示,包括以下步骤:
[0101] 步骤1收集钢种成分、凝固条件及边界条件;本实施例中,采用的是38CrMoAl钢。其成分如表1所示,凝固条件包括冷速、过冷度及浇铸温度;
[0102] 表1钢种主要成分
[0103]
[0104] 步骤2基于金属凝固理论,采用元胞自动机模型计算界面生长方向、液相溶质场与温度场分布以及固相区溶质场分布;同时耦合格子Boltzmann方法中经典模型D2Q9模型去计算流场速度分布;具体方法为:
[0105] 步骤2.1为使元胞自动机模型具有一定物理意义,引入金属凝固理论,采用尖锐界面模型去计算界面生长方向、曲率过冷;
[0106] 首先假设固液界面处于热力学平衡状态,根据界面处溶质守恒定律,界面生长方向由如下公式计算:
[0107]
[0108] 其中 表示界面生长方向,曲率过冷K由如下公式计算:
[0109]
[0110] 其中, 和 分别表示固相率在x轴和y轴上的一阶偏导; 为固相率在先在x轴上求偏导后对y轴求二阶偏导; 和 分别为固相率在x轴和y轴上的二阶偏导;
[0111] 步骤2.2确定流体流动边界,采用D2Q9模型去计算计算流场的速度分布;确定流体流动边界及其边界条件。在二维数值模拟计算中,通常采用D2Q9模型计算流场分布,同时忽略外力作用,其液相流体粒子计算表达式和平衡态分布函数采用如下公式计算:
[0112]
[0113]
[0114] 其中τ为无量纲单步松弛时间,fi(x+ciΔt,t+Δt)、fi(x,t)为液相流体粒子分布函数, 为液相流体粒子平衡分布函数;其中ci为液相流体粒子在格子i方向的迁移速度,Δt为时间步长;由于采用D2Q9 模型,其中ci计算公式和为权重系数wi如下公式计算:
[0115]
[0116]
[0117]
[0118] 其中,c与cs分别为格子速度、格子声速;宏观密度ρ、宏观速度u以及流体动力学粘度v由如下公式计算:
[0119]
[0120]
[0121]
[0122] 其中fi为入口节点分布函数,i(取值为0,1,2,3,4,5,6,7,8)为角步节点序数,入口密度ρin、出口密度ρout、流体粒子在x轴方向速度ux和流体粒子在y轴方向速度uy计算公式由如下式计算:
[0123]
[0124]
[0125]
[0126]
[0127]
[0128]
[0129]
[0130]
[0131] f1=f3,f2=f4,f8=f6  (19)
[0132]
[0133] 步骤2.3结合流场速度分布,采用元胞自动机模型去计算液相温度场;具体方法由如下方法计算:
[0134] 首先液相温度场的计算是传热过程中的计算,故液相区温度场采用如下公式计算:
[0135]
[0136] 其中,t为时间,ρ为基体密度,λ为导热系数,T为元胞温度, cp为基体比热容,fs为元胞固相率,L为潜热;qw,w、qw,e、qw,n、qw,s分别为左边界、右边界、上边界和下边界四个方向的热流密度;为了简化计算,假设基体在x轴和y轴方向上的导热系数λ相等;
[0137] 固液界面处过冷度由成分过冷、曲率过冷、热过冷度和动力学过冷度组成,由于动力学过冷度在凝固过程中相对于其他过冷度较小,故忽略不计;因此,过冷度由如下公式计算:
[0138] ΔT=ΔTc+ΔTr+(Tbulk-TL)  (22)
[0139] 其中,ΔTr为曲率过冷、Tbulk为凝固前沿的液相温度、TL为液相线温度、ΔTc为成分过冷;液相线温度、曲率过冷和成分过冷分别由如下公式求解:
[0140] TL=1535-65[%C]-2.7[%Al]-90[%N]  (23)
[0141]
[0142]
[0143] 其中,Γ为Gibbs-Thomson系数;mL,i为元素i的液相线斜率; cL,i表示元素i的液相浓度; 为元素i的固液界面前沿浓度;
[0144] 步骤2.4结合流场速度分布,采用元胞自动机模型去计算液相区和固相区的溶质分布;具体方法如下所示:
[0145] 液相中考虑了溶质之间的相互作用对枝晶的生长,由于固相中的溶质传输比液相中小几个量级,故忽略固相中的溶质相互作用,固相和液相的溶质传输分别采用如下公式进行计算:
[0146]
[0147]
[0148] 其中,CL,i、CS,i分别表示液相和固相中元素i的浓度;CS,i为元素 i在固相基体中的传输系数;i=1,2,3...n-1,n为钢液中元素的总个数,第n个元素表示溶剂; 表示液相中Darken系数矩阵,为简化计算假设基体中x轴和y轴方向 取值相同,利用下式求解:
[0149]
[0150]
[0151]
[0152] 其中,R表示气体常数,ak表示元素k的活度,xk表示元素k的摩尔分数,xj表示元素j的摩尔分数,δki表示Kronecter delta函数,当k=i时,δki取1,否则δki取0;Mk表示元素k在体系中的迁移率,其根据爱因斯坦公式求解,如下公式所示:
[0153]
[0154] 其中, 表示跟踪元素k的扩散系数。假设固液界面处于热力学平衡状态,界面处满足溶质分配定律:
[0155]
[0156] 根据界面处溶质守恒定律,在热力学平衡状态下界面胞的生长速度按下式求解:
[0157]
[0158] 其中,右边的 为溶质传输项,其只考虑液相溶质间相互作用对枝晶生长的影响; 表示界面胞凝固前沿生长速度;通过求解出的界面胞凝固前沿生长速度,单位时间步长内界面胞的生长通过固相率增加来计算,如下公式所示:
[0159]
[0160]
[0161]
[0162] 其中, 和 分别为上一时刻和这一时刻AlN的元胞固相率Δfs,AlN的固相率的增加;Δt为单位时间步长; 表示沿 方向穿过元胞中心的单位长度;Δl=1μm为网格单元长度;θ表示枝晶生长方向和x轴方向的夹角。
[0163] 固液界面法相与生长方向的函数 由如下公式计算:
[0164]
[0165] 其中,θ和 分别为界面生长法向和择优生长方向与x轴正方向的夹角。ε表示各项异性参数,其中由择优生长方向与x轴正方向的夹角 如下公式计算:
[0166]
[0167] 固液两相区生成AlN的反应公式和活度与活度系数计算公式由如下所示:
[0168] [Al]+[N]=(AlN)  (39)
[0169] a[M]=[%M]f[M]  (40)
[0170]
[0171] 其中,a[M]表示元素M的Henry活度;M表示元素Al或N,[M] 表示元素M溶解在钢中;[%M]表示元素M基于质量1%标准的浓度; f[Al]和f[N]分别表示元素Al和N基于质量1%标准的活度系数; 表示元素j对元素i的相互作用系数;(如表2所示);
[0172] 表2溶质相互作用系数
[0173]
[0174] 当计算区域内某一元胞满足AlN析出条件时,标记该元胞并对该元胞进行3×3的网格细化分,同时根据动态化学平衡对其生长进行计算,如下公式所示:
[0175]
[0176] MAlN表示AlN的相对分子质量;Δx表示反应量;
[0177] 对于复杂边界条件,其特征在于采用Mei修正F-H格式,沿x 方向速度,在右侧添加自由边界条件,结合格子Boltzmann方法建立 Fe-C-Al-N四元合金流场下枝晶受力生长模型。
[0178] 所以根据Mei修正F-H格式,首先对固相边界点xb,假设其有一个虚拟态平衡分布函数如下:
[0179]
[0180] 其中uf是xf处的流体速度,ubf为待定虚拟速度,为了求解ubf,构建一插值因子α,α与xw位置有关,ubf由如下公式求得:
[0181]
[0182]
[0183]
[0184] 其中,β是为了计算xf反弹分布构造的插值因子,uw是xw处界面移动速度;根据上述公式可求解虚拟态平衡分布函数;故Mei修正 F-H格式下溶质场分布以及固相区溶质分布由如下公式求得:
[0185]
[0186]
[0187] 所述根据动态化学平衡对元胞生长进行计算时,为减少由时间步长偏大带来的计算误差,对元胞进行空间上的3×3细化分时同时对其时间上进行细化分,在计算强制对流对AlN枝晶受力计算生长时在一个时间步长内的生长进行多次循环计算,以减少时间步长带来的计算误差;当基体元胞中AlN体积增加且接触到邻近液相胞时,邻近液相胞将被细划分为3×3的F-H格式下边界AlN析出元胞,继续AlN 生长。利用数据和分析和可视化处理软件图像显示AlN枝晶的形状、尺寸和受力情况。
[0188] 本实施例中,计算过程中涉及到的模型参数值如表3所示:
[0189] 表3计算过程中的模型参数
[0190]
[0191]
[0192] 实施例基于Visual Studio 2015平台运用C++语言对强制对流对 AlN枝晶受力计算数学模型编写如图2所示的数值模拟程序实现,得到如图5所示的钢液凝固流场中枝晶形貌演变和Al溶质分布模拟图;数值模拟得到流场下枝晶生长模型,为优化凝固技术,对枝晶在强制对流下受力分析进行预测,提高铸坯质量提供理论指导。
[0193] 最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
相关专利内容
标题 发布/更新时间 阅读量
粉末冶金方法 2020-05-12 3
粉末冶金工艺 2020-05-12 352
一种粉末冶金模具 2020-05-13 595
一种粉末冶金烧结炉 2020-05-13 294
粉末冶金方法 2020-05-11 446
粉末冶金工艺 2020-05-11 106
粉末冶金方法 2020-05-11 594
粉末冶金方法 2020-05-12 842
粉末冶金方法 2020-05-11 410
冶金设备 2020-05-11 883
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈