真空熱處理爐瞬態(tài)傳熱過程的數(shù)值模擬計(jì)算
建立了一個(gè)三維真空熱處理爐的非線性有限元模型，該模型的熱源是電流通過加熱體產(chǎn)生的焦耳熱，傳熱途徑主要考慮了加熱體與內(nèi)屏蔽層間、各屏蔽層之間的非線性輻射傳熱。利用有限元軟件ANSYS 對(duì)真空熱處理爐加熱過程中的瞬態(tài)溫度場(chǎng)進(jìn)行了模擬計(jì)算，得到了熱電偶溫度隨時(shí)間的變化曲線，并與實(shí)驗(yàn)結(jié)果進(jìn)行了比較，討論了進(jìn)一步提高計(jì)算準(zhǔn)確性的途徑。
真空焊接是在真空條件下對(duì)工件進(jìn)行加熱，此過程以輻射傳熱方式為主，真空加熱具有在別的介質(zhì)(大氣、可控氣氛、鹽浴)中加熱時(shí)不具備的特點(diǎn)。因?yàn)�，一般的加熱�?shí)際上不能夠在廣泛的溫度范圍和普遍的條件下都完全保持金屬熾熱表面與氣氛碳勢(shì)嚴(yán)格平衡和不起任何化學(xué)反應(yīng)(氧化、脫碳、增碳、浸蝕等)，而真空加熱是在極稀薄的氣氛中進(jìn)行的，避免了上述一般加熱的弊病。
零件真空焊接后質(zhì)量好壞取決于真空焊接工藝參數(shù)的選擇，如加熱溫度、保溫時(shí)間和真空度等。傳統(tǒng)的研究方法是采用實(shí)驗(yàn)來探索合適的真空焊接工藝參數(shù)，這種方法雖然簡(jiǎn)單易行，但卻需要進(jìn)行大量的實(shí)驗(yàn)，具有一定的盲目性，且耗費(fèi)大量的人力、物力和時(shí)間。通過數(shù)值模擬，我們可以建立真空焊接過程溫度場(chǎng)和真空焊接工藝參數(shù)的定量關(guān)系，并通過少量的驗(yàn)證性實(shí)驗(yàn)證明數(shù)值模擬方法在真空焊接過程中的適用性，大量的熱處理工藝參數(shù)篩選工作就可以在計(jì)算機(jī)上完成，它能夠顯著地提高效率，降低生產(chǎn)成本。這對(duì)正確地選擇真空焊接工藝參數(shù)，指導(dǎo)實(shí)際生產(chǎn)具有重要的現(xiàn)實(shí)意義。
本文將建立一個(gè)真空熱處理爐的三維有限元模型，并通過ANSYS 軟件對(duì)其輻射傳熱熱過程進(jìn)行數(shù)值模擬。
1、計(jì)算模型和方法
1.1、熱量傳遞的基本公式
熱量傳遞有三種基本方式，即熱傳導(dǎo)、對(duì)流和熱輻射。
考慮物體溫度隨時(shí)間而變(非穩(wěn)態(tài)問題)，并且內(nèi)部具有熱源。在直角坐標(biā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ì)流邊界條件為：
式中h———對(duì)流換熱系數(shù)，W/m2·℃；Tf———環(huán)境溫度，℃
熱輻射是指由物體表面直接向外界發(fā)射可見和不可見射線，在空間傳遞能量的現(xiàn)象。熱輻射與導(dǎo)熱和對(duì)流不同，在傳遞能量時(shí)，不需要界面的相互接觸，并且是物體之間相互輻射。輻射傳熱的計(jì)算相對(duì)于傳導(dǎo)和對(duì)流來說要更復(fù)雜。熱輻射的一個(gè)最重要的基本定律是斯蒂芬—波爾茲曼(Stefan—Boltzmann)定律：
式中T———物體表面溫度，K；Ta———周圍環(huán)境溫度，K；qr———輻射換熱的熱流密度，W/m2；ε———工件表面輻射率；σ———Stefan—Boltzmann常數(shù)
在兩個(gè)或多個(gè)表面之間的輻射換熱過程中，它不僅取決于參加輻射換熱物體的表面溫度、輻射性質(zhì)(如發(fā)射率、吸收率、反射率和投射率)、輻射的方向和波長(zhǎng)特性，而且還取決于表面的幾何形狀和表面間的空間位置。
分析表面熱輻射最為關(guān)鍵的任務(wù)之一是計(jì)算視角系數(shù)。所謂的視角系數(shù)是指由一個(gè)物體表面輻射和反射出去的全部熱射線有多大比例投射到另一物體的表面。角度系數(shù)決定于兩個(gè)面的形狀、大小及其相對(duì)位置和距離�？紤]如圖1 所示的表面1 和表面2 之間的輻射傳熱。
圖1 視角系數(shù)定義
從表面1 向表面2 傳遞的輻射熱流強(qiáng)度q12可表示成：
其中T1 和T2 分別為表面1 和表面2 的溫度。F12 的物理意義是表面1 上被表面2 所看見的那一部分占表面1 總面積的比例。所以也被形象地稱為視角系數(shù)，其表達(dá)式為：
從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é)犭娕迹糜诒O(jiān)控爐溫。真空熱處理爐溫定義為監(jiān)控熱電偶的溫度。在實(shí)際真空焊接過程的各個(gè)階段，工件溫度、監(jiān)控?zé)犭娕紲囟取⒓訜嵩䴗囟�、各層隔熱屏溫度和爐殼溫度都是不同的。在整個(gè)過程中真空熱處理爐體內(nèi)主要通過熱輻射方式進(jìn)行傳熱。
ANSYS 有限元軟件已經(jīng)被廣泛應(yīng)用于傳熱過程的數(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ì)的；
●爐中熱源是電流通過加熱鉬絲產(chǎn)生的焦耳熱，計(jì)算中用測(cè)量得到加熱元件的體熱產(chǎn)生率(用電流和電壓折算)方法來控制熱量產(chǎn)生。
真空熱處理爐有限元計(jì)算模型具體參數(shù)為：爐體反射屏高為650 mm，內(nèi)徑470 mm，外徑550 mm，每層間距約15 mm，屏厚1 mm。加熱絲高500 mm，鼠籠直徑360 mm。
經(jīng)過簡(jiǎn)化，建立的有限元模型及網(wǎng)格劃分如圖2 所示。
圖2 有限元模型及網(wǎng)格劃分
2、計(jì)算結(jié)果
在加熱鉬絲上加載如圖3 的加熱功率曲線后，得到金屬塊(用于模擬熱電偶)上計(jì)算的溫度隨時(shí)間的變化，作為對(duì)比，圖上給出了測(cè)量得到的溫度隨時(shí)間的變化曲線。
圖3 計(jì)算和測(cè)量的溫度隨時(shí)間變化曲線
圖4 顯示了在保溫階段加熱鉬絲、金屬隔熱屏各層間的溫度分布情況。從圖中可以看出，由于熱傳導(dǎo)，加熱鉬絲、各金屬隔熱屏上的溫度是均勻的，但鉬絲上溫度比其它各部分的溫度高得多，這也是熱輻射的需要。
圖4 有限元法計(jì)算的溫度分布云圖
3、結(jié)論
計(jì)算和測(cè)量溫度的變化范圍和趨勢(shì)基本保持一致，說明計(jì)算模型和計(jì)算過程基本是正確的，可以為實(shí)際生產(chǎn)提供理論指導(dǎo)。
下一步的工作是提高計(jì)算的準(zhǔn)確度和建立更符合實(shí)際的計(jì)算模型，其中包括：
● 平滑加熱絲兩端的電壓電流和加熱絲功率的準(zhǔn)確確定。
● 各種材料輻射率的準(zhǔn)確確定。
● 隔熱屏的厚度會(huì)影響熱容量，導(dǎo)致溫度延遲，要求有準(zhǔn)確值。
● 考慮其它形式的熱損失(如連接部分的熱傳導(dǎo)引起的熱損失，約10%~15%左右)。
● 數(shù)值計(jì)算的精度和收斂性。
