楊恭勇，周小龍，李家飛，梁秀霞

(1.東北電力大學(xué) 工程訓(xùn)練教學(xué)中心 吉林 吉林 132012；2.河南信宇石油機(jī)械制造股份有限公司 河南 濮陽 457001)

?

局部Hilbert邊際能量譜在滾動軸承故障診斷中的應(yīng)用

楊恭勇1，周小龍1，李家飛2，梁秀霞2

(1.東北電力大學(xué) 工程訓(xùn)練教學(xué)中心 吉林 吉林 132012；2.河南信宇石油機(jī)械制造股份有限公司 河南 濮陽 457001)

針對滾動軸承振動信號非線性的特點(diǎn)，提出一種基于局部希爾伯特邊際能量譜和馬氏距離判別法相結(jié)合的故障診斷方法。首先，采用鏡像延拓方法抑制經(jīng)驗(yàn)?zāi)B(tài)方法將待分析信號分解成多階固有模態(tài)函數(shù)和的形式，并根據(jù)相關(guān)性判別算法選取含有主要故障信息的IMF分量；其次，利用局部Hilbert邊際能量譜提取故障信號能量特征參數(shù)；最后，利用馬氏距離算法對滾動軸承的工作狀態(tài)進(jìn)行判別。通過相關(guān)試驗(yàn)證明了此方法的有效性，具有一定的工程應(yīng)用價(jià)值。

經(jīng)驗(yàn)?zāi)B(tài)分解；局部Hilbert邊際能量譜；馬氏距離；滾動軸承；故障診斷

滾動軸承為各種機(jī)械設(shè)備中不可缺少的用于支撐轉(zhuǎn)軸的重要部件[1]。對于存在故障的滾動軸承故障信號而言，其常常具有非平穩(wěn)的時(shí)變特性[2]。對于此類信號，常采用時(shí)頻分析方法對其進(jìn)行故障診斷研究。在眾多時(shí)頻分析方法中，小波變換的應(yīng)用最為廣泛。然而該方法具有自身的局限性[3]，在分析非平穩(wěn)信號時(shí)往往難以得到滿意的效果。希爾伯特-黃變換(Hilbert-Huang Transform，HHT)[4]方法是近年來得到研究人員廣泛關(guān)注的一種新的信號處理和分析方法，HHT自身的自適應(yīng)性使得該方法在旋轉(zhuǎn)機(jī)械故障診斷領(lǐng)域得到廣泛應(yīng)用[5-6]，但該HHT方法的端點(diǎn)效應(yīng)和虛假模態(tài)等問題嚴(yán)重影響了此方法在該領(lǐng)域的發(fā)展。

基于上述分析，并結(jié)合滾動軸承振動信號的特性，本文提出一種基于局部Hilbert邊際能量譜和馬氏距離判別算法相結(jié)合的故障診斷方法。并通過試驗(yàn)研究，驗(yàn)證了所提方法的有效性。

1 改進(jìn)的HHT

HHT方法在處理旋轉(zhuǎn)機(jī)械故障診斷問題方面的不足主要表現(xiàn)為：端點(diǎn)效應(yīng)[7]和虛假固有模態(tài)函數(shù)(I)IMF分量[8]等問題。這些問題嚴(yán)重限制了HHT方法診斷的準(zhǔn)確性。在此，采用邊界局部特征尺度延拓法[9]與相關(guān)性分析法相結(jié)合的改進(jìn)HHT方法解決上述問題。

1.1 端點(diǎn)優(yōu)化對稱延拓算法

EMD方法是根據(jù)信號的局部特性將其分解成多階表征信號特性的IMF，據(jù)此所采用的端點(diǎn)優(yōu)化對稱延拓算法的步驟如下：

(1)

(2)以b、a為中心將所延拓的信號序列側(cè)進(jìn)行對稱延拓以保證信號序列的結(jié)構(gòu)不變性，則經(jīng)對稱延拓后的信號序列h：

(2)

式中：i=0,1,…,n-1。并采用三次樣條插值法獲得此信號序列的包絡(luò)線S。當(dāng)S與hi偏差最小時(shí)的信號序列所延拓的端點(diǎn)為最優(yōu)延拓點(diǎn)，此延拓點(diǎn)可由號hi與其相應(yīng)包絡(luò)Si的偏差評價(jià)函數(shù)來確定：

(3)

在?ε/?b=0且?ε/?a=0時(shí)，求出b、a的值。

由EMD算法可知，EMD出現(xiàn)端點(diǎn)效應(yīng)現(xiàn)象的主要原因是在信號端點(diǎn)處未做約束，而此方法正是根據(jù)信號端點(diǎn)處的自身特性來解決此問題，所以此方法能夠有效解決EMD分解過程中所存在的端點(diǎn)效應(yīng)問題。

1.2 虛假IMF判別算法

由于EMD算法具有正交性，所以經(jīng)EMD分解得到的各IMF分量與原信號間應(yīng)具有較好的相關(guān)性[10]，由此，可以求取各IMF分量同所分解信號間的相關(guān)系數(shù)來選取對信號信息敏感的IMF分量。從而有效避免虛假模態(tài)分量問題對于分析準(zhǔn)確性的影響。

設(shè)信號x(t)經(jīng)改進(jìn)EMD分解后得到n階IMF分量，即c1(t)，c2(t)，…，cn(t)，則其相關(guān)系數(shù)ρi為

(4)

分別計(jì)算各階IMF分量同原始信號間的相關(guān)系數(shù)，并由信號自身特性選取門限閾值判別出包含主要故障信息的IMF。

2 基于局部Hilbert邊際能量譜的特征能量法

對于滾動軸承而言，當(dāng)它自身出現(xiàn)故障時(shí)，會使?jié)L動軸承系統(tǒng)在固有頻率段的信號能量發(fā)生變化[11]，由之前分析可知，經(jīng)HHT方法所獲得的信號的Hilbert譜是信號能量的一種完整的時(shí)-頻分布，這為以Hilbert譜為基礎(chǔ)對滾動軸承進(jìn)行故障診斷提供了理論依據(jù)。

EMD分解出的每一階IMF分量的頻率和能量(幅值)都是時(shí)間的函數(shù)，所以，信號的Hilbert譜H(ω,t)準(zhǔn)確地描述了所分析信號的幅值同其時(shí)間和頻率的變化及分布情況。若將x(t)2作為信號的能量密度，經(jīng)HHT分析后，信號的Hilbert譜H2(ω,t)同樣具有此物理意義，故特將H2(ω,t)定義為信號的Hilbert能量譜。它在時(shí)-頻區(qū)間上表述了信號的能量變化情況，略去信號殘差后，不同工況下滾動軸承的振動信號在HHT分析前后應(yīng)滿足能量守恒，則有：

(5)

由此可定義E(ω)為Hilbert邊際能量譜，即：

(6)

由式(6)可知，信號能量隨頻率的變換情況可在E(ω)中得到準(zhǔn)確反映。而當(dāng)滾動軸承出現(xiàn)故障時(shí)，其固有頻率段的能量會有有所增大，所以可選擇對敏感IMF分量進(jìn)行Hilbert分析，求出表征信號主要狀態(tài)特性的局部Hilbert譜H′(ω,t)：

(7)

將E′(ω)定義為信號x(t)的邊際能量譜，它反應(yīng)了在固有頻率段的能量隨頻率的變化情況，所以可進(jìn)一步定義在固有頻率區(qū)間ω1～ω2的能量為特征向量：

(8)

式中：ω1～ω2為固有頻率段，S為相應(yīng)工作狀態(tài)下滾動軸承振動信號的能量。因此，可將S作為滾動軸承故障的特征值進(jìn)行后續(xù)診斷研究。

3 基于M距離的滾動軸承故障診斷

馬氏距離是多元數(shù)理統(tǒng)計(jì)理論中常用的判別方法之一，由于該方法不受量綱的影響，同時(shí)算法簡單，計(jì)算速度快，因此，非常適用于小樣本問題的分析[12]。

選擇之前計(jì)算的滾動軸承各不同狀態(tài)下在固有頻率段的能量S作為訓(xùn)練樣本，進(jìn)行馬氏距離分析，具體步驟如下：

(1) 按一定的采樣頻率fs對不同工況下滾動軸承的故障信號進(jìn)行采樣，在此，滾動軸承的共有正常狀態(tài)、內(nèi)圈故障狀態(tài)和外圈故障狀態(tài)等三種工況，每種工況的采用次數(shù)為30次；

(2) 對每個(gè)所采集到的振動信號進(jìn)行改進(jìn)EMD分解，獲得其IMF分量；

(3) 采用基于相關(guān)性的虛假IMF判別算法，選取敏感IMF分量進(jìn)行Hilbert分析，獲取其局部Hilbert譜H′(ω,t)；

(4) 根據(jù)公式(7)和公式(8)求出每個(gè)樣本信號的特征能量函數(shù)S(特征值)；

(6)按上述步驟計(jì)算出待檢信號x(t)的特征值Sx，并求出其與標(biāo)準(zhǔn)特征值之間的M距離：

(9)

(7)比較di的大小，其中di的最小值所對應(yīng)的狀態(tài)即為被檢測滾動軸承振動信號所處的工況。

4 試驗(yàn)研究

為驗(yàn)證此方法的正確性，在滾動軸承實(shí)驗(yàn)臺上對不同工況(正常狀態(tài)、內(nèi)圈故障狀態(tài)、外圈故障狀態(tài))滾動軸承的振動信號進(jìn)行采集。實(shí)驗(yàn)臺主要由信號發(fā)生系統(tǒng)、信號采集系統(tǒng)和信號分析系統(tǒng)等部分組成。

試驗(yàn)過程中，采樣頻率為 12kHz，電動機(jī)負(fù)荷為0，轉(zhuǎn)速為1 797r/min，采樣時(shí)間1s。滾動軸承內(nèi)、外圈的故障原因?yàn)辄c(diǎn)火花加工所造成單點(diǎn)損傷。采集到的不同工況下滾動軸承加速度信號，如圖1所示。

圖1 不同狀態(tài)下滾動軸承加速度信號

運(yùn)用改進(jìn)EMD對上述信號進(jìn)行分解，其中，正常狀態(tài)下滾動軸承信號分解得到5階IMF分量，內(nèi)圈和外圈故障信號經(jīng)改進(jìn)EMD分解都得到6階IMF分量，其中不同狀態(tài)下滾動軸承信號的各IMF分量與其原信號間的相關(guān)系數(shù)如表1所示。

表1 不同狀態(tài)下各IMF分量與原信號相關(guān)系數(shù)

由表1可知，滾動軸承的主要信息集中在其前4階IMF分量中，所以，對不同工況的滾動軸承振動信號，在此取IMF1～I(xiàn)MF4作為敏感IMF進(jìn)行后續(xù)分析。不同故障狀態(tài)下敏感IMF分量所求得的局部Hilbert譜示于圖2。

圖2 不同故障狀態(tài)下滾動軸承加速度信號的局部Hilbert譜

由圖2可知，當(dāng)6205-2RS JEM SKF型滾動軸承出現(xiàn)故障時(shí)，激起的固有頻率在2 000 Hz-3 000 Hz，因此，取ω1=2 000 Hz，ω2=3 000 Hz，求出的各狀態(tài)下相應(yīng)的標(biāo)準(zhǔn)特征值和方差var(Si)如表2所示。

表2 標(biāo)準(zhǔn)特征值和其方差

按照馬氏距離判別方法的步驟對不同工況下滾動軸承的60組振動信號進(jìn)行故障診斷研究，診斷結(jié)果表明：在正常狀態(tài)和外圈故障狀態(tài)下各出現(xiàn)1組誤判，內(nèi)圈故障狀態(tài)下診斷結(jié)果全部正確，總體診斷正確率為96.67%，證明該方法在處理此類問題方面的有效性。而出現(xiàn)誤判的原因可能是測量過程中的測量誤差所造成的，部分診斷結(jié)果如表3所示。

表3 部分馬氏距離函數(shù)診斷結(jié)果

5 結(jié) 論

(1)相較于傳統(tǒng)HHT方法，將端點(diǎn)優(yōu)化對稱延拓的EMD方法同基于相關(guān)性分析的虛假IMF判別算法相結(jié)合的改進(jìn)HHT方法可有效解決EMD分解過程中的端點(diǎn)效應(yīng)和虛假模態(tài)問題，保證信號特征提取的準(zhǔn)確性與有效性；

(2)通過M距離判別方法對局部Hilbert邊際能量譜所提取的各IMF分量的能量特征參數(shù)進(jìn)行滾動軸承的工作狀態(tài)和故障類型識別，試驗(yàn)結(jié)果表明，該方法可有效識別出滾動軸承的故障類型和工作狀態(tài)，為該類問題的解決提供了一種有效的方法。

[1] 褚福磊，彭志科，馮志鵬，等.機(jī)械故障診斷中的現(xiàn)代信號處理方法[M].北京：科學(xué)技術(shù)出版社,2009:159-160.

[2] 葛明濤，董素鴿.LMD及馬氏距離敏感閾值的滾動軸承故障診斷[J].機(jī)械設(shè)計(jì)與制造,2015,2:210-213.

[3] 周小龍.希爾伯特-黃變換在故障診斷中的應(yīng)用[D].長春：長春工業(yè)大學(xué),2014.

[4] N.E.Huang,Z.Shen,S.R.Long,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.

[5] 游子躍，王寧，李明明，等.基于EEMD和BP神經(jīng)網(wǎng)絡(luò)的風(fēng)機(jī)齒輪箱故障診斷方法[J].東北電力大學(xué)學(xué)報(bào),2015,35(1):64-72.

[6] 周小龍，姜振海，馬風(fēng)雷.基于改進(jìn)HHT的微弱故障信號特征提取方法[J].東北電力大學(xué)學(xué)報(bào),2016,36(5):52-56.

[7] 王傳菲，安鋼，王凱，等.基于鏡像延拓和神經(jīng)網(wǎng)絡(luò)的EMD端點(diǎn)效應(yīng)改進(jìn)方法 [J].裝甲兵工程學(xué)院學(xué)報(bào),2010,24(2):62-65.

[8] 雷亞國.基于改進(jìn)Hilbert-Huang 變換的機(jī)械故障診斷[J].機(jī)械工程學(xué)報(bào),2011,47(5):71-77.

[9] 黃誠惕.希爾伯特-黃變換及其應(yīng)用研究[D].成都：西南交通大學(xué),2006.

[10] 徐斌，徐德城，朱衛(wèi)平，等.希爾伯特-黃變換方法的改進(jìn) [J].西北工業(yè)大學(xué)學(xué)報(bào),2011,29(2):268-272.

[11] 于德介，程軍圣，楊宇.機(jī)械故障診斷的Hilbert-Huang變換方法[M].北京：科學(xué)技術(shù)出版社,2006:154-166.

[12] 楊德友，蔡國偉.基于因散經(jīng)驗(yàn)?zāi)Ｊ椒纸馀c最小二乘支持向量機(jī)的風(fēng)電場短期風(fēng)速預(yù)測 [J].東北電力大學(xué)學(xué)報(bào),2015,35(3):44-49.

A Study of Rolling Bearing Fault Diagnosis Based on Local Hilbert Marginal Energy spectrum

Yang Gongyong1,Zhou Xiaolong1,Liang Xiuxia2,Li Jiafei2

(1.The Engineering Training Teaching Center,Northeast Electric Power University,Jilin Jilin 132012;2.Henan Xinyu Petroleum Machinery Manufacturing Company,Henan Puyang 457001)

Aiming at the non-linear feature of the rolling bearing vibration signal,a rolling bearing fault diagnosis method based on local Hilbert marginal energy spectrum and Mahalanobis distance is proposed.Firstly,the vibration signal is decomposed by improved empirical mode decomposition,and a set of the intrinsic mode function is obtained and the mode functions which containing the main fault information are selected by the false evaluation method based on correlation analysis technique.Then,energy characteristic parameters of fault signal are extracted from each IMF component by using local Hilbert marginal energy spectrum.Finally,Mahalanobis distance is used to identify the rolling bearing fault pattern.The experiment results show that this method can identify rolling bearing fault patterns effectively and offer a practical method for its fault diagnosis.

Empirical mode decomposition;Local Hilbert marginal energy spectrum;Mahalanobis distance;Rolling bearing;Fault diagnosis

2017-03-12

楊恭勇(1987-)，男，碩士，助理實(shí)驗(yàn)師，主要研究方向：機(jī)械制造及故障診斷.

1005-2992(2017)02-0077-05

TH17

A

電子郵箱： 76025858@qq.com(楊恭勇)；196389679@qq.com(周小龍)；2598037262@qq.com(梁秀霞)；2241430075@qq.com(李家飛)