進氣道//發動機相容性試飛中總壓畸變數據的采集與處理

汪濤，姜健，史建邦

進氣道//發動機相容性試飛中總壓畸變數據的采集與處理

汪濤，姜健，史建邦

(中國飛行試驗研究院發動機所，陜西西安710089)

研究了進氣道/發動機相容性試飛中總壓畸變數據的采集與處理方法，詳細介紹了傳感器選型、測點布局、受感部設計、測試及采集系統構建的要點，開發了通用動態數據處理軟件，并重點研究了動態總壓數據處理方法和參數選擇問題。通過某型飛機/發動機的進氣道/發動機相容性試飛應用及其數據處理表明：提出的測試與采集方案合理可行，工程應用情況良好；動態數據處理中應合理選擇截止頻率、濾波器階次以及濾波函數。

進氣道/發動機相容性；總壓畸變；紊流度；濾波

1 引言

航空發動機的氣動研究理論上都是在均勻進氣條件下進行，即認為進氣道出口氣流參數均勻。但在實際情況中，航空發動機無一例外地必須在非均勻進氣（即進氣畸變）條件下工作。在早期的亞聲速、小機動飛行條件下，設計人員主要重視其性能的好壞，追求高的性能指標，而對進氣畸變的影響，即對進氣道/發動機相容性問題未給予足夠重視。但隨著飛機飛行速度和高度的不斷提高、機動性的不斷增加，進氣畸變的影響越來越突出，在實際飛行中暴露出來的相容性問題也越來越嚴重。在此背景下，自20世紀70年代以來，美、英及蘇聯等國大力研究了進氣流場畸變引起的進氣道/發動機相容性問題[1]。由于進氣流場壓力畸變（包括穩態壓力畸變和動態壓力畸變）主要由飛機側滑、大攻角飛行時引起，這些都會影響發動機穩定工作，因此，在進氣道/發動機相容性試飛中，對進氣總壓畸變的測量尤為關鍵[2～9]。

在進氣道/發動機相容性試飛中，通過測量發動機氣動界面上氣流的總壓分布，來描述發動機進口流場的畸變程度?？倝夯兊臏y量和采集，測量數據尤其是動態總壓畸變數據的處理，是其核心問題。其中涵蓋了傳感器選型、測點布局、受感部設計、測試及采集方法、穩態及動態數據處理方法等技術細節問題。

本文結合某型飛機/發動機的進氣道/發動機相容性試飛，從測試需求、系統構建直至數據處理，逐一進行了研究，尤其是對動態數據的采集和處理方法進行了較為深入的分析，并開發了通用數據處理軟件。

2 傳感器選型及受感部設計

總壓畸變測量，是通過進氣道/發動機氣動界面上氣流總壓分布來描述發動機進口流場。實際中，由于無法測量出整個氣動界面上的氣流參數，故采用抽樣測量幾十個參數的方法來代替。目前，俄標多采用“水”字型測量耙，布置30點穩態總壓和6點動態總壓進行測量；而美標多采用“米”字型測量耙，布置40/48點穩態總壓及40/48點動態總壓進行測量。我國在此領域對俄羅斯的合作和繼承較多，一般采用俄羅斯所使用的方法。為減小測量誤差，測量布局時應參考早期數值計算/風洞試驗結果，盡可能使耙布置在最低/最高壓力區[10]。

在傳感器選型上，穩態總壓測量選用絕對壓力傳感器，測量范圍根據飛機的高度-速度包線確定。動態總壓測量選用固有頻率較高的傳感器，且當氣流脈動頻率高于發動機轉子頻率4倍時，上游氣流畸變對發動機穩定性幾乎無影響[2]，為此全尺寸動態總壓探頭的頻響應為1 000 Hz。因此，在傳感器選擇上，傳感器的固有頻率至少要比對發動機穩定性有影響的最高頻率高2～3倍，以避免信號失真，保證傳感器的平坦響應特性。動態傳感器應靠近測點安裝，以避免由于安裝問題而影響所需頻率的測量。

根據上述原則，本文折中采用了“水”字型測量耙(如圖1所示)，6支耙等周向間距布置，每支耙上的5點穩態總壓和5點動態總壓沿半徑方向等環面布置。穩態壓力傳感器選用了KYCO8A型絕對壓力傳感器，其測量范圍為0～170 kPa。動態總壓傳感器選用了KULITE公司的XCEL-70-152-25A型壓電式傳感器，其測量范圍為0～170 kPa絕壓，固有頻率大于3 000 Hz，適用環境溫度范圍為-50～80℃。為滿足安裝要求，對動態壓力傳感器進行了封裝，如圖2所示。

圖1 總壓畸變測量耙Fig.1 Total pressure distortion measurement rake

圖2 動態壓力傳感器Fig.2 Dynamic total pressure sensor

測量耙上包含了穩態受感部和封裝過的動態壓力傳感器，測量耙為懸臂支板結構[11]，內側整體穿有支撐環，安全性較好，堵塞比小于5%。制造過程中，首先加工一套振動試驗件，并按文獻[12]中有關規定進行振動試驗，振動嚴酷等級選取G級。振動試驗及X射線檢查合格后，按與試驗件相同的技術條件制造裝機用測量耙。

3 測試及采集系統架構

試飛中，測試及采集系統架構非常龐大，一架飛機的測試參數一般在四千個左右。由于進氣道總壓畸變參數有其特殊性，按照分布式測試系統構架，本文專門研制了進氣道畸變子系統進行測試和采集，如圖3所示。該系統由穩態、動態壓力傳感器，兩臺中等增益快變參數調節器，兩臺慢變參數調節器，進氣道畸變采集器和數字記錄器組成。

圖3 測試及采集系統架構圖Fig.3 Measuring and acquisition system

壓力耙上穩態壓力特性為：0～180 kPa，采樣率取32 Hz，穩態壓力從受感部后端引出管處通過紫銅管或氟塑料管引至壓力傳感器處。動態壓力特性為：0～180 kPa，0～400 Hz(即有影響的最高頻率為400 Hz，超過400 Hz的脈動氣流對本文發動機穩定性無影響，該頻率由發動機前期地面試驗得到)，采樣率取1 000 Hz。該系統加裝30個KYCO8A絕對壓力傳感器，傳感器輸出經慢變參數調節器進行放大、濾波處理；快變參數調節器對30個KULITE動態壓力傳感器提供橋壓，并對輸出動態壓力信號進行放大、濾波，然后與慢變參數調節器輸出的信號一并送至進氣道畸變采集器進行采集、編碼，以PCM串行數據流的形式送至記錄器進行記錄。

在采集系統研制過程中還驗證了采集系統與動、穩態壓力傳感器的匹配性，動態壓力信號采集的準確性，以及各通道信號采集的相位一致性(要求400 Hz時相位差小于3°)，以上驗證均滿足要求。

4 飛行試驗數據處理及研究

應用上述方法，本文構建了測試及采集系統，對某型飛機/發動機進行了進氣道/發動機相容性試飛，主要的10種飛行工況如表1所示。表中H為飛行高度，Ma為飛行馬赫數、α為飛行攻角。

表1 飛行工況編號Table 1 Fight test conditions

試飛獲得了連續記錄的穩態和動態總壓數據，進氣道出口綜合畸變指數W由動態紊流度ε和穩態周向畸變指數δO相加得到。即：

穩態周向畸變由30點穩態總壓數據繪制圖譜(如圖4所示)并通過插值計算指數，其計算方法相對固定，在此不贅述。

圖4 穩態總壓畸變云圖Fig.4 Static total pressure distortion contour

動態數據處理較為復雜。動態壓力信號在采集中不可避免地混入了高頻成分，因此首先必須對數字信號進行濾波，濾波后再根據下式進行計算。

式中：n為動態總壓測量點數，m為一段時間內每個測量點所得到的數據樣本量，Ptd_ij為第i支動態壓力測量點第j時刻的動態壓力，Ptave_i為第i支穩態壓力測量點的平均壓力。

動態數據處理關鍵在于濾波，濾波后用式(2)計算紊流度。由于數字濾波方法并不固定，往往不同方法處理結果相差較大，因此本文編制了數字濾波軟件，軟件界面見圖5。濾波函數有五種可選，截止頻率、濾波器階次均可調。輸入、輸出功能完善，可輸入全部30列數據或其中幾列，可輸出濾波后的動態數據及每個測點對應的紊流度，也可根據式(2)直接給出面平均紊流度。本文利用該軟件重點研究了數字濾波方法，給出了數字濾波合理的參數選取。另外，對于隨機脈動數據，為獲得誤差不大于2%的壓力均方差估算值，一般要持續1～2 s的樣本容量。本文統一取飛行姿態穩定后的2 s時間樣本。

4.1 截止頻率的影響

對于采集到的動態總壓數據，選取合適的低通截止頻率十分重要，頻率過大或過小都得不到準確結果。對于進氣道采集到的動態數據，截止頻率的選取與兩因素有關：一是實際應用時發動機的轉速(即擾動頻率)，二是滿足采樣定律的要求。由于目前并無準確的理論公式可供參考，根據文獻[13]～[16]，截止頻率fc的工程估算值可采用下式計算：

圖5 開發的數字濾波軟件Fig.5 Software for dynamic data processing

式中：Hz為發動機擾動頻率(Hz)，Vmax為發動機進口處最大質量流速(m/s)，D1為發動機進口截面直徑(m)。

利用譜分析工具對動態數據信號進行分析表明：信號頻率成分主要集中在300 Hz以下，與截止頻率理論計算結果150 Hz≤fc≤380 Hz比較吻合。

確定信號的主要頻率成分后，利用本文軟件分別以50、100、150、200、250、300、350、400、450 Hz為低通截止頻率，以4階Blackman窗函數進行濾波，各工況的面平均紊流度如圖6所示。從圖中看，隨著低通截止頻率的提高，紊流度逐步上升，到400 Hz后趨于平坦，說明低通截止頻率取400 Hz合適。

4.2 濾波器階次的影響

圖7示出了濾波器階次對紊流度的影響。從圖中看，濾波器選取1階和3階時結果明顯偏小，2階時結果普遍偏大，只有在4階以上時結果才趨于平穩。因此，對于本文對象的試驗數據，濾波器階次應不小于4階。

4.3 濾波窗函數的影響

圖6 截止頻率對紊流度的影響Fig.6 Cut-off frequency vs turbulence

圖7 濾波器階次對紊流度的影響Fig.7 Filter order vs turbulence

圖8 濾波窗函數對紊流度的影響Fig.8 Filtering function vs turbulence

本文的濾波函數選用了窗函數，圖8給出了五種窗函數(Hanning窗、Hamming窗、Blackman窗、Tri?ang窗、Kaise窗)對濾波結果的影響。從圖中看，濾波窗函數對紊流度計算結果幾乎沒影響。

5 結論

(1)本文提出的測試與采集方案合理可行，工程應用情況良好。

(2)濾波參數選擇對動態數據處理非常關鍵，濾波函數對結果影響較小，但低通截止頻率和濾波器階次對結果影響很大。對于本文試驗，低通截止頻率應不低于400 Hz，濾波器階次應不小于4階。

[1]姜健，于芳芳，趙海剛，等.進氣道/發動機相容性評價體系的完善與發展[J].科學技術與工程，2009，9(21)：6474—6483.

[2]劉大響，葉培梁，胡俊，等.航空燃氣渦輪發動機穩定性設計與評定技術[M].北京：航空工業出版社，2004：3—13.

[3]吳虎，廉小純，陳輔群.動態壓力畸變下壓縮系統穩定性分析[J].推進技術，1996，17(4)：32—36.

[4]廉小純，陳輔群，王驪安.進氣總壓總溫畸變對壓氣機穩定性的影響[J].航空動力學報，1993，8(2)：105—108.

[5]吳虎，廉小純，陳輔群.穩態聯合畸變下多級軸流壓氣機的氣動穩定性分析[J].航空學報，1996，17(4)：452—455.

[6]Bowditch D N，Coltrin R E.A Survey of Inlet/Engine Dis?tortion Capability[R].AIAA 83-1166，1983.

[7]Willam T C.History,Philosophy,Physics,and Future Di?rections of Aircraft Propulsion System/Inlet Integration[R]. ASME GT2004-54210，2004.

[8]Rabe D，Bolcs A，Russler P.Influence of Inlet Distortion onTransonicCompressorBladeLoading[R].AIAA 95-2461，1995.

[9]Walsh K R，Yuhas A J，Williams J G，et al.Inlet Distor?tion for an F/A-18A Aircraft during Stabilized Maneuvers up to 60°Angle of Attack[R].NASA TM-104329，1997.

[10]Frederick S，Jeffrey W T，Henry H A.Effect of Number of Probes and Their Orientation on the Calculation of Several CompressorFaceDistortionDescriptors[R].NASA TM-72859，1979.

[11]Andrew J Y，Ronald J R，Richard R B，et al.Design and Development of an F/A-18 Inlet Distortion Rake-A Cost and Time Saving Solution[R].AIAA 94-2132-CP，1994.

[12]GJB150.16-86，軍用設備環境試驗方法——振動試驗[S].

[13]Beale D.An Acceptance Process for the Evaluation of In?let Distortion[R].AIAA 92-39182，1992.

[14]惠增宏，張理，竹朝霞，等.進氣道實驗中動態數據采集及處理方法研究[J].實驗力學，2002，17(1)：48—54.

[15]吳森林，顧衛群，任銘林.壓力畸變動態數據的數字濾波處理方法[J].燃氣渦輪試驗與研究，2007，20(1)：48—51.

[16]Norby W P，Ladd J A，Yuhas A J.Dynamic Inlet Distor?tion Prediction with a Combined Computational Fluid Dy?namics andDistortionSynthesis Approach[R].NASA CR-198053，1996.

Acquisition and Processing of Total Pressure Distortion Data in Inlet/Engine Compatibility Flight Test

WANG Tao，JIANG Jian，SHI Jian-bang
(
Engine Department of Chinese Flight Test Establishment，Xi’an 710089，China)

The methods on acquisition and processing of total pressure distortion data in inlet/engine com?patibility flight test were researched.The sensor selection,station distribution,the rake design,and data ac?quisition system construction were specified.General software for dynamic data processing was developed, data processing method and parameter selection were emphasized.Applying the methods to an inlet/engine compatibility flight test project,final results showed that the measuring and processing methods were rea?sonable and feasible,cut-off frequency and filter order and filtering function should be carefully selected in dynamic data processing.

inlet/engine compatibility；total pressure distortion；turbulivity；filtering

V211.48

A

1672-2620(2012)02-0054-05

2011-08-05；

2011-09-26

汪濤(1983-)，男，陜西安康人，工程師，碩士研究生，研究方向為進氣道/發動機相容性試飛。