介紹了風洞熱風閥在風洞系統(tǒng)中承受高溫及高壓的工況條件及其相關(guān)的技術(shù)要求。建立和描述了風洞閥門的物理及數(shù)學模型，確定了模型的邊界條件。分別采用有限容積法及有限元法對風洞閥的流場、溫度場及其應力應變場進行了耦合計算。給出了閥內(nèi)空氣最大流速、閥座表面平均溫度、閥體局部點最大應力的計算過程和結(jié)果及其閥體加工時所采用的工藝規(guī)則。

1、概述

　　風洞設(shè)備及其相關(guān)的設(shè)計制造技術(shù)是國防尖端科技的基礎(chǔ)，也是世界各國國防工業(yè)的核心技術(shù)。風洞熱風閥安裝在風洞儲能裝置加熱室的出口，是一種可承受高溫、高壓的截止性質(zhì)的閥門。當風洞儲能裝置內(nèi)空氣壓力和溫度達到一定數(shù)值時，風洞閥快速開啟，并在閥后的拉伐爾噴管內(nèi)形成高超聲速氣流，以此沖擊發(fā)動機進行實驗。在風洞實驗過程中，高溫高壓空氣流經(jīng)閥體，閥體承受著高壓及溫差應力的作用，對閥體結(jié)構(gòu)的合理設(shè)計關(guān)系著風洞系統(tǒng)的安全性。

　　傳統(tǒng)的高參數(shù)閥門多采用經(jīng)驗公式法進行設(shè)計，由于經(jīng)驗公式法涉及到的變量較少，很難對閥體結(jié)構(gòu)進行定量優(yōu)化。數(shù)值計算法可對閥門的流態(tài)、熱態(tài)、應力態(tài)等多種狀態(tài)進行模擬分析，從而得到優(yōu)化的結(jié)構(gòu)，但涉及到的流固熱多物理場耦合如何實現(xiàn)還有待探究。

2、多物理場耦合數(shù)學模型

　　熱風閥主體結(jié)構(gòu)為鋼材，為了隔熱在閥體內(nèi)壁襯有隔熱材料。閥座密封面處采用熱噴涂工藝，基材為閥體用鋼材，用水冷方式降低閥座溫度。

　　分析閥體受熱過程可知，閥內(nèi)流動的是潔凈的高溫空氣，可不考慮氣體向閥體內(nèi)壁的熱輻射。高溫空氣以對流形式，將熱量傳遞給隔熱材料，隔熱材料和閥體以熱傳導的形式將熱量傳至閥體表面，最終以對流方式將熱量傳給環(huán)境大氣。由于熱風閥閥體內(nèi)空氣流速不同，導致閥體各部位換熱程度不同，從而引起溫差應力。因而閥體的安全性問題是流動、傳熱及應力等物理場耦合的結(jié)果。準確計算流動狀態(tài)及合理處理邊界層傳熱是關(guān)鍵。

4、計算結(jié)果分析

4.1、流場及溫度場分布

　　圖2 給出了閥體內(nèi)流場分布的狀態(tài)圖，圖中可以看出，空氣在閥腔內(nèi)速度分布并不均勻，最大流速達到9.5m /s，最大流速出現(xiàn)在拐角流線發(fā)生急變處。

圖2 穩(wěn)態(tài)下閥體內(nèi)流速分布

　　圖3 與圖4 分別給出了閥體表面及閥中截面的溫度分布，從閥體表面溫度分布圖中可看出閥表面溫度分布是不均勻的，最高溫度為71.6℃，出現(xiàn)在閥體的拐角處，此處的高溫區(qū)與流場的高速區(qū)相對應，產(chǎn)生這一現(xiàn)象的原因在于流體的正面沖擊邊界層厚度減薄，閥內(nèi)壁換熱增強，繼而引起閥外壁溫度跟著升高。

圖3 閥體外表面溫度場分布

圖4 閥體中截面溫度分布

　　圖5 給出了閥座水冷處的局部溫度分布圖，在密封面處平均溫度約為110℃，但在這一區(qū)域出現(xiàn)了2個高溫點( 圖6) ，這個小區(qū)域溫度與周邊溫度相差很大，即在這一區(qū)域有很大的溫度梯度。這2 個高溫點均位于保溫材料、鋼材、高溫流體相接觸的位置。產(chǎn)生這一現(xiàn)象的原因在于，空氣在流動時先流經(jīng)保溫層，由于保溫材料熱導率很小，在保溫層內(nèi)表面的溫度很高( 接近900℃) ，與保溫層相接觸的流體邊界層內(nèi)的溫度也很高，這部分高溫空氣與鋼材相接觸時，導致鋼材的溫度也升高。隨著流體繼續(xù)向前流動，鋼材帶走了邊界層內(nèi)的熱量，邊界層內(nèi)流體溫度快速下降，鋼材的溫度也快速下降，由于空氣的比熱很小，而且熱導率也很小，因而在這點溫度較高，且產(chǎn)生很大溫度梯度。這一計算結(jié)果也提示我們，在設(shè)計流道時應盡量避免流道內(nèi)凸臺等結(jié)構(gòu)，盡可能使邊界層穩(wěn)定，以提高邊界層的隔熱作用。

4.2、應力與變形分析

　　將溫度場各點數(shù)據(jù)導入閥體結(jié)構(gòu)單元中，對閥體進行應力分析，得到閥體等效分布如圖7 所示。閥體上最大等效應力為284.67MPa，分布在閥體內(nèi)部凸臺上側(cè)。將應力最大點的分布位置局部放大( 圖8) ，可見閥體的高應力分布區(qū)域很小，這是由于溫差和結(jié)構(gòu)突變的共同作用導致該位置的應力值突增，如果將該位置的銳角倒鈍，則該處的應力集中情況將會大幅改善。

圖5 穩(wěn)態(tài)下閥體水冷面溫度分布

圖6 高溫點產(chǎn)生位置

圖7 閥體等效應力分布

圖8 閥體最大等效應力分布位置

　　閥體外壁等效應力分布如圖9 所示，外壁平面的中間位置等效應力相對較大，可達105.02MPa。最小等效應力分布在閥體下角，其值為1.028MPa。閥體總變形如圖10 所示，最大總變形量為0.5mm，位于閥體下方直邊棱處，最小總變形量為0.03mm，位于上端法蘭端面處。

圖9 閥體外壁等效應力分布

圖10 閥體總變形分布

5、結(jié)語

　　(1) 建立了用于閥體流場、溫度場、應力應變場的物理數(shù)學模型，并將有限容積法與有限元法有效結(jié)合在一起，求解了風洞閥內(nèi)流體與固體區(qū)域的各物理場分布。

　　(2) 流場與溫度場計算表明，閥體外表溫度小于75℃，閥座表面平均溫度110℃，能夠達到設(shè)計要求。但設(shè)計中應注意流體通道的光滑，避免突臺等不規(guī)則的結(jié)構(gòu)導致邊界層破壞引起局部高溫。

　　(3) 由于隔熱材料起到了作用，閥體總體的應力與應變均不大，但由于結(jié)構(gòu)突變在局部點有較高的應力值。對于這類小區(qū)域，尖角引起的應力集中問題，可在閥體加工時采用銳角倒鈍的工藝措施。