張 鑫，范 明

（中國航發(fā)控制系統(tǒng)研究所，江蘇無錫 214000）

0 引言

隨著飛機對發(fā)動機性能要求的不斷提高，其工作包線不斷擴大，導(dǎo)致發(fā)動機喘振裕度不足；飛機的載彈量增加，發(fā)射時廢氣尾流對發(fā)動機影響更為強烈；飛機機動性要求更高，發(fā)動機工作條件更為惡劣。以上變化都對發(fā)動機防喘/消喘控制提出了更高要求。

為了提高發(fā)動機防喘/消喘控制的性能，國內(nèi)外均持續(xù)開展了各種防喘/消喘控制系統(tǒng)的研制與設(shè)計。Gaston提出了一種新的判喘方法，在保證防喘/消喘功能的同時盡量減小系統(tǒng)對發(fā)動機性能的影響；Gogoi等建立了一種基于CFD的軸流壓氣機旋轉(zhuǎn)失速和喘振等氣動不穩(wěn)定性預(yù)測模型，可作為初步研究壓縮系統(tǒng)失穩(wěn)的有效工具；張紹基針對戰(zhàn)斗機發(fā)射導(dǎo)彈時的熱擾動因素開展了詳細(xì)的防喘/消喘控制系統(tǒng)設(shè)計與優(yōu)化，并進(jìn)行了臺架驗證；邢家瑞等對發(fā)射武器的防喘/消喘控制開展了分類細(xì)化設(shè)計，進(jìn)行了地面、飛行臺吞噬火藥等多種驗證試驗；秦海波等提出根據(jù)高壓壓氣機轉(zhuǎn)速和發(fā)動機進(jìn)口壓力對判喘依據(jù)進(jìn)行動態(tài)修正，以提高判喘的準(zhǔn)確度；葉巍等對喘振過程中發(fā)動機高壓壓氣機總壓與靜壓的壓差變化情況及其脈動頻率進(jìn)行了研究與分析，并給出了判喘系統(tǒng)的評估方法；孫健國等研究了喘振裕度自適應(yīng)控制方法以及該控制方法對發(fā)動機性能的提升作用；張海波等建立了某型渦扇發(fā)動機喘振實時模型，提出了一種主動穩(wěn)定性控制方法，提高了發(fā)動機過渡態(tài)的動態(tài)性能；李小彪等根據(jù)渦扇發(fā)動機工作特點，設(shè)計了一種基于壓氣機出口總壓的新判喘方法；吳虎等基于氣動耦合原理，發(fā)展了某型渦扇發(fā)動機風(fēng)扇出口分流環(huán)計算模型，建立了喘振及其喘振消除的數(shù)值模擬方法。為了更加深入地研究發(fā)動機的喘振特性，鄒龍政、屈霽云等、唐狄毅等分別采用進(jìn)氣口安裝擾流板、引氣口接入高壓氣源和燃油階躍等多種方式開展了逼喘試驗與數(shù)據(jù)分析；曠桂蘭等開展了某型渦軸發(fā)動機整機逼喘試驗，獲得了該型發(fā)動機的喘振邊界線；江勇等采用擾流插板方式研究了2種不同轉(zhuǎn)速下發(fā)動機進(jìn)氣道進(jìn)口截面的喘振壓力波。但是，基于工程應(yīng)用的判喘參數(shù)優(yōu)化工作開展相對較少。張紹基等、劉世官等和郭貴喜等分別從發(fā)動機判喘總體試驗設(shè)計、試驗詳細(xì)過程以及基于虛擬儀器的仿真系統(tǒng)三大方面闡述了某型發(fā)動機的判喘試驗與優(yōu)化設(shè)計，但未對所有相關(guān)參數(shù)開展優(yōu)化。

本文針對某型發(fā)動機在喘振過程中存在漏判的問題，設(shè)計了用于參數(shù)優(yōu)化的喘振信號機載處理仿真系統(tǒng)，開展了對判喘系統(tǒng)所有判喘參數(shù)的仿真分析與驗證。

1 判喘系統(tǒng)原理

某型渦扇發(fā)動機的判喘系統(tǒng)由高壓壓氣機出口總靜壓引氣管路、喘振壓差傳感器、信號處理電路和現(xiàn)場可編程門陣列（Field Programmable Gate Array,FPGA）處理算法組成，處理邏輯如圖1所示。通過安裝于高壓壓氣機出口末端的引壓裝置將該截面的總壓和靜壓氣源引入喘振壓差傳感器的兩端，傳感器輸出總靜壓差信號。根據(jù)高壓壓氣機出口的總靜壓差信號的脈動幅值與頻率可以有效評估發(fā)動機是否發(fā)生喘振。

圖1 判喘系統(tǒng)處理邏輯

機載信號處理電路對傳感器的信號經(jīng)全波整流后，進(jìn)行低通和帶通濾波處理，得到壓差信號的波動量與均值。其中，波動量反映了壓差信號瞬時的變化，當(dāng)發(fā)動機處于穩(wěn)態(tài)工作時，約為0；反映了壓差信號在一段時間內(nèi)的平均值，當(dāng)發(fā)動機處于穩(wěn)態(tài)工作時，等于當(dāng)前的壓差值。在FPGA端算法處理（運行速度為31.25 kHz）中，將與的比值與閾值進(jìn)行比較，當(dāng)＞且持續(xù)一定時間，判斷發(fā)動機發(fā)生喘振。具體算法如下

式中：為壓差信號直流分量高選后結(jié)果；為壓差信號交流分量計算結(jié)果；為壓差信號直流分量高選值；（）為當(dāng)前＞的持續(xù)時間；（+1）為FPGA端運行下一個周期時，＞的持續(xù)時間；（）為當(dāng)前消喘指令持續(xù)時間。

式（1）的作用是為了對進(jìn)行高選，防止發(fā)動機在慢車及以下狀態(tài)時壓差信號過小而導(dǎo)致誤判為喘振，具體分析見第4.2節(jié)。

式（2）的作用是當(dāng)壓差信號由穩(wěn)態(tài)瞬時增大時，為0；當(dāng)壓差信號由穩(wěn)態(tài)瞬時變小時，取為負(fù)數(shù)的的絕對值來衡量壓差信號波動的劇烈程度。由于當(dāng)發(fā)動機發(fā)生喘振時壓差信號均出現(xiàn)先減小的過程，式（2）可以避免當(dāng)發(fā)動機狀態(tài)正常變化時，壓差信號增大過程被誤認(rèn)為是由喘振引起的，從而降低判喘系統(tǒng)誤判的可能性。

式（3）的作用是記錄發(fā)動機喘振狀態(tài)的持續(xù)時間，其中（0）=0，（）≥0。當(dāng)/≥（即壓差信號瞬態(tài)波動量除以平均值大于一定閾值）時，=4；當(dāng)/＜時，=-1。假設(shè)/≥持續(xù)的時間=20 ms時，認(rèn)為發(fā)動機發(fā)生喘振。由于程序在FPGA中每運行1次的時間為1/31250 s，當(dāng)（）≥0.02×31250×4時，認(rèn)為發(fā)動機發(fā)生喘振。

式（4）的作用是確定發(fā)動機執(zhí)行消喘動作的持續(xù)時間，其中（0）=0，（）≥0。當(dāng)確定發(fā)動機發(fā)生喘振后，令=14開始進(jìn)行累加計算。當(dāng)（）≥×14×31250時，發(fā)出消喘指令。在執(zhí)行消喘指令后，當(dāng)（）＜×31250×4時，即判喘系統(tǒng)認(rèn)為發(fā)動機不在喘振狀態(tài)時，令=-1，對（）進(jìn)行減法，當(dāng)（）＜×14×31250時，停止發(fā)出消喘指令。

其中，低通和帶通的濾波參數(shù)、的高選值、判喘閾值、喘振確認(rèn)持續(xù)時間以及消喘指令持續(xù)時間對判喘的有效性、實時性和可靠性均起到不同的作用。不同發(fā)動機的喘振信號特征存在一定的差異，根據(jù)實際發(fā)動機的喘振信號對上述相關(guān)參數(shù)進(jìn)行優(yōu)化，可以明顯提高判喘系統(tǒng)的性能，從而有效地保障發(fā)動機工作的安全性。

2 機載處理仿真系統(tǒng)

為了實現(xiàn)對判喘系統(tǒng)的優(yōu)化，機載處理仿真系統(tǒng)包括了高頻多通道同步信號采集與記錄、機載處理數(shù)字化實現(xiàn)以及喘振信號模擬仿真3個模塊，上位機及系統(tǒng)架構(gòu)分別如圖2、3所示。

圖2 喘振機載處理仿真系統(tǒng)上位機

圖3 仿真系統(tǒng)架構(gòu)

2.1 信號采集與記錄

發(fā)動機臺架的試車數(shù)據(jù)采集與記錄存在多平臺采集、數(shù)據(jù)保存時間不同步、未采集喘振壓差傳感器原始信號等問題，不利于對喘振信號特征進(jìn)行分析與仿真。本文借助PXI計算機的高性能采集板卡，采用高速大數(shù)據(jù)量專用的tdms格式保存數(shù)據(jù)，實現(xiàn)采樣率為100 kHz的6通道同步信號采集與記錄，采集信號包括喘振壓差傳感器原始信號、高壓壓氣機出口總壓信號和高壓壓氣機轉(zhuǎn)速信號等。

本文系統(tǒng)的信號采集與記錄模塊的優(yōu)勢有以下2點：

（1）數(shù)據(jù)采集的采樣率，國內(nèi)其他試驗低于10 kHz，本文達(dá)到100 kHz。在對壓差信號數(shù)據(jù)進(jìn)行離線分析時發(fā)現(xiàn)發(fā)動機在發(fā)生喘振的瞬間存在2~3個周期100 Hz左右的短時波動。另外，對喘振過程前后的數(shù)據(jù)進(jìn)行頻率分析時發(fā)現(xiàn)存在1 kHz以上的高頻信號。目前暫不確定這些高頻成分是否與喘振相關(guān)。但是高采樣率為后續(xù)進(jìn)行更深入的喘振特征分析提供了可能性，特別是對預(yù)測喘振等復(fù)雜算法的研究。

（2）強調(diào)了數(shù)據(jù)采集的同步性要求。信號采集通過2張板卡PXI 4462和PXI 6281完成，2張板卡之間通過共用采樣時鐘保證同步性。對高壓壓氣機出口總壓、靜壓和壓差等信號同步采集，對研究管路對信號的相位延遲等問題有著重要作用。另外，系統(tǒng)可以增加對發(fā)動機其他位置的壓力測量，以了解喘振的始發(fā)位置。

2.2 機載處理數(shù)字化模塊

在PXI計算中，基于LabVIEW圖形化編程語言采用NI PXI-7852R板卡，實現(xiàn)了機載判喘系統(tǒng)的數(shù)字化。通過相同喘振壓差原始信號的回放模擬，對數(shù)字化模塊與實際發(fā)動機控制器的處理模塊進(jìn)行了測試對比，驗證了二者在信號處理結(jié)果和判喘時間上具有較好的一致性。

該模塊采用上位機+FPGA端的架構(gòu)，如圖4所示。其中，上位機主要功能包括波形顯示、參數(shù)配置和數(shù)據(jù)記錄等；FPGA端實現(xiàn)了全波整流、低通濾波、帶通濾波和其他邏輯算法的功能，即實現(xiàn)了喘振壓差信號的完整處理過程。模塊對低通和帶通的濾波參數(shù)、的高選值、判喘閾值以及喘振確認(rèn)持續(xù)時間等判喘參數(shù)設(shè)定了對應(yīng)的變量，實現(xiàn)了上位機對FPGA端參數(shù)的在線修改功能。該功能可以方便試驗人員對相關(guān)參數(shù)進(jìn)行快速修正，極大地縮短了參數(shù)優(yōu)化與驗證所需要的時間，更加深刻地認(rèn)識相關(guān)參數(shù)對判喘的可靠性和實效性的影響。

圖4 機載處理數(shù)字化程序架構(gòu)

2.3 喘振信號模擬仿真

要得到1組性能優(yōu)良的發(fā)動機消喘系統(tǒng)，必須通過整機逼喘試驗來深入研究和確定消喘控制器的對應(yīng)參數(shù)。俄羅斯研制1套發(fā)動機消喘系統(tǒng)通常要進(jìn)行500余次的整機逼喘試驗來驗證，不僅需要巨額的科研經(jīng)費，而且項目研制周期也很長。為了提高研發(fā)效率、降低試驗成本，采用NI PXI-4461板卡設(shè)計了喘振信號模擬仿真模塊。

該模塊的信號模擬輸出包括真實臺架逼喘試驗的喘振壓差信號和根據(jù)喘振特征設(shè)計的定制波形2類。真實的壓差信號模擬可實現(xiàn)對判喘參數(shù)優(yōu)化設(shè)計的快速驗證，且避免了重復(fù)試驗的資源耗費。定制波形可通過軟件配置實現(xiàn)不同頻率和幅值以及噪聲疊加的信號模擬。定制波形的模擬輸出能夠更加深入全面地驗證判喘系統(tǒng)的特性，還可用于實際控制器生產(chǎn)交付時判喘系統(tǒng)的測試驗收。另外，由于高性能板卡的使用，信號模擬輸出的采樣率可達(dá)到100 kHz，輸出誤差小于±10 mV。

3 逼喘試驗與數(shù)據(jù)統(tǒng)計

本文判喘參數(shù)的優(yōu)化設(shè)計基于某型發(fā)動機的逼喘試驗開展。在試驗前，對發(fā)動機的試車臺架進(jìn)行了適應(yīng)性改造，主要包括：在喘振壓差傳感器的引氣端口處增加絕壓傳感器，分別測量高壓壓氣機出口總壓和靜壓；增加喘振壓差傳感器信號輸出，保證在不影響機載控制器正常工作的情況下進(jìn)行壓差原始信號采集；增加各類傳感器信號輸出和臺架消喘指令、機載消喘指令輸出的線纜，并與機載處理仿真系統(tǒng)連接。

在試驗過程中，采用進(jìn)口安裝插板的方法逼喘。在發(fā)動機進(jìn)口安裝可移動的扇形插板，通過專用的地面設(shè)備控制，采用步進(jìn)的方式推進(jìn)插板。當(dāng)發(fā)動機進(jìn)入喘振狀態(tài)時，通過地面設(shè)備能夠快速推出插板，避免發(fā)動機的喘振狀態(tài)持續(xù)時間過長。

逼喘試驗分別在發(fā)動機的小狀態(tài)和大狀態(tài)下開展，相應(yīng)的高頻同步采集數(shù)據(jù)如圖5所示。

圖5 發(fā)動機喘振數(shù)據(jù)

4 參數(shù)分析與優(yōu)化設(shè)計

通過機載處理仿真系統(tǒng)，利用真實數(shù)據(jù)和定制數(shù)據(jù)分別對所有判喘系統(tǒng)的參數(shù)進(jìn)行仿真分析，評估不同參數(shù)對判喘準(zhǔn)確性和時效性的影響，總結(jié)出1套完整的參數(shù)優(yōu)化方法。該方法包括低通和帶通的濾波參數(shù)、的高選值、判喘閾值、喘振確認(rèn)持續(xù)時間以及消喘指令持續(xù)時間等6個參數(shù)的影響分析與優(yōu)化設(shè)計。

4.1 濾波參數(shù)

當(dāng)喘振引起的壓力脈動頻率不在或部分不在帶通濾波范圍內(nèi)時，會由于喘振造成的壓力脈動的測量幅值衰減而導(dǎo)致判喘失??；當(dāng)?shù)屯V波參數(shù)過大或過小時，會導(dǎo)致無法相對準(zhǔn)確地得到壓力的平均值，影響判喘的有效性。因此，要確定機載判喘的低通和帶通的濾波參數(shù)，需要對發(fā)動機喘振前后分別進(jìn)行快速傅里葉變換（Fast Fourier Tranform，F(xiàn)FT），了解核心機和整機在喘振過程中的頻率成分。

對于低通濾波參數(shù)，分別設(shè)定截止頻率為0.50、0.72、1.00、2.00 Hz，濾波效果如圖6所示。從圖中可見，小于1.00 Hz的截止頻率得到的可以較好地反映壓差的平均值；在0.72 Hz時反映平均值的同時，延遲較小。結(jié)合實際電路設(shè)計，保持原有低通截止頻率為0.72 Hz不變。

圖6 喘振信號不同低通濾波參數(shù)濾波效果

對于帶通濾波參數(shù)，發(fā)動機在喘振過程中的壓差信號的脈動頻率的變化范圍主要為1~10 Hz，如圖7所示。因此，將原帶通濾波范圍3~166 Hz更改為1~166 Hz。

圖7 發(fā)動機大狀態(tài)喘振壓差數(shù)據(jù)FFT結(jié)果

4.2 DPCKDC的高選值DPCKDC0

判喘系統(tǒng)中設(shè)定高選值的目的，是為了避免發(fā)動機在慢車及以下狀態(tài)時由于壓差信號較小，發(fā)動機狀態(tài)的正常波動可能引起判喘系統(tǒng)的誤判。這是因為與的比值作為衡量發(fā)動機波動狀態(tài)的值，當(dāng)分母值較小時，分子小范圍變化也會引起比值發(fā)生較大的變化。原則上，的高選值不能太小，否則存在喘振誤判的風(fēng)險；同時，高選值不能太大，否則可能存在對喘振無法判斷或延遲判斷的可能性。

對多臺份多批次發(fā)動機在正常工作過程中的值進(jìn)行統(tǒng)計分析發(fā)現(xiàn)，在慢車過程中值在1.13±0.04范圍內(nèi)變化。參考前期設(shè)計值，保持為原有設(shè)定值，即為0.983。

4.3 判喘閾值A(chǔ)0

判喘閾值是最直接參與判斷發(fā)動機是否發(fā)生喘振的參數(shù)，是判喘系統(tǒng)相關(guān)參數(shù)中最為重要的設(shè)定值。在完成濾波參數(shù)和的高選值設(shè)定后，將所有實際發(fā)動機逼喘數(shù)據(jù)進(jìn)行模擬回放，通過仿真系統(tǒng)的信號處理數(shù)字化模塊查看實際波動量與均值的比值（值）。在綜合考慮漏報和虛警的情況下，設(shè)定合理的判斷閾值。

邢家瑞和秦海波等提出喘振閾值根據(jù)轉(zhuǎn)速關(guān)系分段設(shè)定。當(dāng)高壓轉(zhuǎn)子轉(zhuǎn)速＜85%時，判喘閾值設(shè)為0.3；當(dāng)高壓轉(zhuǎn)子轉(zhuǎn)速≥85%時，判喘閾值設(shè)為0.6。但在本文發(fā)動機不同狀態(tài)下的喘振試驗中，值變化均小于0.6。根據(jù)該型發(fā)動機的喘振試驗的實際情況，本文的判喘閾值設(shè)置為單一值，不根據(jù)轉(zhuǎn)速進(jìn)行分段。

值仿真結(jié)果如圖8所示。從圖中可見，值在發(fā)動機正常工作過程中為0.30以下，在喘振時最大值為0.45。因此，為了減小系統(tǒng)誤判的風(fēng)險，將判喘閾值由0.30更改為0.35。

圖8 A值仿真結(jié)果

4.4 喘振確認(rèn)持續(xù)時間T0

從實際逼喘數(shù)據(jù)中可見，在喘振過程中壓差信號的波動類似于正弦波，設(shè)定合理的喘振確認(rèn)持續(xù)時間可以避免由于信號干擾和噪聲導(dǎo)致的虛警。而該值的設(shè)定是在確定上述判喘參數(shù)后，經(jīng)過數(shù)據(jù)模擬回放與仿真試驗得到的。

通過對發(fā)動機在小狀態(tài)和大狀態(tài)下的逼喘試驗數(shù)據(jù)的仿真可知，當(dāng)判喘閾值為0.35時，喘振確認(rèn)持續(xù)時間分別需要至少18.3 ms和23.5 ms。為了提高判喘可靠性和實效性，將喘振確認(rèn)持續(xù)時間由20 ms更改為15 ms。

4.5 消喘指令持續(xù)時間T1

設(shè)定消喘指令持續(xù)時間的目的，是為了確定當(dāng)判喘成功并且發(fā)動機退出喘振后，延遲多久退出消喘動作。這樣的操作可以最大程度保證發(fā)動機完全退出喘振狀態(tài)，同時又不會過度執(zhí)行消喘。該值對判喘影響不大，可根據(jù)實際發(fā)動機狀態(tài)設(shè)定。本文在優(yōu)化設(shè)計中保持原值不變。

5 參數(shù)驗證

在參數(shù)優(yōu)化前，判喘系統(tǒng)在發(fā)動機不同狀態(tài)下的逼喘試驗中，存在小狀態(tài)消喘指令發(fā)出延遲時間過長和大狀態(tài)漏判的問題。這樣的結(jié)果無法滿足發(fā)動機正常工作的控制要求，必須保證在全狀態(tài)下判喘的有效性和實時性。

在參數(shù)優(yōu)化后，依次進(jìn)行了控制器在回路的仿真試驗驗證和實際發(fā)動機的逼喘試驗驗證。試驗結(jié)果表明，本文的參數(shù)優(yōu)化有效地縮短了消喘指令的發(fā)出時間，將小狀態(tài)的判喘時間由86 ms縮短到53 ms，發(fā)動機小狀態(tài)判喘優(yōu)化仿真結(jié)果如圖9所示；同時，在發(fā)動機大狀態(tài)的逼喘試驗中，成功判斷出發(fā)動機喘振，參數(shù)優(yōu)化后發(fā)動機大狀態(tài)判喘結(jié)果如圖10所示。其中，當(dāng)消喘指令值為1時，表示發(fā)出消喘指令；當(dāng)消喘指令值為0時，表示不發(fā)出消喘指令。

圖9 發(fā)動機小狀態(tài)判喘優(yōu)化仿真結(jié)果

圖10 參數(shù)優(yōu)化后發(fā)動機大狀態(tài)判喘結(jié)果

6 結(jié)論

（1）在仿真系統(tǒng)應(yīng)用的基礎(chǔ)上，完整地了解了相關(guān)參數(shù)對判喘的影響和作用，形成了完整的參數(shù)優(yōu)化設(shè)計流程；

（2）在喘振過程中的高壓壓氣機總壓與靜壓的壓差信號脈動頻率變化范圍主要為1~10 Hz，將原帶通濾波范圍3~166 Hz更改為1~166 Hz；

（3）通過對喘振數(shù)據(jù)的仿真與分析，考慮判喘可靠性和實效性，將喘振確認(rèn)持續(xù)時間由20 ms更改為15 ms;

（4）經(jīng)參數(shù)優(yōu)化后的判喘系統(tǒng)，將小狀態(tài)的判喘時間由86 ms縮短到53 ms，并完成了大狀態(tài)的喘振判喘。

本文通過對相關(guān)參數(shù)的優(yōu)化設(shè)計，有效地提高了某型發(fā)動機判喘系統(tǒng)的準(zhǔn)確性和時效性，驗證了喘振信號處理仿真系統(tǒng)的有效性，為其他型號航空發(fā)動機的判喘系統(tǒng)設(shè)計與優(yōu)化提供了技術(shù)支持。