基于變分模態(tài)分解近似熵和支持向量機的軸承故障診斷方法

岳應(yīng)娟，孫鋼，蔡艷平

(火箭軍工程大學(xué),西安 710025)

文獻[1]通過仿真及試驗證明了變分模態(tài)分解(Variational Mode Decomposition，VMD)對軸承故障特征頻率的提取能力并提出了基于譜峭度和VMD算法的包絡(luò)譜軸承故障診斷方法。但是由于包絡(luò)分析法的固有缺陷，在形成包絡(luò)信號時需依靠經(jīng)驗進行解調(diào)頻帶參數(shù)的選取[2]，故障特征頻率并不與理論值嚴(yán)格地一一對應(yīng)，人為主觀因素對診斷結(jié)果造成很大影響且難以實現(xiàn)自動化診斷。

近似熵(Approximate Entropy，AE)是從衡量時間序列復(fù)雜性的角度來度量信號中產(chǎn)生新模式的概率大小的物理量，其所需數(shù)據(jù)長度短、抗噪能力強，同時適用于隨機信號和確定性信號[3]。軸承發(fā)生不同的故障時，在不同頻帶內(nèi)信號的近似熵值會發(fā)生改變，故可以通過VMD[4]將信號分解成預(yù)設(shè)的幾個頻帶，然后通過計算每個頻帶的近似熵來判斷是否發(fā)生故障以及故障類型。因此，將VMD近似熵與支持向量機(Support Vector Machine，SVM)相結(jié)合，嘗試用于軸承振動信號的處理，以更好的實現(xiàn)軸承故障的識別。

1 故障診斷方法

1.1 近似熵

近似熵是一種新的序列復(fù)雜性的度量方法，其用一個非負(fù)數(shù)表示某時間序列的復(fù)雜性，時間序列越復(fù)雜，所對應(yīng)的近似熵越大，信號越趨于非平穩(wěn)狀態(tài)，包含的頻率成分也更加豐富、系統(tǒng)更加復(fù)雜;而近似熵越低則表示信號越趨于周期性、信號包含的頻譜越窄[5-6]。當(dāng)軸承處于不同的工作狀態(tài)時，振動信號的頻率分布會發(fā)生不同的變化，同樣近似熵在振動信號不同頻帶的分布也會發(fā)生相應(yīng)的變化，以此作為特征量便可實現(xiàn)軸承不同工作狀態(tài)的判別。

已知X(i)為包含N個數(shù)據(jù)點的時間序列{x(1),x(2),…,x(N)}，其近似熵算法如下：

1)按照設(shè)定模式維數(shù)m重構(gòu)相空間X(i)，X(i)=[x(i),x(i+q),…,x(i+m-1)]，其中i=1,2,…,N-m+1。

2)定義矢量距離d[X(i),X(j)]為X(i)與X(j)對應(yīng)元素中的最大差值，即

。(1)

i,j=1,…,N-m+1,i≠j。

(2)

(3)

5)增加模式維數(shù)m=m+1，重復(fù)1)～4)，求得Φm+1(r)。

6)定義近似熵為

ApEn(m,r,N)=Φm(r)-Φm+1(r)。

(4)

近似熵與m，r和N的取值有關(guān)，但對N的依賴程度最小，一般情況下取m=2，r=0～0.25SD， (SD為時間序列標(biāo)準(zhǔn)差)。從近似熵的計算可以看出,近似熵反映的是時間序列的自相似程度。

1.2 支持向量機

支持向量機是基于結(jié)構(gòu)風(fēng)險最小原理和統(tǒng)計學(xué)習(xí)理論的學(xué)習(xí)機器[7]，其將樣本空間映射到高維空間，通過解決特征空間的線性分類問題實現(xiàn)樣本空間中的非線性分類。相比于傳統(tǒng)學(xué)習(xí)機器，SVM能夠適應(yīng)非常有限的學(xué)習(xí)樣本信息，且對數(shù)據(jù)維數(shù)不敏感，具有很好的推廣性[8-9]。

SVM可以看作一個二分類的問題，其核心思想是尋找一個最優(yōu)分類的超平面，該超平面在保證分類精度的同時，應(yīng)當(dāng)使分類間隔(margin)最大，如圖1所示。

圖1 支持向量機原理示意圖

圖1中H為分類線，在H的兩側(cè)各有一類分類樣本，H1和H2分別表示距離分類線H最近的樣本，H1和H2到H的距離相等且與H保持平行。分類線H的方程為w·x+b=0，對其歸一化，使線性可分樣本集(xi,yi),i=1,…,n,x∈Rd,y∈{-1,1} ,且滿足

yi[(w·xi)+b]-1≥0;i=1,…,l。

(5)

此時的分類間隔為2/‖ω‖，若使間隔最大，則‖ω‖應(yīng)當(dāng)最小，易證‖ω‖2/2最小的分類面為最優(yōu)分類面，位于H1和H2上的訓(xùn)練樣本點即為支持向量。使用SVM實現(xiàn)分類的具體步驟為：

1)有M類的分類問題，給定其訓(xùn)練集

T={(x1,y1),(x2,y2),…,(xl,yl)}；

xi∈Rd,i= 1,2,…,l,y∈{1,2,…,M} 。

(6)

2)對于j=1,2,…,M，把其中一類看作正類，其余M-1類看作負(fù)類，用兩類支持向量機求出決策函數(shù)f(x)，即

fj(x)=sgn(gj(x)),

(7)

式中：g(x)為分類平面函數(shù)。

3)判斷輸入x屬于第J類，其中J是g1(x),g2(x),…,gM(x)中最大者的上標(biāo)。

1.3 故障診斷流程

如圖2所示，基于VMD近似熵和支持向量機的故障診斷方法共分以下幾個步驟：1)采集軸承典型故障下的振動信號，用VMD進行初步分析；2)統(tǒng)計各IMF分量的近似熵，并將得到的近似熵向量作為訓(xùn)練樣本對SVM分類器進行訓(xùn)練；3)將在線采集的信號或待診斷故障信號按照上述提取特征參數(shù)的方法構(gòu)造待診斷故障樣本，并將其輸入訓(xùn)練好的分類器，完成在線故障診斷。

圖2 基于VMD近似熵和支持向量機的軸承故障診斷流程

2 試驗分析

以某變速箱軸承為研究對象，變速箱裝置簡圖如圖3所示[1]。軸承型號為KOYO1205，利用電火花在軸承的外圈溝道、內(nèi)圈溝道和鋼球表面各設(shè)置面積約為3 mm2的點蝕，以模擬外圈故障、內(nèi)圈故障及鋼球故障。變速箱運行時的負(fù)載為25 N·m，軸的轉(zhuǎn)速約600 r/min(fr=10 Hz)，通過B&K3560數(shù)據(jù)采集儀拾取變速箱軸承座的垂直振動信號。采集軸承在正常狀態(tài)及3種故障狀態(tài)下各50個振動信號樣本，每個樣本長度為3 000。

圖3 變速箱裝置簡圖

不同故障軸承實測振動信號的時域波形如圖4所示，從圖中很難對具體的故障類型進行判斷，需對時域信號作進一步處理。

對采集到的200個試驗數(shù)據(jù)分別進行VMD處理，為確保頻帶分解適宜，避免模態(tài)混疊現(xiàn)象，有利于計算各個頻段的近似熵，保證實際分解信號的保真度，設(shè)定分解層數(shù)K=4，懲罰參數(shù)α=2 000，偶權(quán)值參數(shù)τ=0.3[1]。計算VMD處理后得到的IMF分量的近似熵，并將其作為特征向量，由于篇幅所限，從4種工況中各選取3個特征向量進行展示，結(jié)果見表1。

圖4 不同工況下的軸承振動信號

表1 不同工況下振動信號IMF分量的近似熵

分別從4種工況中隨機選出25個近似熵向量，共100個組成訓(xùn)練樣本用于訓(xùn)練SVM。采用5重交叉檢驗法，選擇最優(yōu)的RBF核函數(shù)參數(shù)σ=2.5和懲罰系數(shù)C=250。訓(xùn)練完成后對剩余的100個向量進行分類測試，為減小誤差，上述試驗重復(fù)10次后取其平均值，用識別率作為指標(biāo)進行評價。

VMD近似熵與EMD[10]近似熵、LMD[11]近似熵經(jīng)SVM處理后的分類結(jié)果見表2。由表可知：在EMD近似熵、LMD近似熵作為樣本進行分類時，識別率不高，這是由于使用EMD和LMD在進行信號分解時，受偽分量和模態(tài)混疊的影響較大[12-13]，分解后各分量包含的頻段近似熵也會受到影響，從而降低識別率；而VMD能很好地對信號進行頻段上不同尺度表征，其近似熵向量更具分類特性，易于分類器的分類，識別率也更高，更適用于滾動軸承的故障診斷。

表2 不同故障診斷方法的精度

3 結(jié)束語

通過對VMD近似熵及SVM的分析，可得出如下結(jié)論：

1)VMD是一種自適應(yīng)的信號處理方法，可精確地應(yīng)用于非線性、非平穩(wěn)信號的處理過程中，較好的在不同頻帶尺度下表征原信號。

2)近似熵具有很強的表征信號不規(guī)則性和復(fù)雜性的能力，VMD近似熵與SVM相結(jié)合的方法可以成功地捕捉到不同故障類型信號在不同頻段尺度下變化的不規(guī)則性，對軸承的工作狀態(tài)和故障類型進行辨識。

3)與EMD和LMD近似熵相比，VMD近似熵與SVM相結(jié)合的識別率更高，性能更好，可用于軸承故障的自動化診斷。