高爐球墨鑄鐵凸臺(tái)冷卻壁的溫度場(chǎng)的有限元仿真

劉 杰,黃曉然,杜 鋼

（1.本溪鋼鐵集團(tuán)技術(shù)中心,遼寧 本溪 117000;2.清華大學(xué),北京 100084;3.東北大學(xué) 材料與冶金學(xué)院,沈陽(yáng) 110004）

高爐球墨鑄鐵凸臺(tái)冷卻壁的溫度場(chǎng)的有限元仿真

劉 杰1,黃曉然2,杜 鋼3

（1.本溪鋼鐵集團(tuán)技術(shù)中心,遼寧 本溪 117000;2.清華大學(xué),北京 100084;3.東北大學(xué) 材料與冶金學(xué)院,沈陽(yáng) 110004）

為探討高爐凸臺(tái)冷卻壁的破壞過(guò)程,根據(jù)高爐的生產(chǎn)狀況,確定凸臺(tái)冷卻壁的熱邊界條件,采用ANSYS方法,計(jì)算凸臺(tái)冷卻壁的溫度分布;計(jì)算結(jié)果表明,凸臺(tái)冷卻壁的凸臺(tái)部位和失去鑲磚部位容易燒損;切斷冷卻循環(huán)水可以加速冷卻壁的燒損.

高爐;鑄鐵冷卻壁;侵蝕;溫度場(chǎng)

在高爐生產(chǎn)中,冷卻壁工作條件惡劣,它要不斷承受溫度的變化.溫度的變化將引起壁體熱應(yīng)力的產(chǎn)生,壁體熱應(yīng)力作為各種內(nèi)力之首應(yīng)引起重視.由于溫度變化產(chǎn)生的應(yīng)力超過(guò)其允許拉伸強(qiáng)度就會(huì)引起裂紋的產(chǎn)生,若超過(guò)疲勞強(qiáng)度,經(jīng)過(guò)多次不穩(wěn)定溫度變化則會(huì)引起疲勞裂紋,最終導(dǎo)致冷卻壁破損.冷卻壁大量損毀的部位是在爐襯侵蝕劇烈的爐腹,爐腰和爐身中下部,而其中破損最為嚴(yán)重的是第三代冷卻壁的凸臺(tái)部位[1].如何提高冷卻壁使用壽命以提高爐齡,確保高爐高效穩(wěn)定地生產(chǎn),已成為鋼鐵企業(yè)和相關(guān)行業(yè)的重要課題.

用數(shù)值仿真方法求解穩(wěn)態(tài)下高爐冷卻壁的溫度場(chǎng)和熱應(yīng)力已較普遍采用.研究成果對(duì)分析冷卻壁的破損原因和優(yōu)化冷卻壁設(shè)計(jì)參數(shù)等問(wèn)題都有一定指導(dǎo)意義.但研究熱點(diǎn)大多集中于銅或鑄鋼等新型材料冷卻壁,對(duì)目前我國(guó)高爐生產(chǎn)中普遍使用且破壞較為嚴(yán)重的三、四代鑲磚帶凸臺(tái)冷卻壁的研究卻比較少.

本文運(yùn)用有限元分析軟件 ANSYS的熱 -力耦合的計(jì)算方法,動(dòng)態(tài)模擬了我國(guó)高爐使用的第三代球墨鑄鐵鑲磚帶凸臺(tái)冷卻壁溫度場(chǎng)的變化狀況,為冷卻壁的破壞過(guò)程和維護(hù)提供依據(jù).

1 凸臺(tái)冷卻壁的計(jì)算模型

首先,以安裝在某高爐爐身下部的第三代球墨鑄鐵凸臺(tái)冷卻壁為實(shí)體,建立有限元計(jì)算模型.為模擬冷卻壁破壞過(guò)程,參照冷卻壁結(jié)構(gòu)參數(shù)及冷卻水布置方式[2,3],布置計(jì)算模型的水冷通路:縱向設(shè) 3條冷卻通道;本體背部增設(shè)一路蛇形管,在模型中簡(jiǎn)化為 4條橫向冷卻通路;凸臺(tái)處有一水平彎管.考慮到數(shù)值計(jì)算的規(guī)模,凸臺(tái)設(shè)在壁體中部,則根據(jù)模型對(duì)稱性,截取冷卻壁整體的右下部四分之一建模 （圖 1）.

圖 1 凸臺(tái)冷卻壁計(jì)算模型示意圖Fig.1 Schem a tic d iag ram o f the ca lcu la tingm ode l o f the coo ling stave

模型各部位尺寸設(shè)計(jì)為:L1—爐殼厚度,40 mm;L2—填充層厚度,40mm;L3+L4—冷卻壁本體厚度,320mm;L4—鑲磚寬度 （高爐半徑方向）,100mm;L5—凸臺(tái)深入爐內(nèi)的長(zhǎng)度,410 mm;L6—冷卻壁本體縱向冷卻水管距本體熱面距離,160 mm;L7—冷卻壁本體橫向冷卻水管距本體熱面的距離,220mm;L8—凸臺(tái)冷卻水管距凸臺(tái)熱面的距離;R—凸臺(tái)與冷卻壁本體銜接過(guò)渡半徑,110mm;冷卻壁本體縱向冷卻水管間距 120mm;本體橫向水管間距 400mm;凸臺(tái)U型冷卻水管半徑 70mm.

2 凸臺(tái)冷卻壁的三維傳熱數(shù)學(xué)模型

根據(jù)熱傳導(dǎo)定義,在沒(méi)有任何材料質(zhì)量純運(yùn)動(dòng)的情況下,熱通過(guò)材料的傳遞,沿 x方向上傳導(dǎo)的熱流速率 q為:

式中:k是材料的導(dǎo)熱系數(shù);A是垂直于 x熱流通過(guò)的面積;T是溫度;x是長(zhǎng)度參數(shù).

在沒(méi)有熱源的穩(wěn)態(tài)條件下,冷卻壁三維導(dǎo)熱微分方程為:

式中:T為溫度,℃;kx,ky,kz為與溫度相關(guān)的材料導(dǎo)熱系數(shù),W/（m·℃）;ρ為連續(xù)介質(zhì)的質(zhì)量密度,（g/cm3）;c為連續(xù)介質(zhì)的 C,（J/（kg·K））

式（2）是控制正交各向異性中熱傳導(dǎo)的微分方程.如果假設(shè) x,y和 z方向的熱傳導(dǎo)率相同,即常數(shù),則可將方程改寫為傅立葉方程:

式中,常數(shù) a=（k/ρc）為放熱系數(shù).

第三代球墨鑄鐵帶凸臺(tái)冷卻壁的傳熱可以看作導(dǎo)熱問(wèn)題來(lái)處理,在穩(wěn)態(tài)的工作條件下,三維導(dǎo)熱方程為:

式中:λ（T）是與溫度有關(guān)的導(dǎo)熱系數(shù) （W·m-1·℃-1）,i=1,2,3,表示 x,y和 z軸.

3 凸臺(tái)冷卻壁三維傳熱的邊界條件

冷卻壁溫度場(chǎng)計(jì)算涉及兩種邊界條件:一種是絕熱邊界條件;另一種是對(duì)流換熱邊界條件.冷卻壁與周圍環(huán)境之間、熱面與高溫煤氣之間、冷卻水與壁體之間屬于對(duì)流換熱邊界.具體的邊界條件方程如下:

（1）冷卻壁冷面、熱面邊界方程:

式中,h1,h2為冷卻壁冷面與周圍環(huán)境之間、冷卻壁熱面與高爐煤氣流之間的對(duì)流換熱系數(shù);TE,TB分別為高爐周圍空氣溫度、高爐爐溫.

（2）冷卻水與壁體之間的對(duì)流換熱邊界方程:

式中,hx,hy,hz是邊界外法線的方向余弦;h是冷卻水與壁體的傳熱系數(shù);Tw是冷卻水平均溫度.

（3）冷卻壁側(cè)面和側(cè)對(duì)稱面的絕熱邊界方程:

冷卻壁底面和上對(duì)稱面的絕熱邊界方程:

上述各式中,溫度 T是與時(shí)間無(wú)關(guān)、與坐標(biāo)有關(guān);導(dǎo)熱系數(shù)λ（T）與溫度有關(guān).

4 凸臺(tái)冷卻壁溫度場(chǎng)建模及數(shù)值計(jì)算

4.1 凸臺(tái)冷卻壁溫度場(chǎng)的有限元模擬

在有限元對(duì)實(shí)體模型進(jìn)行網(wǎng)格劃分過(guò)程中,為提高計(jì)算精度和滿足耐火材料建模的合理性,通常的劃分規(guī)則是在結(jié)構(gòu)平緩部位網(wǎng)格粗化,而在結(jié)構(gòu)復(fù)雜的部位網(wǎng)格適當(dāng)細(xì)化.考慮到構(gòu)成凸臺(tái)球墨鑄鐵冷卻壁材料的特殊性,本文采用混合建模方式,對(duì)形狀比較復(fù)雜的冷卻壁壁體,采用ANSYS SOL ID 98有限元四面體單元進(jìn)行網(wǎng)格劃分 （圖 2）;而對(duì)冷卻壁的鑲磚部分,則采用六面體單元 ANSYS SOL ID 5有限元單元 （圖 2）,總單元數(shù)為 53 314個(gè).

ANSYS SOL ID 5單元是六面體單元、八節(jié)點(diǎn)定義,每個(gè)節(jié)點(diǎn)有六個(gè)自由度,計(jì)算精度較高,具有三維溫度場(chǎng)、磁場(chǎng)、電場(chǎng)、壓電場(chǎng)與結(jié)構(gòu)場(chǎng)之間有限耦合和大變形的應(yīng)力剛化功能.ANSYS SOL ID 98十節(jié)點(diǎn)四面體單元是 ANSYS SOL ID 5單元的退化單元,其性質(zhì)同 ANSYSSOL ID 5單元.

將爐殼,壁體,填料,鑲磚視為連續(xù)體考慮,在ANSYS中設(shè)置各種材料隨溫度變化的物性參數(shù)（冷卻水管規(guī)格=60mm×5mm）.

圖 2 有限元網(wǎng)格劃分Fig.2 Schem a tic d iag ram o f the finitee lem en t d ivision

表 1 鑄鐵、鑲磚、爐殼及填料的物性參數(shù)Tab le 1 M a te ria lp rope rties o f cast iron,fireb rick,fu rnacem an tle and fillingm a te ria l

表 2 鑄鐵冷卻壁熱邊界條件Tab le 2 The rm a lbounda ry cond ition o f the coo ling stave

有限元計(jì)算中各材質(zhì)物性參數(shù)的選擇如表 1所示[3,4,5].鑄鐵冷卻壁的熱邊界條件列于表 2[4]:

4.2 凸臺(tái)冷卻壁溫度場(chǎng)的有限元模擬

在對(duì)凸臺(tái)冷卻壁的熱模型分析計(jì)算中,需做如下假設(shè):

（1）在模型的寬度 （z向）和高度 （y向）范圍內(nèi),冷卻壁熱面附近爐溫均勻一致;

（2）忽略冷卻壁與周圍耐火材料的傳熱,即將二者視作絕熱關(guān)系;

（3）用綜合換熱系數(shù)體現(xiàn)水管管壁厚度、水管與壁體間氣隙厚度等的影響;

（4）假設(shè)材料為各向同性,即各向?qū)嵯禂?shù)相等,遵守 VonM ises屈服準(zhǔn)則.

考慮到凸臺(tái)冷卻壁通常安裝在高爐的爐身下部、爐腰和爐腹部位,這些部位已經(jīng)有初渣形成,爐內(nèi)的溫度在 1 100℃以上,同時(shí)又考慮到鑄鐵冷卻壁的工作環(huán)境在 760℃以上時(shí),容易發(fā)生相變,致使自身強(qiáng)度下降,造成燒損.通過(guò)初步計(jì)算和文獻(xiàn)[4]的數(shù)據(jù),確定凸臺(tái)冷卻壁溫度場(chǎng)的熱邊界條件 （見表 2）.

凸臺(tái)冷卻壁溫度場(chǎng)的模擬計(jì)算分為兩部分:（1）凸臺(tái)冷卻壁破壞前平衡溫度場(chǎng);（2）凸臺(tái)冷卻壁破壞后平衡溫度場(chǎng).

4.2.1 凸臺(tái)冷卻壁破壞前的溫度場(chǎng)

凸臺(tái)冷卻壁在破損之前,外形是完整的,其溫度場(chǎng)的計(jì)算結(jié)果見圖 3和圖 4.

從圖 3可以看到,當(dāng)爐溫為 1 158℃,冷卻水流速 1.14m/s時(shí),凸臺(tái)冷卻壁的溫度范圍在 76℃到 988℃之間,最高溫度出現(xiàn)在凸臺(tái)冷卻壁的內(nèi)側(cè)面處,溫度達(dá)到了 988℃,在此高溫下,球墨鑄鐵的抗拉強(qiáng)度急劇下降,性能變壞,極易發(fā)生開裂,剝蝕等永久性破壞.

凸臺(tái)下端冷卻壁壁體與鑲磚接觸處的熱面溫度為 761℃（圖 4）,這是因?yàn)殍偞u部分導(dǎo)熱系數(shù)小,導(dǎo)熱能力差,造成鑲磚和鑲磚附近部位的鑄鐵溫度偏高.這個(gè)溫度也超過(guò)了球墨鑄鐵相變轉(zhuǎn)變溫度 760℃,達(dá)到了破損的溫度極限.凸臺(tái)上冷卻水管流經(jīng)處的溫度較低,不會(huì)被燒損.計(jì)算結(jié)果說(shuō)明了第三代凸臺(tái)冷卻壁的凸臺(tái)部位和與鑲磚接觸部位是該冷卻壁的薄弱環(huán)節(jié),極易被燒蝕破損.

4.2.2 凸臺(tái)冷卻壁破壞后的溫度場(chǎng)

隨著凸臺(tái)冷卻壁凸臺(tái)部位的破損和鑲磚的融蝕脫落,冷卻壁的形狀已殘缺不全.通過(guò)對(duì)凸臺(tái)冷卻壁施加熱負(fù)荷和力學(xué)邊界條件,對(duì)冷卻壁的破損形狀進(jìn)行了預(yù)測(cè),在殘余冷卻壁上重新進(jìn)行網(wǎng)格劃分、計(jì)算,得到殘余凸臺(tái)冷卻壁的溫度場(chǎng) （圖5）.由于凸臺(tái)冷卻壁破損狀況不同,溫度場(chǎng)的分布形態(tài)也略有區(qū)別.從圖 5中可以看出,最高溫度依然出現(xiàn)在凸臺(tái)前端,達(dá)到 932℃,該溫度遠(yuǎn)高于球墨鑄鐵相變溫度 760℃,極易造成燒損.

圖 5 殘缺凸臺(tái)冷卻壁的溫度分布圖Fig.5 The tem pe ra tu re fie ld o f thein com p le te coo ling stave

凸臺(tái)冷卻壁的凸臺(tái)部分處在高爐高溫高負(fù)荷區(qū),往往先被燒壞,但凸臺(tái)部分的燒蝕速度,不僅受爐內(nèi)高溫和高壓的影響,還取決于凸臺(tái)部位冷卻水管完好與否.冷卻水管的冷卻作用,可以減緩?fù)古_(tái)的燒損進(jìn)度.當(dāng)凸臺(tái)處冷卻水管部分或完全裸露在高溫煤氣流中時(shí),一旦斷水,冷卻水管很快損毀,導(dǎo)致凸臺(tái)處失去冷卻作用.加快凸臺(tái)部位的燒蝕.

冷卻壁的鑲磚部位的燒損也有類似的現(xiàn)象.冷卻壁的鑲磚部分一旦融蝕,冷卻壁熱面處的壁體失去耐火材料的保護(hù)而暴露在高溫煤氣流中,局部溫度不斷升高,加速冷卻壁的破損.一旦斷水,后果不可想象.因此,當(dāng)冷卻壁漏水不得不關(guān)閉循環(huán)水的情況下,應(yīng)采取在冷卻壁上開孔,加裝簡(jiǎn)易圓筒冷卻器,延緩整塊冷卻壁的燒損速度,避免溫度過(guò)高,爐墻燒漏、爐體燒穿等重大事故的發(fā)生.

5 結(jié) 論

（1）以第三代球墨鑄鐵凸臺(tái)冷卻壁為實(shí)體模型,建立有限元溫度場(chǎng)的數(shù)值模型;

（2）以球墨鑄鐵冷卻壁抗拉強(qiáng)度曲線作為判定依據(jù),確定凸臺(tái)冷卻壁溫度場(chǎng)的熱邊界條件,采用ANSYS計(jì)算凸臺(tái)冷卻壁的溫度分布.

（3）計(jì)算結(jié)果表明,凸臺(tái)冷卻壁的凸臺(tái)部位和鑲磚的脫落部位容易造成冷卻壁的燒損;而凸臺(tái)冷卻水管和冷卻壁內(nèi)冷卻水管的失效,最終導(dǎo)致鑄鐵凸臺(tái)冷卻壁的完全破壞.

[1]曹傳根,周渝生,曹進(jìn),等.寶鋼 3號(hào)高爐冷卻壁破損機(jī)理的熱態(tài)實(shí)驗(yàn)研究.鋼鐵,1999,34（4）:5.

（CAO Chuan-gen,ZHOU Yu-sheng,CAO Jin,etal.The ho t test study on the stuve of baosteels no3 b last fu rnace[J].Iron and Steel,1999,34（4）:5.）

[2]周傳典.高爐煉鐵生產(chǎn)技術(shù)手冊(cè).北京:冶金工業(yè)出版社,2004.

（ZHOU Chuan-d ian.H andbook of b last fu rnace ironm ak ing p roduction techno logy[M]. Beijing:M etallu rgical Industry Press,2004.）

[3]吳俐俊,程惠爾,錢中.冷卻水管管形變化下的高爐冷卻壁傳熱分析[J].鋼鐵,2005,40（5）:15.

（W U L i-jun;CHENG H ui-er;Q IAN Zhong.Therm al stressand therm ald isto rtion of cast iron coo ling stavew ith su rface-alloyed coo lingw aterp ipes[J].Iron and Steel,2005,40（5）:15.）

[4]石琳,程素森,馮力,等.冷卻水管表面合金化球墨鑄鐵冷卻壁的熱應(yīng)力和熱變形 [J].北京科技大學(xué)學(xué)報(bào),2007,29（9）:942.

（SH IL in,CH ENG Su-sen,FENG L i,et al.Therm al stress and therm al d isto rtion of cast iron coo ling stave w ith su rfacealloyed coo lingw ater p ipes[J].Jou rnalofU niversity of Science and Techno logy Beijing,2007,29（9）:942.）

[5]H ebelR,Steiger R,Streuber C.A dvanced m aterials fo r the b last fu rnace lining in regard to long serv ice life,reasonab le cost and their theo retical qualification[C]//Ironm ak ing Conf Proc.1994,53:397.

F in ite elem en tana lysis of ductile cast iron coo ling stave fa ilure process in b last furnace

L IU Jie1,HUAN G X iao-ran2,D u Gang3
（1.Techno logy Cen ter of Benx i Iron and SteelCo rp.,Benx i117000,China;2.TsinghuaU niversity.,Beijing,100084,Ch ina;3.Schoo lofM aterials and M etallu rgy,N o rtheastern U niversity,Shenyang 110004,China）

In o rder to research the failu re p rocess the ductile iron coo ling stave,calcu late the tem peratu re field and attain the therm al boundary cond ition,u tilize the AN SYS therm alm echan ical fin ite elem en t sim u lation acco rd ing to b last fu rnace p roduction status. The resu lts of the calcu lation show ed that the coo ling stave is easily bu rned up on convex p late and the bu lge w h ich lose the shedding firebricks;. C u tting off the coo ling w ater w ill accelerate the bu rning speed of the coo ling stave.

b last fu rnace;nodu lar cast iron coo ling stave;erosion; tem peratu re field

TB 332

A

1671-6620（2010）01-0007-05

2009-05-15.

劉杰 （1951—）,男,遼寧營(yíng)口人,本溪鋼鐵集團(tuán)工程師,Email:lucyking@vip.sana.com;杜鋼 （1954—）男,遼寧鞍山人,東北大學(xué)教授.