單級(jí)離心壓氣機(jī)葉片表面非定常壓力分布研究

摘 要：針對(duì)某單級(jí)離心壓氣機(jī)，開展了全通道非定常數(shù)值模擬計(jì)算，研究了離心葉輪葉片表面的非定常壓力分布。結(jié)果顯示，離心葉輪葉片表面壓力面上的非定常脈動(dòng)特征比吸力面?zhèn)让黠@，且葉片根部壓力脈動(dòng)最大;整體而言，葉輪葉片表面壓力分布呈現(xiàn)一定的規(guī)律，該規(guī)律對(duì)于離心壓氣機(jī)的振動(dòng)預(yù)測(cè)和改善氣動(dòng)性能及使用壽命具有較好的指導(dǎo)意義。

關(guān)鍵詞：離心壓氣機(jī);數(shù)值模擬;非定常;葉片壓力分布

離心壓氣機(jī)因其單級(jí)壓比高、結(jié)構(gòu)簡(jiǎn)單的優(yōu)點(diǎn)，是中小流量渦軸、渦槳航空發(fā)動(dòng)機(jī)和輔助動(dòng)力裝置最重要的壓縮系統(tǒng)部件。目前，在我國(guó)中小流量航空發(fā)動(dòng)機(jī)研制過程中，出現(xiàn)了不少離心壓氣機(jī)葉片裂紋和掉塊故障，而根據(jù)故障結(jié)果分析，大部分是由于高周疲勞導(dǎo)致的。引起葉輪疲勞破壞一方面是由于材料本身缺陷和加工制造原因，另一方面則是由于葉片交變載荷的長(zhǎng)期作用。針對(duì)離心壓氣機(jī)的葉片疲勞故障，較多的研究主要集中在加工制造和材料缺陷等方面。而對(duì)于離心壓氣機(jī)內(nèi)部非定常流動(dòng)影響葉輪振動(dòng)、強(qiáng)度的研究，也主要是采用理論計(jì)算和定性分析。Ju 等對(duì)一離心式壓縮機(jī)進(jìn)行了非定常流場(chǎng)計(jì)算，考察不同工況下的葉片表面壓力脈動(dòng)、應(yīng)力變化，并將結(jié)果用于壽命評(píng)估[1] 。毛義軍等計(jì)算了壓縮機(jī)首級(jí)閉式葉輪在非定常氣動(dòng)力作用下的葉片等效應(yīng)力，對(duì)疲勞斷裂進(jìn)行了分析，發(fā)現(xiàn)交變應(yīng)力集中位置和葉片斷裂位置的一致性[2]。

為實(shí)現(xiàn)壓氣機(jī)轉(zhuǎn)子葉片振動(dòng)特性、尾流激振響應(yīng)等多方面的分析預(yù)估，需將葉輪葉片在氣動(dòng)載荷激勵(lì)下的振動(dòng)應(yīng)力水平及分布進(jìn)行綜合研究。若需預(yù)測(cè)完整的振動(dòng)響應(yīng)，必須獲取真實(shí)的非定常氣動(dòng)載荷。在真實(shí)的進(jìn)出口環(huán)境下，研究離心葉輪葉片表面壓力脈動(dòng)特征，探索葉片所受交變載荷的脈動(dòng)規(guī)律，支撐離心壓氣機(jī)的氣動(dòng)改進(jìn)和振動(dòng)預(yù)測(cè)，對(duì)于提高離心壓氣機(jī)的性能、可靠性和使用壽命都具有非常重要的現(xiàn)實(shí)意義。

1 研究對(duì)象和計(jì)算方法

研究對(duì)象為單級(jí)離心壓氣機(jī)，包括徑向進(jìn)氣的進(jìn)口支板、離心葉輪、徑向擴(kuò)壓器和回流器共4排葉片，計(jì)算域如圖1示意圖所示。

數(shù)值計(jì)算采用NUMECA軟件包進(jìn)行，基于全通道的離心壓氣機(jī)網(wǎng)格模型，首先開展的定常計(jì)算[3-5]，再以此狀態(tài)點(diǎn)的定常計(jì)算收斂結(jié)果作為初場(chǎng)，保持壓氣機(jī)邊界條件不變，進(jìn)行非定常計(jì)算。

非定常計(jì)算中一個(gè)計(jì)算周期內(nèi)設(shè)置了180個(gè)時(shí)間步長(zhǎng)，即轉(zhuǎn)子葉柵每旋轉(zhuǎn)2°為一個(gè)物理時(shí)間步，時(shí)間步長(zhǎng)為8.889μs。相當(dāng)于非定常計(jì)算的采樣頻率為112.5kHz，足以滿足捕捉非定常流場(chǎng)特征的要求。同時(shí)，每個(gè)物理時(shí)間步長(zhǎng)內(nèi)取50步虛擬迭代。若將全部結(jié)果保存下來，對(duì)計(jì)算時(shí)間和存儲(chǔ)量要求巨大，故設(shè)置每20個(gè)時(shí)間步長(zhǎng)保存一個(gè)詳細(xì)流場(chǎng)計(jì)算結(jié)果，一個(gè)周期共保存9個(gè)時(shí)刻結(jié)果以供詳細(xì)分析，各個(gè)時(shí)刻分別命名為T1、T2、T3、T4、T5、T6、T7、T8和T9。

2 結(jié)果分析

圖3分別給出了各個(gè)時(shí)刻下葉片在50%葉高處的壓力分布對(duì)比圖。從圖中可以看出，離心葉輪出口附近表現(xiàn)出了非常明顯的非定常特性，即葉片近出口表面壓力曲線隨時(shí)間的變化劇烈;葉輪流道內(nèi)部的出口位置附近相對(duì)進(jìn)口表現(xiàn)出更明顯的非定常特性，而壓力面?zhèn)热~片表面壓力隨時(shí)間的波動(dòng)幅度較吸力面?zhèn)雀蟆ｋS著相對(duì)弦長(zhǎng)的增大，葉片表面的壓力脈動(dòng)也越來越大，葉片上的壓力脈動(dòng)主要存在于壓力面?zhèn)龋壹邪l(fā)生在80%相對(duì)弦長(zhǎng)位置到葉輪出口之間。主葉片表面壓力最大值出現(xiàn)在T2時(shí)刻，而分流葉片表面壓力最大值出現(xiàn)較主葉片晚兩個(gè)時(shí)刻，即出現(xiàn)在T4時(shí)刻，這可能由于葉輪出口的低速團(tuán)在周向輸運(yùn)產(chǎn)生時(shí)間滯后有關(guān)。

分析產(chǎn)生上述現(xiàn)象的原因可能是：離心葉輪出口尾跡和葉輪通道內(nèi)主流區(qū)內(nèi)的有勢(shì)流體及間隙泄漏流在遇到擴(kuò)壓器葉片前緣時(shí)產(chǎn)生干涉，且這種干涉是隨著葉輪與擴(kuò)壓器周向位置的不同而變化的，并且在葉輪出口與擴(kuò)壓器進(jìn)口間的無葉段內(nèi)氣流的摻混最為劇烈，無葉段內(nèi)的擾動(dòng)向上游傳播，從而引起了葉輪出口的壓力脈動(dòng)現(xiàn)象。

圖4給出了主葉片、分流葉片在壓力面?zhèn)炔煌~高處出現(xiàn)的壓力最大值和最大壓力脈動(dòng)值。隨著葉高的增加，葉片表面出現(xiàn)的壓力最大值和壓力脈動(dòng)最大值均呈均勻減小的趨勢(shì)，且分流葉片減小的梯度大于主葉片，說明分流葉片表面的壓力變化沿葉高更為劇烈。整個(gè)葉片表面出現(xiàn)最大壓力值和壓力脈動(dòng)最大值的地方均位于葉片的10%葉高位置，而且分流葉片根部的壓力變化較主葉片更明顯，說明葉片根部受到的氣動(dòng)激振力較大。

出現(xiàn)壓力最大值及壓力脈動(dòng)最大值對(duì)應(yīng)的相對(duì)弦長(zhǎng)位置如圖5所示?？梢?，隨著葉高增高，主葉片出現(xiàn)壓力最大值和壓力脈動(dòng)最大值的位置往葉輪出口方向緩慢后移，而分流葉片則往前緣方向均勻前移，且主葉片出現(xiàn)壓力最大值和壓力脈動(dòng)最大值的位置在不同葉高均比分流葉片靠后。

圖6給出了各個(gè)時(shí)刻下各葉片50%葉高處的葉片載荷分布對(duì)比圖，葉片載荷定義為相同弦長(zhǎng)位置點(diǎn)，壓力面靜壓與吸力面靜壓之差。從圖中可見，葉片的載荷加載量隨著弦長(zhǎng)的增大，其非定常特性越顯著。同時(shí)，由于分流葉片前緣位置相對(duì)靠后，其載荷脈動(dòng)的相對(duì)弦長(zhǎng)范圍比主葉片更寬廣。

圖7給出了主葉片和分流葉片不同葉高的相對(duì)載荷在一個(gè)周期內(nèi)不同時(shí)刻的變化示意圖，相對(duì)載荷定義為該葉高位置葉片載荷在整個(gè)相對(duì)弦長(zhǎng)范圍內(nèi)的積分值與最大值的比值，詳見公式（1）、（2）。

其中，上標(biāo)L表示弦長(zhǎng)長(zhǎng)度，Δp表示葉片表面壓力面與吸力面壓力差，ds為單位弦長(zhǎng)，為某葉高位置葉片載荷量，為的最大值。

從圖中可見，各葉片不同葉高的載荷量不同，同時(shí)也隨著時(shí)刻的變化呈現(xiàn)一定的規(guī)律變化。從葉高的變化來看，分流葉片載荷加載量隨葉高增大均勻增加，而主葉片則隨著葉高的升高加載量增長(zhǎng)的趨勢(shì)則不同，葉中和葉尖的載荷量基本一致。從時(shí)間的變化來看，同一個(gè)葉片不同葉高上的載荷加載量隨著時(shí)刻變化的趨勢(shì)基本一致;葉中和葉尖的載荷量隨時(shí)間的變化最大幅值均在20%左右，而葉根的載荷量隨時(shí)間變化幅度較大，接近30%。

3 結(jié)論

本文對(duì)某單級(jí)離心壓氣機(jī)的主葉片和分流葉片表面壓力分布進(jìn)行詳細(xì)研究，分析了不同葉高表面壓力脈動(dòng)特征和載荷變化規(guī)律，得到以下結(jié)論：

1）離心葉輪葉片表面壓力面上的非定常脈動(dòng)特征比吸力面?zhèn)让黠@，靠近葉輪出口后半弦長(zhǎng)范圍內(nèi)的壓力波動(dòng)幅度相對(duì)前半葉片弦長(zhǎng)范圍內(nèi)更大。

2）葉片表面壓力最大值和脈動(dòng)最大值均出現(xiàn)在10%葉高的葉根位置，葉片各葉高的載荷量隨時(shí)間變化的規(guī)律基本一致，但葉根位置隨時(shí)間變化的幅度更大。

參考文獻(xiàn)：

[1] Ju Y， Liu H， Yao Z， Xing P， Zhang C. Fluid-structure Interaction Analysis and Lifetime Estimation of a Natural Gas Pipeline Centrifugal Compressor under Near-choke and Near-surge Conditions[J]. Chinese Journal of Mechanical Engineering， 2015， 28（6）： 1261-1268.

[2] 毛義軍， 祁大同， 許慶余. 離心壓縮機(jī)葉輪葉片疲勞斷裂故障的數(shù)值分析[J]. 西安交通大學(xué)學(xué)報(bào)， 2008， 42（11）：1336-1339.

[3] 劉濤，李杜，賀丹，吳俊峰.湍流模型和差分格式在小流量離心壓氣機(jī)數(shù)值計(jì)算的研究[C]. 中國(guó)航空動(dòng)力分會(huì)第十一屆小型發(fā)動(dòng)機(jī)學(xué)術(shù)研討會(huì)，2017.

[4] 賀丹，李杜，劉濤，張錦綸.某單級(jí)離心壓氣機(jī)的網(wǎng)格無關(guān)性研究[C].中國(guó)航空動(dòng)力分會(huì)第十一屆小型發(fā)動(dòng)機(jī)學(xué)術(shù)研討會(huì)，2017.

[5] 吳俊峰，李杜，劉濤，賀丹.單級(jí)離心壓氣機(jī)三維全通道數(shù)值模擬計(jì)算技術(shù)研究[C].中國(guó)航空動(dòng)力分會(huì)第十一屆小型發(fā)動(dòng)機(jī)學(xué)術(shù)研討會(huì)，2017.

作者簡(jiǎn)介：

賀丹（1988-），女，土家族，湖南株洲人，碩士，工程師，現(xiàn)就職于中國(guó)航發(fā)湖南動(dòng)力機(jī)械研究所，研究方向：壓氣機(jī)氣動(dòng)設(shè)計(jì)。