牛立群　袁　旗　王旭明　楊德生　馮廷彪

(1.蘭州蘭石能源裝備工程研究院有限公司，甘肅730314；2.甘肅省高端鑄鍛件工程技術(shù)研究中心，甘肅730314；3.蘭州蘭石鑄鍛有限責(zé)任公司，甘肅730314)

?

錠模斜度及澆注工藝對(duì)鋼錠質(zhì)量影響的數(shù)值模擬研究

牛立群1,2袁旗1,2王旭明1,2楊德生1,2馮廷彪3

(1.蘭州蘭石能源裝備工程研究院有限公司，甘肅730314；2.甘肅省高端鑄鍛件工程技術(shù)研究中心，甘肅730314；3.蘭州蘭石鑄鍛有限責(zé)任公司，甘肅730314)

采用THERCAST模擬軟件，模擬35 t鋼錠的澆注和凝固過(guò)程，分析鋼錠內(nèi)部缺陷形成的原因。模擬計(jì)算錠模斜度、澆注溫度和澆注速度對(duì)缺陷分布的影響，得出材料為30CrNiMo8的35 t鋼錠最佳澆注工藝方案。

THERCAST；錠模斜度；澆注溫度；澆注速度

1　有限元模型建立

鋼錠屬于軸對(duì)稱幾何體，為提高效率，節(jié)省計(jì)算時(shí)間，通過(guò)軟件的對(duì)稱功能，選取需要模擬的十二棱鋼錠體的1/12作為計(jì)算基體。有限元模型分為七個(gè)部分：鋼錠，錠模，底盤，底磚，冒口套，發(fā)熱劑，碳化稻殼。鋼錠單元數(shù)是45 844個(gè)。

2　初始條件和邊界條件

2.1加載材料

模擬材料的成分見表1。

表1　化學(xué)成分(質(zhì)量分?jǐn)?shù)，%)

2.2計(jì)算模式選擇

本次仿真分析選擇熱力學(xué)計(jì)算。熱交換及摩擦設(shè)置：金屬液與環(huán)境之間的熱交換選擇adiabatic(絕熱)，模具與模具之間的熱交換及摩擦選擇mold_mold_resistance和Constant_Coulomb，金屬與模具之間的熱交換及摩擦選擇air_gap_dependent_resistance和Sliding。

2.3澆注屬性設(shè)置

澆注屬性設(shè)置的主要參數(shù)有澆注溫度和澆注速度。澆注過(guò)程和冷卻過(guò)程中熱學(xué)條件設(shè)置為：環(huán)境溫度20℃，散熱系數(shù)12。澆注完成后的法向速度設(shè)定為0 mm/s。

3　縮松形成分析

鋼水由液態(tài)變成固態(tài)的過(guò)程中，需要經(jīng)歷三個(gè)階段，即完全液態(tài)、液態(tài)與固態(tài)共存、完全固態(tài)[3]。鋼錠在這三個(gè)階段中，都會(huì)發(fā)生體積收縮的現(xiàn)象。如果鋼液最后凝固區(qū)域無(wú)法得到外部鋼液的補(bǔ)充，則會(huì)產(chǎn)生孔洞。Campbell提出凝固過(guò)程中的五種補(bǔ)縮方式，包括固體補(bǔ)縮方式、爆炸充填方式、質(zhì)量補(bǔ)縮方式、枝晶補(bǔ)縮方式以及液體補(bǔ)縮方式。這些補(bǔ)縮方式以液體補(bǔ)縮作用尤為重要，所以需要模擬液相分?jǐn)?shù)的變化規(guī)律。

鋼錠凝固過(guò)程中液相分布如圖1所示?？梢钥闯觯挥袖撳V中心軸線部位凝固最慢，距離冒口約1/3處，液相分?jǐn)?shù)大約0.9的區(qū)域，被完全孤立出來(lái)兩個(gè)區(qū)域，這就可能導(dǎo)致鋼錠在此處產(chǎn)生兩處縮松缺陷。當(dāng)凝固至885步時(shí)，可以清楚觀察到，此時(shí)在距離冒口約2/3處，液相分?jǐn)?shù)大于0.9的黑色區(qū)域，也分離出來(lái)2個(gè)區(qū)域，表明在這個(gè)位置產(chǎn)生縮松缺陷可能性較大。當(dāng)凝固至923步時(shí)，可以清楚的看到，液相分?jǐn)?shù)大于0.9的黑色區(qū)域呈現(xiàn)倒放漏斗狀，隨著凝固過(guò)程推進(jìn)，此時(shí)冒口就不能達(dá)到對(duì)鋼錠的補(bǔ)縮作用。

圖2給出了模擬凝固時(shí)Niyama分布圖。可以看出，鋼錠的縮孔出現(xiàn)在軸線處，尤其是冒口中心位置和距離冒口頂端軸線1/3處相對(duì)嚴(yán)重。圖3為解剖鋼錠疏松缺陷分布圖，模擬疏松缺陷的位置與實(shí)際解剖鋼錠的結(jié)果基本一致。

3.2不同錠模斜度模擬結(jié)果對(duì)比

鋼錠錠身的剖面圖如圖4所示。

5.鄉(xiāng)鎮(zhèn)居民基本養(yǎng)老保險(xiǎn)和基本醫(yī)療保險(xiǎn)的繳費(fèi)。針對(duì)稅務(wù)部門收縮鄉(xiāng)鎮(zhèn)征收機(jī)構(gòu)，對(duì)鄉(xiāng)鎮(zhèn)居民基本養(yǎng)老保險(xiǎn)和基本醫(yī)療保險(xiǎn)的征繳帶來(lái)一定程度影響的現(xiàn)實(shí)情況，可考慮平移現(xiàn)有的征繳措施，如繼續(xù)委托村委會(huì)或街道辦進(jìn)行繳費(fèi)審核和確認(rèn)，繼續(xù)委托原來(lái)的商業(yè)銀行代扣繳社會(huì)保險(xiǎn)費(fèi)。商業(yè)銀行扣款成功后將扣款明細(xì)數(shù)據(jù)返回社會(huì)保險(xiǎn)經(jīng)辦機(jī)構(gòu)記賬，同時(shí)將各險(xiǎn)種的扣款總額向稅務(wù)部門申報(bào)入庫(kù)。

式中，D為鋼錠上端直徑；d為鋼錠下端直徑；H為錠身高度；φ為錠模斜度。

圖1　鋼錠凝固過(guò)程中液相分布

圖2　鋼錠凝固缺陷分布

圖3　35 t鋼錠局部疏松圖

圖4　鋼錠錠身的剖面圖

tan(φ)比值越小，鋼錠越“矮胖”，使得鋼錠變表面的傳熱面積變小，從而導(dǎo)致鋼錠的結(jié)晶速度慢，因此會(huì)產(chǎn)生較為嚴(yán)重的偏析現(xiàn)象。然而，此比值也不能過(guò)大，否則鋼液中的氣體和夾雜物比較難以排除，縮孔縮松會(huì)形成在錠身較深的部位，而且在這種情況下，由于鋼錠內(nèi)部存在有較大的靜壓力，會(huì)增大鋼錠縱向裂紋產(chǎn)生的幾率[4]。

據(jù)此，通過(guò)三維建模軟件SolidWorks建立幾何模型，表2列出模型模擬用的基本信息。此外建立了不同拔模斜度下的鋼錠模、底盤、冒口、絕熱板等幾何模型，并用SolidWorks將各個(gè)零件裝配起來(lái)。裝配完畢后，將每個(gè)零件單獨(dú)另存為STL格式文件。

表2　模擬模型基本信息

3.2.1溫度場(chǎng)分布

圖5為采用澆注溫度1 555℃，模具溫度為60℃，凝固時(shí)間分別是2 800 s、8 000 s、12 000 s，采用不同錠模斜度值下的溫度場(chǎng)分布對(duì)比切割圖?？梢钥闯觯撳V的凝固趨勢(shì)是由軸向從鋼錠底部向冒口推進(jìn)，徑向上從鋼錠四周模壁處向鋼錠的中心逐步進(jìn)行。

從圖5(a)可以看出，在凝固初期，隨著錠模斜度的增大，鋼錠內(nèi)部溫度場(chǎng)的分布差別不是很大，黑色區(qū)域類似于漏斗形狀。

圖5　不同錠模斜度值下的鋼錠溫度場(chǎng)切片圖

從圖5(b)可以看出，隨著凝固時(shí)間推移，鋼錠內(nèi)部溫度場(chǎng)有變化，隨著錠模斜度增大，白色區(qū)域面積明顯增加，說(shuō)明增大拔模斜度則內(nèi)部的溫度的推進(jìn)速度也減慢，鋼錠的凝固速度降低。尤其是φ=1°時(shí)，白色區(qū)域發(fā)生了分叉，說(shuō)明錠模斜度改變了凝固時(shí)溫度場(chǎng)的分布。

從圖5(c)可以看出，φ=1°時(shí)白色區(qū)域消失；φ=2°時(shí)白色區(qū)發(fā)生了明顯頸縮現(xiàn)象；φ=3°時(shí)白色區(qū)域也出現(xiàn)頸縮，相比沒(méi)有那么明顯。

3.2.2縮松縮孔分布

由于鋼錠存在液態(tài)收縮，因此在鋼錠的最后凝固區(qū)域如果無(wú)法得到外部鋼液的補(bǔ)充，就會(huì)產(chǎn)生縮孔。從橫切面角度觀察，軸線縮松在中心處，而從縱向剖面角度觀察則出現(xiàn)在軸線部位。根據(jù)溫度場(chǎng)分布情況，縮松出現(xiàn)在軸線處，與實(shí)際情況完全吻合。鋼液從上到下的溫度梯度小，在凝固后期，由于鋼錠軸心部等不到冒口及上部鋼液的補(bǔ)縮，冒口的補(bǔ)縮通道被過(guò)早地中斷，降低了冒口的補(bǔ)縮能力，導(dǎo)致疏松缺陷在此區(qū)域產(chǎn)生。不同φ值的鋼錠縮松縮孔切片圖如圖6所示?？梢钥闯觯S著拔模斜度值增大，鋼錠軸線區(qū)域內(nèi)縮松縮孔產(chǎn)生趨勢(shì)減小。

圖6　不同φ值的鋼錠縮松縮孔切片圖

綜合分析，鋼錠模斜度過(guò)大或過(guò)小均會(huì)使鋼錠產(chǎn)生嚴(yán)重且難以鍛合的縮松缺陷。35 t鋼錠模斜度為2°～3°之間比較合適。

3.3不同澆注速度模擬結(jié)果對(duì)比

表3列出模型模擬用的基本信息。錠模斜度為2°，澆注溫度為1 555℃，澆注速度分別為411 838 mm3/s、348 680 mm3/s、247 103 mm3/s。

不同澆注速度下縮松切片圖如圖7所示?？梢杂^察到，隨著澆注速度由411 838 mm3/s變?yōu)?48680mm3/s，縮松缺陷區(qū)域由鋼錠的底部上移；澆注速度為247 103 mm3/s時(shí)，鋼錠軸線區(qū)域的白色區(qū)域變大。這是由于澆注溫度一定時(shí)，鋼錠由下到上散熱逐漸加快，鋼錠底部橫向凝固速度快于縱向，使鋼錠下部的鋼液過(guò)早的被孤立，澆注速度過(guò)快，鋼液沒(méi)有及時(shí)起到補(bǔ)縮作用，鋼液下部就已經(jīng)達(dá)到固液兩相區(qū)，鋼液起不到補(bǔ)縮作用。如果澆注溫度一定，澆注速度過(guò)慢，此時(shí)就會(huì)導(dǎo)致整個(gè)凝固過(guò)程中補(bǔ)縮效果都不是很理想。

表3　模擬模型基本信息

圖7　不同澆注速度縮松切片圖

3.4不同澆注溫度模擬結(jié)果對(duì)比

表4列出模型模擬用的基本信息。錠模斜度為2°，澆注速度348 680 mm3/s，澆注溫度分別為1 545℃、1 555℃、1 565℃。

表4　模擬模型基本信息

在不同的澆注溫度下，鋼錠的凝固順序都是從側(cè)模內(nèi)壁和鋼錠低模向中心部位和上方移動(dòng)，隨著時(shí)間的推進(jìn)，鋼錠模側(cè)壁的冷卻能力增加的同時(shí)，底盤的冷卻能力慢慢減弱，從而致使未凝固區(qū)呈V字形[5]。當(dāng)澆注溫度過(guò)低，冒口區(qū)域的發(fā)熱劑與保溫磚不能及時(shí)發(fā)揮出保溫的作用，而低溫導(dǎo)致的鋼錠凝固過(guò)快，樹枝狀晶體組織在鋼液凝固中容易出現(xiàn)，從而使得鋼液補(bǔ)縮困難，在鋼錠的中軸線上易形成相關(guān)缺陷。此外，由于溫度過(guò)低致使鋼液粘度增大，阻礙了夾雜物和氣體的上浮，也易產(chǎn)生相應(yīng)缺陷[6]。相反，澆注溫度必須控制到恰到好處，如果澆注溫度過(guò)高，則鋼錠收縮量過(guò)大，鋼錠內(nèi)部易出現(xiàn)縮孔縮松，如圖8所示。

圖8　不同澆注溫度下縮松切片圖

4　結(jié)論

大型鍛件質(zhì)量好壞與鋼錠的質(zhì)量密切相關(guān)。利用有限元軟件THERCAST進(jìn)行數(shù)值模擬計(jì)算分析，為35 t鋼錠的錠模選擇和30CrNiMo8材料澆注工藝優(yōu)化提供了理論依據(jù)。

(1)有限元軟件THERCAST模擬計(jì)算缺陷分布的情況與實(shí)際較好吻合。

(2)錠模斜度。錠模斜度不是越大越好，也不是越小越好。通過(guò)模擬認(rèn)為35 t鋼錠錠模斜度為2°～3°之間比較合適。

(3)澆注工藝。采用低溫快澆，高溫慢澆。30NiCrMo8鋼錠合理的澆注溫度為1 555℃，澆注速度348 680 mm3/s。

(4)鋼錠設(shè)計(jì)。為使鋼錠所產(chǎn)生的縮松缺陷最小，且在鍛造工藝中容易鍛合，工藝生產(chǎn)過(guò)程中采用增大冒口的高度和直徑來(lái)增大補(bǔ)縮量，但利用率會(huì)降低。

[1]李振彥. 大型鍛造用鋼錠疏松缺陷的模擬與實(shí)驗(yàn)研究[D]. 太原: 太原科技大學(xué),2012:1-2.

[2]黃希祜. 鋼鐵冶金原理[M].北京:冶金工業(yè)出版社, 2011.

[3]張春宇. 72噸鋼錠的凝固過(guò)程數(shù)值模擬與優(yōu)化[D].鞍山：遼寧科技大學(xué), 2013:18-20.

[4]陳進(jìn)，余連權(quán).大型鋼錠模的設(shè)計(jì)[J].重慶科技學(xué)院學(xué)報(bào)(自然科學(xué)版)，2010,12(5)：116-118.

[5]更健，劉莊. 鋼模結(jié)構(gòu)及澆注工藝對(duì)鋼錠內(nèi)部縮松的影響[J]. 大型鑄鍛件，1997(3):15-17.

[6]金楊, 安紅萍.大型鋼錠凝固特性的初步研究[J]. 大型鑄鍛件，2011(1):5-8.

編輯杜青泉

Research of Numerical Simulation on Influence of Mold Angle and Pouring Process on Quality of Steel Ingot

Niu Liqun, Yuan Qi, Wang Xuming, Yang Desheng, Feng Tingbiao

By adopting THERCAST simulation software, the casting and solidification process of 35 t ingot has been simulated, so as to analyze the causes of internal defects of steel ingot. The optimal pouring process of 35 t 30CrNiMo8 steel ingot has been obtained by simulating and calculating the influence of mold angle, pouring temperature and pouring speed on defect distribution.

THERCAST; mold angle; pouring temperature; pouring speed

2016—07—01

甘肅省科技計(jì)劃資助項(xiàng)目(2015GS05896)。

牛立群(1986—)，男，助理工程師，研究生，主要從事鑄鍛工藝過(guò)程模擬及應(yīng)用研究。

TG316.1+92

B