民航發動機拆發日期預測研究

彭鴻博　劉孟萌　王悅閣

(中國民航大學航空工程學院,天津　300300)

?

民航發動機拆發日期預測研究

彭鴻博劉孟萌王悅閣

(中國民航大學航空工程學院,天津300300)

針對民航發動機拆發送修業務的需求，采用Daubechines的3階小波對民航發動機歷史性能數據進行降噪平滑處理，建立民航發動機性能衰退模型。結合影響民航發動機拆發日期的時壽件的軟時間、適航指令/服務通告限制等因素，對民航在翼發動機進行拆發日期預測。試驗結果驗證了該方法在預測民航在翼發動機拆發日期的可行性，可以滿足民航發動機拆發送修業務的需求。

民航發動機時壽件拆發日期預測狀態監控小波變換排氣溫度裕度AD/SB

0　引言

民航發動機是飛機的“心臟”，其健康狀況與飛行安全息息相關。對民航發動機采取視情維修的策略。當發動機性能衰退到較差的狀況時，就要拆下發動機送修，以保證其安全性和經濟性[1]。送修時，如果過早拆下發動機，就會造成發動機資源的浪費，增加航空公司運營成本；過晚拆下，就會對飛行安全構成威脅。這兩種情況都是不允許的[2]。民航發動機性能衰退的主要表征是排氣溫度裕度(exhausted gas temperature margin，EGTM)不斷減小。

小波變換在時域、頻域均具有良好的局部性質[3]，本文利用小波變換對民航發動機EGTM的歷史數據進行處理并分析，得出民航發動機的EGTM衰退模型，并結合發動機時壽件(life limited parts，LLP)的軟時間、適航指令/服務通告(airworthiness dictate/service bulletin，AD/SB)限制等因素，預測民航在翼發動機的拆發日期。

1　發動機拆發日期影響因素分析

1.1性能衰退的影響

性能衰退情況是指隨著發動機使用循環次數的增加，發動機的EGTM、燃油流量、壓氣機轉子轉速等氣路參數出現衰退的趨勢[4]。本文采用EGTM作為預測發動機拆發日期的性能參數。在民用發動機投入使用后，隨著發動機使用時間的增加，EGTM呈現一定規律的衰退趨勢。當EGTM接近警戒值時，就必須拆下發動機送修，以恢復發動機性能。

1.2時壽件的影響

時壽件是指民航發動機結構中有明確使用壽命限制的零部件，如發動機中的風扇盤、風扇軸、低壓渦輪軸、中壓渦輪盤等。由于受材料、結構及工作環境等因素的影響，必須在時壽件使用到規定的時間前，將發動機拆下送修，進行時壽件的修理或更換[5]。時壽件是影響送修目標的硬性條件。利用式(1)計算該臺發動機的剩余最低循環數，再折算成下發日期。

C1=min(CLLP,1-C0.1,CLLP,2-C0.2,…，

CLLP,i-C0,i)

(1)

式中：CLLP,i為該臺發動機第i個時壽件的設計循環數；C0,i為第i個時壽件使用的循環數。

1.3AD/SB的影響

AD/SB的影響是指適航部門和制造商規定發動機使用一定循環數后，必須拆下發動機進行某些工作。例如，某發動機規定在運行14 035次循環后，應拆下發動機進行繞線管的孔探檢查[6]。所以，可以通過這個循環數來計算發動機的下發日期。

下發日期的計算公式為：

C2=min(CAD/SB,1-C0,CAD/SB,2-C0,…，

CAD/SB,i-C0)

(2)

式中：CAD/SB,i為第i個AD/SB中規定的循環數；C0為已經使用的循環數。

2　發動機拆發日期的預測

2.1發動機拆發日期預測流程

發動機從性能衰退、時壽件和AD/SB限制等角度，均可以得到一個拆發日期。取最靠前的那個日期就是該臺發動機的最終拆發日期。綜上所述，民航發動機拆發日期的預測流程如圖1所示。

2.2粗大誤差剔除

圖2　去除粗大誤差前后的EGTM樣本圖

2.3小波變換去噪平滑

為了精確地提取EGTM的衰退趨勢，利用小波變換，對去除粗大誤差后的樣本數據作進一步的降噪平滑。小波變換是在傳統的傅里葉變換的基礎上發展演變而來的，因其具有良好的時頻局部特性，得到了廣泛的應用。一般而言，數據源f(n)都是包含噪聲的，我們可以用式(3)表示一個數據源的噪聲模型，e(n)為高斯白噪聲，小波變換就是要抑制e(n)，以恢復f(n)，從而達到去除噪聲的目的[9]。

s(n)=f(n)+σe(n)

(3)

式中：e(n)為噪聲；σ為噪聲強度。

EGTM采樣信號是一維信號，運用小波變換進行降噪平滑的具體步驟為：

①對去除粗大誤差的EGTM采樣信號進行小波分解；

②設定各層細節的閾值，對分解的小波進行過濾；

③將得到的EGTM信號趨勢進行小波重構。

本文采用Daubechines的3階濾波器對EGTM衰退的趨勢進行提取，小波分析趨勢對比如圖3所示。

圖3　小波分析趨勢對比圖

2.4EGTM衰退趨勢線的擬合

本文采用一元線性回歸的相關方法，求解民航發動機EGTM的衰退率。一元線性回歸模型可以表示為[10]：

yi=β0+β1xi+εi

(4)

式中：yi為EGTM估計值；β0為EGTM初始衰退值；β1為EGTM衰退率；xi為飛機飛行循環數；εi為隨機誤差。

直接運用Matlab中的Curve Fitting Tool對式(4)求解，得到的EGTM衰退趨勢如圖4所示。

圖4　EGTM衰退趨勢線

將式(4)求得的結果代入一元線性回歸模型，得到EGTM衰退模型，即：

y=-0.017 5x+114.3

(5)

為保證衰退模型完全可靠，排除小波降噪平滑后出現的異常值、周期性因素等干擾，運用殘差圖診斷回歸效果?！皻埐顖D”是指以飛行循環次數為橫坐標，以殘差為縱坐標，將每一次飛行循環的殘差描繪在該平面坐標上所形成的圖形。當描繪的點圍繞殘差等于0的直線上下隨機散布，說明回歸直線對原觀測值的擬合情況良好；否則，說明回歸直線對原觀測值的擬合不理想。利用Matlab得到的殘差分布圖如圖5所示。

圖5　殘差分布圖

從圖5可以直觀看出，只有很小一部分點的置信區間不包含零點，說明所得到的一元線性模型還是十分可靠的。EGTM衰退率為0.017 5 ℃/循環。

航空實際運營表明：全壽命周期內發動機的EGTM衰退率并不是一成不變的，而是在裝機初期衰退較快，中期比較平穩，末期衰退加速。根據不同階段的EGTM數據，采用本文的方法，可求解各階段的EGTM衰退率。

3　發動機拆發日期預測應用驗證

對某航空公司一臺CF56-22發動機拆發日期的預測方法進行驗證，編號為67889的發動機已經飛行了28 717次飛行循環。當前EGTM為105.55℃，EGTM的警戒值(當EGTM降為警戒值時，發動機必須拆下送修)為22 ℃。該臺發動機的部分時壽件清單如表1所示。AD/SB規定：該發動機在40 000次循環后要拆下進行SPOOL檢查。

表1　部分時壽件清單

利用本文方法對該案例進行分析。

①從性能衰退角度考慮，當前EGTM為105.55 ℃，警戒值為22 ℃，按照當前的EGTM衰退率0.017 5，求得再經過4 775次飛行循環，EGTM就達到了警戒值。在實際使用過程中，航空公司會對發動機的核心機進行幾次水洗，假設為3次，每次可以獲得10 ℃左右的EGTM恢復；將這些因素考慮在內，再經過4 775+2 057=6 832次飛行循環，發動機的EGTM就達到警戒值，必須拆下送修。

②從時壽件角度考慮，由表1和式(1)可知，壽命件Fan Disk再用1 283次循環，就必須更換。

③從AD/SB限制角度考慮，再飛行11 283次循環，就必須拆下發動機進行SPOOL檢查。

綜上所述，該航空公司編號為67 889的發動機再

使用1 283次循環，就要送修。發動機的日使用率約為4次循環，所以大約再使用320天后必須送修。

4　結束語

發動機拆發日期的預測可以有效減少發動機維修業務中的維修不足或者過度維修。本文結合性能衰退、時壽件到壽、AD/SB限制等參數，提出了一種拆發日期的預測方法。該方法可以直觀、有效地預測民航在翼發動機短、中期的拆發日期，為發動機維修工程師提供指導。

但該方法仍存在一些不足，如EGTM的線性擬合只適用于短、中期預測。發動機拆發日期預測精度的提高，還需收集更多的性能參數，以建立更準確的狀態變化趨勢模型。另外，不能僅局限于本文提出的影響發動機拆發日期的這三個因素，還應根據航空公司的實際情況，考慮其他的影響因素，如故障因素等。因此，該方法還需進一步完善，以得到更精確的預測日期。

[1] 雷達,鐘詩勝.面向MRO的民航發動機拆發日期預測系統[J].計算機集成制造系統,2013,19(7):1715-1720.

[2] 付旭云,鐘詩勝.民用航空發動機維修計劃啟發式算法[J].計算機集成制造系統,2010,16(7): 52-57.

[3] 王佰行.基于小波分析的數據平滑處理算法研究與應用[J].科技創新導報,2009(31):142-142.

[4] 曹惠玲,林達錦,曾學鋒.CFM56-7B型發動機衰退模型建立與研究[J].中國民航大學學報,2010,28(4):9-12.

[5] 王悅.CFM56-7B發動機送修時機選擇[J].航空維修與工程,2012(2):94-96.

[6] 左洪福，蔡景，陳志雄.航空維修工程學[M].北京：科學出版社，2011:211-219.

[7] 郝英,孫健國,白杰.航空燃氣渦輪發動機氣路故障診斷現狀與展望[J].航空動力學報,2003,18(6):753-760.

[8] 彭云飛.航空發動機狀態參數處理技術及其應用研究[D].哈爾濱:哈爾濱工業大學,2009.

[9] 吳偉,蔡培升.基于MATLAB的小波去噪仿真[J].信息與電子工程,2008,6(3):220-222.

[10]隋永志,李書明,倪繼良,等.基于 MATLAB的航空發動機試車數據的線性回歸研究[J].航天制造技術,2009(6): 45-48.

Research on Prediction of the Demolition Date for Repairing of Civil Aviation Engine

In accordance with the demands for maintenance of civil aviation engines,the demolition date for repairing has to be determined.By adopting Daubechines 3 wavelet,the historical performance data of civil aviation engines are smoothed and de-noised to establish the performance deterioration model.In addition,combining with some of the limiting factors that influencing the demolition date for repairing,such as the soft time of life limited parts,and airworthiness directives / service bulletin,etc.,the prediction of the in-service engine is conducted.The practical examples verify the feasibility of this method of predicting the demolition date for repairing; it meets the requirement of business of demolition for repairing of civil aviation engines.

Civil aviation engineLife limited partRemoval datePredictionStatus monitoring Wavelet transform Exhausted gas temperature marginAirworthiness dictate/service bulletin(AD/SB)

彭鴻博(1972—)，男，2007年畢業于法國國立民航大學航空發動機專業，獲碩士學位，副教授；主要從事航空發動機健康管理方向的研究。

TH707;TP274+.2

A

10.16086/j.cnki.issn 1000-0380.201608003

修改稿收到日期：2015-12-30。