應(yīng)用單元體參數(shù)的航空發(fā)動機性能預(yù)警方法

黃燕曉, 郝紅勛, 郭家琛

(1. 中國民航大學(xué) 通用航空學(xué)院， 天津 300300；2. 中國民航大學(xué) 飛行技術(shù)學(xué)院， 天津 300300；3. 南京航空航天大學(xué) 民航學(xué)院， 江蘇 南京 210016)

?

應(yīng)用單元體參數(shù)的航空發(fā)動機性能預(yù)警方法

黃燕曉1, 郝紅勛2, 郭家琛3

(1. 中國民航大學(xué) 通用航空學(xué)院， 天津 300300；2. 中國民航大學(xué) 飛行技術(shù)學(xué)院， 天津 300300；3. 南京航空航天大學(xué) 民航學(xué)院， 江蘇 南京 210016)

基于安全關(guān)口前移的性能預(yù)警要求，提出以單元體參數(shù)的航空發(fā)動機性能預(yù)警方法,提取發(fā)動機原理分析性能參數(shù)相關(guān)的可測參數(shù)與參數(shù)組合.以發(fā)動機機隊為研究對象，運用數(shù)據(jù)統(tǒng)計的σ準(zhǔn)則構(gòu)建性能指數(shù)動態(tài)閾值模型.以PW4077D型機隊5臺發(fā)動機1 440～1 980循環(huán)數(shù)據(jù)驗證模型，并與主成分性能閾值模型比較.結(jié)果表明：構(gòu)建的模型計算結(jié)果與主成分模型結(jié)果一致，即性能動態(tài)閾值隨循環(huán)數(shù)增加而復(fù)雜下降，在1 500,1 560和1 840等循環(huán)附近閾值有所恢復(fù).

航空發(fā)動機； 單元體參數(shù)； 預(yù)警方法； 動態(tài)閾值；σ準(zhǔn)則

民航干線飛機裝配的渦輪風(fēng)扇發(fā)動機，具有推力大、耗油率低特點，因長期工作在高速、高壓和高溫等惡劣環(huán)境下[1-3]，容易引起性能衰退或突變，帶來安全隱患.因此，航空發(fā)動機性能需要可靠監(jiān)控預(yù)警[4-5].傳統(tǒng)的監(jiān)控預(yù)警方法常以氣路參數(shù)或振動參數(shù)構(gòu)建性能超限或趨勢模型，預(yù)測突發(fā)故障或分析性能衰退趨勢等[5-6].然而，以氣路或振動參數(shù)而不是單元體參數(shù)進行性能監(jiān)控預(yù)警，將使單元體結(jié)構(gòu)設(shè)計優(yōu)勢無法充分應(yīng)用，也不能將發(fā)動機安全關(guān)口前移至單元體層面[7]，采用經(jīng)驗確定的靜態(tài)閾值不能滿足性能精確監(jiān)控預(yù)警要求.從單元體層面構(gòu)建性能參數(shù)邊界或閾值建立性能預(yù)警監(jiān)控模型，國內(nèi)外學(xué)者都做了相關(guān)工作[1,8,9].這些研究工作更注重以經(jīng)驗監(jiān)控發(fā)動機瞬態(tài)性能，判斷精度不高且沒有預(yù)警能力.鑒于此，本文提出單元體參數(shù)表征發(fā)動機性能的預(yù)警方法.

1　單元體及表征參數(shù)

民航飛機裝配的渦輪風(fēng)扇發(fā)動機均是以單元體結(jié)構(gòu)設(shè)計為原則，常用機型為PW4077D發(fā)動機[10].目前，使用機型為GE90系列發(fā)動機，文中以二者為例進行比較分析.PW4077D發(fā)動機有風(fēng)扇、核心機、低壓渦輪和附件傳動裝置4個核心單元體，以13個占位表示著不同的發(fā)動機位置，為使表述簡潔，不計各占位的小數(shù)點，如12.5占位以125表述.航空發(fā)動機單元體性能參數(shù)與可測參數(shù)，如表1所示.

表1　航空發(fā)動機單元體性能參數(shù)與可測參數(shù)

2　航空發(fā)動機性能預(yù)警方法

2.1發(fā)動機性能單元體參數(shù)權(quán)值計算方法

(1)

各參數(shù)的相關(guān)系數(shù)αi，j為

(2)

2.2機隊各循環(huán)性能動態(tài)閾值計算方法

以機隊性綜合閾值替代各臺發(fā)動機性能閾值作為準(zhǔn)則，依據(jù)發(fā)動機單元體參數(shù)值及各參數(shù)權(quán)值，提出基于單元體參數(shù)的發(fā)動機性能指數(shù).不同循環(huán)t下的性能指數(shù)計算式為

(3)

3　性能預(yù)警定模型算例驗證

3.1定量計算單元體參數(shù)性能權(quán)值

發(fā)動機性能的主成分計算結(jié)果，如表2所示.

表2　發(fā)動機性能的主成分計算結(jié)果

由表2可知：主成分1,2對應(yīng)特征根大于1，且方差貢獻率分別為59.154 0%和23.118 0%，累計貢獻率為82.272 0%，在80%～90%準(zhǔn)則內(nèi).故主成分1,2能夠包含選取的7個參數(shù)表征的發(fā)動機性能，選取主成分k=2，特征根λk={(4.414 0, 1.168 0)|k=1,2}.計算各主成分對應(yīng)于7個參數(shù)的載荷系數(shù)和線性組合系數(shù)，然后，計算相關(guān)系數(shù)矩陣中發(fā)動機排氣溫度特征值和權(quán)值，所得結(jié)果如表3所示.

表3　發(fā)動機單元體參數(shù)相關(guān)系數(shù)矩陣特征值與權(quán)值

圖1　5臺發(fā)動機性能指數(shù)趨勢Fig.1　performance index trends of five engines

3.2機隊各臺發(fā)動機的性能指數(shù)

由式(3)可計算出機隊各臺發(fā)動機的性能指數(shù).機隊中5臺發(fā)動機性能指數(shù)趨勢，如圖1所示.由圖1可知：PW4077D機隊中5臺發(fā)動機在該循環(huán)區(qū)間內(nèi)性能總體處于降趨勢，但過程變化復(fù)雜.

3.3發(fā)動機機隊性能指數(shù)動態(tài)閾值

圖2　機隊發(fā)動機性能閾值趨勢Fig.2　Fleet performance threshold values trend

由于5臺發(fā)動機在相同循環(huán)的發(fā)動機性能指數(shù)值視作正態(tài)分布，按照σ準(zhǔn)則分析該機隊的性能指數(shù)動態(tài)閾值,如圖2所示.圖2中：顯著性水平α=0.05;n為工作循環(huán)數(shù).

依據(jù)主成分分析法，計算第1,2主成分特征值對應(yīng)的發(fā)動機性能指數(shù).計算機隊各循環(huán)性能的閾值分別為：6.305 0， 1.764 4， 4.531 0，1.339 3， 2.179 9，0.850 1，-0.725 2，-2.579 0，-1.492 8，-3.442 2，5.277 7，-0.200 1，-1.209 6，-1.382 4.由圖2可知：在1 440～1 980循環(huán)區(qū)間的性能指數(shù)閾值變化與定權(quán)值法得到的結(jié)果從趨勢上是一致的，故基于定權(quán)值法和數(shù)據(jù)統(tǒng)計結(jié)合計算出的機隊性能指數(shù)閾值曲線能夠替代發(fā)動機性能基線監(jiān)控發(fā)動機無故障狀態(tài)下的性能預(yù)警監(jiān)控.

4　結(jié)論

1) 結(jié)合航空發(fā)動機單元體結(jié)構(gòu)設(shè)計的主要特點，提出以單元體參數(shù)及組合替代氣路參數(shù)和振動參數(shù)表征發(fā)動機性能[11].

2) 提出以機隊而不是單臺發(fā)動機為預(yù)警計算的基準(zhǔn)，定權(quán)值和主成分法得到的性能閾值趨勢變化均隨工作循環(huán)數(shù)增加而單調(diào)下降.下降過程比較復(fù)雜，在維護維修等工作后性能閾值會顯著恢復(fù)，但不會恢復(fù)到初始循環(huán)水平.

[1]SIMONDL,ARMSTONGJB.Anintegratedapproachforaircraftengineperformanceestimationandfaultdiagnostics[J].JournalofEngineeringforGasTurbinesandPower,2013,135(7):1-10.

[2]HOCHGREBS,DENNISD,AYRANCII,etal.Forcedandself-excitedinstabilitiesfromleanpremixed,liquid-fuelledaero-engineinjectorsathighpressuresandtemperature[C]∥ASMETurboExpTurbineTechnicalConferenceandExposition.SanAntonio:AmericanSocietyofMechanicalEngineers,2013:1-12.

[3]GREITZEREM.SomeaerodynamicproblemsofaircraftEngines:fiftyyearsafter-the2007IGTIscholarlecture[J].JournalofTurboMachinery,2009,131(3):1-13.

[4]孫見忠,左洪福,劉鵬鵬,等.基于動態(tài)線性模型的民航發(fā)動機性能監(jiān)控模型與在翼壽命預(yù)測[J].系統(tǒng)工程理論與實踐,2013,33(2):3243-3250.

[5]VOLPONIAJ.Gasturbineenginehealthmanagement:Past,present,andfuturetrends[J].JournalofEngineeringforGASturbinesandPower,2014,136(5):1-20.

[6]單曉明,宋云峰,黃金泉,等.基于神經(jīng)網(wǎng)絡(luò)和模糊邏輯的航空發(fā)動機狀態(tài)監(jiān)視[J].航空動力學(xué)報,2009,24(10):2356-2361.

[7]趙威,張德志.單元體結(jié)構(gòu)發(fā)動機總體結(jié)構(gòu)尺寸控制與分析[J].航空發(fā)動機,2014,40(5):45-49.

[8]劉志榮,朱睿,梁忠生,等.發(fā)動機健康基線及評估準(zhǔn)則研究[J].廈門大學(xué)學(xué)報(自然科學(xué)版),2010,49(4):520-525.

[9]周俊,劉俊華,楊濤,等.航空發(fā)動機在線振動檢測系統(tǒng)的開發(fā)[J].計算機測量與控制,2015,23(11):3599-3602.

[10]黃燕曉,李書明,王凌云,等.面向核心單元體的航空發(fā)動機性能評估研究[J].制造業(yè)自動化,2015,37(9):86-92,105.

[11]李東民，張慧，李玉善.不同氧空比下發(fā)動機動力性能與 NOx排放仿真分析[J].華僑大學(xué)學(xué)報(自然科學(xué)版),2016,37(3):281-286.

(責(zé)任編輯： 陳志賢 英文審校： 崔長彩)

Research on Aero-Engine Performance Early Alerting Method Using Module Parameters

HUANG Yanxiao1, HAO Hongxun2, GUO Jiachen2

(1. College of General Aviation, Civil Aviation University of China, Tianjin 300300， China;2. College of Flight Technology, Civil Aviation University of China, Tianjin 300300， China;3. College of Civil Aviation, Nanjing University of Aeronautics and Astronautics, Nanjing 210016， China)

To satisfy the requirements caused by the pre crossing safety for aero-engine performance early alerting, the paper put forward to take out the measurable parameters and their combinations related to the performance parameters by engine operating principle with aero-engine performance early alerting method based on module parameters. With the engine fleet as the research objects, then proposed performance index dynamic threshold model using the statistical criteria with the fleet as the research objects. The performance data of 5 engines with 1 440-1 980 cycles in PW4077D fleet were imputed into the model to verify the performance index dynamic threshold model. Then the paper compared the results generated by combining the fixed weighted value and the statistical criteria with that generated by the PCA method. The comparison showed that the results of two methods were in agreement with each other, i.e., when the performance index dynamic threshold values decrease complexly as the cycle number within the 1 440-1 980 cycles increases, and the threshold values of the 1 500, 1 560 and 1 840 cycles in their neighborhood had recovered somewhat.

aero-engine; module parameters; alerting method; dynamic threshold;σcriteria

10.11830/ISSN.1000-5013.201605001

2016-07-23

黃燕曉(1980-)，男，講師，博士，主要從事飛機發(fā)動機健康管理技術(shù)的研究.E-mail:oldsea0592@163.com.

國家安全生產(chǎn)監(jiān)督管理局安全生產(chǎn)重大事故防治關(guān)鍵技術(shù)科技項目(tianjin-0001-2015AQ)

V 239

A

1000-5013(2016)05-0527-04