冒口對鋼錠凝固過程固相轉變影響的數值模擬

袁林華，于景坤，李志強

（1沈陽東大冶金科技股份有限公司，遼寧沈陽 110004；2東北大學材料與冶金學院，遼寧沈陽 110004）

試驗研究

冒口對鋼錠凝固過程固相轉變影響的數值模擬

袁林華1，于景坤2，李志強2

（1沈陽東大冶金科技股份有限公司，遼寧沈陽 110004；2東北大學材料與冶金學院，遼寧沈陽 110004）

通過自建模型模擬研究了冒口材質和高度對大鋼錠凝固過程中固相轉變的影響。結果表明：鋼錠凝固過程中，錠模及底座附近的鋼液發生相轉變的速率較快且固液共存區寬度較窄，而冒口附近的鋼液發生固相轉變的速率較慢且糊狀區較寬。相比于傳統的黏土磚，輕質磚和絕熱板均有助于改善冒口的保溫性能，且均可減緩冒口外表面溫度上升的速率。冒口高于鋼液上表面部分高度的變化對于鋼錠凝固過程中冒口附近鋼液相轉變速率的影響較小。

冒口；鋼錠凝固；數值模擬

隨著電力工業的發展和發電機組的大型化，低壓轉子鍛件的尺寸也隨之增加。尺寸增加給低壓轉子大鍛件的煉鋼、鑄錠、鍛造和熱處理都帶來了非常大的挑戰［1-3］。由于實驗研究代價高，數值模擬成為認識這些缺陷和優化鋼錠設計的重要途徑［4］。

通過數值模擬方法可以得到鋼錠凝固過程中鋼錠本身及錠模的溫度場變化狀況［5-6］。錠模冒口的使用可以有效改善鋼錠凝固過程中錠頭凝固的固相率轉變狀況［7］。大型鋼錠凝固過程的三維數值模擬結果表明冒口頂部附近鋼液固相轉變較快［8］，但鮮有關于冒口內襯材料及高度對錠頭固相率轉變影響的報道。

本研究通過自建模型模擬探討某48 t鋼錠凝固過程中固相轉變的狀況及冒口內襯耐火材料種類和高度對鋼錠凝固過程中固相轉變的影響。

1 研究對象

以國內某重型設備制造公司的48 t轉子鍛造鋼錠為研究對象（化學成分見表1，其中J系數為10），錠模內壁為多邊形。為簡化計算過程，本文選取此鋼錠及錠模縱截面的一半為研究對象（如圖1所示）。此外，為考察冒口內襯耐火材料對鋼錠凝固的影響，選取輕質磚和絕熱板與傳統的黏土磚進行對比研究。各材質的物性參數見表2。

2 模型建立

2.1 假設條件

假設鑄型瞬間充滿，液態金屬的初始溫度即為澆注溫度；液態金屬內部無對流作用，無能量與質量傳輸；鋼錠及錠模的材料均質且各向同性。

表1 鋼種化學成分（質量分數）%

圖1 模型計算區域示意圖

表2 材料的熱物性參數

2.2 傳熱控制方程及潛熱處理

在鋼錠（包括固相和液相）和錠模內部，溫度變化主要通過熱傳導的方式進行。在材料均質且各向同性的條件下，有內熱源的非穩態傳熱過程可通過傅立導熱方程來描述：

式中：ρ為密度，kg/m3；T為溫度，K；cp為定壓比熱容，J/（kg·K）；t為時間，s；Q為熱源項，J；x、y、z均為方向分量。

忽略鋼液凝固過程中液態金屬的流動，鋼錠凝固過程溫度變化及分布即可由式（1）求得。式（1）中Q為鋼液凝固發生相變時產生的熱量，即潛熱。

本模型采用等價比熱法和溫度回升法相結合的方法處理鋼液凝固過程中潛熱的釋放。在一個時間步內，當計算區域的溫度變化范圍包含液相線（或固相線）溫度時，潛熱處理采用溫度回升法；當計算區域溫度變化在固液相線溫度之間時，潛熱處理采用等價比熱法。

2.3 體積收縮處理

在t到t+dt時間步長內，所有固液兩相區單元格的體積收縮總量為：

式中：n為兩相區單元格總數，Vi為單元i的體積。將△V平均分配到鋼液上表面單元格內。當這些單元格的體積收縮達到單元格體積時，將此單元格的屬性向下移動到鄰近的單元格，并將該單元格狀態更新為空白，不再參與下一個步進時間的計算。

2.4 初始條件和邊界條件

鋼液凝固初始條件主要包括鋼液溫度和錠模溫度，其中，設鋼液初始溫度為1 540℃，錠模溫度為室溫（30℃），冒口預熱溫度為200℃。

鋼液凝固熱分析邊界條件主要包括4類：1）中軸線和鋼液上表面設為絕熱邊界條件；2）錠模底部和側部設為綜合考慮輻射和對流的等效對流換熱邊界條件；3）冒口內襯與外側鋼板設為緊密接觸傳熱邊界條件；4）鋼錠-錠模邊界設為輻射換熱等效對流傳熱邊界條件。

2.5 模型計算方法

采用顯式有限差分法對模型進行離散處理并求解，步進時間為ds2/（4DL）min。

3 模擬結果與討論

3.1 鋼錠凝固過程鋼液狀態及溫度場變化

鋼液凝固過程中，伴隨著溫度的降低，鋼液向固態轉變，模擬結果如圖2所示。

圖2 鋼錠凝固過程狀態變化

從圖2a～c可以看出，鋼錠凝固初期，錠模及底座附近的鋼液發生相轉變的速率較快且糊狀區寬度較窄，而冒口附近的鋼液發生相轉變的速率較慢且糊狀區較寬，在冒口與鋼液的交界面上，鋼液上表面交界處的鋼液發生相轉變的時間比其它位置更早。而圖2d～f表明，液相消失后，固液兩相區內熔融鋼液向固相鋼錠轉變是從其上表面與冒口交界處和軸向向上兩個方向進行的。鋼液上面與冒口交界處鋼液的快速冷卻凝固有可能導致最終產品產生疏松缺陷，應盡量避免。

3.2 冒口內襯材料對鋼錠凝固的影響

3.1 節結果和分析表明冒口內襯材料對鋼錠凝固過程的影響較大。與傳統的黏土磚相比，若改為輕質磚或絕熱板，鋼錠凝固過程冒口附近錠頭凝固的固相轉變對比結果及冒口外壁的溫度變化分別如圖3、圖4所示。

圖3 冒口內襯材料對鋼錠凝固過程的影響

圖4 內襯材料對錠模外表面溫度的影響

由圖3可知，在相同時刻，相比于傳統的黏土磚，選用輕質磚有利于改善錠頭上表面與冒口內襯交界處固相轉變的狀態變化；若選用絕熱板，則錠頭上表面與冒口內襯交界處固相轉變的狀況會和冒口內襯邊界其他位置的狀態非常接近，更有利于改善鋼液凝固的過程，可以為鋼錠的最終質量提供保證。

由圖4可知，輕質磚和絕熱板均可減緩冒口外壁溫度上升的速率。例如，當凝固進行了4 h時，黏土磚內襯冒口的外表面溫度為900 K，輕質磚內襯冒口的外表面溫度為750 K，而絕熱板內襯冒口的外表面溫度僅為500 K。保溫性能的改善有利于延長錠頭的凝固時間，從而給夾雜物的上浮提供充足的時間。

3.3 冒口高度對鋼錠凝固的影響

由于冒口高出鋼液（錠）上表面部分溫度升高需要的熱量只能來源于鋼液（錠）上表面與冒口內表面的交界處的鋼液（錠），這部分的固相轉變速率較快。當選用黏土磚為冒口內襯時，冒口高度對于鋼錠凝固的影響結果如圖5所示。

由圖5可知，在同一時刻，冒口高于鋼液（錠）上表面部分高度的變化對于鋼錠凝固過程中冒口部分固相轉變速率的影響較小。

圖5 冒口高度對鋼錠凝固的影響

4 結論

4.1 鋼錠凝固過程中，錠模及底座附近的鋼液發生相轉變的速率較快且糊狀區寬度較窄，而冒口附近的鋼液發生相轉變的速率較慢且糊狀區較寬。

4.2 相比于傳統的黏土磚，輕質磚和絕熱板均有助于改善冒口的保溫性能，且均可減緩冒口外表面溫度上升的速率。

4.3 冒口高于鋼液上表面部分的高度變化對于鋼錠凝固過程中冒口附近鋼液相轉變速率的影響比較小。

［1］劉顯惠，林錦棠.國內外汽輪機大型轉子鍛件材料的技術進展（一）［J］.國外金屬熱處理，1999（3）：5-7.

［2］雷雪，徐薇平，韓利戰，等.低壓轉子加熱過程溫度場與相變的計算機模擬［J］.金屬熱處理，2009，34（9）：48-51.

［3］劉鑫，鐘約先，馬慶賢，等.核電汽輪機低壓轉子技術的發展［J］鍛壓裝備與制造技術，2009（3）：13-18.

［4］Radovic Z，Lalovic M.Numerical simulation of steel ingot solidification process［J］.Journal of Materials processing Technology，2005（160）：156-159.

［5］王芹，袁守謙，鄧林濤，等.鍛造用鋼錠凝固過程溫度場數值模擬及其應用［J］.大型鑄鍛件，2005（1）：10-12.

［6］董潔，袁守謙，鄧林濤，等.鍛造用鋼錠凝固過程溫度場數值模擬［J］.鑄造技術，2007，28（2）：268-270.

［7］馬薇，趙建華，宋剛.冒口保溫條件對大型鋼錠凝固過程影響的數值模擬［J］.熱加工工藝，2011，40（21）：41-43.

［8］李文勝，沈丙振，周翔，等.大型鋼錠凝固過程三維數值模擬［J］.大型鑄鍛件，2010（3）：1-4.

Numerical Simulation on theEffect of Riser on Solid PhaseTransformation in Solidification Processof Ingot

YUAN Linhua1,YU Jingkun2,LI Zhiqiang2

（1 Northeastern University Metallurgical Technology Institute Co.,Ltd.,Shenyang 110004,China;
2 School of Material and Metallurgy,Northeastern University,Shenyang 110004,China）

The effects of the materials and height of riser on the solid phase transformation of steel in the solidification process of an steel ingot was simulated through a self designed model.The results showed that during the solidification process of the ingot,the transforming rate of steel from liquid to solid was fast and the width of the liquid/solid coexisting area was narrow nearby the side and bottom of mould,while the transforming rate of steel was slow and the width of the area was wide nearby the riser.Compared with the traditional clay brick,the applications of the light weight brick and the heat insulating shield not only would improve the heat insulating property of the riser but also slow down the increasing rate of temperature outside the riser.The effect of the height change of riser which was higher than the top surface of ingot on the transforming rate of steel nearby the riser in the solidification process was tiny.

riser;ingot solidification;numerical simulation

TP391；TF771+.1

A

1004-4620（2014）03-0026-03

2014-04-08

袁林華，男，1962年生，1983年畢業于東北大學冶金工程專業。現為沈陽東大冶金科技股份有限公司工程師，從事冶金輔料研究工作。