周 浩，賈民平

(東南大學(xué)機械工程學(xué)院，江蘇南京211189)

0 引 言

據(jù)統(tǒng)計，旋轉(zhuǎn)機械的故障有30%是由滾動軸承故障引起的，它的好壞對機器的工作狀況影響極大。當(dāng)滾動軸承局部存在損傷或缺陷時，其振動信號多表現(xiàn)為非平穩(wěn)、非線性、調(diào)制等特征［1］。利用傳統(tǒng)單一的時域參數(shù)分析法或者頻域分析法等方法很難從該類振動信號中提取故障特征信息。

EMD是將信號分解為有限個高頻到低頻的IMF的和，具有很多獨特優(yōu)點:基函數(shù)的自動產(chǎn)生;自適應(yīng)的濾波特性;自適應(yīng)的多分辨率;正交性和完備性;IMF分量的調(diào)制特性等。上述優(yōu)點使得EMD已經(jīng)被應(yīng)用到許多領(lǐng)域，它非常適合對非線性、非平穩(wěn)信號進行分析［2］。峭度系數(shù)是一個重要的數(shù)值統(tǒng)計量，它可以用來描述信號圖形頂峰的凸平度。當(dāng)峭度為3時，分布曲線具有正常峰度(即零峭度);當(dāng)峭度大于3時，分布曲線具有正峭度，表示信號中沖擊成分較多。同時，Hilbert包絡(luò)解調(diào)通過對故障軸承的振動信號進行幅值和頻率解調(diào)，放大和分離了故障特征信息，極大地提高了信噪比［3］。

基于此，本研究提出一種基于EMD和峭度指標(biāo)的Hilbert包絡(luò)解調(diào)分析的信號處理方法，通過對仿真信號和實際信號分析，表明該方法可以應(yīng)用于滾動軸承的故障診斷。

1 EMD方法和Hilbert包絡(luò)解調(diào)

1.1 EMD基本理論

經(jīng)驗?zāi)B(tài)分解法從本質(zhì)上來講就是通過對非平穩(wěn)信號的分解，來獲得一系列表征信號特征時間尺度本征模函數(shù)的一種信號平穩(wěn)化處理方法［4-5］。其中，每個IMF分量必須滿足如下2個條件:

(1)本征模函數(shù)在整個時間范圍內(nèi)，局部極值點和過零點的個數(shù)必須相等，或最多相差1個;

(2)在時間范圍內(nèi)的任意時刻，由局部極大值點形成的上包絡(luò)線和由局部極小值點形成的下包絡(luò)線的平均值為零，即上、下包絡(luò)線關(guān)于時間軸呈局部對稱。

EMD方法的分解過程如下:

(1)找出原始信號x(t)所有的極大值點和極小值點，并用三次樣條插值函數(shù)分別擬合這些數(shù)據(jù)點，形成上、下包絡(luò)線。

(2)計算上包絡(luò)線和下包絡(luò)線的均值m1，將原始信號減去該平均包絡(luò)得到一個新的數(shù)據(jù)序列h1:

如果h1全部滿足本征模函數(shù)的兩個條件，則h1就是原始信號x(t)分解出來的第一個IMF分量c1。否則說明h1還不是一個本征模函數(shù)，還需將h1作為一個原始信號重復(fù)步驟(1)，繼續(xù)進行篩選，直至產(chǎn)生第一個IMF分量c1。

(3)計算:

再將r1作為原始信號重復(fù)步驟(1)～(3)，以該類推得原始信號x(t)的n個滿足條件的分量和一個不能再分解的單調(diào)函數(shù)，即殘余分量rn為止［6］，可得:

式中:rn—信號的平均趨勢。

1.2 Hilbert包絡(luò)解調(diào)

利用希爾伯特變換進行信號包絡(luò)的原理是讓測試信號產(chǎn)生一個90°的相移，從而與原信號構(gòu)成一個解析信號，該解析信號即構(gòu)成包絡(luò)信號［7］。

一個實信號x(t)的希爾伯特變換定義為:

希爾伯特變換的包絡(luò)A(t)定義為:

瞬時相位φ(t)定義為:

經(jīng)過希爾伯特變換檢波，除去高頻的振動分量，然后用含有缺陷激振分量的A(t)代替原始信號x(t)進行頻譜分析，便可以進行缺陷的故障診斷［8］。

2 診斷方法

2.1 峭度指標(biāo)

四階中心矩也是一個重要的統(tǒng)計估計量。它可以用來表示樣本的函數(shù)圖形頂峰的凸平度，即峭度系數(shù)，其數(shù)學(xué)描述為:

式中:K—原始信號x的峭度指標(biāo);N—信號的長度;μ—信號x的均值;σ—信號x的標(biāo)準差。

當(dāng)峭度值為3時，分布曲線具有正常峰度(即零峭度)。當(dāng)標(biāo)準差小于正常狀態(tài)下的標(biāo)準差，即觀測值的分散程度較小沖擊成分較多時，峭度值增大，此時曲線峰頂?shù)母叨雀哂谡Ｕ龖B(tài)分布曲線，稱為正峭度。當(dāng)標(biāo)準差大于正常狀態(tài)下的標(biāo)準差，即觀測值的分散程度較大時，峭度系數(shù)減小，此時曲線峰頂?shù)母叨鹊陀谡Ｕ龖B(tài)分布曲線，稱為負峭度。

2.2 融合EMD、峭度和Hilbert解調(diào)的診斷方法

滾動軸承故障振動信號中沖擊成分所占比重較多，頻率成分很豐富，在不同頻帶內(nèi)所包含的軸承故障信息不同。當(dāng)研究者分析滾動軸承振動信號時，經(jīng)過EMD分解的各階IMF分量會包含由故障引起的不同頻段的固有振動成分［9-10］。依據(jù)峭度系數(shù)理論，峭度值較大的IMF分量中周期性沖擊成分較多，即包含滾動軸承的故障信息越多［11］。因此，本研究選取原始振動信號經(jīng)過EMD分解后峭度值較大幾階IMF分量重構(gòu)信號，并將其作為Hilbert包絡(luò)解調(diào)的對象，提取故障特征頻率，與滾動軸承各個部件的理論故障頻率對比，找出故障部位。因此，本研究采用下述分析步驟對滾動軸承進行故障診斷:

(1)對滾動軸承故障信號進行EMD分解，得到多個頻率從高到低排列的IMF分量;

(2)根據(jù)式(10)計算各階IMF分量的峭度系數(shù);

(4)對重構(gòu)信號進行Hilbert包絡(luò)解調(diào)分析，得到包絡(luò)信號;

(5)對包絡(luò)信號進行FFT，獲得信號的包絡(luò)譜;

(6)參照滾動軸承各個部件的理論故障特征頻率，對包絡(luò)譜進行分析，根據(jù)譜上特征頻率來判斷軸承故障類型。

3 仿真信號分析

滾動軸承運行時，當(dāng)外圈、內(nèi)圈或者滾動體上發(fā)生磨損、點蝕、剝落等局部損傷時會產(chǎn)生突變的沖擊脈沖力，該沖擊作用會誘發(fā)軸承系統(tǒng)的高頻固有振動成分，如軸承內(nèi)外圈的徑向彎曲固有振動、滾動體的固有振動或者測振傳感器的固有振動等，導(dǎo)致在時域表現(xiàn)為簡潔的周期性沖擊信號被各種高頻信號調(diào)制，形成較為復(fù)雜的調(diào)制信號。因此，本研究對該類信號進行包絡(luò)解調(diào)處理，分析其特征頻率將能夠準確地辨識出滾動軸承的故障。

可以用如下數(shù)學(xué)模型描述滾動軸承內(nèi)圈點蝕故障:

式中:fi—滾動軸承的內(nèi)圈故障特征頻率;fz—載波頻率;Bcos(2πfi)—測試信號;A—一個直流分量，使偏置后的信號［A+Bcos(2πfi)］都具有正值;cos(2πfzt)—載波信號。

仿真程序隨機生成滾動軸承內(nèi)圈的故障特征頻率和對應(yīng)的高頻載波頻率，考慮實際工況下振動信號含有噪聲，此處添加高斯白噪聲信號n(t)。仿真形成滾動軸承內(nèi)圈故障的振動信號為:

采用SPSS 21.00統(tǒng)計學(xué)軟件對數(shù)據(jù)進行處理，計量資料以“±s”表示，采用t檢驗（血液流變學(xué)指標(biāo)），計數(shù)資料以百分數(shù)（%）表示，采用x2檢驗（患者滿意度），以P＜0.05為差異有統(tǒng)計學(xué)意義。

本研究設(shè)置采樣頻率為12 000 Hz，取持續(xù)1 s的仿真數(shù)據(jù)分析，對應(yīng)的時域波形和頻域波形如圖1所示。由圖1可以看出存在邊頻信號，雖然兩者相減就是特征信號頻率，但是在實際應(yīng)用中，振動信號包含大量的干擾信號，其頻譜圖譜線較多，很難清晰地辨別出轉(zhuǎn)頻及其變頻帶，從而無法準確辨識出滾動軸承故障發(fā)生的部位。

圖1 調(diào)幅仿真信號時域和頻域波形圖

本研究對x(t)進行EMD分解，得到12階IMF分量。由于高斯隨機白噪聲的存在，仿真信號每次模擬時都會變動，導(dǎo)致各階IMF分量的峭度系數(shù)會發(fā)生改變。但經(jīng)過多次實驗，發(fā)現(xiàn)前4階高頻IMF分量的峭度系數(shù)一直較大，現(xiàn)取前4階IMF分量重構(gòu)原信號，并對其進行Hilbert包絡(luò)解調(diào)分析，解調(diào)后的譜圖如圖2所示。

圖2 經(jīng)Hilbert包絡(luò)解調(diào)的重構(gòu)信號頻譜圖

由圖2可知，含有噪聲的仿真信號經(jīng)過EMD分解、選擇性重組和Hilbert包絡(luò)解調(diào)后，可以清晰地看出其頻譜圖有一條非常突出的譜線，且該譜線對應(yīng)的頻率剛好是滾動軸承內(nèi)圈故障特征頻率105 Hz。可見，融合EMD分解、峭度指標(biāo)和Hilbert包絡(luò)解調(diào)的分析方法能夠很好地提取滾動軸承的故障特征頻率，具有一定的實際應(yīng)用價值。

4 滾動軸承故障診斷實例分析

本次使用的故障數(shù)據(jù)取自美國Western Reserve University電氣工程實驗室的滾動軸承故障模擬實驗臺的滾動軸承數(shù)據(jù)。該實驗臺由一個1.5 kW的電動機，一個扭矩傳感器/譯碼器，一個功率測試計，還有電子控制器構(gòu)成。待檢測的軸承支撐著電動機的轉(zhuǎn)軸，驅(qū)動端軸承為SKF6205深溝球軸承。軸承的局部損傷為電火花加工的單點損傷。采樣頻率為12 000 Hz，軸承轉(zhuǎn)速為r=1 750 r/min［12］。

查閱手冊可得SKF6205深溝球軸承的標(biāo)準參數(shù)，通過理論公式計算可得軸承內(nèi)圈的理論故障特征頻率為157.9 Hz。驅(qū)動端滾動軸承內(nèi)圈故障的時域波形如圖3所示，采樣點數(shù)12 000。

圖3 滾動軸承內(nèi)圈故障的時域波形

本研究利用融合EMD分解、峭度指標(biāo)和Hilbert包絡(luò)解調(diào)的分析方法對原始振動信號進行處理。EMD分解后得到的12階IMF分量c1～c12和殘余分量res如圖4所示。

圖4 軸承故障振動信號的第1～12階IMF和res分量

各階IMF分量的峭度值如表1所示。

表1 各階IMF峭度值

從表1中可以看出，前4階IMF分量的峭度值較大，相對而言其包含的故障信息最大，對前4階IMF分量進行疊加重構(gòu)原信號，并對重構(gòu)信號進行Hilbert包絡(luò)解調(diào)處理，其包絡(luò)譜如圖5所示。從圖5中可以較為清晰地看出，包絡(luò)譜中出現(xiàn)了滾動軸承內(nèi)圈故障特征頻率的1倍頻(158 Hz)、2倍頻(316 Hz)、3倍頻(473 Hz)等，其中譜線在頻率29 Hz處為軸承的轉(zhuǎn)頻。從上述分析可以確認該滾動軸承內(nèi)圈發(fā)生故障，與實際相符。因此，本研究所提出的基于EMD和峭度指標(biāo)的Hilbert包絡(luò)解調(diào)分析方法能夠應(yīng)用在滾動軸承故障中。

圖5 重構(gòu)的包絡(luò)譜(局部放大圖)

5 結(jié)束語

EMD分解能將非線性非平穩(wěn)振動信號分解為多個能夠反映原始信號局部特征的固有模態(tài)函數(shù)IMF，通過對部分IMF分量重構(gòu)分析，能夠清楚地揭示信號中包含的故障信息。Hilbert包絡(luò)解調(diào)能有效地將故障特征信息從復(fù)雜的調(diào)制信號中提取出來，有效地排除了各種高頻干擾、提高信噪比。

本研究提出了基于EMD和峭度指標(biāo)的Hilbert包絡(luò)解調(diào)分析的信號處理方法，并將其應(yīng)用于滾動軸承故障診斷。筆者通過對仿真信號與實測信號的應(yīng)用，清晰而準確地得到了滾動軸承的故障特征信息，表明了該方法可以較好地用于滾動軸承故障信號處理，為滾動軸承的故障診斷提供了一種有效的方法。

［1］鐘秉林，黃 仁.機械故障診斷學(xué)［M］.3版.北京:機械工業(yè)出版社，2007.

［2］李 力.機械信號處理及其應(yīng)用［M］.武漢:華中科技大學(xué)出版社，2007.

［3］胡清忠，褚福磊.基于Hilbert解調(diào)的滾動軸承故障診斷分析［J］.振動與沖擊，2010，29(S):1-3.

［4］焦衛(wèi)東，朱有劍.基于EMD與倒譜分析的軸承故障診斷［J］.機電工程，2009，26(2):18-21.

［5］張梅軍，曹 勤，王 闖.EMD分解，分形理論和RBF神經(jīng)網(wǎng)絡(luò)相結(jié)合的軸承智能故障診斷研究［J］.機械，2012(11):16-21.

［6］HUANG N E，SHEN Z，LONG S R，et al.The empirical mode decomposition and the Hilbert spectrum for non-linear non-stationary time series analysis［J］.Proceedings of the Royal Society，1998(454):903-995.

［7］宋曉美，孟繁超，張 玉.基于包絡(luò)解調(diào)分析的滾動軸承故障診斷研究［J］.儀表儀器與分析監(jiān)測，2012(1):1-4.

［8］LI Hui，ZHANG Yu-ping，ZHENG Hai-qi.Hilbert-Huang transform and maginal spectrum for detection and diagnosis of localized defects in roller bearing［J］.Journal of Mechanical Science and Technology，2009，23(2):291-301.

［9］朱 瑜，王 殿，王海洋.基于EMD和信息熵的滾動軸承故障診斷［J］.軸承，2012，15(6):50-53.

［10］施國友.滾動軸承對螺桿壓縮機運行的影響［J］.流體機械，2012，40(7):35-36.

［11］YANG YU.Roller bearing fault diagnosis method based on EMD［J］.China Mechanical Engineering，2004，15(10):908-920.

［12］The Case Westren Reserve University Bearing Data Center Website.Bearing data center seeded fault test data［EB/OL］.［2007-01-01］.http://www.eecs.case.edu/laboratory/bearing.