石向陽
(濱州學院 飛行學院,山東 濱州 256603)
由于航空發(fā)動機在飛機系統(tǒng)中的特殊重要性,維修過程中不允許出現(xiàn)差錯,一旦出現(xiàn)故障就可能導致飛機出現(xiàn)空中停車[1]等危險事故,嚴重影響飛機的正常使用,甚至釀成航空事故,從而造成大的經(jīng)濟損失。因此,航空發(fā)動機的故障診斷和狀態(tài)監(jiān)控格外重要。
近年來,有關神經(jīng)網(wǎng)絡的算法研究成為人工智能領域熱點問題。神經(jīng)元模型來源于生物的大腦,大腦里有非常多的腦細胞,每一個腦細胞可以看作是一個神經(jīng)元,神經(jīng)元之間相互連接構成一個復雜的神經(jīng)網(wǎng)絡[2]。神經(jīng)網(wǎng)絡包含有多種模型,例如基于梯度算法的BP神經(jīng)網(wǎng)絡、Elman神經(jīng)網(wǎng)絡和RBF徑向基神經(jīng)網(wǎng)絡等[3],每種算法的特點不同,適用范圍也有所區(qū)別,其中,BP神經(jīng)網(wǎng)絡應用范圍最廣。BP神經(jīng)網(wǎng)絡屬于一種依據(jù)誤差逆?zhèn)鞑ニ惴ㄓ柧毜奶厥舛鄬忧梆伨W(wǎng)絡,該方法能夠解決一些復雜問題,且能夠對此展開推測、聯(lián)想與記憶等系列環(huán)節(jié)[4],使得故障診斷率得以提高,所以BP神經(jīng)網(wǎng)絡很適合應用于現(xiàn)代大型儀器設備,例如汽車、飛機等的故障診斷。
國內(nèi)外許多學者針對BP神經(jīng)網(wǎng)絡應用展開了研究,例如FENG提出了基于BP神經(jīng)網(wǎng)絡的航空活塞發(fā)動機機油系統(tǒng)故障診斷[5],喬文生等提出了將BP神經(jīng)網(wǎng)絡應用于立磨齒輪箱專家診斷系統(tǒng)[6],牛華提出了將BP神經(jīng)網(wǎng)絡應用于汽車發(fā)動機故障診斷[7],朱濤提出將ABC-BP神經(jīng)網(wǎng)絡應用于航空發(fā)動機故障診斷[8],均為該領域研究奠定了一定基礎。本文采用目前應用最廣泛的BP神經(jīng)網(wǎng)絡算法,針對PW4000航空發(fā)動機氣路故障進行診斷研究。
BP神經(jīng)網(wǎng)絡可以包含大量的輸入與輸出關系,它的網(wǎng)絡結構包含三層:輸入層、隱含層和輸出層[3]。在隱含層的傳輸輸入后,如果樣本輸出的期望值與實際輸出值相差比較大,則轉向反向傳播,系統(tǒng)調(diào)整各層的連接權值以減少誤差,直到滿足精度要求。一個典型的BP神經(jīng)網(wǎng)絡結構如圖1所示[7]。
圖1 BP神經(jīng)網(wǎng)絡圖
因為航空發(fā)動機的數(shù)據(jù)很多,根據(jù)相關的文獻資料,確定出幾個可以反映出發(fā)動機故障類型的數(shù)據(jù)。本文所選取的數(shù)據(jù)來自某航空公司PW4000發(fā)動機的巡航狀態(tài)時的歷史監(jiān)測數(shù)據(jù),一共討論了4種氣路故障,分別是:發(fā)動機高壓壓氣機故障、低壓壓氣機故障、高壓渦輪故障和低壓渦輪故障[9]。為簡便計算,數(shù)據(jù)的采集只有4種,分別是發(fā)動機排氣溫度TEGT(K)、燃油流量mf(kg/s)、高壓轉子轉速N2(rpm)和低壓轉子轉速N1(rpm)。實驗中,每種故障各收集了50組數(shù)據(jù),故障總數(shù)為200組。將40組數(shù)據(jù)用做測試,10組數(shù)據(jù)用來驗證。
線性函數(shù)歸一化在處理數(shù)據(jù)中很常見,通過線性函數(shù)歸一化可以把數(shù)據(jù)放置于單位為一的區(qū)間內(nèi),可以更快更方便的計算數(shù)據(jù)。公式如下:
x=(xi-xmin)/(xmax-xmin),
式中,x表示線性函數(shù)歸一化處理后得到的數(shù)值,xi表示每個發(fā)動機性能參數(shù)第i個樣本值,xmin表示數(shù)據(jù)樣本中的最小值,xmax表示數(shù)據(jù)樣本中的最大值。
因為所處理數(shù)據(jù)太多,本文僅選擇部分典型數(shù)值,以作為線性函數(shù)歸一化數(shù)據(jù),見表1。
表1 部分線性函數(shù)歸一化數(shù)據(jù)
神經(jīng)網(wǎng)絡一般會有很多個輸入變量,每個變量的特征都很難預先判定。其中會有很多無關的輸入變量,這些變量會增加模型運算的復雜程度,使模型的誤差變大。所以,對輸入變量進一步篩選后再入模型成為很重要的一步。本文采用的是平均影響值算法,它用來表示輸入對輸出影響的大小,運算結果的符號代表是正相關還是負相關,絕對值的大小表示產(chǎn)生影響的大小。將訓練的樣本同時增加10%,構成一組新的數(shù)據(jù),把這組數(shù)據(jù)作為模型的輸入向量。在一般的情況下,TEGT、mf、N1、N2這4個性能參數(shù)為主要的監(jiān)測對象。運算的公式如下:
Miv=|Y-Yout|。
式中,Yout表示輸出向量,Y表示輸入向量增加10%后的結果,Miv值表示無量綱。
經(jīng)過計算,平均影響值如表2所示。
表2 故障平均影響值表
由表2可以得出,N2的總和最小,說明N2對網(wǎng)絡的影響程度最小,當要簡化時,可以首先考慮減少N2的數(shù)據(jù)。
將上述確定的4種故障類型,分別用Y1~Y4表示,數(shù)字1表示該部件發(fā)生了故障,數(shù)字0表示該部件完好。故障模式Y=(Y1,Y2,Y3,Y4)為神經(jīng)網(wǎng)絡的輸出量,并將它們的輸出量用阿拉伯數(shù)字1、2、3、4表示,得表3故障模式。綜上,本文確定的故障模式為Y=(TEGT,mf,N1,N2)。
表3 故障模式
根據(jù)建立的故障模型,確定該BP神經(jīng)網(wǎng)絡結構輸入層是4個神經(jīng)元,輸出層是1個神經(jīng)元。由于BP神經(jīng)網(wǎng)絡的自身局限性,無法確定隱含層的神經(jīng)元數(shù)目??紤]到該數(shù)目將會直接影響輸出層結果,甚至影響準確度。因此,本文依據(jù)試驗數(shù)據(jù)與經(jīng)驗總結的公式相結合的方法,公式為
式中,m表示輸入節(jié)點數(shù)目,n表示輸出節(jié)點數(shù)目,a為1~10 之間的整數(shù),l表示取值范疇的區(qū)間(1,12)。因為l的取值為1~12,所以從1至12逐個進行實驗驗證,算出每個節(jié)點數(shù)對應的均方誤差值,均方誤差表達式為
根據(jù)表4所示,當l取值為7時均方誤差的值最小,說明當隱含層為7時,運算的誤差最小。所以確定的網(wǎng)絡結構為(4-7-1)型。
表4 不同隱含層的均方誤差表
根據(jù)上述信息,基本確定網(wǎng)絡結構模型,然后采用BP神經(jīng)網(wǎng)絡進行訓練和測試。按照前面數(shù)據(jù)進行訓練,網(wǎng)絡迭代次數(shù)的精確度定為0.001,上限設為50 000步,只要程序符合其中任一條件,就會終止訓練。該程序運行結果能夠反映網(wǎng)絡的運行情況,結果如圖2所示。圖3為BP網(wǎng)絡訓練時的誤差收斂曲線,變化平穩(wěn),無較大波動,收斂情況良好。
圖2 BP神經(jīng)網(wǎng)絡界面
圖3 BP神經(jīng)網(wǎng)絡誤差收斂曲線
根據(jù)BP神經(jīng)網(wǎng)絡的數(shù)據(jù),判斷是否達到了預期目標。為了證明所設計的BP網(wǎng)絡結構模型準確度,采用航空公司更為成熟的Elman網(wǎng)絡數(shù)據(jù)進行對比驗證。將所采集的PW4000發(fā)動機故障數(shù)據(jù),按照前文4個公式依次計算,然后建立相應的故障數(shù)據(jù)模型,并對數(shù)據(jù)展開訓練,最后通過MATLAB軟件仿真輸出曲線圖和預測絕對誤差對比圖,具體如圖4、圖5所示。
由圖4和圖5可知,BP神經(jīng)網(wǎng)絡和Elman神經(jīng)網(wǎng)絡的實際輸出與期望輸出曲線值很接近,兩者的絕對誤差也都很接近。但是,考慮到Elman神經(jīng)網(wǎng)絡會出現(xiàn)訓練速度慢和容易陷入局部極小點的缺點,對神經(jīng)網(wǎng)絡的訓練較難達到全局最優(yōu)。BP神經(jīng)網(wǎng)絡的誤差不僅保持在10%內(nèi),而且具備速度更快、誤差率更低等優(yōu)點。因此,可以驗證BP神經(jīng)網(wǎng)絡更能夠確定氣路故障的位置。
圖4 BP神經(jīng)網(wǎng)絡輸出曲線圖
圖5 BP網(wǎng)絡和Elman網(wǎng)絡的預測絕對誤差對比圖
本文以航空公司常用的PW4000航空發(fā)動機為研究對象,以該型發(fā)動機的多種氣路故障作為案例,借助MATLAB軟件,對有關故障數(shù)據(jù)進行訓練,通過BP神經(jīng)網(wǎng)絡和Elman神經(jīng)網(wǎng)絡的實際輸出與期望輸出曲線進行比較,驗證了BP神經(jīng)網(wǎng)絡算法診斷速度更快、診斷結果更精確、誤差率更低,更能夠確定氣路故障的位置,并在一定程度上彌補了其他診斷方法的不足。