CN107229803A - 一种预测连铸生产的铸坯微观组织的方法 - Google Patents

一种预测连铸生产的铸坯微观组织的方法 Download PDF

Info

Publication number
CN107229803A
CN107229803A CN201710459388.5A CN201710459388A CN107229803A CN 107229803 A CN107229803 A CN 107229803A CN 201710459388 A CN201710459388 A CN 201710459388A CN 107229803 A CN107229803 A CN 107229803A
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
steel
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710459388.5A
Other languages
English (en)
Inventor
罗森
王朝辉
陈耀
赵群
王卫领
朱苗勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN201710459388.5A priority Critical patent/CN107229803A/zh
Publication of CN107229803A publication Critical patent/CN107229803A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Continuous Casting (AREA)

Abstract

本发明提供一种预测连铸生产的铸坯微观组织的方法,方法包括:收集使用待分析钢种生产铸坯的连铸设备的设备参数和工艺参数;针对待分析的钢种参数,确定该钢种的热物性参数;利用预先建立的宏微观多尺度数学模型结合所述设备参数、工艺参数和所述钢热物性参数,预测当前连铸设备对应的铸坯的组织信息。上述方法能够较好的预测连铸坯凝固组织,为优化连铸结晶器工艺条件提供直观可靠的信息。

Description

一种预测连铸生产的铸坯微观组织的方法
技术领域
本发明涉及的是冶金连铸技术,具体涉及一种预测连铸生产的铸坯微观组织的方法。
背景技术
钢连铸过程实质上是钢液连续冷却凝固成型过程,当钢液在结晶器水冷作用下,钢液与结晶器表面接触位置大量形核,并逐渐向铸坯中心生长。铸坯凝固组织生长的同时伴随着钢中溶质在不同相之间的重新分配,由于钢中溶质元素在固相中的溶解度低于在液相中的溶解度,所以凝固过程中,枝晶从铸坯表面向中心生长的同时伴随着溶质元素向枝晶臂壁间富集,如果富集在枝晶臂间的溶质元素在凝固末期不能有效扩散,将会造成短距离(10~100μm)溶质成分波动形成枝晶间微观偏析。
连铸过程中钢液流动促使枝晶间溶质元素发生长距离的迁徙,将会引起铸坯较大范围内的溶质成分波动进而形成铸坯宏观偏析。同时,铸坯凝固组织生长同时伴随着凝固收缩,如果枝晶过度生长,造成枝晶搭桥,使得凝固收缩得不到钢液补充,将会形成疏松、缩孔等缺陷。这些连铸坯凝固缺陷在后续热加工处理过程得不到有效消除,将会影响最终钢产品性能。
由此可见,连铸坯凝固组织与铸坯缺陷之间有直接关系,且对最终钢产品质量产生重要影响。为此,连铸坯凝固微观组织预测对于优化连铸工艺,控制连铸坯凝固组织,提高铸坯质量至关重要。
现有技术中的连铸坯凝固组织检验方法是判断连铸坯质量的常用手段,主要有硫印检验方法、热酸蚀检验方法、冷酸蚀检验方法、电解腐蚀检验方法和枝晶腐蚀检验方法等。这些连铸坯凝固组织检测方法主要针对已经生产连铸坯进行取样检测分析,需要一定的试样加工和分析时间,检测周期较长,成本较大,劳动强度大,环境污染,对人体健康有伤害。为此,现有的连铸坯凝固组织检测大多作为连铸坯生产过程随检手段,检测样本少,难以真实反映大生产过程不同连铸工艺条件下连铸坯凝固组织,且检测周期长,滞后性严重,难以用于实时指导连铸工业生产。
另外,现有技术中的数值模拟作为一种高效低成本手段,越来越受到冶金工作者的重视。数值模拟预测法基于钢连铸凝固传热机理,结合凝固形核原理,建立钢连铸凝固组织长数学模型,预测连铸钢液凝固过程。该方法较试验测量法具有成本低、效率高、分析全面等优点。由于钢连铸过程是个复杂的高温钢液凝固过程,涉及到宏微观多物理场的耦合计算,现有的钢连铸坯凝固组织数学模型,大多基于ProCAST商业软件的模块所建立,该模型移植性较差,且后期二次开发较困难,难以满足钢连铸生产过程工艺优化和凝固组织控制方面的要求。
发明内容
为解决现有技术中的问题,本发明提供一种预测连铸生产的铸坯微观组织的方法。
第一方面,本发明提供一种预测连铸生产的铸坯微观组织的方法,包括:
收集使用待分析钢种生产铸坯的连铸设备的设备参数和工艺参数;
针对待分析的钢种参数,确定该钢种的热物性参数;
利用预先建立的宏微观多尺度数学模型结合所述设备参数、工艺参数和所述钢热物性参数,预测当前连铸设备对应的铸坯的组织信息。
可选地,针对待分析的钢种参数,确定该钢种的热物性参数,包括:
依据枝晶间溶质守恒原理和枝晶间溶质扩散原理建立枝晶间溶质微观偏析数学模型;
根据待分析的钢种成分、冷却强度,通过所述枝晶间溶质微观偏析数学模型,获取钢液凝固过程中枝晶间溶质偏析和相变信息;
根据所述枝晶间溶质偏析和相变信息,获取钢种的热物性参数;
所述热物性参数包括下述的一种或多种:液相线温度、固相线温度、导热系数、密度、比热和凝固潜热数据。
可选地,根据所述枝晶间溶质偏析和相变信息,获取钢种的热物性参数,包括:
根据公式一,确定钢液的液相线温度Tl
公式一:
根据公式二,确定钢液凝固过程中δ/γ相转变开始温度TAr4
公式二:
根据公式三,确定钢液的导热系数k;
公式三:k=kαfα+kδfδ+kγfγ+klfl
根据公式四,确定钢液的钢密度ρ;
公式四:ρ=ργfγδfδlfl
根据公式五,确定钢液的凝固潜热数据c;
公式五:c=cαfα+cδfδ+cγfγ+clfl
其中,Cδ l,i为溶质元素i在钢液的δ/l相界面处液相中的浓度;mi和ni分别为伪二元Fe-i相图中液相线和Ar4线斜率;C0 l,i为元素i在钢液中的初始浓度,%;ki δ/l为溶质元素i在δ/l界面平衡分配系数;
kα,kδ,kγ,kl分别为α相,δ相,γ相和l相导热系数;fα,fδ,fγ,fl分别为α相,δ相,γ相和l相所占的分率;cα,cδ,cγ,cl分别为α相,δ相,γ相和l相比热容;ρα,ρδ,ργ,ρl分别为α相,δ相,γ相和l相密度。
可选地,所述工艺参数包括下述的一种或多种:
浇铸温度、拉速、钢种成分、结晶器水量、温升和二冷各区水量;
和/或,
所述设备参数包括下述的一种或多种:
结晶器长度、结晶器的有效高度、结晶器水槽尺寸与分布、铜板厚度、二冷区长度、二冷喷嘴布置;
和/或,
所述铸坯的组织信息包括下述的一种或多种:
定量化地表征晶粒尺寸、定量化地表征柱状晶向等轴晶转变位置、定量化地表征等轴晶率和定量化地表征铸坯不同区域的坯壳厚度。
可选地,利用预先建立的宏微观多尺度数学模型结合所述设备参数、工艺参数和所述钢热物性参数,预测当前连铸设备对应的铸坯的组织信息的步骤之前,所述方法还包括:
依据钢连铸过程的传热原理及生产铸坯的各连铸设备的传热信息,建立宏微观多尺度数学模型。
可选地,所述宏微观多尺度数学模型包括:
依据钢连铸过程的传热原理及生产铸坯的各连铸设备的传热信息建立的宏观传热模型;
在预测过程中使用宏观传热模型的预测结果对铸坯组织信息预测的微观枝晶生长模型。
可选地,所述宏观传热模型包括:
定义钢连铸过程的二维凝固传热控制方式为:
式中:T为温度,℃;ρ为密度,kg/m3;cp为热容,J/(kg·℃);keff为导热系数,W/(m·℃);L为凝固潜热,J/kg;fs为固相分率;t为时间,s;x和y分别为铸坯横截面宽面和窄面方向长度,m;
定义结晶器冷却区的结晶器传热为:
式中:q为热流密度,W/m2;z为离弯月面的距离,m;vc为拉速,m/s;
定义连铸二冷区的传热表示对流传热和辐射传热之和,具体为:
q=h(Tsurf-Tamb)+σε[(Tsurf+273.15)4-(Tamb+273.15)4];
式中:h为综合对流换热系数,W/(m2·℃);Tsurf为铸坯表面温度,℃;Tamb为环境温度,℃;σ为Stefan-Boltzman常数5.67×10-8W/(m2·K4);ε为铸坯表面黑度系数;
定义空冷区辐射换热为:
q=σε[(Tsurf+273.15)4-(Tamb+273.15)4];
定义铸坯不同区域/位置处冷却速率CRi为:
式中:Tl为液相线温度,℃;Ts为固相线温度,℃;△t为凝固时间。
可选地,所述微观枝晶生长模型包括:
固液两相中的溶质扩散控制过程表示为:
式中:t为时间,Dl和Ds分别为固液两相中溶质扩散系数,fs为固相分率,k为溶质平衡分配系数;
枝晶尖端过冷度ΔT表示为:
式中:C0为合金初始浓度,Cl *为固液界面处液相溶质浓度,m0为液相线斜率,Г为Gibbs-Thomson系数,κ为固液界面曲率,为界面能各向异性函数,θ为枝晶择优生长角度,为固液界面法相与x轴的夹角;
固液界面曲率κ由固液界面处固相率梯度表示:
固液相界面张力的各向异性表示为:
枝晶尖端生产速度受固液相界面处溶质浓度控制,表示为:
式中:(fs)x和(fs)y为固相率的一阶偏导,(fs)xx,(fs)xy和(fs)yy为固相率的二阶偏导,ε为界面各向异性参数;
vn为界面法相移动速度,k为溶质平衡分配系数,Ds为溶质元素在固相中扩散系数,Dl为溶质元素在液相中扩散系数。
本发明具有的有益效果如下:
本发明的方法为了保证计算效率,在连铸宏观凝固传热的基础上,耦合钢液内部形核和晶粒生长在内的微观传输现象,建立连铸坯凝固组织的宏微观多尺度数学模型,实现建立的模型预测不同连铸工艺条件对连铸坯凝固组织的影响,为优化连铸工艺,控制连铸坯凝固组织提供了理论指导。
此外,宏微观多尺度数学模型内所涉及钢种物性参数均由溶质微观偏析模型计算所得,避免了现有的传统模型采用常数和简单的经验回归公式处理钢种热物性参数(例如液相线温度、固相线温度、导热系数、密度、比热、凝固潜热等),无法表征钢液凝固过程溶质偏析和相变对钢种热物性参数的影响的弊端,显著提高了模型计算精度。
附图说明
图1为本发明提供的预测连铸生产的铸坯微观组织的方法的功能结构图;
图2为本发明提供的宏微观多尺度数学模型中使用的枝晶形貌示意图;
图3为本发明提供的宏微观多尺度数学模型中连铸宏观传热计算的示意图;
图4为本发明提供的预测连铸生产的铸坯微观组织的方法的流程图;
图5(a)为实测重轨钢U75V连铸坯的组织信息的示意图;
图5(b)为预测重轨钢U75V连铸坯的组织信息的示意图。
具体实施方式
为了更好的解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
结合图1至图3所示,本实施例的预测方法可由热物性参数计算、宏观凝固传热计算、凝固组织生长计算、结果输出四部分内容组成。
(1)热物性参数计算
图2示出了钢连铸坯枝晶生长示意图。在图2中,假设枝晶形貌横截面为正六边形,首先依据枝晶间溶质守恒原理和溶质扩散原理建立溶质微观偏析数学模型;
然后根据钢种成分和冷却强度,结合钢液凝固过程中枝晶间溶质偏析数学模型,计算钢液凝固过程中枝晶间溶质偏析和相变信息,并在此基础上获取钢种的热物性参数,例如计算液相线温度、固相线温度、导热系数、密度、比热、凝固潜热、等钢种热物性参数,从而为宏微观多尺度数学模型提供准确地钢种物性参数。
溶质微观偏析数学模型中溶质扩散采用如下公式计算:
初始条件:当t=0时,
边界条件:当x=0,λ/2时,
以上各式中,λ为二次枝晶间距,m;C0 l,i,Cl,i,Cs,i分别为元素i在钢液中的初始浓度,液相l和固相s中溶质浓度,%;Ds,i(T)为溶质元素i在固相s中的扩散系数,m2/s;t为时间,s;ki s/l,ki δ/l,ki γ/l分别为溶质元素i在s/l,δ/l,γ/l界面平衡分配系数,具体物性参数见表1。
钢液相线温度TL和凝固过程中δ/γ相转变开始温度TAr4分别采用如下表达式:
式中:Cδ l,i为元素i在δ/l相界面处液相中的浓度;mi和ni分别为伪二元Fe-i相图中液相线和Ar4线斜率,具体见表1。
随着钢凝固过程的推进,钢液中溶质元素Mn、S在固液界面前沿富集。当溶质元素浓度超过钢液中MnS平衡溶解度时,将以[Mn]+[S]=(MnS)反应析出MnS,其标准吉布斯自由能变ΔGΘ(J/mol)采用下式计算:
液相中:ΔGΘ=-165248.81+90.90T (6)
δ相中:ΔGΘ=-202598.52+81.77T (7)
γ相中:ΔGΘ=-176892.30+57.95T (8)
此时,固液界面前沿剩余溶质元素浓度为平衡浓度:
式中:fMn和fS分别为钢液中溶质元素Mn,S的活度系数,可由下式计算所得:
式中:ei i,ei j为活度相互作用系数,如表2所示。
表1元素的平衡分配系数和扩散系数
注:R=1.987cal/mol·K,T为开尔文温度。
表2 1873K时活度相互作用系数
钢导热系数与温度T(℃),碳含量wC(%)以及相分率之间的如下关系式:
k=kαfα+kδfδ+kγfγ+klfl (12)
式中:fα,fδ,fγ,fl分别为α相,δ相,γ相和l相分别所占的分率。
kγ=21.6-8.35·10-3T (15)
kl=39.0 (16)
a1=0.425-4.385·10-4T (17)
a2=0.209+1.09·10-3T (18)
钢凝固过程中比热与钢中相组成之间满足如下关系式:
式中:
钢密度采用如下公式计算:
ρ=ργfγδfδ+ρ1lfl (24)
ρl=7100-73wC-(0.8-0.09wC)(T-1550) (27)
式中:ρ为钢的密度,kg/m3;fγ,fδ和fl分别为钢中γ-奥氏体相、δ-铁素体相和液相所占的分率;T为温度,℃;wC为C含量,%。
(2)宏观凝固传热计算
本实施例中为节约计算时间,选择铸坯横截面的1/4作为研究对象,如图3所示,钢连铸过程的二维凝固传热控制方程如下:
式中:T为温度,℃;ρ为密度,kg/m3;cp为热容,J/(kg·℃);keff为导热系数,W/(m·℃);L为凝固潜热,J/kg;fs为固相分率;t为时间,s;x和y分别为铸坯横截面宽面和窄面方向长度,m。
为了计算整个铸流铸坯温度分布,将整个铸流划分为连续单元切片,切片产生于连铸结晶器弯月面处,整个切片初始温度与浇铸温度相等。随后单元切片以拉速相同的运动速度从结晶器向下运动,直到消失于连铸坯火焰切割处。单元切片温度采用有效体积法并结合不同位置处边界条件所确定:
(a)结晶器
在结晶器冷却区,结晶器传热采用Davies等提出的经验公式
式中:q为热流密度,W/m2;z为离弯月面的距离,m;vc为拉速,m/s。
(b)二冷区
连铸二冷区传热边界条件较为复杂,二冷区铸坯表面热量传递包括:铸坯表面辐射散热、铸坯表面与冷却水雾间的强制对流换热、冷却水加热蒸发传热、铸坯与支撑辊接触导热等。为简化处理,将二冷区冷却水加热蒸发传热、铸坯与支撑辊接触导热等传热方式带走的热量,采用增大对流系数的方式来考虑。因此连铸二冷区的传热可简化表示为对流传热和辐射传热之和。
q=h(Tsurf-Tamb)+σε[(Tsurf+273.15)4-(Tamb+273.15)4] (30)
式中:h为综合对流换热系数,W/(m2·℃);Tsurf为铸坯表面温度,℃;Tamb为环境温度,℃;σ为Stefan-Boltzman常数5.67×10-8W/(m2·K4);ε为铸坯表面黑度系数,取0.8。
(c)空冷区
在空冷区内,铸坯表面主要以辐射的方式向外散热,辐射换热的热流公式如下:
q=σε[(Tsurf+273.15)4-(Tamb+273.15)4] (31)
结合单元切片所处位置和经历的时间关系,采用显示差分法求得单元切片各节点温度,从而求得铸坯不同位置处冷却速率CRi如下:
式中:Tl为液相线温度,℃;Ts为固相线温度,℃;△t为凝固时间,s。
(3)凝固组织生长计算
要完整地描述枝晶的生长行为,必须综合考虑热扩散、溶质扩散、流动、界面能以及高速生长时的动力学效应。由于金属合金的热扩散能力是溶质扩散能力的103~104倍,金属合金凝固过程微观组织演变动力学主要受溶质传输控制。在不考虑液相中强制对流和自然对流的情况下,固液两相中的溶质扩散控制方程可以用如下方程表示:
式中:t为时间,Dl和Ds分别为固液两相中溶质扩散系数,fs为固相分率,k为溶质平衡分配系数。
枝晶尖端过冷度ΔT可用如下公式表示:
ΔT=ΔTc+ΔTt+ΔTk+ΔTr (35)
式中:ΔTc、ΔTt、ΔTk和ΔTr分别为溶质过冷、热过冷、动力学过冷和曲率过冷。
由于动力学过冷只在非常高的凝固速率(快速凝固)下才对界面平衡温度产生影响,所以钢连铸凝固过程可以忽略动力学过冷的影响,只考虑热过冷、溶质过冷和曲率过冷对枝晶尖端过冷度的影响。因此,枝晶尖端过冷度如下公式表示:
式中:C0为合金初始浓度,Cl *为固液界面处液相溶质浓度,m0为液相线斜率,Г为Gibbs-Thomson系数,κ为平均固液界面曲率,为界面能各向异性函数,θ为枝晶择优生长角度,为固液界面法相与x轴的夹角。
固液界面曲率由界面处固相率梯度确定,可如下公式计算:
固液界面张力的各向异性可由如下公式计算:
式中:(fs)x和(fs)y为固相率的一阶偏导,(fs)xx,(fs)xy和(fs)yy为固相率的二阶偏导,ε为界面各向异性参数。
枝晶尖端生长速度受界面处溶质浓度控制,可由如下公式计算:
式中:vn为界面法相移动速度,k为溶质平衡分配系数,Ds为溶质元素在固相中扩散系数,Dl为溶质元素在液相中扩散系数。
(4)结果输出
主要是根据浇铸温度、拉速、钢种成分、冷却水量等连铸工艺条件,采用宏微观多尺度数学模型,预测连铸过程钢液凝固过程坯壳生长规律,并图像化地显示连铸坯凝固组织演变规律,以及定量化地表征连铸坯不同位置处坯壳厚度、晶粒尺寸、CET转变位置、等轴晶率等信息,为连铸工艺优化提供定量化的理论支撑。
从图4中可以直接看出,本实施例的预测方法可包括:
步骤一,收集连铸设备与工艺参数,设备参数主要包括:结晶器长度、有效高度、结晶器水槽尺寸与分布、铜板厚度、二冷区长度、二冷喷嘴布置等;工艺参数主要包括:钢种成分、浇铸温度、拉速、结晶器水量、温升、二冷各区水量、水流密度等;
步骤二,根据浇铸钢种成分和实测枝晶间距,采用溶质微观偏析模型,获得钢热物性参数数据。
应说明的是,本步骤中的实测枝晶间距是实现本发明之前利用上述的连铸设备制备在铸坯后实测该铸坯的枝晶间距。
步骤三,将实际连铸工艺参数和钢热物性参数计算数据作为宏微观多尺度数学模型的输入值,采用宏微观多尺度数学模型,预测不同连铸工艺条件连铸坯凝固组织形貌,如定量化地表征晶粒尺寸、定量化地表征柱状晶向等轴晶(CET)转变位置、定量化地表征等轴晶率和定量化地表征铸坯不同区域的坯壳厚度。
实验证明
下面以具体实施案例,针对某钢厂连铸生产重轨钢U75V连铸坯,方坯断面尺寸为1830×230mm,其结晶器有效高度为850mm。方坯连铸机总长度为20m,其二冷区包括5个喷水冷却区(总长度10.36m)和1个空冷区(长度9.64m),见表2。浇注过热度为25℃,拉速为0.68m/min。重轨钢U75V成分如表3所示。
表2方坯铸机各区冷却长度和水量
表3重轨钢U75V钢主要成分,wt.%
具体过程可描述如下:
第一步、收集该钢种连铸工艺与设备条件,为预测连铸工况条件下连铸坯凝固组织做准备,具体收集材料见表2和表3所示。
第二步、将重轨钢U75V和冷却速率作为溶质微观偏析模型的输入值,计算重轨钢U75V钢种热物性参数。
第三步、将重轨钢U75V浇铸工艺和设备条件,以及钢种热物性参数计算值作为宏微观多尺度数学模型输入值,计算重轨钢U75V连铸凝固组织生长。
第四步、图5(a)和图5(b)为实测重轨钢U75V连铸坯凝固组织与预测组织的比较,可以看出本发明提出的一种钢连铸坯凝固组织预测方法能够较好的预测连铸坯凝固组织,能够为优化连铸结晶器工艺条件提供直观可靠的信息。
最后应说明的是:以上所述的各实施例仅用于说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或全部技术特征进行等同替换;而这些修改或替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (8)

1.一种预测连铸生产的铸坯微观组织的方法,其特征在于,包括:
收集使用待分析钢种生产铸坯的连铸设备的设备参数和工艺参数;
针对待分析的钢种参数,确定该钢种的热物性参数;
利用预先建立的宏微观多尺度数学模型结合所述设备参数、工艺参数和所述钢热物性参数,预测当前连铸设备对应的铸坯的组织信息。
2.根据权利要求1所述的方法,其特征在于,针对待分析的钢种参数,确定该钢种的热物性参数,包括:
依据枝晶间溶质守恒原理和枝晶间溶质扩散原理建立枝晶间溶质微观偏析数学模型;
根据待分析的钢种成分和冷却强度,通过所述枝晶间溶质微观偏析数学模型,获取钢液凝固过程中枝晶间溶质偏析和相变信息;
根据所述枝晶间溶质偏析和相变信息,获取钢种的热物性参数;
所述热物性参数包括下述的一种或多种:液相线温度、固相线温度、导热系数、密度、比热和凝固潜热数据。
3.根据权利要求2所述的方法,其特征在于,根据所述枝晶间溶质偏析和相变信息,获取钢种的热物性参数,包括:
根据公式一,确定钢液的液相线温度Tl
公式一:
根据公式二,确定钢液凝固过程中δ/γ相转变开始温度TAr4
公式二:
根据公式三,确定钢液的导热系数k;
公式三:k=kαfα+kδfδ+kγfγ+klfl
根据公式四,确定钢液的钢密度ρ;
公式四:ρ=ργfγδfδlfl
根据公式五,确定钢液的凝固潜热数据c;
公式五:c=cαfα+cδfδ+cγfγ+clfl
其中,Cδ l,i为溶质元素i在钢液的δ/l相界面处液相中的浓度;mi和ni分别为伪二元Fe-i相图中液相线和Ar4线斜率;C0 l,i为元素i在钢液中的初始浓度,%;ki δ/l为溶质元素i在δ/l界面平衡分配系数;
kα,kδ,kγ,kl分别为α相,δ相,γ相和l相导热系数;fα,fδ,fγ,fl分别为α相,δ相,γ相和l相所占的分率;cα,cδ,cγ,cl分别为α相,δ相,γ相和l相比热容;ρα,ρδ,ργ,ρl分别为α相,δ相,γ相和l相密度。
4.根据权利要求1或3所述的方法,其特征在于,所述工艺参数包括下述的一种或多种:
浇铸温度、拉速、钢种成分、结晶器水量、温升和二冷各区水量;
和/或,
所述设备参数包括下述的一种或多种:
结晶器长度、结晶器的有效高度、结晶器水槽尺寸与分布、铜板厚度、二冷区长度、二冷喷嘴布置;
和/或,
所述铸坯的组织信息包括下述的一种或多种:
定量化地表征晶粒尺寸、定量化地表征柱状晶向等轴晶转变位置、定量化地表征等轴晶率和定量化地表征铸坯不同区域的坯壳厚度。
5.根据权利要求4所述的方法,其特征在于,利用预先建立的宏微观多尺度数学模型结合所述设备参数、工艺参数和所述钢热物性参数,预测当前连铸设备对应的铸坯的组织信息的步骤之前,所述方法还包括:
依据钢连铸过程的传热原理及生产铸坯的各连铸设备的传热信息,建立宏微观多尺度数学模型。
6.根据权利要求5所述的方法,其特征在于,所述宏微观多尺度数学模型包括:
依据钢连铸过程的传热原理及生产铸坯的各连铸设备的传热信息建立的宏观传热模型;
在预测过程中使用宏观传热模型的预测结果对铸坯组织信息预测的微观枝晶生长模型。
7.根据权利要求6所述的方法,其特征在于,所述宏观传热模型包括:
定义钢连铸过程的二维凝固传热控制方式为:
<mrow> <msub> <mi>&amp;rho;c</mi> <mi>p</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>T</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mrow> <mi>e</mi> <mi>f</mi> <mi>f</mi> </mrow> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>T</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>k</mi> <mrow> <mi>e</mi> <mi>f</mi> <mi>f</mi> </mrow> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>T</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;rho;</mi> <mi>L</mi> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>;</mo> </mrow>
式中:T为温度,℃;ρ为密度,kg/m3;cp为热容,J/(kg·℃);keff为导热系数,W/(m·℃);L为凝固潜热,J/kg;fs为固相分率;t为时间,s;x和y分别为铸坯横截面宽面和窄面方向长度,m;
定义结晶器冷却区的结晶器传热为:
<mrow> <mi>q</mi> <mo>=</mo> <mo>{</mo> <mn>2.64</mn> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>z</mi> <msub> <mi>v</mi> <mi>c</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <mn>0.91</mn> <mi>exp</mi> <mo>{</mo> <mo>-</mo> <mrow> <mo>(</mo> <mfrac> <mi>z</mi> <mrow> <mn>11.3</mn> <msub> <mi>v</mi> <mi>c</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>}</mo> <mo>+</mo> <mn>0.93</mn> <mo>}</mo> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mn>6</mn> </msup> <mo>;</mo> </mrow>
式中:q为热流密度,W/m2;z为离弯月面的距离,m;vc为拉速,m/s;
定义连铸二冷区的传热表示对流传热和辐射传热之和,具体为:
q=h(Tsurf-Tamb)+σε[(Tsurf+273.15)4-(Tamb+273.15)4];
式中:h为综合对流换热系数,W/(m2·℃);Tsurf为铸坯表面温度,℃;Tamb为环境温度,℃;σ为Stefan-Boltzman常数5.67×10-8W/(m2·K4);ε为铸坯表面黑度系数;
定义空冷区辐射换热为:
q=σε[(Tsurf+273.15)4-(Tamb+273.15)4];
定义铸坯不同区域/位置处冷却速率CRi为:
<mrow> <msub> <mi>CR</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mi>l</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>s</mi> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> </mrow>
式中:Tl为液相线温度,℃;Ts为固相线温度,℃;△t为凝固时间。
8.根据权利要求6所述的方法,其特征在于,所述微观枝晶生长模型包括:
固液两相中的溶质扩散控制过程表示为:
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>l</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msub> <mi>D</mi> <mi>l</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>l</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>+</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msub> <mi>D</mi> <mi>l</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>l</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>+</mo> <msub> <mi>C</mi> <mi>l</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>k</mi> <mo>)</mo> </mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>;</mo> </mrow>
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msub> <mi>D</mi> <mi>s</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>+</mo> <mfrac> <mo>&amp;part;</mo> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <msub> <mi>D</mi> <mi>s</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>&amp;rsqb;</mo> <mo>;</mo> </mrow>
式中:t为时间,Dl和Ds分别为固液两相中溶质扩散系数,fs为固相分率,k为溶质平衡分配系数;
枝晶尖端过冷度ΔT表示为:
式中:C0为合金初始浓度,为固液界面处液相溶质浓度,m0为液相线斜率,Г为Gibbs-Thomson系数,κ为固液界面曲率,为界面能各向异性函数,θ为枝晶择优生长角度,为固液界面法相与x轴的夹角;
固液界面曲率κ由固液界面处固相率梯度表示:
<mrow> <mi>&amp;kappa;</mi> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> </msub> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>y</mi> </msub> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mn>2</mn> </msubsup> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>y</mi> <mn>2</mn> </msubsup> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>x</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mi>y</mi> <mn>2</mn> </msubsup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>;</mo> </mrow>
固液相界面张力的各向异性表示为:
枝晶尖端生产速度受固液相界面处溶质浓度控制,表示为:
<mrow> <msub> <mi>v</mi> <mi>n</mi> </msub> <msubsup> <mi>C</mi> <mn>1</mn> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>D</mi> <mi>s</mi> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mi>s</mi> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>n</mi> </mrow> </mfrac> <msup> <mo>|</mo> <mo>*</mo> </msup> <mo>-</mo> <msub> <mi>D</mi> <mn>1</mn> </msub> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>C</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>&amp;part;</mo> <mi>n</mi> </mrow> </mfrac> <msup> <mo>|</mo> <mo>*</mo> </msup> <mo>;</mo> </mrow>
式中:(fs)x和(fs)y为固相率的一阶偏导,(fs)xx,(fs)xy和(fs)yy为固相率的二阶偏导,ε为界面各向异性参数;
vn为界面法相移动速度,k为溶质平衡分配系数,Ds为溶质元素在固相中扩散系数,Dl为溶质元素在液相中扩散系数。
CN201710459388.5A 2017-06-16 2017-06-16 一种预测连铸生产的铸坯微观组织的方法 Pending CN107229803A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710459388.5A CN107229803A (zh) 2017-06-16 2017-06-16 一种预测连铸生产的铸坯微观组织的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710459388.5A CN107229803A (zh) 2017-06-16 2017-06-16 一种预测连铸生产的铸坯微观组织的方法

Publications (1)

Publication Number Publication Date
CN107229803A true CN107229803A (zh) 2017-10-03

Family

ID=59936136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710459388.5A Pending CN107229803A (zh) 2017-06-16 2017-06-16 一种预测连铸生产的铸坯微观组织的方法

Country Status (1)

Country Link
CN (1) CN107229803A (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108038342A (zh) * 2018-01-24 2018-05-15 东北大学 一种预测包晶钢凝固过程相变的相场模拟方法
CN108213369A (zh) * 2018-01-05 2018-06-29 东北大学 一种降低重轨钢中a类夹杂物评级的连铸坯凝固组织控制方法
CN108296460A (zh) * 2018-04-24 2018-07-20 东北大学 一种利用增材稀释消除连铸坯中心线偏析的方法
CN108897972A (zh) * 2018-07-20 2018-11-27 辽宁石油化工大学 一种电渣重熔铸锭凝固微观组织的预测方法
CN109187610A (zh) * 2018-07-20 2019-01-11 中冶连铸技术工程有限责任公司 铸坯微观组织模拟方法
CN110008631A (zh) * 2019-04-18 2019-07-12 燕山大学 一种铜管轧制和冷却过程中参数模拟方法及其验证方法
CN110008561A (zh) * 2019-03-28 2019-07-12 东北大学 一种钢液中颗粒移动的预测方法
CN110263418A (zh) * 2019-06-17 2019-09-20 哈尔滨理工大学 一种体心立方合金微观偏析数值预测方法
CN110625079A (zh) * 2019-10-21 2019-12-31 北京科技大学 一种智能型连铸电磁搅拌在线控制系统及方法
CN111027172A (zh) * 2019-10-31 2020-04-17 东北大学 一种钢液对流情况下枝晶生长的预测方法
CN112185488A (zh) * 2020-09-30 2021-01-05 中南大学 基于数据驱动多组元镍基高温合金γ`相演化的预测方法
CN115007820A (zh) * 2022-07-28 2022-09-06 西北工业大学 一种消减球状石墨组织壁厚敏感性的超声连铸方法及系统
CN115229150A (zh) * 2022-07-04 2022-10-25 攀钢集团攀枝花钢铁研究院有限公司 钢轨夹杂物的控制方法
CN116187103A (zh) * 2023-04-26 2023-05-30 北京科技大学 在线预测连铸凝固终点位置和溶质元素分布的方法及系统
CN116384159A (zh) * 2023-05-29 2023-07-04 北京科技大学 一种连铸工艺温度仿真和宏观组织预测的方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122324A (zh) * 2011-03-03 2011-07-13 宣化钢铁集团有限责任公司 热轧棒线材力学性能预报方法
CN104384469A (zh) * 2014-12-16 2015-03-04 东北大学 一种钢连铸结晶器内初凝坯壳厚度的预测系统及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102122324A (zh) * 2011-03-03 2011-07-13 宣化钢铁集团有限责任公司 热轧棒线材力学性能预报方法
CN104384469A (zh) * 2014-12-16 2015-03-04 东北大学 一种钢连铸结晶器内初凝坯壳厚度的预测系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗森: "连铸坯凝固过程微观偏析与组织模拟及轻压下理论研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108213369B (zh) * 2018-01-05 2020-07-17 东北大学 一种降低重轨钢中a类夹杂物评级的连铸坯凝固组织控制方法
CN108213369A (zh) * 2018-01-05 2018-06-29 东北大学 一种降低重轨钢中a类夹杂物评级的连铸坯凝固组织控制方法
CN108038342B (zh) * 2018-01-24 2020-01-24 东北大学 一种预测包晶钢凝固过程相变的相场模拟方法
CN108038342A (zh) * 2018-01-24 2018-05-15 东北大学 一种预测包晶钢凝固过程相变的相场模拟方法
CN108296460A (zh) * 2018-04-24 2018-07-20 东北大学 一种利用增材稀释消除连铸坯中心线偏析的方法
CN108897972A (zh) * 2018-07-20 2018-11-27 辽宁石油化工大学 一种电渣重熔铸锭凝固微观组织的预测方法
CN109187610A (zh) * 2018-07-20 2019-01-11 中冶连铸技术工程有限责任公司 铸坯微观组织模拟方法
CN110008561A (zh) * 2019-03-28 2019-07-12 东北大学 一种钢液中颗粒移动的预测方法
CN110008631A (zh) * 2019-04-18 2019-07-12 燕山大学 一种铜管轧制和冷却过程中参数模拟方法及其验证方法
CN110263418A (zh) * 2019-06-17 2019-09-20 哈尔滨理工大学 一种体心立方合金微观偏析数值预测方法
CN110263418B (zh) * 2019-06-17 2022-10-21 哈尔滨理工大学 一种体心立方合金微观偏析数值预测方法
CN110625079B (zh) * 2019-10-21 2020-12-15 北京科技大学 一种智能型连铸电磁搅拌在线控制系统及方法
CN110625079A (zh) * 2019-10-21 2019-12-31 北京科技大学 一种智能型连铸电磁搅拌在线控制系统及方法
CN111027172A (zh) * 2019-10-31 2020-04-17 东北大学 一种钢液对流情况下枝晶生长的预测方法
CN112185488A (zh) * 2020-09-30 2021-01-05 中南大学 基于数据驱动多组元镍基高温合金γ`相演化的预测方法
CN112185488B (zh) * 2020-09-30 2023-08-18 中南大学 基于数据驱动多组元镍基高温合金γ`相演化的预测方法
CN115229150A (zh) * 2022-07-04 2022-10-25 攀钢集团攀枝花钢铁研究院有限公司 钢轨夹杂物的控制方法
CN115229150B (zh) * 2022-07-04 2024-05-14 攀钢集团攀枝花钢铁研究院有限公司 钢轨夹杂物的控制方法
CN115007820A (zh) * 2022-07-28 2022-09-06 西北工业大学 一种消减球状石墨组织壁厚敏感性的超声连铸方法及系统
CN115007820B (zh) * 2022-07-28 2023-03-21 西北工业大学 一种消减球状石墨组织壁厚敏感性的超声连铸方法及系统
CN116187103A (zh) * 2023-04-26 2023-05-30 北京科技大学 在线预测连铸凝固终点位置和溶质元素分布的方法及系统
CN116187103B (zh) * 2023-04-26 2023-06-20 北京科技大学 在线预测连铸凝固终点位置和溶质元素分布的方法及系统
CN116384159A (zh) * 2023-05-29 2023-07-04 北京科技大学 一种连铸工艺温度仿真和宏观组织预测的方法及系统
CN116384159B (zh) * 2023-05-29 2023-08-22 北京科技大学 一种连铸工艺温度仿真和宏观组织预测的方法及系统

Similar Documents

Publication Publication Date Title
CN107229803A (zh) 一种预测连铸生产的铸坯微观组织的方法
Louhenkilpi Continuous casting of steel
Beckermann Modelling of macrosegregation: applications and future needs
CN104384469B (zh) 一种钢连铸结晶器内初凝坯壳厚度的预测系统及方法
Pickering Macrosegregation in steel ingots: the applicability of modelling and characterisation techniques
CN108038342A (zh) 一种预测包晶钢凝固过程相变的相场模拟方法
Liu et al. Physical and numerical simulation of mixed columnar-equiaxed solidification during cold strip feeding in continuous casting
CN109785907A (zh) 一种钢液凝固过程中TiN夹杂物析出情况的预测方法
El-Bealy On the formation of macrosegregation and interdendritic cracks during dendritic solidification of continuous casting of steel
CN107052292A (zh) 一种基于热物性参数分布计算的连铸坯热跟踪计算方法
Li et al. Modeling of species transport and macrosegregation in heavy steel ingots
CN104540617B (zh) 凝固结束位置控制方法以及凝固结束位置控制装置
Shen et al. Dendrite growth behavior in directionally solidified Fe–C–Mn–Al alloys
Canté et al. Experimental analysis of the columnar-to-equiaxed transition in directionally solidified Al–Ni and Al–Sn alloys
CN110765599A (zh) 一种钢液凝固过程AlN夹杂物析出情况的预测方法
CN113591242B (zh) 一种连铸坯坯壳厚度的动态软测量方法
Wang et al. Mathematical modeling of thermal resistances of mold flux and air gap in continuous casting mold based on an inverse problem
Duan et al. A numerical study of the effect of multiple pouring on macrosegregation in a 438-ton steel ingot
Ramírez-López et al. Simulation of heat transfer in steel billets during continuous casting
CN110910967A (zh) 一种钢液凝固过程中CrN析出情况的预测方法
Brito et al. Cellular growth during the transient directional solidification of Zn-rich Zn–Cu monophasic and peritectic alloys
Šarler et al. A multiscale slice model for continuous casting of steel
Bai et al. Modeling effect of cooling conditions on solidification process during thermal cycle of rollers in twin-roll strip casting
Hu et al. Effect of different parameters on solidification structure of multi-crystalline silicon produced by continuous casting
Ma et al. Numerical simulation on the centre macrosegregation during Billet continuous casting with the soft reduction

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20171003

RJ01 Rejection of invention patent application after publication