黃剛 劉婕妤 王瑤 梁仍康 楊思帆

摘要： 文章以E3高壓渦輪為研究對(duì)象，研究侵入式探針對(duì)單級(jí)高壓渦輪的氣動(dòng)性能影響機(jī)理，分析在高溫高速氣體沖擊下探針自身所產(chǎn)生的形變。結(jié)果表明探針在轉(zhuǎn)子出口測(cè)試截面，會(huì)影響渦輪的總體性能，使渦輪等熵效率降低，探針尾流區(qū)域熵增，尾跡與主流產(chǎn)生摻混損失，出口氣流角發(fā)生改變。探針在安全范圍內(nèi)發(fā)生形變，且與發(fā)動(dòng)機(jī)無(wú)共振隱患。

Abstract： Taking E3 high-pressure turbine as the research object， this paper studies the influence mechanism of intrusive probe on the aerodynamic performance of single-stage high-pressure turbine， and analyzes the deformation of probe itself under the impact of high-temperature and high-speed gas. The results show that the test section of the probe at the rotor outlet will still affect the overall performance of the turbine， reduce the isentropic efficiency and flow of the turbine， increase the entropy of the probe wake area， produce mixing loss between the wake and the mainstream， and change the outlet flow angle. The probe deforms within the safe range and has no potential resonance with the engine.

關(guān)鍵詞： 侵入式探針;高壓渦輪;等熵效率;共振

Key words： invasive probe;high pressure turbine;isentropic efficiency;resonance

中圖分類號(hào)：V231.3 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文章編號(hào)：1674-957X（2022）03-0065-04

0 ?引言

渦輪是航空發(fā)動(dòng)機(jī)的重要部件，其處于高溫高壓的惡劣環(huán)境，是有著強(qiáng)壓力和溫度梯度的強(qiáng)三維性、非定常性和強(qiáng)剪切的復(fù)雜流動(dòng)[1]。為測(cè)得渦輪的性能參數(shù)，一般采用總體性能測(cè)量方法，在渦輪的進(jìn)出口布置穩(wěn)態(tài)測(cè)試梳或測(cè)試耙，獲取截面氣流平均總參數(shù)分析得渦輪性能。

試驗(yàn)時(shí)安裝侵入式探針會(huì)浸入到渦輪流場(chǎng)中，造成流動(dòng)阻塞，當(dāng)探針測(cè)試截面位于渦輪進(jìn)口或級(jí)間時(shí)，探針探頭和支桿產(chǎn)生的尾跡會(huì)進(jìn)入渦輪流場(chǎng)中而對(duì)渦輪性能產(chǎn)生影響。甚至在某三級(jí)軸流壓氣機(jī)性能試驗(yàn)時(shí)，探針支桿的尾跡脫落渦誘導(dǎo)了轉(zhuǎn)子葉片產(chǎn)生超過(guò)安全閾值的振動(dòng)應(yīng)變[2-4]。

國(guó)內(nèi)外研究人員對(duì)試驗(yàn)過(guò)程中探針支桿對(duì)壓氣機(jī)和渦輪氣動(dòng)性能的影響已有一些研究。Wyler[5]研究了在封閉通道中的圓柱形探針支桿的阻塞效應(yīng)，使得流場(chǎng)靜壓降低，馬赫數(shù)增大，為減小探針影響，應(yīng)該盡可能的減小探針阻塞比。羅鉅[6]采用了一種直徑8mm等環(huán)面分布6點(diǎn)的梳狀總壓探針來(lái)獲取兩級(jí)軸流壓氣機(jī)進(jìn)出口的總壓。朱銘敏[7]在四級(jí)低速大尺寸壓氣機(jī)進(jìn)出口截面分別布置了4支5點(diǎn)和4支7點(diǎn)總壓探針。按照航標(biāo)建議，進(jìn)口總溫總壓測(cè)量截面可在導(dǎo)葉前緣1.5倍葉中弦長(zhǎng);出口總溫和總壓探針可布置于0.5～1.0倍轉(zhuǎn)子葉中弦長(zhǎng)，在探針尺寸限制下，測(cè)點(diǎn)數(shù)應(yīng)盡量多的布置[8]。Frederick等[9]研究了探針的測(cè)試布局對(duì)流場(chǎng)測(cè)量結(jié)果的影響，探針數(shù)量少于在8梳40點(diǎn)的測(cè)試布局對(duì)導(dǎo)致對(duì)總壓的測(cè)量誤差增大。付少林[10]等人利用近壁邊界層模型改進(jìn)了渦輪測(cè)試效率的精度，認(rèn)為高精度測(cè)試效率需要7～10個(gè)徑向測(cè)點(diǎn)和5個(gè)周向測(cè)點(diǎn)。張小慶[11]采用熱流耦合模擬篩選出合適的總溫探針，并用其測(cè)量了60mm脈沖燃燒風(fēng)洞流場(chǎng)。

本文工作主要是采用數(shù)值模擬方法，研究探針支桿的物理堵塞和局部干擾，關(guān)注侵入式探針改變?cè)辛鲌?chǎng)的情況。首次研究在渦輪出口測(cè)試截面帶有實(shí)際探針支桿的測(cè)試布局，分析其對(duì)單級(jí)高壓渦輪的性能影響。

1 ?計(jì)算模型及方法

本文采用有詳細(xì)試驗(yàn)數(shù)據(jù)的普惠E3單級(jí)高壓渦輪為研究對(duì)象，有24個(gè)鈍前緣導(dǎo)葉片和54個(gè)高扭轉(zhuǎn)角轉(zhuǎn)子葉片。采用UG軟件對(duì)E3高壓渦輪葉片進(jìn)行逆向建模，得到較為準(zhǔn)確的CAD模型;探針模型采用直徑5mm的圓柱支桿，如圖1、圖2所示。

本次計(jì)算采用單通道模型驗(yàn)證CFD計(jì)算模型的精度和逆向建模的準(zhǔn)確度。計(jì)算含有探針的渦輪流場(chǎng)時(shí)采用非定常方法，4：9轉(zhuǎn)靜通道，探針支桿周向均布6個(gè)，呈“水”字形分布，軸向位置為距轉(zhuǎn)子0.5倍轉(zhuǎn)子葉中弦長(zhǎng)（記為b），即測(cè)試截面位置，考慮真實(shí)試驗(yàn)情況，探針支桿距輪轂面有2mm間距。計(jì)算域進(jìn)口至導(dǎo)葉距離為2.5倍導(dǎo)葉高度（記為h），為消除計(jì)算中的回流影響，出口延長(zhǎng)至距轉(zhuǎn)子后緣5b位置。本文采用Numeca AG5對(duì)渦輪葉片區(qū)域生成O4H型結(jié)構(gòu)網(wǎng)格，采用Icem對(duì)探針計(jì)算域劃分非結(jié)構(gòu)網(wǎng)格，探針尾流方向沿出氣角方向加密，其中第一層邊界層網(wǎng)格高度為2e-6m，保證流場(chǎng)y+均值<2，以滿足湍流模型的使用條件。單通道網(wǎng)格（不含探針）數(shù)量約100萬(wàn)，多通道網(wǎng)格（含探針）數(shù)量為600萬(wàn)。

1.1 邊界條件

數(shù)值計(jì)算方面分別采用了ANSYS CFX定常及非定常變比熱模塊，采用SST Kω湍流模型，周期性邊界，絕熱無(wú)滑移壁面，導(dǎo)熱系數(shù)采用溫度相關(guān)函數(shù)，采用Sutherland方法計(jì)算粘性系數(shù)。渦輪進(jìn)口總壓為1324491Pa，總溫1633k，出口靜壓270kPa，轉(zhuǎn)子轉(zhuǎn)速13232Rpm。

1.2 無(wú)關(guān)性驗(yàn)證

采用定常單通道無(wú)探針的計(jì)算條件，對(duì)普惠E3單級(jí)高壓渦輪進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證，以消除網(wǎng)格數(shù)量對(duì)數(shù)值計(jì)算結(jié)果的影響。圖3為網(wǎng)格量73萬(wàn)、100萬(wàn)和120萬(wàn)網(wǎng)格量的徑向分布渦輪落壓比對(duì)比，可見73萬(wàn)網(wǎng)格與后兩者有誤差，因此在保證精度的條件下，本文采用的100w單通道渦輪計(jì)算域網(wǎng)格比較合適。

2 ?計(jì)算結(jié)果及討論

2.1 流場(chǎng)計(jì)算分析

圖4為渦輪導(dǎo)葉中截面靜壓計(jì)算值與E3渦輪試驗(yàn)值對(duì)比，其壓力面計(jì)算結(jié)果符合程度好，而導(dǎo)葉吸力面受多種激波影響，在吸力面中徑至近尾緣處符合程度較差。

表1為有無(wú)探針情況下的數(shù)值仿真結(jié)果與試驗(yàn)值的對(duì)比，可見無(wú)探針計(jì)算條件下，等熵效率、落壓比等流場(chǎng)參數(shù)均與試驗(yàn)值相差不大，可見數(shù)值計(jì)算的精度很好。而在渦輪出口測(cè)試截面加入了探針后，等熵效率降低了0.15%，流量降低了0.6%，其余參數(shù)影響程度較低?？梢娫跍u輪出口的探針對(duì)上游的導(dǎo)葉和轉(zhuǎn)子的氣動(dòng)影響較小，而對(duì)總的性能參數(shù)如效率和流量有著一定的影響，若試驗(yàn)時(shí)不修正探針的影響，將導(dǎo)致渦輪試驗(yàn)的精度降低。

表1已經(jīng)顯示侵入式探針對(duì)渦輪總體性能的影響，因此下面詳細(xì)對(duì)比下流場(chǎng)參數(shù)的變化位置和影響情況。對(duì)于渦輪氣動(dòng)性能來(lái)說(shuō)，總壓、氣流角等氣動(dòng)參數(shù)及它們的分布起到關(guān)鍵作用。對(duì)于整個(gè)轉(zhuǎn)子葉片來(lái)講，流場(chǎng)參數(shù)的周向分布并不均勻，相鄰葉片間的荷載和氣動(dòng)參數(shù)均有差異，因此采用多通道和非定常計(jì)算方法能更準(zhǔn)確的捕捉實(shí)際渦輪葉片的真實(shí)流動(dòng)情況。

圖5為渦輪導(dǎo)葉中徑葉高的絕對(duì)馬赫數(shù)分布對(duì)比云圖。可見探針上游的絕對(duì)馬赫數(shù)分布沒(méi)有大的變化，而探針下游區(qū)域有明顯的低馬赫數(shù)區(qū)域，受雷洛數(shù)影響，5mm探針支桿尾跡脫渦程度較輕，但低數(shù)脫渦區(qū)影響了軸向4倍直徑以上的區(qū)域，尾跡有一定程度地?cái)U(kuò)散。

圖6為渦輪導(dǎo)葉中徑葉高的靜熵分布對(duì)比云圖，由于整個(gè)數(shù)值計(jì)算過(guò)程為絕熱流動(dòng)，渦輪進(jìn)口截面至導(dǎo)葉區(qū)域不做功，因此該區(qū)域靜熵?zé)o變化。而導(dǎo)葉吸力面、尾跡及轉(zhuǎn)子固壁、尾跡區(qū)域受膨脹波、正激波、尾緣脫渦激波、燕尾波及剪切層等波系影響，其熵增很大，即總壓損失增加。明顯可見，探針的加入，尾跡與主流摻混，導(dǎo)致探針尾流區(qū)域熵增，產(chǎn)生了較大的損失。圖7為渦輪導(dǎo)葉中徑葉高的絕對(duì)出口氣流角分布對(duì)比云圖。探針的存在使得氣流角周向分布不一致，探針的脫渦改變了原本的氣流方向角。

圖8為測(cè)試截面絕對(duì)氣流角沿徑向分布曲線圖，對(duì)比了有無(wú)探針情況下在測(cè)試截面位置（0.5b）處的沿徑向平均分布的絕對(duì)氣流角曲線?？梢娞结樀募尤胧沟迷镜臍饬鹘钦w減小了2-5°，在探針根部位置，受探針尖部間隙及繞流的影響，氣流角偏轉(zhuǎn)程度最大。倘若氣流角偏轉(zhuǎn)值超過(guò)了下游低壓渦輪的不敏感攻角，將對(duì)發(fā)動(dòng)機(jī)性能產(chǎn)生影響。

圖9和圖10為測(cè)試截面絕對(duì)總壓和等熵效率沿徑向分布曲線圖，探針支桿的存在使渦輪流場(chǎng)的總壓減少，及流場(chǎng)可利用能量的損耗，在0.05-0.95葉高的范圍內(nèi)，總壓損失較大。同樣的，探針支桿的存在使，在探針截面位置的等熵效率同樣減少，在0.2-0.4和0.65-0.95葉高范圍內(nèi)的等熵效率均有所損失。

2.2 流熱固耦合分析

為了研究探針在測(cè)試截面位置的實(shí)際受力變形程度，其靜強(qiáng)度和動(dòng)強(qiáng)度是否達(dá)標(biāo)，由于渦扇發(fā)動(dòng)機(jī)實(shí)際運(yùn)行中不會(huì)有大的結(jié)構(gòu)形變，因此采用單向流熱固耦合模擬。同時(shí)采用探針支桿和帶6個(gè)等環(huán)分布延伸測(cè)孔的探針進(jìn)行模態(tài)分析，校驗(yàn)安全裕度。

為方便研究，渦輪導(dǎo)葉及轉(zhuǎn)子均不考慮冷卻裝置。圖11是導(dǎo)葉葉片及轉(zhuǎn)子葉片在氣流沖擊和對(duì)流換熱情況下的總應(yīng)變?cè)茍D，其中最大應(yīng)變發(fā)生在導(dǎo)葉尾緣中徑位置，約有1.3mm的形變，可得其撓度系數(shù)=2.2%<4%，而圖12是兩種探針的應(yīng)變?cè)茍D，可見有實(shí)際的延伸測(cè)點(diǎn)對(duì)于探針總的變形影響不大，均滿足靜強(qiáng)度的使用要求。

根據(jù)探針實(shí)際裝配狀態(tài)，對(duì)探針梳進(jìn)行模態(tài)分析，設(shè)置探針頂部為固定約束?？捎上率叫：颂结樀膭?dòng)強(qiáng)度，即

式中：nf為安全裕度，f為激振頻率，fn為自振頻率（n取1，2，3）。

由于探針承受渦輪轉(zhuǎn)子的出口氣流沖擊，因此激振頻率可等效為轉(zhuǎn)子的轉(zhuǎn)動(dòng)頻率為220.533Hz。而經(jīng)計(jì)算探針前三階自振頻率分別為938.11、939.21、5729.8Hz，顯而易見探針安全裕度均大于15%的要求，滿足使用的條件。

3 ?結(jié)論

①在轉(zhuǎn)子出口測(cè)試截面的探針不會(huì)影響上游的導(dǎo)葉及轉(zhuǎn)子氣動(dòng)參數(shù)，但會(huì)使渦輪總體的等熵效率減小，流量減小。②探針的存在使流場(chǎng)總壓損失增加，絕對(duì)出口氣流角偏轉(zhuǎn)2-5°。探針產(chǎn)生低馬赫數(shù)尾跡區(qū)域，探針尾跡脫渦與主流混合有摻混損失，且使得轉(zhuǎn)子出口流場(chǎng)參數(shù)周向分布不均勻。③在氣流沖擊和高溫環(huán)境下，導(dǎo)葉在尾緣中徑位置產(chǎn)生最大應(yīng)變，探針的動(dòng)靜強(qiáng)度均符合一般要求，無(wú)共振隱患。

參考文獻(xiàn)：

[1]鄒正平.航空燃?xì)廨啓C(jī)渦輪氣體動(dòng)力學(xué)：流動(dòng)機(jī)理及氣動(dòng)設(shè)計(jì)：Turbine aerodynamics for aero-engine： flow analysis and aerodynamics design[M].上海交通大學(xué)出版社，2014.

[2]向宏輝，高杰，楊榮菲，劉氦旭，吳森林.探針支桿尾跡對(duì)壓氣機(jī)轉(zhuǎn)子葉片振動(dòng)特性的影響研究[J].實(shí)驗(yàn)流體力學(xué)，2021，35（02）：58-66.

[3]楊榮菲，劉氦旭，向宏輝，高杰，王靖宇.進(jìn)口探針支桿尾跡誘發(fā)壓氣機(jī)轉(zhuǎn)子葉片共振的數(shù)值研究[J].推進(jìn)技術(shù)，2021，42（05）：1002-1012.

[4]劉氦旭.進(jìn)口探針誘發(fā)壓氣機(jī)葉片振動(dòng)現(xiàn)象的數(shù)值研究[D].南京航空航天大學(xué)，2020.

[5]Wyler J S. Probe blockage effects in free jets and closed tunnels[J]. 1975.

[6]羅距.高性能風(fēng)扇/壓氣機(jī)三維葉片氣動(dòng)設(shè)計(jì)與實(shí)驗(yàn)研究[D].南京航空航天大學(xué)，2013.

[7]朱銘敏，羌曉青，居振州，馬宇晨，滕金芳.四級(jí)低速大尺寸壓氣機(jī)的試驗(yàn)與數(shù)值模擬研究[J].工程熱物理學(xué)報(bào)，2021，42（05）：1150-1160.

[8]HB 7081-2012，航空燃?xì)鉁u輪發(fā)動(dòng)機(jī)軸流渦輪氣動(dòng)性能試驗(yàn)方法[S].

[9]Frederick S，Jeffrey W T，Henry H A. Effect of Number of Probes and Their Orientation on the Calculation of Several Compressor Face Distortion Descriptors[R]. NASATM-72859，1979.

[10]付少林，楊榮菲，劉長(zhǎng)青，等.測(cè)試布局對(duì)渦輪效率的影響研究[J].推進(jìn)技術(shù)，2019，40（11）：2464-2472.

[11]張小慶，曾來(lái)榮，劉偉雄.“總溫探針在脈沖燃燒風(fēng)洞中的應(yīng)用”.第十五屆全國(guó)激波與激波管學(xué)術(shù)會(huì)議論文集（下冊(cè)）.Ed.， 2012，46-50.