桑則林，白楚楓，劉作宏，閆衛(wèi)青，柳國(guó)印，徐志峰

（中國(guó)航發(fā)沈陽(yáng)發(fā)動(dòng)機(jī)研究所，遼寧沈陽(yáng) 110015）

隨著現(xiàn)代戰(zhàn)斗機(jī)對(duì)高機(jī)動(dòng)性能要求的提高，惡劣的飛行環(huán)境和大迎角飛行動(dòng)作等都會(huì)造成發(fā)動(dòng)機(jī)進(jìn)口氣動(dòng)測(cè)量截面（Aerodynamic Interface Plane，AIP）的流場(chǎng)壓力參數(shù)嚴(yán)重不均，降低了發(fā)動(dòng)機(jī)的可用穩(wěn)定余量，既會(huì)對(duì)發(fā)動(dòng)機(jī)的工作穩(wěn)定性造成破壞，又可能對(duì)發(fā)動(dòng)機(jī)結(jié)構(gòu)的完整性和壽命造成嚴(yán)重影響［1-2］。進(jìn)氣壓力畸變對(duì)航空發(fā)動(dòng)機(jī)穩(wěn)定性的影響是整個(gè)發(fā)動(dòng)機(jī)研制過(guò)程中一個(gè)關(guān)鍵的問(wèn)題，其會(huì)造成氣流參數(shù)徑向分布的變化，影響模態(tài)波擾動(dòng)和失速團(tuán)在動(dòng)態(tài)失速的三維擴(kuò)散過(guò)程，造成風(fēng)扇和壓氣機(jī)因葉片受力不均和顫動(dòng)而引起氣動(dòng)性能和結(jié)構(gòu)整體性下降［3-4］。進(jìn)氣畸變條件下壓氣機(jī)葉片的激振源更豐富，同一共振轉(zhuǎn)速下葉片的振動(dòng)幅值比自然進(jìn)氣條件下大［5］，因此研究進(jìn)氣畸變對(duì)發(fā)動(dòng)機(jī)的影響具有重要意義?？赏ㄟ^(guò)數(shù)值模擬［6］和建立模擬試驗(yàn)裝置［7］等方法開(kāi)展航空發(fā)動(dòng)機(jī)進(jìn)氣畸變研究，Hale等［8］基于非定常歐拉方程數(shù)值計(jì)算方法完成了單轉(zhuǎn)子均勻進(jìn)氣條件下壓氣機(jī)級(jí)的性能分析，發(fā)展了進(jìn)氣畸變對(duì)發(fā)動(dòng)機(jī)穩(wěn)定性影響的理論模型；芮長(zhǎng)勝等［9］通過(guò)對(duì)軸流壓氣機(jī)總壓畸變?cè)囼?yàn)數(shù)據(jù)進(jìn)行模擬計(jì)算，預(yù)測(cè)出該軸流壓氣機(jī)畸變條件下的失速首發(fā)級(jí)；張興發(fā)等［10］對(duì)軸流壓氣機(jī)與插板畸變發(fā)生器的耦合進(jìn)行數(shù)值仿真，得到了不同轉(zhuǎn)速下影響壓氣機(jī)穩(wěn)定邊界的主要因素；桑增產(chǎn)等［11］利用插板式畸變發(fā)生器試驗(yàn)研究了進(jìn)氣總壓畸變對(duì)某型渦噴發(fā)動(dòng)機(jī)氣動(dòng)穩(wěn)定性的影響；胡駿等［12］對(duì)某型大涵道比渦扇發(fā)動(dòng)機(jī)穩(wěn)定性的影響進(jìn)行了總壓畸變和總溫畸變?cè)u(píng)估，得到了臨界畸變指標(biāo)和發(fā)動(dòng)機(jī)首發(fā)失穩(wěn)級(jí)組；張亞維等［13］通過(guò)設(shè)計(jì)進(jìn)氣道綜合監(jiān)測(cè)系統(tǒng)，分析了進(jìn)氣畸變與發(fā)動(dòng)機(jī)實(shí)際工況的匹配關(guān)系；孫鵬等［14］通過(guò)數(shù)值分析研究了插板后非均勻流場(chǎng)的變化，得出規(guī)律性的壓氣機(jī)前總壓分布形式。從工程應(yīng)用實(shí)際出發(fā)，可移動(dòng)式擾流板因其設(shè)備簡(jiǎn)單、實(shí)用性強(qiáng)等特點(diǎn)成為被廣泛應(yīng)用的畸變發(fā)生裝置［15-16］。在理論體系研究方面，我國(guó)已形成了以綜合壓力畸變指數(shù)來(lái)預(yù)測(cè)并建立進(jìn)口壓力畸變對(duì)發(fā)動(dòng)機(jī)影響的評(píng)定體系，其由穩(wěn)態(tài)周向總壓畸變指數(shù)和動(dòng)態(tài)總壓畸變指數(shù)組合而成［17］。

依據(jù)大量畸變?cè)囼?yàn)數(shù)據(jù)的積累，隨發(fā)動(dòng)機(jī)狀態(tài)和插板深度動(dòng)態(tài)變化的畸變指數(shù)成為研究畸變?cè)囼?yàn)的重要參數(shù)，擋板類(lèi)畸變發(fā)生裝置產(chǎn)生的畸變指數(shù)的準(zhǔn)確性是研究畸變?cè)囼?yàn)對(duì)發(fā)動(dòng)機(jī)影響的重要前提，本文采用插板式壓力畸變發(fā)生器，對(duì)某航空發(fā)動(dòng)機(jī)壓力畸變?cè)囼?yàn)過(guò)程中AIP截面參數(shù)進(jìn)行分析總結(jié)，探究壓力畸變的影響因素，為航空發(fā)動(dòng)機(jī)進(jìn)氣壓力畸變?cè)囼?yàn)提供借鑒。

1 畸變發(fā)生器

在發(fā)動(dòng)機(jī)進(jìn)氣口前通過(guò)安裝插板式畸變發(fā)生器產(chǎn)生壓力畸變。插板裝置通常分為液壓驅(qū)動(dòng)和伺服電機(jī)驅(qū)動(dòng)兩種，當(dāng)進(jìn)氣道尺寸較小時(shí)一般選用伺服電機(jī)驅(qū)動(dòng)插板裝置作為畸變發(fā)生器，其結(jié)構(gòu)組成如圖1 所示。采用伺服電機(jī)、直線(xiàn)推桿驅(qū)動(dòng)方式，可控制插板停留到中間的任何位置且可控制插板按指定速度移動(dòng)，實(shí)現(xiàn)插板對(duì)發(fā)動(dòng)機(jī)進(jìn)氣口的部分遮擋進(jìn)而引起發(fā)動(dòng)機(jī)進(jìn)氣口壓力畸變。作動(dòng)筒內(nèi)設(shè)有線(xiàn)位移傳感器，可向測(cè)控系統(tǒng)實(shí)時(shí)反饋插板的位置，借此對(duì)插板位置實(shí)現(xiàn)閉環(huán)反饋控制。

圖1 伺服電動(dòng)驅(qū)動(dòng)插板裝置結(jié)構(gòu)組成

為保證試驗(yàn)安全，可將機(jī)載壓差判喘系統(tǒng)及執(zhí)行機(jī)構(gòu)、臺(tái)架判喘系統(tǒng)與臺(tái)架畸變發(fā)生器的插板裝置進(jìn)行連鎖設(shè)置，畸變?cè)囼?yàn)聯(lián)合防喘系統(tǒng)示意圖如圖2 所示。機(jī)載喘振壓差傳感器識(shí)別到喘振信號(hào)后，機(jī)載控制器執(zhí)行消喘邏輯，同時(shí)將喘振信號(hào)發(fā)送至畸變發(fā)生器，插板裝置執(zhí)行應(yīng)急退板；臺(tái)架判喘系統(tǒng)識(shí)別到喘振信號(hào)后，將喘振信號(hào)發(fā)送給機(jī)載控制器和畸變發(fā)生器，機(jī)載控制器執(zhí)行消喘邏輯，插板裝置執(zhí)行應(yīng)急退板。

圖2 畸變?cè)囼?yàn)聯(lián)合防喘系統(tǒng)示意圖

2 參數(shù)測(cè)量和處理方法

發(fā)動(dòng)機(jī)進(jìn)氣壓力畸變?cè)囼?yàn)測(cè)量截面示意圖如圖3所示，進(jìn)氣測(cè)量截面包括：0 -0 截面、1 -1 截面和AIP截面，各截面的位置如下。

圖3 發(fā)動(dòng)機(jī)進(jìn)氣壓力畸變?cè)囼?yàn)測(cè)量截面示意圖

①0 -0 截面主要用于測(cè)量進(jìn)氣流量，位于插板前2D處，布置有4 支×4 點(diǎn)周向均布的總壓受感部和4 個(gè)周向均布的壁面靜壓測(cè)孔，0 -0 截面受感部分布示意圖（逆航向）如圖4 所示。

圖4 0 -0截面受感部分布示意圖（逆航向）

②A(yíng)IP截面用于測(cè)量和計(jì)算進(jìn)氣穩(wěn)態(tài)周向總壓畸變指數(shù)，位于發(fā)動(dòng)機(jī)進(jìn)氣機(jī)匣前安裝邊前150 mm處，布置有6 支×5 點(diǎn)周向均布穩(wěn)態(tài)總壓測(cè)量受感部和6 個(gè)周向均布的壁面靜壓測(cè)壓孔，穩(wěn)態(tài)總壓測(cè)量受感部測(cè)試精度為±0.3%。

③1 -1 截面用于測(cè)量進(jìn)氣紊流度，在距AIP測(cè)量截面前85 mm處，1 -1 截面在0.9R 環(huán)面周向均布6支高頻總壓脈動(dòng)測(cè)量傳感器，其頻率響應(yīng)＞3.5 kHz，測(cè)試誤差＜1%，整個(gè)系統(tǒng)的計(jì)算誤差＜10%，1 -1 截面總壓脈動(dòng)測(cè)耙周向分布示意圖（逆航向）如圖5所示。

圖5 1 -1截面總壓脈動(dòng)測(cè)耙周向分布示意圖（逆航向）

④可移動(dòng)插板位于A(yíng)IP 截面前3D 處，插板式壓力畸變發(fā)生器簡(jiǎn)圖（逆航向）如圖6 所示。

圖6 插板式壓力畸變發(fā)生器簡(jiǎn)圖（逆航向）

進(jìn)氣壓力畸變相關(guān)氣動(dòng)參數(shù)計(jì)算公式及定義如下。AIP截面上單個(gè)測(cè)點(diǎn)總壓恢復(fù)系數(shù)σ（，θ）為

式中：r- ＝r／rti為測(cè)點(diǎn)半徑與輪緣半徑之比。

周向位置的徑向平均總壓恢復(fù)系數(shù)σr，av（θ）為

低壓區(qū)內(nèi)平均總壓恢復(fù)系數(shù)σ0為

式中：θ1～θ2為周向低壓區(qū)范圍，θ1和θ2分別為低壓區(qū)的開(kāi)始角和結(jié)束角。

面平均總壓恢復(fù)系數(shù)σav為

動(dòng)態(tài)總壓脈動(dòng)均方根值（ΔP1-1D）RMS為

周向θ位置的紊流度ε（θ）為

面平均紊流度εav為

綜合壓力畸變指數(shù)W為

3 壓力畸變?cè)囼?yàn)

航空發(fā)動(dòng)機(jī)進(jìn)氣壓力畸變?cè)囼?yàn)表明：總壓恢復(fù)系數(shù)的量值、不均勻度與插板深度密切相關(guān)，擾動(dòng)氣流參數(shù)隨著插板相對(duì)深度的增加而急劇增大［18］。發(fā)動(dòng)機(jī)在不同狀態(tài)下穩(wěn)態(tài)周向壓力畸變指數(shù)、面平均紊流度（又稱(chēng)“動(dòng)態(tài)畸變指數(shù)”）、綜合壓力畸變指數(shù)、穩(wěn)態(tài)周向壓力畸變指數(shù)與動(dòng)態(tài)畸變指數(shù)的比值（以下簡(jiǎn)稱(chēng)“穩(wěn)動(dòng)態(tài)比例”）與插板深度的關(guān)系如圖7 所示，其中狀態(tài)1 ～狀態(tài)3 為利用CFD軟件進(jìn)行穩(wěn)態(tài)和動(dòng)態(tài)數(shù)值模擬計(jì)算的結(jié)果。按照實(shí)際工作條件建立模型，通過(guò)給定的描述穩(wěn)態(tài)和動(dòng)態(tài)總壓畸變的邊界條件進(jìn)行模擬計(jì)算，測(cè)量截面沿周向均布的30 個(gè)點(diǎn)，徑向按照實(shí)際測(cè)點(diǎn)位置布置的5 個(gè)點(diǎn)，按式（1）～式（5）計(jì)算獲得不同插板深度下穩(wěn)態(tài)畸變指數(shù)；動(dòng)態(tài)數(shù)據(jù)以獲得的總壓脈動(dòng)數(shù)據(jù)為基礎(chǔ)，結(jié)合穩(wěn)動(dòng)態(tài)比例，通過(guò)試驗(yàn)數(shù)據(jù)校核，得出計(jì)算模擬下的動(dòng)態(tài)畸變指數(shù)，其中狀態(tài)1 為0.8 轉(zhuǎn)速狀態(tài)，不同狀態(tài)之間為相差2%的遞增穩(wěn)態(tài)轉(zhuǎn)速臺(tái)階。在相同插板深度下，隨著發(fā)動(dòng)機(jī)狀態(tài)的提高，穩(wěn)態(tài)畸變和動(dòng)態(tài)畸變指數(shù)變大；在相同發(fā)動(dòng)機(jī)狀態(tài)下，隨著插板深度的增加，綜合畸變指數(shù)按照一定的比例關(guān)系逐漸增大。

圖7 畸變指數(shù)與插板深度的關(guān)系

與相同技術(shù)狀態(tài)和進(jìn)氣口條件下的數(shù)值模擬計(jì)算結(jié)果相比，基于整機(jī)畸變調(diào)試試驗(yàn)過(guò)程穩(wěn)態(tài)數(shù)據(jù)融合處理所獲得的AIP 截面參數(shù)［19］計(jì)算得到的綜合畸變指數(shù)偏低（以圖7 所示的壓力畸變?cè)囼?yàn)結(jié)果為例，試驗(yàn)轉(zhuǎn)速與狀態(tài)2 計(jì)算轉(zhuǎn)速一致），相同插板深度下調(diào)試階段穩(wěn)態(tài)壓力畸變指數(shù)低于狀態(tài)1 下的數(shù)據(jù)，動(dòng)態(tài)畸變指數(shù)高于狀態(tài)1 的數(shù)據(jù)，但仍低于相同發(fā)動(dòng)機(jī)狀態(tài)2 數(shù)值模擬預(yù)估數(shù)值，因此導(dǎo)致調(diào)試階段綜合畸變指數(shù)與預(yù)估值之間的偏差最大可達(dá)18%。由穩(wěn)動(dòng)態(tài)比例可見(jiàn)，調(diào)試過(guò)程中隨著插板深度逐漸增大，穩(wěn)動(dòng)態(tài)比例偏離預(yù)估值也越高。選取插板相對(duì)深度H ＝25%時(shí)測(cè)量截面參數(shù)進(jìn)行詳細(xì)分析，由圖7 可見(jiàn)畸變調(diào)試過(guò)程中各畸變特征指數(shù)和穩(wěn)動(dòng)態(tài)比例與數(shù)值模擬獲得的預(yù)估數(shù)值均存在較大偏差，須對(duì)試驗(yàn)數(shù)據(jù)和試驗(yàn)結(jié)果進(jìn)行詳細(xì)分析，探究進(jìn)氣壓力畸變?cè)囼?yàn)的影響因素。

3.1 畸變?cè)O(shè)備整體結(jié)構(gòu)的影響

在相同狀態(tài)下，AIP 截面上總壓恢復(fù)系數(shù)沿徑向分布（同一半徑上6 點(diǎn)測(cè)量平均值）和AIP 截面上總壓恢復(fù)系數(shù)沿周向分布（同一角度位置5 點(diǎn)測(cè)量平均值）如圖8 所示。由圖8 可見(jiàn)，由于平板結(jié)構(gòu)相對(duì)于水平直徑的對(duì)稱(chēng)性，插板后面壓力場(chǎng)也是周向?qū)ΨQ(chēng)的，進(jìn)而證實(shí)了在不同板位下壓力場(chǎng)分布相應(yīng)地于180°位置周向?qū)ΨQ(chēng)。結(jié)合徑向和周向壓力分布可見(jiàn)，兩組壓力參數(shù)低壓區(qū)位置和徑向分布形式基本一致，綜合考慮測(cè)量和計(jì)算誤差等因素，低壓區(qū)域壓力測(cè)量值偏高是導(dǎo)致整體畸變指數(shù)出現(xiàn)較大差異的主要原因。

圖8 AIP截面上總壓恢復(fù)系數(shù)分布情況

依據(jù)AIP截面總壓云圖（圖9）可以更加明顯地看到AIP截面上壓力分布沿直徑方向上下對(duì)稱(chēng)分布，對(duì)比畸變?cè)囼?yàn)和數(shù)值模擬可以看出相同截面下高壓區(qū)域范圍基本一致，低壓區(qū)域壓力水平和影響范圍不同造成畸變量的差異。從而可以判斷，畸變發(fā)生設(shè)備整體結(jié)構(gòu)對(duì)畸變結(jié)果的影響較大，需對(duì)插板位置和畸變發(fā)生器結(jié)構(gòu)進(jìn)行梳理和檢查。

圖9 AIP截面總壓云圖

通過(guò)重新標(biāo)定插板位置，可以確認(rèn)兩次對(duì)比數(shù)據(jù)插板達(dá)到了相同的位置。對(duì)畸變發(fā)生器和進(jìn)氣道之間的結(jié)構(gòu)進(jìn)行檢查：如圖10 所示，進(jìn)氣道與插板前轉(zhuǎn)接段之間、轉(zhuǎn)接段和測(cè)量段之間均采用螺栓壓緊結(jié)構(gòu)，此結(jié)構(gòu)可保證兩側(cè)連接的緊實(shí)度。插板后轉(zhuǎn)接段與轉(zhuǎn)接段之間采用密封轉(zhuǎn)接段的形式進(jìn)行連接，密封轉(zhuǎn)接段為兩段獨(dú)立的轉(zhuǎn)接段采用插接方式進(jìn)行連接，插接后以塑膠進(jìn)行密封填補(bǔ)，以便調(diào)節(jié)插板距離、AIP截面的距離和適應(yīng)發(fā)動(dòng)機(jī)軸向的移動(dòng)。與螺栓連接相比，插接結(jié)構(gòu)更容易造成進(jìn)氣轉(zhuǎn)接段漏氣，使得低壓區(qū)壓力水平偏離預(yù)估值，通過(guò)檢查發(fā)現(xiàn)密封轉(zhuǎn)接段氣密性差是造成整體畸變指數(shù)偏低的重要因素。

圖10 畸變發(fā)生器與發(fā)動(dòng)機(jī)連接結(jié)構(gòu)圖

3.2 測(cè)點(diǎn)位置及計(jì)算方法的影響

經(jīng)過(guò)調(diào)整發(fā)動(dòng)機(jī)進(jìn)氣道密封轉(zhuǎn)接段，檢查確認(rèn)密封性后在發(fā)動(dòng)機(jī)達(dá)到狀態(tài)2 的轉(zhuǎn)速和插板深度下進(jìn)行畸變?cè)囼?yàn)?；冎笖?shù)相對(duì)比例值如表1 所示，可見(jiàn)在相同插板深度下動(dòng)態(tài)畸變指數(shù)達(dá)到預(yù)估水平，穩(wěn)態(tài)周向壓力畸變指數(shù)與預(yù)估數(shù)值相比偏低，導(dǎo)致綜合畸變指數(shù)存在接近9%的偏差。通過(guò)對(duì)比圖11 中試驗(yàn)總壓恢復(fù)系數(shù)的分布情況，可見(jiàn)不同方向壓力分布情況相近，須對(duì)兩次試驗(yàn)過(guò)程的穩(wěn)態(tài)畸變算法進(jìn)行校驗(yàn)。

表1 畸變指數(shù)相對(duì)值對(duì)比

圖11 狀態(tài)2下AIP截面上總壓恢復(fù)系數(shù)分布情況

在臺(tái)架畸變?cè)囼?yàn)中，為準(zhǔn)確掌握畸變指數(shù)隨插板深度和發(fā)動(dòng)機(jī)狀態(tài)變化的關(guān)系，同時(shí)避免測(cè)量受感部過(guò)多造成的誤差干擾，AIP截面測(cè)量?jī)H采用6 支×5 點(diǎn)周向均布的“水”字耙進(jìn)行穩(wěn)態(tài)總壓參數(shù)測(cè)量，如圖12所示。由1.2 節(jié)計(jì)算公式可見(jiàn)，穩(wěn)態(tài)周向壓力畸變指數(shù)取決于低壓區(qū)平均總壓恢復(fù)系數(shù)和面平均總壓恢復(fù)系數(shù)，其中主要因素為周向θ 位置對(duì)應(yīng)的總壓恢復(fù)系數(shù)及周向測(cè)點(diǎn)個(gè)數(shù)；周向測(cè)點(diǎn)個(gè)數(shù)的影響主要在于隨著低壓區(qū)周向測(cè)點(diǎn)個(gè)數(shù)的增多，計(jì)算得到的低壓區(qū)平均總壓恢復(fù)系數(shù)減小，在相同的面平均總壓系數(shù)下，穩(wěn)態(tài)周向壓力畸變指數(shù)增大。在數(shù)值模擬時(shí)，AIP 截面上提取的數(shù)據(jù)點(diǎn)個(gè)數(shù)遠(yuǎn)多于試驗(yàn)時(shí)的30 個(gè)穩(wěn)態(tài)數(shù)據(jù)點(diǎn)，因此處理試驗(yàn)數(shù)據(jù)時(shí)需細(xì)化整個(gè)測(cè)量截面壓力參數(shù)，才能更好地反映低壓區(qū)的影響范圍。根據(jù)總壓恢復(fù)系數(shù)沿周向分布及總壓云圖數(shù)據(jù)，AIP 截面上壓力參數(shù)相應(yīng)地于180°徑向位置對(duì)稱(chēng)分布，同時(shí)測(cè)點(diǎn)布置沿180°徑向位置非對(duì)稱(chēng)分布，因此可將測(cè)量參數(shù)進(jìn)行上下鏡像對(duì)稱(chēng)處理，這對(duì)低壓區(qū)域平均總壓恢復(fù)系數(shù)的計(jì)算起到重要作用。經(jīng)過(guò)鏡像測(cè)點(diǎn)處理后，測(cè)點(diǎn)位置如表2 所示，畸變指數(shù)可達(dá)到預(yù)估水平，能夠更好地反映AIP截面上的畸變水平。

表2 兩種數(shù)據(jù)分析方法下周向測(cè)點(diǎn)分布

圖12 AIP截面受感部布局示意圖（逆航向）

4 結(jié)束語(yǔ)

本文對(duì)影響航空發(fā)動(dòng)機(jī)進(jìn)氣畸變?cè)囼?yàn)的因素進(jìn)行了研究，通過(guò)對(duì)比分析不同數(shù)據(jù)下壓力參數(shù)和畸變指數(shù)，可得到以下結(jié)論。

①畸變強(qiáng)度的大小可由低壓區(qū)域壓力水平和影響面積間接反映，畸變發(fā)生器的連接氣密性可導(dǎo)致試驗(yàn)結(jié)果與理論值的綜合畸變指數(shù)最大偏差達(dá)到18%。

②周向測(cè)點(diǎn)個(gè)數(shù)影響穩(wěn)態(tài)周向壓力畸變指數(shù)的計(jì)算值大小，采用鏡像對(duì)稱(chēng)處理的方法細(xì)化測(cè)量截面參數(shù)，能夠使試驗(yàn)環(huán)境下的畸變指數(shù)計(jì)算值提高接近9%，更好地反映壓力畸變水平。