郭政波,魏智輝,潘鵬飛
(中國飛行試驗(yàn)研究院發(fā)動(dòng)機(jī)所,陜西?西安?710089)
提高渦輪效率是提高航空發(fā)動(dòng)機(jī)效率、降低飛行成本和減少航空排氣污染的重要策略之一。研究表明,高壓渦輪的間隙每減小0.254 mm,可以使得燃油消耗率降低1%,從而顯著提高發(fā)動(dòng)機(jī)壽命。此外,還能減少NOx,CO和CO2的排放[1]。渦輪葉尖間隙主動(dòng)控制技術(shù)的應(yīng)用可以提高發(fā)動(dòng)機(jī)效率、減小燃油消耗率、增加部件壽命[2,3]。出于減少燃油消耗、降低飛行使用成本和環(huán)保的需求,渦輪葉尖間隙控制技術(shù)成為研究熱點(diǎn)。在控制過程中,渦輪葉尖間隙過大會(huì)導(dǎo)致效率低下,渦輪葉尖間隙過小又容易引起擦碰等情況,影響發(fā)動(dòng)機(jī)安全。鑒于渦輪葉尖間隙控制技術(shù)可以帶來巨大的收益,國內(nèi)外學(xué)者對(duì)該技術(shù)進(jìn)行了深入研究。
Kypuros等研究了預(yù)測(cè)渦輪葉尖間隙變化的動(dòng)態(tài)模型,該模型考慮發(fā)動(dòng)機(jī)工作時(shí)熱應(yīng)力、轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)產(chǎn)生的離心力和機(jī)匣內(nèi)外壓差等因素對(duì)間隙變化的影響,并應(yīng)用于飛機(jī)起飛時(shí)渦輪葉尖間隙的計(jì)算[4]。Lattime等設(shè)計(jì)了液壓驅(qū)動(dòng)的機(jī)械系統(tǒng),并開展了試驗(yàn)評(píng)估它的響應(yīng)和定位準(zhǔn)確性,該系統(tǒng)可以用于進(jìn)一步開發(fā)渦輪葉尖間隙主動(dòng)控制系統(tǒng)的作動(dòng)器[5]。賈丙輝研究了轉(zhuǎn)子對(duì)高壓渦輪葉尖間隙變化規(guī)律的影響,討論了轉(zhuǎn)子在飛行器機(jī)動(dòng)飛行情況下的振動(dòng)幅值對(duì)葉尖間隙的影響[6]。楊家禮等設(shè)計(jì)了一種冷氣自動(dòng)調(diào)節(jié)機(jī)構(gòu),研究了一種渦輪葉尖間隙控制技術(shù),證明了該方法對(duì)渦輪葉尖間隙可有效控制[7]。
綜上可以看到,國內(nèi)外對(duì)渦輪葉尖間隙的研究主要集中于對(duì)渦輪葉尖間隙形成的機(jī)理、葉尖間隙的仿真計(jì)算、葉尖間隙的控制技術(shù)、葉尖間隙對(duì)渦輪性能影響的研究。本文從渦輪葉尖間隙的變化機(jī)理出發(fā),建立了機(jī)匣、葉片和轉(zhuǎn)子的簡化模型,并應(yīng)用某型發(fā)動(dòng)機(jī)的試車數(shù)據(jù)計(jì)算了發(fā)動(dòng)機(jī)渦輪葉尖間隙,研究地面起動(dòng)時(shí)發(fā)動(dòng)機(jī)渦輪葉尖間隙動(dòng)態(tài)變化規(guī)律,為發(fā)動(dòng)機(jī)地面起動(dòng)時(shí)渦輪葉尖間隙控制策略的設(shè)計(jì)提供參考。
葉尖間隙變化的主要原因是發(fā)動(dòng)機(jī)靜止或轉(zhuǎn)動(dòng)部件由于大量負(fù)載而發(fā)生位移或者變形。負(fù)載主要包括兩種:發(fā)動(dòng)機(jī)負(fù)載和飛行負(fù)載。發(fā)動(dòng)機(jī)負(fù)載包括離心力、熱應(yīng)力、發(fā)動(dòng)機(jī)內(nèi)壓和推力負(fù)載。飛行負(fù)載包括慣性力(重力)、空氣動(dòng)力學(xué)力(外部環(huán)境壓力)和陀螺負(fù)載。如圖1所示,發(fā)動(dòng)機(jī)負(fù)載可產(chǎn)生對(duì)稱的和非對(duì)稱的間隙變化,飛行負(fù)載引起非對(duì)稱的間隙變化。
圖1 ??兩種葉尖間隙變化
一般最小間隙發(fā)生在飛機(jī)起飛和發(fā)動(dòng)機(jī)起動(dòng)過程中,在這種情況下,由于離心力和葉片被快速加熱,轉(zhuǎn)子組合體膨脹很快。與此同時(shí),轉(zhuǎn)子周圍的機(jī)匣結(jié)構(gòu)由于熱效應(yīng)膨脹的速度相對(duì)較慢,這就導(dǎo)致了間隙的迅速減小。最終,機(jī)匣的增長速率超過轉(zhuǎn)子組件的增長,間隙稍有增大。在這一過程中,為了避免摩擦,渦輪機(jī)設(shè)計(jì)時(shí)必須預(yù)留間隙。這一額外的間隙導(dǎo)致非最佳間隙的產(chǎn)生,并使得發(fā)動(dòng)機(jī)在大部分正常工作狀態(tài)下燃油油耗增大。
渦輪葉尖間隙模型包含3個(gè)基本單元:機(jī)匣、渦輪轉(zhuǎn)盤、渦輪葉片。由于熱應(yīng)力、壓力分布所帶來的位置偏移動(dòng)態(tài)疊加,導(dǎo)致渦輪葉尖間隙的變化,所以有:
式中,rshroud(t)、rrotor(t)、lblade(t)分別為機(jī)匣內(nèi)徑、轉(zhuǎn)子盤外徑、渦輪葉片長度。ra、ro、L是機(jī)匣、轉(zhuǎn)子盤、渦輪葉片的初始長度,u為由于熱和機(jī)械力產(chǎn)生的動(dòng)態(tài)形變量。下面對(duì)3個(gè)單元分別建模計(jì)算。
機(jī)匣單元:為方便計(jì)算,可將復(fù)雜的機(jī)匣結(jié)構(gòu)簡化為圓環(huán)結(jié)構(gòu),如圖2所示。其內(nèi)表面涂著一層耐磨材料,同時(shí)也起著隔熱層的作用。襯套內(nèi)徑暴露在加熱的燃?xì)庵校摐囟燃s等于渦輪進(jìn)口溫度Tturbine。外表面暴露在壓氣機(jī)出口氣流中,假定該溫度等于壓氣機(jī)放氣溫度Tcompressor,假定渦輪葉尖間隙區(qū)域渦輪進(jìn)口和出口壓力是線性的,壓氣機(jī)放氣流過了襯套和機(jī)匣的空間。忽略壓氣機(jī)出口和機(jī)匣外表面溫度差異。由于傳導(dǎo)特性,耐磨層消散了很多渦輪熱量,降低了合金層rc所要承受的溫度Tc。Tc是熱邊界狀態(tài)合金的內(nèi)表面用來計(jì)算由于熱力的偏差。為簡化模型,耐熱層內(nèi)表面和合金層外表面假定是準(zhǔn)穩(wěn)態(tài)的。耐磨層內(nèi)表面壁溫Ta熱轉(zhuǎn)換效率ha,機(jī)匣外表面溫度Tb熱轉(zhuǎn)換效率hb。經(jīng)過以上簡化,經(jīng)過推導(dǎo)可得由于熱力導(dǎo)致的機(jī)匣徑向變化量[8]:
式中,α為熱膨脹系數(shù)。由于壓力分布導(dǎo)致的機(jī)匣在結(jié)合處的徑向rc變化量為:
式中:E為楊氏模量, 為泊松比。
圖2 ??機(jī)匣單元結(jié)構(gòu)簡圖
轉(zhuǎn)子盤:轉(zhuǎn)子盤模型結(jié)構(gòu)如圖3所示,在建模時(shí)做了簡化,視轉(zhuǎn)子盤為薄厚均勻的圓盤,只考慮離心力和壓入空氣溫度導(dǎo)致的變化。壓氣機(jī)放氣通常用來清洗和冷卻圍繞轉(zhuǎn)子旁邊的空腔,只有很小的尖部表面暴露在高溫核心氣流中。假定熱膨脹只是由轉(zhuǎn)子盤和壓氣機(jī)放氣溫度的熱交換造成的,為了簡化分析,假定轉(zhuǎn)子與葉片無熱交換。壓氣機(jī)放氣流過大多是轉(zhuǎn)子盤表面積。由于轉(zhuǎn)子盤的表面暴露在壓氣機(jī)放氣溫度中,假定轉(zhuǎn)子盤的溫度等于壓氣機(jī)放氣溫度。轉(zhuǎn)子由于熱應(yīng)力導(dǎo)致的變化量為:
式中,α為熱膨脹系數(shù),r0為轉(zhuǎn)子無熱應(yīng)力時(shí)的半徑。
離心力作用下的轉(zhuǎn)子徑向形變量可估算為:
圖3 ??轉(zhuǎn)子盤模型結(jié)構(gòu)
渦輪葉片:渦輪葉片模型結(jié)構(gòu)如圖4所示。葉片在工作中既受到離心力作用,又受到熱膨脹的作用。渦輪葉片相對(duì)較薄,因此可忽略葉片厚度上面的溫度梯度。假定葉片有一個(gè)統(tǒng)一的溫度,葉片溫度Tblade。因此,葉片變化量和溫度有線性關(guān)系:
同轉(zhuǎn)子一樣,葉片受到離心力的作用。為簡化計(jì)算,葉片和轉(zhuǎn)子增長是非常小的,因此可以近似離心力是葉片的角速度的函數(shù)。由于離心力導(dǎo)致葉片的變化量可表示為:
圖4 ??渦輪葉片模型結(jié)構(gòu)
本文中計(jì)算數(shù)據(jù)為某型發(fā)動(dòng)機(jī)試車數(shù)據(jù),測(cè)試參數(shù)有壓氣機(jī)前溫度T2,壓力P2,壓氣機(jī)后壓力P3,渦輪后溫度T5,壓力P5。
根據(jù)變比熱方法[9]可計(jì)算出壓氣機(jī)后溫度T3,渦輪前溫度T4。由于發(fā)動(dòng)機(jī)材料特性及幾何尺寸未知,本文采用公開文獻(xiàn)[4]的發(fā)動(dòng)機(jī)數(shù)據(jù),圖5為發(fā)動(dòng)機(jī)地面起動(dòng)時(shí)主要參數(shù)變化曲線,從開始推油門桿至慢車狀態(tài)穩(wěn)定,發(fā)動(dòng)機(jī)轉(zhuǎn)速、壓氣機(jī)后溫度和渦輪前溫度增大。從油門桿起動(dòng)開始至轉(zhuǎn)速穩(wěn)定大約需85 s,渦輪前溫度達(dá)到最大值后緩慢下降達(dá)到穩(wěn)定。
圖5 ??地面起動(dòng)時(shí)發(fā)動(dòng)機(jī)參數(shù)歷程
利用上述計(jì)算方法求取了起動(dòng)過程中發(fā)動(dòng)機(jī)機(jī)匣徑向變化量、轉(zhuǎn)子徑向變化量、渦輪葉片變化量,代入公式(1)中得到發(fā)動(dòng)機(jī)渦輪葉尖間隙總變化量,結(jié)果如圖6所示。
圖6 ??地面起動(dòng)發(fā)動(dòng)機(jī)渦輪葉尖間隙變化量計(jì)算結(jié)果
結(jié)合圖5和圖6可以看出,由于離心力與熱應(yīng)力的作用,發(fā)動(dòng)機(jī)轉(zhuǎn)子的變化量和葉片變化量隨著轉(zhuǎn)速和渦輪進(jìn)口溫度增加較快,由于葉片體積相對(duì)轉(zhuǎn)子較小,葉片的形變量小于轉(zhuǎn)子。機(jī)匣由于熱應(yīng)力的作用隨著溫度的增加而形變?cè)黾樱覚C(jī)匣的形變速率快于轉(zhuǎn)子和葉片的形變速率。在發(fā)動(dòng)機(jī)地面起動(dòng)過程中,渦輪葉尖間隙呈現(xiàn)減小趨勢(shì),在渦輪后溫度T5達(dá)最大值時(shí),渦輪葉尖間隙有最小值,隨后緩慢上升趨于穩(wěn)定,整個(gè)過程渦輪葉尖間隙變化率約為25.9%。
(1)發(fā)動(dòng)機(jī)地面起動(dòng)過程渦輪葉尖間隙呈減小趨勢(shì),機(jī)匣形變量上升速率最大,葉片變化量對(duì)葉尖間隙變化影響最小。
(2)本文建立了發(fā)動(dòng)機(jī)葉尖間隙模型,并用試車數(shù)據(jù)計(jì)算了發(fā)動(dòng)機(jī)地面起動(dòng)渦輪葉尖間隙變化規(guī)律,計(jì)算結(jié)果符合預(yù)期趨勢(shì)。