變分模態(tài)分解和改進的自適應(yīng)共振技術(shù)在軸承故障特征提取中的應(yīng)用

李華 伍星 劉韜

摘要： 針對滾動軸承早期故障特征提取困難的問題，提出了基于變分模態(tài)分解（Variational Mode Decomposition， VMD）和改進的自適應(yīng)共振技術(shù)的滾動軸承故障特征提取方法。針對軸承故障信號所在頻帶難以選擇的問題，提出了基于改進的自適應(yīng)共振技術(shù)（Improved Adaptive Resonance Technology， IART）的IMF選取方法。首先，確定模態(tài)數(shù)，提出了峭度最大值的模態(tài)數(shù)確定方法；然后，對原始振動信號進行VMD分解，獲得既定數(shù)目的本征模態(tài)分量（Intrinsic Mode Function， IMF）；其次，利用IART選取包含豐富故障信息的IMF分量；最后，（如有需要）對選取的IMF分量進行基于IART的帶通濾波，并進行包絡(luò)解調(diào)分析提取故障特征頻率。將該方法應(yīng)用到軸承仿真數(shù)據(jù)和實際數(shù)據(jù)中，能夠?qū)崿F(xiàn)軸承故障特征的精確診斷，證明了該方法的有效性。

關(guān)鍵詞： 故障診斷； 滾動軸承； 變分模態(tài)分解； 峭度； 改進的自適應(yīng)共振

中圖分類號：TH165+.3； TN911.7文獻標志碼： A文章編號： 1004-4523（2018）04-0718-09

DOI：10.16385/j.cnki.issn.1004-4523.2018.04.020

引言

對軸承的狀態(tài)監(jiān)測與故障診斷一直是機械設(shè)備故障診斷的熱點[1]。當滾動軸承發(fā)生故障時，其振動信號包含了大量的運行狀態(tài)信息，表現(xiàn)為非平穩(wěn)性和多分量性的調(diào)制信號，特別在故障早期，由于調(diào)制源弱，早期故障信號微弱，并且受周圍設(shè)備、環(huán)境的噪聲干擾，導(dǎo)致故障特征頻率難以提取、識別[2]。

故障診斷的關(guān)鍵是從原始信號中提取故障特征信號（故障特征頻率）。Huang等[3]提出了經(jīng)驗?zāi)Ｊ椒纸?，此方法具有自適應(yīng)分解特性，對非平穩(wěn)和非線性信號的處理具有較高的效率。因此，基于EMD的故障診斷方法層出不窮。但EMD存在模態(tài)混疊、端點效應(yīng)、受采樣頻率影響較大等不足[4-5]?；诖?，Dragomiretskiy等[6]提出了一種自適應(yīng)信號處理新方法——變分模態(tài)分解（VMD），該方法通過迭代搜尋變分模型最優(yōu)解來確定每個分量的頻率中心及帶寬，相比于EMD，VMD將信號分解轉(zhuǎn)化為非遞歸、變分模態(tài)分解模態(tài)并且具有堅實的理論基礎(chǔ)，有效地彌補了EMD和EEMD等的不足。很多學(xué)者已將VMD應(yīng)用于軸承故障診斷研究。Mohanty等[7]將VMD引入軸承故障診斷，并取得了比EMD更好的效果。劉長良等[8]提出了基于變分模態(tài)分解的故障特征提取方法，利用觀察中心頻率的方法確定模態(tài)K，并采用模糊C均值聚類進行故障識別。趙洪山等[9]提出了基于變分模態(tài)分解和奇異值分解的故障特征提取方法，利用峭度指標選取敏感IMF分量，對選取的IMF進行奇異值分解重構(gòu)，取得良好效果。劉尚坤等[10]以互信息為判據(jù)對原VMD方法進行了迭代停止條件的改進，并結(jié)合Teager能量算子的優(yōu)點，提出了Teager-VMD時頻分析新方法。針對故障滾動軸承振動信號的非平穩(wěn)和非高斯特性，Te Han等[11]提出了一種基于變分模態(tài)分解（VMD）和自回歸（AR）模型參數(shù)的故障特征提取方法。并且在滾動軸承故障診斷領(lǐng)域提出了一種新型隨機森林分類器用于模式識別。夏均忠等[12]提出了一種用最大相關(guān)峭度解卷積（Maximum Correlated Kurtosis Deconvolution， MCKD）和變分模態(tài)分解相結(jié)合的方法提取滾動軸承故障特征。Zipeng Li等[13]針對背景噪聲嚴重影響模式個數(shù)的設(shè)置，提出了一種基于相關(guān)性分析的面向獨立的VMD方法，用于自適應(yīng)提取輪對軸承的弱故障和復(fù)合故障特征。

當滾動軸承發(fā)生故障時，由于振動信號表現(xiàn)出的幅值調(diào)制特性，通過包絡(luò)分析可以得到軸承的故障特征頻率[14]。鑒于VMD的模態(tài)數(shù)需事先確定的問題及其敏感IMF分量的選取問題，本文提出VMD和改進的自適應(yīng)共振技術(shù)相結(jié)合的診斷方法。

1基礎(chǔ)理論

1.1變分模態(tài)分解（VMD）變分模態(tài)分解（VMD）[6]是一種基于維納濾波的自適應(yīng)信號處理新方法。在VMD算法中，每一個IMF均被稱為AM—FM信號，即ukt=Aktcoskt（1）式中ωkt=′kt=dktdt，Akt為信號的瞬時幅值，ωkt為瞬時頻率。在t-δ，t+δδ=2π'kt內(nèi)，可將ukt視為幅值為Akt，頻率為ωkt的諧波信號。

第4期李華，等：變分模態(tài)分解和改進的自適應(yīng)共振技術(shù)在軸承故障特征提取中的應(yīng)用振 動 工 程 學(xué) 報第31卷VMD的求解過程是一個變分問題的求解過程。假設(shè)每個模態(tài)是具有中心頻率的有限帶寬，中心頻率和帶寬在分解過程中不斷更新，則，變分問題可表示為尋求k個模態(tài)函數(shù)ukt，并使得所有模態(tài)函數(shù)的估計帶寬和最小，模態(tài)之和為輸入信號f。具體步驟如下：

1.2峭度

峭度（Kurtosis）反映振動信號分布特性的數(shù)值統(tǒng)計量，是描述波形尖峰度的無量綱參數(shù)[15]。

其數(shù)學(xué)描述為K=Ex-μ4σ4（6）式中x為所分析的振動信號，μ為信號x的均值，σ為信號x的標準差。

峭度對沖擊信號敏感，適合軸承早期故障診斷。滾動軸承無故障時，峭度指標值約等于3。當滾動軸承發(fā)生局部故障時，峭度值會增大。

1.3改進的自適應(yīng)共振技術(shù)（IART）

由于傳統(tǒng)共振解調(diào)技術(shù)的缺點，共振頻率和帶通濾波器的中心頻率都是固定的，王平等[16]提出了一種自適應(yīng)諧振解調(diào)技術(shù)。在此基礎(chǔ)上，本文提出了一種基于功率譜的共振頻率選取方法：直接對信號進行功率譜分析，依靠功率譜分析來自動識別原始振動信號和各IMF的共振頻率（功率譜反映了能量的集中情況）。并在此基礎(chǔ)上確定帶寬，設(shè)計帶通濾波器。由于實際情況下共振頻率一般在故障頻率的10倍以上，因此設(shè)定一個閾值，當振幅峰值出現(xiàn)在故障頻率的10倍以下時，將被放棄。

選取共振頻率存在如下兩種情況：1）當滿足閾值條件的功率譜中只有一個幅度峰值時，可以直接選擇為中心頻率，即共振頻率；2）當滿足閾值條件的功率譜中存在多條譜線時，選擇幅值最大的功率譜線作為中心頻率，即共振頻率。

然后確定帶寬，這里參考文獻[12]的設(shè)定方法，將帶通濾波器帶寬設(shè)定為3倍故障頻率。

2基于VMD和IART的故障特征提取針對滾動軸承早期故障診斷中故障特征微弱難以提取的問題，本文提出VMD和IART相結(jié)合的方法。然后，對選取的IMF進行基于IART的帶通濾波，進一步提高信噪比，并對濾波后的信號進行包絡(luò)解調(diào)分析，提取故障特征頻率。

2.1基于峭度最大值的模態(tài)優(yōu)化

因為峭度指標是故障的一種表征，當存在某個模態(tài)的峭度值具有最大值時，可認為此時具有良好的分解效果（在給定的搜索范圍內(nèi)）?；诖耍疚奶岢隽嘶谇投茸畲笾祪?yōu)化模態(tài)數(shù)K的方法。具體步驟及流程在2.3節(jié)中給出。

2.2基于IART選取敏感IMF分量

軸承振動信號經(jīng)過VMD分解為一系列IMF分量，若要進行故障特征提取，最好的方法是選取包含故障信息最多的IMF分量進行分析，為此，本文提出了基于IART的敏感IMF選取方法。

對信號進行功率譜分析，選取共振頻率，并以此為中心頻率設(shè)計帶通濾波器，將此頻帶范圍作為特征頻帶區(qū)間，利用區(qū)間從屬關(guān)系選取有效IMF分量。用雙值特征區(qū)間系數(shù)代替一般的單值特征參考系數(shù)，有效地避免了單值參考系數(shù)的不易選取的問題。

首先，對原信號和VMD分解獲得的若干IMF分量分別求取功率譜，選取共振頻率；然后，以此共振頻率作為中心頻率設(shè)計帶通濾波器，并將此濾波頻帶作為其特征頻帶區(qū)間；其次，比較各個IMF分量的特征頻帶區(qū)間與原信號特征頻帶區(qū)間的從屬關(guān)系，若滿足，則選取為敏感IMF，否則舍棄；如果存在多個敏感IMF，則對其進行疊加、重構(gòu)。

2.3算法步驟及流程

基于VMD和IART的故障特征提取方法的具體步驟及診斷流程圖如圖1所示。

（1）獲取原始振動信號，初始化模態(tài)數(shù)K=2，懲罰因子α和帶寬τ使用默認值：α=2000，τ=0；

（2）對軸承原始振動信號進行VMD分解，計算各個模態(tài)的峭度值，并通過比較獲得此模態(tài)數(shù)下的峭度最大值，然后K=K+1繼續(xù)以上分析，直到取到K=16為止；比較在各個模態(tài)數(shù)下所取得的峭度最大值的大小，將最大的峭度值所對應(yīng)的模態(tài)數(shù)K判定為最優(yōu)模態(tài)數(shù)K。為了確定模態(tài)數(shù)K的搜索范圍和步長，本文借鑒文獻[17]給出的模態(tài)數(shù)K的討論范圍，由于K過大，造成效率低，計算負擔(dān)重；K過小，易引入噪聲，所以搜索范圍設(shè)為[2，16]。

（3）利用IART選取敏感IMF分量。計算原始振動信號與各個IMF分量的功率譜，獲得共振頻率，并以此共振頻率作為中心頻率設(shè)計帶通濾波器。選取濾波頻帶從屬于原始信號濾波頻帶的IMF分量為敏感IMF進行信號重構(gòu)；

（4）對重構(gòu)信號進行基于IART的帶通濾波，此處的帶通濾波器由步驟（3）獲得，進一步提高SNR；

（5）對步驟（4）中濾波后的信號進行包絡(luò)解調(diào)分析，判斷故障發(fā)生部位。

3信號仿真分析

對原信號進行11個模態(tài)的VMD分解，獲得11個IMF分量，由于篇幅限制，僅給出前6個IMF分量的時域波形，如圖4所示。由于時域圖差別不大，無法直觀地選擇敏感IMF分量進行分析。因此，需要合適的方法選取敏感IMF分量。

3.2有效IMF選取及帶通濾波

利用本文提出的IART選取有效的IMF分量，對原始信號和各IMF分量求取功率譜，按照第1.3節(jié)獲得共振頻率，以此作為中心頻率設(shè)計帶通濾波器，獲得頻帶區(qū)間。同樣地，由于篇幅限制，只給出了原信號和IMF5-IMF7的功率譜圖，如圖5所示。

由圖5可知，原信號的共振頻率選取為3000 Hz。以此作為中心頻率設(shè)計帶通濾波器，根據(jù)文獻[12]，取濾波帶寬為3倍的故障特征頻率，即Bw=3fin=237 Hz，進一步獲得濾波頻帶為[2882，3119]，將其作為原信號的特征頻帶區(qū)間。從各IMF的功率譜可知，IMF6的共振頻率在3000 Hz。同樣地，設(shè)計帶通濾波器，Bw=3fin=237 Hz，其特征頻帶區(qū)間為[2882，3119]?？芍?，IMF6的濾波頻帶從屬于原信號的濾波頻帶，故選取IMF6為敏感IMF分量進行分析。原信號和IMF5-IMF7的特征頻帶區(qū)間如表1所示。

對IMF6進行包絡(luò)解調(diào)分析，如圖6所示。從圖中可以清晰地提取的故障特征頻率（78.52 Hz）、轉(zhuǎn)頻（28.13 Hz）以及二倍頻成分，均與理論值很接近，由此可判別故障類型為內(nèi)圈故障。這充分說明了本文提出的方法的有效性。而且，此時已經(jīng)幾乎沒有噪聲的影響，只有少許微弱噪聲毛刺，所以最后一步的帶通濾波可以舍棄（此處仍給出了帶通濾波的包絡(luò)譜，如圖7所示。由圖可知，殘余的微弱噪聲被有效的去除，效果良好）。

4實驗驗證及分析

為了證明本文提出方法的有效性及優(yōu)越性，將其用于實際軸承數(shù)據(jù)分析。本文采用來源于美國西儲大學(xué)的軸承數(shù)據(jù)[19]，軸承的型號為6205RS JEM SKF，采樣頻率fs=12000 Hz，試驗數(shù)據(jù)選擇轉(zhuǎn)速為1750 r/min，負載為2 hp（1492 W），故障尺寸為0.021″，驅(qū)動端軸承在滾動體故障狀態(tài)下的數(shù)據(jù)。

4.1VMD分析

如圖8（a）和（b）分別是軸承滾動體故障數(shù)據(jù)的原始信號時域圖和包絡(luò)譜。時域圖中由于噪聲的嚴重影響，無法提取沖擊特征；包絡(luò)譜中，雖然能夠提取故障特征頻率，但受到很嚴重的寬頻帶噪聲影響，如圖中 2000～3000 Hz范圍內(nèi)的波峰尤為突出（尤其紅色橢圓框內(nèi)），無法判別。因此，需要對信號進行預(yù)處理，提高信噪比。

為此，對原始振動信號進行VMD分解，初始化模態(tài)數(shù)K=2，設(shè)定K的搜索范圍為[2，16]（由文獻[13]可知，K過大，造成效率低，計算負擔(dān)重；K過小，易引入噪聲），利用本文提出的峭度最大值原則確定最優(yōu)模態(tài)K。經(jīng)過分析，當K=6時，取得峭度最大值為3.4620，所以取最優(yōu)模態(tài)數(shù)為K=6。模態(tài)數(shù)與峭度最大值關(guān)系圖如圖9所示。

對原信號進行6個模態(tài)的VMD分解，獲得6個IMF分量，各個IMF分量的時域波形如圖9所示。同樣地，由于時域圖各個IMF分量差別不大，無法直觀地選擇敏感IMF分量進行分析。因此，利用本文提出的方法選取敏感IMF分量。

4.2有效IMF選取及帶通濾波

利用本文提出的基于IART的方法選取有效IMF分量。由于篇幅限制，只給出了原信號和IMF3-IMF6的功率譜圖，如圖11所示。

如圖11（a）所示，根據(jù)1.3節(jié)所述，可取原始振動信號的共振頻率為3345 Hz。以此作為中心頻率設(shè)計帶通濾波器，按文獻[12]取濾波帶寬為3倍的故障特征頻率，即Bw=3fb=412 Hz，進一步獲得濾波頻帶為[3139，3551]，并將其作為原信號的特征頻帶區(qū)間。

然后，獲取各個IMF分量的共振頻率，并設(shè)計帶通濾波器，獲得特征頻帶區(qū)間。經(jīng)分析，可知IMF5的共振頻率為3345 Hz，與原始振動信號的一致。同樣地，設(shè)計帶通濾波器，Bw=3fin=412 Hz，其特征頻帶區(qū)間為[3139，3551]?？芍?，IMF5的濾波頻帶從屬于原信號的濾波頻帶，故選取IMF5為敏感分量進行分析。原信號和IMF4-IMF6的特征頻帶區(qū)間如表3所示。

為了證明本文提出的IMF選取方法的優(yōu)越性，這里利用峭度最大值原則對有效IMF進行選取。對各個IMF分量進行峭度計算，經(jīng)比較，IMF4具有峭度最大值，因此，選取IMF4為有效分量。對其進行包絡(luò)解調(diào)分析。如圖13所示，從圖中雖然能夠提取轉(zhuǎn)頻，但故障特征頻率幾乎無法提取，而且存在很多未知的波峰干擾。顯然，本文圖12具有更好的效果。

從圖12可看到，雖然已經(jīng)能夠比較清晰地提取故障特征頻率，但存在一定的噪聲影響，比如在310 Hz附近的波峰。這里，筆者提出利用選取IMF時設(shè)計的帶通濾波器對選取的IMF進行帶通濾波，進一步提高信號的信噪比。如圖14所示，經(jīng)過濾波后，有效地去除了干擾噪聲（310 Hz附近的波峰被去除）。證明了基于VMD和IART的故障診斷方法的有效性。

為了驗證基于峭度最大值優(yōu)化模態(tài)數(shù)K的優(yōu)勢，取K=11，利用本文所述的IART方法選取的敏感IMF分量（經(jīng)分析，只有IMF7的特征頻帶與原信號的特征頻帶滿足從屬關(guān)系，故選取的敏感IMF分量為IMF7）的包絡(luò)譜如圖15所示。容易看出，雖然可以提取出轉(zhuǎn)頻，但故障頻率幾乎淹沒在噪聲中，由圖12所呈現(xiàn)的效果明顯的更為理想。由此表明，相比于人為主觀決策，基于峭度最大值原則的模態(tài)數(shù)確定方法提供了理論依據(jù)，且更為可靠。

4.3與EEMD對比

采用EEMD算法對上述滾動體故障信號進行分解，得到14個IMF分量，取前3個IMF分量求取頻譜，如圖16所示。由圖可以看到EEMD分解存在頻率混疊，不同分量之間的頻帶差不明顯，且互相影響。

4.4基于VMD和IART的弱故障分析

為了能說明本文方法的適用性以及對微弱故障的處理能力，在這里分析了文[19]中軸承的型號為6205RS JEM SKF，故障尺寸為0.007″，轉(zhuǎn)速為1，730 r/min，負載為3 hp（2238 W），采樣頻率fs=12000 Hz的滾動體故障數(shù)據(jù)。計算可知，轉(zhuǎn)頻fr=28.7 Hz，滾動體故障特征頻率fb=134.89 Hz。

4.4.1VMD分析

如圖19（a）和（b）分別是軸承滾動體弱故障數(shù)據(jù)的原始信號時域圖和包絡(luò)譜。時域圖中由于噪聲的嚴重影響，無法提取沖擊特征；包絡(luò)譜中由于受到很嚴重的噪聲影響，只能提取到轉(zhuǎn)頻，而無法提取故障特征頻率。故利用本文方法度其進行處理，提高信噪比。

利用峭度最大值原則確定最優(yōu)模態(tài)K。經(jīng)分析，當K=11時，具有峭度最大值（3.762），所以最優(yōu)模態(tài)為K=11。對原信號進行11個模態(tài)的VMD分解，獲得11個IMF分量。

4.4.2有效IMF選取

根據(jù)IART選取有效IMF分量，可知，IMF8為有效IMF分量，故對其進行包絡(luò)解調(diào)分析，如圖20（a）所示。同樣地，利用峭度最大值原則選取IMF3為敏感分量，其包絡(luò)譜如圖20（b）。由圖可知，圖（a）可以清晰地提取故障特征頻率（131.3 Hz）和轉(zhuǎn)頻（29.3 Hz），均與理論計算值很接近，而圖20（b）只能提取2倍轉(zhuǎn)頻，無法提取故障特征頻率。由此可知，本文提出的選取方法具有明顯優(yōu)勢。證明了本文方法可以有效的應(yīng)用于微弱特征的提取。由圖20（a）可知，此時已經(jīng)幾乎沒有噪聲影響，所以無需進行二次濾波。

同樣地，利用EEMD對信號進行分析，可得14個IMF分量?；贗ART選取敏感分量， IMF1和IMF2滿足從屬關(guān)系，將其疊加重構(gòu)，包絡(luò)譜如圖21所示。由圖可知，包絡(luò)譜僅能提取轉(zhuǎn)頻，無法提取故障特征頻率。且存在很嚴重的噪聲干擾。故也說明了本文方法具有明顯優(yōu)勢。

5結(jié)論

1）首先，本文針對模態(tài)數(shù)K，提出了峭度最大值原則的優(yōu)化方法，能夠有效地確定最優(yōu)模態(tài)數(shù)K。并且比人為主觀決策更有理論依據(jù)更可靠。

2）針對軸承故障信號所在頻帶難以選擇的問題，提出了基于改進的自適應(yīng)共振技術(shù)（IART）的IMF選取方法。在對軸承滾動體故障的兩個實例分析中，均取得了比峭度最大值原則選取IMF更好的效果。

3）針對VMD分解后，若幾乎不存在噪聲影響，則完成診斷過程；若存在的殘余噪聲對故障特征提取的影響，提出了將VMD與IART相結(jié)合的故障特征提取方法。并在對實際軸承故障的分析中證明了其有效性。

4）VMD還存在缺陷，K，α等需事先確定，而且參數(shù)的搜索范圍建立在前人討論的基礎(chǔ)上，缺乏理論依據(jù)，有待進一步完善。而自適應(yīng)共振的中心頻率及帶寬的選取也需進一步研究。

參考文獻：

[1]Wang H Q， Hou W， Tang G， et al. Fault Detection Enhancement in Rolling Element Bearings via Peak-Based Multiscale Decomposition and Envelope Demodulation[J]. Mathematical Problems in Engineering， 2014， 2014（1）：135—142.

[2]丁康， 黃志東， 林慧斌. 一種譜峭度和Morlet 小波的滾動軸承微弱故障診斷方法[J]. 振動工程學(xué)報， 2014， 27（1）：128—135.

DING Kang， HUANG Zhidong， LIN Huibin. A weak fault diagnosis method for rolling element bearings based on morlet wavelet and spectral kurtosis [J]. Journal of Vibration Engineering， 2014， 27（1）：128—135.

[3]Huang N E， Shen Z， Long S R， et al. The empirical mode decomposition and the Hilbert spectrum for non-linear and non-stationary time series analysis[C]. Proceedings of the Royal Society Lond，1998， 454：903—995.

[4]Wu Zhaohua， 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]Rilling G， Flandrin P. On the influence of sampling on the empirical Mode Decomposition[C]. IEEE International Conference on Acoustics， Speech and Signal Processing. ICASSP 2006 Proceedings，IEEE， 2006：444—447.

[6]Dragomiretskiy K， Zosso D. Variational Mode Decomposition[J]. IEEE Transactions on Signal Processing， 2014， 62（3）：531—544.

[7]Mohanty， Gupta K K， Raju K S. Bearing fault analysis using variational mode decomposition[C].International Conference on Industrial and Information Systems. IEEE， 2014：1—6.

[8]劉長良， 武英杰， 甄成剛. 基于變分模態(tài)分解和模糊C均值聚類的滾動軸承故障診斷[J]. 中國電機工程學(xué)報， 2015， 35（13）：3358—3365.

LIU Changliang， WU Yingjie， ZHENG Chenggang. Rolling bearing fault diagnosis based on variational mode decomposition and fuzzy C means clustering [J]. Proceeding of the CSEE， 2015， 35 （13）： 3358—3365.

[9]趙洪山， 郭雙偉， 高奪. 基于奇異值分解和變分模態(tài)分解的軸承故障特征提取[J]. 振動與沖擊， 2016， 35（22）：183—188.

ZHAO Hongshan， GUO Shuangwei， GAO Duo. Fault feature extraction of bearing faults based on singular value decomposition and variational modal decomposition [J]. Journal of Vibration and Shock， 2016， 35 （22）： 183-188.

[10]劉尚坤， 唐貴基， 王曉龍. 基于改進變分模態(tài)分解的旋轉(zhuǎn)機械故障時頻分析方法[J]. 振動工程學(xué)報， 2016， 29（6）：1119—1126.

LIUShangkun， TANG Guiji， WANG Xiaolong. Time frequency analysis method for rotary mechanical fault based on improved variational mode decomposition[J]. Journal of Vibration Engineering， 2016， 29（6）：1119—1126.

[11]Han T， Jiang D. Rolling bearing fault diagnostic method based on VMD-AR model and random forest classifier[J]. Shock and Vibration， 2016， 2016：1—11.

[12]夏均忠， 趙磊， 白云川，等. 基于MCKD和VMD的滾動軸承微弱故障特征提取[J]. 振動與沖擊， 2017， 36（20）：78—83.

XIAJunzhong， ZHAO Lei， BAI Yunchuan， et al. Feature extraction for rolling element bearing weak fault based on MCKD and VMD[J]. Journal of Vibration and Shock， 2017， 36（20）：78—83.

[13]Li Z， Chen J， Zi Y， et al. Independence-oriented VMD to identify fault feature for wheel set bearing fault diagnosis of high speed locomotive[J]. Mechanical Systems & Signal Processing， 2017， 85：512—529.

[14]Yang Y， Yu D， Cheng J. A fault diagnosis approach for roller bearing based on IMF envelope spectrum and SVM[J]. Measurement， 2007， 40（9-10）：943-950.

[15]胡愛軍， 馬萬里， 唐貴基. 基于集成經(jīng)驗?zāi)B(tài)分解和峭度準則的滾動軸承故障特征提取方法[J]. 中國電機工程學(xué)報， 2012， 32（11）：106-111.

Hu Aijun， Ma Wanli， Tang Guiji. Rolling bearing fault feature extraction method based on ensemble empirical mode decomposition and kurtosis criterion[J]. Proceedings of the CSEE， 2012， 32（11）：106—111.

[16]王平， 廖明夫. 滾動軸承故障診斷的自適應(yīng)共振解調(diào)技術(shù)[J]. 航空動力學(xué)報， 2005， 20（4）：606—612.

WANG Ping， LIAO Ming-fu. Adaptive Demodulated Resonance Technique for the Rolling Bearing Fault Diagnosis[J]. Journal of Aerospace Power， 2005， 20（4）： 606—612.

[17]唐貴基， 王曉龍. 變分模態(tài)分解方法及其在滾動軸承早期故障診斷中的應(yīng)用[J]. 振動工程學(xué)報， 2016， 29（4）：638—648.

TANG Guiji， WANG Xiaolong. Variational mode decomposition method and its application on incipient fault diagnosis of rolling bearing [J]. Journal of Vibration Engineering， 2016， 29 （4）： 638—648.

[18]Liu T， Chen J， Dong G， et al. The fault detection and diagnosis in rolling element bearings using frequency band entropy[J]. Proceedings of Institution of Mechanical Engineers Part C Journal of Mechanical Engineering Science， 2013， 227（1）：87—99.

[19]http：//www.cwru.edu/laboratory/bearing/welcome_over view. htm.

Abstract： According to the difficult problem that the fault features extraction of rolling bearings in early failure duration， an incipient fault diagnosis method for rolling bearings based on the variational mode decomposition （VMD） and the improved adaptive resonance technology （IART） is proposed. According to the problem that the frequency band of the bearing fault signal is difficult to choose， an intrinsic mode function （IMF） selection method based on IART is proposed. Firstly， the mode number determination method based on kurtosis maximum value is proposed here to determine the mode number. Then， the original vibration signal is decomposed by VMD to obtain IMF， and the IMF component with abundant fault information is selected by IART. Finally， the selected IMF component is subjected to band-pass filtering based on IART if necessary， and the fault characteristic frequency is extracted by envelope demodulation analysis. The method can be applied to the bearing simulation data and the actual data， which can realize the accurate diagnosis of bearing fault characteristics and prove the effectiveness of the method.

Key words： fault diagnosis； rolling bearing； variational mode decomposition； kurtosis； improved adaptive resonance technology