周小龍， 劉薇娜， 姜振海； 馬風(fēng)雷

(1. 長春理工大學(xué) 機(jī)電工程學(xué)院， 長春 130022； 2. 長春工業(yè)大學(xué) 機(jī)電工程學(xué)院， 長春 130012)

基于改進(jìn)HHT和馬氏距離的齒輪故障診斷

周小龍1， 劉薇娜1， 姜振海2； 馬風(fēng)雷2

(1. 長春理工大學(xué) 機(jī)電工程學(xué)院， 長春 130022； 2. 長春工業(yè)大學(xué) 機(jī)電工程學(xué)院， 長春 130012)

針對齒輪振動信號非線性和非平穩(wěn)的特點(diǎn)，提出一種基于改進(jìn)希爾伯特-黃變換與馬氏距離相結(jié)合的故障診斷方法。利用自適應(yīng)白噪聲的完備經(jīng)驗(yàn)?zāi)B(tài)分解將齒輪振動信號分解成一系列固有模態(tài)函數(shù)，并采用敏感固有模態(tài)函數(shù)判別算法判斷出對故障信息敏感的模態(tài)函數(shù)；通過對敏感固有模態(tài)分量的局部希爾伯特瞬時(shí)能量譜的分析，得出信號能量隨時(shí)間變化的精確表達(dá)；以不同故障信號局部希爾伯特瞬時(shí)能量譜的最大峰值作為特征向量，采用馬氏距離對齒輪故障進(jìn)行狀態(tài)識別。試驗(yàn)結(jié)果表明，該方法可有效識提取齒輪故障特征，實(shí)現(xiàn)不同故障狀態(tài)識別。

齒輪；自適應(yīng)白噪聲完備經(jīng)驗(yàn)?zāi)B(tài)分解；瞬時(shí)能量譜；馬氏距離；故障診斷

齒輪是機(jī)械設(shè)備中的常用部件，而齒輪傳動也是機(jī)械設(shè)備中常見的傳動方式。通常，齒輪故障是導(dǎo)致機(jī)械設(shè)備失效的主要原因，因此，對齒輪進(jìn)行故障診斷具有重要意義[1]。

當(dāng)齒輪發(fā)生故障時(shí)，其振動信號往往是非線性、非平穩(wěn)的時(shí)變信號，對于該類信號時(shí)頻分析是有效的方法[2]。以經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)為核心的希爾伯特-黃變換(Hilbert-Huang Transform，HHT)方法是近年來由Huang等[3]提出的一種處理非線性、非平穩(wěn)信號的有效方法。然而在EMD分解過程中易產(chǎn)生模態(tài)混疊現(xiàn)象，這將影響EMD分解效果。為解決該問題，Wu等[4]提出了平均總體經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition，EEMD)方法，并被廣泛應(yīng)用于旋轉(zhuǎn)機(jī)械的故障診斷中。李輝等[5]提出一種基于EEMD和Teager-Huang變換的齒輪箱故障診斷方法；張超等[6]采用EEMD和支持向量機(jī)相結(jié)合的方法對齒輪故障進(jìn)行了有效診斷；吳小濤等[7]針對軸承振動信號的特點(diǎn)，提出了基于峭度準(zhǔn)則EEMD和改進(jìn)形態(tài)濾波方法的軸承故障診斷方法。雖然基于EEMD的HHT方法在機(jī)械故障診斷中取得了一定成果，但該方法仍然存在不足，并限制該方法在故障診斷領(lǐng)域的發(fā)展。其主要表現(xiàn)為：①EEMD分解后的固有模態(tài)函數(shù)(Intrinsic Mode Function，IMF)內(nèi)仍有一定幅值的噪聲殘留，運(yùn)算過程中選擇加入白噪聲的幅值和分解次數(shù)將直接影響著分解結(jié)果，如果參數(shù)選擇不合理，將會降低分解效率，并使分解結(jié)果中會產(chǎn)生較多虛假IMF[8-9]。②EEMD分解出的IMF中，只有部分IMF能夠反映故障特征，而其它的為干擾成分，如果干擾成分不剔除，將會影響HHT的故障診斷精度。

針對上述問題，提出一種基于自適應(yīng)白噪聲的完備經(jīng)驗(yàn)?zāi)B(tài)分解(Complete Ensemble Empirical Mode Decomposition with Adaptive Noise，CEEMDAN)局部Hilbert瞬時(shí)能量譜和馬氏距離相結(jié)合的故障診斷方法。首先，利用CEEMDAN將非平穩(wěn)的齒輪振動信號分解成一系列IMF，并采用敏感IMF判別算法判斷出對故障信息敏感的IMF；其次，根據(jù)HHT方法的基本原理構(gòu)建出敏感IMF分量的局部Hilbert瞬時(shí)能量譜，用于提取故障特征；最后，利用馬氏距離對齒輪故障狀態(tài)進(jìn)行識別。試驗(yàn)結(jié)果表明：該方法可有效識別齒輪故障狀態(tài)。

1 改進(jìn)的HHT

1.1 CEEMDAN

相較于EEMD方法，CEEMDAN在分解過程中，采用對加入高斯白噪聲進(jìn)行EMD分解的方法代替EEMD分解過程中每次加入的高斯白噪聲，可實(shí)現(xiàn)自適應(yīng)加噪[10]，并能在較少平均次數(shù)下，實(shí)現(xiàn)重構(gòu)誤差幾乎為零[11]。因此，CEEMDAN不僅能夠克服EMD的模態(tài)混疊問題，同時(shí)能夠提高分解效率，減少虛假IMF的產(chǎn)生。

步驟1與EEMD方法相同，在待分解信號x(t)中加入不同幅值的高斯白噪聲ωi(t)，并進(jìn)行I次分解。由EEMD對信號x(t)+ε0ωi(t)分解出的第1個(gè)IMF分量即為CEEMDAN的第1個(gè)模態(tài)分量。

步驟2在第一階段(k=1)，計(jì)算第1個(gè)殘余信號

(1)

步驟3對信號r1(t)+ε1E1(ωi(t)) (i=1,2,…,I)進(jìn)行i次分解，直到得到第1個(gè)EMD模態(tài)分量為止，并計(jì)算第2個(gè)模態(tài)分量

(2)

步驟4對其余每個(gè)階段(k=2,…,K)(i=1, 2,…,I)計(jì)算第k個(gè)殘余信號

(3)

步驟5分解x(t)+εkEk(ωi(t))(i=1,2,…,I)～第1個(gè)EMD模態(tài)分量，并定義第k+1個(gè)模態(tài)分量為

(4)

步驟6執(zhí)行步驟4，并將k+1。當(dāng)殘余信號的極值點(diǎn)個(gè)數(shù)不超過2個(gè)的時(shí)候分解結(jié)束。

分解結(jié)束時(shí)，設(shè)所有模態(tài)分量的數(shù)量為K，則最終的殘差信號為

(5)

因此，原信號經(jīng)CEEMDAN分解后可表示為

(6)

由CEEMDAN分解過程可知，其分解過程相對于EEMD更加完整，能夠?qū)Υ纸膺M(jìn)行精確重構(gòu)。在每一個(gè)模態(tài)分解階段，根據(jù)系數(shù)εk可實(shí)現(xiàn)信噪比的自主選擇。

為驗(yàn)證CEEMDAN的有效性，仿真一個(gè)由正弦波和沖擊成分組成的信號。分別使用EEMD和CEEMDAN對仿真信號進(jìn)行分解，分解過程中加入的白噪聲標(biāo)準(zhǔn)差為0.2，分解次數(shù)500次，分解結(jié)果如圖1和圖2所示。由圖1和圖2可知，EEMD和CEEMDAN都將正弦信號和沖擊信號分解出來，但由于參數(shù)選擇問題，EEMD分解結(jié)果中出現(xiàn)了較多虛假IMF分量。而CEEMDAN有效減少了虛假IMF分量的產(chǎn)生，提高了分解準(zhǔn)確性。

(a)仿真信號 (b)IMF1 (c)IMF2

(d)IMF3 (e)IMF4 (f)IMF5圖1 EEMD分解結(jié)果Fig.1 EEMD decomposition results

(a)仿真信號 (b)IMF1 (c)IMF2 (d) IMF3圖2 CEEMDAN分解結(jié)果Fig.2 CEEMDAN decomposition results

圖3為EEMD和CEEMDAN兩種分解方法的重構(gòu)誤差。由圖3可知，CEEMDAN的重構(gòu)誤差小于10-15，而EEMD由于受到參數(shù)選擇的影響，出現(xiàn)較大重構(gòu)誤差。因此，EEMD在一定程度上影響了信號分解的完整性，若想減小重構(gòu)誤差，需要提高分解次數(shù)，但這將導(dǎo)致計(jì)算規(guī)模的增加并降低計(jì)算速率。

(a)EEMD重構(gòu)誤差 (b) CEEMDAN重構(gòu)誤差圖3 不同分解方法的重構(gòu)誤差Fig.3 Reconstruction error of different decomposition methods

1.2 敏感IMF判別算法

信號經(jīng)CEEMDAN分解最終得到一系列IMF，雖然CEEMDAN能夠有效減少虛假IMF分量的產(chǎn)生，但對于機(jī)械故障診斷而言，虛假IMF會影響故障診斷的準(zhǔn)確性，同時(shí)，真實(shí)IMF中也只有部分IMF是包含故障特征或?qū)收咸卣髅舾械摹Ｒ虼?，在對IMF進(jìn)行Hilbert變換前，需判別各IMF分量與所診斷故障的相關(guān)程度，保證故障特征提取的精度與診斷的有效性。在此采用一種基于信號相關(guān)性的敏感IMF判斷法[12]。具體算法如下:

(1)計(jì)算故障信號x(t)的各IMF分量c1(t),…，cn(t)與原信號x(t)的相關(guān)系數(shù)αi；

(2)計(jì)算故障信號x(t)的各IMF分量c1(t),…，cn(t)與正常狀態(tài)信號y(t)的相關(guān)系數(shù)βi；

(3)由兩個(gè)相關(guān)系數(shù)計(jì)算得到各IMF分量包含故障信息的敏感相關(guān)系數(shù)γi，記作

γi=|αi-βi| (i=1,2,…,n)

(7)

(4)計(jì)算各IMF的故障敏感系數(shù)

(8)

(9)

(6)判斷出敏感系數(shù)差值最大值所對應(yīng)的序號k，那么排序后的前k個(gè)IMF即為故障敏感的IMF。

對于該方法，其主要原理是：與故障信號x(t)相關(guān)系數(shù)越大的IMF分量，所包含信號x(t)的故障信息越多；計(jì)算與正常狀態(tài)信號y(t)的相關(guān)系數(shù)βi，可知每個(gè)IMF分量中與故障信息無關(guān)的正常信息。敏感IMF判別方法將各IMF與故障信號的相關(guān)性同正常信號的相關(guān)性相結(jié)合，配合敏感系數(shù)，能突顯出故障信息，削弱正常信息的影響。

對于旋轉(zhuǎn)機(jī)械系統(tǒng)而言，當(dāng)系統(tǒng)構(gòu)件出現(xiàn)故障時(shí)，會產(chǎn)生瞬態(tài)沖擊。圖1中的仿真信號可以看成一故障信號，經(jīng)CEEMDAN分解后，IMF1為故障特征成分，而IMF2為虛假分量，IMF3為正常狀態(tài)信息。采用敏感IMF判別算法計(jì)算各相關(guān)系數(shù)和敏感系數(shù)，結(jié)果如圖4所示。

圖4 仿真信號的敏感IMF判別結(jié)果Fig.4 The sensitive correlation coefficient results of simulation signal

由圖4可知，故障特征成分IMF1被準(zhǔn)確判別出來，而虛假成分IMF2和含有正常狀態(tài)信息成分的IMF3被剔除，表明該算法的準(zhǔn)確性。

1.3 局部Hilbert瞬時(shí)能量譜

對CEEMDAD分解后的各敏感IMF分量進(jìn)行Hilbert變換

(10)

構(gòu)造解析信號

(11)

進(jìn)一步可求出瞬時(shí)頻率

(12)

將上述敏感IMF的瞬時(shí)頻率和瞬時(shí)幅值組合，可得信號的局部Hilbert譜

(13)

(14)

由此可以定義

(15)

式(15)為局部Hilbert瞬時(shí)能量譜(簡稱：局部瞬時(shí)能量譜)，它反映了信號能量隨時(shí)間分布的情況。

常用的齒輪故障診斷方法一般基于振動能量的變化，但由于齒輪故障多為發(fā)生在各別齒上，所引起的振動能量變化在整個(gè)齒輪振動中所占比例較小，因此對于這種變化檢測比較困難。當(dāng)齒輪存在局部故障時(shí)，故障齒輪會產(chǎn)生較大的動載荷，并會影響到相鄰齒，當(dāng)故障程度增加時(shí)，動載荷增加，并會使振動能量發(fā)生較大變化。由敏感IMF求得的局部瞬時(shí)能量譜是對故障能量隨時(shí)間變化的精確表示，可有效提取相關(guān)故障信息。

由之前分析可知，仿真信號的IMF1是敏感IMF，圖5為其局部瞬時(shí)能量譜。在發(fā)生故障時(shí)，信號中的動載荷增加、振動加強(qiáng)，使相關(guān)時(shí)間段內(nèi)信號的振動能量出現(xiàn)較大變化，結(jié)合時(shí)域信號特征，在振動開始和結(jié)束的階段沖擊最強(qiáng)。在圖5中這些現(xiàn)象得到準(zhǔn)確反映，顯示出局部瞬時(shí)能量譜在故障特征提取方面的有效性。

圖5 仿真信號的局部瞬時(shí)能量譜Fig.5 The local instantaneous energy spectrum of simulation signal

2 試驗(yàn)研究

為驗(yàn)證所提方法在故障特征提取方面的有效性，在齒輪試驗(yàn)臺上進(jìn)行了齒輪具有斷齒、磨損和點(diǎn)蝕等三種故障的試驗(yàn)，其中斷齒和點(diǎn)蝕故障發(fā)生在從動齒輪，磨損故障發(fā)生在驅(qū)動齒輪。試驗(yàn)所用的齒輪試驗(yàn)臺型號為PQZZ-II型，圖6為試驗(yàn)臺裝置圖。

在齒輪故障模擬時(shí)，使用鐵錘將齒輪故障實(shí)驗(yàn)臺的從動齒輪的一個(gè)齒敲掉，折斷線位于分度圓附近，用于模擬斷齒故障；對另一從動齒輪的一個(gè)齒面進(jìn)行點(diǎn)焊，并采用油石磨蝕使焊接殘留的焊渣剝落，在其齒面上形成坑點(diǎn)，點(diǎn)蝕面積率為2%，用來模擬中等程度的點(diǎn)蝕故障；為模擬齒輪磨損故障，設(shè)置驅(qū)動齒輪的磨損為單齒齒面嚴(yán)重磨損，齒面磨損原因?yàn)槿藶槟ヤS，磨損程度為齒厚的15%。

圖6 齒輪試驗(yàn)臺Fig.6 Gear tester

試驗(yàn)過程中，電動機(jī)轉(zhuǎn)頻為50 Hz，對旋轉(zhuǎn)機(jī)械而言，軸承座是設(shè)備易產(chǎn)生劣化現(xiàn)象的易損點(diǎn)，為保證所測數(shù)據(jù)的準(zhǔn)確性，同時(shí)為全面監(jiān)測齒輪箱的振動狀態(tài)，加速度傳感器安裝在齒輪箱軸承座上。數(shù)據(jù)采集卡采用ADA16-8/2(LPCI)型高速多功能采集卡，設(shè)置采樣頻率為5 120 Hz，采樣時(shí)間為1 s，采樣精度為16 bt。平行軸齒輪箱中各齒輪均為標(biāo)準(zhǔn)直齒輪，其主要參數(shù)如表1所示。

表1 平行軸齒輪箱參數(shù)

試驗(yàn)采集到的各種狀態(tài)條件下齒輪加速度信號如圖7所示。從圖7可知，在信號采集過程中，背景噪聲較強(qiáng)，而由于沒有消噪室，使信號的信噪比降低，這將影響改進(jìn)HHT分析的準(zhǔn)確性，同時(shí)增加故障檢測難度。為此，采用快速獨(dú)立分量分析(Fast Independent Component Analysis，F(xiàn)astICA)法對時(shí)域信號進(jìn)行降噪處理，降噪后信號如圖8所示。

(a)齒輪正常狀態(tài)信號 (b)齒輪斷齒故障信號 (c)齒輪磨損故障信號 (d)齒輪點(diǎn)蝕故障信號圖7 齒輪振動信號Fig.7 Gear vibration signals

(a)齒輪正常狀態(tài)降噪信號 (b)齒輪斷齒降噪信號 (c)齒輪磨損降噪信號 (d)齒輪點(diǎn)蝕降噪信號圖8 降噪后齒輪振動信號Fig.8 Noise reduction gear vibration signals

從圖8可知，環(huán)境噪聲等因素對于信號的干擾已得到明顯改善。大部分的無用高頻分量被濾去，降噪信號的波形更加突顯出原信號的信息。各種齒輪故障信號在時(shí)域中均表現(xiàn)為脈沖沖擊，同高斯幅值調(diào)制信號相類似。雖然不同故障的時(shí)域信號存在差異，但是以此為依據(jù)不能準(zhǔn)確區(qū)分故障類型。采用基于CEEMDAN和敏感IMF的改進(jìn)HHT方法對降噪后齒輪故障信號進(jìn)行分析，3類故障信號經(jīng)CEEMDAN分解后，斷齒信號得到6個(gè)IMF分量，齒輪磨損和點(diǎn)蝕信號分別得到5個(gè)IMF分量。根據(jù)敏感IMF判別算法，計(jì)算各故障類型IMF分量的相關(guān)系數(shù)(αn，βn)和敏感系數(shù)(λn),如圖9所示，通過敏感IMF判別算法，選擇斷齒信號的IMF1、磨損信號的IMF1和點(diǎn)蝕信號的IMF1 ～ IMF3為敏感IMF分量。各故障信號中其它IMF分量中所包含的是正常信號與故障信號的共同信息，所以將它們剔除。

(a)齒輪斷齒信號的敏感IMF判別結(jié)果 (b)齒輪磨損信號的敏感IMF判別結(jié)果 (c)齒輪點(diǎn)蝕信號的敏感IMF判別結(jié)果圖9 故障信號的敏感IMF判別結(jié)果Fig.9 The sensitive correlation coefficient results of fault signals

對齒輪各狀態(tài)下的敏感IMF分量進(jìn)行Hilbert變換，相應(yīng)的局部瞬時(shí)能量譜,如圖10所示。從圖10可知，相較于齒輪正常嚙合信號，各故障信號的局部瞬時(shí)能量譜具有明顯的沖擊特征，由于齒輪存在故障，這將導(dǎo)致齒輪在嚙合過程中產(chǎn)生瞬態(tài)沖擊的較大能量，不同故障類型所產(chǎn)生的瞬態(tài)沖擊能量的大小也不相同。對各狀態(tài)下信號的局部瞬時(shí)能量譜分析可知，當(dāng)齒輪正常嚙合時(shí)，局部瞬時(shí)能量譜有4個(gè)主要峰值，峰值均值為55.75，其中最大峰值為70；齒輪斷齒狀態(tài)下有4個(gè)主要峰值，其平均值為568.25，其中最大峰值為762； 齒輪磨損狀態(tài)下有3個(gè)主要峰值，其平均值為2 430.47，其中最大峰值為3 752；齒輪點(diǎn)蝕狀態(tài)下有4個(gè)主要峰值，其平均值為1 024.21，其中最大峰值為1 315。由上述分析可知，對于齒輪各種狀態(tài)信號而言，瞬時(shí)能量譜峰值的最大值和均值存在明顯差別，可將其作為故障特征進(jìn)行診斷研究。

(a)齒輪正常信號局部 (b)齒輪斷齒信號局部 (c)齒輪磨損信號局部 (d)齒輪點(diǎn)蝕信號局部 瞬時(shí)能量譜 瞬時(shí)能量譜 瞬時(shí)能量譜 瞬時(shí)能量譜圖10 齒輪不同狀態(tài)的局部瞬時(shí)能量譜Fig.10 Local instantaneous energy spectrum of different gear conditions

3 馬氏距離的故障診斷分類

馬氏距離是多元數(shù)理統(tǒng)計(jì)理論中常用的判別方法之一，用于表示數(shù)據(jù)的協(xié)方差距離。該方法的計(jì)算是建立在總體樣本的基礎(chǔ)上，排除了變量之間相關(guān)性的干擾，是一種衡量兩個(gè)未知樣本集相似程度的有效方法[13]。由于算法簡單，并適用于小樣本問題的處理，因此，被廣泛應(yīng)用于故障診斷[14]領(lǐng)域。

根據(jù)上節(jié)分析，選擇不同狀態(tài)下齒輪瞬時(shí)能量譜的最大值作為訓(xùn)練樣本S進(jìn)行馬氏距離分析，具體步驟如下：

步驟1在齒輪故障綜合診斷實(shí)驗(yàn)臺上分別獲取四種工況各30組數(shù)據(jù)；

步驟2為消除背景噪聲對于故障診斷準(zhǔn)確性的影響，采用FastICA降噪；

步驟3對降噪后各個(gè)樣本信號進(jìn)行改進(jìn)HHT分析，得到瞬時(shí)能量譜；

步驟4以瞬時(shí)能量譜的最大值作為訓(xùn)練樣本S；

步驟6對待測齒輪故障信號x(t)按步驟2和步驟3求出其待測特征值Sx;

(16)

式中,k=1,2,3,4；分別代表齒輪正常、磨損、斷齒和點(diǎn)蝕四種工作狀態(tài)。

步驟8比較d1、d2、d3和d4的大小，并將最小判別距離所對應(yīng)的狀態(tài)作為待測信號的故障類型。

基于局部瞬時(shí)能量譜和馬氏距離的齒輪故障診斷法流程圖，如圖11所示。

圖11 齒輪故障診斷方法流程Fig.11 The technology process of gear fault diagnosis method

試驗(yàn)表明：通過小樣本即可實(shí)現(xiàn)對齒輪故障的識別，對于未知樣本的分類速度和識別準(zhǔn)確率都較高，40組數(shù)據(jù)中出現(xiàn)了2組誤判，總體識別正確率為95%。而誤判數(shù)據(jù)的產(chǎn)生可能是由于測量過程中的測量誤差和各種干擾因素的影響。部分識別結(jié)果如表2所示。表2部分識別結(jié)果

表2 部分識別結(jié)果

按上述方法，分別對制動力矩的輸出電流在0.1 A和0.2 A工況條件下的60組數(shù)據(jù)進(jìn)行識別，識別結(jié)果與實(shí)際故障類型基本一致，其中正常狀態(tài)全部識別正確，磨損和斷齒故障各有1組識別錯(cuò)誤，齒輪點(diǎn)蝕有2組識別錯(cuò)誤，總體診斷正確率為93.33%，進(jìn)一步表明該方法的有效性。

4 結(jié) 論

(1)利用CEEMDAN分析信號，不但能夠克服EMD分解過程中出現(xiàn)的模態(tài)混疊現(xiàn)象，同時(shí)可以解決EEMD分解效率低的問題。

(2)通過敏感IMF算法可以有效剔除與故障信息無關(guān)的IMF，增加故障特征提取的準(zhǔn)確性。

(3)Hilbert瞬時(shí)能量譜能夠準(zhǔn)確反映出齒輪故障振動信號的能量隨時(shí)間的分布情況，利用Hilbert瞬時(shí)能量譜可以較好的提取出齒輪故障信號的特征。

(4)基于改進(jìn)HHT瞬時(shí)能量譜和馬氏距離的故障診斷方法可以有效識別齒輪的故障狀態(tài)，滿足齒輪故障診斷要求，同時(shí)，實(shí)際應(yīng)用時(shí)應(yīng)盡量排除各種干擾因素，保證數(shù)據(jù)樣本的數(shù)量。

[ 1 ] 于德介，程軍圣，楊宇. Hilbert-Huang 變換在齒輪故障診斷中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào)， 2005，41(6)： 102-107.

YU Dejie，CHENG Junsheng，YANG Yu. Application of Hilbert-Huang transform method to gear fault diagnosis[J]. Chinese Journal of Mechanical Engineering, 2005，41(6)： 102-107.

[ 2 ] 程軍圣，楊宇，于德介. 局部均值分解方法及其在齒輪故障診斷中的應(yīng)用[J]. 振動工程學(xué)報(bào)，2009，22(1)： 76-84.

CHENG Junsheng，YANG Yu，YU Dejie. The local mean decomposition method and its application to gear fault diagnosis[J]. Journal of Vibration Engineering, 2009，22(1)： 76-84.

[ 3 ] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis [J]. Proceeding of the Royal Society of London-Series A: Mathematical, Physical and Engineering Sciences, 1998, 454(1971): 903-995.

[ 4 ] WU Z, HUANG N E. Ensemble empirical mode decomposition: a noise-assisted data analysis method [J]. Advances in Adaptive Data Analysis, 2009, 1(1): 1-41.

[ 5 ] 李輝，鄭海起，唐力偉. 基于EEMD 和THT 的齒輪故障診斷方法[J]. 振動、測試與診斷，2011，31(4)： 496-500.

LI Hui，ZHENG Haiqi，TANG Liwei. Gear fault diagnosis method based on EEMD and THT [J]. Journal of Vibration，Measurement amp; Diagnosis, 2011，31(4)： 496-500.

[ 6 ] 張超，陳建軍，郭迅. 基于EEMD能量熵和支持向量機(jī)的齒輪故障診斷方法[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012，43(3)： 932-939.

ZHANG Chao，CHEN Jianjun，GUO Xun. Gear fault diagnosis method based on ensemble empirical mode decomposition energy entropy and support vector machine[J]. Journal of Central South University (Science and Technology), 2012，43(3)： 932-939.

[ 7 ] 吳小濤，楊錳，袁曉輝，等. 基于峭度準(zhǔn)則EEMD及改進(jìn)形態(tài)濾波方法的軸承故障診斷[J]. 振動與沖擊，2015，34(2)： 38-44.

WU Xiaotao, YANG Meng, YUAN Xiaohui, et al. Bearing fault diagnosis using EEMD and improved morphological filtering method based on kurtosis criterion [J]. Journal of Vibration and Shock, 2015，34(2)： 38-44.

[ 8 ] TORRES M E, COLOMINAS M A, SCHLOTTHAUER G, et al. A complete ensemble empirical mode decomposition with adaptive noise [C] // IEEE International Conference on Speech and Signal Processing. Prague: ICASSP, 2011: 4144-4147.

[ 9 ] COLOMINAS M A, SCHLOTTHAUER G, TORRES M E. Improved complete ensemble EMD: a suitable tool for biomedical signal processing [J]. Biomedical Signal Processing amp; Control, 2014, 14(1): 19-29.

[10] 韓慶陽，孫強(qiáng)，王曉東，等. CEEMDAN去噪在拉曼光譜中的應(yīng)用研究[J]. 激光與光電子學(xué)進(jìn)展，2015，52(11)： 267-274.

HAN Qingyang，SUN Qiang，WANG Xiaodong，et al. Application of CEEMDAN in Raman spectroscopy denoising[J]. Laser amp; Optoelectronics Progress, 2015，52(11)： 268-274.

[11] 李軍，李青. 基于CEEMDAN-排列熵和泄漏積分ESN的中期電力負(fù)荷預(yù)測研究[J]. 電機(jī)與控制學(xué)報(bào)，2015，19(8)： 70-80.

LI Jun，LI Qing. Medium term electricity load forecasting based on CEEMDAN-permutation entropy and ESN with leaky integrator neurons[J]. Electric Machines and Control, 2015，19(8)： 70-80.

[12] 雷亞國. 基于Hilbert-Huang變換的機(jī)械故障診斷[J]. 機(jī)械工程學(xué)報(bào)，2011，47(5)： 71-77.

LEI Yaguo. Machinery fault diagnosis based on improved Hilbert-Huang transform[J]. Journal of Mechanical Engineering, 2011，47(5)： 71-77.

[13] 劉昱，張俊紅，畢鳳榮，等. 基于LMD邊際譜的柴油機(jī)氣門故障診斷[J]. 內(nèi)燃機(jī)學(xué)報(bào)，2014，35(6)： 96-100.

LIU Yu，ZHANG Junhong，BI Fengrong，et al. Diesel engine valve fault diagnosis based on LMD marginal spectrum[J]. Chinese Internal Combustion Engine Engineering, 2014，35(6)： 96-100.

[14] 駱志高，李旭東，趙俊麗，等. 利用馬氏距離判別法準(zhǔn)確實(shí)現(xiàn)對裂紋的識別[J]. 振動與沖擊，2013，32(21)：186-188.

LUO Zhigao, LI Xudong, ZHAO Junli, et al. Crack identification with mahalanobis distance discrimination method [J]. Journal of Vibration and Shock, 2013，32(21)：186-188.

GearfaultdiagnosisbasedonimprovedHHTandMahalanobisdistance

ZHOU Xiaolong1， LIU Weina1， JIANG Zhenhai2， MA Fenglei2

(1. College of Mechanical and Electric Engineering, Changchun University of Science and Technology, Changchun 130022, China; 2.School of Mechatronic Engineering, Changchun University of Technology, Changchun 130012, China)

In view of nonlinear and non-stationary characteristics of gear vibration signals, a fault diagnosis method based on improved Hilbert-Huang transform and Mahalanobis distance was proposed. The gear vibration signals were decomposed by complete ensemble empirical mode decomposition with adaptive noise, the intrinsic mode functions were obtained and sensitive intrinsic mode functions were selected by the sensitivity evaluation method. Then, the local Hilbert instantaneous energy spectrum of the sensitive intrinsic mode components was analyzed, and the fault information can be extracted from the distribution of the energy of the gear vibration signal with the change of time. Finally, the maximum peak value of the local Hilbert instantaneous energy spectrum was treated as the fault features and the Mahalanobis distance method was used for judging the gear fault. Experimental results show that the method can effectively extract gear fault features and apply for different fault identification.

gear; complete ensemble empirical mode decomposition with adaptive noise; instantaneous energy spectrum; Mahalanobis distance; fault diagnosis

國家自然科學(xué)基金資助項(xiàng)目(51075041)；吉林省教育廳科技發(fā)展項(xiàng)目(2014124)

2016-05-11 修改稿收到日期： 2016-09-20

周小龍 男，博士生，1987年生

劉薇娜 女，博士，教授，1956年生

TH17

A

10.13465/j.cnki.jvs.2017.22.034