30Cr2Ni4MoV低壓轉(zhuǎn)子加熱過程的計算機模擬

王子榮　黨淑娥　侯　微　劉海龍　秦尚武

(太原科技大學材料科學與工程學院，山西030024)

?

30Cr2Ni4MoV低壓轉(zhuǎn)子加熱過程的計算機模擬

王子榮黨淑娥侯微劉海龍秦尚武

(太原科技大學材料科學與工程學院，山西030024)

摘要：運用DEFORM-HT模塊以30Cr2Ni4MoV低壓轉(zhuǎn)子為對象研究其升溫過程中內(nèi)部產(chǎn)生的溫度場和熱應力，獲得了加熱過程中低壓轉(zhuǎn)子截面幾個關鍵位置處的升溫曲線及熱應力曲線，分析了最大截面上的溫度分布以及溫差對熱應力的影響，得出了加熱速度和保溫時間與溫差和熱應力之間的關系。

關鍵詞：30Cr2Ni4MoV鋼；低壓轉(zhuǎn)子；溫度場；熱應力；計算機模擬

大型低壓轉(zhuǎn)子的熱加工因其存在尺寸效應，所以在加熱過程中會受到傳熱和導熱等因素的影響。當加熱工藝參數(shù)設定不當時，轉(zhuǎn)子內(nèi)部溫度分布極不均勻，由此會引起轉(zhuǎn)子內(nèi)部不同部位在加熱過程奧氏體化的程度不一樣，即奧氏體晶粒大小分布不均勻，表層一定深度可能產(chǎn)生過熱甚至過燒，而心部可能未發(fā)生完全奧氏體化，會導致高溫狀態(tài)下轉(zhuǎn)子內(nèi)部力學性能不均勻。尤其轉(zhuǎn)子的半徑方向會產(chǎn)生明顯的溫度差，這直接導致了熱應力及組織應力的產(chǎn)生和累計，且溫差越大，兩者相互影響疊加后的綜合應力就越大。當綜合應力大于低壓轉(zhuǎn)子鋼在設定溫度下的屈服強度時就會造成工件變形，而大于抗拉強度則會造成斷裂[1～5]。由于近年來計算機熱處理模擬技術(shù)飛速發(fā)展產(chǎn)生了質(zhì)的飛躍，融入了奧氏體化相變動力學和傳熱學等，所以不用實驗也可以求出零件在熱處理過程中任一時刻的溫度場分布及熱應力分布[6～11]。所以，若工件結(jié)構(gòu)、加熱介質(zhì)和各種參數(shù)等詳細而準確，就能實現(xiàn)實際生產(chǎn)過程的精準計算[12]。

本文運用DEFORM-HT模塊以30Cr2Ni4MoV低壓轉(zhuǎn)子為對象研究其升溫過程中內(nèi)部產(chǎn)生的溫度場和熱應力。對加熱過程轉(zhuǎn)子截面幾個關鍵位置處的升溫曲線、最大截面上的溫度分布、轉(zhuǎn)子截面溫差對熱應力的影響進行了分析，從而根據(jù)模擬結(jié)果對轉(zhuǎn)子加熱過程所使用的工藝的合理性進行分析，為接下來改進加熱工藝提供了理論指導。

1加熱過程數(shù)值模型的建立

本文選取的低壓轉(zhuǎn)子見圖1，加熱工藝曲線如圖2，以兩者為基礎進行溫度場模擬，在此模擬基礎上將其模擬結(jié)果作為輸入?yún)?shù)進行瞬時應力-應變模擬。

圖1　工件尺寸圖

圖2　熱處理工藝曲線

汽輪機低壓轉(zhuǎn)子作為非常典型的軸對稱工件，其在柱狀坐標系下的溫度場方程為：

(1)

式中，ρ為材料密度；Cp為材料比定壓熱容；T為工件瞬態(tài)溫度；λ為材料導熱系數(shù)；t為過程進行的時間；q′為相變潛熱的熱流密度；r、z分別為沿徑向和軸向的坐標。

對瞬態(tài)溫度場進行求解時，既需要相應的邊界條件，也必須知道工件的初始條件。初始條件指的是所求解的非穩(wěn)態(tài)傳熱問題在初始時刻區(qū)域中各點的溫度，即工件裝爐時的溫度[13]。常用的換熱邊界條件為：

(2)

式中，H為邊界綜合換熱系數(shù)；TC為環(huán)境溫度，TW為低壓轉(zhuǎn)子與環(huán)境接觸邊界上溫度。

2DEFORM-HT有限元模擬

2.1材料熱物性參數(shù)的確定

低壓轉(zhuǎn)子鋼30Cr2Ni4MoV導熱系數(shù)λ與比熱容Cp均為溫度函數(shù)[14]，即

λ=33.44295+0.00224×T-1.57612E-5×T2

(馬氏體，室溫～700℃)

(3)

λ=118.23001-0.16649×T+745954E-5×T2

(奧氏體，800～1050℃)

(4)

Cp=0.5602-4.902×10-4T+1.407×10-6T2

(馬氏體，室溫～725℃)

(5)

Cp=2.746-3.350×10-3T+1.067×10-6T2

(奧氏體，770～900℃)

(6)

由上述公式可得出鋼的比熱容和導熱系數(shù)，見表1和表2。

表1　30Cr2Ni4MoV鋼不同溫度的導熱系數(shù)

表2　30Cr2Ni4MoV鋼比熱容隨溫度的變化

2.2模型的建立

圖1所示的低壓轉(zhuǎn)子為典型軸對稱階梯軸，所以選取最大直徑的1/4處為模擬對象，如圖3所示。圖3中P1、P2、P3分別表示工件徑向不同距離處的節(jié)點，P4點表示轉(zhuǎn)子棱邊，P1-P3表示心部至表面的三等劃分。

圖3　工件網(wǎng)格劃分和取點位置

2.3生產(chǎn)工藝和所需各類條件確定

大型鍛件由于尺寸效應，加熱過程中會產(chǎn)生很大的溫差，所以采用階梯緩慢加熱的方式來減小表層和心部的溫差。其加熱工藝曲線如圖2所示，即室溫(本文選取20℃)時將轉(zhuǎn)子裝入爐中，以30℃/h的加熱速度加熱至300℃、保溫10 h，再以17.5℃/h的加熱速度升溫到650℃、保溫35 h，接下來以16℃/h的加熱速度升溫到890℃、保溫75 h。

由于工件軸徑比大，體積大，所以一般情況下利用井式電爐吊掛加熱，工件的初始溫度就是其剛裝爐時的溫度(本文選取20℃)。邊界條件主要考慮工件和外界之間的熱交換過程，本文換熱邊界條件采用第三類邊界條件，主要是工件和加熱爐氣體的輻射與對流，換熱系數(shù)經(jīng)驗公式如下所示：

h=2.56(TW-TC)0.25+4.6×10-8×

(7)

h=3.26(TW-TC)0.25+4.6×10-8×

(8)

3模擬結(jié)果與分析

3.1溫度場模擬結(jié)果及分析

本文運用DEFORM-HT模塊，模擬可得低壓轉(zhuǎn)子徑向截面從表層至軸心不同深度幾個位置的升溫曲線，見圖4。從圖4可知，加熱過程中轉(zhuǎn)子表面及棱邊處溫度升高較快，且接近原設定工藝，棱邊和表層最大溫差為34℃，但棱邊達到指定溫度要比表層早54 h，表明越接近設定工藝溫度溫升越慢，尤其當從885℃升高到890℃時耗時34 h，所以設定工藝時在溫差允許范圍內(nèi)可適當提高溫度上限以大大縮減加熱時間。越靠近心部溫度升高越慢，其升溫過程已近似于一條斜線，其與設定工藝偏差較大，且越靠近表面處溫差越大，轉(zhuǎn)子表面與R/2處的溫度差值最高值為216℃，但是R/2與軸心處溫度差值最高值只有53℃，由此可知低壓轉(zhuǎn)子加熱過程中同一徑向截面不同位置的溫度差值并非等距縮減，從工件表層到工件軸心處溫度分布差異明顯。

圖4　轉(zhuǎn)子最大截面處特殊部位的升溫曲線

圖5為轉(zhuǎn)子加熱過程中棱邊、R/2、心部與表層的溫差。可以看出，隨著加熱溫度增大，轉(zhuǎn)子徑向溫度差值飛速加大，且停止加熱時其溫度差值達到巔峰。保溫時間越久轉(zhuǎn)子徑向溫度差值漸漸變小，停止保溫時溫度差值變?yōu)樽钚　Ｈ绠敿訜嶂?50℃時，表面與心部溫差達到最大值268℃，當650℃保溫35 h后，溫差減小至62℃，保溫前后共減小溫差達206℃。由此表明，階梯式加熱會有效減小大尺寸低壓轉(zhuǎn)子表面與心部溫差，且加熱過程階梯越多，表面與心部溫差越小。為了預測該轉(zhuǎn)子在以上加熱工藝條件下的熱應力，以溫度場模擬結(jié)果作為輸入?yún)?shù)，對轉(zhuǎn)子在加熱過程中的瞬時應力分布進行簡單分析，結(jié)果如圖6所示。將圖6與圖5對比可以看出，熱應力變化趨勢與溫差變化趨勢相對一致，但不完全一致。轉(zhuǎn)子棱邊、表層、R/2、心部的最大熱應力分別為15.3 MPa、21.6 MPa、23.6 MPa、23.4 MPa。棱邊和表層的最大熱應力發(fā)生在35 600 s處，由圖4可知，該處是棱邊和表層實際升溫速率最大值的終了處，而R/2點極大熱應力發(fā)生在159 000 s時，心部極大熱應力發(fā)生在61 400 s時，兩者均是實際加熱速率極大值的終了處。由此可得，最大熱應力并不是發(fā)生在最大溫差處而是發(fā)生在最大升溫速率終了處。由于低壓轉(zhuǎn)子一般要求其屈服應力為690 MPa～820 MPa，斷裂應力約為830 MPa，如果不考慮其加熱過程中的組織應力，單獨來看熱應力并不會致使其產(chǎn)生裂紋。

圖5　工件升溫過程中不同位置和表層之間溫差

圖6　轉(zhuǎn)子軸向應力與時間變化的關系曲線

4結(jié)論

(1)低壓轉(zhuǎn)子最大截面處即使以0.29℃/min速度加熱，也依然會造成表面與心部268℃的溫差。所以在實際生產(chǎn)中一味降低加熱速度并不能很好控制轉(zhuǎn)子截面溫度分布的均勻性。

(2)采用階梯式加熱可以有效地減小大截面低壓轉(zhuǎn)子的表面與心部溫差，得到較均勻的溫度分布。

(3)低壓轉(zhuǎn)子表面溫度在升溫過程中越接近設定溫度，升溫越緩慢，所以為了使工件盡快到達設定溫度，可將爐溫在誤差允許范圍內(nèi)上調(diào)，這樣可節(jié)省大量時間，縮短生產(chǎn)周期。

(4)低壓轉(zhuǎn)子最大截面上溫差造成的熱應力遠遠小于工件允許的變形及斷裂應力，以給定的加熱工藝曲線進行加熱，其最大截面上的溫差引起的熱應力并不會對轉(zhuǎn)子造成嚴重的變形及開裂。

參考文獻

[1]孫福民,徐濤,王年臣,等.汽輪機熱加工技術(shù)手冊[M].哈爾濱:哈爾濱汽輪機廠有限責任公司, 2006.

[2]康大韜,葉國斌.大型鍛件材料及熱處理[M].北京:龍門書局,1998.

[3]雷雪,徐薇平,韓利戰(zhàn).低壓轉(zhuǎn)子加熱過程模擬及工藝優(yōu)化[D]. 上海交通大學，2009(9):48-51.

[4]張振國.大直徑曲軸熱處理后的殘余應力分析[J].金屬熱處理，2005(08).

[5]羅利強.汽車輪轂軸承淬火殘余應力的建模與仿真[D].浙江工業(yè)大學， 2010.

[6]Inoue T,A rim o to K. Implementation of program for simulating heat treatment processes and some results[J]. Industrial Heating, 1995(4): 41-45.

[7]William D Disler. Using on-line predictive computer modeling to optimize heat treat processing[J].Industrial Heating, 1996(9): 51-56.

[8]趙雙群,林富生. 長期服役后高中壓轉(zhuǎn)子30Cr1Mo1V鋼的性能研究[J]. 動力工程學報， 2012(05).

[9]朱偉，鄧志成，申秀蘭，等.核電汽輪低壓轉(zhuǎn)子強度有限元計算[J].發(fā)電設備，2014(01).

[10]李世鍵，孫明月，李殿中，孫嫘，孫海燕. 百萬千瓦級核電轉(zhuǎn)子中心壓實過程的有限元模擬[J]. 熱加工工藝， 2011(21).

[11]任武,張趙威,吳運新,李益華.多體列車車軸淬火溫度場模擬和殘余應力研究[J].熱加工工藝，2014(18).

[12]吳淑珍，陳敬超，馮晶，等.鋼鐵淬火過程模擬的研究進展[J].材料導報，2008(6).

[13]匡琦，潘建生，葉建松.三維非線性有限元在熱處理爐CAD中的應用研究[J].金屬熱處理，2001(11).

[14]王曉燕，顧劍鋒，陳睿愷等.30Cr2Ni4MoV低壓轉(zhuǎn)子鋼中不同微觀組織的比熱和導熱系數(shù)的測定[J].金屬熱處理， 2011(4):1-3.

編輯杜青泉

Computer Simulation of Heating Process of

30Cr2Ni4MoV Low Pressure Rotor

Wang Zirong, Dang Shu′e, Hou Wei, Liu Hailong, Qin Shangwu

Abstract：By adopting DEFORM-HT module, the internally generated temperature field and thermal stress during the temperature rise period of 30Cr2Ni4MoV low pressure rotor have been studied. Then, the temperature rising curve and the thermal stress curve at several key points of cross section of rotor during the heating process have been obtained. Meanwhile, the distribution of temperature on the maximum cross section and the influence of thermal stress on temperature difference have been analyzed. Finally, the relationships between heating rate and holding time, temperature difference and thermal stress have been reached.

Key words：30Cr2Ni4MoV steel; low pressure rotor; temperature field; thermal stress; computer simulation

作者簡介：王子榮(1990—)，男，碩士研究生，主要研究方向：大型鍛造理論與新技術(shù)，熱處理數(shù)值模擬。電話：15035123077，E-mail：244706565@qq.com

基金項目：國家自然科學基金資助項目(51275330)；山西省自然科學基金資助項目(2012011022-4)。

收稿日期：2015—09—14

中圖分類號：O242.1

文獻標志碼：A