王 強(qiáng),吳 偉,劉 東,婁華語,王良模
(1.航空工業(yè)沈陽飛機(jī)設(shè)計研究所, 沈陽 110000;2.南京理工大學(xué) 機(jī)械工程學(xué)院, 南京 210014)
艙門是民機(jī)的運(yùn)動功能部件,其功能、壽命、安全性、維修性和可靠性對民機(jī)安全有直接影響[1]。因此,針對飛機(jī)艙門收放系統(tǒng)的故障診斷十分重要。目前,對飛機(jī)液壓系統(tǒng)的故障檢測工作主要依靠實(shí)驗(yàn)平臺進(jìn)行實(shí)物實(shí)驗(yàn)。但是,單純依靠實(shí)驗(yàn)獲取故障數(shù)據(jù)集存在工作繁瑣、數(shù)據(jù)量不足等問題。因此,大量研究人員采用仿真方法對液壓系統(tǒng)進(jìn)行故障研究。王明錦等[2]結(jié)合 AMESim+Simulink平臺,仿真分析了閥控馬達(dá)液壓系統(tǒng)控制效果與工作特性。Liu[3]分析了民機(jī)主起落架艙門伸縮原理和故障原因,并深入研究了上位鎖(Uplock)的控制過程和邏輯。Zhu等[4]基于AMESim和LMSlian對飛機(jī)艙門鎖機(jī)構(gòu)進(jìn)行聯(lián)合建模仿真,并通過向模型中注入故障信息為故障診斷和分析提供了一種獲取故障數(shù)據(jù)的方法。
在故障診斷方法研究方面,儲瑄等[5]采用故障樹系統(tǒng)分析理論對飛機(jī)某艙門收放系統(tǒng)進(jìn)行診斷分析,并進(jìn)行了可靠性試驗(yàn)。近年來,機(jī)器學(xué)習(xí)和神經(jīng)網(wǎng)絡(luò)逐漸發(fā)展成系統(tǒng)故障診斷的主流方向。王帥星等[6]在建模基礎(chǔ)上通過故障注入獲取內(nèi)泄漏故障數(shù)據(jù),并采用粒子群(PSO)優(yōu)化BP神經(jīng)網(wǎng)絡(luò)進(jìn)行故障診斷研究。Liu等[7]利用變分模式分解(VMD)、多重分形去趨勢波動分析(MFDFA)和概率神經(jīng)網(wǎng)絡(luò)(PNN)對機(jī)艙門執(zhí)行機(jī)構(gòu)進(jìn)行故障診斷,有效解決了由非線性和非平穩(wěn)執(zhí)行機(jī)構(gòu)振動信號引起的故障診斷困難。Cui等[8]在DBN模型的結(jié)構(gòu)基礎(chǔ)上,提出利用SSA算法對優(yōu)化DBN模型隱藏層的節(jié)點(diǎn)數(shù)量,并建立SSA-DBN網(wǎng)絡(luò)模型用于飛機(jī)艙門液壓系統(tǒng)的故障診斷。Yang等[9]在飛機(jī)艙門執(zhí)行機(jī)構(gòu)故障診斷方面,采用稀疏自動編碼器(SAE)算法避免簡單的直接篩選過程,采用雙向長短期內(nèi)存RNN算法(BiLSTM-RNN)應(yīng)用于故障診斷,減少了傳感器數(shù)量。
綜上,民機(jī)艙門收放系統(tǒng)故障模擬存在故障數(shù)據(jù)少、故障診斷精度低的問題,因此,提出基于CPSO-BP神經(jīng)網(wǎng)絡(luò)的飛機(jī)艙門收放系統(tǒng)故障診斷模型。首先,根據(jù)飛機(jī)艙門系統(tǒng)工作特性及高發(fā)故障,分析確定系統(tǒng)典型故障模式;其次,建立飛機(jī)艙門收放系統(tǒng)仿真模型,通過典型故障的仿真分析獲取故障數(shù)據(jù)集;最后,分別采用BP神經(jīng)網(wǎng)絡(luò)、混沌粒子群算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)進(jìn)行故障診斷。
飛機(jī)艙門收放系統(tǒng)驅(qū)動方式主要分直線式和旋轉(zhuǎn)式。由于旋轉(zhuǎn)式驅(qū)動方式占用空間小、艙門開閉時間短,已成為先進(jìn)飛機(jī)的首選方案。
圖1所示的旋轉(zhuǎn)式飛機(jī)艙門收放系統(tǒng)采用閥控液壓馬達(dá)驅(qū)動飛機(jī)艙門完成開閉運(yùn)動。液壓油從液壓源出發(fā)經(jīng)過電液比例閥,經(jīng)由流量控制閥驅(qū)動控制液壓馬達(dá)旋轉(zhuǎn),最后經(jīng)減速器實(shí)現(xiàn)艙門開閉。控制時采用經(jīng)典PID控制。
圖1 旋轉(zhuǎn)式飛機(jī)艙門收放系統(tǒng)原理示意圖
飛機(jī)艙門收放系統(tǒng)正常開閉時,通過外場統(tǒng)計艙門收放系統(tǒng)故障現(xiàn)象主要為艙門不動作、艙門動作不到位和艙門動作緩慢。根據(jù)外場高發(fā)故障及飛機(jī)艙門系統(tǒng)工作特性,確定4類典型故障模式:
1) 液壓馬達(dá)外泄漏:液壓馬達(dá)隨著使用時間的增加,管嘴保護(hù)圈容易從管嘴端面與殼體端面間隙中擠出,出現(xiàn)外部油液滲漏現(xiàn)象[10-11]。隨著泄漏量增大,導(dǎo)致液壓馬達(dá)進(jìn)出口兩端壓差降低,進(jìn)而導(dǎo)致艙門動作緩慢甚至動作不到位。
2) 液壓油污染:液壓油混入空氣后,其彈性系數(shù)降低,變得更易被壓縮,會導(dǎo)致作動元件出現(xiàn)爬行現(xiàn)象等問題,造成的故障占液壓系統(tǒng)故障的75%以上[12]。
3) 節(jié)流閥阻塞:固體顆粒污染液壓油最常見且危害最嚴(yán)重,主要包括硬質(zhì)顆粒和軟質(zhì)顆粒2種污染物。這些物質(zhì)極易堵塞如節(jié)流閥等元件的閥口,造成元件的閥口流通面積減小,對系統(tǒng)產(chǎn)生負(fù)面影響。
4) 流量控制閥閥芯磨損:流量控制閥是閥控系統(tǒng)控制信號直接作用元件,頻繁變化閥口開度會導(dǎo)致閥芯與閥體產(chǎn)生磨損,液壓油液污染物更是加劇了這一現(xiàn)象[13],造成流量控制閥內(nèi)泄漏故障,降低控制效果。
依照艙門收放系統(tǒng)結(jié)構(gòu),在AMESim中選取相應(yīng)元器件,建立旋轉(zhuǎn)式飛機(jī)艙門驅(qū)動系統(tǒng)仿真模型如圖2所示。系統(tǒng)采用閥控液壓馬達(dá)驅(qū)動艙門運(yùn)動,控制中采用位置PID和速度PID相結(jié)合的方法,通過位置閉環(huán)和速度閉環(huán)保證到位精度和艙門運(yùn)行速度控制效果。為簡化結(jié)構(gòu)將部分控制信號用超級元件“MOVE_SPEED”、“TIP_SPEED”、“OPEN_DECELERATION_ ANGLE”、“CLOSE_DECELERATION_ANGLE”及“LOCATION”代替。
設(shè)置各部件相關(guān)參數(shù),仿真時間為5 s,時間間隔為0.01 s。運(yùn)行艙門模型,結(jié)束后查看系統(tǒng)中的各種參數(shù)。飛機(jī)艙門的位移和速度曲線如圖3所示。
圖3 艙門位移和速度曲線
從圖3中可以看出,飛機(jī)艙門收放系統(tǒng)在仿真周期內(nèi)完成了打開和關(guān)閉兩階段動作,所用時間分別為1.84 s和1.87 s,低于設(shè)計許用要求的2 s;且艙門在開閉末端有明顯減速動作,末端速度保持在15(°)/s以內(nèi),滿足系統(tǒng)設(shè)計要求。即搭建的艙門收放系統(tǒng)模型可用。
如圖2所示,根據(jù)飛機(jī)艙門收放系統(tǒng)工作特性及典型故障,進(jìn)行故障注入:① 液壓馬達(dá)進(jìn)出口兩端并聯(lián)節(jié)流閥(見圖中1處),并調(diào)節(jié)等效孔徑大小實(shí)現(xiàn)故障模擬;② 修改液壓油空氣含量參數(shù)實(shí)現(xiàn)液壓油污染故障注入(見圖中2處);③ 修改節(jié)流閥等效孔徑實(shí)現(xiàn)節(jié)流閥阻塞故障注入(見圖中3處);④ 流量控制閥進(jìn)出口兩端并聯(lián)節(jié)流閥(見圖中4處),實(shí)現(xiàn)流量控制閥閥芯磨損故障注入。
基于建立的飛機(jī)艙門收放系統(tǒng)故障仿真模型,進(jìn)行飛機(jī)艙門收放系統(tǒng)正常工作及液壓馬達(dá)泄漏、液壓油污染、節(jié)流閥阻塞和流量控制閥磨損等典型故障的仿真分析。飛機(jī)艙門收放系統(tǒng)仿真參數(shù)見表1所示,在正常值和故障值區(qū)間內(nèi)分別等距取10個和20個數(shù)值,其他參數(shù)取正常值,每種故障分別進(jìn)行30次仿真,共獲取120組仿真數(shù)據(jù)。原始數(shù)據(jù)存在重復(fù)等問題,故需要對原始數(shù)據(jù)進(jìn)行預(yù)處理:對重復(fù)數(shù)據(jù)手動進(jìn)行直接刪除。
表1 飛機(jī)艙門收放系統(tǒng)仿真參數(shù)
設(shè)置仿真時間為5 s,采樣頻率為100 Hz,提取A點(diǎn)流量和壓力信號、B點(diǎn)位置信號形成一組故障信號,每組信號共3×500個數(shù)據(jù)點(diǎn)。從正常值和故障值區(qū)間隨機(jī)各抽取一組數(shù)據(jù)構(gòu)成測試集,其余構(gòu)成訓(xùn)練集。
圖4所示為不同狀態(tài)下的流量、壓力和位移變化曲線。相比于正常工況,隨著系統(tǒng)的故障注入,艙門到位速度變慢,系統(tǒng)消耗流量增加,系統(tǒng)壓力降低。此外,壓力、流量和艙門位移在不同故障模式下變化。因此,如果系統(tǒng)中的鎖定機(jī)構(gòu)發(fā)生故障,可以通過系統(tǒng)的工作壓力、流量和位移信號來推斷故障原因。
圖4 飛機(jī)艙門收放系統(tǒng)故障仿真結(jié)果
為解決故障數(shù)據(jù)較少穩(wěn)態(tài),在仿真獲取故障數(shù)據(jù)的基礎(chǔ)上,采用滑動窗口方法對故障信號進(jìn)行重疊采樣[14-15]。為保持重構(gòu)信號原有信號的特征,考慮艙門開閉最長時間為2.5 s,設(shè)置窗口長度為256,步長為1。經(jīng)滑動窗口重疊采樣后,每組故障數(shù)據(jù)可獲得245個樣本,隨后按4∶1劃分訓(xùn)練集和測試集。
采集壓力、流量、速度和位移多種傳感器信號組成故障信號,在數(shù)據(jù)增強(qiáng)基礎(chǔ)上,采用橫向拼接的時序信號特征融合[16]方法,將原有格式為3×256的樣本數(shù)據(jù)重構(gòu)為長度為768的一維信號樣本。橫向拼接特征融合原理如圖5所示。橫向拼接是將特征直接進(jìn)行拼接,將其從頭到尾串接起來,雖然這樣會增加模型的運(yùn)算量和運(yùn)算時間,但能保留最完備的信息,不會造成信息的損失。
對數(shù)據(jù)進(jìn)行抽稀與歸一化操作。最終建立故障數(shù)據(jù)集,包含測試集樣本5 880個,訓(xùn)練集樣本23 535個,每個樣本格式為長度256、數(shù)據(jù)分布在[0,1]間的一維信號。
BP神經(jīng)網(wǎng)絡(luò)是一種按誤差逆向傳播數(shù)據(jù)訓(xùn)練的多層前饋網(wǎng)絡(luò),目前其相關(guān)應(yīng)用最為廣泛。將飛機(jī)艙門收放系統(tǒng)故障數(shù)據(jù)集作為神經(jīng)網(wǎng)絡(luò)輸入,電磁閥磨損、液壓馬達(dá)外泄漏、液壓油混入空氣和節(jié)流閥阻塞作為輸出,建立神經(jīng)網(wǎng)絡(luò)。該BP神經(jīng)網(wǎng)絡(luò)為三層典型網(wǎng)絡(luò)結(jié)構(gòu),含1個隱層,隱層神經(jīng)元個數(shù)設(shè)置為8,最大訓(xùn)練次數(shù)1 000,學(xué)習(xí)率為0.01,目標(biāo)為訓(xùn)練均方差1×10-6。進(jìn)行5次診斷,最終的故障診斷結(jié)果如表2所示??梢钥闯?使用傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)進(jìn)行故障診斷時,平均診斷正確率可達(dá)85.36%,多次診斷結(jié)果的魯棒性較強(qiáng),但診斷結(jié)果仍不太理想。
表2 BP神經(jīng)網(wǎng)絡(luò)綜合故障診斷結(jié)果
BP神經(jīng)網(wǎng)絡(luò)綜合故障診斷結(jié)果不太理想的原因是網(wǎng)絡(luò)訓(xùn)練時BP網(wǎng)絡(luò)中的初始權(quán)重和閾值為隨機(jī)選取產(chǎn)生,容易出現(xiàn)局部收斂極小點(diǎn),會降低擬合效果,故需要對神經(jīng)網(wǎng)絡(luò)初始權(quán)重和閾值進(jìn)行優(yōu)化。本文選取混沌粒子群算法(CPSO)優(yōu)化BP神經(jīng)網(wǎng)絡(luò)。
3.2.1混沌粒子群(CPSO)優(yōu)化算法
粒子群優(yōu)化算法是一種采用速度-位置搜索模型的群體智能尋優(yōu)算法,具有編碼簡單、收斂速度快的特點(diǎn)[17]。粒子群算法中更新粒子位置和速度的公式如下:
vid=wvid+c1r1(piBest-xid)+
c2r2(pgBest-xid)
(1)
xid=xid+vid
(2)
其中:xid表示為粒子位置;vid表示粒子速度;w為慣性權(quán)重;c1和c2為學(xué)習(xí)因子;r1和r2為[0,1]間的隨機(jī)數(shù);piBest和pgBest分別為個體極值和群體極值。
基本粒子群算法尋優(yōu)搜索過程有2個不足:① 隨機(jī)初始化和進(jìn)化過程使得個體極值和群體極值的更新有一定盲目性,影響收斂速度;② 利用式(1)和式(2)更新粒子位置,本質(zhì)是采用正向反饋方式,通過當(dāng)前粒子信息、個體極值和群體極值3個信息獲取下一步迭代位置,容易陷入局部最優(yōu)解。
針對粒子群優(yōu)化算法的不足,提出基于混沌優(yōu)化思想的混沌粒子群優(yōu)化(CPSO)算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)。一般將混沌定義為由確定性方程得到且具有隨機(jī)性運(yùn)動的運(yùn)動狀態(tài)?;煦鐒恿W(xué)中比較典型的混沌模型是Logistic映射[18],其表達(dá)式為
zn+1=μzn(1+zn)
(3)
式中:μ為控制參數(shù),一般取值為4。
針對粒子群算法后期收斂速度慢、易陷入局部最優(yōu)點(diǎn)而出現(xiàn)早熟的問題,通過粒子群的方差來判斷是否早熟,方差表達(dá)式為
(4)
式中:Fi為粒子的適應(yīng)度;Favg為粒子的平均適應(yīng)度。
當(dāng)方差小于給定值時判定粒子進(jìn)入早熟狀態(tài),此時需要重新設(shè)定粒子速度和位置以跳出局部最優(yōu)。
采用混沌理論改進(jìn)粒子群算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)故障診斷流程見圖6。
3.2.2故障診斷結(jié)果的對比分析
為驗(yàn)證CPSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的故障診斷性能,基于已建立的故障數(shù)據(jù)集,分別采用CPSO-BP和BP兩種不同方法對比進(jìn)行故障診斷。CPSO算法種群大小設(shè)置為20,迭代尋優(yōu)100次,其他初始化超參數(shù)如表3所示。
圖6 CPSO-BP神經(jīng)網(wǎng)絡(luò)流程
表3 CPSO算法超參數(shù)
圖7所示為CPSO算法尋優(yōu)適應(yīng)度進(jìn)化曲線。分析可知,CPSO-BP算法在迭代尋優(yōu)44次后,適應(yīng)度不再減小,即尋優(yōu)獲取到最優(yōu)目標(biāo)值。
圖7 CPSO-BP適應(yīng)度進(jìn)化曲線
各神經(jīng)網(wǎng)絡(luò)對測試集故障診斷結(jié)果的混淆矩陣元素如圖8所示。其中0—4分別代表正常工況、液壓馬達(dá)泄漏、液壓油污染、節(jié)流閥阻塞及流量控制閥磨損五類情況。可以清楚得出:相比BP神經(jīng)網(wǎng)絡(luò),CPSO-BP神經(jīng)網(wǎng)絡(luò)針對液壓馬達(dá)泄漏、液壓油污染、節(jié)流閥及流量控制閥磨損四類故障的診斷正確率均有明顯提升,其中液壓油污染和流量控制閥磨損兩類故障診斷的正確率分別提高了27%和34%。
圖8 測試集故障診斷混淆矩陣元素示意圖
此外,在CPSO-BP算法基礎(chǔ)上,研究神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點(diǎn)個數(shù)對CPSO-BP及BP診斷結(jié)果的影響。隨著網(wǎng)絡(luò)節(jié)點(diǎn)個數(shù)的增加,故障診斷模型的診斷精度和診斷時間如圖9所示。可知優(yōu)化后,8節(jié)點(diǎn)時,CPSO-BP診斷精度較其他節(jié)點(diǎn)更高,正確率達(dá)到93%,且診斷時間相比7、9和10節(jié)點(diǎn)更低,綜合計算效率更高。
圖9 不同節(jié)點(diǎn)故障診斷正確率與計算時間
1) 基于AMESim搭建飛機(jī)艙門收放系統(tǒng)故障模型以代替實(shí)際設(shè)備進(jìn)行故障診斷研究,彌補(bǔ)了閥控液壓馬達(dá)驅(qū)動艙門系統(tǒng)建模研究較少的不足,具有仿真精度高、故障數(shù)據(jù)采集方便和數(shù)據(jù)量豐富等優(yōu)點(diǎn);
2) 采用滑動窗口方式進(jìn)行故障數(shù)據(jù)增強(qiáng),極大豐富了故障數(shù)據(jù)集,多傳感數(shù)據(jù)橫向拼接的特征融合方式提高了神經(jīng)網(wǎng)絡(luò)訓(xùn)練與診斷精度;
3) 采用BP神經(jīng)網(wǎng)絡(luò)進(jìn)行故障初步診斷的平均診斷正確率僅為85.37%。為提升診斷正確率,采用混沌粒子群算法(CPSO)優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的初始權(quán)重和閾值,CPSO-BP神經(jīng)網(wǎng)絡(luò)的系統(tǒng)故障診斷正確率可以達(dá)到93%,提高了液壓油污染和流量控制閥磨損兩類故障的診斷正確率,能較準(zhǔn)確識別故障類型,可用于飛機(jī)艙門系統(tǒng)故障診斷。