冷凝型礦熱爐多物理場數(shù)值模擬研究

摘 """""要： 為了研究冷凝型礦熱爐內(nèi)應(yīng)力分布情況，基于多相流、多物理場以及冷凝型礦熱爐爐襯結(jié)構(gòu)的特點，構(gòu)建冷凝型礦熱爐多物理場耦合的傳熱計算模型。利用CFD仿真手段合理地預(yù)測冷凝型礦熱爐爐襯的溫度場分布情況，進行溫度場的分析?；跓?結(jié)構(gòu)耦合理論和冷凝型礦熱爐爐襯結(jié)構(gòu)特點，構(gòu)建冷凝型礦熱爐工作狀態(tài)下爐襯熱應(yīng)力仿真計算模型，分析出冷卻水溫度與爐襯應(yīng)力和變形間的關(guān)系。

關(guān) "鍵 "詞：礦熱爐；冷凝型礦熱爐；數(shù)值模擬

中圖分類號：TF806.7""""""文獻標(biāo)志碼：A """"文章編號：1004-0935（2024）10-1574-04

目前，回轉(zhuǎn)窯-礦熱爐（RKEF）工藝因所冶煉的鎳鐵合金具有品位高、有害元素少、生產(chǎn)效率高和工藝成熟等優(yōu)點^[1]已被廣泛應(yīng)用于鎳鐵合金冶煉工藝中。礦熱爐是生產(chǎn)多種鐵合金的重要設(shè)備。根據(jù)中國鐵合金行業(yè)發(fā)展研究報告，中國具有世界上產(chǎn)量最大、品種最全的鐵合金生產(chǎn)體系，近年來鐵合金總體產(chǎn)量一直維持在3 000萬t以上。由于鐵合金生產(chǎn)屬于高耗能高排放行業(yè)，為了實現(xiàn)鐵合金生產(chǎn)的節(jié)能與高質(zhì)量發(fā)展，政府通過多項政策推動礦熱爐的大型化發(fā)展。大型礦熱爐具有熱效率更高、產(chǎn)量更高、操作更穩(wěn)定等優(yōu)勢，也利于煙塵的凈化回收和余熱利用。但在礦熱爐大型化后，爐襯的砌筑與維修成本明顯增加，所以對礦熱爐爐襯進行研究勢在必行。

由于礦熱爐內(nèi)部條件復(fù)雜，難以通過實驗的方法進行研究，多數(shù)學(xué)者都是采用數(shù)值模擬的方法展開研究。劉奇等^[2]用數(shù)值模擬方法對銅鋼復(fù)合冷卻壁和銅冷卻壁進行了比較，發(fā)現(xiàn)銅鋼復(fù)合銅冷卻壁抗變形能力優(yōu)于銅冷卻壁。孟凡東等^[3]利用ANSYS熱-結(jié)構(gòu)耦合仿真計算功能進行爐墻熱應(yīng)力計算和分析。計算分析不同內(nèi)襯剩余厚度和爐殼開裂長度、寬度大小對內(nèi)襯外緣環(huán)向應(yīng)力的影響。GAN等^[4]通過建立熱風(fēng)爐的三維熱力學(xué)模型計算了運行過程中的溫度和應(yīng)力分布，分析出最大主應(yīng)力主要發(fā)生在耐火襯里外，而最小主應(yīng)力主要發(fā)生在耐火襯里內(nèi)。ANIL等^[5]通過數(shù)值模擬分析了碳化硅磚和高鋁磚在傳熱過程中的性能，分析出性能較好的襯里是碳化硅磚。WEI等^[6]利用ANSYS模擬了氣化爐運行中的應(yīng)力場，研究了熱應(yīng)力對耐火材料和爐渣層的破壞。CHEN等^[7]采用等效對流傳熱系數(shù)簡化冷卻爐結(jié)構(gòu)，建立了侵蝕狀態(tài)下高爐爐襯熱應(yīng)力軸對稱計算模型，分析了常見侵蝕狀態(tài)下的熱應(yīng)力。ZHONG等^[8]利用有限元軟件ANSYS建立了銅冷卻槽的熱結(jié)構(gòu)耦合模型，計算出銅冷卻爐能滿足高爐最壞條件下的要求，在高溫負荷下不會產(chǎn)生疲勞裂紋。孫昊等^[9]研究了不同電極插入深度對于爐內(nèi)溫度場分布的影響，并對熔池內(nèi)的熔化凝固情況進行了分析。

以鎳鐵冶煉工藝中48 MW冷凝型礦熱爐為研究對象，利用計算流體動力學(xué)（CFD）對冷凝型礦熱爐溫度場與應(yīng)力場進行模擬，預(yù)測了爐襯的溫度場與應(yīng)力場的分布規(guī)律，為礦熱爐的安全高效生產(chǎn)提供理論依據(jù)。

1 "模型建立與網(wǎng)格劃分

1.1 "模型建立

冷凝型礦熱爐結(jié)構(gòu)如圖1所示。圖1（a）為冷凝型礦熱爐物理模型結(jié)構(gòu)。爐殼、爐蓋以及懸臂梁均采用Q235鋼，爐襯結(jié)構(gòu)由高鋁磚、鎂磚與耐火顆粒組成，爐底由鎂磚、鎂制搗打料和高鋁磚構(gòu)成。外層建有50"mm的冷卻水層用來模擬實際工況下的噴淋式冷卻。為了提高數(shù)值模擬的精確度，采用ICEM軟件劃分精度更高的六面體結(jié)構(gòu)化網(wǎng)格，網(wǎng)格總量為163萬，相較于四面體網(wǎng)格，六面體網(wǎng)格具有較高的精確度，并且在達到計算精度的要求下，所需要的網(wǎng)格數(shù)要比四面體網(wǎng)格少，從而達到節(jié)省計算時間和節(jié)約計算機資源的目的。礦熱爐爐體六面體網(wǎng)格如圖1（b）所示。

2 "數(shù)學(xué)模型

2.1 "熱應(yīng)力數(shù)學(xué)模型

根據(jù)熱彈性理論，爐膛的熱變形行為可以用以下方程^[10]來描述。方程（1）、方程（2）、方程（3）分別為平衡方程、運動學(xué)方程和本構(gòu)方程。

（1）

（2）

（3）

式中：T₀—初始溫度；

σ_i—平行于x、y、z軸的熱應(yīng)力分量，i=x，y，z；

ε_i—平行于x、y、z軸的熱應(yīng)變分量，i=x，y，z；

τ_xy、τ_yz、τ_zx和γ_xy、γ_yz、γ_zx—剪應(yīng)力的分量；

u、v、w—位移分量；

E—拉伸彈性模量和壓縮彈性模量；

μ—泊松比；

α—熱膨脹系數(shù)。

2.2 "氧化鎳反應(yīng)方程

研究的是礦熱爐在生產(chǎn)鎳鐵時的冶煉狀態(tài)，為更好地還原在冶煉過程中各物質(zhì)的參數(shù)，在計算時加入了化學(xué)反應(yīng)機理。氧化鎳在碳熱還原時發(fā)生的主要化學(xué)反應(yīng)由方程（4）至方程（6）來表示。

3Fe₂O₃（s）+C（s）=2Fe₃O₄（s）+CO（g） """"（4）

Fe₃O₄（s）+C（s）=3FeO（s）+CO（g） """""（5）

NiO（s）+C（s）=Ni（s）+CO（g） """""""（6）

3 "結(jié)果與討論

3.1 "溫度場分析

礦熱爐爐體溫度分布如圖2所示。由圖2可以看出，爐膛底部高溫區(qū)出現(xiàn)在坩堝區(qū)正下方，最高溫度為2 334 K。由于耐火材料高熱阻特性，爐體溫度梯度明顯高于熔池區(qū)域。爐襯內(nèi)壁與外爐殼間的溫度則由1 744 K降至308 K，較高的溫度梯度是爐體產(chǎn)生熱應(yīng)力熱變形的主要原因。

3.2""不同冷卻水溫度時的變形量分析

不同冷卻水溫度下爐殼變形量對比曲線如圖3所示。由圖3可以看出，隨著冷卻水溫度的升高，爐殼的變形量在逐漸增大。最大變形位置出現(xiàn)在z=2.57～2.84 m處，這是受爐內(nèi)高溫渣-金層的影響，此處爐襯溫度最大所。當(dāng)冷卻水溫度T=50 ℃時爐殼變形量最大，為13.45 mm；冷卻水溫度T=25 ℃時最大變形量降為13.12 mm。

3.2 "不同冷卻水溫度時的應(yīng)力分析

不同冷卻水溫度下爐殼Von Mises等效應(yīng)力對比曲線如圖4所示。

由圖4可以看出，隨著冷卻水溫度的升高，外爐殼的應(yīng)力在逐漸降低。爐殼底部則因為外爐殼與爐底的連接工藝是焊接，約束爐體的熱膨脹而產(chǎn)生了較大的擠壓應(yīng)力，這一應(yīng)力達到了160 MPa。在外爐殼z=2.73～3.03 m處，爐殼應(yīng)力最大，在水溫T=25 ℃時，爐殼最大應(yīng)力為203 MPa；T=50 ℃時，爐殼最大應(yīng)力降為181 MPa。這是因為：結(jié)構(gòu)體部件的熱膨脹或擠壓收縮受到外界約束；結(jié)構(gòu)體內(nèi)某一部件中的溫度梯度過大；不同線膨脹系數(shù)和爐襯材料組合的影響。在外爐殼z=6.2～6.8 m處出現(xiàn)應(yīng)力集中現(xiàn)象，這是由于在爐殼上部區(qū)域存在支撐平臺承載著爐蓋懸臂梁，約束外爐殼形變。可以發(fā)現(xiàn)，冷卻水進口溫度的變化對外爐殼應(yīng)力的影響較為明顯，提高水溫可以降低外爐殼應(yīng)力，但有使外爐殼形變過大的風(fēng)險。

4 "結(jié) 論

1）水冷型礦熱爐外爐殼熱應(yīng)力主要集中于爐殼底部和z=2.73～3.03 m高度。外爐殼底部應(yīng)力集中是由于爐殼與爐底是處于焊接狀態(tài)，阻礙了耐火材料的進一步膨脹，故應(yīng)力隨著高度的減小而變大；而z=2.73～3.03 m處應(yīng)力集中是因為：爐襯結(jié)構(gòu)的熱膨脹或擠壓收縮受到外界約束；爐襯結(jié)構(gòu)中溫度梯度過大；不同線膨脹系數(shù)爐襯材料組合的影響。

2）隨著冷卻水溫度的升高，外爐殼熱應(yīng)力隨之降低，爐殼變形量逐漸增大，這是由于冷卻水溫度越高，爐殼的膨脹量越大，爐襯內(nèi)部溫度梯度降低，所以使得爐殼應(yīng)力逐漸減小。

參考文獻：

[1] LIU P， LI B K， CHEUNG S C P， et al. Material and energy flows in rotary kiln-electric furnace smelting of ferronickel alloy with energy saving[J]. Applied Thermal Engineering， 2016， 109： 542-559.

[2]"劉奇，程樹森，趙宏博，等.銅鋼復(fù)合冷卻壁熱變形分析[J]. 工程科學(xué)學(xué)報，2016，38（1）：108-117.

[3] 孟凡東. 高爐爐缸爐殼開裂的三維應(yīng)力分析[D].沈陽：東北大學(xué)，2014.

[4] GAN Y F， JANG J Y， WU T Y. 3D dynamic thermal and thermo-"mechanical stress analysis of a hot blast stove[J]. Ironmaking amp; Steelmaking， 2020， 47（9）： 959-972.

[5] ANIL K， NANDAN B S， RITURAJ C. Computational modeling of blast furnace cooling stave based on heat transfer analysis[J]. Materials Physics and Mechanics， 2012， 15（1）： 46-65.

[6] WEI Y，YANG W G，QIN H B， et al. Numerical simulation of thermal stress field of water-cooled-wall pulverized-coal gasifiers[J]. China's Refractories， 2019，28（4）： 17-21.

[7] CHEN L Y， LI Y， GUI J H. Thermal stress of blast furnace hearth linings under erosion state[J]. Applied Mechanics and Materials， 2009， 16-19： 1101-1105.

[8] ZHONG Y， LU X Z， FAN X M， et al. Thermal simulation study on copper cooling stave of blast furnace[J]. Advanced Materials Research， 2015， 1088： 739-743.

[9]"孫昊，劉鵬，張新杰，等.電極插入深度對礦熱爐加熱影響的CFD分析[J].遼寧化工，2022，51（7）：991-993.

[10] WANG L， CHEN L Y， ZHAO L， et al. Thermofluid-solid coupling numerical simulation model of blast furnace hearth protection measures[J]. IEEE Access， 2022， 10： 65080-65091.

Numerical Simulation Study on Multiple Physical

Fields of Condensing Submerged Arc Furnace

LI Yaqi， SUN Hao， LIU Peng， KOU"Liping， DU Jinqi， LAN Xiaoyu

（School of Mechanical and Power Engineering， Shenyang University of Chemical Technology， Shenyang Liaoning 113142， China）

Abstract：""In order to study the stress distribution in the condensing submerged arc furnace， based on the characteristics of multiphase flow， multiple physical fields and the lining structure of the condensing submerged arc furnace， a coupled heat transfer calculation model of multiple physical fields of the condensing submerged arc furnace was constructed. CFD simulation was"used"to reasonably predict the temperature field distribution of the condensing submerged arc furnace lining and analyze the temperature field. Based on the thermal structure coupling theory and the structural characteristics of the lining of a condensing submerged arc furnace， a simulation calculation model for the thermal stress of the lining under working conditions of the condensing submerged arc furnace was"constructed， and the relationship between the cooling water temperature and the stress and deformation of the lining was"analyzed.

Key words：""Submerged arc furnace; Condensing submerged arc furnace; Numerical simulation