賴生智, 吳亞東, 田 杰, 歐陽華
(上海交通大學(xué) a.航空航天學(xué)院;b.機(jī)械與動力工程學(xué)院;c.燃?xì)廨啓C(jī)與民用航空發(fā)動機(jī)教育部工程研究中心, 上海 200240)
由于葉頂間隙的存在,葉頂泄漏流普遍存在于壓氣機(jī)中,并引發(fā)非定常壓力脈動,導(dǎo)致葉片受非定常力而產(chǎn)生一系列振動問題[1-2],流固耦合又加劇了這種流動不穩(wěn)定性.旋轉(zhuǎn)不穩(wěn)定性(RI)是一種典型的葉頂區(qū)域非定常流動,靠近葉頂區(qū)域壁面壓力頻譜上,在小于葉片通過頻率(BPF)的某個頻率范圍內(nèi)出現(xiàn)連續(xù)的等間隔尖峰組成的寬頻帶凸起.大量研究表明,RI現(xiàn)象通常發(fā)生在小流量、大間隙工況下,隨之產(chǎn)生的葉片非同步振動[3-4]及葉頂間隙噪聲[5]都對壓氣機(jī)的工作性能及安全性產(chǎn)生較大的危害.
RI現(xiàn)象最早是由Baumgartner等[1]發(fā)現(xiàn)的,他們認(rèn)為RI現(xiàn)象會導(dǎo)致葉片非同步振動.之后,Kameier等[6]在對一臺軸流風(fēng)扇進(jìn)行噪音測試時發(fā)現(xiàn)在大葉頂間隙情況下,聲壓頻譜上小于BPF的范圍內(nèi)出現(xiàn)了寬頻凸起,葉頂區(qū)域的壓力譜也出現(xiàn)了相同的特征.隨后,Kameier等[7]對RI現(xiàn)象中的一些頻譜特性做了實驗研究及理論推導(dǎo),并得到了RI現(xiàn)象的周向模態(tài)特性.Mailach等[8]在軸流壓氣機(jī)實驗中發(fā)現(xiàn)當(dāng)發(fā)生RI現(xiàn)象時,葉頂泄漏渦的強(qiáng)度是高度脈動的,同時具有周向傳播的特征,并認(rèn)為葉頂泄漏渦的非定常性是造成RI現(xiàn)象的主要原因.Marz等[9]發(fā)現(xiàn)在壓氣機(jī)中葉頂泄漏流與來流在葉片前緣附近發(fā)生干涉,形成了一個從同一葉片通道內(nèi)的吸力面向壓力面移動的渦結(jié)構(gòu),并認(rèn)為這是RI現(xiàn)象產(chǎn)生的主要原因.Pardowitz等[10-12]則在對環(huán)形葉柵的數(shù)值與實驗研究中對周向模態(tài)隨時間的變化等進(jìn)行了探索,并提出了新的觀點,認(rèn)為RI的成因是葉頂區(qū)域剪切層的不穩(wěn)定性.陳智洋等[13]用多通道數(shù)值模擬發(fā)現(xiàn)通道中泄漏流自誘導(dǎo)的非定常性及周向傳播導(dǎo)致了RI現(xiàn)象的產(chǎn)生.王昊等[14-16]通過對壓力頻譜的分析發(fā)現(xiàn)脈動葉頂泄漏流產(chǎn)生的旋轉(zhuǎn)壓力波與旋轉(zhuǎn)不穩(wěn)定性具有一致的頻率特性.
上述觀點表明,葉頂間隙流動與RI現(xiàn)象關(guān)系緊密,不同間隙下的流動特性和頻譜特征等還需細(xì)致的闡述.本文在前期試驗捕捉到RI現(xiàn)象特征的基礎(chǔ)上,進(jìn)一步采用數(shù)值計算方法分析不同間隙情況下RI的流動特征和周向模態(tài)等特性,討論與RI相關(guān)的葉頂區(qū)域流動結(jié)構(gòu)、湍流特性和通道渦的產(chǎn)生及傳播特性.周向模態(tài)和相關(guān)性分析解釋了RI模態(tài)及頻率分布特點.對比不同間隙下葉頂區(qū)域流動及頻譜特征,并得到了RI現(xiàn)象隨葉頂間隙減小而變化規(guī)律.
研究對象為上海交通大學(xué)的一臺單級低速軸流壓氣機(jī),結(jié)構(gòu)如圖1所示,其動葉葉型為NACA65,葉片數(shù)為21個,設(shè)計轉(zhuǎn)速為 3 000 r/min,葉頂間隙與葉高的比值為2.2%,其他裝置細(xì)節(jié)描述詳見文獻(xiàn)[14-16].其氣動性能曲線和本文選取的計算工況點及某流量系數(shù)下機(jī)匣壁面的壓力功率譜(PSD)如圖2所示,圖中φ為流量系數(shù),f為頻率,k為總壓系數(shù),σ為功率譜密度.可以看到,在該間隙下,功率譜上BPF以下存在明顯的RI頻率特征,表現(xiàn)為多個寬頻帶凸起,每個凸起由一系列等頻率間隔的尖峰組成.
數(shù)值計算使用商業(yè)軟件ANSYS CFX,采用的湍流模型為模型k-ε.由于動靜干涉效應(yīng)的存在,導(dǎo)葉的尾跡對動葉葉頂區(qū)域的流場脈動存在干擾,影響RI特性的識別.本文中僅對動葉區(qū)域進(jìn)行計算,此簡化在文獻(xiàn)[14]中已有討論.同時,為了捕捉到非葉片周期的周向非定常流動特性,采用全通道模型.計算選用的網(wǎng)格如圖3所示,其無關(guān)性驗證后的網(wǎng)格總數(shù)約為1.29×107.邊界條件為入口給定總溫總壓,出口給定靜壓.時間分辨率為 0.000 1 s.計算穩(wěn)定后,取 4 096 個時間步的流場結(jié)果,對應(yīng)旋轉(zhuǎn)周期超過20個,用于后文的模態(tài)分解,模態(tài)分解的頻率分辨率約為2.44 Hz.
圖1 試驗臺結(jié)構(gòu)圖
圖2 壓氣機(jī)試驗數(shù)據(jù)
圖3 計算網(wǎng)格細(xì)節(jié)
圖4 監(jiān)測點布置
計算中在周向和軸向布置一系列監(jiān)測點以獲取葉頂區(qū)域在絕對坐標(biāo)系和與轉(zhuǎn)子相對靜止的相對坐標(biāo)系下的壓力時域信號.監(jiān)測點布置方式如圖4所示.在軸向方向設(shè)置9個截面,分布在從葉頂前緣前方到尾緣后方,記為Z1~Z9.圖中橫向為葉頂區(qū)域一周的展開方向.在每個截面內(nèi)近機(jī)匣壁面處周向均布120個監(jiān)測點.
(1)
(2)
式中:N為葉片通道數(shù);p為各個葉片通道的相同周向位置點的壓力;θ為監(jiān)測點所在方位角;α為一個葉片通道所占周向角度.pRMS表征了流場的脈動強(qiáng)度.
對空間位置為θj測點的壓力時域信號p(θj,t)進(jìn)行傅里葉變換,可以得到壓力頻域分布p(θj,ω),進(jìn)一步可以得到模態(tài)m的幅值的復(fù)數(shù)值A(chǔ)(m,ω),主要計算公式為
(3)
(4)
(5)
式中:ω為監(jiān)測點壓力頻率.根據(jù)恩奎斯特采樣定律,模態(tài)分解得到的m的范圍為[-59,60].本文對4種葉頂間隙與葉高比值即1.98%、1.69%、1.35%和0.84%時RI現(xiàn)象的特性進(jìn)行分析,為了討論方便,這4種間隙條件簡寫為1.98%TC、1.69%TC、1.35%TC以及0.84%TC.
圖5 不同時刻機(jī)匣壁面靜壓云圖
圖6 流道截面壓力云圖
2.1.1葉頂區(qū)域流動結(jié)構(gòu) 間隙為1.98%TC時,相對坐標(biāo)系下的機(jī)匣壁面靜壓分布如圖5所示,圖中展示了12個葉片范圍內(nèi)機(jī)匣壁面靜壓(相對壓力,pr)隨時間(t)的變化.數(shù)字1~12代表葉片所在位置,每張橫圖表示轉(zhuǎn)子中12個葉片的葉頂間隙區(qū)域的展開,其軸向范圍為葉片向前后各延展約一倍軸向弦長長度(c).縱向13張橫圖對應(yīng)連續(xù)的13個時間步.在t時刻,以1號葉片下游虛線橢圓標(biāo)出的兩處較小的高壓區(qū)為特征,可以識別到另一周期數(shù)的周向流動結(jié)構(gòu).根據(jù)橫向規(guī)律,此特征在7、12以及圖中未展示的17號葉片的下游再次出現(xiàn),即可以看作是1個周期數(shù)為4的周向流動結(jié)構(gòu).4個區(qū)域所包含的葉片數(shù)不一致的原因為21個轉(zhuǎn)子葉片引起的基本流動很強(qiáng),從而影響了較弱流動的周向可分辨性.沿圖中箭頭方向,1號葉片下游的此特征在13t時刻再次出現(xiàn),即此周向流動結(jié)構(gòu)對同一葉片的擾動頻率約為1/(12t)=833 Hz.由此可以看出,由于葉頂間隙的存在,周向上出現(xiàn)了葉片周期以外的流動結(jié)構(gòu).
t時刻的葉片尾緣截面靜壓和葉片前緣截面總壓分布如圖6所示.圖6(a)中,在葉片尾緣能清晰地看到該周向流動的周期性.以圖中靠近機(jī)匣處的2個藍(lán)色低壓區(qū)和2個紅色高壓區(qū)為標(biāo)志,21個葉片通道可由圖中黑色虛線劃分為4個區(qū)域,即1個6通道區(qū)和3個5通道區(qū),表征1個周期數(shù)為4的周向流動結(jié)構(gòu).圖6(b)中,在葉片前緣葉尖高度有25個高壓區(qū),表征1個周期數(shù)為25的周向流動結(jié)構(gòu),是除葉片數(shù)周期流動之外最明顯的一個結(jié)構(gòu).機(jī)匣壁面處的對應(yīng)位置存在25個低壓區(qū),這些低壓區(qū)的強(qiáng)度在周向上也呈現(xiàn)出明顯的周期性.圖中1、2,7、8,13、14,以及20、21這4組區(qū)域機(jī)匣壁面處出現(xiàn)了幅值相近的低壓區(qū),同樣如圖中黑色虛線所劃,3個6低壓區(qū)和1個7低壓區(qū)表征的周期數(shù)為4的周向流動結(jié)構(gòu),與圖6(a)中觀察到結(jié)果一致.
綜上所述,葉頂區(qū)域存在一個周期數(shù)為25的周向流動,與葉片發(fā)生干涉,產(chǎn)生了1個強(qiáng)度相近、周期數(shù)為4的周向流動結(jié)構(gòu),該結(jié)構(gòu)在葉片前緣和尾緣均清晰可見,而周期數(shù)為25的周向流動結(jié)構(gòu)僅在前緣可見,說明這種流動結(jié)構(gòu)起始于葉片前緣,在達(dá)到尾緣前已幾乎消失,且在周向和軸向都具有傳播特性.
圖7 機(jī)匣壁面和葉片表面壓力云圖
圖8 葉頂區(qū)域流線與通道渦結(jié)構(gòu)
2.1.3通道渦形成和傳播 經(jīng)過處理得到葉頂區(qū)域流線與通道渦結(jié)構(gòu),如圖8所示,圖中p′為相對總壓.圖8(a)顯示了前緣渦和葉頂泄漏渦合成一個通道渦的過程.前緣渦起始于葉片前緣并在周向傳播的過程中,葉頂泄漏流在葉片中部從壓力面向吸力面運動,加入到前緣渦中并一起傳播.兩者合成了通道渦,跨過葉片通道到后一葉片的前緣附近,如此周向傳播.
圖8(b)和(c)所示的葉頂區(qū)域流線圖更清楚地展示了通道渦合成的細(xì)節(jié).圖8(b)中,前一葉片的葉頂泄漏流線翻過葉頂區(qū)域形成通道渦,遇到下一葉片的壓力面,部分從通道中流出,同時形成了下一葉片的泄漏流動,但實際上應(yīng)將前一葉片的葉頂泄漏流分為兩部分看待.圖8(c)中,藍(lán)色為葉片前緣處流線,綠色為泄漏流線的前部分,紅色為泄漏流線的后部分.藍(lán)色流線和綠色流線在葉片吸力面附近發(fā)生干涉,合成了通道渦,隨后部分流出通道,部分成為下一葉片的泄漏流;而紅色流線基本未參與這個渦的形成,而是跨過葉片通道,組成了下一葉片的前緣流動,故其也是前緣渦的一個重要來源.綜上,通道渦的形成是前緣渦和葉頂泄漏渦的前部分與軸向主流干涉形成的.
圖9所示為通道渦的傳播過程中13個連續(xù)時間步的計算結(jié)果.圖中可以看到通道渦像一道波浪,打在下一葉片的前緣,破碎成兩部分,一部分為下一葉片的前緣渦,一部分向壓力面?zhèn)鞑?同時可以看到,t時刻和13t時刻的通道渦分布幾乎相同,渦結(jié)構(gòu)的發(fā)展過程清晰可辨,可以認(rèn)為通道渦的波動周期是12t,即頻率為833 Hz,與上文的結(jié)果一致.分析表明,由前緣渦和葉頂泄漏渦合成的通道渦是周期數(shù)為25且擾動頻率約為833 Hz的周向流動結(jié)構(gòu).
圖9 通道渦的發(fā)展過程
對相對坐標(biāo)系下相鄰葉片通道內(nèi)同角度位置點P1和P2的壓力信號相干性和相位關(guān)系作以下分析.以相干函數(shù)表示兩組信號的相干性:
(6)
式中:G(f)表示功率譜密度函數(shù)的復(fù)數(shù)值,下標(biāo)數(shù)字相同表示一組信號的自功率譜,下標(biāo)數(shù)字不同表示兩組信號的互功率譜.相干函數(shù)值在0~1之間,函數(shù)值越大,表示相干性越高.
圖10 P1和P2信號的相關(guān)分析
兩組信號的互功率譜(C)及相干程度如圖10所示,圖中δ為相位.在RI所在頻率RIF與2RIF發(fā)生的頻率內(nèi)兩組信號的相干函數(shù)接近于1,相干程度明顯高于其他頻率范圍.同時,兩組信號的互功率譜相位在此頻率范圍內(nèi)呈現(xiàn)出良好的線性關(guān)系.這證明了造成RI現(xiàn)象的流動在周向上具有傳播特性,同時與葉片數(shù)也有很大關(guān)系.
圖中線性關(guān)系的斜率可以作如下解釋.相對坐標(biāo)系下,模態(tài)m旋轉(zhuǎn)一周,監(jiān)測點壓力信號相位總共變化2πm,所以P1、P2點的相位差為
(7)
式中:B為葉片數(shù);i=1,2,…,即考慮到圓周的周期性,互功率譜相位中的相位差限制在-2π到2π內(nèi),在圖10中可以得到驗證,例如832.51 Hz處的 -25 模態(tài),依照上式計算值為 -7.5,取i=1得到相位差為 -1.2(π取3.14),與圖中讀到的數(shù)值相符.
考慮相鄰模態(tài)的相位差值與互功率譜斜率的關(guān)系,即
(8)
葉片前緣即Z2面上120個監(jiān)測點在絕對和相對坐標(biāo)系下的周向模態(tài)分解結(jié)果如圖11所示.圖中S為聲壓級,數(shù)字為模態(tài)階數(shù).圖中,每個頻率尖峰都對應(yīng)于一個周向模態(tài).其中幅值最大的±25模態(tài)是RI現(xiàn)象的主模態(tài),相對坐標(biāo)系中頻率為832.51 Hz,與2.1節(jié)中流場分析發(fā)現(xiàn)的周期數(shù)為25、擾動頻率為833 Hz的通道渦結(jié)構(gòu)一致,說明RI現(xiàn)象是由通道渦的周向傳播引起的.
圖11 Z2面上監(jiān)測點周向模態(tài)分解
在絕對坐標(biāo)系中,可以清晰地看到在BPF以下存在4個頻率凸起帶,每個凸起帶都由一系列頻率間隔相近的尖峰組成,頻率間隔約為21.0 Hz.RI現(xiàn)象對應(yīng)于圖中的第2個凸起帶,其發(fā)生的頻率范圍約為340~480 Hz,對應(yīng)的模態(tài)階數(shù)為21~28;而第3個凸起帶在各個尖峰的頻率上恰好為BPF與RIF之差,其范圍約為570~710 Hz,尖峰的模態(tài)階數(shù)也為BPF與RIF對應(yīng)模態(tài)之差,即 -7~0.另外兩個凸起帶對應(yīng)于2RIF和BPF-2RIF,模態(tài)階數(shù)和頻率的和差關(guān)系也同樣存在.各凸起帶的頻率、模態(tài)階數(shù)的和差關(guān)系說明了造成RI現(xiàn)象的通道渦與葉片之間存在干涉現(xiàn)象.
而在相對坐標(biāo)系中,RIF連續(xù)峰值的頻率間隔約為 -29.0 Hz.RIF對應(yīng)于連續(xù)模態(tài) -15~-31,其中-21~-28模態(tài)與絕對坐標(biāo)系下標(biāo)出的21~28模態(tài)一一對應(yīng),負(fù)號代表在兩個參考系下,模態(tài)波的旋轉(zhuǎn)方向相反.同一個模態(tài)在兩個坐標(biāo)系下對應(yīng)的頻率不同,所有模態(tài)組合起來,產(chǎn)生了RI頻率范圍偏移現(xiàn)象,且兩種坐標(biāo)系下連續(xù)模態(tài)的角頻率間隔之差(Δω)恰好為轉(zhuǎn)子旋轉(zhuǎn)角速度(Ω),即
ΔωF-ΔωR=Ω
(9)
式中:上標(biāo)F和R分別代表絕對坐標(biāo)系和相對坐標(biāo)系.其頻率間隔符合ΔfF-ΔfR=50 Hz.
對應(yīng)RI的模態(tài)m在頻譜上的頻率與旋轉(zhuǎn)角速度有關(guān),在不考慮RI自身振動頻率的情況下,其在絕對坐標(biāo)系和相對坐標(biāo)系下的旋轉(zhuǎn)角速度分別為
(10)
(11)
由式(10)與(11)計算得到圖11中RI在絕對和相對坐標(biāo)系下不同模態(tài)下的旋轉(zhuǎn)角速度,并列出模態(tài)和頻率的關(guān)系,如圖12所示.圖中,絕對和相對坐標(biāo)系下頻率與模態(tài)m呈線性關(guān)系,線形擬合得到斜率分別為21.0和 -29.0 Hz,兩坐標(biāo)系下擬合結(jié)果的縱軸截距均為 -106.8 Hz.但是不同模態(tài)的旋轉(zhuǎn)角速度并不相同,隨模態(tài)數(shù)增加呈現(xiàn)遞增或遞減趨勢.根據(jù)2.2節(jié)的流動分析,RI現(xiàn)象是由周向傳播的通道渦導(dǎo)致的,其旋轉(zhuǎn)角速度在RIF上各尖峰模態(tài)應(yīng)當(dāng)為一定值.
根據(jù)Kameier[7]提出的旋轉(zhuǎn)源機(jī)理假設(shè),旋轉(zhuǎn)角速度不同的現(xiàn)象是對應(yīng)RI的流動結(jié)構(gòu)自身振動的結(jié)果.Kameier在絕對坐標(biāo)系和相對坐標(biāo)系之外,在RI流動結(jié)構(gòu)(本文中為周向傳播的通道渦)上建立旋轉(zhuǎn)源坐標(biāo)系Q,并得到以下關(guān)系式:
圖12 f,Ω與m的線性擬合
式(14)顯示了RI現(xiàn)象在絕對坐標(biāo)系和相對坐標(biāo)系下同一模態(tài)下頻率的偏移,兩者絕對數(shù)值的加和等于轉(zhuǎn)子頻率與模態(tài)數(shù)的乘積.圖11中所顯示的RI在絕對和相對坐標(biāo)系下對應(yīng)相同模態(tài)數(shù)下的頻率關(guān)系驗證了該式.
以上結(jié)果從機(jī)理上指出了RI現(xiàn)象的脈動特征,并解釋了頻譜中呈現(xiàn)出等頻率間隔且模態(tài)數(shù)連續(xù)尖峰的現(xiàn)象以及其在絕對和相對坐標(biāo)系下的頻率之間的關(guān)系.同時,表明了RI自身擾動頻率的存在及其在絕對和相對坐標(biāo)系下的旋轉(zhuǎn)角速度與頻率間隔的關(guān)系.
本節(jié)針對不同間隙情況下與RI相關(guān)的流動分布進(jìn)行討論.圖13為0.84%TC、1.35%TC、1.69%TC和1.98%TC中的葉片前緣流道截面總壓云圖,1.98%TC的結(jié)果已在前文給出,為了對比方便,此處再次列出,下同.1.69%TC和1.35%TC中,葉頂區(qū)域可以分別識別出周期數(shù)為26和30的周向流動結(jié)構(gòu),與這兩種間隙條件下的RI主模態(tài)階數(shù)一致,而0.84%TC中則看不到周期流動結(jié)構(gòu).圖中僅定性地展示壓力區(qū)的周期性,因此未給出圖例.
圖14是不同間隙條件下葉頂區(qū)域的渦量和總壓分布.隨著葉頂間隙減小,葉頂區(qū)域流動明顯減弱.尤其最小間隙0.84%TC中,葉頂流動強(qiáng)度小且無明顯周期性變化,未觀察到明顯的渦合成以及通道渦在前緣破碎的過程.
圖15是0.84%TC下Z2截面的壓力信號頻譜.可以看出,在絕對坐標(biāo)系下,小于BPF的范圍內(nèi)未發(fā)現(xiàn)明顯RI凸起帶特征,取而代之的是全頻率范圍內(nèi)的連續(xù)尖峰,相鄰間隔為50 Hz,為明顯的轉(zhuǎn)子頻率的特征,可以認(rèn)為在這種間隙下RI現(xiàn)象幾乎消失,壓氣機(jī)中最顯著的流動是轉(zhuǎn)子旋轉(zhuǎn)引起的.然而,在500 Hz左右,即大間隙下RI出現(xiàn)的位置,尖峰的幅值仍稍高于其他尖峰,表明此時仍有微弱的RI現(xiàn)象存在,這一點在相對坐標(biāo)系下更為清晰,在原RI出現(xiàn)的頻率段仍有小幅的頻率凸起帶,但其高階倍頻則不可見.
圖13 葉片前緣流道截面總壓云圖
圖14 不同間隙條件葉頂流動
圖15 0.84%TC中Z2截面葉頂壓力頻譜
0.84%TC間隙下模態(tài)分解結(jié)果如圖16所示.圖中顯示該間隙下的周向模態(tài)特性與大間隙情況有很大的區(qū)別.絕對坐標(biāo)系下,連續(xù)尖峰對應(yīng)的模態(tài)不連續(xù),而是以3組模態(tài)交替出現(xiàn),且每組模態(tài)的間隔為3,即圖中所示的9,12,15…;-30,-27,-24…;51,54,57…,整個頻率段內(nèi)都具有此規(guī)律.同時,在旋轉(zhuǎn)坐標(biāo)系下的模態(tài)分解圖中幾乎觀察不到顯著的周向模態(tài)特征,每一個單獨模態(tài)都表現(xiàn)出近似的頻譜特性.這些特點都表明了小間隙下RI流動的消失.
交替連續(xù)的模態(tài)特性的形成可通過圖17進(jìn)行解釋.圖17(a)為絕對坐標(biāo)系下各頻率下的所有模態(tài)分布圖.可以看到,圖中有數(shù)條斜率為1/(50 Hz)的亮帶,但觀察不到各階RIF.只保留每個頻率下的最大模態(tài),去掉其他背景后得到圖17(b).可以看到頻率逐漸增大時,最大模態(tài)交替出現(xiàn)在3條亮帶中,即按照圖16(a)中的9,-30,51,12,-27,54,15,-24,57…的順序出現(xiàn).這些亮帶都是轉(zhuǎn)子旋轉(zhuǎn)的表現(xiàn),在各種間隙情況下都存在,只是當(dāng)小間隙下RI強(qiáng)度很小的時候,才在頻譜圖和模態(tài)分解圖中顯現(xiàn).這說明了在此小間隙下RI現(xiàn)象基本消失.而轉(zhuǎn)子旋轉(zhuǎn)的模態(tài)為何表現(xiàn)出3條平行的亮帶,且只觀察到階數(shù)為3的倍數(shù)的模態(tài)則需要進(jìn)一步探究.
上文討論了在0.84%TC和1.98%TC小間隙與大間隙的周向模態(tài)特性,下文分析1.35%TC的模態(tài)特性.圖18是絕對坐標(biāo)系下的模態(tài)分解圖.圖18(a)是所有模態(tài)結(jié)果,圖18(b)~18(d)分別是轉(zhuǎn)子模態(tài)、BPF-RIF及RIF對應(yīng)的模態(tài)細(xì)節(jié).圖18(a)中可以觀察到雜亂的模態(tài)分布,既有0.84%TC中間隔為3的三組模態(tài)交替連續(xù)的特征(圖18(b)),也有1.98%TC中典型的RI連續(xù)模態(tài)特征(圖18(d)),且兩種特征的幅值相近,在模態(tài)圖上都能顯現(xiàn)出來,因此可以看作是兩種特征的疊加.這說明了1.35%TC的模態(tài)特征介于大小間隙之間,RI現(xiàn)象強(qiáng)度相較于大間隙情況有所減弱,與轉(zhuǎn)子轉(zhuǎn)動模態(tài)相當(dāng).另外,圖18(c)和(d)顯示RIF與BPF-RIF的頻率范圍很接近,約為0.5BPF,也造成了模態(tài)圖的雜亂.
圖16 0.84%TC中Z2截面葉頂壓力頻譜
圖17 0.84%TC中絕對坐標(biāo)系下模態(tài)階數(shù)與頻率的關(guān)系
圖18 絕對坐標(biāo)系下模態(tài)細(xì)節(jié)
對一臺單級低速軸流壓氣機(jī)RI特性進(jìn)行分析.在前期大間隙情況下試驗和數(shù)值計算的基礎(chǔ)上,進(jìn)一步采用全通道數(shù)值計算方法分析4種葉頂間隙條件下旋轉(zhuǎn)不穩(wěn)定性的頻率特性、模態(tài)特性和流場特性的變化及原因,主要結(jié)論如下.
(1)大葉頂間隙1.98%TC情況下,RI現(xiàn)象發(fā)生時,葉片前緣存在一個周期數(shù)約為25、擾動頻率約為833 Hz的周向流動結(jié)構(gòu),與RI的主模態(tài)階數(shù)和頻率相對應(yīng).在葉頂前緣吸力面位置渦結(jié)構(gòu)產(chǎn)生的區(qū)域壓力脈動強(qiáng)烈.葉片前緣高幅值的前緣渦與葉片較低幅值葉頂泄漏渦的中前部合成一個通道渦,傳播到下一葉片前緣時破碎成下一葉片的前緣渦和葉頂泄漏渦兩部分.頻率和模態(tài)的一致性證明了該通道渦是造成RI現(xiàn)象的主要原因.
(2)相鄰葉片通道內(nèi)同角度位置壓力信號相干性分析顯示,在RIF與2RIF發(fā)生的頻率內(nèi)兩組信號的相干函數(shù)接近于1,相干程度明顯高于其他頻率范圍.同時,兩組信號的互功率譜相位在此頻率范圍內(nèi)呈現(xiàn)出良好的線性關(guān)系.證明產(chǎn)生RI流動的通道渦結(jié)構(gòu)在周向上具有傳播特性.
(3)RI的頻率特性表現(xiàn)在頻譜上小于BPF的范圍內(nèi),存在多個由一系列等頻率間隔尖峰組成的寬頻凸起帶,對應(yīng)尖峰頻率的空間模態(tài)數(shù)連續(xù)分布.采用旋轉(zhuǎn)源機(jī)理解釋了RI現(xiàn)象是由一種自身具有脈動角頻率,同時在絕對和相對坐標(biāo)系下以一定旋轉(zhuǎn)角速度轉(zhuǎn)動的通道渦流動結(jié)構(gòu)造成的.同一坐標(biāo)系下該通道渦的旋轉(zhuǎn)角速度相同,使得不同模態(tài)下RI在絕對和相對坐標(biāo)下表現(xiàn)出等頻率間隔分布.同一模態(tài)在絕對坐標(biāo)系和相對坐標(biāo)系下的頻率存在偏移,兩者絕對數(shù)值的加和等于轉(zhuǎn)子頻率與模態(tài)數(shù)的乘積.
(4)葉頂間隙較大時,出現(xiàn)典型的RI現(xiàn)象.隨間隙的減小,葉頂流動強(qiáng)度減弱,RI主模態(tài)階數(shù)增大,且RI的強(qiáng)度逐漸減小,而強(qiáng)度相對較小的轉(zhuǎn)子轉(zhuǎn)動的影響逐漸顯露,兩坐標(biāo)系下頻率間隔之和保持在約50 Hz.0.84%TC中,葉頂泄漏渦幾乎不具有脈動特性,RI現(xiàn)象幾乎消失,模態(tài)分解結(jié)果顯示出3組交替連續(xù)的轉(zhuǎn)子模態(tài)特征.1.35%TC模態(tài)結(jié)果表現(xiàn)為0.84%TC和1.98%TC的特征疊加的雜亂現(xiàn)象.證明葉頂流動中通道渦的脈動強(qiáng)度和周向傳播與RI現(xiàn)象強(qiáng)度間具有重要聯(lián)系.
致謝本研究得到了北京先進(jìn)航空發(fā)動機(jī)協(xié)同創(chuàng)新中心的支持.