朱繼安，劉義

(1.重慶機(jī)電職業(yè)技術(shù)大學(xué) 電氣電子學(xué)院， 重慶 400036；2.重慶紅巖內(nèi)燃機(jī)有限責(zé)任公司， 重慶 400712)

0 引言

配氣機(jī)構(gòu)是柴油機(jī)的重要組成部件，其中配氣機(jī)構(gòu)的氣門組件工作環(huán)境最為惡劣，受頻繁沖擊負(fù)荷作用，產(chǎn)生變形；受高溫高壓燃?xì)鉀_刷，產(chǎn)生高的熱應(yīng)力；潤(rùn)滑條件差，是柴油機(jī)的主要故障點(diǎn)之一。柴油機(jī)工作時(shí)，配氣機(jī)構(gòu)的氣門、推桿、挺柱等零件因溫度升高而伸長(zhǎng)，為了保證在熱態(tài)下氣門的密封性，必須在氣門與搖臂等傳動(dòng)件之間留有緩沖配氣機(jī)構(gòu)熱膨脹的冷間隙和防止凸輪直接作用于氣門的熱間隙。一般在常溫下，進(jìn)氣門冷間隙為0.20～0.35 mm，排氣門冷間隙在0.30～0.40 mm。在柴油機(jī)使用中，由于零件的磨損與變形，氣門間隙會(huì)逐漸增大，造成進(jìn)、排氣門遲開(kāi)、早關(guān)，導(dǎo)致進(jìn)、排氣時(shí)間變短，進(jìn)氣不足，排氣不凈，致使柴油機(jī)的動(dòng)力性與經(jīng)濟(jì)性下降，同時(shí)使各零件之間撞擊和磨損加劇，噪聲增大。若氣門間隙過(guò)小，則會(huì)引起氣門密封不嚴(yán)而漏氣，導(dǎo)致柴油機(jī)功率下降，油耗增加，甚至燒壞氣門組件。因此，對(duì)氣門間隙的檢測(cè)和故障診斷引起了大量專家學(xué)者理論研究和實(shí)踐，如費(fèi)紅姿[1]等在柴油機(jī)上模擬氣門間隙異常、氣門漏氣等故障，采用階比跟蹤方法重采樣柴油機(jī)缸蓋振動(dòng)角度域信號(hào)并進(jìn)行EMD經(jīng)驗(yàn)?zāi)B(tài)分解和Hilbert-Huang變換等算法提取反映柴油機(jī)氣門故障的特征參數(shù)，實(shí)現(xiàn)柴油機(jī)氣門間隙故障診斷；江志農(nóng)[2]等利用Teager能量算子提出基于能量算子梯度鄰域的振動(dòng)自適應(yīng)提取柴油機(jī)氣門間隙故障在線診斷方法等，取得了一些成果，但存在模式混淆或干擾造成振動(dòng)時(shí)刻檢測(cè)誤差大的問(wèn)題。本文采集汽缸蓋上的振動(dòng)信號(hào)，避免模式混淆而采用基于變分模態(tài)VMD算法對(duì)振動(dòng)信號(hào)進(jìn)行分解，計(jì)算多模態(tài)分量IMF的相對(duì)頻譜能量矩得到柴油機(jī)氣門間隙故障特征參數(shù)，采用概率神經(jīng)網(wǎng)絡(luò)技術(shù)對(duì)進(jìn)排氣門間隙等進(jìn)行故障診斷方法的研究。

1 柴油機(jī)振動(dòng)信號(hào)采集與處理

1.1 柴油機(jī)振動(dòng)信號(hào)采集

汽缸蓋上的振動(dòng)信號(hào)源非常豐富，但主要成分是進(jìn)排氣門開(kāi)、關(guān)產(chǎn)生的振動(dòng)和燃油燃燒氣壓帶來(lái)的振動(dòng)，振動(dòng)信號(hào)能夠反映柴油機(jī)內(nèi)部零部件的工作狀況。通過(guò)安裝在柴油機(jī)缸蓋上部的壓電加速度傳感器和數(shù)據(jù)采集系統(tǒng)采樣振動(dòng)信號(hào)，得到原始的振動(dòng)離散信號(hào)X(n)。該振動(dòng)信號(hào)很明顯地表現(xiàn)出進(jìn)排氣門開(kāi)/關(guān)和燃油燃燒氣壓帶來(lái)的振動(dòng)[3]。

1.2 原始采樣數(shù)據(jù)的處理

柴油機(jī)缸蓋振動(dòng)信號(hào)是一種非線性、非平穩(wěn)信號(hào)。針對(duì)這類信號(hào)處理，學(xué)者們提出了多種方法，如小波變換、EMD分解、EEMD分解、LMD分解等算法，小波變換存在基函數(shù)和噪音閾值先驗(yàn)選擇問(wèn)題；經(jīng)驗(yàn)?zāi)Ｊ椒纸釫MD方法可根據(jù)信號(hào)本身的時(shí)間尺度特征自適應(yīng)地將原始采樣信號(hào)分解為一系列本征模式分量IMF的線性和，比小波變換有了進(jìn)步，但EMD方法存在模式混淆和端點(diǎn)效應(yīng)等問(wèn)題，為克服模式混淆問(wèn)題，集成經(jīng)驗(yàn)?zāi)Ｊ椒纸釫EMD算法應(yīng)運(yùn)而生。該方法從噪聲利用的角度，對(duì)信號(hào)添加白噪聲來(lái)改善傳統(tǒng)EMD方法的模式混淆問(wèn)題，然而EEMD的應(yīng)用效果依賴于添加的白噪聲大小和總體平均次數(shù)，缺乏可靠通用的噪聲協(xié)助準(zhǔn)則。為避免上述各種信號(hào)分解算法的不足，本文采用變分模態(tài)VMD信號(hào)分解算法。

1.2.1 變分模態(tài)分解VMD算法

變分模態(tài)分解VMD算法是一種完全非遞歸自適應(yīng)信號(hào)分解方法，它可有效避免EMD等算法產(chǎn)生的模態(tài)混淆和端點(diǎn)效應(yīng)問(wèn)題，其整體框架是變分問(wèn)題[4]。根據(jù)預(yù)設(shè)模態(tài)分量個(gè)數(shù)K對(duì)信號(hào)進(jìn)行分解，將原始信號(hào)f(x)分解為K個(gè)中心頻率為ωk模態(tài)函數(shù)uk，輸入信號(hào)被分解為K個(gè)IMF分量。

VMD算法重新定義本征模態(tài)分量IMF為一調(diào)幅調(diào)頻分量：

uk(t)=Ak(t)cos (Ф(t))

(1)

式中：Ф(t)為非遞減函數(shù)，Ak(t)≥0。

(2)

式中：{uk}={u1，u2,…,uK}為分解得到的K個(gè)模態(tài)分量；{ωk}={ω1,ω2,…,ωK}為各模態(tài)分類的頻率中心；δ(t)為脈沖函數(shù)；?t(·)表示函數(shù)對(duì)t求偏導(dǎo)。

L({uk}，{ωk}，λ)=

(3)

(4)

(5)

(6)

式中：n為迭代次數(shù)；g∈[1,K];Γ為更新因子；ε為一正數(shù)，代表精度。

(7)

利用Parseval/Plancheral傅里葉等距變換，解式(4)、(5)、(6)得：

(8)

(9)

(10)

1.2.2 VMD分解多模態(tài)分量IMF的能量矩分析

VMD具有較好的信號(hào)分解能力，能將信號(hào)不同頻帶成分進(jìn)行有效地分離并且隨著柴油機(jī)故障的發(fā)展，振動(dòng)信號(hào)不同分量的能量會(huì)發(fā)生變化，其頻率成分也會(huì)隨之變化。在此，通過(guò)計(jì)算各IMF分量的相對(duì)頻譜能量矩，提取柴油機(jī)缸蓋振動(dòng)信號(hào)特征值，進(jìn)一步對(duì)柴油機(jī)氣門故障進(jìn)行診斷[5]。

頻譜能量矩定義如下：

令缸蓋振動(dòng)信號(hào)x(t)=u0(t),則缸蓋振動(dòng)信號(hào)及其IMF分量的頻譜能量矩Mk為

(11)

式中：ω為頻率變量，Uk(ω)為uk(t)對(duì)應(yīng)的頻域函數(shù)，k=0,1,2，…,K。

進(jìn)一步定義相對(duì)頻譜能量矩mk為：

(12)

式中：X(ω)為x(t)的頻譜函數(shù)，Mmax為式(11)中Mk的最大值。

1.2.3 振動(dòng)信號(hào)的特征提取

相對(duì)頻譜能量矩特征不僅能反映信號(hào)自身及其IMF分量的頻率信息，而且反映了它們的能量特征，刻畫了信號(hào)能量和頻率的綜合特征，全面描述了柴油機(jī)缸蓋振動(dòng)信號(hào)特征[5]。

基于VMD的柴油機(jī)氣門間隙故障診斷特征參數(shù)提取步驟：

1) 采樣不同氣門間隙下的缸蓋振動(dòng)信號(hào)，得到原始振動(dòng)信號(hào)數(shù)據(jù)樣本；

2) 對(duì)采集的振動(dòng)信號(hào)進(jìn)行VMD分解，得到不同間隙下的多模態(tài)分量IMF；

3) 計(jì)算每個(gè)氣門間隙IFM的相對(duì)頻譜能量矩，得到相應(yīng)柴油機(jī)氣門間隙故障診斷特征參數(shù)。

2 柴油機(jī)實(shí)驗(yàn)實(shí)測(cè)信號(hào)分析和故障診斷樣本的獲取

本實(shí)驗(yàn)在4135四缸四沖程柴油機(jī)上進(jìn)行，發(fā)動(dòng)機(jī)轉(zhuǎn)速為1 500 r/min,輸出功率為5 kW。測(cè)試平臺(tái)如圖1所示。

圖1 柴油機(jī)故障診斷實(shí)驗(yàn)測(cè)試平臺(tái)

由于柴油機(jī)做周期性旋轉(zhuǎn)，檢測(cè)凸輪軸旋轉(zhuǎn)的判缸信號(hào)和曲軸轉(zhuǎn)速信號(hào)相結(jié)合作為采樣缸蓋振動(dòng)信號(hào)的同步信號(hào)，使采樣的缸蓋振動(dòng)信號(hào)為整周期內(nèi)振動(dòng)信號(hào)，缸蓋振動(dòng)信號(hào)采樣頻率為50 kHz。實(shí)驗(yàn)工況如表1所示。

表1 實(shí)驗(yàn)間隙工況

在氣門間隙偏小、氣門間隙偏大、氣門間隙大和正常間隙運(yùn)行4種狀態(tài)下進(jìn)行。

實(shí)驗(yàn)按下列步驟進(jìn)行：

1) 同步等時(shí)間間隔采樣缸蓋振動(dòng)信號(hào)；

2) 完成零均值化、剔除異常值、去除趨勢(shì)項(xiàng)[6]；

3) 階比跟蹤等角度樣條插值重采樣得到原始的缸蓋振動(dòng)采樣角域信號(hào)；

4) 接著對(duì)得到的角域信號(hào)完成VMD分解(K=4、α=2 000)[7]，如圖2所示；

圖2 VMD分解IMF分量時(shí)域圖

5) 對(duì)原始振動(dòng)信號(hào)及其4個(gè)IMF分量進(jìn)行能量矩計(jì)算；

6) 得到每種運(yùn)行狀態(tài)下缸蓋振動(dòng)信號(hào)的特征數(shù)據(jù)mk。

每種運(yùn)行狀態(tài)完成振動(dòng)信號(hào)5周期采樣樣本，得到的樣本相對(duì)頻譜能量矩見(jiàn)圖3，其中樣本號(hào)1～5為A工況運(yùn)行的缸蓋振動(dòng)信號(hào)5個(gè)樣本的特征數(shù)據(jù)，6～10為B工況的，11～15為C工況的，16～20為D工況的。如圖3所示，只有m3和m4的相對(duì)頻譜能量矩可以作為特征參數(shù)用于建立概率神經(jīng)網(wǎng)絡(luò)柴油機(jī)氣門間隙故障診斷的訓(xùn)練樣本，如表2所示。

圖3 樣本相對(duì)頻譜能量矩

3 概率神經(jīng)網(wǎng)絡(luò)

概率神經(jīng)網(wǎng)絡(luò)PNN是以Parazen窗口函數(shù)為激活函數(shù)的一種前向傳播網(wǎng)絡(luò)，不需要反向傳播優(yōu)化參數(shù)，PNN結(jié)合貝葉斯決策來(lái)判斷測(cè)試樣本的類別。PNN吸收徑向基神經(jīng)網(wǎng)絡(luò)和經(jīng)典概率密度估計(jì)原理的優(yōu)點(diǎn)，其網(wǎng)絡(luò)結(jié)構(gòu)與神經(jīng)網(wǎng)絡(luò)RBF的相似特別適合用于模式分類場(chǎng)合，故障診斷就是模式分類的一種特例。PNN網(wǎng)絡(luò)結(jié)構(gòu)由輸入層、隱含層、求和層和輸出層等四層組成，輸入層神經(jīng)元數(shù)量為樣本向量維數(shù)，隱含層的神經(jīng)元數(shù)量為訓(xùn)練樣本數(shù)，求和層神經(jīng)元數(shù)量為分類類型數(shù)量，輸出層只有一個(gè)神經(jīng)元。

3.1 概率神經(jīng)網(wǎng)絡(luò)PNN的訓(xùn)練過(guò)程

PNN的訓(xùn)練過(guò)程就是確定訓(xùn)練樣本和優(yōu)化平滑因子σ的過(guò)程。

3.1.1 確定訓(xùn)練樣本

設(shè)某類別Ci有m個(gè)訓(xùn)練樣本,再進(jìn)行以下操作：

① 將監(jiān)測(cè)對(duì)象運(yùn)行于類別Ci狀態(tài)并采集h個(gè)樣本的原始數(shù)據(jù)，h>m；

② 對(duì)采集的h個(gè)樣本數(shù)據(jù)采用VMD分解得到相應(yīng)的IMF分量；

③ 計(jì)算IMF分量相對(duì)頻譜能量矩作為故障診斷的特征參數(shù)并制成采樣樣本；

④ 采用K聚類算法將h個(gè)采樣樣本聚類為m個(gè)訓(xùn)練樣本。

3.1.2 優(yōu)化平滑因子σ

采用粒子群尋優(yōu)算法優(yōu)化平滑因子σ。將已知類別樣本作為優(yōu)化樣本，改變?chǔ)抑?，?yàn)證優(yōu)化樣本的分類正確性，追求最高的分類正確率，得到優(yōu)化的平滑因子σ。

4 柴油機(jī)氣門間隙故障診斷用概率神經(jīng)網(wǎng)絡(luò)的構(gòu)成

如表2所示，將氣門間隙小、氣門間隙偏大、氣門間隙大、氣門間隙正常等4個(gè)運(yùn)行工況下VMD分解的IMF3分量和IMF4分量的相對(duì)頻譜能量矩m3和m4作為柴油機(jī)氣門間隙特征參數(shù)，也作為概率神經(jīng)網(wǎng)絡(luò)的輸入層的2個(gè)輸入神經(jīng)元；8個(gè)訓(xùn)練樣本作為隱含層的8神經(jīng)元，表2中行數(shù)據(jù)作為隱含層的聚類中心Xi向量；求和層采用4個(gè)神經(jīng)元，進(jìn)行歸類平均概率的計(jì)算；輸出層為一個(gè)神經(jīng)元，得出故障判別結(jié)果。即本柴油機(jī)燃油噴射系統(tǒng)故障診斷概率神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為2-8-4-1。

表2 柴油機(jī)氣門間隙故障診斷訓(xùn)練樣本

5 柴油機(jī)氣門間隙故障診斷用概率神經(jīng)網(wǎng)絡(luò)的檢驗(yàn)

柴油機(jī)按轉(zhuǎn)速1 500 r/min和5 kW負(fù)載分別運(yùn)行在氣門間隙小、氣門間隙偏大、氣門間隙大、氣門間隙正常等4種工況下，采用優(yōu)化平滑因子和未優(yōu)化平滑因子兩種概率神經(jīng)網(wǎng)絡(luò)進(jìn)行工況判斷，每種網(wǎng)絡(luò)每種工況連續(xù)4 h，故障診斷100次，記錄運(yùn)行工況誤判的次數(shù)。檢驗(yàn)記錄如表3所示。

表3 故障診斷檢驗(yàn)表

由表3可知，平滑因子優(yōu)化后，工況診斷正確率有明顯提高。

6 結(jié)論

1) 柴油機(jī)氣門間隙故障造成柴油機(jī)動(dòng)力性和經(jīng)濟(jì)性下降，是柴油機(jī)常見(jiàn)主要故障之一，需要及時(shí)對(duì)該類故障進(jìn)行診斷和維護(hù)維修。

2) 實(shí)驗(yàn)數(shù)據(jù)表明采用變分模態(tài)VMD對(duì)柴油機(jī)缸蓋振動(dòng)信號(hào)進(jìn)行分解并計(jì)算多模態(tài)分量IMF相對(duì)頻譜能量矩可有效提取柴油機(jī)氣門間隙故障特征參數(shù)，采用粒子算法對(duì)概率神經(jīng)網(wǎng)絡(luò)平滑因子尋優(yōu)能明顯提高基于概率神經(jīng)網(wǎng)絡(luò)的柴油機(jī)氣門間隙故障診斷正確率。

3) 對(duì)柴油機(jī)缸蓋振動(dòng)信號(hào)采用的變分模態(tài)VMD分解算法是一種完全非遞歸信號(hào)分解算法，它可有效避免EMD等算法產(chǎn)生的模態(tài)混淆和端點(diǎn)效應(yīng)問(wèn)題，得到清晰中心頻率和IMF分量，是比較先進(jìn)的信號(hào)分解算法。

4) 本文提出的基于缸蓋振動(dòng)信號(hào)相對(duì)頻譜能量矩和優(yōu)化概率神經(jīng)網(wǎng)絡(luò)的柴油機(jī)氣門間隙故障診斷方法能夠有效提高故障診斷正確率，是柴油機(jī)故障診斷領(lǐng)域的有益嘗試。本方法僅在4135四缸四沖程柴油機(jī)進(jìn)行了實(shí)驗(yàn)驗(yàn)證，在其他型號(hào)規(guī)格柴油機(jī)的適用泛化性上還需要進(jìn)一步的實(shí)驗(yàn)驗(yàn)證。