甄 冬，田少寧，郭俊超，孟召宗，谷豐收，4

（1.河北工業(yè)大學(xué)機(jī)械工程學(xué)院，天津 300130；2.天津理工大學(xué)天津市先進(jìn)機(jī)電系統(tǒng)設(shè)計(jì)與智能控制重點(diǎn)實(shí)驗(yàn)室，天津 300384；3.天津理工大學(xué)機(jī)電工程國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心，天津 300384；4.Centre for Efficiency and Performance Engineering，University of Huddersfield，Huddersfield HD1 3DH）

引言

滾動(dòng)軸承作為機(jī)械設(shè)備的重要零件之一，在現(xiàn)代工業(yè)中得到了廣泛的應(yīng)用。然而，由于常年在惡劣環(huán)境、變速和交變載荷下運(yùn)行，滾動(dòng)軸承不可避免地會(huì)產(chǎn)生各種失效故障［1-2］。這些早期故障會(huì)導(dǎo)致設(shè)備損壞，并造成重大經(jīng)濟(jì)損失甚至致命的事故。因此，尋求精確的滾動(dòng)軸承早期故障診斷方法和有效的維護(hù)策略是當(dāng)務(wù)之急。在實(shí)際工程中，當(dāng)滾動(dòng)軸承存在早期局部故障時(shí)，將產(chǎn)生一系列微弱的脈沖信號(hào)。然而，這些微弱脈沖信號(hào)通常被強(qiáng)烈的背景噪聲和復(fù)雜的干擾信號(hào)所淹沒(méi)，使得滾動(dòng)軸承的故障診斷變得相當(dāng)復(fù)雜和困難［3-4］。因此，探索一種準(zhǔn)確有效的滾動(dòng)軸承早期微弱故障診斷方法具有重要現(xiàn)實(shí)意義和工程應(yīng)用價(jià)值。

目前，許多典型的故障特征提取方法已被廣泛應(yīng)用于滾動(dòng)軸承故障診斷中，包括隨機(jī)共振（SR）［5］、Wigner-Ville 分布（WVD）［6］、奇異值分解（SVD）［7］和稀疏分解［8］等。盡管這些方法對(duì)于滾動(dòng)軸承故障檢測(cè)具有一定的效果，但是其固有局限性也影響到實(shí)際的應(yīng)用。例如，WVD 具有高的時(shí)頻分辨率和能量集中性等優(yōu)點(diǎn)，但存在嚴(yán)重的交叉項(xiàng)干擾，阻礙了其對(duì)信號(hào)的有效分析和各分量參數(shù)的提取。SR 是一種提取瞬態(tài)特征的有效工具，但是其系統(tǒng)參數(shù)難以確定，因此檢測(cè)能力受到限制。稀疏分解具有良好的信號(hào)分解性能，但過(guò)多依賴于原子庫(kù)的設(shè)計(jì)和分解算法。SVD 是一種有效的降噪方法，但是當(dāng)測(cè)量信號(hào)具有較低的信噪比（SNR）時(shí)，其有效性無(wú)法得到保證。實(shí)際上，復(fù)雜的振動(dòng)信號(hào)包含非平穩(wěn)和非線性分量，非自適應(yīng)信號(hào)處理方法對(duì)于實(shí)際信號(hào)的分析效率并不高。經(jīng)驗(yàn)?zāi)B(tài)分解（EMD）適合分析非線性、非平穩(wěn)信號(hào)［9］，能夠?qū)?fù)雜的信號(hào)分解為具有幅度和頻率調(diào)制特性的本征模態(tài)函數(shù)（IMFs），但是EMD 的端點(diǎn)效應(yīng)和模態(tài)混疊問(wèn)題容易導(dǎo)致IMFs 失去其特定的物理意義。

集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）是由Wu 等［10］提出的一種改進(jìn)型EMD 方法，該方法通過(guò)向測(cè)量信號(hào)添加有限的白色噪聲以消除EMD 端點(diǎn)效應(yīng)和模態(tài)混疊的問(wèn)題。基于上述優(yōu)點(diǎn)，EEMD 在旋轉(zhuǎn)機(jī)械故障診斷領(lǐng)域得到了廣泛的關(guān)注［11-13］。張琛等［11］提出了一種基于EEMD 奇異值熵的故障診斷方法，對(duì)滾動(dòng)軸承的微弱故障進(jìn)行了有效的診斷。胡蔦慶等［12］提出了一種基于EEMD 和深度卷積神經(jīng)網(wǎng)絡(luò)（DCNN）的行星齒輪故障檢測(cè)方法，并對(duì)行星齒輪的故障類型進(jìn)行了分類。Wang 等［13］提出一種基于EEMD 和可調(diào)Q因子小波變換（TQWT）的故障特征提取方法，用于檢測(cè)滾動(dòng)軸承的早期故障。上述研究表明，基于EEMD 分解選擇故障敏感的信號(hào)成分，能夠有效地提高振動(dòng)信號(hào)的信噪比，改善故障特征提取的性能。但是，如何選擇故障敏感的IMFs仍然是一個(gè)棘手的問(wèn)題。因此，國(guó)內(nèi)外學(xué)者對(duì)故障敏感IMFs 的選擇開(kāi)展了大量的研究工作。Chen等［14］開(kāi)發(fā)了一種使用信噪比選擇敏感IMFs 的方法。劉永強(qiáng)等［15］提出應(yīng)用自相關(guān)函數(shù)標(biāo)準(zhǔn)來(lái)確定最佳IMF。Xue 等［16］提出使用峭度準(zhǔn)則來(lái)選擇故障敏感IMF。Osman 等［17］開(kāi)發(fā)了一種基于結(jié)構(gòu)元素尺度來(lái)選擇敏感IMFs 的方法。但是，這些研究?jī)H專注于分析單個(gè)或者幾個(gè)獨(dú)立的IMF 來(lái)提取故障相關(guān)特征，而沒(méi)有考慮到不同IMFs 在揭示故障時(shí)表現(xiàn)出不同程度的有效性。郭俊超等［18］提出一種加權(quán)平均的EEMD（WEEMD）來(lái)揭示每個(gè)IMF 在故障特征提取中的重要性，但是，所獲取的WEEMD 濾波信號(hào)仍然存在隨機(jī)噪聲成分，這將影響故障特征提取的準(zhǔn)確度。鑒于此，本文提出了一種基于累計(jì)均值（MSAM）和小波閾值降噪的EEMD（即IEEMD）分析方法，不僅可以有效地避免遺漏含有重要故障信息的IMFs 分量，而且能夠有效地減少IMFs 中的隨機(jī)噪聲。但是，IEEMD 濾波信號(hào)中仍然存在耦合頻率和線性調(diào)制分量。

調(diào)制信號(hào)雙譜（MSB）被廣泛地應(yīng)用于故障特征提取中，能夠有效地解調(diào)線性調(diào)制成分和檢測(cè)耦合頻率［19-22］。Gu 等［19］利用MSB 對(duì)電動(dòng)機(jī)的轉(zhuǎn)子故障進(jìn)行診斷研究，仿真和實(shí)驗(yàn)分析證明該方法能對(duì)電動(dòng)機(jī)的轉(zhuǎn)子故障進(jìn)行有效的診斷。陳峙等［20］提出使用MSB 檢測(cè)齒輪的故障，并對(duì)齒輪的不同故障類型進(jìn)行了分類。Guo 等［21］將MSB 運(yùn)用到行星齒輪箱的故障檢測(cè)中，通過(guò)仿真和實(shí)驗(yàn)研究驗(yàn)證了MSB在行星齒輪箱狀態(tài)監(jiān)測(cè)和故障診斷方面的可行性。Guo 等［22］提出一種增強(qiáng)MSB 的軸承故障檢測(cè)方法，該方法充分利用AR 模型降低非高斯噪聲和MSB抑制高斯噪聲的特性來(lái)提取故障特征，為滾動(dòng)軸承故障診斷提供了新思路。上述研究表明，MSB 能夠從復(fù)雜的振動(dòng)信號(hào)中有效地解調(diào)信號(hào)中固有的調(diào)制成分和檢測(cè)耦合頻率分量。然而，這些工作假設(shè)振動(dòng)信號(hào)是平穩(wěn)的，在處理非平穩(wěn)信號(hào)時(shí)會(huì)產(chǎn)生某些干擾頻率分量。另外，由于頻譜泄露的影響，MSB的應(yīng)用會(huì)限制其噪聲抑制能力。

綜上所述，本文提出了一種基于IEEMD 和MSB 的故障診斷方法，用于滾動(dòng)軸承微弱故障特征提取。首先，使用EEMD 將滾動(dòng)軸承的振動(dòng)信號(hào)分解為一系列IMFs。鑒于不同的IMFs 在揭示故障特征方面顯示出不同程度的有效性，利用MSAM 準(zhǔn)則將IMFs 自適應(yīng)地分為低頻IMFs 和高頻IMFs，其中高頻IMFs 采用小波閾值降噪進(jìn)行處理。隨后，將降噪后的高頻IMFs 與低頻IMFs 進(jìn)行重構(gòu)以獲取瞬態(tài)脈沖信號(hào)。最后，利用MSB 進(jìn)一步抑制剩余隨機(jī)噪聲和確定性干擾分量，以精確地提取調(diào)制信號(hào)特征。通過(guò)仿真和故障軸承實(shí)驗(yàn)信號(hào)對(duì)IEEMD-MSB 進(jìn)行性能驗(yàn)證。結(jié)果表明，在故障特征提取方面，IEEMD-MSB 的性能優(yōu)于譜峭度（SK）和WEEMD-MSB。

1 改進(jìn)集成經(jīng)驗(yàn)?zāi)B(tài)分解

1.1 IEEMD 算法

集成經(jīng)驗(yàn)?zāi)B(tài)分解（EEMD）是一種自適應(yīng)分析方法，適合處理非平穩(wěn)信號(hào)，能夠?qū)?fù)雜的振動(dòng)信號(hào)自適應(yīng)地分解為一系列含有不同頻率成分的IMFs。

式中y(t)表示測(cè)量信號(hào)；IMFi(i=1，2，…，N)表示第i個(gè)IMF；N表示IMF的數(shù)目；rN(t) 表示殘差分量。

由式（1）可以得出，EEMD本質(zhì)上能夠?qū)?fù)雜的信號(hào)分解為一系列頻率由高到低的幅值調(diào)制和頻率調(diào)制（AM/FM）且零均值的IMFs，實(shí)現(xiàn)了由滾動(dòng)軸承故障引起的調(diào)制信號(hào)成分與干擾信號(hào)成分及隨機(jī)噪聲的分離。為了選擇與故障調(diào)制信號(hào)相關(guān)性強(qiáng)的IMFs，利用基于累計(jì)均值（MSAM）的尺度選擇標(biāo)準(zhǔn)，將高頻IMFs與所有IMFs分開(kāi)。MSAM定義表示為［23］：

綜上所述，IEEMD 方法的具體流程如圖1所示。

圖1 基于IEEMD 的流程圖Fig.1 The flowchart based on IEEMD

1.2 仿真對(duì)比

為了驗(yàn)證IEEMD 在改善降噪效果方面的性能，構(gòu)造滾動(dòng)軸承外圈局部缺陷的仿真信號(hào)進(jìn)行分析［27］：

式中K和Ak分別表示故障脈沖的數(shù)量和第k個(gè)故障脈沖的幅值；ti=t-表示滾動(dòng)軸承的隨機(jī)滑移，其中，T為故障脈沖間隔，τi為零均值均勻分布隨機(jī)變量的第i個(gè)實(shí)現(xiàn)，標(biāo)準(zhǔn)偏差為0.01T～0.02T；ζ和fn分別為阻尼特性因子和故障共振頻率；u(ti)為單位階躍函數(shù)；n(t)表示信噪比為-15 dB 的高斯白噪聲。

仿真信號(hào)的采樣頻率為96 kHz，采樣點(diǎn)數(shù)為960000。圖2 為仿真信號(hào)x(t)的時(shí)域波形及其頻譜，可以看出，周期性脈沖已經(jīng)被噪聲淹沒(méi)。

圖2 仿真信號(hào)Fig.2 Simulation signal

仿真信號(hào)由EEMD 進(jìn)一步處理，該信號(hào)被EEMD 自適應(yīng)地分解為18 個(gè)IMFs。為了有效地獲取故障特征信息，利用MSAM 準(zhǔn)則將分解的IMFs劃分為低頻IMFs 和高頻IMFs，如圖3 所示。根據(jù)式（2）可以得出，區(qū)分高頻IMFs 和低頻IMFs 的尺度為9。其中，高頻IMFs（IMF1～I(xiàn)MF8）使用小波閾值降噪算法進(jìn)行降噪處理，然后將其與低頻IMFs（IMF9～I(xiàn)MF18）重構(gòu)以獲取故障相關(guān)的重構(gòu)信號(hào)，如圖4 所示。

圖3 仿真信號(hào)MSAM 值的計(jì)算結(jié)果Fig.3 Calculation results of MSAM of simulation signal

圖4 IEEMD 濾波結(jié)果Fig.4 The results processed by IEEMD

為體現(xiàn)IEEMD 在降噪性能方面的優(yōu)越性，將IEEMD 的降噪效果與基于相關(guān)系數(shù)準(zhǔn)則的EEMD［15］、基于峭度準(zhǔn)則的EEMD［16］和加權(quán)平均的EEMD（WEEMD）［18］的降噪效果進(jìn)行了對(duì)比。首先應(yīng)用相關(guān)系數(shù)準(zhǔn)則對(duì)EEMD 分解的IMFs 進(jìn)行分析，結(jié)果如表1 所示。由文獻(xiàn)［15］可知，相關(guān)系數(shù)值越大表明含有越多的故障特征信息。因此，選擇IMF1作為最敏感的IMF 進(jìn)行分析，結(jié)果如圖5（a）所示?？梢钥闯?，雖然濾除了一部分噪聲，但在信號(hào)降噪的同時(shí)故障脈沖也在很大程度上被削弱。隨后，應(yīng)用峭度準(zhǔn)則對(duì)EEMD 分解后的IMFs 進(jìn)行分析，結(jié)果如表2 所示。由文獻(xiàn)［16］可知，峭度能夠表征瞬時(shí)故障脈沖的強(qiáng)度，峭度值越大則表明瞬態(tài)脈沖的沖擊越大。因此，選擇IMF2作為最敏感的IMF進(jìn)行分析，結(jié)果如圖5（b）所示?？梢钥闯觯?jīng)濾波后信號(hào)的故障脈沖不明顯，仍然含有大量的背景噪聲。最后，將WEEMD 應(yīng)用于仿真信號(hào)，對(duì)WEEMD 分解后的IMFs 進(jìn)行加權(quán)分析，結(jié)果如圖5（c）所示。從圖5（c）中可以看出，經(jīng)WEEMD 降噪后的信號(hào)仍然包含相當(dāng)成分的噪聲，盡管能夠識(shí)別出仿真信號(hào)的瞬態(tài)脈沖成分，但是其降噪效果不理想。因此，對(duì)比分析結(jié)果表明，IEEMD 對(duì)噪聲的抑制能力明顯優(yōu)于以上三種方法。

表1 IMFs 的相關(guān)系數(shù)值Tab.1 Correlation coefficients of IMFs

表2 IMFs 的峭度值Tab.2 Kurtosis of IMFs

圖5 對(duì)比算法的濾波結(jié)果Fig.5 Process results of comparison algorithms

為了進(jìn)一步評(píng)估IEEMD 的降噪效果，利用信噪比（SNR）和均方根誤差（RMSE）評(píng)測(cè)IEEMD 和上述三種算法的性能，其定義如下：

式中s(l)為原始信號(hào)為降噪信號(hào)；L為采樣長(zhǎng)度。

表3 給出了四種算法降噪性能的評(píng)測(cè)結(jié)果，可以看出，IEEMD 的SNR最大且RMSE最低，降噪效果明顯優(yōu)于其余三種算法，進(jìn)一步驗(yàn)證了IEEMD降噪性能的優(yōu)越性。

表3 四種降噪算法的評(píng)測(cè)結(jié)果Tab.3 Evaluation results of four denoising algorithms

2 調(diào)制信號(hào)雙譜分析

2.1 MSB 算法

調(diào)制信號(hào)雙譜（MSB）分析是一種在常規(guī)雙譜分析的基礎(chǔ)上改進(jìn)的先進(jìn)解調(diào)方法。其原理是通過(guò)抑制振動(dòng)信號(hào)中的隨機(jī)噪聲和干擾分量來(lái)解調(diào)信號(hào)中固有的調(diào)制成分。對(duì)于離散時(shí)間信號(hào)x(t)進(jìn)行相應(yīng)的離散傅里葉變換X(f)，MSB 在頻域中表示為［19］：

式中 上標(biāo)“*”表示復(fù)共軛；BMS(fc，fx)和E<>分別表示信號(hào)x(t)的調(diào)制信號(hào)雙譜和期望算子；fc表示載波頻率；fx表示調(diào)制頻率；(fc+fx)和(fc-fx)分別表示較高和較低的邊帶頻率。

為了得到更精確量化的邊帶幅值，利用幅值歸一化消除fc的影響來(lái)改善MSB 的性能。歸一化MSB 邊帶估計(jì)器（MSB-SE）定義如下［20］：

式中BMS(fc，0)表示fx=0 時(shí)的功率譜平方估計(jì)。MSB-SE 的典型分析結(jié)果如圖6 所示。

圖6 MSB-SE 的分析結(jié)果Fig.6 The analysis results of MSB-SE

2.2 MSB 檢測(cè)器

圖6 所示的MSB-SE 分析結(jié)果表明，用于檢測(cè)故障的最佳頻率處于fc特定值范圍內(nèi)，且最大峰值出現(xiàn)在該頻段內(nèi)。通過(guò)對(duì)fx增量方向上有效的MSB 幅值進(jìn)行平均獲得fc切片，結(jié)果如圖7 所示，其B(fc)定義為［21］：

圖7 MSB 切片的結(jié)果Fig.7 The results of MSB slice

式中 ?f表示fx方向的頻率分辨率；M為某個(gè)特定載波頻率中調(diào)制頻率（故障特征頻率）的倍頻數(shù)。

為了使獲得的結(jié)果更可靠，對(duì)圖7 中標(biāo)有“*”的幾個(gè)MSB 切片進(jìn)行平均來(lái)進(jìn)一步得出MSB 檢測(cè)器，如下式所示：

MSB 檢測(cè)器的結(jié)果如圖8 所示。

圖8 MSB 檢測(cè)器的結(jié)果Fig.8 The results of MSB detector

3 故障診斷方法的流程

本文所提出的基于IEEMD 和MSB 的軸承微弱故障特征提取方法，其具體的診斷步驟如下，流程圖如圖9 所示。

圖9 IEEMD-MSB 的故障特征提取流程圖Fig.9 The fault feayure extraction flowchart of IEEMD-MSB

（1）應(yīng)用EEMD 將原始信號(hào)分解為本征模態(tài)函數(shù)（IMFs）集。

（2）計(jì)算每個(gè)IMF 的累計(jì)均值（MSAM）。

（3）對(duì)步驟（2）中得出的MSAM 值進(jìn)行判斷。如果MSAM 值在第m層顯示偏離零，則第m層之前的IMFs 為高頻IMFs，其余IMFs 為低頻IMFs。

（4）采用小波閾值降噪對(duì)高頻IMFs 進(jìn)行降噪處理，然后將其與低頻IMFs 結(jié)合以獲得重構(gòu)信號(hào)。

（5）利用MSB 對(duì)重構(gòu)信號(hào)進(jìn)行解調(diào)并提取故障特征。

4 實(shí)驗(yàn)研究

為了驗(yàn)證所提方法的有效性，開(kāi)展實(shí)驗(yàn)研究。滾動(dòng)軸承的試驗(yàn)臺(tái)如圖10 所示。試驗(yàn)臺(tái)由交流電動(dòng)機(jī)、兩個(gè)振動(dòng)傳感器、兩個(gè)支撐滾動(dòng)軸承以及三個(gè)撓性聯(lián)軸器所組成。實(shí)驗(yàn)裝置布置如下：一個(gè)振動(dòng)傳感器安裝在電機(jī)驅(qū)動(dòng)端滾動(dòng)軸承座上；另一個(gè)振動(dòng)傳感器放置在支撐滾動(dòng)軸承座上。實(shí)驗(yàn)過(guò)程中，交流電動(dòng)機(jī)轉(zhuǎn)速為1470 r/min，實(shí)驗(yàn)數(shù)據(jù)的采樣頻率和采樣長(zhǎng)度分別為96 kHz 和1920000。實(shí)驗(yàn)軸承故障模式為人工模擬的局部點(diǎn)蝕故障，包括電動(dòng)機(jī)滾動(dòng)軸承外圈故障和支撐滾動(dòng)軸承內(nèi)圈故障，如圖11 所示。表4 為故障滾動(dòng)軸承的主要參數(shù)；表5 為滾動(dòng)軸承的故障特征頻率。

表4 滾動(dòng)軸承的主要參數(shù)Tab.4 Main parameters of rolling bearings

表5 滾動(dòng)軸承的故障特征頻率Tab.5 Fault characteristic frequency of rolling bearings

圖10 滾動(dòng)軸承試驗(yàn)臺(tái)Fig.10 Rolling bearing experiment platform

圖11 故障模式Fig.11 The fault modes

4.1 電機(jī)軸承外圈故障分析

圖12 為實(shí)驗(yàn)采集的電機(jī)軸承外圈故障振動(dòng)信號(hào)的波形及其頻譜。顯然，外圈故障特征頻率的高次諧波被準(zhǔn)確地識(shí)別，但故障特征頻率fo被強(qiáng)背景噪聲和干擾分量所淹沒(méi)，如圖12（b）所示。

圖12 電機(jī)軸承外圈故障信號(hào)Fig.12 Outer race fault signal of motor bearing

為了準(zhǔn)確地提取故障特征頻率及其諧波，首先應(yīng)用EEMD 將測(cè)得的電機(jī)軸承外圈故障信號(hào)自適應(yīng)地分解為19 個(gè)IMFs。圖13 表示每個(gè)IMF 的MSAM 值。可以看出，區(qū)分高頻IMFs 和低頻IMFs的合適尺度為11。隨后，應(yīng)用小波閾值降噪對(duì)高頻IMFs（IMF1～I(xiàn)MF10）進(jìn)行降噪處理，并將降噪后的高頻IMFs 與低頻IMFs（IMF11～I(xiàn)MF19）組合以獲取重構(gòu)信號(hào)。同時(shí)，SNR 和RMSE 指標(biāo)被用來(lái)評(píng)估IEEMD 和第1.2 節(jié)仿真分析中所提及三種算法的降噪性能。表6 為四種算法降噪性能的評(píng)估結(jié)果，可以看出，IEEMD 的SNR最大且RMSE最低，降噪性能明顯高于其余三種算法，進(jìn)一步驗(yàn)證了IEEMD 降噪性能的優(yōu)越性。

表6 四種降噪算法的評(píng)測(cè)結(jié)果Tab.6 Evaluation results of four denoising algorithms

圖13 外圈故障信號(hào)MSAM 值的計(jì)算結(jié)果Fig.13 Calculation results of MSAM for outer race fault signal

最后，利用MSB 對(duì)重構(gòu)信號(hào)進(jìn)行處理，提取與故障相關(guān)的特征頻率，其分析結(jié)果如圖14 所示?？梢钥闯?，IEEMD-MSB 方法可以清楚地識(shí)別與故障相關(guān)的特征頻率及其諧波。

圖14 IEEMD-MSB 的分析結(jié)果Fig.14 Analysis results of IEEMD-MSB

為了進(jìn)一步證明所提方法的優(yōu)越性，本文采用譜峭度（SK）［28］和WEEMD-MSB［18］分析方法對(duì)測(cè)得的電機(jī)軸承外圈故障的振動(dòng)信號(hào)進(jìn)行分析，結(jié)果如圖15，16 所示。從圖15（b）中雖然可以識(shí)別故障特征頻率fo及其諧波，但是SK 算法的分析結(jié)果卻被大量的背景噪聲和干擾頻率所干擾。圖16（b）中也存在一些干擾頻率成分，尤其是在故障特征頻率的高次諧波附近。綜上分析得出，相比SK 和WEEMDMSB，IEEMD-MSB 在電機(jī)軸承外圈故障診斷中能夠取得更為準(zhǔn)確的結(jié)果。

圖15 SK 的分析結(jié)果Fig.15 Analysis results of SK

圖16 WEEMD-MSB 的分析結(jié)果Fig.16 Analysis results of WEEMD-MSB

4.2 支撐軸承內(nèi)圈故障分析

圖17 為實(shí)驗(yàn)采集的支撐軸承內(nèi)圈故障振動(dòng)信號(hào)的波形及其頻譜。顯然，故障特征頻率fi及其諧波被背景噪聲和干擾諧波所淹沒(méi)，如圖17（b）所示，無(wú)法準(zhǔn)確、有效地識(shí)別故障特征信息。

圖17 支撐軸承內(nèi)圈故障信號(hào)Fig.17 Inner race fault signal of supporting bearing

為了準(zhǔn)確地提取故障特征頻率及其諧波，首先利用EEMD 將測(cè)得的支撐軸承內(nèi)圈故障信號(hào)自適應(yīng)地分解為18 個(gè)IMFs。圖18 表示每個(gè)IMF 的MSAM值?？梢园l(fā)現(xiàn)，區(qū)分高頻IMFs和低頻IMFs的合適尺度為7。隨后，應(yīng)用小波閾值降噪對(duì)高頻IMFs（IMF1～I(xiàn)MF6）進(jìn)行降噪處理，并將降噪后的高頻IMFs 與低頻IMFs（IMF7～I(xiàn)MF18）進(jìn)行組合以獲取重構(gòu)信號(hào)。同樣，計(jì)算IEEMD 和第1.2節(jié)仿真分析中所提及三種算法的SNR 和RMSE 指標(biāo)，結(jié)果如表7所示。從表7中可以看出，IEEMD 的降噪效果優(yōu)于其余三種算法，進(jìn)一步證明了IEEMD 的降噪性能。

表7 四種降噪算法的評(píng)測(cè)結(jié)果Tab.7 Evaluation results of four denoising algorithms

圖18 內(nèi)圈故障信號(hào)MSAM 值的計(jì)算結(jié)果Fig.18 Calculation results of MSAM for inner race fault signal

最后，使用MSB 對(duì)重構(gòu)信號(hào)進(jìn)行解調(diào)分析，其結(jié)果如圖19 所示。從圖19（b）中可以看出，IEEMDMSB 可以清晰地識(shí)別與故障相關(guān)的特征頻率及其諧波。

圖19 IEEMD-MSB 的分析結(jié)果Fig.19 Analysis results of IEEMD-MSB

為了更好地對(duì)比，采用SK 和WEEMD-MSB 對(duì)所采集的振動(dòng)信號(hào)進(jìn)行分析，結(jié)果如圖20 和21 所示。從圖20（b）中雖然可以提取故障特征頻率fi及其諧波，但是混有較多的背景噪聲和干擾諧波。對(duì)于圖21（b）WEEMD-MSB 的分析結(jié)果，盡管能夠識(shí)別出故障特征頻率，但故障特征頻率幅值和干擾成分幅值的對(duì)比度較差。因此，IEEMD-MSB 比SK和WEEMD-MSB 更能有效、準(zhǔn)確地檢測(cè)支撐軸承內(nèi)圈故障。

圖20 SK 的分析結(jié)果Fig.20 Analysis results of SK

圖21 WEEMD-MSB 的分析結(jié)果Fig.21 Analysis results of WEEMD-MSB

5 結(jié)論

針對(duì)滾動(dòng)軸承振動(dòng)信號(hào)的強(qiáng)非線性和非平穩(wěn)特性，提出了一種基于IEEMD 和MSB 的軸承微弱故障特征提取方法。通過(guò)對(duì)電機(jī)軸承外圈和支撐軸承內(nèi)圈故障診斷的研究，得出以下結(jié)論：

（1）IEEMD 能夠有效地解決MSB 處理非平穩(wěn)信號(hào)時(shí)的不足，并進(jìn)一步突出振動(dòng)信號(hào)的調(diào)制分量。

（2）通過(guò)MSB 對(duì)IEEMD 濾波信號(hào)進(jìn)行殘余噪聲抑制和調(diào)制成分解調(diào)，可以有效地改善故障特征提取精度。

（3）通過(guò)對(duì)不同類型實(shí)驗(yàn)軸承故障的診斷，驗(yàn)證了所提方法的有效性，且該方法的故障特征提取能力優(yōu)于WEEMD-MSB 和SK，能夠更有效地提高滾動(dòng)軸承故障診斷的精度。

（4）從上述分析可以得出，基于IEEMD 和MSB的滾動(dòng)軸承微弱故障特征提取方法在實(shí)際應(yīng)用中的潛力值得進(jìn)一步挖掘和研究。