CN108984918A - 一种电渣重熔自耗电极熔化速率的预测方法 - Google Patents

一种电渣重熔自耗电极熔化速率的预测方法 Download PDF

Info

Publication number
CN108984918A
CN108984918A CN201810805418.8A CN201810805418A CN108984918A CN 108984918 A CN108984918 A CN 108984918A CN 201810805418 A CN201810805418 A CN 201810805418A CN 108984918 A CN108984918 A CN 108984918A
Authority
CN
China
Prior art keywords
electroslag remelting
slag
electrode
heat
equation
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.)
Granted
Application number
CN201810805418.8A
Other languages
English (en)
Other versions
CN108984918B (zh
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.)
Liaoning Shihua University
Original Assignee
Liaoning Shihua University
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 Liaoning Shihua University filed Critical Liaoning Shihua University
Priority to CN201810805418.8A priority Critical patent/CN108984918B/zh
Publication of CN108984918A publication Critical patent/CN108984918A/zh
Application granted granted Critical
Publication of CN108984918B publication Critical patent/CN108984918B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • CCHEMISTRY; METALLURGY
    • C22METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22BPRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
    • C22B9/00General processes of refining or remelting of metals; Apparatus for electroslag or arc remelting of metals
    • C22B9/16Remelting metals
    • C22B9/18Electroslag remelting
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00Technologies related to metal processing
    • Y02P10/25Process efficiency

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Manufacturing & Machinery (AREA)
  • Computer Hardware Design (AREA)
  • Materials Engineering (AREA)
  • Mechanical Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Manufacture And Refinement Of Metals (AREA)

Abstract

本发明提供一种电渣重熔自耗电极熔化速率的预测方法,涉及冶金熔炼技术领域。该方法首先收集电渣重熔设备的设备参数和工艺参数,并采用ANSYS的EMAG电磁模块对电渣重熔体系内的电极、渣池、铸锭和周围空气进行电磁场分析,确定电渣重熔体系内电磁力和焦耳热分布;然后采用ANSYS FLUENT软件确定电渣重熔体系内温度场和流场分布;最后根据电渣重熔过程电极端部温度分布,确定自耗电极熔化速度。本发明提供的电渣重熔自耗电极熔化速率的预测方法,能够定量预测电渣重熔过程电极熔化速率,避免了实验测定时操作困难、精度较差、成本较高的缺点,为优化连铸工艺,控制电渣重熔过程熔化速率提供了理论指导。

Description

一种电渣重熔自耗电极熔化速率的预测方法
技术领域
本发明涉及冶金熔炼技术领域,尤其涉及一种电渣重熔自耗电极熔化速率的预测方法。
背景技术
电渣重熔(Electroslag Remelting,简称ESR)是一种具精炼和定向凝固于一体的特殊冶炼方法。电流通过自耗电极进入渣池的过程中,由于在整个供电回路中熔渣电导率远大于自耗电极、铸锭和外部短网,渣池将占据变压器二次电压的二次压降,在渣池中产生大量的焦耳热,使得插入渣池内的自耗电极温度升高至熔化温度,电极端部逐渐熔化并滴落,金属熔滴穿越渣池汇聚在金属熔池,在铜质水冷结晶器的冷却作用下由下到上逐渐凝固,形成自下而上的定向凝固组织。此外,电渣重熔过程电磁力使得渣池内湍流明显,渣金混合充分,金属熔滴形成和穿越渣池的过程中与熔渣之间发生剧烈的化学反应,有效去除金属熔体内有害元素,从而达到净化金属熔体的作用。由此可见,电渣重熔作为一种具金属熔化、精炼和凝固成型于一体的特种熔炼技术,能够有效地去除有害杂质元素和非金属夹杂物、提高金属纯净度和改善凝固结晶组织,获得高洁净度、高均匀性、超细组织和高性能的优质重熔铸锭,从而为船舶、核电、军工等特殊行业生产提供优质毛坯铸件。
在电渣重熔生产过程中,一般把金属熔池形貌控制作为控制电渣重熔铸锭凝固组织和质量的关键参数。这是因为熔池形状对凝固质量有着重要的影响,主要表现为枝晶的生长方向,其生长方向垂直于等温面,这样液态金属的结晶沿着熔池底部曲面的法线方向生长,则熔池的形状和深度直接决定着结晶方向。同时,金属熔池内枝晶生长也会影响溶质元素的偏析。为了保证电渣重熔锭的高质量,金属熔池形貌浅平为佳。然而,由于电渣重熔体系为高温非透明,现有手段很难在线监测电渣重熔过程金属熔池形貌,并加以控制。研究表明,电渣重熔金属熔池形貌不仅受电磁重熔电流强度影响,还受熔渣特性、厚度、填充比、插入深度等多因素影响。这些电渣重熔过程工艺参数将影响电渣重熔过程自耗电极的熔化速度,从而决定电渣重熔过程金属熔池深度和形貌。为此,控制电渣重熔过程自耗电极熔速对于获得高质量的电渣重熔铸锭至关重要。
电渣重熔熔速控制传统上主要有恒功率控制和恒熔速控制两种。恒功率控制在熔炼过程中由于在熔炼过程中由于热条件的变化,锭结构发生变化,熔池深度随时间的推移而增加,由于电极缩短,电阻损失降低以及感应环闭合,感应损失降低,有可能使熔速增加,从而无法保证熔炼过程稳定熔速的要求。采用恒熔速控制需要熔速检测系统,根据电极称重(文献:邓鑫,姜周华,董艳武.东北大学学报(自然科学版),2011,32(1):94-97.)和电极位移检测(文献:Li W Z,Wang WY,Hu Y C,Chen Y X.Metallurgical and MaterialsTransactions B,43B(2):276-289.)控制,但由于电渣重熔过程检测精度较差,很难满足整个电渣重熔过程熔速控制要求。电渣重熔过程熔速影响最终铸锭凝固组织形貌,根据凝固组织形貌和二次枝晶臂间距可定量预测电渣重熔电极熔化速率,但需要对电渣重熔铸锭进行解剖分析,检测周期较长,成本较大,劳动强度大,难以适用于电渣重熔工业生产过程。
发明内容
针对现有技术的缺陷,本发明提供一种电渣重熔自耗电极熔化速率的预测方法,实现对电渣重熔过程中自耗电极熔化速率的预测。
一种电渣重熔自耗电极熔化速率的预测方法,包括以下步骤:
步骤1、收集电渣重熔设备的设备参数和工艺参数;
所述设备参数包括:熔渣直径、熔渣高度、电极直径、电极长度以及铸锭直径和长度;所述工艺参数包括:钢种成分、熔化速度、冷却水流量;
步骤2、采用ANSYS软件中的EMAG电磁模块对电渣重熔体系内的电极、渣池、铸锭和周围空气进行电磁场分析,确定电渣重熔体系内电磁力和焦耳热分布,具体方法为:
步骤2.1、计算电渣重熔体系内的电磁场;
交变电流能够在电渣重熔体系内产生大量的焦耳热和磁场,能够为电渣重熔提供足够的热量熔化自耗电极,和产生电磁力搅拌渣池和均匀化温度;由于电渣重熔体系内存在渣金两相,交变电流经过渣相和金属相时将产生不同的磁场分布,进而影响电渣重熔体系内宏观流动、传热和传质传输现象;所述电渣重熔体系内的宏观流动、传热和传质传输电磁现象由Maxwell方程组描述,对于磁准静态场,Maxwell方程组如下:
M·B=0 (3)
其中,E为电渣重熔体系内的电场强度,B为电渣重熔体系内的磁感应强度,H为电渣重熔体系内的磁场强度,J为电渣重熔体系内的电流密度,D为电渣重熔体系内电位移矢量,ρe为电渣重熔体系内电荷体密度;
为表征在电磁场作用下媒质的宏观特性,除以上方程外,还需描述材料电磁特性的本构方程:
D=ε0εrE (5)
B=μ0μrH (6)
J=σE (7)
其中,ε0为真空介电常数,εr为相对介电常数,μ0为真空磁导率,μr为相对磁导率,σ为媒质电导率;
根据Lorentz定律,电渣重熔体系内电流密度J和磁感应强度B满足以下关系式:
F=J×B (8)
其中,F为电磁力;
根据Joule定律,电渣重熔体系内的电场强度E和电流密度J满足以下关系式:
QJ=E·J (9)
其中,QJ为电磁搅拌作用区域的发热密度;
求解电磁特性的本构方程所需设定的边界条件如下:
①电极顶部和铸锭底部满足的边界条件如下公式所示:
其中,Hdz为电极顶部表面和铸锭底部表面的磁场强度,n为电极顶部表面和铸锭底部表面的法向向量方向;
②电极侧面满足的边界条件,如下公式所示:
其中,Hdc为电极侧面的磁场强度,Re为电极半径,I为流过电极的电流强度;
③渣池表面满足的边界条件,如下公式所示:
其中,Hzb为渣池表面磁场强度x,y为渣池表面上的点的横纵坐标值,且满足Rm为铸锭半径;
④渣池侧面满足的边界条件,如下公式所示:
其中,Hzc为渣池侧面磁场强度,Rm为铸锭半径;
⑤铸锭侧面满足的边界条件,如下公式所示:
其中,Hzdc为铸锭侧面磁场强度;
⑥渣金界面处满足的边界条件,如下公式所示:
其中,Hzj为渣金界面处磁场强度,σs为渣的电导率,σm为金属的电导率,n′为渣金界面法向向量方向;
步骤2.2、计算电渣重熔体系内的流场;
电渣重熔体系为渣金两相共存体系,电渣重熔过程中电磁力和热浮力会促使渣金两相流动,并且自耗电极熔化和熔滴滴落也会对电渣重熔体系内流动行为产生影响;同时电渣重熔体系内的渣池和液体金属流动也会影响电磁场和温度场分布,进而又影响熔体流动行为;电渣重熔体系内渣金两相流动行为通过连续性方程、动量守恒方程,湍流模型和VOF模型描述;
所述连续性方程为单位时间内流体微元体中质量的增加,等于同一时间间隔内流入该微元体的净质量,其张量形式如下公式所示:
其中,t为时间,ρ为流体密度,ui为流体流动速度,xi表示坐标轴方向;
所述动量守恒方程为微元体中流体的动量对时间的变化率等于外界作用在该微元体上各种力之和,其张量形式如下公式所示:
其中,F为外力源项,p为压力,μl为层流动力黏度,μt为湍流黏度,g为重力加速度;
所述湍流黏度μt如下公式所示:
其中,fμ=exp[-3.4/(1+Ret/50)2],k为湍动能,ε为湍动能耗散率,Cμ为经验常数;
所述湍流模型又称为低Re数k-ε双方程模型,包括湍动能k方程和湍动能耗散率ε方程;
所述湍动能k方程如下公式所示:
其中,σk为经验常数;
所述湍动能耗散率ε方程如下公式所示:
其中,f1=1.0,f2=1-0.3exp(-Ret 2),C1、C2均为经验常数;
所述多相流VOF模型为多种互不相溶流体共用一组动量方程,在全部计算区域内追踪流体体积分数;电渣重熔体系存在熔渣和钢液两相,对于不可压缩流体,钢液体积分数αst满足如下传输方程:
其中,U为电渣重熔体系内速度矢量,当αst=1时代表钢液,αst=0=0时为熔渣,αst在0~1之间为钢/渣界面;
所述流场满足的边界条件分别为:
①电极熔化端部以及速度入口处满足的边界条件如下公式所示:
其中,Udd为电极熔化端部以及速度入口处速度,Me为自耗电极重熔速度;
②渣池表面满足自由滑移边界条件,如下公式所示:
其中,uk为z轴方向流动速度,n″为渣池表面法向向量方向;
③渣池侧面满足无滑移边界条件,如下公式所示:
Uzc=0 (24)
其中,Uzc为渣池侧面速度;
④铸锭侧面满足无滑移边界条件,如下公式所示:
Udc=0 (25)
其中,Udc为铸锭侧面速度;
⑤铸锭底部满足充分流动出口边界条件,如下公式所示:
其中,Uzd为铸锭侧面速度,n″′为铸锭底部法向向量方向;
步骤3、采用ANSYS FLUENT软件的UDF接口通过自编程的方法将ANSYS EMAG对电渣重熔过程电磁场分析所得到的稳态电磁场结果文件,分别以动量源项和能量源项的形式导入体积电磁力和体积焦耳热,从而耦合考虑电渣重熔过程交变电流通过电渣重熔体系所产生的电磁力和焦耳热对电渣重熔体系内多相流流动和热量传输的影响,确定电渣重熔体系内温度场和流场分布;
电渣重熔体系由自耗电极、渣池、金属熔池和铸锭组成,各部分之间相互传热;电渣重熔过程热量传输由如下方程描述:
其中,H为热焓,Cp为等压比热容,λl和λt分别为层流导热系数和湍流导热系数;
温度场边界条件如下:
①电极熔化端部采用温度边界条件;
电极端部在渣池焦耳热作用下熔化,为了简化处理,将电极端部温度设置为金属液相线温度Tl,m,如下公式所示:
T1=Tl,m (28)
其中,T1为电极端部温度;
②渣池表面采用混合传热边界条件;
渣池表面换热比较复杂,不仅存在与大气的对流和辐射两种传热方式,还存在与结晶器和电极之间的辐射传热;为简化处理,只考虑渣池与大气之间的对流和辐射传热;因此,渣池表面与大气交换的总热量为:
Q=Qc+Qr=h(Ts-Tamb)A (29)
其中,Qc为对流换热热量,Qr为辐射换热热量,h为综合换热系数,Ts为渣池表面温度,Tamb为环境温度,A为辐射界面面积;
③渣池侧面采用温度边界条件;
渣池表面与结晶器接触处会形成一层薄的渣壳,为了简化处理,渣池侧面温度与熔渣固相线温度Ts,slag相同,如下公式所示:
T2=Ts,slag (30)
其中,T2为渣池侧面温度;
④铸锭侧面采用混合传热边界条件;
由于铸锭凝固收缩,使得铸锭之间会形成不规则的气隙分布;为了简化处理,将渣膜、气隙、结晶器铜板和冷却水之间的复杂传热,采用综合换热系数简化处理成对流传热;
⑤铸锭底部采用混合传热边界条件;
铸锭底部与冷却水箱之间的传热同样十分复杂,为了简化处理,将采用综合传热系数类比为对流传热;
步骤4、根据电渣重熔过程电极端部温度分布,确定自耗电极熔化速度;
电渣重熔过程,电流流经渣池所产生的热量,通过渣池传递到电极端部,加热熔化电极;根据渣/电极界面处热量守恒推导出了电极熔化速度与热量传递之间的如下关系式:
其中,me为自耗电极的熔化速度,Cp,e为电极比热,Tl,m为金属液相线温度,L为金属潜热,Re为电极半径,qse为渣/电极界面处渣池向电极传递的热通量,r为距离电极横截面中心的距离。
由上述技术方案可知,本发明的有益效果在于:本发明提供的一种电渣重熔自耗电极熔化速率的预测方法,全面考虑了电渣重熔体系内电磁场、焦耳热和流动的影响,建立了电渣重熔体系内多物理场耦合的宏观传输数学模型,能够定量预测电渣重熔过程电极熔化速率,该方法具有简单、快捷、易实现、结果精确地特点,避免了实验测定时操作困难、精度较差、成本较高的缺点,为优化连铸工艺,控制电渣重熔过程熔化速率提供了理论指导。
附图说明
图1为本发明实施例提供的一种电渣重熔自耗电极熔化速率的预测方法的流程图;
图2为本发明实施例提供的电渣重熔体系的示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例以某特殊钢厂电渣重熔生产转子钢30Cr1Mo1V过程为例,使用本发明的一种电渣重熔自耗电极熔化速率的预测方法,对电渣重溶体系内的自耗电极的熔化速率进行预测。
一种电渣重熔自耗电极熔化速率的预测方法,如图1所示,包括以下步骤:
步骤1、收集电渣重熔设备的设备参数和工艺参数;
设备参数包括:熔渣直径、熔渣高度、电极直径、电极长度以及铸锭直径和长度;工艺参数包括:钢种成分、熔化速度、冷却水流量。
本实施例中,电渣重熔的设备参数如表1所示,转子钢30Cr1Mo1V的钢种成分如表2所示。表1电渣重熔设备明细
表2转子钢30Cr1Mo1V成分(wt%)
C Si Mn S P Cr Mo V
0.3 0.25 0.8 0.01 0.01 1.2 1.2 0.25
步骤2、采用ANSYS软件中的EMAG电磁模块对如图2所示的电渣重熔体系内的电极、渣池、铸锭和周围空气进行电磁场分析,确定电渣重熔体系内电磁力和焦耳热分布,具体方法为:
步骤2.1、计算电渣重熔体系内的电磁场;
交变电流能够在电渣重熔体系内产生大量的焦耳热和磁场,能够为电渣重熔提供足够的热量熔化自耗电极,和产生电磁力搅拌渣池和均匀化温度;由于电渣重熔体系内存在渣金两相,交变电流经过渣相和金属相时将产生不同的磁场分布,进而影响电渣重熔体系内宏观流动、传热和传质传输现象;所述电渣重熔体系内的宏观流动、传热和传质传输电磁现象由Maxwell方程组描述,对于磁准静态场,Maxwell方程组如下:
其中,E为电渣重熔体系内的电场强度,B为电渣重熔体系内的磁感应强度,H为电渣重熔体系内的磁场强度,J为电渣重熔体系内的电流密度,D为电渣重熔体系内电位移矢量,ρe为电渣重熔体系内电荷体密度;
为表征在电磁场作用下媒质的宏观特性,除以上方程外,还需描述材料电磁特性的本构方程:
D=ε0εrE (5)
B=μ0μrH (6)
J=σE (7)
其中,ε0为真空介电常数,εr为相对介电常数,μ0为真空磁导率,μr为相对磁导率,σ为媒质电导率;
根据Lorentz定律,电渣重熔体系内电流密度J和磁感应强度B满足以下关系式:
F=J×B (8)
其中,F为电磁力;
根据Joule定律,电渣重熔体系内电场强度E和电流密度J满足以下关系式:
QJ=E·J (9)
其中,QJ为电磁搅拌作用区域的发热密度;
求解电磁特性的本构方程所需设定的边界条件如下:
①电极顶部和铸锭底部满足的边界条件如下公式所示:
其中,Hdz为电极顶部表面和铸锭底部表面的磁场强度,n为电极顶部表面和铸锭底部表面的法向向量方向;
②电极侧面满足的边界条件,如下公式所示:
其中,Hdc为电极侧面的磁场强度,Re为电极半径,I为流过电极的电流强度;
③渣池表面满足的边界条件,如下公式所示:
其中,Hzb为渣池表面磁场强度x,y为渣池表面上的点的横纵坐标值,且满足
④渣池侧面满足的边界条件,如下公式所示:
其中,Hzc为渣池侧面磁场强度,Rm为铸锭半径;
⑤铸锭侧面满足的边界条件,如下公式所示:
其中,Hzdc为铸锭侧面磁场强度;
⑥渣金界面处满足的边界条件,如下公式所示:
其中,Hzj为渣金界面处磁场强度,σs为渣的电导率,σm为金属的电导率,n′为渣金界面法向向量方向。
步骤2.2、计算电渣重熔体系内的流场;
电渣重熔体系为渣金两相共存体系,电渣重熔过程中电磁力和热浮力会促使渣金两相流动,并且自耗电极熔化和熔滴滴落也会对电渣重熔体系内流动行为产生影响;同时电渣重熔体系内的渣池和液体金属流动也会影响电磁场和温度场分布,进而又影响熔体流动行为;电渣重熔体系内渣金两相流动行为通过连续性方程、动量守恒方程,湍流模型和VOF模型描述;
所述连续性方程为单位时间内流体微元体中质量的增加,等于同一时间间隔内流入该微元体的净质量,其张量形式如下公式所示:
其中,t为时间,s;ρ为流体密度,ui为流体流动速度,xi表示坐标轴方向;
所述动量守恒方程为微元体中流体的动量对时间的变化率等于外界作用在该微元体上各种力之和,其张量形式如下公式所示:
其中,F为外力源项,p为压力,μl为层流动力黏度,μt为湍流黏度,g为重力加速度;
所述湍流黏度μt如下公式所示:
其中,fμ=exp[-3.4/(1+Ret/50)2],k为湍动能,ε为湍动能耗散率,Cμ为经验常数;
所述湍流模型又称为低Re数k-ε双方程模型,包括湍动能k方程和湍动能耗散率ε方程;
所述湍动能k方程如下公式所示:
其中,σk为经验常数;
所述湍动能耗散率ε方程如下公式所示:
其中,f1=1.0,f2=1-0.3exp(-Ret 2),C1、C2均为经验常数;
所述多相流VOF模型为多种互不相溶流体共用一组动量方程,在全部计算区域内追踪流体体积分数;电渣重熔体系存在熔渣和钢液两相,对于不可压缩流体,钢液体积分数αst满足如下传输方程:
其中,U为电渣重熔体系内速度矢量,当αst=1时代表钢液,αst=0=0时为熔渣,αst在0~1之间为钢/渣界面;
所述流场满足的边界条件分别为:
①电极熔化端部以及速度入口处满足的边界条件如下公式所示:
其中,Udd为电极熔化端部以及速度入口处速度,Me为自耗电极重熔速度;
②渣池表面满足自由滑移边界条件,如下公式所示:
其中,uk为z轴方向流动速度,n″为渣池表面法向向量方向;
③渣池侧面满足无滑移边界条件,如下公式所示:
Uzc=0 (24)
其中,Uzc为渣池侧面速度;
④铸锭侧面满足无滑移边界条件,如下公式所示:
Udc=0 (25)
其中,Udc为铸锭侧面速度;
⑤铸锭底部满足充分流动出口边界条件,如下公式所示:
其中,Uzd为铸锭侧面速度,n″′为铸锭底部法向向量方向。
步骤3、采用ANSYS FLUENT软件的UDF接口通过自编程的方法将ANSYS EMAG对电渣重熔过程电磁场分析所得到的稳态电磁场结果文件,分别以动量源项和能量源项的形式导入体积电磁力和体积焦耳热,从而耦合考虑电渣重熔过程交变电流通过电渣重熔体系所产生的电磁力和焦耳热对电渣重熔体系内多相流流动和热量传输的影响,确定电渣重熔体系内温度场和流场分布;
电渣重熔体系由自耗电极、渣池、金属熔池和铸锭组成,各部分之间相互传热;电渣重熔过程热量传输由如下方程描述:
其中,H为热焓,Cp为等压比热容,λl和λt分别为层流导热系数和湍流导热系数;
温度场边界条件如下:
①电极熔化端部采用温度边界条件;
电极端部在渣池焦耳热作用下熔化,为了简化处理,将电极端部温度设置为金属液相线温度Tl,m,如下公式所示:
T1=Tl,m (28)
其中,T1为电极端部温度;
②渣池表面采用混合传热边界条件;
渣池表面换热比较复杂,不仅存在与大气的对流和辐射两种传热方式,还存在与结晶器和电极之间的辐射传热;为简化处理,只考虑渣池与大气之间的对流和辐射传热;因此,渣池表面与大气交换的总热量为:
Q=Qc+Qr=h(Ts-Tamb)A (29)
其中,Qc为对流换热热量,Qr为辐射换热热量,h为综合换热系数,Ts为渣池表面温度,Tamb为环境温度,A为辐射界面面积;
③渣池侧面采用温度边界条件;
渣池表面与结晶器接触处会形成一层薄的渣壳,为了简化处理,渣池侧面温度与熔渣固相线温度Ts,slag相同,如下公式所示:
T2=Ts,slag (30)
其中,T2为渣池侧面温度;
④铸锭侧面采用混合传热边界条件;
由于铸锭凝固收缩,使得铸锭之间会形成不规则的气隙分布;为了简化处理,将渣膜、气隙、结晶器铜板和冷却水之间的复杂传热,采用综合换热系数简化处理成对流传热;
⑤铸锭底部采用混合传热边界条件;
铸锭底部与冷却水箱之间的传热同样十分复杂,为了简化处理,将采用综合传热系数类比为对流传热;
步骤4、根据电渣重熔过程电极端部温度分布,确定自耗电极熔化速度;
电渣重熔过程,电流流经渣池所产生的热量,通过渣池传递到电极端部,加热熔化电极;根据渣/电极界面处热量守恒推导出了电极熔化速度与热量传递之间的如下关系式:
其中,me为自耗电极的熔化速度,Cp,e为电极比热,Tl,m为金属液相线温度,L为金属潜热,Re为电极半径,qse为渣/电极界面处渣池向电极传递的热通量,r为距离电极横截面中心的距离。
本实施例中,电渣重熔过程电极熔速计算值与实测值比较如表3所示,从表中可以看出,本发明的电渣重熔自耗电极熔化速率的预测方法,能够较好的预测电渣重熔过程中自耗电极的熔化速率。
表3电渣重熔过程电极熔速计算值与实测值比较
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (6)

1.一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:包括以下步骤:
步骤1、收集电渣重熔设备的设备参数和工艺参数;
步骤2、采用ANSYS软件中的EMAG电磁模块对电渣重熔体系内的电极、渣池、铸锭和周围空气进行电磁场分析,确定电渣重熔体系内电磁力和焦耳热分布;
步骤2.1、计算电渣重熔体系内的电磁场;
步骤2.2、计算电渣重熔体系内的流场;
步骤3、采用ANSYS FLUENT软件的UDF接口,通过自编程的方法将ANSYS EMAG对电渣重熔过程电磁场分析所得到的稳态电磁场结果文件,分别以动量源项和能量源项的形式导入体积电磁力和体积焦耳热,从而耦合考虑电渣重熔过程交变电流通过电渣重熔体系所产生的电磁力和焦耳热对电渣重熔体系内多相流流动和热量传输的影响,确定电渣重熔体系内温度场和流场分布;
步骤4、根据电渣重熔过程电极端部温度分布,确定自耗电极熔化速度。
2.根据权利要求1所述的一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:步骤1所述设备参数包括:熔渣直径、熔渣高度、电极直径、电极长度以及铸锭直径和长度;所述工艺参数包括:钢种成分、熔化速度、冷却水流量。
3.根据权利要求1所述的一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:所述步骤2.1的具体方法为:
交变电流能够在电渣重熔体系内产生大量的焦耳热和磁场,能够为电渣重熔提供足够的热量熔化自耗电极,和产生电磁力搅拌渣池和均匀化温度;由于电渣重熔体系内存在渣金两相,交变电流经过渣相和金属相时将产生不同的磁场分布,进而影响电渣重熔体系内宏观流动、传热和传质传输现象;所述电渣重熔体系内的宏观流动、传热和传质传输电磁现象由Maxwell方程组描述,对于磁准静态场,Maxwell方程组如下:
其中,E为电渣重熔体系内的电场强度,B为电渣重熔体系内的磁感应强度,H为电渣重熔体系内的磁场强度,J为电渣重熔体系内的电流密度,D为电渣重熔体系内电位移矢量,ρe为电渣重熔体系内电荷体密度;
为表征在电磁场作用下媒质的宏观特性,除以上方程外,还需描述材料电磁特性的本构方程:
D=ε0εrE (5)
B=μ0μrH (6)
J=σE (7)
其中,ε0为真空介电常数,εr为相对介电常数,μ0为真空磁导率,μr为相对磁导率,σ为媒质电导率;
根据Lorentz定律,电渣重熔体系内电流密度J和磁感应强度B满足以下关系式:
F=J×B (8)
其中,F为电磁力;
根据Joule定律,电渣重熔体系内电场强度E和电流密度J满足以下关系式:
QJ=E·J (9)
其中,QJ为电磁搅拌作用区域的发热密度;
求解电磁特性的本构方程所需设定的边界条件如下:
①电极顶部和铸锭底部满足的边界条件如下公式所示:
其中,Hdz为电极顶部表面和铸锭底部表面的磁场强度,n为电极顶部表面和铸锭底部表面的法向向量方向;
②电极侧面满足的边界条件,如下公式所示:
其中,Hdc为电极侧面的磁场强度,Re为电极半径,I为流过电极的电流强度;
③渣池表面满足的边界条件,如下公式所示:
其中,Hzb为渣池表面磁场强度x,y为渣池表面上的点的横纵坐标值,且满足
④渣池侧面满足的边界条件,如下公式所示:
其中,Hzc为渣池侧面磁场强度,Rm为铸锭半径;
⑤铸锭侧面满足的边界条件,如下公式所示:
其中,Hzdc为铸锭侧面磁场强度;
⑥渣金界面处满足的边界条件,如下公式所示:
其中,Hzj为渣金界面处磁场强度,σs为渣的电导率,σm为金属的电导率,n′为渣金界面法向向量方向。
4.根据权利要求3所述的一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:所述步骤2.2的具体方法为:
电渣重熔体系内渣金两相流动行为通过连续性方程、动量守恒方程,湍流模型和VOF模型描述;
所述连续性方程为单位时间内流体微元体中质量的增加,等于同一时间间隔内流入该微元体的净质量,其张量形式如下公式所示:
其中,t为时间,ρ为流体密度,ui为流体流动速度,xi表示坐标轴方向;
所述动量守恒方程为微元体中流体的动量对时间的变化率等于外界作用在该微元体上各种力之和,其张量形式如下公式所示:
其中,F为外力源项,p为压力,μl为层流动力黏度,μt为湍流黏度,g为重力加速度;
所述湍流黏度μt如下公式所示:
其中,fμ=exp[-3.4/(1+Ret/50)2],k为湍动能,ε为湍动能耗散率,Cμ为经验常数;
所述湍流模型又称为低Re数k-ε双方程模型,包括湍动能k方程和湍动能耗散率ε方程;
所述湍动能k方程如下公式所示:
其中,σk为经验常数;
所述湍动能耗散率ε方程如下公式所示:
其中,f1=1.0,f2=1-0.3exp(-Ret 2),C1、C2均为经验常数;
所述多相流VOF模型为多种互不相溶流体共用一组动量方程,在全部计算区域内追踪流体体积分数;电渣重熔体系存在熔渣和钢液两相,对于不可压缩流体,钢液体积分数αst满足如下传输方程:
其中,U为电渣重熔体系内速度矢量,当αst=1时代表钢液,αst=0=0时为熔渣,αst在0~1之间为钢/渣界面;
所述流场满足的边界条件分别为:
①电极熔化端部以及速度入口处满足的边界条件如下公式所示:
其中,Udd为电极熔化端部以及速度入口处速度,Me为自耗电极重熔速度;
②渣池表面满足自由滑移边界条件,如下公式所示:
其中,uk为z轴方向流动速度,n″为渣池表面法向向量方向;
③渣池侧面满足无滑移边界条件,如下公式所示:
Uzc=0 (24)
其中,Uzc为渣池侧面速度;
④铸锭侧面满足无滑移边界条件,如下公式所示:
Udc=0 (25)
其中,Udc为铸锭侧面速度;
⑤铸锭底部满足充分流动出口边界条件,如下公式所示:
其中,Uzd为铸锭侧面速度,n″′为铸锭底部法向向量方向。
5.根据权利要求4所述的一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:所述步骤3的具体方法为:
电渣重熔体系由自耗电极、渣池、金属熔池和铸锭组成,各部分之间相互传热;电渣重熔过程热量传输由如下方程描述:
其中,H为热焓,Cp为等压比热容,λl和λt分别为层流导热系数和湍流导热系数;
温度场边界条件如下:
①电极熔化端部采用温度边界条件;
电极端部在渣池焦耳热作用下熔化,为了简化处理,将电极端部温度设置为金属液相线温度Tl,m,如下公式所示:
T1=Tl,m (28)
其中,T1为电极端部温度;
②渣池表面采用混合传热边界条件;
渣池表面换热比较复杂,不仅存在与大气的对流和辐射两种传热方式,还存在与结晶器和电极之间的辐射传热;为简化处理,只考虑渣池与大气之间的对流和辐射传热;因此,渣池表面与大气交换的总热量为:
Q=Qc+Qr=h(Ts-Tamb)A (29)
其中,Qc为对流换热热量,Qr为辐射换热热量,h为综合换热系数,Ts为渣池表面温度,Tamb为环境温度,A为辐射界面面积;
③渣池侧面采用温度边界条件;
渣池表面与结晶器接触处会形成一层薄的渣壳,为了简化处理,渣池侧面温度与熔渣固相线温度Ts,slag相同,如下公式所示:
T2=Ts,slag (30)
其中,T2为渣池侧面温度;
④铸锭侧面采用混合传热边界条件;
由于铸锭凝固收缩,使得铸锭之间会形成不规则的气隙分布;为了简化处理,将渣膜、气隙、结晶器铜板和冷却水之间的复杂传热,采用综合换热系数简化处理成对流传热;
⑤铸锭底部采用混合传热边界条件;
铸锭底部与冷却水箱之间的传热同样十分复杂,为了简化处理,将采用综合传热系数类比为对流传热。
6.根据权利要求5所述的一种电渣重熔自耗电极熔化速率的预测方法,其特征在于:所述步骤4的具体方法为:
电渣重熔过程,电流流经渣池所产生的热量,通过渣池传递到电极端部,加热熔化电极;根据渣/电极界面处热量守恒推导出了电极熔化速度与热量传递之间的如下关系式:
其中,me为自耗电极的熔化速度,Cp,e为电极比热,Tl,m为金属液相线温度,L为金属潜热,Re为电极半径,qse为渣/电极界面处渣池向电极传递的热通量,r为距离电极横截面中心的距离。
CN201810805418.8A 2018-07-20 2018-07-20 一种电渣重熔自耗电极熔化速率的预测方法 Active CN108984918B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810805418.8A CN108984918B (zh) 2018-07-20 2018-07-20 一种电渣重熔自耗电极熔化速率的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810805418.8A CN108984918B (zh) 2018-07-20 2018-07-20 一种电渣重熔自耗电极熔化速率的预测方法

Publications (2)

Publication Number Publication Date
CN108984918A true CN108984918A (zh) 2018-12-11
CN108984918B CN108984918B (zh) 2023-04-18

Family

ID=64549443

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810805418.8A Active CN108984918B (zh) 2018-07-20 2018-07-20 一种电渣重熔自耗电极熔化速率的预测方法

Country Status (1)

Country Link
CN (1) CN108984918B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113934139A (zh) * 2020-06-29 2022-01-14 宝武特种冶金有限公司 一种基于在线仿真模型的真空电弧重熔过程熔化速度控制方法
CN113987892A (zh) * 2021-12-28 2022-01-28 北京钢研高纳科技股份有限公司 一种控制高温合金偏析的真空电弧重熔3d模型及控制方法
WO2023142422A1 (zh) * 2022-01-28 2023-08-03 苏州大学 单进多出式水冷结晶器电渣重熔装置及电渣重熔的方法

Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1370721A (en) * 1971-03-17 1974-10-16 Inst Elektroswarki Patona Method of electroslag remelting of consumable electrodes
WO1980001809A1 (en) * 1979-02-26 1980-09-04 Inst Elektroswarki Patona Method and device for controlling the process of electroslag remelting of a plurality of consumable electrodes
US4450007A (en) * 1982-12-13 1984-05-22 Cabot Corporation Process for electroslag remelting of manganese-base alloys
GB9216573D0 (en) * 1991-12-10 1992-09-16 Leybold Durferrit Gmbh System for re-melting an electrode
US5686199A (en) * 1996-05-07 1997-11-11 Alliedsignal Inc. Flow field plate for use in a proton exchange membrane fuel cell
US6115404A (en) * 1999-02-03 2000-09-05 Sandia Corporation Dynamic control of remelting processes
TW449820B (en) * 1996-02-15 2001-08-11 Tokai Carbon Kk Plasma-etching electrode plate
US20050071039A1 (en) * 2003-09-30 2005-03-31 Tokyo Electron Limited System and method for using first-principles simulation to provide virtual sensors that facilitate a semiconductor manufacturing process
CN1751834A (zh) * 2005-10-20 2006-03-29 武汉理工大学 一种磁控熔化极焊接方法及其拓展应用和通用设备
CN101032788A (zh) * 2007-04-28 2007-09-12 罗键 一种电磁复合场熔化极堆焊方法、设备及其拓展应用
CN201793688U (zh) * 2010-09-27 2011-04-13 东北大学 一种降低氢、氧含量的电渣炉保护罩
CN201999975U (zh) * 2011-02-22 2011-10-05 沈阳东大材料先进制备技术工程研究中心 一种具有电极称重功能的横臂移动式电渣炉
CN102297669A (zh) * 2011-01-25 2011-12-28 东北大学 一种电渣炉熔化电极剩余长度的软测量方法
JP2013111505A (ja) * 2011-11-25 2013-06-10 Gastar Corp 加圧容器およびその加圧容器を用いた空気溶存装置
CN104330411A (zh) * 2014-11-10 2015-02-04 东北大学 一种动态观测模拟电渣重熔过程熔滴行为的装置及方法
CN104388689A (zh) * 2014-11-28 2015-03-04 中国科学院金属研究所 一种电磁复合控制电渣重熔细晶铸造方法及装置
CN104451183A (zh) * 2013-09-16 2015-03-25 西安扩力机电科技有限公司 一种真空自耗电弧炉熔速智能控制系统
CN104713370A (zh) * 2015-02-13 2015-06-17 中钢集团邢台机械轧辊有限公司 电渣炉自耗电极插入深度及熔速的控制方法
CN105364269A (zh) * 2015-12-08 2016-03-02 辽宁石油化工大学 一种铜板及其合金的熔化焊接方法
CN107301291A (zh) * 2017-06-21 2017-10-27 河北工业大学 一种基于电渣重熔过程的电磁场、温度场和流动场耦合有限元计算方法
CN108278901A (zh) * 2018-01-30 2018-07-13 上海荣丰科技发展有限公司 一种用于冶炼熔渣成分调整的工业炉及生产工艺

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1370721A (en) * 1971-03-17 1974-10-16 Inst Elektroswarki Patona Method of electroslag remelting of consumable electrodes
WO1980001809A1 (en) * 1979-02-26 1980-09-04 Inst Elektroswarki Patona Method and device for controlling the process of electroslag remelting of a plurality of consumable electrodes
US4450007A (en) * 1982-12-13 1984-05-22 Cabot Corporation Process for electroslag remelting of manganese-base alloys
GB9216573D0 (en) * 1991-12-10 1992-09-16 Leybold Durferrit Gmbh System for re-melting an electrode
TW449820B (en) * 1996-02-15 2001-08-11 Tokai Carbon Kk Plasma-etching electrode plate
US5686199A (en) * 1996-05-07 1997-11-11 Alliedsignal Inc. Flow field plate for use in a proton exchange membrane fuel cell
US6115404A (en) * 1999-02-03 2000-09-05 Sandia Corporation Dynamic control of remelting processes
US20050071039A1 (en) * 2003-09-30 2005-03-31 Tokyo Electron Limited System and method for using first-principles simulation to provide virtual sensors that facilitate a semiconductor manufacturing process
CN1751834A (zh) * 2005-10-20 2006-03-29 武汉理工大学 一种磁控熔化极焊接方法及其拓展应用和通用设备
CN101032788A (zh) * 2007-04-28 2007-09-12 罗键 一种电磁复合场熔化极堆焊方法、设备及其拓展应用
CN201793688U (zh) * 2010-09-27 2011-04-13 东北大学 一种降低氢、氧含量的电渣炉保护罩
CN102297669A (zh) * 2011-01-25 2011-12-28 东北大学 一种电渣炉熔化电极剩余长度的软测量方法
CN201999975U (zh) * 2011-02-22 2011-10-05 沈阳东大材料先进制备技术工程研究中心 一种具有电极称重功能的横臂移动式电渣炉
JP2013111505A (ja) * 2011-11-25 2013-06-10 Gastar Corp 加圧容器およびその加圧容器を用いた空気溶存装置
CN104451183A (zh) * 2013-09-16 2015-03-25 西安扩力机电科技有限公司 一种真空自耗电弧炉熔速智能控制系统
CN104330411A (zh) * 2014-11-10 2015-02-04 东北大学 一种动态观测模拟电渣重熔过程熔滴行为的装置及方法
CN104388689A (zh) * 2014-11-28 2015-03-04 中国科学院金属研究所 一种电磁复合控制电渣重熔细晶铸造方法及装置
CN104713370A (zh) * 2015-02-13 2015-06-17 中钢集团邢台机械轧辊有限公司 电渣炉自耗电极插入深度及熔速的控制方法
CN105364269A (zh) * 2015-12-08 2016-03-02 辽宁石油化工大学 一种铜板及其合金的熔化焊接方法
CN107301291A (zh) * 2017-06-21 2017-10-27 河北工业大学 一种基于电渣重熔过程的电磁场、温度场和流动场耦合有限元计算方法
CN108278901A (zh) * 2018-01-30 2018-07-13 上海荣丰科技发展有限公司 一种用于冶炼熔渣成分调整的工业炉及生产工艺

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JARDY A: "Magnetohydrodynamic and thermal behavior of electroslag remelting slags", 《METALLURGICAL TRANSACTIONS B》 *
KHARICHA A: "Influence of an imposed vertical current on the droplet formation during a melting process", 《THE SEVENTH INTERNATIONAL CONFERENCE ON MULTIPHASE FLOW TAMPA》 *
ZHONGJUN YANG: "Prediction of melt rate of vibrating-electrode Electroslag Remelting process using artificial neural network", 《2015 5TH INTERNATIONAL CONFERENCE ON INFORMATION SCIENCE AND TECHNOLOGY (ICIST)》 *
张华: "电渣熔铸异形件温度场数值模拟技术", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 *
王晓花: "高中压转子钢电渣重熔过程宏微观多尺度数值模拟及应用研究", 《中国博士学位论文全文数据库工程科技I辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113934139A (zh) * 2020-06-29 2022-01-14 宝武特种冶金有限公司 一种基于在线仿真模型的真空电弧重熔过程熔化速度控制方法
CN113934139B (zh) * 2020-06-29 2023-06-16 宝武特种冶金有限公司 基于在线仿真模型的真空电弧重熔过程熔化速度控制方法
CN113987892A (zh) * 2021-12-28 2022-01-28 北京钢研高纳科技股份有限公司 一种控制高温合金偏析的真空电弧重熔3d模型及控制方法
WO2023142422A1 (zh) * 2022-01-28 2023-08-03 苏州大学 单进多出式水冷结晶器电渣重熔装置及电渣重熔的方法

Also Published As

Publication number Publication date
CN108984918B (zh) 2023-04-18

Similar Documents

Publication Publication Date Title
CN108984918B (zh) 一种电渣重熔自耗电极熔化速率的预测方法
CN108897972A (zh) 一种电渣重熔铸锭凝固微观组织的预测方法
Yang et al. Numerical analysis for electromagnetic field influence on heat transfer behaviors in cold crucible used for directional solidification
Du et al. Effect of operating conditions on inclusion of die steel during electroslag remelting
Yao et al. Effects of mold electromagnetic stirring on heat transfer, species transfer and solidification characteristics of continuous casting round billet
Cui et al. Comparative investigation on ingot evolution and product quality under different arc distributions during vacuum arc remelting process
Zuo et al. Study on the sump and temperature field during low frequency electromagnetic casting a superhigh strength Al–Zn–Mg–Cu alloy
LIU et al. Modeling progress of high-temperature melt multiphase flow in continuous casting mold
Qi et al. Effect of directional solidification in electroslag remelting on the microstructure and cleanliness of an austenitic hot-work die steel
Kelkar et al. Computational modeling of electroslag remelting processes
Tang et al. Effects of annulus gap on flow and temperature field in electromagnetic direct chill casting process
Wen et al. Numerical Simulation on the oxidation of lanthanum during the Electroslag remelting process
Silva et al. The columnar to equiaxed transition during the horizontal directional solidification of Sn–Pb alloys
Bojarevics et al. Travelling magnetic field mixing for particle dispersion in liquid metal
Li et al. Numerical simulation of a new electroslag remelting technology with current conductive stationary mold
CN115125396B (zh) 采用var工艺对特殊钢枝晶细化及夹杂物去除的控制方法
CHEN et al. Temperature field calculation on cold crucible continuous melting and directional solidifying Ti50Al alloys
Chen et al. Hot test and simulation of ESR hollow ingots formation in current supplying mould with electrodes change
Malinowski et al. Modelling of the strand and mold temperature in the continuous steel caster
Hua et al. Numerical simulation of the solidification process of Cu-0.45% Sn alloy in upward continuous casting
Kharicha et al. Selected numerical investigations on ESR process
Zhang et al. Effect of pull-down rate and power on electromagnetic separation of hypereutectic Ti–Si alloy under vacuum
Bojarevics et al. Cold crucible melting with bottom pouring nozzle
Ludwig et al. Important Key Process Simulations in the Field of Steel Metallurgy
Protopopov et al. Increasing the efficiency of refining rail steel in intermediate CCM ladles based on rational organization of hydrodynamic processes

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
GR01 Patent grant
GR01 Patent grant