真空熱處理爐瞬態(tài)傳熱過(guò)程的數(shù)值模擬計(jì)算

2014-09-10 石磊解永蓉

　　建立了一個(gè)三維真空熱處理爐的非線性有限元模型，該模型的熱源是電流通過(guò)加熱體產(chǎn)生的焦耳熱，傳熱途徑主要考慮了加熱體與內(nèi)屏蔽層間、各屏蔽層之間的非線性輻射傳熱。利用有限元軟件ANSYS 對(duì)真空熱處理爐加熱過(guò)程中的瞬態(tài)溫度場(chǎng)進(jìn)行了模擬計(jì)算，得到了熱電偶溫度隨時(shí)間的變化曲線，并與實(shí)驗(yàn)結(jié)果進(jìn)行了比較，討論了進(jìn)一步提高計(jì)算準(zhǔn)確性的途徑。

　　真空焊接是在真空條件下對(duì)工件進(jìn)行加熱，此過(guò)程以輻射傳熱方式為主，真空加熱具有在別的介質(zhì)(大氣、可控氣氛、鹽浴)中加熱時(shí)不具備的特點(diǎn)。因?yàn)�，一般的加熱�?shí)際上不能夠在廣泛的溫度范圍和普遍的條件下都完全保持金屬熾熱表面與氣氛碳勢(shì)嚴(yán)格平衡和不起任何化學(xué)反應(yīng)(氧化、脫碳、增碳、浸蝕等)，而真空加熱是在極稀薄的氣氛中進(jìn)行的，避免了上述一般加熱的弊病。

　　零件真空焊接后質(zhì)量好壞取決于真空焊接工藝參數(shù)的選擇，如加熱溫度、保溫時(shí)間和真空度等。傳統(tǒng)的研究方法是采用實(shí)驗(yàn)來(lái)探索合適的真空焊接工藝參數(shù)，這種方法雖然簡(jiǎn)單易行，但卻需要進(jìn)行大量的實(shí)驗(yàn)，具有一定的盲目性，且耗費(fèi)大量的人力、物力和時(shí)間。通過(guò)數(shù)值模擬，我們可以建立真空焊接過(guò)程溫度場(chǎng)和真空焊接工藝參數(shù)的定量關(guān)系，并通過(guò)少量的驗(yàn)證性實(shí)驗(yàn)證明數(shù)值模擬方法在真空焊接過(guò)程中的適用性，大量的熱處理工藝參數(shù)篩選工作就可以在計(jì)算機(jī)上完成，它能夠顯著地提高效率，降低生產(chǎn)成本。這對(duì)正確地選擇真空焊接工藝參數(shù)，指導(dǎo)實(shí)際生產(chǎn)具有重要的現(xiàn)實(shí)意義。

　　本文將建立一個(gè)真空熱處理爐的三維有限元模型，并通過(guò)ANSYS 軟件對(duì)其輻射傳熱熱過(guò)程進(jìn)行數(shù)值模擬。

1、計(jì)算模型和方法

　　1.1、熱量傳遞的基本公式

　　熱量傳遞有三種基本方式，即熱傳導(dǎo)、對(duì)流和熱輻射。

　　考慮物體溫度隨時(shí)間而變(非穩(wěn)態(tài)問(wèn)題)，并且內(nèi)部具有熱源。在直角坐標(biāo)系下，熱傳導(dǎo)方程為：

熱傳導(dǎo)方程

　　式中T———物體的瞬態(tài)溫度，℃；k———材料的導(dǎo)熱系數(shù)，W/m·℃；c———材料的比熱容，J/kg·℃；ρ———材料的密度，kg/m3；qv———內(nèi)熱源強(qiáng)度，W/m3

　　對(duì)流邊界條件為：

對(duì)流邊界條件

　　式中h———對(duì)流換熱系數(shù)，W/m2·℃；Tf———環(huán)境溫度，℃

　　熱輻射是指由物體表面直接向外界發(fā)射可見(jiàn)和不可見(jiàn)射線，在空間傳遞能量的現(xiàn)象。熱輻射與導(dǎo)熱和對(duì)流不同，在傳遞能量時(shí)，不需要界面的相互接觸，并且是物體之間相互輻射。輻射傳熱的計(jì)算相對(duì)于傳導(dǎo)和對(duì)流來(lái)說(shuō)要更復(fù)雜。熱輻射的一個(gè)最重要的基本定律是斯蒂芬—波爾茲曼(Stefan—Boltzmann)定律：

　　式中T———物體表面溫度，K；Ta———周圍環(huán)境溫度，K；qr———輻射換熱的熱流密度，W/m2；ε———工件表面輻射率；σ———Stefan—Boltzmann常數(shù)

　　在兩個(gè)或多個(gè)表面之間的輻射換熱過(guò)程中，它不僅取決于參加輻射換熱物體的表面溫度、輻射性質(zhì)(如發(fā)射率、吸收率、反射率和投射率)、輻射的方向和波長(zhǎng)特性，而且還取決于表面的幾何形狀和表面間的空間位置。

　　分析表面熱輻射最為關(guān)鍵的任務(wù)之一是計(jì)算視角系數(shù)。所謂的視角系數(shù)是指由一個(gè)物體表面輻射和反射出去的全部熱射線有多大比例投射到另一物體的表面。角度系數(shù)決定于兩個(gè)面的形狀、大小及其相對(duì)位置和距離�？紤]如圖1 所示的表面1 和表面2 之間的輻射傳熱。

視角系數(shù)定義

圖1 視角系數(shù)定義

　　從表面1 向表面2 傳遞的輻射熱流強(qiáng)度q12可表示成：

　　其中T1 和T2 分別為表面1 和表面2 的溫度。F12 的物理意義是表面1 上被表面2 所看見(jiàn)的那一部分占表面1 總面積的比例。所以也被形象地稱為視角系數(shù)，其表達(dá)式為：

視角系數(shù)

　　從F12 的表達(dá)式中可以看到F12 完全取決于兩個(gè)進(jìn)行熱輻射表面的幾何和空間相對(duì)位置。目前常用Monte Calo 方法(蒙特卡洛方法)計(jì)算熱輻射視角系數(shù)。其基本思想就是基于光線發(fā)射和截?cái)嘣�，將從表�? 上隨機(jī)發(fā)射并被表面2 截?cái)嗟墓饩€占總光線的比例等效成視角系數(shù)F12。

　　1.2、有限元模型的基本假設(shè)

　　實(shí)際的真空熱處理爐結(jié)構(gòu)非常復(fù)雜，通常主要由爐殼、爐體、爐蓋、爐底和加熱元件等組成，其中爐殼是采用雙層不銹鋼板焊接而成的殼體，中間通水冷卻。爐體采用全金屬隔熱反射屏，有多層，最內(nèi)層為鉬板，次內(nèi)層為高溫不銹鋼，其余為普通不銹鋼。加熱元件采用鉬絲制成的“U”形狀電熱元件。在爐壁中部安裝1 支鎧裝監(jiān)控?zé)犭娕�，用于監(jiān)控爐溫。真空熱處理爐溫定義為監(jiān)控熱電偶的溫度。在實(shí)際真空焊接過(guò)程的各個(gè)階段，工件溫度、監(jiān)控?zé)犭娕紲囟取⒓訜嵩䴗囟�、各層隔熱屏溫度和爐殼溫度都是不同的。在整個(gè)過(guò)程中真空熱處理爐體內(nèi)主要通過(guò)熱輻射方式進(jìn)行傳熱。

　　ANSYS 有限元軟件已經(jīng)被廣泛應(yīng)用于傳熱過(guò)程的數(shù)值模擬，其技術(shù)也日趨成熟和完善，ANSYS 中關(guān)于輻射的重要假設(shè)和方法：

　　● 認(rèn)為輻射是平面現(xiàn)象，因此適合用不透明平面建模。

　　● 不直接計(jì)入平面反射率。考慮到效率，假設(shè)平面吸收率和發(fā)射率相等(a=ε)，因此，只有發(fā)射率特性需要在ANSYS 輻射分析中定義。

　　● ANSYS 不自動(dòng)計(jì)入發(fā)射率的方向特性，也不允許發(fā)射率定義隨波長(zhǎng)變化。

　　● 所有分隔輻射面的介質(zhì)在計(jì)算輻射能量交換時(shí)都看作是非參與的(不吸收也不發(fā)射能量) 。在現(xiàn)階段對(duì)真空熱處理爐的完整建模是困難的，必須抓住重點(diǎn)，對(duì)復(fù)雜的真空熱處理爐進(jìn)行簡(jiǎn)化，建立的有限元計(jì)算模型主要簡(jiǎn)化和假設(shè)條件如下：

　　● 真空熱處理爐的模型由金屬外殼，5 層隔熱屏，24 根加熱鉬絲組成，用小的金屬塊模擬為熱電偶行為；

　　● 外殼上采用水冷，因此用水對(duì)流邊界；

　　● 隔熱屏及加熱區(qū)靠熱輻射傳熱；

　　●真空熱處理爐各種材料是各向同性且均質(zhì)的；

　　●爐中熱源是電流通過(guò)加熱鉬絲產(chǎn)生的焦耳熱，計(jì)算中用測(cè)量得到加熱元件的體熱產(chǎn)生率(用電流和電壓折算)方法來(lái)控制熱量產(chǎn)生。

　　真空熱處理爐有限元計(jì)算模型具體參數(shù)為：爐體反射屏高為650 mm，內(nèi)徑470 mm，外徑550 mm，每層間距約15 mm，屏厚1 mm。加熱絲高500 mm，鼠籠直徑360 mm。

　　經(jīng)過(guò)簡(jiǎn)化，建立的有限元模型及網(wǎng)格劃分如圖2 所示。

有限元模型及網(wǎng)格劃分