基于VMD和MED的滾動軸承故障診斷方法

郭家宇，熊炘，劉浩

（上海大學 機電工程與自動化學院，上海 200072）

滾動軸承的健康狀況影響著整個機械系統的工作狀態[1]，通過狀態監測和故障特征提取技術了解滾動軸承的運行狀態，及早發現故障問題，可以有效提高機械設備的運行水平并降低維修成本。然而，從機械設備現場采集的滾動軸承振動信號由于受嚴重的背景噪聲和其他隨機信號干擾，所獲信號往往具有非線性、非平穩、非高斯的特征，信噪比較低；特別是齒輪箱中的軸承故障特征信息往往淹沒于齒輪嚙合頻率及背景噪聲，嚴重影響了故障特征信號的識別，用常規的時域分析方法（時域波形分析、自相關分析、峭度指標和波形指標等）與頻域分析方法（幅值譜、功率譜、倒譜和包絡分析等）難以有效地進行共振頻帶的選取，無法準確提取故障特征。

時頻分析方法[2-6]能夠同時再現信號的時、頻域特征，對處理非線性、非平穩信號非常有效，但均具有各自的局限性：小波變換可以有效抑制白噪聲，但小波基函數和閾值的選擇不具有自適應性；經驗模態分解（Empirical Mode Decomposition，EMD）方法分析非平穩、非線性信號的能力突出，但其數學基礎不嚴格，模態混疊和端點效應問題表現明顯[7]；集合經驗模態分解（Ensemble Empirical Mode Decomposition，EEMD）對 EMD加以改進，在一定程度上抑制了模態混疊的產生，但算法的計算量過大[8]；局部均值分解（Local Mean De-composition，LMD）方法亦存在著模態混疊、迭代計算量大以及端點效應等問題。

變分模態分解（Variational Mode Deconvolution，VMD）可以自適應地分析信號，具有堅實的理論基礎和更好的噪聲魯棒性，且相比于EMD，EEMD和LMD等算法，其分解層數少、精度高，不存在模態混疊現象[9-11]。最小熵解卷積（Minimum Entropy Deconvolution，MED）應用于混合多源干擾信號中沖擊信號的提取，能有效減弱信號受采集路徑而衰減的影響，可以進一步突出振動信號的沖擊特征，在滾動軸承故障特征提取中取得了較好的分析效果[12-14]。因此，將VMD與MED相結合，用于提取滾動軸承故障特征，同時與EMD進行比較，通過仿真信號和實測信號驗證此方法的有效性及優越性。

1 變分模態分解

變分模態分解是建立在Wiener濾波、Hilbert變換、解析信號、頻率混合和外差解調等概念基礎上的非遞歸、自適應的信號處理方法，具有嚴格的數學基礎。其本質是通過迭代搜尋變分模型最優解來確定每個分量的頻率中心及帶寬，從而能自適應地實現信號頻域的有效分離[15]。VMD算法可解釋為變分問題的構造和求解，對輸入信號的具體處理步驟可分為變分模型的建立和變分問題的求解2個步驟。

1.1 變分模型的建立

假設輸入信號yp（t）是由有限個中心頻率不同、帶寬有限的離散模態分量 uk（t）（k＝1，2，…，K）組成。通過Hilbert變換將各模態分量uk（t）轉換成解析信號，獲得各分量單邊頻譜。接著將各模態分量的解析信號與e-jωkt相乘，各分量頻譜被調制到相應基頻帶。對分量解析信號計算梯度的平方L2范數，估計出中心脈沖頻率ωk的帶寬。引入約束條件，構建受約束的變分模型為

式中：f為原始信號；uk為模態函數；ωk為各個模態的中心頻率；K為模態分量數目。

1.2 變分問題的求解

在上述變分模型中引入α和構造求解無約束性變分問題的增廣拉格朗日函數 L（｛uk｝，｛ωk｝，λ），即

式中：α為二次懲罰因子，用于保證處于Gauss噪聲下信號的重構精度；λ（t）為Lagrange乘法算子，使得約束條件保持嚴格性。

利用乘法算子交替方向法（Alternate Direction Method of Multipliers，ADMM）對無約束性變分問題求解，交替迭代更新un+1k，ωn+1k，λn+1，尋找增廣Lagrange函數的鞍點，即為（2）式約束變分模型的最優解，從而將原始信號分解為K個窄帶本征模態分量[16]。最終求得的更新表達式為

解得中心頻率的更新方法為

1.3 VMD算法流程

2 最小熵解卷積

滾動軸承振動信號經VMD處理后得到離散模態分量uk（t），為提取信號中較大的尖脈沖成分，恢復原始的輸入沖擊信號，篩選出和原信號相關系數大的分量進行重構，需要對得到的重構信號進行MED處理。

3 實際驗證

3.1 故障診斷流程

綜上所述，將VMD和MED相結合，首先，將原始故障采集信號yp（t）進行VMD處理，選取與原信號相關系數大的分量重構；然后，將重構信號y（t）進行MED處理，去除采集路徑因素的影響，提高信號信噪比，凸顯故障沖擊成分；最后，對VMD-MED處理后的濾波信號x（t）采用包絡分析，找出信號特征頻率，對比軸承故障頻率可獲得診斷結果。故障診斷流程如圖1所示。

圖1 VMD＋MED的故障診斷流程圖Fig.1 Flow chart of fault diagnosis of VMD＋MED

3.2 試驗條件

滾動軸承故障信號采集自如圖2所示的二級齒輪減速器平臺，其選用標準直齒輪，齒輪齒數分別為Z1＝32，Z2＝64，Z3＝96，齒輪模數為 1.5。交流電動機輸出軸通過聯軸器與減速器輸入軸連接。振動信號來源于固定在位于減速器輸入軸處左側故障軸承軸承座上端的振動加速度傳感器，采樣頻率為25.6 kHz。試驗采用6004-2RSH型深溝球軸承（具體參數見表1），利用線切割方式在軸承內圈溝道上加工寬0.3 mm、深0.2 mm的裂紋故障，從而模擬軸承內圈早期故障缺陷。

圖2 二級減速器試驗平臺Fig.2 Secondary reducer test platform

表1 試驗軸承基本參數Tab.1 Basic parameters for test bearing

3.3 軸承故障診斷

試驗中，電動機轉速為2 400 r/min，即轉頻fr＝40 Hz，計算可得軸承內圈故障特征頻率為fi＝215.6 Hz。選取采樣點數為16 384，實際測得內圈微弱故障原始信號及其頻譜如圖3所示，從時域圖中雖然可見部分脈沖沖擊，但是其被強烈背景噪聲所干擾，無法獲取有效故障特征信息；頻譜中的高頻成分較多，故障特征頻率同樣難以觀察。

圖3 內圈故障原始信號及其頻譜Fig.3 Original fault signal of inner ring and its spectrum

考慮到特征頻率范圍和計算效率，選擇分解層數為5層，α采用VMD默認值2 000；τ選為0.3以保證實際信號分解的保真度。內圈故障信號VMD處理結果如圖4所示。

圖4 VMD分解所得分量Fig.4 Components acquired by VMD

計算各分量與原信號之間的相關系數，計算結果見表2?；谙嚓P系數較大的原則選取V1和V2分量對信號進行重構，降噪后信號的頻譜圖如圖5所示，從圖中可以看到信號中的主要高頻干擾被濾除。

表2 各模態分量和原信號的相關系數Tab.2 Correlation coefficient between each modal component and original signal

圖5 內圈故障重構信號的頻譜Fig.5 Spectrum of reconstructed fault singal inner ring

運用MED對經過VMD分解并重構的信號進行處理，所得信號的時域波形如圖6所示，從圖中可以看到故障沖擊信息明顯加強。

圖6 MED處理后信號時域圖Fig.6 Time domain signal after MED processing

為驗證VMD＋MED算法的優越性，對內圈故障原始信號采用不同算法進行處理，結果如圖7所示。從圖中可以看出：圖7a和圖7b均能提取出內圈故障特征頻率及其1～4倍頻，而圖7c中則清晰可見內圈故障特征頻率及其1～9倍頻，說明經過VMD和MED方法處理后的信號包絡譜能夠凸顯更多的故障特征頻率倍頻信息。

圖7 內圈故障原始信號的處理結果Fig.7 Processing results of original fault signal of inner ring

進一步對軸承外圈及鋼球故障信號進行分析，VMD＋MED算法均可以較好的提取出強背景噪聲下滾動軸承的故障信息，且由于VMD的分解層數更少，于滾動軸承故障在線診斷中更具優勢。

4 結束語

針對齒輪箱中工作的滾動軸承受強背景噪聲和齒輪嚙合頻率的影響，難以用常規分析方法提取故障特征頻率，將VMD和MED相結合，用于提取軸承故障特征信息，經過對齒輪箱故障軸承的振動信號分析表明：VMD降噪時分解層數少且分解精度較高，具有良好的分解效率，信號中的背景噪聲能有效地被濾除掉；基于VMD和MED的方法可以有效提取強背景噪聲下的滾動軸承故障特征，且相對于基于EMD和MED的特征提取方法提取效果更優。