李華 伍星 劉韜

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

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

中圖分類號(hào)：TH165+.3； TN911.7文獻(xiàn)標(biāo)志碼： A文章編號(hào)： 1004-4523（2018）04-0718-09

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

引言

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

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

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

1基礎(chǔ)理論

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

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

1.2峭度

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

其數(shù)學(xué)描述為K=Ex-μ4σ4（6）式中x為所分析的振動(dòng)信號(hào)，μ為信號(hào)x的均值，σ為信號(hào)x的標(biāo)準(zhǔn)差。

峭度對(duì)沖擊信號(hào)敏感，適合軸承早期故障診斷。滾動(dòng)軸承無故障時(shí)，峭度指標(biāo)值約等于3。當(dāng)滾動(dòng)軸承發(fā)生局部故障時(shí)，峭度值會(huì)增大。

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

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

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

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

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

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

因?yàn)榍投戎笜?biāo)是故障的一種表征，當(dāng)存在某個(gè)模態(tài)的峭度值具有最大值時(shí)，可認(rèn)為此時(shí)具有良好的分解效果（在給定的搜索范圍內(nèi)）?；诖?，本文提出了基于峭度最大值優(yōu)化模態(tài)數(shù)K的方法。具體步驟及流程在2.3節(jié)中給出。

2.2基于IART選取敏感IMF分量

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

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

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

2.3算法步驟及流程

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

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

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

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

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

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

3信號(hào)仿真分析

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

3.2有效IMF選取及帶通濾波

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

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

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

4實(shí)驗(yàn)驗(yàn)證及分析

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

4.1VMD分析

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

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

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

4.2有效IMF選取及帶通濾波

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

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

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

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

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

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

4.3與EEMD對(duì)比

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

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

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

4.4.1VMD分析

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

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

4.4.2有效IMF選取

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

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

5結(jié)論

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

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

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

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

參考文獻(xiàn)：

[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 小波的滾動(dòng)軸承微弱故障診斷方法[J]. 振動(dòng)工程學(xué)報(bào)， 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均值聚類的滾動(dòng)軸承故障診斷[J]. 中國電機(jī)工程學(xué)報(bào)， 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]. 振動(dòng)與沖擊， 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]劉尚坤， 唐貴基， 王曉龍. 基于改進(jìn)變分模態(tài)分解的旋轉(zhuǎn)機(jī)械故障時(shí)頻分析方法[J]. 振動(dòng)工程學(xué)報(bào)， 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的滾動(dòng)軸承微弱故障特征提取[J]. 振動(dòng)與沖擊， 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)驗(yàn)?zāi)B(tài)分解和峭度準(zhǔn)則的滾動(dòng)軸承故障特征提取方法[J]. 中國電機(jī)工程學(xué)報(bào)， 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]王平， 廖明夫. 滾動(dòng)軸承故障診斷的自適應(yīng)共振解調(diào)技術(shù)[J]. 航空動(dòng)力學(xué)報(bào)， 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)分解方法及其在滾動(dòng)軸承早期故障診斷中的應(yīng)用[J]. 振動(dòng)工程學(xué)報(bào)， 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