第40卷 第3期 上 海 金 属 Vol.40,No.390 2018年5月 SHANGHAI METALS May ,2018基金项目:省部共建高品质特殊钢冶金与制备国家重点实验室开放课题(SKLASS 2015⁃04)作者简介:赵静,女,副教授,博士,研究方向为金属凝固理论及控制HPMO 作用下铸锭内部电磁场㊁流场和温度场的数值模拟赵 静1,2 徐智帅2(1.唐山学院机电工程系,河北唐山 063000;2.上海大学省部共建高品质特殊钢冶金与制备国家重点实验室,上海 200072) 【摘要】 建立了冒口脉冲磁致振荡(HPMO )作用下18t 铸锭凝固过程的数学模型,运用商业软件ANSYS 对铸锭凝固过程的电磁场㊁流场和温度场的分布进行了数值模拟㊂研究结果表明,HPMO 在铸锭内部产生的电磁效应主要集中在铸锭冒口区域,并随时间和空间不断变化,该结果为HPMO 在大铸锭工业生产中的应用提供了重要的参考依据㊂【关键词】 冒口脉冲磁致振荡技术 电磁场 流场 温度场 数值模拟Numerical Simulation of Electromagnetic ,Flow and Temperature Distribution during Ingot Solidification Process under HPMOZhao Jing 1,2 Xu Zhishuai 2(1.Mechanical and Electronic Engineering ,Tangshan College ,Tangshan Hebei 063000,China ;2.State Key Laboratory of Advanced Special Steel ,Shanghai University ,Shanghai 200072,China ) 【Abstract 】 The mathematical model was built to describe the solidification process of an 18tingot under HPMO.The distributions of electromagnetic field ,flow field and temperature field were simulated with the commercial ANSYS finite element software.The simulation results showed that the electromagnetic effects of HPMO concentrated mainly in the hot⁃top of ingot and changed with the time and space ,which would provide important references for the application of HPMO in ingot production.【Key Words 】 hot⁃top pulsed⁃magneto oscillation (HPMO ),electromagnetic filed ,flowfiled ,temperature ,numerical simulation 由于大铸锭散热缓慢㊁凝固时间长,凝固过程中容易形成严重的晶粒粗大㊁组织不均匀㊁宏观偏析等缺陷㊂这些缺陷不仅为后续加工造成很大困难,还会成为产品使用中的潜在隐患[1],但是人们至今没有找到改善大铸锭铸造缺陷的理想方法㊂上海大学翟启杰等开发的脉冲磁致振荡技术(pulsed⁃magneto oscillation,简称PMO)通过电磁效应促进固液界面和熔体表面附近形核并形成 结晶雨”,从而细化整个铸坯和铸锭的凝固组织[2⁃4]㊂目前,该技术已经应用于连铸生产,但在模铸中的应用仍处于探索阶段㊂在此基础上,他们又提出了适用于铸锭生产的PMO 技术 冒口脉冲磁致振荡技术(hot⁃top pulsed⁃magnetooscillation,简称HPMO)㊂大铸锭的凝固过程十分复杂,无法采用试验研究的方法探明铸锭内部电磁场㊁流场和温度场的分布情况㊂本文采用数值模拟的方法研究了HPMO 作用下18t 铸锭内部电磁场㊁流场和温度场的分布,其结果可为HPMO 技术在大铸锭生产中的应用提供一定的理论依据㊂1 数学模型的建立(1)几何模型的建立数学模拟采用的几何模型以某钢厂18t 铸锭尺寸为原型,在计算模拟过程中,共分为铸锭㊁铸型㊁线圈㊁硅钢片㊁空气5个区域㊂单元类型为ANSYS 中的PLANE53[7],如图1所示㊂第3期 赵 静等:HPMO 作用下铸锭内部电磁场㊁流场和温度场的数值模拟91 图1 几何模型Fig.1 Geometric model (2)计算过程首先,根据实际生产设备提取物理模型,在一定的假设条件下进行合理的简化,获得进行数值模拟的几何模型;其次,运用ANSYS 软件对HPMO 作用下熔体内部的电磁场进行数值模拟,获取熔体内部电磁力和焦耳热等电磁分量的分布;最后,将电磁力和焦耳热作为动量方程和能量方程的源项以此研究电磁场对流场和温度场的影响㊂(3)控制方程麦克斯韦方程组能够描述一切宏观电磁场学问题㊂本次模拟以麦克斯韦方程组(式(1))为出发点,脉冲电流是由RLC 振荡回路产生的,电流随时间的变化关系如式(2)所示[8]㊂▽×→E =-∂→B ∂t,▽×→H =∂→D ∂t +→J ,▽㊃→D =0,▽㊃→B =0(1)I (t )=I 0exp(-ξt )sin(ωt )(2)描述流体湍流流动的方程为连续性方程㊁动量方程和K ⁃ε方程,可以用式(3)~式(6)进行描述[9]㊂∂ρ∂t +ρ∂u i ∂x i=0(3)∂(ρu i )∂t +∂(ρu i u j )∂x j =-∂P ∂x j +∂∂x j μ∂u i ∂x æèçöø÷j +∂∂x j μ∂u j ∂x æèçöø÷j +S(4)∂∂t (ρK )+∂∂x j (ρu j K )=∂∂x j μ+μt σk ㊃∂K ∂x æèçöø÷j +G K -ρε(5)∂∂t (ρε)+∂∂x j (ρu j ε)=∂∂x j μ+μt σε㊃∂ε∂x æèçöø÷j +εK (c 1G k-c 2ρε)(6)式中:ρ为钢液的密度,u i 为沿坐标轴方向的速度矢量,X j 为坐标轴的方向,P 为压力,S 为动量源,K 为湍动能,ε为湍动能耗散率㊂将电磁力F 以动量源的形式加载到动量方程中进行计算,由式(7)计算得到㊂F =→J ×→B(7) 运用ANSYS 程序中的FLOTRAN CFD 模块进行流体计算,所用计算单元为FLUID141,其结构特征与PLANE13(电磁场计算单元)类似㊂将电磁场模型中的钢液部分赋予FLUID141单元属性和流体属性,其他部分赋予零单元属性即得到流场计算的有限元模型,从而实现电磁场与流场的耦合计算㊂HPMO 作用下钢液内部的温度分布可以看成是具有内热源的热传导过程,可以用式(8)进行描述㊂ρC p ∂T ∂t +v x ∂T ∂x +v y ∂T ∂y +v z ∂T ∂()z=∂∂x k x∂T ∂()x +∂∂y k y∂T ∂()y +∂∂z k z∂T ∂()z+̇q V(8) 将焦耳热导入温度场计算模块中,从而得到HPMO 对铸锭内部温度分布的影响㊂ q =∫J (t )2σd t(9)式中:J 为电流密度,σ为电阻率㊂(4)初始条件和边界条件1)初始条件电磁场:本模拟中电磁场的计算周期为一个电脉冲的放电周期㊂计算开始时,各电磁量的初始值均为0㊂流场和温度场:由于充型过程相对于整个凝固过程而言,时间很短,且HPMO 为充型结束后对熔体进行处理㊂模拟中忽略充型过程对凝固过程的影响,计算开始时,设定熔体内各点的速度为0,温度为1550℃㊂2)边界条件电磁场:空气区域的外边界条件设为无限远边界条件,即A =0,φ=0㊂中心对称轴为超导边界条92 上 海 金 属 第40卷件,▽×A =0,∂A n∂n=0㊂不同介质之间,A 1=A 2,φ1=φ2,(σ1-σ2)∂A n ∂t =σ2∂φ∂æèçöø÷n 2-σ1∂φ∂æèçöø÷n 1㊂流场:液面处,法向流速V y =0㊂型壁处,V x =0,V y =0㊂对称轴,V x =0㊂温度场:液面处㊁型壁处只考虑对流换热,均采用第3类边界条件,即采用q =h (T -T a )的形式给出,其中T a 为环境温度,取25℃㊂液面处对流换热系数h =10w㊃(m 2㊃K)-1,型壁处h =1000w㊃(m 2㊃K)-1㊂对称面为绝热边界条件,q =0㊂(5)物性参数模拟所用电磁参数如表1所示㊂ 表1 电磁参数Table 1 Electromagnetic parameters相对磁导率电阻率熔体11.4e⁃6铸型1005e⁃6线圈12e⁃7硅钢片20000空气1 模拟所用材料为P 91钢,其热物性参数由ProCAST 的Level 模型计算得到㊂所有参数均随温度变化,变化曲线如图2所示㊂图2 P 91钢的热物性参数Fig.2 Thermal parameters of P 91steel2 电磁场分布2.1 磁感应强度的分布图3为一个周期内磁感应强度的分布规律㊂从图中可以看出,HPMO 作用下熔体内部磁感应强度大小和方向随时间和空间不断发生变化㊂在一个放电周期内的T /4和3T /4时刻,磁感应强度的最大值基本相同,但方向相反㊂当t >T 1(T 1为放电脉宽)时,熔体内部的磁场强度逐渐减弱,其大小不及放电期间的十分之一㊂从空间分布来看,由于趋肤效应的影响,HPMO作用下产生的感图3 一个周期内磁感应强度分布Fig.3 Distributions of magnetic field in a period应磁场主要集中在铸锭的上部,靠近冒口型壁线圈位置处磁感应强度的数值最大㊂2.2 焦耳热的分布变化的磁场会产生感应电流,感应电流会在熔体内部产生一定的热量㊂图4为一个周期内熔体内部焦耳热的分布情况㊂HPMO 产生的焦耳热主要集中在靠近线圈位置的熔体内部,相对整个铸锭,其作用区域很小㊂但其数值很大,在当前模拟条件下,其最大值达到2.2×107W /m 3㊂图4 一个周期焦耳热分布Fig.4 Distributions of Joule heat in a period第3期 赵 静等:HPMO 作用下铸锭内部电磁场㊁流场和温度场的数值模拟93 2.3 电磁力的分布电磁力为磁感应强度B 与感应电流I 的乘积,而磁感应强度与感应电流之间存在一定的相位差,因此电磁力的变化频率与磁感应强度不同㊂下面将一个放电周期内T /8㊁T /4㊁3T /8㊁T /2㊁5T /8㊁3T /4㊁7T /8㊁T 等8个时刻的电磁力分布进行分析,如图5所示㊂从整个放电周期来看,电磁力的大小和方向不断发生变化,但大部分时间指向熔体内部,在t =3T /4时,电磁力的数值最大,在t =T 时,电磁力数值最小㊂图5 一个周期内电磁力分布Fig.5 Distributions of electromagnetic force in a period3 流场分布HPMO 作用产生的电磁力会对熔体的流动产生一定的影响,该流动为强制流动㊂由图5可见,电磁力的作用区域主要分布在铸锭的冒口区域,因此HPMO 产生的强制流动也主要集中在冒口区域,在铸锭的下半部分流动很弱㊂图6为自降温开始后100㊁400㊁800㊁1200s 内铸锭内部的钢液流动情况㊂可以看出,HPMO 作用下在冒口区域形成了很强的回流区,靠近型壁处流速最大㊂由于电磁力的周期性变化,又因在0~1200s 内冒口处熔体温度一直处于液相线温度以上,因此铸锭内上半部分钢液的流动状态在0~1200s 内几乎保持不变㊂由于铸锭尺寸较大,内外散热不均匀,在熔体内部存在一定的温度场,引起密度的不同,会导致熔体内部产生一定的自然对流㊂因此,降温一段时间后,在熔体的下半部分形成了另一个回流区㊂图6 铸锭内部钢液流动随时间的变化Fig.6 Flow change over time in ingot4 温度场分布图7为自降温后100㊁400㊁800㊁1200s 内熔体内部温度变化,与无HPMO 处理时钢锭内部温度分布相比发现,HPMO 未改变铸锭的凝固顺序,仍然为 U”形凝固[10]㊂但HPMO 产生的焦耳热对冒口处钢液的温度产生了一定的影响,延长了冒口处钢液的凝固时间,增强了冒口的保温作用,有助于减少铸锭内部的缩孔,并利于熔体内部夹杂物的上浮㊂图7 铸锭内部温度分布Fig.7 Distributions of temperature in ingot 5 结论(1)HPMO 作用下电磁场分布主要集中在铸锭冒口区域,磁感应强度㊁电磁力和焦耳热均随输入电流的周期性变化而变化㊂(2)HPMO 产生的焦耳热延长了冒口处的凝固时间,增强了冒口的保温效果㊂(3)HPMO 产生的强制对流在铸锭上半部分94 上 海 金 属 第40卷形成一个回流区㊂参考文献[1]黄天佑,刘小刚,康进武,等.我国大型铸钢件生产的现状与关键技术[J].铸造,2007,56(9):899⁃904.[2]LI B,YIN Z,GONG Y Y,et al.Effect of temperature field onsolidification structure of pure Al under pulse magneto⁃oscillation [J].China Foundry,2011,8(2):172⁃176.[3]LIANG D,LIANG Z,ZHAI Q J,et al.Nucleation and grainformation of pure Al under Pulsed Magneto⁃Oscillation treatment[J].Materials Letters,2014,130:48⁃50.[4]徐智帅,李祺欣,梁柱元,等.脉冲磁致振荡下Al⁃4.5wt%Cu合金微观组织形态[J].上海金属,2015,37(2):31⁃35.[5]李祺欣,俞基浩,梁冬,等.脉冲磁致振荡细化工业纯铝机制研究[J].上海金属,2015,37(4):48⁃51.[6]EDRY I,MORDECHAI T,FRAGE N.Effects of treatmentduration and cooling rate on pure aluminum solidification upon pulse magneto⁃oscillation treatment [J ].Metallurgical and Materials Transactions A,2016,47(3):1261⁃1267.[7]Inc ANSYS.ANSYS 耦合场分析指南[M].SASI Inc,1998.[8]简柏敦.导电与导磁物质中的电磁场[M].北京:人民教育出版社,1981.[9]郭鸿志,张欣欣,刘向军,等.传输过程数值模拟[M].北京:冶金工业出版社,1998.[10]赵静,周湛,张捷宇,等.钢锭凝固过程中温度场和流场的数值模拟研究[J].上海金属,2014,36(1):55⁃58.收修改稿日期:2017⁃08⁃14(上接第89页)[4]杨天钧,段国锦,周渝生,等.高炉无料钟布料炉料分布预测模型的开发研究[J].钢铁,1991,26(11):10⁃14,28.[5]陈令坤.武钢新3号高炉无钟布料料面预测模型的开发[D].武汉:武汉科技大学,1993.[6]高征凯,杨天钧,赵永福,等.无钟布料仿真模型在美钢联Fairfield 厂8号高炉的应用[C]//全国炼铁生产技术暨炼铁年会文集.北京:中国金属学会,2002:646⁃653.[7]经文波,陈小雷.高炉应用数学模型研究[J].江西冶金,2003,23(6):72⁃74,83.[8]ZHAO H T ,ZHU M H,DU P,et al.Uneven distribution ofburden materials at blast burnace top in bell⁃less top with parallel bunkers[J].ISIJ International,2012,52(12):2177⁃2185.[9]ZHAO H T,DU P,REN L Q,et al.Simultaneouslyoverdeveloped central and peripheral gas flow of a blast furnace [J].ISIJ International,2015,55(10):2064⁃2073.[10]杜鹤桂,余艾冰.高炉无钟炉顶布料的模拟试验研究[J].钢铁,1986,21(11):1⁃9.[11]张建良,张雪松,张汝旺,等.无钟炉顶多环布料料面形状的模拟[J].过程工程学报,2008,8(S 1):135⁃139.[12]李志全,贾文君,毕耜友,等.邯钢高炉无钟布料实物模型研制及应用[C]//全国炼铁生产技术会议暨炼铁学术年会文集(下).北京:中国金属学会,2012:563⁃569.[13]JIMENEZ J,MOCHON J,FORMOSO A,et al.Burdendistribution analysis by digital image processing in a scale modelof a blast furnace shaft[J].ISIJ International,2000,40(2):114⁃120.[14]郭可信.宝钢2号高炉无料钟布料实践[J].宝钢技术,1993(1):5⁃10,34.[15]杜鹤桂.高炉无钟布料的重要环节 平台的形成[J].炼铁,1995,14(3):33⁃36.收修改稿日期:2017⁃09⁃22。
