基于層次模糊熵和改進(jìn)支持向量機(jī)的軸承診斷方法研究*

李永波， 徐敏強(qiáng)， 趙海洋, 黃文虎

(哈爾濱工業(yè)大學(xué)深空探測基礎(chǔ)研究中心， 黑龍江 哈爾濱 150001)

?

基于層次模糊熵和改進(jìn)支持向量機(jī)的軸承診斷方法研究*

李永波， 徐敏強(qiáng)， 趙海洋, 黃文虎

(哈爾濱工業(yè)大學(xué)深空探測基礎(chǔ)研究中心， 黑龍江 哈爾濱 150001)

提出一種新的軸承故障特征提取方法——層次模糊熵 (Hierarchical Fuzzy Entropy ,HFE)。層次模糊熵包括層次分析和模糊熵計(jì)算。與多尺度模糊熵相比, 層次模糊熵既分析信號的低頻分量又分析信號的高頻分量，因而能提取更全面、準(zhǔn)確的故障信息。改進(jìn)支持向量機(jī)(Improved support vector machine based binary tree，ISVM-BT)相比其他多分類器具有識別率更高的優(yōu)勢，因此提出了一種基于層次模糊熵和改進(jìn)支持向量機(jī)的軸承故障診斷方法。首先將HFE作為故障特征提取工具，然后將所得的特征向量輸入到改進(jìn)支持向量機(jī)進(jìn)行模式識別。通過軸承故障診斷的工程應(yīng)用，表明該方法可以有效提取軸承故障特征，實(shí)現(xiàn)軸承不同故障類型和故障程度的準(zhǔn)確識別。

故障診斷； 層次模糊熵(HFE)； 改進(jìn)支持向量機(jī)(ISVM-BT)； 滾動軸承

引　言

滾動軸承作為一種基礎(chǔ)的機(jī)械構(gòu)件，在工業(yè)領(lǐng)域中發(fā)揮著極其重要的作用。人們對滾動軸承的故障診斷進(jìn)行了大量的研究[1-3]。對于軸承的故障診斷通常采用振動分析方法，通過振動信號采集、故障特征提取和模式識別進(jìn)行軸承故障診斷。針對滾動軸承故障信號的非平穩(wěn)、非線性特征，如何有效提取軸承故障信息是最重要的環(huán)節(jié)。近年來，已有多種非線性分析方法如小波分析、神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)、經(jīng)驗(yàn)?zāi)B(tài)分解、分形維數(shù)等。

Pincus[4-5]最早提出近似熵的概念，通過比較數(shù)據(jù)和其自身來反映時(shí)間序列的復(fù)雜度。Richman[6]在近似熵的基礎(chǔ)上提出了樣本熵，克服了近似熵自身匹配的缺點(diǎn)。由于近似熵和樣本熵都采用階躍函數(shù)定義相似性，這與實(shí)際樣本類邊緣比較模糊不符。為克服上述不足，陳偉婷[7]提出了模糊熵(Fuzzy Entropy)的概念。然而近似熵、樣本熵和模糊熵只能從單一尺度上度量時(shí)間序列復(fù)雜性，與復(fù)雜度之間沒有直接對應(yīng)關(guān)系，Costa[8-9]等提出了多尺度熵概念，首先對時(shí)間序列進(jìn)行空間尺度分割，然后計(jì)算每個(gè)尺度下樣本的熵值，進(jìn)而得到時(shí)間序列在不同尺度下的復(fù)雜度。鄭近德[10]將粗粒分割與模糊熵相結(jié)合提出多尺度模糊熵，并將其應(yīng)用于軸承不同類型故障的診斷中。

文獻(xiàn)[11-12]指出多尺度熵只能分析時(shí)間序列的低頻部分，忽略了高頻成分，進(jìn)而提出層次熵，它既可反映信號高頻分量復(fù)雜度又可反映低頻分量復(fù)雜度。本文在模糊熵的基礎(chǔ)上，借鑒層次分析的思想，提出了一種新的熵值計(jì)算方法——層次模糊熵。由于測得的軸承振動加速度信號的故障信息頻帶分布比較豐富，只考慮信號中的低頻成分不能完全提取隱藏故障特征。為了能夠更加全面分析非線性、非平穩(wěn)信號在各個(gè)頻段的特征，本文將層次模糊熵作為故障特征提取的工具應(yīng)用到滾動軸承故障診斷中。

二叉樹的層次結(jié)構(gòu)在二叉樹支持向量機(jī)(SVM-BT)中發(fā)揮著重要作用。目前二叉樹的層次結(jié)構(gòu)的設(shè)計(jì)方法主要包括：類間樣本距離法[13]和類內(nèi)樣本分布法[14]。其中類間樣本距離是先分出距離最遠(yuǎn)的類，而類內(nèi)樣本分布法是先分出分布最廣的類?；趦煞N方法的優(yōu)點(diǎn)，設(shè)計(jì)一種新的可分性測度，稱為ISVM-BT，可有效提高識別率[15]。

本文利用層次模糊熵提取故障特征，結(jié)合ISVM-BT提出了一種新的軸承診斷方法。在此基礎(chǔ)上，將此方法與多尺度模糊熵進(jìn)行比較。實(shí)驗(yàn)結(jié)果表明，該方法可有效完成軸承不同故障類型和程度的準(zhǔn)確識別診斷。

1　多尺度模糊熵和層次模糊熵

1.1多尺度模糊熵

樣本熵和模糊熵的定義參考文獻(xiàn)[7,16], 為了更準(zhǔn)確地描述時(shí)間序列的復(fù)雜性，Costa提出多尺度概念[8]。鄭近德借鑒多尺度分析與模糊熵的概念提出了多尺度模糊熵(multi-scale fuzzy entropy，MFE)并將其應(yīng)用到軸承故障診斷中[10]。MFE方法的計(jì)算流程如下：

(1)

式中τ為尺度因子，τ=1,2,…，n(n為正整數(shù))。

(2) 對于每一個(gè)得到的粗粒向量求取模糊熵，得到n個(gè)粗粒向量的模糊熵值，將其表示成尺度因子τ的函數(shù)，稱為多尺度模糊熵分析。

(2)

式中相似容限r(nóng)=(0.1～0.25)*SD，SD為原始序列的標(biāo)準(zhǔn)差。

從上述多尺度模糊熵的計(jì)算流程可以看出，多尺度模糊熵是從不同尺度因子上度量時(shí)間序列的復(fù)雜度。實(shí)際上，在多尺度模糊熵中，粗粒化的過程是求平均化的過程，即只考慮了原始序列的低頻分量，忽略高頻分量。對于故障信息主要分布較豐富的時(shí)間序列，多尺度模糊熵不能滿足要求。為了提取信號中高頻分量的故障信息，本文引入了層次模糊熵的概念，相比于多尺度模糊熵，層次模糊熵同時(shí)考慮了信號中低頻分量和高頻分量，從而能提供更加全面、準(zhǔn)確的時(shí)間模式信息。

1.2層次模糊熵

江英提出了層次熵的概念，用于度量時(shí)間序列在不同節(jié)點(diǎn)處的復(fù)雜性，并將其成功應(yīng)用到生物學(xué)信號分析中[11]。借鑒層次熵中層次分割的優(yōu)勢，結(jié)合模糊熵的定義，提出了層次模糊熵(hierarchical fuzzy entropy，HFE)的概念。與層次熵類似，層次模糊熵的計(jì)算流程如下：

(1) 給定長度為N的時(shí)間序列{u(i),i=1,2,…，N}，定義平均數(shù)算子Q0和Q1如下：

(3)

(4)

注：N=2n，n是正整數(shù)。算子Q0和算子Q1的長度為2n-1。根據(jù)算子Q0和Q1原始序列可重構(gòu)為

(5)

當(dāng)j=0或j=1,定義矩陣Qj算子如下

(6)

(2) 構(gòu)造一個(gè)n維向量[γ1,γ2,…,γn]∈{0,1}, 則整數(shù)e可表示為

(7)

式中正整數(shù)e對應(yīng)的向量為[γ1,γ2,…,γn]。

(3) 基于向量[γ1,γ2,…,γn]，定義時(shí)間序列u(i)每一層分解的節(jié)點(diǎn)分量如下

(8)

式中k表示層次分割中的第k層，如圖1所示，原始時(shí)間序列u(i)在第k+1層的低頻和高頻部分分別用uk,0,uk,1表示。

(4) 對所得每一個(gè)層次分量求其模糊熵，得到2k個(gè)層次分量的模糊熵值，即為層次模糊熵分析，層次模糊熵可表示為：

HFE(u,k,e,m,r)=FuzzyEn(uk,e,m,r)

(9)

圖1　時(shí)間序列u(i)的層次分割示意圖(k=3)Fig.1　Hierarchical tree diagram of time series u(i) with k=3

綜上所述，Q0和Q1算子是低頻部分和高頻部分，與Haar小波的低通和高通濾波的原理相一致[11]。在圖1中最左側(cè)的分解節(jié)點(diǎn)U1, 0,U2, 0和U3, 0的模糊熵值分別對應(yīng)多尺度分析中尺度所得的模糊熵值，即分解節(jié)點(diǎn)Uk, 0對應(yīng)多尺度分析中尺度的模糊熵值。以上分析說明多尺度模糊熵只分析時(shí)間序列低頻部分的模式信息(層次模糊熵的最左側(cè)的分解節(jié)點(diǎn))，忽略了高頻部分的模式信息，而層次模糊熵在分析時(shí)間序列高頻部分的同時(shí)還分析低頻部分的模式信息。實(shí)際測得軸承振動信號在其高頻部分也存在重要的故障信息，信號低頻部分的信息并不能完全反映軸承故障的本質(zhì)特征，說明了對信號進(jìn)行層次分析的重要性。

1.3參數(shù)選擇

根據(jù)層次模糊熵的定義，在進(jìn)行層次模糊熵計(jì)算前需要設(shè)置4個(gè)參數(shù)：嵌入維數(shù)m，模糊函數(shù)的相似容限r(nóng)，模糊函數(shù)的邊界梯度n和層次分解的層數(shù)k。(1)嵌入維數(shù)m的取值會影響動態(tài)重構(gòu)時(shí)信息量的多少。m取值過小會造成重構(gòu)時(shí)信息的丟失，m取值過大會有很多的詳細(xì)信息，但m越大，計(jì)算所需要的數(shù)據(jù)長度N就越大(N=10m～30m)，因此綜合考慮設(shè)置m=2。(2)相似容限r(nóng)的取值對統(tǒng)計(jì)特性有影響，r過大會造成統(tǒng)計(jì)信息的丟失；r過小估計(jì)出的統(tǒng)計(jì)特性不準(zhǔn)確。因此一般取r為0.1～0.25倍原始數(shù)據(jù)的標(biāo)準(zhǔn)差，即r=(0.1～0.25) SD。(3)模糊函數(shù)的邊界梯度n在向量的相似性計(jì)算中發(fā)揮權(quán)重作用，n>1時(shí)，更多地計(jì)入較近的向量對其相似度的貢獻(xiàn)，而更少地計(jì)入較遠(yuǎn)的向量的相似度貢獻(xiàn)；n過大會導(dǎo)致細(xì)節(jié)信息的喪失，而n<1時(shí)則相反。為了獲取盡量多的細(xì)節(jié)信息，根據(jù)文獻(xiàn)[7]，取較小的整數(shù)值n=2。(4)分解層數(shù)k，k值過大會影響計(jì)算效率并且會導(dǎo)致參與每一個(gè)層次分量計(jì)算的點(diǎn)減少，同時(shí)k值過小會導(dǎo)致原始序列頻帶劃分不夠詳細(xì)，從而不能獲得足夠的從低頻到高頻的層次分量，綜合考慮取k=3。

1.4仿真數(shù)據(jù)試驗(yàn)

為了比較多尺度模糊熵和層次模糊熵，根據(jù)參考文獻(xiàn)[11]定義一個(gè)白噪聲信號，數(shù)據(jù)長度N=5000 如圖2(a)所示，對其進(jìn)行傅里葉變換和歸一化處理，其結(jié)果如圖2(b)所示。由圖2可知，白噪聲信號在各個(gè)頻帶上的復(fù)雜度近乎不變。

然后用層數(shù)k=3的HFE和尺度τ=8的MFE分別求解白噪聲信號的模糊熵值，同時(shí)為了觀察HFE和MFE算法的穩(wěn)定性，在每一個(gè)尺度上分別求取100組獨(dú)立白噪聲的熵值的誤差棒如圖3所示。

圖2　白噪聲的時(shí)域波形和頻譜圖Fig.2　Waveform and FT( Fourier transform) spectrum of white noise

圖3　白噪聲的層次模糊熵和多尺度模糊熵 Fig.3　HFE and MFE curves of white noise

由圖3可以看出，對于白噪聲信號，由HFE計(jì)算得出的層次熵值隨分解節(jié)點(diǎn)的變化平穩(wěn)，熵值變化近似獨(dú)立于節(jié)點(diǎn)的變化, 此現(xiàn)象與不同頻帶上白噪聲的熵值近似不變的結(jié)論保持一致[11]。同時(shí)觀察MFE計(jì)算的白噪聲的模糊熵值變化，發(fā)現(xiàn)模糊熵值隨著尺度增加單調(diào)下降，這與白噪聲在各個(gè)頻帶復(fù)雜度近似不變結(jié)論不符。由于MFE只分析了仿真信號的低頻成分，忽略其高頻的成分，而HFE能同時(shí)分析信號中的低頻和高頻成分，因而能提取更加豐富的信息。 最后觀察HFE和MFE的計(jì)算100組獨(dú)立白噪聲的熵值的誤差棒，通過圖3發(fā)現(xiàn)，尺度較小時(shí)HFE和MFE都具有很好穩(wěn)定性，但隨著尺度增大，HFE相比于MFE的誤差棒更小，這說明HFE相比于MFE的穩(wěn)定更好。

2　改進(jìn)二叉樹支持向量機(jī)(ISVM-BT)

二叉樹支持向量機(jī)(ISVM-BT)識別能力主要依賴于二叉樹層次結(jié)構(gòu)，目前的二叉樹生成算法主要有：類間樣本距離[13]和類內(nèi)樣本分布[14]。兩種算法從不同的角度描述了樣本的可分性，且具有較好的推廣能力。文獻(xiàn)[15]在結(jié)合兩種算法優(yōu)勢的基礎(chǔ)上提出了一種新的生成二叉樹層次結(jié)構(gòu)的算法，理論上能夠?qū)㈩愰g樣本距離較大且類內(nèi)樣本分布越廣的類優(yōu)先分離出來。

為實(shí)現(xiàn)上述目的，采用歐氏距離來分別描述類間樣本距離和類內(nèi)樣本分布的情況。給定同類特征樣本集合{xi，i=1，2，…，k}，類內(nèi)樣本距離的定義如下：

類內(nèi)樣本歐氏距離為

(10)

樣本xi到其他類內(nèi)樣本的平均距離

(11)

定義類內(nèi)所有樣本到其他類內(nèi)樣本的平均距離

(12)

同理定義類間樣本距離，首先給定兩類樣本集合{xi，i=1，2，…，k}和{yj，j=1，2，…，k}，其中xi∈類A，yj∈類B，定義類A與類B的類間樣本平均距離如下：

類A與類B間的歐氏距離

(13)

計(jì)算類A樣本xi對到類B各樣本間的平均距離

(14)

定義類A和類B間的平均距離

(15)

將類內(nèi)平均距離和類間平均距離結(jié)合起來定義類A和類B的可分性IA,B為

(16)

式中AVIxy為類A與類B的平均距離，AVx和AVy為類內(nèi)樣本平均距離，K為權(quán)重值?？煞中訧A,B表達(dá)式中含有類內(nèi)樣本距離和類間樣本距離，因此能同時(shí)反映類A和類B的類間樣本距離和類內(nèi)樣本分布的情況。通過調(diào)節(jié)權(quán)重系統(tǒng)K(K取值為[2-4，2-3，…，23，24])，可生成最優(yōu)的二叉樹層次結(jié)構(gòu)。

ISVM-BT的計(jì)算流程圖如圖4所示。

圖4　ISVM-BT的計(jì)算流程Fig.4　Flowchart of ISVM-BT procedure

具體算法流程如下：

(1) 給定權(quán)值K的范圍(K=2n，n=-4，-3，…，4)，計(jì)算訓(xùn)練樣本的AV和AVI；

(2)選取某一權(quán)值K，計(jì)算當(dāng)前K下的可分性矩陣SI=[Ii,j],i,j=1,2,…,N,i≠j(SI為對稱矩陣)；

(17)

(3)對SI的每一行分別求和，并按值大小進(jìn)行排序，若數(shù)值相等，序號小的在前；

(4)子分類器的構(gòu)建：將可分性最大的類作為子分類器的第一個(gè)正類，余下的N-1類樣本集合作為子分類器的負(fù)類；

(5)重復(fù)步驟(4)的過程，直至建立N-1個(gè)SVM子分類器，生成二叉樹層次結(jié)構(gòu)。

(6)將測試樣本輸入當(dāng)前層次結(jié)構(gòu)下支持向量機(jī), 得出測試樣本的識別準(zhǔn)確率，并記錄；

(7) 令n=n+1，得到不同的權(quán)值K，重復(fù)步驟(2)～(6)，通過最高的測試樣本識別率確定最優(yōu)權(quán)值K。

3　層次模糊熵和改進(jìn)支持向量機(jī)的算法流程

基于層次模糊熵故障特征提取和ISVM-BT模式識別的軸承故障診斷方法流程如下：

(1)采用層次模糊熵對采集不同狀態(tài)下的軸承信號進(jìn)行層次k=3的模糊熵計(jì)算，得到8個(gè)節(jié)點(diǎn)的熵值，對每一個(gè)節(jié)點(diǎn)進(jìn)行模糊熵計(jì)算時(shí)設(shè)置嵌入維數(shù)m=2，相似容限r(nóng)=0.15*SD；

(2)將得到的8個(gè)節(jié)點(diǎn)熵值作為故障特征向量，同時(shí)將特征向量分為訓(xùn)練樣本集和測試樣本集；

(3)將訓(xùn)練樣本集輸入到ISVM-BT中訓(xùn)練，根據(jù)最高識別準(zhǔn)確率，選出最優(yōu)二叉樹層次結(jié)構(gòu)；

(4)將測試樣本輸入到由步驟(3)得到的最優(yōu)二叉樹層次結(jié)構(gòu)中做故障模式識別。

需要注意的是經(jīng)遺傳算法計(jì)算的罰函數(shù)C和核函數(shù)γ會略有不同，因此以5次優(yōu)化過程中最高識別率作為當(dāng)前權(quán)重K值對應(yīng)的識別準(zhǔn)確率。

4　工程應(yīng)用實(shí)例

4.1試驗(yàn)數(shù)據(jù)介紹

以美國西儲大學(xué)(Case Western Reserve University)軸承實(shí)驗(yàn)室滾動軸承為研究對象[17]，試驗(yàn)臺包括一個(gè)2馬力電機(jī)，一個(gè)轉(zhuǎn)矩傳感器和一個(gè)功率計(jì)，如圖5所示。本次試驗(yàn)采用驅(qū)動端的深溝型球軸承，其型號為6205-2RS JEM SKF。采用電火花機(jī)分別對軸承做3種故障的損傷：軸承內(nèi)圈故障(IRF)，外圈故障(ORF)和滾動體故障(BF)，并做4種不同程度損傷，損傷直徑分別為0.1778 mm, 0.3556 mm, 0.5334 m和0.7112 mm。振動加速度傳感器固定在驅(qū)動電機(jī)支撐軸的軸承座上方，并在四種不同工況下采集振動信號。采樣頻率為12 kHz，采樣時(shí)長為1 s，驅(qū)動電機(jī)的轉(zhuǎn)速為1730 r/min，載荷為2206.50 W。

圖5　軸承試驗(yàn)臺結(jié)構(gòu)簡圖Fig.5　The sketch of rolling bearing experiment system

4.2振動信號分析與故障特征提取

本次試驗(yàn)采用振動數(shù)據(jù)由3種故障狀態(tài)振動信號和正常狀態(tài)振動信號組成，每種故障狀態(tài)又包含不同程度的故障，每種狀態(tài)取60組數(shù)據(jù)，其中10組為訓(xùn)練樣本，剩余50組為測試樣本，數(shù)據(jù)長度N=2048，詳細(xì)信息如表1所示。8種不同故障類型和不同故障程度的軸承故障振動加速度信號的時(shí)域波形如圖6所示。

表1　試驗(yàn)數(shù)據(jù)列表

圖6　8種不同狀態(tài)下的振動信號時(shí)域波形圖Fig.6　The waveforms of rolling bearing vibration signal under eight different conditions

由圖6 可知，由于背景噪聲干擾和故障類型較多，難以從時(shí)域波形上對不同類型和不同程度的軸承故障進(jìn)行區(qū)分。首先采用MFE對不同狀態(tài)下的振動數(shù)據(jù)進(jìn)行處理，為了避免偶然性，每種狀態(tài)下選取50組數(shù)據(jù)進(jìn)行多尺度模糊熵分析，然后求均值，結(jié)果如圖7所示。

圖7　8種不同狀態(tài)下的振動信號多尺度模糊熵Fig.7　MFE values over 8 scales for analyzing 8 health bearing conditions

由圖7 可得出以下結(jié)論：首先正常狀態(tài)下，滾動軸承的振動信號在大尺度上熵值較大，故障狀態(tài)下熵值較小。這是由于正常滾動軸承振動信號是隨機(jī)振動信號[10]，在大尺度上信號無規(guī)則性更高，自相似性較低，因此熵值較大。相比于正常狀態(tài)，滾動軸承出現(xiàn)故障時(shí)，采集的振動信號由于受到周期性沖擊，規(guī)則性增加，自相似性降低，因此熵值偏小。其次3種軸承輕度故障的多尺度模糊在大尺度上排序?yàn)椋狠S承外圈故障<軸承內(nèi)圈故障<軸承滾動體故障。由于滾動軸承的外圈在工作時(shí)固定在基座上，不隨軸旋轉(zhuǎn)，出現(xiàn)故障時(shí)，測得振動信號自相似性最強(qiáng)，熵值最低；內(nèi)圈是隨著軸一起轉(zhuǎn)動，且相比于外圈故障，內(nèi)圈故障的傳遞路徑較長，測得振動信號雖具有一定沖擊性，但也具有干擾。因此相比于外圈故障，內(nèi)圈故障復(fù)雜性較高，熵值較大；滾動體在隨軸轉(zhuǎn)動同時(shí)還有自身的轉(zhuǎn)動，且傳遞路徑最遠(yuǎn)，因此理論上滾動體故障相比于其他兩種故障具有最高復(fù)雜度，熵值最大。

雖然通過多尺度模糊熵能夠區(qū)分軸承正常工況和3種不同類型輕度故障工況，然而對于同一故障類型不同故障程度(如內(nèi)圈輕度和中度故障)以及不同故障類型和故障程度(如內(nèi)圈輕度故障、外圈中度故障和滾動體輕度故障)的故障信號，采用多尺度模糊熵進(jìn)行計(jì)算時(shí)，在每個(gè)尺度上，模糊熵值很接近不能進(jìn)行有效區(qū)分。

圖8　8種不同狀態(tài)下的振動信號層次模糊熵Fig.8　HFE values of the 8 hierarchical decomposition nodes for analyzing 8 health bearing conditions

針對多尺度模糊熵在分析軸承不同故障類型和程度存在的不足，提出了采用多層次模糊熵進(jìn)行軸承故障特征的提取，設(shè)置分割層數(shù)k=3，軸承在不同狀態(tài)下層次模糊熵曲線如圖8所示。由圖8可得出如下結(jié)論。首先軸承正常狀態(tài)的層次模糊熵值隨分解節(jié)點(diǎn)數(shù)增加呈單調(diào)下降趨勢，說明振動信號的低頻成分包含主要的時(shí)間模式信息。其次當(dāng)軸承處于故障狀態(tài)時(shí)，振動信號在低頻部分(節(jié)點(diǎn)=3)和高頻部分(節(jié)點(diǎn)=7)熵值較大，說明故障信息不僅隱藏在振動信號的低頻部分，振動信號的高頻部分也包含重要的故障信息。第三每種狀態(tài)下選取50組數(shù)據(jù)進(jìn)行層次模糊熵分析，每一個(gè)分解節(jié)點(diǎn)處的不同狀態(tài)模糊熵值排序保持不變，說明層次模糊熵分析具有很好的穩(wěn)定性。

綜上所述，滾動軸承的振動信號的高頻部分也包含重要的故障信息，只考慮低頻部分的信息不能完全刻畫軸承故障的本質(zhì)，從而驗(yàn)證了層次模糊熵的優(yōu)越性，適合進(jìn)行軸承故障特征的提取。

4.3軸承故障診斷與結(jié)果分析

通過4.2節(jié)的分析，利用層次模糊熵提取軸承的故障特征。選取4.2節(jié)中8種工況各60組特征向量作為樣本集，其中以每種工況10 樣本作為訓(xùn)練樣本，剩余50組為測試樣本。權(quán)重K的取值范圍設(shè)為[2-4，2-3，…，23，24]，對于每一個(gè)K值對應(yīng)的二叉樹層次結(jié)構(gòu)，采用LibSVM工具箱對剩余的50組樣本進(jìn)行測試，得出在各個(gè)權(quán)值K下識別準(zhǔn)確率，識別結(jié)果如表2所示。從表2中可以看出K=21時(shí)，識別率達(dá)到100%，從而確定出最優(yōu)二叉樹層次結(jié)構(gòu)如圖9所示。注：N 表示正常工況；I1, I2, I3表示內(nèi)圈的輕度，中度和嚴(yán)重故障；O1, O2表示外圈的輕度和中度故障；B1, B2表示滾動體的輕度和重度故障。

表2　不同權(quán)重值K下的層次結(jié)構(gòu)和識別準(zhǔn)確率

圖9　最優(yōu)層次結(jié)構(gòu)下子分類器示意圖Fig.9　The order of multi-fault classifier under optimum hierarchy structure

為了評價(jià)ISVM-BT的識別性能，與單獨(dú)的類間樣本距離法(Inter-BT)類內(nèi)樣本分布法(Intra-BT)進(jìn)行對比，識別結(jié)果如表3所示。從表3可知ISVM-BT具有高的識別率，驗(yàn)證了其識別性能的優(yōu)越性。

為了驗(yàn)證層次模糊熵在故障特征提取能力上的優(yōu)越性，采用多尺度模糊熵對相同樣本進(jìn)行計(jì)算，為了使得輸入特征向量維數(shù)與層次模糊熵保持一致，設(shè)置尺度, 采用相同步驟對多尺度模糊熵進(jìn)行模式識別，K取2-1對應(yīng)最高識別率為98.75%，共出現(xiàn)5個(gè)測試樣本被誤分。其中3個(gè)內(nèi)圈輕度故障錯(cuò)分為滾動體輕度故障, 2個(gè)外圈中度故障錯(cuò)分為內(nèi)圈輕度故障，與上述的理論分析一致。

表3　3種不同分類器下的層次結(jié)構(gòu)和識別準(zhǔn)確率

同時(shí)為進(jìn)一步說明層次模糊熵的優(yōu)越性，將其與小波包分解的方法進(jìn)行對比。為了保證兩種方法在輸入特征向量維數(shù)上的一致性，使用3層小波包分解[18]。小波包分解的方法如下：首先對軸承振動信號進(jìn)行小波包分解，然后計(jì)算各頻帶分量的模糊熵值作為特征向量，最后輸入ISVM-BT進(jìn)行模式識別。當(dāng)K=20對應(yīng)最高識別率為99.25%，出現(xiàn)3個(gè)測試樣本被誤分。其中2個(gè)軸承內(nèi)圈輕度故障錯(cuò)分為內(nèi)圈中度故障和1個(gè)滾動體輕度故障錯(cuò)分為內(nèi)圈輕度故障。對比結(jié)果進(jìn)一步驗(yàn)證了層次模糊熵能有效提取故障特征，實(shí)現(xiàn)軸承不同故障類型和故障程度的診斷。同時(shí)對層次模糊熵和小波包分解的方法在計(jì)算效率方面進(jìn)行對比。采用長度N=2018的軸承振動信號作為測試信號，分別采用兩種方法進(jìn)行求解。層次模糊熵的計(jì)算時(shí)間為2.16 s，而小波包分解方法的計(jì)算時(shí)間為50.42 s，結(jié)果表明層次模糊熵在計(jì)算效率上優(yōu)于小波包分解的方法。

為避免采用單一分類器出現(xiàn)的偶然性，本文采用一對一方法(OAO)[19]、一對多方法(OAA)[20]、Inter-BT方法 和Intra-BT方法等不同多分類器分別對層次模糊熵法和多尺度模糊熵法進(jìn)行模式識別，識別結(jié)果如表4所示。由表4可得出以下結(jié)論：首先在每種多分類器下，相比于多尺度模糊熵，層次模糊熵都擁有較高的識別準(zhǔn)確率，排除了單一分類器識別準(zhǔn)確率高的偶然性，進(jìn)一步驗(yàn)證了層次模糊熵能夠更全面、更深刻的提取軸承故障信息。其次相比于其他多分類器，ISVM-BT對于層次模糊熵和多尺度模糊熵都具有最高識別準(zhǔn)確率，也說明了ISVM-BT在識別能力上具有優(yōu)越性。

表4　基于HFE和MFE的最優(yōu)參數(shù)和診斷結(jié)果

5　結(jié)　論

(1)在模糊熵的基礎(chǔ)上，結(jié)合層次分析思想提出一種新的特征提取方法——層次模糊熵。與多尺度模糊熵相比，層次模糊熵同時(shí)分析信號中的低頻分量和高頻分量的故障信息，對故障特征描述更加全面。

(2)給出了層次模糊熵的算法步驟，并利用仿真信號驗(yàn)證了該算法的優(yōu)越性。

(3)基于層次模糊熵提取滾動軸承的故障特征，應(yīng)用ISVM-BT實(shí)現(xiàn)了軸承不同故障類型和故障程度的準(zhǔn)確診斷，為軸承故障診斷提供了一種新思路。

[1]李輝，鄭海起，楊紹普. 基于EMD和Teager能量算子的軸承故障診斷研究[J]. 振動與沖擊，2008，27(10)：15—17,22.

Li Hui,Zheng Haiqi,Yang Shaopu. Bearing fault diagnosis based on emd and teager kaiser energy operator[J]. Journal of Vibration and Shock, 2008，27(10)：15—17,22.

[2]Qiu H, Lee J, Lin J, et al. Wavelet filter-based weak signature detection method and its application on rolling element bearing prognostics[J]. Journal of Sound and Vibration, 2006, 289: 1066—1090.

[3]張超, 陳建軍, 徐亞蘭. 基于 EMD 分解和奇異值差分譜理論的軸承故障診斷方法[J]. 振動工程學(xué)報(bào), 2011, 24(5): 539—545.

Zhao Chao, Chen Jianjun, Xu Yalan. A bearing fault diagnosis method based on EMD and difference spectrum theory of singular value[J]. Journal of Vibration Engineering, 2011, 24(5): 539—545.

[4]Pincus S M. Approximate entropy as a measure of system complexity [J]. Proceedings of the National. Academy of Sciences, 1991, 88: 2297—2301.

[5]Pincus S M. Approximate entropy as a complexity measure [J]. Chaos, 1995, 5(1): 110—117.

[6]Richman J S, Mo orman J R. Physiological time-series analysis using approximate entropy and sample entropy [J]. American Journal of Physiology-Heart Circulatory Physiology, 2000, 278: 2039—2049.

[7]陳偉婷. 基于熵的表面肌電信號特征提取研究 [D]. 上海:上海交通大學(xué), 2008.

Chen Weiting. A study of feature extraction from SEMG signal based on entropy [D]. Shanghai: Shanghai Jiao Tong University, 2008.

[8]Costa M, Goldberger A L, Peng C K. Multiscale entropy analysis of physiologic time series [J]. Phys. Rev Lett., 2002, 89: 062102.

[9]Costa M, Goldberger A L, Peng C K. Multiscale entropy analysis of biological signals [J]. Phys. Rev. E., 2005, 71: 021906.

[10]鄭近德, 陳敏均, 程軍圣, 等. 多尺度模糊熵及其在滾動軸承故障診斷中的應(yīng)用[J]. 振動工程學(xué)報(bào), 2014, 27(1):145—151.

Zheng Jinde,Chen Minjun,Cheng Junsheng.Multiscale fuzzy entropy and its application in rolling bearing fault diagnosis. [J].Journal of Vibration Engineering, 2014, 27(1):145—151.

[11]Jiang Ying, Peng C K, Xu Yuesheng. Hierarchical entropy analysis for biological signals [J]. Journal of Computational and Applied Mathematics, 2011, 236:728—742.

[12]朱可恒, 宋希庚, 薛冬新. 一種基于層次熵和 SVM 的軸承故障診斷方法[J]. 噪聲與振動控制, 2013, 33(6): 163—167.

Zhu Keheng, Song Xigeng, Xue Dongxin. Roller bearing fault diagnosis based on hierarchical entropy and SVM[J].Noise and Vibration Control, 2013, 33(6): 163—167.

[13]Xia Shiyu, Li Jiuxian, Xia Liangzheng, et al. Tree-structured support vector machines for multi-class classification[C]. Lecture Notes in Computer Science, Berlin, Heidelberg: Springer-Verlag, 2007,4493: 392—398.

[14]袁勝發(fā),褚福磊.次序二叉樹支持向量機(jī)多類故障診斷算法研究[J]. 振動與沖擊, 2009, 28(3):51—55.

Yuan Sheng-fa, Chu Fu-lei. Multi-class fault diagnosis based on support vector machineswith sequenced binary tree archtecture[J]. Journal of Vibration and Shock, 2009, 28(3):51—55.

[15]趙海洋, 徐敏強(qiáng), 王金東. 改進(jìn)二叉樹支持向量機(jī)及其故障診斷方法研究[J]. 振動工程學(xué)報(bào), 2013, 26(5): 764—770.

Zhao Hai-yang, Xu Min-qiang, Wang Jin-dong. An improved binary tree SVM and application for fault diagnosis[J]. Journal of Vibration Engineering, 2013, 26(5): 764—770.

[16]Richman J S, Moorman J R. Physiological time-series analysis using approximate entropy and sample entropy[J]. American Journal of Physiology-Heart and Circulatory Physiology, 2000,278:H2039-H2049.

[17]Bearing Data Center, Case Western Reserve University[WB/OL]. http://csegroups.case.edu/bearingdatacenter/pages/download-data-file.

[18]趙志宏, 楊紹普. 基于小波包變換與樣本熵的滾動軸承故障診斷[J]. 振動、測試與診斷, 2012, 32(4): 640—644.

Zhao Zhi-hong, Yang shao-pu. Roller bearing fault diagnosis based on wavelet packet transform and sample entropy[J]. Journal of Vibration, Measurement & Diagnosis, 2012, 32(4): 640—644.

[19]Knerr S, Personnaz L, Dreyfus G. Single-layer learning revisited: a stepwise procedure for building and training a neural network[C]//Fogelman-Souli F, Herault J ed. Neurocomputing: Algorithms, Architectures and Applications. Springer, Berlin, 1990.

[20]Vapnik V. The Nature of Statistical Learning Theory[M]. New York: Springer, 1995.

A study on rolling bearing fault diagnosis method based on hierarchical fuzzy entropy and ISVM-BT

LIYong-bo,XUMin-qiang,ZHAOHai-yang,HUANGWen-hu

(School of Astronautics, Harbin Institute of Technology, Harbin 150080, China)

A new rolling bearing fault feature extractor called hierarchical fuzzy entropy (HFE) is proposed in this paper, which is composedcomprises the of hierarchical procedure and the fuzzy entropy calculation. Compared with multi-scale fuzzy entropy (MFE) method, HFE method considers both the low and high frequency components of the vibration signals, which can provide a much more accurate estimation of entropy. Besides, improved support vector machine based binary tree SVM (ISVM-BT) has the priority of high recognition accuracy compared with other classifiers. HenceTherefore, in this paper we proposed a novel rolling bearing fault diagnosis method based on HFE and ISVM-BT is proposed in this paper. Firstly, HFE is utilized to extract fault features and then the fault features are fed into the ISVM-BT to automatically fulfill the fault patterns identifications. The experimental results show the proposed method is effective in recognizing the different categories and severities of rolling bearings.

fault diagnosis; hierarchical fuzzy entropy (HFE); improved support vector machine based binary tree(ISVM-BT); rolling bearing

2014-09-01；

2015-04-27

國家自然科學(xué)基金資助項(xiàng)目(10772061)

TH165+.3

A

1004-4523(2016)01-0184-09

10.16385/j.cnki.issn.1004-4523.2016.01.023

李永波(1986—),男，博士研究生。E-mail：liyongbo0532@126.com

徐敏強(qiáng)(1960—)，男，博士生導(dǎo)師，教授。E-mail：xumq@hit.edu.cn