彭其春,彭 勝,楊進(jìn)玲,劉紅軍,曹建新,楊秀芝
(1.武漢科技大學(xué)鋼鐵冶金及資源利用省部共建教育部重點(diǎn)實(shí)驗(yàn)室,湖北 武漢,4300812.湖南華菱漣源鋼鐵集團(tuán)有限公司,湖南 婁底,417009)
應(yīng)用二維非穩(wěn)態(tài)傳熱數(shù)學(xué)模型建立板坯連鑄凝固傳熱模型對(duì)鑄坯溫度場(chǎng)進(jìn)行仿真研究是當(dāng)前連鑄技術(shù)研究的重要手段。有關(guān)該類(lèi)模型用于連鑄結(jié)晶器冷卻、二冷配水設(shè)計(jì)分析等取得了很好的效果[1-6]。生產(chǎn)實(shí)際表明,結(jié)晶器寬、窄面的冷卻水量以及結(jié)晶器等設(shè)備的運(yùn)行狀況均可影響連鑄坯角部裂紋的形成[5]。
本文結(jié)合湖南華菱漣源鋼鐵集團(tuán)有限公司(以下簡(jiǎn)稱漣鋼)210轉(zhuǎn)爐廠現(xiàn)場(chǎng)實(shí)際,建立板坯凝固傳熱模型,利用該模型對(duì)連鑄結(jié)晶器冷卻(一次冷卻)進(jìn)行模擬分析,提出了使鑄坯角部裂紋發(fā)生率最小的結(jié)晶器冷卻水控制方案。
基本假設(shè)如下:①鑄坯傳熱為二維傳熱;②鑄坯厚度和寬度方向的傳熱為軸對(duì)稱;③注入結(jié)晶器的鋼液溫度恒定;④忽略鑄坯與夾輥的傳熱;⑤引入等效導(dǎo)熱系數(shù)來(lái)考慮鑄坯液芯的對(duì)流傳熱;⑥忽略鑄坯密度的變化,鋼的導(dǎo)熱系數(shù)是溫度的函數(shù);⑦凝固潛熱的釋放按等效比熱容處理。
根據(jù)微元體內(nèi)能量守衡定律,板坯凝固傳熱微分方程為

式中:T為鑄坯中心鋼液溫度,℃;ρ為鋼的密度,kg·m-3;ce為等效比熱容,J/(kg·K);λ為導(dǎo)熱系數(shù),W/(m·K)。
板坯從結(jié)晶器內(nèi)鋼水彎月面向下以一定速度移動(dòng),熱量從鑄坯中心向表面?zhèn)鬟f,其熱量的傳遞受鑄坯表面邊界條件和金屬熱物理性能的影響。以板坯左下角為原點(diǎn)、寬度方向?yàn)閤軸和厚度方向?yàn)閥軸建立坐標(biāo)系。板坯凝固傳熱方程計(jì)算區(qū)域如圖1所示。圖1中,A為鑄坯斷面寬度中心;C為鑄坯斷面厚度中心;L為鑄機(jī)冶金長(zhǎng)度。

圖1 凝固傳熱方程計(jì)算區(qū)域Fig.1 Computation region for solidification heat transfer equation
1.3.1 初始條件

式中:t為鋼水在結(jié)晶器內(nèi)停留時(shí)間,min;T0為鋼水澆注溫度,℃。
1.3.2 邊界條件
鑄坯中心

式中:q為熱流密度,W/m2;a、b為實(shí)驗(yàn)常數(shù)。
在靜止的水所冷卻的結(jié)晶器內(nèi)測(cè)定,熱流密度與鋼水在結(jié)晶器內(nèi)停留時(shí)間的關(guān)系[7]為

式中:c、d為實(shí)驗(yàn)常數(shù);l為結(jié)晶器有效長(zhǎng)度,m;V為鋼的拉速,m/min。
1.4.1 鋼的液、固相線溫度計(jì)算[7-9]

1.4.2 凝固潛熱Lf
用數(shù)學(xué)模型計(jì)算凝固潛熱Lf時(shí),可通過(guò)液固兩相等效比熱容求得,其計(jì)算式為

式中:ceff為液固兩相等效比熱容,J/(kg·K);c為比熱容,J/(kg·K);Tl、Ts分別為液相線、固相線溫度,℃。
1.4.3 導(dǎo)熱系數(shù)λ
采用激光熱導(dǎo)儀測(cè)定不同鋼種固相、液相的的導(dǎo)熱系數(shù)。液相區(qū)流動(dòng)鋼液導(dǎo)熱系數(shù)相當(dāng)于靜止鋼液導(dǎo)熱系數(shù)的4~6倍,預(yù)取5倍。液固兩相區(qū)導(dǎo)熱系數(shù)λe采用下式計(jì)算:

式中:λs、λl分別為固相區(qū)、液相區(qū)導(dǎo)熱系數(shù),W/(m·K)。
1.4.4 比熱容c

利用有限元編程來(lái)模擬結(jié)晶器內(nèi)溫度場(chǎng)的變化。以漣鋼210轉(zhuǎn)爐廠1#連鑄機(jī)組冷卻水寬面水量4200L/min、窄面水量400L/min為冷卻條件,研究結(jié)晶器內(nèi)寬、窄面的坯殼厚度變化以及角部溫度和中心溫度的變化。結(jié)晶器內(nèi)寬、窄面鑄坯坯殼厚度隨時(shí)間的變化分別如圖2、圖3所示。由圖中可以看出,拉速一定,結(jié)晶器內(nèi)鑄坯寬、窄面的坯殼厚度隨時(shí)間的推移逐漸增大。

圖2 寬面坯殼厚度隨時(shí)間的變化Fig.2 Variation of the thickness of broadside shell with time

圖3 窄面坯殼厚度隨時(shí)間的變化Fig.3 Variation of the thickness of narrow side shell with time
結(jié)晶器鑄坯坯殼厚度隨拉速的變化如圖4所示,結(jié)晶器鑄坯角部溫度隨時(shí)間的變化如圖5所示。從圖4中可看出,隨著拉速的增大,鑄坯出結(jié)晶器的寬、窄面坯殼厚度逐漸減小。這是因?yàn)殡S著拉速增大,結(jié)晶器內(nèi)的冷卻強(qiáng)度降低,從而減小了鑄坯出結(jié)晶器時(shí)的坯殼厚度。從圖5中可看出,初始階段,隨著時(shí)間的推移,鑄坯角部溫度迅速降低,在11~15s時(shí),鑄坯角部溫度出現(xiàn)了小量回升。

圖4 不同拉速下出結(jié)晶器鑄坯坯殼厚度變化Fig.4 Variation of the thickness of slab shell with casting speed

圖5 鑄坯角部溫度隨時(shí)間的變化Fig.5 Variation of slab temperature in corner portion with time
結(jié)晶器鑄坯寬、窄面中心溫度隨時(shí)間的變化分別如圖6、圖7所示。從圖中可看出,當(dāng)結(jié)晶器水量為4200/400L/min時(shí),鑄坯寬、窄面中心溫度隨時(shí)間推移迅速降低,在7s左右,降勢(shì)趨緩。

圖6 鑄坯窄面中心溫度隨時(shí)間的變化Fig.6 Variation of the core temperature of slab narrow side with time
基于上述模擬分析,以結(jié)晶器水量為4200/400L/min、典型拉速1.00m/min和實(shí)際生產(chǎn)中最大拉速1.15m/min為條件,模擬鑄坯出結(jié)晶器坯殼厚度如表1所示。
按照表1中模擬鑄坯出結(jié)晶器窄面坯殼厚度數(shù)據(jù),利用公式(7),可計(jì)算出理論出結(jié)晶器坯殼厚度如表2所示。

表1 模擬鑄坯出結(jié)晶器坯殼厚度Table1 Simulated slab shell thickness

式中:em為坯殼厚度,mm;t為凝固時(shí)間,min(t=H/V,H 為結(jié)晶器有效高度,m;V 為拉速,m/min);K 為凝固系數(shù),mm/min。
對(duì)比表1和表2可知,模擬出結(jié)晶器坯殼厚度僅高于理論出結(jié)晶器坯殼厚度約2mm,表明該廠結(jié)晶器弱冷卻控制較好,但仍有下調(diào)空間。
結(jié)合漣鋼201轉(zhuǎn)爐廠現(xiàn)場(chǎng)條件,在不同結(jié)晶器冷卻水量和拉速條件下,對(duì)結(jié)晶器內(nèi)鑄坯坯殼厚度隨結(jié)晶器有效高度的變化進(jìn)行模擬,結(jié)果分別如圖8、圖9所示。由圖8、圖9中可以看出,隨著結(jié)晶器有效高度的增大,鑄坯坯殼的厚度逐漸增大;隨著結(jié)晶器冷卻水水量減小,鑄坯坯殼厚度逐漸減小,這是由于結(jié)晶器冷卻水量減小,導(dǎo)致結(jié)晶器的冷卻強(qiáng)度下降,當(dāng)結(jié)晶器有效高度大于0.5m時(shí),這種趨勢(shì)更為明顯。

圖8 坯殼厚度隨結(jié)晶器有效高度的變化(拉速為1.00 m/min)Fig.8 Variation of the thickness of slab shell with the effective height of crystallizer(casting speed:1.00m/min)

圖9 坯殼厚度隨結(jié)晶器有效高度的變化(拉速為1.15 m/min)Fig.9 Variation of the thickness of slab shell with the effective height of crystallizer(casting speed:1.15m/min)
結(jié)晶器冷卻水量調(diào)整后的鑄坯出結(jié)晶器坯殼厚度如表3所示。由表3中可看出,拉速一定,鑄坯出結(jié)晶器寬、窄面坯殼厚度隨水量的減小逐漸減小。對(duì)比表2可知,當(dāng)結(jié)晶器水量為3600/340 L/min時(shí),鑄坯出結(jié)晶器坯殼厚度達(dá)到了理論下限值,即結(jié)晶器冷卻水量大于3600/340L/min時(shí),鑄坯出結(jié)晶器坯殼厚度理論上均能滿足要求。實(shí)際生產(chǎn)中發(fā)現(xiàn),當(dāng)結(jié)晶器冷卻水量從4200/400 L/min降至4000/380L/min時(shí),鑄坯角裂情況得到明顯改善;當(dāng)結(jié)晶器冷卻水量降至3800/360 L/min和3600/340L/min時(shí),角裂又有不同程度的加劇。原因在于結(jié)晶器水量過(guò)大引起坯殼冷卻不均,使得結(jié)晶器彎月面處的初生坯殼凝固不均勻,在鑄坯表面形成角裂;結(jié)晶器水量過(guò)低出結(jié)晶器的鑄坯強(qiáng)度相對(duì)較低,拉坯過(guò)程中易產(chǎn)生應(yīng)力集中,使坯殼角部出現(xiàn)細(xì)小的橫裂紋,其細(xì)小裂紋在二冷區(qū)擴(kuò)展,最終形成角裂。

表3 結(jié)晶器冷卻水量調(diào)整后的鑄坯出結(jié)晶器坯殼厚度Table3 Adjusted slab shell thickness

圖10 鑄坯各部位溫度隨結(jié)晶器有效高度的變化Fig.10 Variation of the temperature of slab shell with the effective height of crystallizer
不同水量條件下,鑄坯各部位溫度隨結(jié)晶器有效高度的變化如圖10所示。從圖10中可看出,鑄坯角部溫度低于寬面中心及窄面中心的溫度,這是因?yàn)榻遣渴嵌S傳熱,而寬面和窄面是一維傳熱,角部比寬面中心和窄面中心冷卻強(qiáng)度大。
基于以上分析可認(rèn)為,實(shí)際生產(chǎn)中結(jié)晶器冷卻水量為4000/380L/min時(shí),可有效降低鑄坯角裂的發(fā)生率。
(1)模擬出結(jié)晶器坯殼厚度高于理論出結(jié)晶器坯殼厚度約2mm,漣鋼210轉(zhuǎn)爐廠結(jié)晶器弱冷卻仍有下調(diào)空間。
(2)結(jié)晶器冷卻水量為4000/380L/min時(shí),鑄坯角裂可得到明顯的改善。
[1]王葉婷,孫猛.板坯連鑄表面溫度值二冷動(dòng)態(tài)控制[J].熱處理,2009,24(1):42-45.
[2]韓志偉,馮科,王勇,等.板坯連鑄二冷計(jì)算機(jī)仿真軟件的應(yīng)用[J].鋼鐵技術(shù),2007(1):23-26.
[3]謝建輝,岳學(xué)軍.石鋼1號(hào)連鑄機(jī)二冷配水設(shè)計(jì)分析[J].連鑄,2001(4):5-6.
[4]李朗.不同拉速下內(nèi)外復(fù)合冷卻結(jié)晶器內(nèi)鋼液流場(chǎng)的數(shù)值模擬[J].冶金叢刊,2008(177):1-3.
[5]譚飛.連鑄坯角部縱裂紋原因分析及對(duì)策[J].四川冶金,2006,28(4):20-23.
[6]張有慧,王宥宏,虞明香,等.碳鋼板坯連鑄凝固模擬[J].特殊鋼,2004,25(6):13-15.
[7]蔡開(kāi)科.連續(xù)鑄鋼[M].北京:冶金工業(yè)出版社,1990.
[8]馮科,徐楚韶,陳登福.方坯連鑄凝固傳熱的復(fù)合數(shù)值模擬[J].特殊鋼,2003,24(6):5-8.
[9]魏勇.板坯連鑄凝固傳熱及壓下位置預(yù)報(bào)模型的研究[D].武漢:武漢科技大學(xué),2007.