張宇航,吳日銘
(上海工程技術(shù)大學(xué) 材料工程學(xué)院,上海 201620)
模具鋼淬火冷卻[1]是一個(gè)復(fù)雜的過程,芯部的溫度變化一直是模具鋼熱處理工藝研究的焦點(diǎn)和難點(diǎn)。為預(yù)測(cè)芯部溫度變化以及由此帶來的組織和殘余應(yīng)力變化,熱處理的數(shù)值模擬也由此逐漸被廣泛應(yīng)用[2-5]。
工件尺寸遠(yuǎn)大于實(shí)驗(yàn)室模擬測(cè)試的試樣尺寸,文獻(xiàn)報(bào)道的模具鋼淬火多建立在小尺寸鋼塊試驗(yàn)研究(如端淬試驗(yàn)),而實(shí)際生產(chǎn)中,隨著工件尺寸增大及熱處理環(huán)境的改變,小尺寸工件試樣數(shù)值模擬研究結(jié)果和規(guī)律已無法與大截面尺寸工件實(shí)際情況相匹配[6]。
從淬火介質(zhì)狀態(tài)考慮,高溫工件淬入介質(zhì),介質(zhì)立即沸騰。以水為例,淬火開始時(shí),工件壁面的溫度遠(yuǎn)高于水的汽化溫度,水首先進(jìn)入膜態(tài)沸騰階段。隨著工件表面溫度的降低,淬火冷卻過程依次進(jìn)入過渡沸騰和核態(tài)沸騰階段;在膜態(tài)沸騰階段,工件表面覆蓋了一層完整氣膜阻止了表面與附近流體的直接接觸,熱量交換主要依靠表面與氣膜間的導(dǎo)熱和輻射進(jìn)行。隨著工件表面溫度的降低,過熱度減小,氣膜逐漸變薄直至破裂[7,8],沸騰程度與工件的表面積和表面粗糙度有關(guān),工件表面積越大,為氣泡的生長提供了越多的成核場(chǎng)所。
為準(zhǔn)確預(yù)測(cè)實(shí)際生產(chǎn)中大截面尺寸H13模具鋼在淬火過程中溫度場(chǎng)的變化情況,從優(yōu)化計(jì)算模型和模擬工業(yè)試驗(yàn)兩方面展開相關(guān)研究,即考慮淬火過程中由于沸騰狀態(tài)引起換熱改變,優(yōu)化計(jì)算模型,進(jìn)行溫度場(chǎng)、應(yīng)力場(chǎng)模擬。選用大截面尺寸鋼塊在生產(chǎn)現(xiàn)場(chǎng)完成淬火,實(shí)時(shí)監(jiān)測(cè)和記錄其內(nèi)部溫度變化,結(jié)合工件組織加以驗(yàn)證,以實(shí)現(xiàn)厚度較厚的模具鋼溫度場(chǎng)的準(zhǔn)確預(yù)測(cè),避免存在殘余應(yīng)力導(dǎo)致其開裂。
試驗(yàn)選用H13鋼,其成分含量如表1所示,選用工件的尺寸為φ300 mm×350 mm。按照?qǐng)D1(a)所示尺寸選取溫度監(jiān)測(cè)點(diǎn),并在測(cè)試點(diǎn)打孔,孔尺寸均為φ10 mm×175 mm。以便于熱電偶直接測(cè)得不同位置芯部溫度,為防止淬火時(shí)水進(jìn)入測(cè)溫孔內(nèi),在每個(gè)孔表面焊接耐高溫的P91無縫鋼管,如圖1(b)所示,保證整個(gè)鋼塊可以浸入淬火液中。
圖1 H13鋼工件
表1 H13模具鋼的化學(xué)成分 質(zhì)量分?jǐn)?shù)
1.2.1 淬火溫度實(shí)時(shí)監(jiān)測(cè)
為防止淬火時(shí)工件開裂,將工件進(jìn)行預(yù)退火處理,退火工藝如圖2(a)所示,退火工藝完成后,按照?qǐng)D2(b)所示工藝進(jìn)行淬火加熱和保溫處理。大型淬火池(12 000 mm×5 000 mm×5 600 mm)如圖3所示,圖3分別展現(xiàn)了入水前、入水過程中和入水后的工件狀態(tài)。從圖3(c)可以看到焊接的無縫鋼管可以保證在淬火冷卻過程中熱電偶正常工作。淬火過程中采用JTWARK-191F型熱電偶對(duì)材料芯部溫度變化進(jìn)行測(cè)量和監(jiān)測(cè),一端按厚度方向插入探測(cè)孔內(nèi)(175 mm),另一端與溫度記錄儀相連,數(shù)據(jù)采集頻率為2 s。
圖2 H13鋼熱處理工藝
圖3 淬火過程中工件狀態(tài)
1.2.2 數(shù)值模擬
利用數(shù)值分析和有限元理論,使用DeForm-3D熱處理模塊分析和模擬工件淬火過程中溫度場(chǎng)和應(yīng)力場(chǎng)變化狀態(tài)。根據(jù)工件實(shí)際尺寸建立模型并進(jìn)行網(wǎng)格劃分,通過試驗(yàn)、計(jì)算和查閱文獻(xiàn)獲取材料相關(guān)參數(shù)。并通過H13鋼相變特性參數(shù),使用JMATPRO軟件獲得H13鋼的CCT(過冷奧氏體連續(xù)冷卻轉(zhuǎn)變)曲線。
導(dǎo)熱微分方程是求取淬火過程中溫度場(chǎng)數(shù)據(jù)的常用方法,根據(jù)傅里葉定律以及能量守恒定律推導(dǎo)三維非穩(wěn)態(tài)導(dǎo)熱方程為[9-11]:
初始邊界條件為:
傳熱過程采用第三類邊界條件,公式為:
其中,λ為熱導(dǎo)率,W/(m?℃);T為溫度,Tw為材料表面溫度,Tf為外界溫度,℃;t為時(shí)間,s;ρ為材料密度,kg/m3;c為比熱容,J/(kg?K);h 為對(duì)流換熱系數(shù),W/(m2?℃);x,y,z為坐標(biāo)位置;n指向物體外法線方向。
由于工件淬火過程中受沸騰因素影響,產(chǎn)生了穿過蒸汽膜的輻射傳熱,并且輻射起著增大蒸汽膜厚度的作用,假定輻射和對(duì)流過程可以簡單相加則是不合理的。BROMLEY通過對(duì)水平表面上的膜態(tài)沸騰進(jìn)行研究,得出以下形式的超越方程計(jì)算總的換熱系數(shù)[12]。
其中,Ts為工件表面溫度,℃;Tsat為相應(yīng)壓力下淬火液的飽和溫度,℃;ε是固體發(fā)射率,其隨溫度變化如表2所示;σ是斯蒂芬-波爾玆曼常數(shù),其值為1.380 649× 10-23J/K。
表2 不同溫度下鋼的固體發(fā)射率ε
圖4所示為計(jì)算后得到總對(duì)流換熱系數(shù)與參考文獻(xiàn)[13]中對(duì)流換熱系數(shù)對(duì)比,改進(jìn)的對(duì)流換熱系數(shù)隨著工件表面溫度的升高先變大后變小,與參考文獻(xiàn)[13]中的趨勢(shì)一致。優(yōu)化后的對(duì)流換熱系數(shù)值比參考文獻(xiàn)[13]略高。
圖4 修正后的對(duì)流換熱系數(shù)與參考文獻(xiàn)[13]中的數(shù)據(jù)比較
通過數(shù)值模擬工件整個(gè)淬火過程,模擬所需的工件網(wǎng)格劃分為20 000個(gè),如圖(5)所示。通過溫度場(chǎng)云圖變化情況可以看到,從淬火開始,工件表面迅速降溫,但工件芯部依然是高溫狀態(tài),根據(jù)H13鋼的CCT曲線(見圖6),工件表面立刻進(jìn)行馬氏體轉(zhuǎn)變,芯部大部分奧氏體化,在大約750 s后工件內(nèi)部開始由外到內(nèi)逐漸進(jìn)行相轉(zhuǎn)變。由于工件尺寸大,在冷卻過程中,工件表面與芯部產(chǎn)生了較大的溫差,工件的表面與芯部在完成馬氏體轉(zhuǎn)變的時(shí)間差增大,且奧氏體向馬氏體轉(zhuǎn)變的體膨脹系數(shù)差別明顯,使淬火后的工件變形較大、殘余應(yīng)力較大,有開裂傾向[1]。
圖5 H13工件網(wǎng)格劃分及淬火過程中工件剖面溫度場(chǎng)變化
圖6 H13鋼CCT曲線
通過比較圖7(a)、(b)可知,由沸騰引起對(duì)流換熱系數(shù)的變化會(huì)改變淬火過程的溫度場(chǎng),通過優(yōu)化對(duì)流換熱系數(shù)后的模擬數(shù)據(jù),圖7(b)中曲線吻合程度明顯高于圖7(a),驗(yàn)證了優(yōu)化后的準(zhǔn)確性,同時(shí)也表明水淬過程熱輻射的作用不容忽視。觀察圖7(b)可以看到,4#位置位于工件邊緣處,其計(jì)算曲線與試驗(yàn)曲線吻合度最好;1#、3#位置分別處于距離工件2/3R和1/3R處,只有在前半段吻合度較好,但分別在淬火1 250 s和750 s后出現(xiàn)相對(duì)較大的誤差,可能與淬火馬氏體相變放熱有關(guān)。2#位置處于工件中心位置,發(fā)現(xiàn)其計(jì)算曲線整段未貼合實(shí)驗(yàn)曲線,這也是溫度預(yù)測(cè)的難點(diǎn)之一。
圖7 工件不同位置的冷卻曲線測(cè)量與計(jì)算結(jié)果比較
在降溫過程中,工件由外到內(nèi)溫度變化越來越平緩,冷卻速度越來越低。從淬火開始,工件芯部與表面的溫差逐漸增大,持續(xù)了約750 s,其中溫差最大達(dá)到860℃,但隨著淬火冷卻時(shí)間的延長,工件表面與芯部的溫差逐漸降低,在2 500 s時(shí),溫差小于100℃。
通過優(yōu)化對(duì)流換熱系數(shù)后,工件芯部的溫度趨勢(shì)預(yù)測(cè)仍存在明顯誤差,說明厚度較厚的H13鋼模塊芯部的溫度場(chǎng)的影響因素仍未完全揭示,如殘余拉/壓應(yīng)力對(duì)未轉(zhuǎn)變組織的溫度場(chǎng)影響、已轉(zhuǎn)變馬氏體組織導(dǎo)熱系數(shù)的變化等,還需要進(jìn)一步試驗(yàn)驗(yàn)證。
淬火過程中由于工件每個(gè)位置存在溫度差異和組織差異,工件內(nèi)部會(huì)產(chǎn)生較大的內(nèi)應(yīng)力。圖8所示為淬火過程中工件內(nèi)部4個(gè)探測(cè)點(diǎn)的應(yīng)力變化曲線,每個(gè)探測(cè)點(diǎn)都出現(xiàn)了應(yīng)力極值,特別是2#、3#、4#位置出現(xiàn)了明顯的極值點(diǎn),在750 s左右時(shí)分別達(dá)到了290、250、-200 MPa。在靠近芯部的2#、3#探測(cè)點(diǎn)先呈壓應(yīng)力后呈拉應(yīng)力,在表面位置上的4#探測(cè)點(diǎn)先呈拉應(yīng)力后呈壓應(yīng)力。這是因?yàn)樵诖慊鸪跗?,工件表面迅速冷卻收縮,而芯部由于溫度過高且散熱緩慢,產(chǎn)生不了相應(yīng)的收縮,從而表面因收縮受阻產(chǎn)生拉應(yīng)力,芯部產(chǎn)生對(duì)應(yīng)的壓應(yīng)力。隨著淬火的持續(xù)進(jìn)行,芯部開始馬氏體轉(zhuǎn)變,此時(shí)相變應(yīng)力起主導(dǎo)作用,并對(duì)外層壓應(yīng)力產(chǎn)生補(bǔ)償,芯部應(yīng)力轉(zhuǎn)變?yōu)槔瓚?yīng)力,表面4#位置因受到芯部的收縮而產(chǎn)生壓應(yīng)力,且逐漸減小,而應(yīng)力曲線中出現(xiàn)的極值是因淬火冷卻過程中組織應(yīng)力和熱應(yīng)力主導(dǎo)作用交替導(dǎo)致[14]。由于最大的拉/壓應(yīng)力產(chǎn)生在淬火750 s前后,可認(rèn)為該時(shí)刻是淬火開裂的關(guān)鍵時(shí)間點(diǎn),縮小在該時(shí)間點(diǎn)的芯部和表面拉壓應(yīng)力峰值差(約490 MPa)是避免淬火開裂的關(guān)鍵。
圖8 φ300mm×350mm H13鋼淬火過程內(nèi)應(yīng)力曲線
對(duì)淬火后的工件進(jìn)行取樣,分別取4#和2#孔底部位置進(jìn)行金相觀察,觀察結(jié)果如圖9(a)、(b)所示,工件表面組織主要為馬氏體,工件芯部組織主要為馬氏體和貝氏體及少量的殘余奧氏體。由表面到芯部馬氏體含量減少,貝氏體含量增加。因不同位置處冷卻速度不同,靠近表面位置冷卻速度快,自由能差值大,相變驅(qū)動(dòng)力大,過冷奧氏體主要轉(zhuǎn)變?yōu)轳R氏體組織;而芯部位置由于冷卻速度下降,自由能差值變小,過冷奧氏體的馬氏體轉(zhuǎn)變量減少,從而貝氏體轉(zhuǎn)變量增加[15]。
采用掃描電子顯微鏡對(duì)試樣芯部進(jìn)行觀察,如圖9(c)所示,發(fā)現(xiàn)在貝氏體周圍存在顆粒狀和塊狀碳化物且分布不均勻,可能由于降溫速率太小,使奧氏體向貝氏體轉(zhuǎn)變,碳化物沿著貝氏體鐵素體相界面析出導(dǎo)致;碳化物不僅量多且部分成排狀分布,這也是明顯降低沖擊韌性的影響因素[16]。
圖9 φ300 mm×350 mm H13鋼工件不同位置顯微組織照片
通過對(duì)尺寸為φ300 mm×350 mm的H13鋼工件進(jìn)行在線淬火監(jiān)測(cè),分別設(shè)置芯部、表面、1/3R和2/3R處的監(jiān)測(cè)點(diǎn),可準(zhǔn)確預(yù)測(cè)厚度較厚的H13鋼熱處理過程中的溫度場(chǎng)和淬火殘余應(yīng)力。在考慮鋼/水對(duì)流換熱系數(shù)受沸騰條件和熱輻射的影響的基礎(chǔ)上,對(duì)流換熱系數(shù)h進(jìn)行了優(yōu)化,得出了如下結(jié)果。
(1)各溫度監(jiān)測(cè)點(diǎn)的試驗(yàn)結(jié)果與計(jì)算結(jié)果相符,表明經(jīng)過優(yōu)化后的對(duì)流換熱系數(shù)可應(yīng)用于大截面尺寸工件溫度場(chǎng)的預(yù)測(cè)。
(2)將優(yōu)化好的數(shù)據(jù)導(dǎo)入計(jì)算機(jī)數(shù)值模擬應(yīng)力場(chǎng),結(jié)合顯微組織觀察發(fā)現(xiàn),淬火在750 s前后出現(xiàn)芯部、表面最大的拉壓峰值應(yīng)力差,證明大尺寸工件淬火存在淬火危險(xiǎn)期,危險(xiǎn)期將導(dǎo)致工件的變形和開裂,若材料的純凈度組織和性能足以抵抗淬火的拉壓峰值應(yīng)力沖擊,則淬火后期沒有開裂風(fēng)險(xiǎn)。