基于实验与3D-CAFE法的高硅钢铸锭凝固行为

工程科学学报，第39卷，第3期：360-368,2017年3月 Chinese Journal of Engineering,Vol.39,No.3:360-368,March 2017 DOI:10.13374/j.issn2095-9389.2017.03.007:http://journals.ustb.edu.cn 基于实验与3D一CAE法的高硅钢铸锭凝固行为 宋炜，张炯明区，王顺玺 北京科技大学钢铁治金新技术国家重点实验室，北京100083 ☒通信作者，E-mail:jmz2233@sina.com 摘要通过空冷和水冷实验研究了高硅钢的铸态组织，发现高硅钢铸态组织主要由粗大的柱状晶构成，水冷铸锭中柱状晶 比例高达90%以上.依据铸锭的化学成分和晶粒统计结果，确定了3 D-CAFE法模拟所需的枝晶生长动力学系数及高斯分布 等参数.采用CAFE法对不同冷却条件下高硅钢的凝固过程进行模拟研究，发现空冷铸锭较水冷铸锭的温度场更均匀，糊状 区更宽阔：空冷铸锭呈“过渡式”凝固，水冷铸锭呈“分层式”凝固：空冷流场较水冷流场更稳定，凝固末期冒口处出现明显的 抽吸现象，而水冷模拟结果中未观察到该现象.组织模拟结果发现，模拟得到的高硅钢凝固组织无论是形貌还是晶粒尺寸都 与实验结果相一致：最后通过改变浇注温度模拟研究了过热度对高硅钢凝固组织的影响，结果表明，随着过热度的降低，铸锭 中心等轴晶率提高，晶粒数量增加，晶粒尺寸变得细小 关键词硅钢：凝固：温度场：流动速率 分类号TG142.71 Solidification behavior of high-silicon steel based on experimental and 3D-CAFE method SONG Wei,ZHANG Jiong-ming,WANG Shun-xi State Key Laboratory of Advanced Metallurgy,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:jmz2233@sina.com ABSTRACT The as-cast structure of high-silicon steel ingots under different cooling conditions was studied in this paper.It is found that the as-cast structure of the ingot is formed mainly by coarse columnar crystals,especially in the water cooling ingot,and the ratio is reaches as 90%.The dendrite tip growth kinetic coefficients and Gauss distribution parameters for 3D-CAFE simulation were deter- mined according to the compositions of high-silicon steel and the results of as-east structure.Then the solidification process of high-sil- icon steel under different cooling conditions was simulated by 3D-CAFE method.The results show that the temperature field under air cooling is more uniform,the mushy zone is broader,and it exhibits a transitional solidification pattern,however,which shows a lay- ered solidification pattern under water cooling.The flow field under air cooling is more stable than that under water cooling,there is a remarkable suction region within the feeder head of the air cooling ingot,and this phenomenon is not observed in the water cooling one.The CAFE results including both morphology and grain size show a good agreement with the results from experiments.Moreover, the influence of superheat on the solidification structures was researched to find that the ratio and quantity of equiaxed structures in- crease with the decrease of the superheat,and the grain size becomes finer KEY WORDS silicon steel:solidification:temperature field:flow rate 硅钢片作为一种最重要的软磁材料被广泛应用于 右田.磁电转换过程中采取高频化可以实现电气设备 电力行业，用量占磁性材料总量的90%~95%左 的高效化和小型轻便化，然而高频化会带来铁损及磁 收稿日期：201605一12

工程科学学报，第 39 卷，第 3 期: 360--368，2017 年 3 月 Chinese Journal of Engineering，Vol． 39，No． 3: 360--368，March 2017 DOI: 10． 13374 /j． issn2095--9389． 2017． 03． 007; http: / /journals． ustb． edu． cn 基于实验与 3D--CAFE 法的高硅钢铸锭凝固行为 宋 炜，张炯明，王顺玺 北京科技大学钢铁冶金新技术国家重点实验室，北京 100083 通信作者，E-mail: jmz2233@ sina． com 摘 要 通过空冷和水冷实验研究了高硅钢的铸态组织，发现高硅钢铸态组织主要由粗大的柱状晶构成，水冷铸锭中柱状晶 比例高达 90% 以上． 依据铸锭的化学成分和晶粒统计结果，确定了 3D--CAFE 法模拟所需的枝晶生长动力学系数及高斯分布 等参数． 采用 CAFE 法对不同冷却条件下高硅钢的凝固过程进行模拟研究，发现空冷铸锭较水冷铸锭的温度场更均匀，糊状 区更宽阔; 空冷铸锭呈“过渡式”凝固，水冷铸锭呈“分层式”凝固; 空冷流场较水冷流场更稳定，凝固末期冒口处出现明显的 抽吸现象，而水冷模拟结果中未观察到该现象． 组织模拟结果发现，模拟得到的高硅钢凝固组织无论是形貌还是晶粒尺寸都 与实验结果相一致; 最后通过改变浇注温度模拟研究了过热度对高硅钢凝固组织的影响，结果表明，随着过热度的降低，铸锭 中心等轴晶率提高，晶粒数量增加，晶粒尺寸变得细小． 关键词 硅钢; 凝固; 温度场; 流动速率 分类号 TG142. 71 Solidification behavior of high-silicon steel based on experimental and 3D--CAFE method SONG Wei，ZHANG Jiong-ming  ，WANG Shun-xi State Key Laboratory of Advanced Metallurgy，University of Science and Technology Beijing，Beijing 100083，China Corresponding author，E-mail: jmz2233@ sina． com ABSTＲACT The as-cast structure of high-silicon steel ingots under different cooling conditions was studied in this paper． It is found that the as-cast structure of the ingot is formed mainly by coarse columnar crystals，especially in the water cooling ingot，and the ratio is reaches as 90% ． The dendrite tip growth kinetic coefficients and Gauss distribution parameters for 3D--CAFE simulation were deter￾mined according to the compositions of high-silicon steel and the results of as-cast structure． Then the solidification process of high-sil￾icon steel under different cooling conditions was simulated by 3D--CAFE method． The results show that the temperature field under air cooling is more uniform，the mushy zone is broader，and it exhibits a transitional solidification pattern，however，which shows a lay￾ered solidification pattern under water cooling． The flow field under air cooling is more stable than that under water cooling，there is a remarkable suction region within the feeder head of the air cooling ingot，and this phenomenon is not observed in the water cooling one． The CAFE results including both morphology and grain size show a good agreement with the results from experiments． Moreover， the influence of superheat on the solidification structures was researched to find that the ratio and quantity of equiaxed structures in￾crease with the decrease of the superheat，and the grain size becomes finer． KEY WOＲDS silicon steel; solidification; temperature field; flow rate 收稿日期: 2016--05--12 硅钢片作为一种最重要的软磁材料被广泛应用于 电力 行 业，用量占磁性材料总量的 90% ～ 95% 左 右［1］． 磁电转换过程中采取高频化可以实现电气设备 的高效化和小型轻便化，然而高频化会带来铁损及磁

宋炜等：基于实验与3D-CAFE法的高硅钢铸锭凝固行为 361 致伸缩大幅增长，从而导致能量的巨大损失和噪声的 (3)能量守恒方程 急剧增大回.实践证明6.5%Si高硅钢在1000~ 10000Hz频率范围内磁化时与传统硅钢片相比0：铁 p+p7+"+ 损P低40%~70%，导磁率μ高几倍至几十倍，磁致 (4) 伸缩入，几乎为零，使用频率越高，其磁学性能优势越 (a)+()+是() 为突出.因此，高硅钢被认为是下一代最具有应用前 其中m: H=∫e,dr+1-f) (5) 景的软磁材料因 式中：u、v和0分别为x、y和z方向的速度分量，m· 然而随着$í含量的提高，高硅钢变得又脆又硬， st为时间，s为液相率：为固相率p为压力， 很难采用传统工艺进行制造.近年来，国内外学者对 Pa;g,为x方向重力分量，ms2:p为密度，kg"m3:μ 高硅钢进行改性处理或改进轧制技术来提高其可轧 性，但均因成本高、产量低、污染环境等缺点没能 为绝对黏度，Pas;K为渗透率，m2;入为热传导率，W· 规模化应用于生产.目前限制高硅钢采用传统工艺生 (mK);c为定压比热容，J·(kgK)1;L为凝固潜 产的一个重要因素是难以获得质量优良的铸锭，铸锭 热，Jkg;T为节点温度，K;H为热焓，Jmol 凝固组织中晶粒粗大、柱状晶发达，在浇注和后续加工 1.2微观模型 过程中容易产生裂纹及瓦垅状缺陷四.因此，研究高 1.2.1异质形核模型 硅钢的凝固特性及其影响因素具有理论和实际意义. 微观组织模拟过程中异质形核包括两种形式▣： 凝固过程模拟采用的微观组织模型一般包括确定 瞬时形核和连续形核.CAFE法采用连续形核模型来 性模型、随机模型和相场模型-，其中，随机模型采 处理形核，假设形核现象发生在一系列不同的形核位 用随机方法处理晶粒形核和长大，能够模拟晶粒竞争 置上，核密度的变化采用连续分布函数dn/d(△T)来 生长机制n.Gandin等2国将有限元FE和元胞自动 描述，即Gauss分布函数： 机CA结合起来建立了宏观-微观相耦合的CAFE法 dn a7-2mA (△T-AT)21 -exp (6) 随机模型，物理机制明确、计算速度快、模拟精度高，被 2A7 国内外学者-阁广泛接受并采用.本文采用实验结合 式中，n为正态分布从0到∞得到的最大形核密度， 模拟的办法对不同冷却条件下高硅钢凝固过程中温度 面形核密度n,m2,体形核密度n,m3;△T为平均形 场、凝固过程、流场凝固组织等凝固特征进行研究，对 核过冷度，K;△T.为形核过冷度标准偏差，K;n为晶粒 制造高硅钢铸锭具有一定的指导意义 密度.过冷度△T的增加时，晶粒密度dn也随之增加， 当给定过冷度△T时，形成的晶核密度n(△T)可由下 1数学物理模型 式求得 1.1宏观模型 n(△)=，dA刀d(A). dn (7) 1.1.1热物性计算 本文采用双混合模型网计算热物性参数，如下式 1.2.2生长动力学模型 描述： 合金凝固过程中，枝晶尖端生长主要受过冷度的 P=∑B+Σ三璃∑g-(① 影响，总过冷度△T可由下式表示： △T=△T.+△T.+△T,+△T (8) 式中：P为相的特性：P为该相中元素i特性：为二 式中，△T。、AT、△T△T,分别为成分过冷度、热力学过 元相互作用参数：xx,分别为元素i寸在该相中的摩尔 冷度、固液界面曲率过冷度、生长动力学过冷度，K.通 分数；v为决定二元相互作用参数的变量. 常情况下△T,、△T,、△T相比△T很小，在计算中可以忽 1.1.2控制方程 略。柱状晶和等轴晶的生长速度采用KGT模型描述四： (1)质量守恒方程. v=a2△T+a3△T. (9) 坐+a(pu+apml+apw=0. (2) ay 式中：a2a是生长动力学系数，m·(sK)l (2)动量守恒方程（只列出x方向，y,:方向与x 1.3FE与CA模型耦合 方向类似) CAFE模型中定义了FE节点和CA元胞之间的插 值因子，并引入凝固潜热的影响，确保了微观组织是温 度场的函数（图1所示）.CA的元胞v与有限元节点 品（片）+品（片）+是（片）-（发） ij小k之间分别具有非零的插值因子中中，中，这些 因子结合E节点的温度就可以确定网格中元胞的温 (3) 度，并对枝晶形核、生长过程释放的潜热求和，更新节

宋 炜等: 基于实验与 3D--CAFE 法的高硅钢铸锭凝固行为 致伸缩大幅增长，从而导致能量的巨大损失和噪声的 急剧 增 大［2］． 实践 证 明 6. 5% Si 高 硅 钢 在 1000 ～ 10000 Hz 频率范围内磁化时与传统硅钢片相比［1］: 铁 损 Pt低 40% ～ 70% ，导磁率 μ 高几倍至几十倍，磁致 伸缩 λs几乎为零，使用频率越高，其磁学性能优势越 为突出． 因此，高硅钢被认为是下一代最具有应用前 景的软磁材料［3］． 然而随着 Si 含量的提高，高硅钢变得又脆又硬， 很难采用传统工艺进行制造． 近年来，国内外学者对 高硅钢进行改性处理或改进轧制技术来提高其可轧 性［4--6］，但均因成本高、产量低、污染环境等缺点没能 规模化应用于生产． 目前限制高硅钢采用传统工艺生 产的一个重要因素是难以获得质量优良的铸锭，铸锭 凝固组织中晶粒粗大、柱状晶发达，在浇注和后续加工 过程中容易产生裂纹及瓦垅状缺陷［1］． 因此，研究高 硅钢的凝固特性及其影响因素具有理论和实际意义． 凝固过程模拟采用的微观组织模型一般包括确定 性模型、随机模型和相场模型［7--10］，其中，随机模型采 用随机方法处理晶粒形核和长大，能够模拟晶粒竞争 生长机制［11］． Gandin 等［12--13］将有限元 FE 和元胞自动 机 CA 结合起来建立了宏观--微观相耦合的 CAFE 法 随机模型，物理机制明确、计算速度快、模拟精度高，被 国内外学者［14--18］广泛接受并采用． 本文采用实验结合 模拟的办法对不同冷却条件下高硅钢凝固过程中温度 场、凝固过程、流场凝固组织等凝固特征进行研究，对 制造高硅钢铸锭具有一定的指导意义． 1 数学物理模型 1. 1 宏观模型 1. 1. 1 热物性计算 本文采用双混合模型［19］计算热物性参数，如下式 描述: P = ∑i xiPi + ∑i ∑ ji xixj ∑υ Ωυ ij( xi － xj ) υ ． ( 1) 式中: P 为相的特性; Pi 为该相中元素 i 特性; Ωυ ij 为二 元相互作用参数; xi、xj分别为元素 i、j 在该相中的摩尔 分数; υ 为决定二元相互作用参数的变量． 1. 1. 2 控制方程 ( 1) 质量守恒方程． ρ t + ( ρu) x + ( ρv) y + ( ρw) z = 0． ( 2) ( 2) 动量守恒方程( 只列出 x 方向，y，z 方向与 x 方向类似) ． ρ fl u t + ρ f 2 ( l u u x + v u y + w u  ) z = － p x + ρgx +   ( x u fl u  ) x +   ( y u fl u  ) y +   ( z u fl u  ) z ( － μ ) K u ． ( 3) ( 3) 能量守恒方程． ρ H t + ρ H  ( T u T x + v T y + w T  ) z =   ( x λ T  ) x +   ( y λ T  ) y +   ( z λ T  ) z ． ( 4) 其中［20］: H = ∫ T 0 cpdT + L( 1 － fs) ． ( 5) 式中: u、v 和 w 分别为 x、y 和 z 方向的速度分量，m· s － 1 ; t 为时间，s; fl 为液相率; fs 为固相率; p 为压力， Pa; gx 为 x 方向重力分量，m·s － 2 ; ρ 为密度，kg·m － 3 ; μ 为绝对黏度，Pa·s; K 为渗透率，m2 ; λ 为热传导率，W· ( m·K) － 1 ; cp为定压比热容，J·( kg·K) － 1 ; L 为凝固潜 热，J·kg － 1 ; T 为节点温度，K; H 为热焓，J·mol． 1. 2 微观模型 1. 2. 1 异质形核模型 微观组织模拟过程中异质形核包括两种形式［11］: 瞬时形核和连续形核． CAFE 法采用连续形核模型来 处理形核，假设形核现象发生在一系列不同的形核位 置上，核密度的变化采用连续分布函数 dn /d( ΔT) 来 描述，即 Gauss 分布函数: dn d( ΔT) = nmax 槡2πΔTσ [ exp － ( ΔT － ΔTn ) 2 2ΔT2 ] σ ． ( 6) 式中，nmax为正态分布从 0 到∞ 得到的最大形核密度， 面形核密度 ns，m － 2，体形核密度 nv，m － 3 ; ΔTn为平均形 核过冷度，K; ΔTσ为形核过冷度标准偏差，K; n 为晶粒 密度． 过冷度 ΔT 的增加时，晶粒密度 dn 也随之增加， 当给定过冷度 ΔT 时，形成的晶核密度 n( ΔT) 可由下 式求得 n( ΔT) = ∫ ΔT 0 dn d( ΔT) d( ΔT) ． ( 7) 1. 2. 2 生长动力学模型 合金凝固过程中，枝晶尖端生长主要受过冷度的 影响，总过冷度 ΔT 可由下式表示: ΔT = ΔTc + ΔTt + ΔTr + ΔTk ( 8) 式中，ΔTc、ΔTt、ΔTr、ΔTk分别为成分过冷度、热力学过 冷度、固液界面曲率过冷度、生长动力学过冷度，K． 通 常情况下 ΔTt、ΔTr、ΔTk相比 ΔTc很小，在计算中可以忽 略． 柱状晶和等轴晶的生长速度采用 KGT 模型描述［21］: v = a2ΔT2 + a3ΔT3 ． ( 9) 式中: a2、a3是生长动力学系数，m·( s·K) － 1 ． 1. 3 FE 与 CA 模型耦合 CAFE 模型中定义了 FE 节点和 CA 元胞之间的插 值因子，并引入凝固潜热的影响，确保了微观组织是温 度场的函数( 图 1 所示) ． CA 的元胞 ν 与有限元节点 i、j、k 之间分别具有非零的插值因子 Φνi、Φνj 、Φνk，这些 因子结合 FE 节点的温度就可以确定网格中元胞的温 度，并对枝晶形核、生长过程释放的潜热求和，更新节 · 163 ·

·362· 工程科学学报，第39卷，第3期 点温度四 2实验过程及结果 实验采用7kg真空感应炉，浇注温度为1773K 单元/ 共开展2炉实验，一炉铸锭在空气中自然冷却，另一炉 铸锭在循环水箱中快速水冷.铸件实物如图2所示， 图2(a)中1为空冷铸件，2为水冷铸件：图2(b)和 (c)分别为两个铸锭冒口处的一次缩孔形貌，图2(d) 和(e)分别为两种冷却条件下铸锭一次缩孔的纵剖面 结果.可以看出，空冷铸锭一次缩孔最深处约为18 mm,水冷铸锭一次缩孔最深达33mm,两个铸锭的一 图1FE节点与CA节点关系示意图四 次缩孔均未延伸到本体.下文中凝固组织的模拟研究 Fig.1 Schematic diagram of the relationship between FE nodes and 不考虑冒口的影响，于是将其切除.铸锭中主要化学 CA nodes 元素含量检测结果见表1. a (b) d 20mm (c) 2 (e) 30mm 20mm 图2铸锭实物图.(a)空冷铸锭和水冷铸锭：(b)和()分别为空冷和水冷铸锭冒口处的一次缩孔形貌：(d)和()分别为空冷和水冷 铸锭一次缩孔的解剖结果 Fig.2 Actual ingot:(a)ingots under air cooling and water cooling:(b,c)morphologies of the cavity in feeder head of the ingots under air cooling and water cooling respectively:(d,e)longitudinal profiles of the cavity under air cooling and water cooling,respectively 表1母料及高硅钢铸锭中的化学成分（质量分数） Table 1 Results of the element composition of raw material and high sil- icon steel 呢 试样 C Mn P S Cr Ni 母料3.140.00250.220.0090.00070.0200.01 铸锭6.570.00250.210.0090.00070.0190.01 将高硅钢铸锭纵向剖开、打磨抛光，用5%硝酸酒 精侵蚀，低倍结果如图3所示.由图3可以看出，高硅 (a (b) 钢铸态组织主要由高比例且尺寸粗大的柱状晶构成， 空冷铸锭中等轴晶比例约为40%，而水冷铸锭中等轴 图3铸锭低倍组织.(a)空冷：(b)水冷 Fig.3 As-cast structure morphology of the ingots:(a)air cooling: 晶比例仅为10%左右.发达的柱状晶有利与最终产品 (b)water cooling 的磁性能3，但也会对后续加工造成不良影响. 从空冷铸锭边部到中心取4块10mm×10mm× 序分布、晶粒随机生长，晶粒内部主要由二次、三次枝 10m微观组织试样，打磨、抛光后采用10%的硝酸酒 晶构成 精侵蚀.微观组织结果如图4所示，图4(a)和(b)分 3数值模拟 别为高硅钢的柱状晶和等轴晶组织，可以看出高硅钢 柱状晶粗大，晶粒间晶界平行，具有明显的方向性，晶 3.1过程及参数 粒内部主要由发达的一次枝晶构成：而等轴晶晶界无 实体建模与网格划分：采用Po/E建立3D几何模

工程科学学报，第 39 卷，第 3 期 点温度［22］． 图 1 FE 节点与 CA 节点关系示意图［22］ Fig． 1 Schematic diagram of the relationship between FE nodes and CA nodes 2 实验过程及结果 实验采用 7 kg 真空感应炉，浇注温度为 1773 K． 共开展 2 炉实验，一炉铸锭在空气中自然冷却，另一炉 铸锭在循环水箱中快速水冷． 铸件实物如图 2 所示， 图 2( a) 中 1# 为空冷铸件，2# 为水冷铸件; 图 2 ( b) 和 ( c) 分别为两个铸锭冒口处的一次缩孔形貌，图 2( d) 和( e) 分别为两种冷却条件下铸锭一次缩孔的纵剖面 结果． 可以看出，空冷铸锭一次缩孔最深处约为 18 mm，水冷铸锭一次缩孔最深达 33 mm，两个铸锭的一 次缩孔均未延伸到本体． 下文中凝固组织的模拟研究 不考虑冒口的影响，于是将其切除． 铸锭中主要化学 元素含量检测结果见表 1． 图 2 铸锭实物图． ( a) 空冷铸锭和水冷铸锭; ( b) 和 ( c) 分别为空冷和水冷铸锭冒口处的一次缩孔形貌; ( d) 和 ( e) 分别为空冷和水冷 铸锭一次缩孔的解剖结果 Fig． 2 Actual ingot: ( a) ingots under air cooling and water cooling; ( b，c) morphologies of the cavity in feeder head of the ingots under air cooling and water cooling respectively; ( d，e) longitudinal profiles of the cavity under air cooling and water cooling，respectively 表 1 母料及高硅钢铸锭中的化学成分 ( 质量分数) Table 1 Ｒesults of the element composition of raw material and high sil￾icon steel % 试样 Si C Mn P S Cr Ni 母料 3. 14 0. 0025 0. 22 0. 009 0. 0007 0. 020 0. 01 铸锭 6. 57 0. 0025 0. 21 0. 009 0. 0007 0. 019 0. 01 将高硅钢铸锭纵向剖开、打磨抛光，用 5% 硝酸酒 精侵蚀，低倍结果如图 3 所示． 由图 3 可以看出，高硅 钢铸态组织主要由高比例且尺寸粗大的柱状晶构成， 空冷铸锭中等轴晶比例约为 40% ，而水冷铸锭中等轴 晶比例仅为 10% 左右． 发达的柱状晶有利与最终产品 的磁性能［23--24］，但也会对后续加工造成不良影响． 从空冷铸锭边部到中心取 4 块 10 mm × 10 mm × 10 mm 微观组织试样，打磨、抛光后采用 10% 的硝酸酒 精侵蚀． 微观组织结果如图 4 所示，图 4( a) 和( b) 分 别为高硅钢的柱状晶和等轴晶组织，可以看出高硅钢 柱状晶粗大，晶粒间晶界平行，具有明显的方向性，晶 粒内部主要由发达的一次枝晶构成; 而等轴晶晶界无 图 3 铸锭低倍组织 . ( a) 空冷; ( b) 水冷 Fig． 3 As-cast structure morphology of the ingots: ( a) air cooling; ( b) water cooling 序分布、晶粒随机生长，晶粒内部主要由二次、三次枝 晶构成． 3 数值模拟 3. 1 过程及参数 实体建模与网格划分: 采用 Pro /E 建立 3D 几何模 · 263 ·

宋炜等：基于实验与3D-CAFE法的高硅钢铸锭凝固行为 ·363· (a) (b) 树枝品 3mm 2mm 图4高硅钢铸锭凝固组织.()柱状品组织：(b)等轴品组织 Fig.4 Macrostructure of the high-silicon steel ingots:(a)columnar crystals:(b)equiaxed crystals 型，模型尺寸与铸锭实物尺寸一致.利用Procast中的 动力学参数计算模块计算出高硅钢生长动力学参数 Meshcast模块对模型进行面网格和体网格划分，结果 a2=2.183×10-7,a3=1.440×10-7 如图5所示. 表2二元铁基合金元素计算参数 Table 2 Calculated parameters of binary Fe-based alloy 质量分数/ 液相线 平衡分配自扩散系数/ 元素 % 斜率 系数 (109m2s1) C 0.0025 -78.00 0.340 2.0 Si 6.6100 -26.50 0.650 2.0 Mn 0.2100 -3.32 0.750 2.0 P 0.0090 -27.10 0.090 2.0 0.0007 -30.40 0.024 2.0 Al 0.9600 -5.00 0.420 2.8 Cr 0.0190 -2.61 0.760 2.0 Ni 0.0100 -1.60 0.940 2.0 图5模拟模型及尺寸参数 Cu 0.0100 -1.70 0.960 2.0 Fig.5 Simulation model and size parameters 0.0100 -6.00 0.230 31 材料设定：将双混合模型计算得到的热物性参数 3.2结果与讨论 代入到凝固计算模型中参与计算，高硅钢的固相线和 3.2.1温度场 液相线温度分别为1681K和1721K. 考虑到模型的对称性，将高硅钢凝固过程中相同 传热参数设定：不考虑充型过程，初始温度设为 时刻的温度场和凝固过程模拟结果耦合在一起，如 1773K与实验浇注温度一致.空冷和水冷的传热系数 图6和图7所示.而图8为两种冷却条件下铸锭不同 设定参考文献14]，分别为100W·(m2·K)-和1500 位置（图6(a))的冷却曲线.从图6和图7可以看 W·(m2.K),环境温度为303K 出，浇注结束后，由于高温金属液与低温的模壁接触产 高斯形核分布参数设定：依据ASTM标准推荐 生很大的温度梯度，温度迅速下降，当温度降至液相线 的计算方法n,=0.8n,而面形核最大密度(n.m)由 温度(T)时结晶开始，此时金属液会释放大量结晶潜 实验统计的晶粒密度确定为1.76×10m2,因此，计 热，当结晶释放的潜热与金属液散失的热量达到平衡 算体形核密度(n.)为5.91×10？m3:面形核平均 时，凝固过程会出现“回温”现象，即图8中冷却曲线 过冷度AT。=0.5K,面形核过冷度标准方差△T。= 中的温度“平台”，而靠近铸锭边部位置(P3)由于冷速 0.1K,体形核平均过冷度△T。=10K,体形核过冷度 大，凝固速度快，冷却过程来不及出现这种“平台”随 标准方差△T。=0.1K. 着温度不断下降，凝固首先从铸锭的角部开始，并沿铸 生长动力学系数设定：采用文献6-17]中数据 壁向中心逐渐推进，铸锭冒口逐渐向下塌陷，一次缩孔 (表2)，其中Si的平衡分配系数k以及液相线斜率m 随着凝固进程逐渐向着中心蔓延，空冷最终的一次缩 根据二元Fe-Si相图测定，应用商业软件CalcoSoft中 孔约20mm,水冷约为35mm,这与实验结果（图2(d)

宋 炜等: 基于实验与 3D--CAFE 法的高硅钢铸锭凝固行为 图 4 高硅钢铸锭凝固组织 . ( a) 柱状晶组织; ( b) 等轴晶组织 Fig． 4 Macrostructure of the high-silicon steel ingots: ( a) columnar crystals; ( b) equiaxed crystals 型，模型尺寸与铸锭实物尺寸一致． 利用 Procast 中的 Meshcast 模块对模型进行面网格和体网格划分，结果 如图 5 所示． 图 5 模拟模型及尺寸参数 Fig． 5 Simulation model and size parameters 材料设定: 将双混合模型计算得到的热物性参数 代入到凝固计算模型中参与计算，高硅钢的固相线和 液相线温度分别为 1681 K 和 1721 K． 传热参数设定: 不考虑充型过程，初始温度设为 1773 K 与实验浇注温度一致． 空冷和水冷的传热系数 设定参考文献［14］，分别为 100 W·( m2 ·K) － 1和 1500 W·( m2 ·K) － 1，环境温度为 303 K． 高斯形核分布参数设定［25］: 依据 ASTM 标准推荐 的计算方法 nv = 0. 8n3 /2 s ，而面形核最大密度 ( ns，max ) 由 实验统计的晶粒密度确定为 1. 76 × 105 m － 2，因此，计 算体形核密度( nv，max ) 为 5. 91 × 107 m － 3 ; 面形核平均 过冷度 ΔTs，n = 0. 5 K，面形核过冷度标准方差 ΔTs，σ = 0. 1 K，体形核平均过冷度 ΔTv，n = 10 K，体形核过冷度 标准方差 ΔTv，σ = 0. 1 K． 生长动力学系数设定: 采用文献［16--17］中数据 ( 表 2) ，其中 Si 的平衡分配系数 k 以及液相线斜率 m 根据二元 Fe--Si 相图测定，应用商业软件 CalcoSoft 中 动力学参数计算模块计算出高硅钢生长动力学参数 a2 = 2. 183 × 10 － 7，a3 = 1. 440 × 10 － 7 ． 表 2 二元铁基合金元素计算参数 Table 2 Calculated parameters of binary Fe-based alloy 元素 质量分数/ % 液相线 斜率 平衡分配 系数 自扩散系数/ ( 10 － 9 m2 ·s － 1 ) C 0. 0025 － 78. 00 0. 340 2. 0 Si 6. 6100 － 26. 50 0. 650 2. 0 Mn 0. 2100 － 3. 32 0. 750 2. 0 P 0. 0090 － 27. 10 0. 090 2. 0 S 0. 0007 － 30. 40 0. 024 2. 0 Al 0. 9600 － 5. 00 0. 420 2. 8 Cr 0. 0190 － 2. 61 0. 760 2. 0 Ni 0. 0100 － 1. 60 0. 940 2. 0 Cu 0. 0100 － 1. 70 0. 960 2. 0 V 0. 0100 － 6. 00 0. 230 3. 1 3. 2 结果与讨论 3. 2. 1 温度场 考虑到模型的对称性，将高硅钢凝固过程中相同 时刻的温度场和凝固过程模拟结果耦合在一起，如 图 6 和图 7 所示． 而图 8 为两种冷却条件下铸锭不同 位置 ( 图 6( a) ) 的冷却曲线． 从图 6 和图 7 可以看 出，浇注结束后，由于高温金属液与低温的模壁接触产 生很大的温度梯度，温度迅速下降，当温度降至液相线 温度( TL ) 时结晶开始，此时金属液会释放大量结晶潜 热，当结晶释放的潜热与金属液散失的热量达到平衡 时，凝固过程会出现“回温”现象，即图 8 中冷却曲线 中的温度“平台”，而靠近铸锭边部位置( P3) 由于冷速 大，凝固速度快，冷却过程来不及出现这种“平台”． 随 着温度不断下降，凝固首先从铸锭的角部开始，并沿铸 壁向中心逐渐推进，铸锭冒口逐渐向下塌陷，一次缩孔 随着凝固进程逐渐向着中心蔓延，空冷最终的一次缩 孔约 20 mm，水冷约为 35 mm，这与实验结果( 图 2( d) · 363 ·

·364 工程科学学报，第39卷，第3期 固相率 温度℃ Tliq ■1.000 1448 1480 0.933 14 0.867 0.800 402 150 0.733 1320 0.667 1280 0.600 1240 0.533 四☒ 1200 0.467 0.400 H20 0333 1080 0267 1040 0.200 1000 0.133 80 0.067 880 0 (a b (e) (d) 图6空冷条件下高硅钢铸锭不同时刻的凝固过程结果(a)和温度场(b~d).(a)99.22s:(b)204.23s:(c)274.23s:(d)339.23s Fig.6 Simulated solidification process (a)and temperature field(b-d)results of the high-silicon steel casting under air cooling condition at differ- ent time:(a）99.22s:(b)204.23s:(c)274.23s:(d)339.23s 固相率 温度℃ 1.000 1520 0.933 1426 0.867 1333 0.800 1240 0.733 7146 0.667 1053 0.600 0.533 0.467 73 0,400 0.333 0.267 493 0.200 8 0.133 213 0.067 120 a b (c) d 图7水冷条件下高硅钢铸锭在不同时刻的凝圆过程结果(a)和温度场（化~d).(a)99.22s:(b)204.23s:(c)274.23s:(d)339.23s Fig.7 Simulated solidification process (a)and temperature field(b-d)results of the high-silicon steel casting under water cooling condition at dif- ferent time:(a)99.22s:(b)204.23s:(c)274.23s:(d)339.23s 1800 1800 1700 1600 6 2 1600 P4 4 1400 1500 1200 400 1000 1300 00 600 1200 100200300400 500 400 600 0 10203040506070 时间s 时间s 图8空冷()及水冷(b)条件下高硅钢铸锭不同位置(P1~P4)的冷却曲线模拟结果 Fig.8 Simulated cooling curves at different locations (PI-P4 )in high-silicon steel casting under air cooling (a)and water cooling condition (b) 和())基本一致.整体对比两种冷却条件下高硅钢 相区较宽，呈“过渡式”凝固：而水冷条件下高硅钢凝 温度场和凝固过程的模拟结果可以发现，空冷条件下 固过程中温度梯度始终很大，两相区范围较为狭狭，呈 的温度场更均匀，温度梯度随着凝固进行逐渐减小，两 “分层式”凝固.空冷条件下宽阔的两相区、均匀的温

工程科学学报，第 39 卷，第 3 期 图 6 空冷条件下高硅钢铸锭不同时刻的凝固过程结果 ( a) 和温度场 ( b ～ d) ． ( a) 99. 22 s; ( b) 204. 23 s; ( c) 274. 23 s; ( d) 339. 23 s Fig． 6 Simulated solidification process ( a) and temperature field ( b--d) results of the high-silicon steel casting under air cooling condition at differ￾ent time: ( a) 99. 22 s; ( b) 204. 23 s; ( c) 274. 23 s; ( d) 339. 23 s 图 7 水冷条件下高硅钢铸锭在不同时刻的凝固过程结果 ( a) 和温度场 ( b ～ d) . ( a) 99. 22 s; ( b) 204. 23 s; ( c) 274. 23 s; ( d) 339. 23 s Fig． 7 Simulated solidification process ( a) and temperature field ( b--d) results of the high-silicon steel casting under water cooling condition at dif￾ferent time: ( a) 99. 22 s; ( b) 204. 23 s; ( c) 274. 23 s; ( d) 339. 23 s 图 8 空冷 ( a) 及水冷 ( b) 条件下高硅钢铸锭不同位置 ( P1 ～ P4) 的冷却曲线模拟结果 Fig． 8 Simulated cooling curves at different locations ( P1--P4 ) in high-silicon steel casting under air cooling ( a) and water cooling condition ( b) 和( e) ) 基本一致． 整体对比两种冷却条件下高硅钢 温度场和凝固过程的模拟结果可以发现，空冷条件下 的温度场更均匀，温度梯度随着凝固进行逐渐减小，两 相区较宽，呈“过渡式”凝固; 而水冷条件下高硅钢凝 固过程中温度梯度始终很大，两相区范围较为狭狭，呈 “分层式”凝固． 空冷条件下宽阔的两相区、均匀的温 · 463 ·

宋炜等：基于实验与3D-CAFE法的高硅钢铸锭凝固行为 ·365 度场有利于扩大等轴晶区，而水冷条件下较大的温度 域（图9(a)所示）：随着凝固的进行，铸锭边部的凝固 梯度使得热流传导加快且方向性增强，有利于柱状晶 坯壳逐渐向中心推进，回旋区域逐渐缩小并向铸锭上 的生长 方移动，最后消失（如图9(b)、()所示）.对比两种 3.2.2流场 冷却条件下的流场模拟结果可以发现，在凝固过程中 流体流动在凝固过程中至关重要，它影响金属液 空冷流场更均匀、更稳定，水冷流场的流动规律性并不 中溶质元素分布，决定铸锭的宏观偏析水平，还可以促 明显，这主要是由于水冷铸锭存在较大的温度梯度，热 进金属液中夹杂物上浮，改善铸锭的洁净度网.图9 对流强烈，流体在流动过程中流股互相干扰，导致流场 和图10分别为空冷和水冷条件下高硅钢凝固过程中 不稳定.另外，在凝固末期，空冷铸锭的冒口下方存在 不同时刻的流场模拟结果.从图中可以看出，凝固开 一个明显的抽吸区域（图9(）所示)，这种现象可以 始时铸锭边部的金属液密度随着温度的降低而增加， 对后凝固的区域进行补缩，减轻或消除二次缩孔等缺 因此边部的金属液在重力的驱动下向下流动，而铸锭 陷，但同样也可能产生中心偏析：而水冷铸锭中未 中心的金属液密度相对较小，在浮力的作用下向上流 观察到这种现象（图10(）所示)，因此推测水冷铸锭 动，于是在铸锭心中两侧位置形成了两个流动回旋区 的冒口下方可能存在微观缩孔或疏松 流速mm·s 流速(mm·s 流速/(mm·s 流速/(mm· 84231 6.695 0.002 83.70 4.812 1.013 0.00568 78.12 4.491 0.945 0.00530 72.54 4.170 0.878 0.00492 66.96 3.850 0.810 0.00454 61.38 529 0.743 0.00417 55.80 3.208 0.675 0.00379 50.22 2.887 0.608 0.00341 44.64 2.566 0.540 0.00303 39.06 2.246 0.473 0.00265 1.925 0.405 0.00227 1.604 0.338 0.00189 2232 1.283 0.270 0.00151 16.74 0.962 0.203 0.00114 1L.16 0.642 0.135 0.00076 5.58 0.321 0.068 0.00038 0 (b) (d) 图9空冷条件下高硅钢铸锭在不同时刻的流场分布.(a)6.75s:(b)72.22s:(c)154.63s:(d)331.23s Fig.9 Simulated flow field results of the high-silicon steel casting at different time under air cooling condition:(a)6.75 s:(b)72.22 s:(c) 154.63s:(d)331.23 流速/mm 流速(mm·g 流速/(mm·s) 流速/mm*s- 36.255 44.294 18.020 0.850 16.00 37.0 14.00 2.200 14.93 34.53 13.07 2.053 13.87 32.07 12.13 1.907 12.80 29W 183? 1.760 11.73 27.13 1.613 10.67 9.33 1.467 9.60 22.20 840 1.320 85 19,73 747 1.173 47 17.27 6.53 1.027 14.80 5.60 0.880 12.33 0.733 427 0.587 2.80 0.440 2.13 1.87 0.293 1.07 2.47 0.93 0.147 (b) ⑨y (d) 图10水冷条件下高硅钢铸锭在不同时刻的流场分布.(a)6.75s:(b)72.22s:(c)154.63s:(d)331.23s Fig.10 Simulated flow field results of the high-silicon steel casting at different time under water cooling condition:(a)6.75s:(b)72.22s:(c) 154.63s:(d)331.23s

宋 炜等: 基于实验与 3D--CAFE 法的高硅钢铸锭凝固行为 度场有利于扩大等轴晶区，而水冷条件下较大的温度 梯度使得热流传导加快且方向性增强，有利于柱状晶 的生长． 3. 2. 2 流场 流体流动在凝固过程中至关重要，它影响金属液 中溶质元素分布，决定铸锭的宏观偏析水平，还可以促 进金属液中夹杂物上浮，改善铸锭的洁净度［26］． 图 9 和图 10 分别为空冷和水冷条件下高硅钢凝固过程中 不同时刻的流场模拟结果． 从图中可以看出，凝固开 始时铸锭边部的金属液密度随着温度的降低而增加， 因此边部的金属液在重力的驱动下向下流动，而铸锭 中心的金属液密度相对较小，在浮力的作用下向上流 动，于是在铸锭心中两侧位置形成了两个流动回旋区 域 ( 图 9( a) 所示) ; 随着凝固的进行，铸锭边部的凝固 坯壳逐渐向中心推进，回旋区域逐渐缩小并向铸锭上 方移动，最后消失 ( 如图 9( b) 、( c) 所示) ． 对比两种 冷却条件下的流场模拟结果可以发现，在凝固过程中 空冷流场更均匀、更稳定，水冷流场的流动规律性并不 明显，这主要是由于水冷铸锭存在较大的温度梯度，热 对流强烈，流体在流动过程中流股互相干扰，导致流场 不稳定． 另外，在凝固末期，空冷铸锭的冒口下方存在 一个明显的抽吸区域 ( 图 9( d) 所示) ，这种现象可以 对后凝固的区域进行补缩，减轻或消除二次缩孔等缺 陷，但同样也可能产生中心偏析［26］; 而水冷铸锭中未 观察到这种现象 ( 图 10( d) 所示) ，因此推测水冷铸锭 的冒口下方可能存在微观缩孔或疏松． 图 9 空冷条件下高硅钢铸锭在不同时刻的流场分布 . ( a) 6. 75 s; ( b) 72. 22 s; ( c) 154. 63 s; ( d) 331. 23 s Fig． 9 Simulated flow field results of the high-silicon steel casting at different time under air cooling condition: ( a) 6. 75 s; ( b) 72. 22 s; ( c) 154. 63 s; ( d) 331. 23 s 图 10 水冷条件下高硅钢铸锭在不同时刻的流场分布 . ( a) 6. 75 s; ( b) 72. 22 s; ( c) 154. 63 s; ( d) 331. 23 s Fig． 10 Simulated flow field results of the high-silicon steel casting at different time under water cooling condition: ( a) 6. 75 s; ( b) 72. 22 s; ( c) 154. 63 s; ( d) 331. 23 s · 563 ·

·366 工程科学学报，第39卷，第3期 图11为两种不同冷却条件下高硅钢铸锭不同 成较强的热对流，所以流速相对其他位置较高.此 位置(P1~P4)的流速随时间的变化曲线，从图11 外，空冷铸锭在凝固前期流速曲线出现波动现象，随 中可以看出，无论是空冷条件还是水冷条件，流速变 着凝固进行波动越来越小逐渐趋于平稳，这是因为 化规律大致都是随着凝固进行先增大后减小，空冷 空冷冷速较小，凝固速度慢，凝固初期释放的结晶潜 铸锭中P3处流速最大、P1处次之：而水冷铸锭中P4 热不断地影响温度场变化进而影响流速变化.对比 处和P1处流速明显高于其他位置的流速，冒口下方 两种冷却条件下的流速变化曲线可以发现，空冷流 流速最大，这主要是因为空冷铸锭冷速慢、流场较为 场整体更稳定，这种稳定的流动状态可以促进柱状 均匀，而P1和P3分别处于流动回旋区的上、下涡流 晶向等轴晶转变(CET),有利于中心等轴晶形 线上，所以流速较高：而水冷铸锭边部位置冷速较 成叨，而水冷铸锭凝固时间较短，温度梯度始终较 快，中心及冒口处熔体温度较高，与附近低温区域形 大，导致流场始终不稳定 55 (a) P1 40b 50 P P2 35 P3 40 P4 3 25 25 10 25 50 75100125150175200225250 05 1015202530354045505560 时间/s 时间/s 图11空冷(a)及水冷(b)条件下高硅钢铸锭不同位置(PI,2,P3和P4)的流速变化曲线 Fig.11 Velocity of flow curves of high-silicon steel casting at different locations (Pl,P2,P3 and P4)under air cooling (a)and water cooling (b) respectively 3.2.3凝固组织 达，晶粒粗大，水冷铸锭几乎由近90%的柱状晶构成， (1)冷却条件的影响 且晶粒平均尺寸和平均面积都要大于空冷铸锭.这主 CAFE模型计算所得到空冷和水冷条件下的高硅 要是由于空冷铸锭凝固前沿的温度梯度较小，两相区 钢凝固组织如图12(b)和(d),从图12中可以发现， 较宽（图6所示），更利于枝晶生长，凝固前沿流体流 CAE模拟结果无论是组织形貌或是晶粒尺寸都与实 动将纤细的枝晶尖端打碎并带到熔体中心成为等轴晶 验结果基本一致.进一步对晶粒尺寸等信息进行统计 的形核核心：而水冷铸锭凝固前沿温度梯度较大，两相 分析，结果见表3.结合表3中的统计结果及图12(b) 区狭窄（图7所示），更有利于柱状晶生长，因此水冷 和()中的组织模拟结果可以发现，高硅钢柱状晶发 铸锭中柱状晶比例高 图12高硅钢铸态组织()实验结果，空冷：(b)模拟结果，空冷：(c)实验结果，水冷：(d)模拟结果，水冷 Fig.12 Structures of high-silicon steel casting:(a)experimental results,air cooling:(b)simulated results,air cooling:(c)experimental results, water cooling:(d)simulated results,water cooling (2)过热度的影响 铸锭质量.利用空冷模型，通过改变初始温度来研究 过热度是浇注过程中重要的工艺参数，直接影响 不同过热度对凝固组织的影响，模拟结果如图13及表4

工程科学学报，第 39 卷，第 3 期 图 11 为两种不同冷却条件下高硅钢铸锭不同 位置 ( P1 ～ P4) 的流速随时间的变化曲线，从图 11 中可以看出，无论是空冷条件还是水冷条件，流速变 化规律大致都是随着凝固进行先增大后减小，空冷 铸锭中 P3 处流速最大、P1 处次之; 而水冷铸锭中 P4 处和 P1 处流速明显高于其他位置的流速，冒口下方 流速最大． 这主要是因为空冷铸锭冷速慢、流场较为 均匀，而 P1 和 P3 分别处于流动回旋区的上、下涡流 线上，所以 流 速 较 高; 而 水 冷 铸 锭 边 部 位 置 冷 速 较 快，中心及冒口处熔体温度较高，与附近低温区域形 成较强的热对流，所 以 流 速 相 对 其 他 位 置 较 高． 此 外，空冷铸锭在凝固前期流速曲线出现波动现象，随 着凝固进行波动越来越小逐渐趋于平稳，这是因为 空冷冷速较小，凝固速度慢，凝固初期释放的结晶潜 热不断地影响温度场变化进而影响流速变化． 对比 两种冷却条件下的流速变化曲线可以发现，空冷流 场整体更稳定，这种稳定的流动状态可以促进柱状 晶向 等 轴 晶 转 变 ( CET) ，有 利 于 中 心 等 轴 晶 形 成［27］，而水冷 铸 锭 凝 固 时 间 较 短，温 度 梯 度 始 终 较 大，导致流场始终不稳定． 图 11 空冷( a) 及水冷( b) 条件下高硅钢铸锭不同位置 ( P1，P2，P3 和 P4) 的流速变化曲线 Fig． 11 Velocity of flow curves of high-silicon steel casting at different locations ( P1，P2，P3 and P4) under air cooling ( a) and water cooling ( b) ， respectively 3. 2. 3 凝固组织 ( 1) 冷却条件的影响． CAFE 模型计算所得到空冷和水冷条件下的高硅 钢凝固组织如图 12 ( b) 和( d) ，从图 12 中可以发现， CAFE 模拟结果无论是组织形貌或是晶粒尺寸都与实 验结果基本一致． 进一步对晶粒尺寸等信息进行统计 分析，结果见表 3． 结合表 3 中的统计结果及图 12( b) 和( d) 中的组织模拟结果可以发现，高硅钢柱状晶发 达，晶粒粗大，水冷铸锭几乎由近 90% 的柱状晶构成， 且晶粒平均尺寸和平均面积都要大于空冷铸锭． 这主 要是由于空冷铸锭凝固前沿的温度梯度较小，两相区 较宽( 图 6 所示) ，更利于枝晶生长，凝固前沿流体流 动将纤细的枝晶尖端打碎并带到熔体中心成为等轴晶 的形核核心; 而水冷铸锭凝固前沿温度梯度较大，两相 区狭窄( 图 7 所示) ，更有利于柱状晶生长，因此水冷 铸锭中柱状晶比例高． 图 12 高硅钢铸态组织 . ( a) 实验结果，空冷; ( b) 模拟结果，空冷; ( c) 实验结果，水冷; ( d) 模拟结果，水冷 Fig． 12 Structures of high-silicon steel casting: ( a) experimental results，air cooling; ( b) simulated results，air cooling; ( c) experimental results， water cooling; ( d) simulated results，water cooling ( 2) 过热度的影响． 过热度是浇注过程中重要的工艺参数，直接影响 铸锭质量． 利用空冷模型，通过改变初始温度来研究 不同过热度对凝固组织的影响，模拟结果如图 13 及表 4 · 663 ·

宋炜等：基于实验与3D-CAFE法的高硅钢铸锭凝固行为 367 表3不同冷却条件下凝固组织数据统计结果 K时，铸锭中等轴晶比例从49.53%增加到69.66%， Table 3 Statistics results of solidification structure under different cool- 平均晶粒直径由1.802mm减小到1.571mm,平均晶粒 ing conditions 面积随过热度降低而有所减小，而晶粒数随过热度降 冷却 等轴晶 计算截面 平均晶粒 平均品粒 低有所增加.这主要是由于高过热度导致凝固前沿温 方式 比例/% 品粒数 直径/mm 面积/mm2 度梯度增加，促进了柱状晶生长：同时较高的过热度也 空冷 49.53 2169 1.802 4.807 会提高熔体中异质核心熔解消失的概率，抑制了等轴 水冷 11.12 1303 4.508 9.443 晶形成：因此，低过热度浇铸是提高等轴晶率、细化晶 所示.由模拟结果可以看出，过热度从52K降低至20 粒的一种有效办法 图13不同过热度下铸锭凝固组织模拟结果.(a)52K:(b)40K:(c)30K:(d)20K Fig.13 Simulation results of solidification structure under the different superheats:(a)52 K:(b)40 K:(c)30 K:(d)20 K 表4不同过热度下凝固组织晶粒数据统计结果 (4)CAFE法模拟得到的凝固组织无论是组织形 Table 4 Statistics results of solidification structure under different su- 貌还是晶粒尺寸都与实验结果基本一致；空冷铸锭中 perheats 无论是等轴晶比例还是晶粒尺寸等凝固特征都优于水 过热度/ 等轴品 计算截面 平均品粒 平均品粒 冷铸锭；随着过热度的降低，高硅钢铸锭中心等轴晶率 K 比例/% 品粒数 直径/mm 面积/mm2 提高，晶粒数增加，晶粒尺寸变得更细小 20 69.66 2325 1.571 4.312 吃 62.56 2244 1.653 4.525 参考文献 40 60.35 2207 1.737 4.692 [He ZZ,Zhao Y,Luo H W.Electrical Sheet.Beijing:Metallurgy 49.53 2169 1.802 4.807 Industry Press,2012 (何忠志，赵宇，罗海文.电工钢.北京：治金工业出版社， 2012) 4结论 2] Haiji H,Okada K,Hiratani T,et al.Magnetic properties and workability of 6.5%Si steel sheet.J MagnMagn Mater,1996, (1)高硅钢铸态组织主要由粗大的柱状晶构成， 160:109 水冷铸锭中柱状晶比例高达90%，空冷铸锭中心有少 3]Fang X S,Liang Y F,Ye F,et al.Effect of rolling reduction on 量等轴晶区：柱状晶晶粒由发达的一次枝晶构成，且具 the texture of 6.5wt%Si electric steel during warm rolling.J 有明显的方向性：等轴晶晶粒由二次、三次枝晶构成， Funct Mater,2012,43(24):3346 取向随机分布 (房现石，梁永锋，叶丰，等.压下率对6.5%Si电工钢温轧 版织构的影响规律.功能材料，2012,43(24)：3346) (2)在高硅钢凝固过程中，空冷铸锭的温度场更 4] Yang JS,Xie JX,Zhou C.Preparation technology and prospeet 均匀，两相区更宽阔，凝固从铸锭角部开始沿铸壁逐渐 of 6.5%Si steel.J Funct Mater,2003,34(3)244 向中心推进，呈“过渡式”凝固；而水冷铸锭温度梯度 (杨劲松，谢建新，周成.6.5%Si高硅钢的制备工艺及发展 始终较大，两相区狭小，呈“分层式”凝固 前景.功能材料，2003,34(3)：244) (3)在凝固前期空冷铸锭的中心两侧位置存在两 5]Zheng X,Yan B.Properties and preparation techniques of Fe- 个较大的流动回旋区，凝固末期在铸锭冒口下方出现 6.5%Si high silicon steel.Mater Rer,2012,26(19)392 了抽吸补缩现象：而水冷铸件在凝固过程中流股相互 (郑鑫，严彪.Fe6.5%Si高硅钢的性能及制备技术.材料导 报，2012,26(19)：392) 干扰，流场不稳定，凝固末期未观察到补缩现象.空冷 6]Lin J P,Ye F,Chen G L,et al.Fabrication technology,micro- 铸锭在凝固开始时流速波动较大，随凝固进行逐渐平 structures and properties of Fe-.5w%Si alloy sheets by cold 稳，而水冷铸锭在凝固过程中流速始终都在波动 rolling.Front Sci,2007(2)13

宋 炜等: 基于实验与 3D--CAFE 法的高硅钢铸锭凝固行为 表 3 不同冷却条件下凝固组织数据统计结果 Table 3 Statistics results of solidification structure under different cool￾ing conditions 冷却 方式 等轴晶 比例/% 计算截面 晶粒数 平均晶粒 直径/mm 平均晶粒 面积/mm2 空冷 49. 53 2169 1. 802 4. 807 水冷 11. 12 1303 4. 508 9. 443 所示． 由模拟结果可以看出，过热度从 52 K 降低至 20 K 时，铸锭中等轴晶比例从 49. 53% 增加到 69. 66% ， 平均晶粒直径由1. 802 mm 减小到1. 571 mm，平均晶粒 面积随过热度降低而有所减小，而晶粒数随过热度降 低有所增加． 这主要是由于高过热度导致凝固前沿温 度梯度增加，促进了柱状晶生长; 同时较高的过热度也 会提高熔体中异质核心熔解消失的概率，抑制了等轴 晶形成; 因此，低过热度浇铸是提高等轴晶率、细化晶 粒的一种有效办法． 图 13 不同过热度下铸锭凝固组织模拟结果 . ( a) 52 K; ( b) 40 K; ( c) 30 K; ( d) 20 K Fig． 13 Simulation results of solidification structure under the different superheats: ( a) 52 K; ( b) 40 K; ( c) 30 K; ( d) 20 K 表 4 不同过热度下凝固组织晶粒数据统计结果 Table 4 Statistics results of solidification structure under different su￾perheats 过热度/ K 等轴晶 比例/% 计算截面 晶粒数 平均晶粒 直径/mm 平均晶粒 面积/mm2 20 69. 66 2325 1. 571 4. 312 30 62. 56 2244 1. 653 4. 525 40 60. 35 2207 1. 737 4. 692 52 49. 53 2169 1. 802 4. 807 4 结论 ( 1) 高硅钢铸态组织主要由粗大的柱状晶构成， 水冷铸锭中柱状晶比例高达 90% ，空冷铸锭中心有少 量等轴晶区; 柱状晶晶粒由发达的一次枝晶构成，且具 有明显的方向性; 等轴晶晶粒由二次、三次枝晶构成， 取向随机分布． ( 2) 在高硅钢凝固过程中，空冷铸锭的温度场更 均匀，两相区更宽阔，凝固从铸锭角部开始沿铸壁逐渐 向中心推进，呈“过渡式”凝固; 而水冷铸锭温度梯度 始终较大，两相区狭小，呈“分层式”凝固． ( 3) 在凝固前期空冷铸锭的中心两侧位置存在两 个较大的流动回旋区，凝固末期在铸锭冒口下方出现 了抽吸补缩现象; 而水冷铸件在凝固过程中流股相互 干扰，流场不稳定，凝固末期未观察到补缩现象． 空冷 铸锭在凝固开始时流速波动较大，随凝固进行逐渐平 稳，而水冷铸锭在凝固过程中流速始终都在波动． ( 4) CAFE 法模拟得到的凝固组织无论是组织形 貌还是晶粒尺寸都与实验结果基本一致; 空冷铸锭中 无论是等轴晶比例还是晶粒尺寸等凝固特征都优于水 冷铸锭; 随着过热度的降低，高硅钢铸锭中心等轴晶率 提高，晶粒数增加，晶粒尺寸变得更细小． 参 考 文 献 ［1］ He Z Z，Zhao Y，Luo H W． Electrical Sheet． Beijing: Metallurgy Industry Press，2012 ( 何忠志，赵宇，罗海文． 电工钢． 北京: 冶金工业出版社， 2012) ［2］ Haiji H，Okada K，Hiratani T，et al． Magnetic properties and workability of 6. 5% Si steel sheet． J Magn Magn Mater，1996， 160: 109 ［3］ Fang X S，Liang Y F，Ye F，et al． Effect of rolling reduction on the texture of 6. 5wt% Si electric steel during warm rolling． J Funct Mater，2012，43( 24) : 3346 ( 房现石，梁永锋，叶丰，等． 压下率对 6. 5% Si 电工钢温轧 板织构的影响规律． 功能材料，2012，43( 24) : 3346) ［4］ Yang J S，Xie J X，Zhou C． Preparation technology and prospect of 6. 5% Si steel． J Funct Mater，2003，34( 3) : 244 ( 杨劲松，谢建新，周成． 6. 5% Si 高硅钢的制备工艺及发展 前景． 功能材料，2003，34( 3) : 244) ［5］ Zheng X，Yan B． Properties and preparation techniques of Fe-- 6. 5% Si high silicon steel． Mater Ｒev，2012，26( 19) : 392 ( 郑鑫，严彪． Fe--6. 5% Si 高硅钢的性能及制备技术． 材料导 报，2012，26( 19) : 392) ［6］ Lin J P，Ye F，Chen G L，et al． Fabrication technology，micro￾structures and properties of Fe--6. 5wt% Si alloy sheets by cold rolling． Front Sci，2007( 2) : 13 · 763 ·

·368· 工程科学学报，第39卷，第3期 (林均品，叶丰，陈国良，等.6.5wt%Si高硅钢冷轧薄板制备 Res1nt,2011,82(10):1173 工艺、结构和性能.前沿科学，2007(2)：13) [17]Wang J L,Wang F M,Zhao YY,et al.Numerical simulation of Spittle JA.Brown SG R.Computer simulation of the effects of 3D-microstructures in solidification processes based on the CAFE alloy variables on the grain structures of castings.Acta Metall, method.Int J Miner Metall Mater,2009,16(6):640 1989,37(7):1803 [18]Hou Z B,Jiang F,Cheng GG.Solidification structure and com- [8]Wang S L.Sekerka R F,Wheeler AA,et al.Thermodynamical- pactness degree of central equiaxed grain zone in continuous cast- ly-consistent phase-field models for solidification.Phys D,1993, ing billet using cellular automaton-finite element method.IS/ 69(12):189 ,2012,52(7):1301 9]Natsume Y,Ohsasa K,Narita T.Phase-field simulation of transi- [19]Kattner U R.The thermodynamic modeling of multicomponen ent liquid phase bonding process of Ni using Ni-P binary filler phase equilibria.JOM,1997,49(12)14 metal.Mater Trans,2003,44(5)819 [20] Thevoz P.Desbiolles JL Rappaz M.Modeling of equiaxed mi- [10]Zhao Y Z,Shi Y W,Shi L F.Current research status of comput- crostructure formation in casting.Metall Trans A.1989.20(2): er simulation of casting solidification structure.Met Form 311 Technol,2002,20(6):53 21] Kurz W,Giovanola B,Trivedi R.Theory of microstructural de- (赵玉珍，史耀武，史立丰.铸件凝固过程组织计算机模拟 velopment during rapid solidification.Acta Metall,1986,34 研究动态.金属成型工艺，2002,20(6)：53) (5):823 [11]Rappaz M,Gandin C A.Probabilistic modelling of microstructure 22] Gandin C A,Rappaz M.A coupled finite element-cellular au- formation in solidification processes.Acta Metall Mater,1993, tomaton model for the prediction of dendritic grain structures in 41(2):345 solidification processes.Acta Metall Mater,1994,42(7):2233 [2]Gandin CA,Rappaz M.A 3D cellular automaton algorithm for 3]Kovac F,Stoyka V,Petryshynets I.Strain-induced grain growth the prediction of dendritic grain growth.Acta Mater,1997,45 in non-oriented electrical steels.J MagnMagn Mater,2008, (5):2187 320(20):c627 [13]Gandin C A,Desbiolles J L,Rappaz M,et al.A three-dimen- 4]Zhang N.Yang P,Mao W M.Influence of columnar grains on sional cellular automation-finite element model for the prediction the recrystallization texture evolution in Fe-3%Si electrical of solidification grain structures.Metall Mater Trans A,1999,30 steel.Acta Metall Sin,2012,48(3):307 (12):3153 (张宁，杨平，毛卫民.柱状品对F©-3%Si电工钢再结晶织 [14]Pang R P,Wang F M,Zhang G Q,et al.Study of solidification 构演变规律的影响.金属学报，2012,48(3)：307) thermal parameters of 430 ferrite stainless steel based on 3D-CA- 25] Song W,Zhang J M,Liu Y,et al.Numerical simulation of so- FE method.Acta Metall Sin,2013,49(10):1234 lidification structure of 6.5 wt%Si steel ingot slab.lronmaking (庞瑞朋，王福明，张国庆，等.基于3D-CAFE法对430铁 Steelmaking,2015,42(9）:656 素体不锈钢凝固热参数的研究.金属学报，2013,49(10)： [26 Pickering EJ.Macrosegregation in steel ingots:the applicability 1234) of modelling and characterisation techniques.IS//Int,2013,53 [15]Lno YZ,Zhang J M,Wei X D,et al.Numerical simulation of (6):935 solidification structure of high carbon SWRH77B billet based on 27]Patil P,Nalawade R,Balachandran G,et al.Analysis of solidi- the CAFE method.Ironmaking Steelmaking,2012,39(1):26 fication behaviour of low alloy steel ingot casting-simulation and [6]Jing CL,Wang X H,Jiang M.Study on solidification structure experimental validation.fronmaking Steelmaking,2015,42 of wheel steel round billet using FE-CA coupling model.Steel (7):512

工程科学学报，第 39 卷，第 3 期 ( 林均品，叶丰，陈国良，等． 6. 5wt% Si 高硅钢冷轧薄板制备 工艺、结构和性能． 前沿科学，2007( 2) : 13) ［7］ Spittle J A，Brown S G Ｒ． Computer simulation of the effects of alloy variables on the grain structures of castings． Acta Metall， 1989，37( 7) : 1803 ［8］ Wang S L，Sekerka Ｒ F，Wheeler A A，et al． Thermodynamical￾ly-consistent phase-field models for solidification． Phys D，1993， 69( 1-2) : 189 ［9］ Natsume Y，Ohsasa K，Narita T． Phase-field simulation of transi￾ent liquid phase bonding process of Ni using Ni--P binary filler metal． Mater Trans，2003，44( 5) : 819 ［10］ Zhao Y Z，Shi Y W，Shi L F． Current research status of comput￾er simulation of casting solidification structure． Met Form Technol，2002，20( 6) : 53 ( 赵玉珍，史耀武，史立丰． 铸件凝固过程组织计算机模拟 研究动态． 金属成型工艺，2002，20( 6) : 53) ［11］ Ｒappaz M，Gandin C A． Probabilistic modelling of microstructure formation in solidification processes． Acta Metall Mater，1993， 41( 2) : 345 ［12］ Gandin C A，Ｒappaz M． A 3D cellular automaton algorithm for the prediction of dendritic grain growth． Acta Mater，1997，45 ( 5) : 2187 ［13］ Gandin C A，Desbiolles J L，Ｒappaz M，et al． A three-dimen￾sional cellular automation-finite element model for the prediction of solidification grain structures． Metall Mater Trans A，1999，30 ( 12) : 3153 ［14］ Pang Ｒ P，Wang F M，Zhang G Q，et al． Study of solidification thermal parameters of 430 ferrite stainless steel based on 3D--CA￾FE method． Acta Metall Sin，2013，49( 10) : 1234 ( 庞瑞朋，王福明，张国庆，等． 基于 3D--CAFE 法对 430 铁 素体不锈钢凝固热参数的研究． 金属学报，2013，49( 10) : 1234) ［15］ Luo Y Z，Zhang J M，Wei X D，et al． Numerical simulation of solidification structure of high carbon SWＲH77B billet based on the CAFE method． Ironmaking Steelmaking，2012，39( 1) : 26 ［16］ Jing C L，Wang X H，Jiang M． Study on solidification structure of wheel steel round billet using FE--CA coupling model． Steel Ｒes Int，2011，82( 10) : 1173 ［17］ Wang J L，Wang F M，Zhao Y Y，et al． Numerical simulation of 3D-microstructures in solidification processes based on the CAFE method． Int J Miner Metall Mater，2009，16( 6) : 640 ［18］ Hou Z B，Jiang F，Cheng G G． Solidification structure and com￾pactness degree of central equiaxed grain zone in continuous cast￾ing billet using cellular automaton-finite element method． ISIJ Int，2012，52( 7) : 1301 ［19］ Kattner U Ｒ． The thermodynamic modeling of multicomponent phase equilibria． JOM，1997，49( 12) : 14 ［20］ Thevoz P，Desbiolles J L，Ｒappaz M． Modeling of equiaxed mi￾crostructure formation in casting． Metall Trans A，1989，20( 2) : 311 ［21］ Kurz W，Giovanola B，Trivedi Ｒ． Theory of microstructural de￾velopment during rapid solidification． Acta Metall，1986，34 ( 5) : 823 ［22］ Gandin C A，Ｒappaz M． A coupled finite element-cellular au￾tomaton model for the prediction of dendritic grain structures in solidification processes． Acta Metall Mater，1994，42( 7) : 2233 ［23］ Kovac F，Stoyka V，Petryshynets I． Strain-induced grain growth in non-oriented electrical steels． J Magn Magn Mater，2008， 320( 20) : e627 ［24］ Zhang N，Yang P，Mao W M． Influence of columnar grains on the recrystallization texture evolution in Fe-- 3% Si electrical steel． Acta Metall Sin，2012，48( 3) : 307 ( 张宁，杨平，毛卫民． 柱状晶对 Fe--3% Si 电工钢再结晶织 构演变规律的影响． 金属学报，2012，48( 3) : 307) ［25］ Song W，Zhang J M，Liu Y，et al． Numerical simulation of so￾lidification structure of 6. 5 wt-% Si steel ingot slab． Ironmaking Steelmaking，2015，42( 9) : 656 ［26］ Pickering E J． Macrosegregation in steel ingots: the applicability of modelling and characterisation techniques． ISIJ Int，2013，53 ( 6) : 935 ［27］ Patil P，Nalawade Ｒ，Balachandran G，et al． Analysis of solidi￾fication behaviour of low alloy steel ingot casting-simulation and experimental validation． Ironmaking Steelmaking， 2015， 42 ( 7) : 512 · 863 ·