楊 斌, 張家瑋, 樊改榮， 王建國

(1.內(nèi)蒙古科技大學(xué) 機(jī)械工程學(xué)院，內(nèi)蒙古 包頭 014010； 2.內(nèi)蒙古一機(jī)集團(tuán)第四分公司，內(nèi)蒙古 包頭 014032)

滾動軸承作為旋轉(zhuǎn)機(jī)械中的重要零部件，其一旦出現(xiàn)單一微弱的故障或者局部損傷，并不會馬上停機(jī)檢查，設(shè)備如果一直在惡劣的環(huán)境中工作，單一微弱故障將會變成強(qiáng)故障，進(jìn)而導(dǎo)致復(fù)合故障的發(fā)生，復(fù)合故障將導(dǎo)致旋轉(zhuǎn)機(jī)械運(yùn)轉(zhuǎn)狀況的惡化，引起相關(guān)部件相繼發(fā)生故障，因此研究滾動軸承復(fù)合故障診斷方法具有重要的意義[1]。當(dāng)故障和干擾成分同時存在時，故障沖擊成分因?yàn)槟芰啃〕１谎蜎]于強(qiáng)烈的背景噪聲中，故對于復(fù)合故障而言，復(fù)合故障中的微弱成分特征提取更加困難，會容易造成漏診和誤診，換言之，復(fù)合故障特征的分離、提取是研究的重點(diǎn)[2]。

最小熵反褶積(Minimum Entropy Deconvolution，MED)方法是一種系統(tǒng)辨識方法，廣泛應(yīng)用于地震信號中的反射參數(shù)信息的提取[3]。基本原理是解反褶積結(jié)果突出少數(shù)大的尖脈沖，對旋轉(zhuǎn)機(jī)械沖擊性故障的降噪處理非常適用，但是MED方法在降噪時未考慮信號的周期性，降噪效果受限，并限制了其故障特征提取能力。為了使軸承周期性沖擊信號更加明顯，方便進(jìn)行后續(xù)處理，Mcdonald在MED的基礎(chǔ)上研究出了最大相關(guān)峭度解卷(Maxim Correlated Kurtosis Deconvolution，MCKD)方法。該方法選取一個有限沖擊響應(yīng)濾波器f使周期已知信號通過濾波后相關(guān)峭度值達(dá)到最大，并且通過解卷積運(yùn)算使信號被噪聲掩蓋的周期沖擊成分變的突出。

眾多學(xué)者對MCKD方法展開了一系列研究，胡勇[4]將MCKD方法應(yīng)用到滾動軸承故障診斷中，利用MCKD清除傳遞路徑的影響，使故障特征成分得到一次增強(qiáng)。呂宏強(qiáng)等[5]將自回歸(AR)模型和MCKD方法相結(jié)合有效提取了齒輪點(diǎn)蝕故障特征，包絡(luò)譜體現(xiàn)了齒輪點(diǎn)蝕過程。任學(xué)平等[6]將MCKD方法和計(jì)算階次跟蹤相結(jié)合提取滾動軸承周期成分，減少現(xiàn)場周圍噪聲的干擾。

但是MCKD方法中的重要參數(shù)，濾波器長度L和沖擊周期T選擇方法尚未確定，上述方法中并未明確。濾波器長度與沖擊周期決定著MCKD方法降噪和故障特征提取性能的好壞。王建東等[7]運(yùn)用粒子群算法對參數(shù)L進(jìn)行自適應(yīng)尋優(yōu)，呂中亮等[8]引入網(wǎng)格搜索法搜尋最優(yōu)濾波器階數(shù)以及周期。但是利用上述方法存在選擇結(jié)果不穩(wěn)定，只對部分振動信號適用，并且需要設(shè)定合理的采樣頻率，以達(dá)到搜尋要求。因此，部分學(xué)者仍然按照Mcdonald等[9]提出的推薦選擇范圍通過實(shí)驗(yàn)選取，朱軍[10]雖然提出利用信噪比選擇濾波器長度L來判斷MCKD的降噪性能，但是對于沖擊周期T并沒有給出合適的選擇方法。

基于以上分析，本文提出了一種最優(yōu)參數(shù)MCKD選擇方法。通過引入包絡(luò)譜稀疏度判定不同沖擊周期值下MCKD的故障特征提取效果；以排列熵值作為選擇最佳的濾波器長度的度量指標(biāo)，使MCKD達(dá)到最佳的降噪和故障特征提取效果。

此外，本文將最優(yōu)參數(shù)MCKD方法與ELMD方法結(jié)合對軸承復(fù)合故障診斷進(jìn)行應(yīng)用研究，對軸承復(fù)合故障信號進(jìn)行有效分離和提取。首先利用最優(yōu)參數(shù)MCKD提取沖擊明顯的故障；其次，運(yùn)用ELMD進(jìn)行信號分解，對各個PF分量進(jìn)行希爾伯特包絡(luò)解調(diào)，判斷出現(xiàn)的頻率特征；再次利用最優(yōu)參數(shù)MCKD方法對其他故障信號進(jìn)行分離和故障提取。

1 MCKD方法原理

相關(guān)峭度是在峭度基礎(chǔ)上提出來的概念，它充分考慮沖擊成分的連續(xù)性和周期性，使得信號被噪聲掩蓋的周期沖擊成分變的突出，具體步驟如下

相關(guān)峭度的定義為

(1)

式中：M為移位數(shù)；T為沖擊信號的周期；yn為降噪后的信號。移位數(shù)M的增加會使脈沖個數(shù)越來越多，使算法的故障檢測能力得到提升。移位數(shù)M根據(jù)經(jīng)驗(yàn)選取為1～7，當(dāng)M>7時會出現(xiàn)迭代方法超出浮點(diǎn)指數(shù)范圍導(dǎo)致計(jì)算精度的降低,本文取M為5。

MCKD的目標(biāo)函數(shù)為

(2)

為了得到使CKM(T)取得最大值的最優(yōu)濾波器求解方程，令：

(3)

k=1,2,…,L

求得結(jié)果以矩陣的形式表述為

(4)

其中xr為移位數(shù)M和沖擊周期T確定的矩陣

r=0,T, 2T,…,mT；

通過迭代方式MCKD的算法流程如下

步驟1 確定濾波器的長度L，移位數(shù)M和沖擊信號的周期T;

步驟3 求得濾波輸出信號y(n)；

步驟4 根據(jù)y(n)計(jì)算αm與β；

步驟5 更新濾波器系數(shù)f；

如果濾波前后信號ΔCKM(T)<ε,則停止迭代，跳回步驟3。

2 最優(yōu)參數(shù)MCKD方法

2.1 確定沖擊周期T

在MCKD方法中，當(dāng)輸入不同的沖擊周期時，沖擊信號將被MCKD算法濾去其他的周期成分，提取期望的故障特征。理論上，故障發(fā)生部位的沖擊間隔與算法中設(shè)定的沖擊周期T相等，MCKD方法將達(dá)到最佳的解卷積效果。但是實(shí)際工程應(yīng)用中，故障部位的脈沖間隔是未知的[11]，因此解卷積周期參數(shù)的選擇關(guān)系著信號中隱含的周期性脈沖序列能否被成功突顯出來，同時也是MCKD算法進(jìn)行軸承復(fù)合故障診斷的關(guān)鍵。當(dāng)信號x(n)的幅值分布較為均勻、差別不大時，信號的稀疏性會呈現(xiàn)比較微弱的狀態(tài)，導(dǎo)致稀疏度值較小，如果信號中出現(xiàn)較明顯的周期性連續(xù)脈沖，信號則呈現(xiàn)出較強(qiáng)的稀疏特性，稀疏度指標(biāo)也隨之增大[12]。稀疏度定義如下

(5)

利用時域信號的稀疏度作為評價標(biāo)準(zhǔn)容易受單個或少量大幅值脈沖影響較大，將信號從時域轉(zhuǎn)換成頻域，可以避免這一問題。因此引入包絡(luò)譜稀疏度指標(biāo)選擇MCKD算法中的沖擊周期T，當(dāng)軸承發(fā)生故障時，采集到的振動信號沖擊是周期性的，通過設(shè)定周期范圍和包絡(luò)譜稀疏度的對應(yīng)關(guān)系，利用觀察到的峰值篩選合適信號的周期值既可以確定MCKD中的參數(shù)亦可以初步判斷出故障部位。最佳周期選擇過程如下：

(1) 為了減少無用運(yùn)算，提高參數(shù)的搜索效率，同時根據(jù)軸承各部位故障頻率計(jì)算得到的沖擊間隔，設(shè)定故障沖擊范圍為1～100，即可滿足實(shí)驗(yàn)要求。

(2) 設(shè)定比較小的初始濾波器長度L，同時為了避免降噪效果不佳對周期沖擊的干擾，本文選取為150。

(3) 利用搜索范圍內(nèi)的各數(shù)值點(diǎn)作為沖擊周期參數(shù)T，利用MCKD算法對故障信號進(jìn)行處理，并計(jì)算包絡(luò)譜稀疏度值，按順序每次增加一個單位沖擊周期T，直到?jīng)_擊周期T到達(dá)100為止，最佳沖擊周期就是包絡(luò)譜稀疏度對應(yīng)達(dá)到最大的那個值。

2.2 確定濾波器長度L

濾波器長度過大會使計(jì)算量增大，導(dǎo)致故障診斷時間過長，濾波器長度過小，MCKD算法不能起到合理的降噪效果。排列熵值可以作為度量時間序列的復(fù)雜性和檢測動力學(xué)突變的指標(biāo)，信號排列熵值越大，信號中的高頻故障沖擊成分越多?；诖吮疚囊肱帕徐睾饬坎煌瑸V波器長度的降噪效果。在選擇合適的濾波器長度之前，固定好已經(jīng)選擇的最佳周期，逐漸增加濾波器長度。根據(jù)Mcdonald的研究結(jié)果初始值設(shè)定為100，考慮計(jì)算成本和濾波器長度過大會導(dǎo)致降噪后的信號失去原信號的特性等因素[13]，將其限定在100～500，根據(jù)信號排列熵值越大高頻故障沖擊成分越多的原理，選擇排列熵峰值對應(yīng)的濾波器長度。

需要注意的是先選擇沖擊周期再選擇濾波器長度是較為合理的順序，如果不能提前確定故障沖擊，MCKD算法即使對該信號進(jìn)行降噪處理，可能會導(dǎo)致沖擊周期設(shè)定錯誤無法提取正確故障特征。

2.3 仿真信號分析

為了驗(yàn)證最優(yōu)參數(shù)MCKD方法可以更有效降噪，突出故障沖擊成分，分別計(jì)算原信號經(jīng)過最優(yōu)參數(shù)MCKD方法和MCKD方法處理后的時域波形圖的峭度值來衡量信號中包含故障信號的強(qiáng)弱程度。仿真信號x(t)由以下信號混合而成，f1=4 000 Hz，fo=30 Hz。

x1(t)=sin(800πt)

x2(t)=2exp(-350t1)sin(2πf1t),t1=mod(t,1/f0)

x3(t)=0.4randn(1,n),n=length(t)

x(t)=x1(t)+x2(t)+x3(t)

仿真波形圖如圖1所示。

圖1 仿真信號的時域波形(峭度值=2.72)

利用本文的方法選取MCKD的參數(shù)，首先固定好濾波器長度L=150，確定最佳沖擊周期T。在沖擊范圍為1～100之間依次對仿真信號進(jìn)行MCKD處理，求取包絡(luò)譜稀疏度，得到的不同沖擊周期下的包絡(luò)譜稀疏度值，如圖2所示，可以看出當(dāng)沖擊周期T為33時，包絡(luò)譜稀疏度值達(dá)到最大。因此沖擊周期選擇為33。圖3為排列熵與濾波器長度的關(guān)系圖，它的選擇范圍較廣，不能像沖擊周期的選擇一樣，每次增加1，會導(dǎo)致計(jì)算效率的損失，故每次增加10。觀察圖3只有當(dāng)濾波器長度為310時，排列熵值達(dá)到頂峰，因此對于仿真信號濾波器長度選擇為310。

圖2 包絡(luò)譜稀疏度與沖擊周期的關(guān)系

圖4為最優(yōu)MCKD參數(shù)下的仿真信號的時域圖與包絡(luò)功率譜圖，峭度值為7.17，說明降噪后的信號沖擊成分比原信號中增多。從圖4可知，時域圖被淹沒的沖擊部分變的突出，干擾成分明顯減少，并且包絡(luò)功率譜提取的頻率成分為仿真信號的基頻與倍頻成分。因此最優(yōu)參數(shù)MCKD方法在降噪和故障特征提取方面具有較好的效果。另從理論計(jì)算的沖擊間隔為33.33 s，接近包絡(luò)譜稀疏度選取的33的沖擊周期，從另一方面也證明了本文方法選擇參數(shù)的有效性。

圖3 排列熵值與濾波器長度的關(guān)系

(a) 時域圖

(b) 包絡(luò)功率譜圖

Fig.4 The time domain graph and envelope power spectrum of the optimal parameter MCKD (Kurtosis value =7.17)

將本文方法與文獻(xiàn)[8]提出的以小波shannon熵為目標(biāo)函數(shù)的最大相關(guān)峭度解卷積參數(shù)尋優(yōu)方法進(jìn)行對比分析，進(jìn)一步驗(yàn)證本文方法的優(yōu)越性，利用小波shannon熵同時確定的最后濾波器長度為313、沖擊周期為28，其包絡(luò)功率譜如圖5所示。從圖5可知，包絡(luò)功率譜圖中干擾成分較多，確定的基頻成分為36和計(jì)算出的故障特征頻率有一定偏差。這是因?yàn)楦鶕?jù)熵值的原理可以知道其對確定濾波器長度效果較好，而將其作為沖擊周期的評定指標(biāo)，確定的沖擊周期將會導(dǎo)致一定偏差，本文利用包絡(luò)譜稀疏度作為沖擊周期的確定指標(biāo)，因?yàn)槠溽槍Φ木褪墙饩矸e效果的分析，恰恰與最大相關(guān)峭度解卷積方法的原理相吻合，故而可以確定出較為準(zhǔn)確的沖擊周期參數(shù)。

圖5 小波Shannon熵確定的包絡(luò)功率譜圖

同時隨著對最大相關(guān)峭度解卷積方法的研究不斷深入，得出結(jié)論為濾波器長度在50以下，其對原信號基本無降噪效果，文獻(xiàn)[8]將濾波器長度選擇范圍確定為[2，500]，步長每次增加1，相比本文確定的范圍為[100，500]，步長每次增加10，存在搜尋效率下降的問題。因此，根據(jù)上述分析說明了本文方法的優(yōu)勢。

為了說明本文選擇最佳沖擊周期時初定的濾波器長度和經(jīng)過排列熵選擇的最佳濾波器長度在降噪效果上的差異，利用[150，33]的組合與最佳參數(shù)組合[310，33]進(jìn)行實(shí)驗(yàn)對比，觀察圖6和圖4可知，[150，33]組合的包絡(luò)功率譜圖相對[310，33]組合的包絡(luò)功率譜圖的干擾成分要多，且對故障信息提取造成了影響，因此進(jìn)一步說明了只有確定最佳的濾波器長度，MCKD方法才能發(fā)揮最佳功用，本文確定的初始濾波器長度只是為了防止部分噪聲對最佳周期的確定造成干擾。本實(shí)驗(yàn)具有普遍性，因篇幅限制實(shí)驗(yàn)信號的分析不再贅述。

圖6 濾波器長度為150的包絡(luò)功率譜圖

2.4 實(shí)驗(yàn)信號分析

利用軸承內(nèi)圈故障信號對最優(yōu)參數(shù)MCKD方法對做進(jìn)一步驗(yàn)證，內(nèi)圈故障特征頻率為fi=124 Hz，根據(jù)上述選擇方法確定合理周期和濾波器長度，圖7為沖擊周期與包絡(luò)譜稀疏度關(guān)系示意圖。沖擊周期T為22時包絡(luò)譜稀疏度達(dá)到最大，因此確定的最佳周期為22。同理，根據(jù)排列熵值和濾波器長度的關(guān)系，確定的濾波器長度為440。

圖8為內(nèi)圈排列熵值與濾波器長度的關(guān)系。

內(nèi)圈振動時域信號如圖9所示。峭度值為6.49。軸承內(nèi)圈信號經(jīng)過最優(yōu)參數(shù)MCKD方法處理后得到時域和包絡(luò)功率譜如圖10所示。峭度值為12.44，小波Shannon熵確定MCKD參數(shù)(確定的沖擊周期為26，濾波器長度為437)方法的包絡(luò)功率譜如圖11所示。對比圖10、11，分析發(fā)現(xiàn)圖10中包絡(luò)功率譜中故障成分有效得到提取。圖11中因沖擊周期選擇不恰當(dāng)，包絡(luò)功率譜中激發(fā)出的沖擊成分與內(nèi)圈故障信號頻率成分存在一定誤差，并且其倍頻周圍存在較多干擾成分，因此最優(yōu)參數(shù)MCKD方法相比小波Shannon熵確定MCKD參數(shù)方法選擇的沖擊周期更加合理，使得特定故障特征提取更加有效。

圖7 內(nèi)圈包絡(luò)譜稀疏度與沖擊周期的關(guān)系

Fig.7 Relationship between envelope spectrum sparsity and impact period of inner-race fault

圖8 內(nèi)圈排列熵值與濾波器長度的關(guān)系

圖9 內(nèi)圈信號的時域波形(峭度值=6.49)

(a) 時域圖

(b) 包絡(luò)功率譜圖

Fig.10 The time domain graph and envelope power spectrum of the optimal parameter MCKD (Kurtosis value=12.44)

圖11 小波Shannon熵確定的包絡(luò)功率譜圖

觀察圖7和圖8可知，在峰值22和410處的附近有相差不大的較高幅值出現(xiàn)，為了驗(yàn)證所出現(xiàn)幅值處的濾波器長度值和沖擊周期值對信號降噪和沖擊成分的提取作用與本文方法提取的最佳參數(shù)組合對信號的作用不一致，利用最佳參數(shù)MCKD方法處理的包絡(luò)功率譜圖與之比較，說明本文方法的有效性和穩(wěn)定性。

觀察圖7和圖8確定的組合值為[410，30]和較遠(yuǎn)處的幅值較高的[320，54]，觀察兩組值確定的MCKD方法作用的包絡(luò)功率譜圖12(a)、圖12(b)可知，雖然兩組參數(shù)都有一定的降噪效果，干擾成分變少，但是相對于最佳濾波器長度作用的包絡(luò)功率譜圖，降噪效果并不佳，同時可以看到因?yàn)榇_定的周期值與內(nèi)圈信號包含的故障沖擊成分無關(guān)，包絡(luò)功率譜圖并未提取到想要的故障信息。

(a)

(b)

因此，通過上述的實(shí)驗(yàn)說明，MCKD方法需要的參數(shù)信息對應(yīng)的是峰值處的濾波器長度和沖擊周期值，與其他較大幅值處對應(yīng)的值無關(guān)，進(jìn)一步說明了本文方法的有效性。

3 最優(yōu)參數(shù)MCKD與ELMD診斷軸承復(fù)合故障

總體局部均值分解(Ensemble Local Mean Decomposition，ELMD)方法[14]能夠有效抑制模態(tài)混疊，是一種噪聲輔助分析方法。即首先將設(shè)定好的有限幅值白噪聲信號加入到待分解的目標(biāo)信號中；然后采用 LMD方法對加入白噪聲后的混合信號進(jìn)行分解，得到多個PF分量和一個殘余分量，以上過程中每次添加不同的白噪聲信號；最后將分解得到的多個相應(yīng)PF分量的均值作為最終的分解結(jié)果。ELMD消除模態(tài)混疊的原理在于白噪聲具有均勻分布的特性，當(dāng)一個信號加入白噪聲后，將被 LMD 方法分解到由白噪聲確定的濾波器組當(dāng)中和其相關(guān)的通頻帶中，這樣便會有效地減少模態(tài)混淆現(xiàn)象。另外由于數(shù)量足夠多的白噪聲其均值接近于零，于是對目標(biāo)信號通過多次添加不同白噪聲進(jìn)行LMD 分解，然后將多次分解結(jié)果的總體平均值作為最終的結(jié)果，這樣可以保留真實(shí)的分解結(jié)果。

當(dāng)滾動軸承發(fā)生復(fù)合故障時，復(fù)合故障中出現(xiàn)不同的故障特征頻率，其故障特征頻譜和時間尺度不同，利用傳統(tǒng)的整體頻譜分析將復(fù)合故障表現(xiàn)在同一頻譜圖中，會因?yàn)槟芰枯^大，信號較強(qiáng)的成分存在，導(dǎo)致能量小，信號弱的故障特征被淹沒[15]。故而，對復(fù)合故障的分析，首先應(yīng)該對原始信號進(jìn)行降噪處理，對不同的信號特征進(jìn)行分離，通過頻域分析有針對的提取某一故障。對比單一故障的診斷方法，對復(fù)合故障診斷方法還需進(jìn)一步研究。

本文以兩種故障并存為例，可把故障特征明顯，容易提取的故障稱為主故障，其次分離出來的故障命名為次故障。因?yàn)檩S承出現(xiàn)復(fù)合故障時，故障沖擊周期不同，輸入不同的沖擊周期成分，即可提取不同的故障特征。根據(jù)前面分析，MCKD方法具有不同沖擊周期對不同故障敏感性不同和有效降噪的特點(diǎn)，能夠逐個將不同的故障特征從復(fù)合故障信號中提取出來，因此可將MCKD方法引入到軸承的復(fù)合故障診斷中，以達(dá)到分離復(fù)合故障的目的。為了達(dá)到最佳的降噪效果，有效提取信號中的全部故障特征，本文提出最優(yōu)參數(shù)MCKD和ELMD方法診滾動軸承復(fù)合故障。

首先將原始信號進(jìn)行最優(yōu)參數(shù)MCKD處理，既能有效降噪，又能利用合理參數(shù)確定能量大，信號強(qiáng)的主故障信息，并運(yùn)用包絡(luò)功率譜提取主故障的故障特征頻率；然后利用ELMD方法將目標(biāo)信號分解為由高到低的信號頻率成分構(gòu)成的多個PF分量，將各本征模態(tài)函數(shù)進(jìn)行頻譜分析，將得到的分量信號中包含主故障信息的某一分量剔除，防止對次故障濾波器長度和沖擊周期兩個參數(shù)的確定造成干擾。

將其他分量重構(gòu)的信號再次利用最優(yōu)參數(shù)MCKD方法進(jìn)行分析，確定次故障的最佳參數(shù)，并再進(jìn)行包絡(luò)功率譜處理，即可識別和診斷次故障的存在情況，基于此復(fù)合故障診斷流程圖，如圖13所示。

圖13 最優(yōu)參數(shù)MCKD與ELMD診斷復(fù)合故障流程

Fig.13 Process of bearing compound fault diagnosis based on the optimal parameter MCKD and ELMD

根據(jù)以上分析，利用最優(yōu)參數(shù)MCKD對復(fù)合故障信號中的故障參數(shù)進(jìn)行選擇，可將復(fù)合故障信號中的故障信息進(jìn)行逐一提取，而將其與ELMD分解方法結(jié)合診斷軸承復(fù)合故障，可以有效解決特征明顯的故障信號掩蓋微弱故障頻譜特征的問題并有效分離復(fù)合故障，以達(dá)到診斷多種故障并存的目的。

4 應(yīng)用分析

為了驗(yàn)證最優(yōu)參數(shù)MCKD與ELMD診斷軸承復(fù)合故障方法的有效性，本文采用軸承與轉(zhuǎn)子故障模擬實(shí)驗(yàn)臺模擬了軸承外圈和內(nèi)圈的復(fù)合故障，以及軸承外圈和滾動體的復(fù)合故障實(shí)驗(yàn)，實(shí)驗(yàn)系統(tǒng)由故障診斷實(shí)驗(yàn)臺、壓電式加速度傳感器、數(shù)據(jù)采集儀、筆記本電腦組成。軸承型號為ER-10K,結(jié)構(gòu)參數(shù)如表1所示。傳感器安裝在軸承座的軸向、垂直徑向和水平徑向這三個測試點(diǎn)上，如圖14所示。設(shè)定實(shí)驗(yàn)電機(jī)轉(zhuǎn)速為1 380 r/min，采樣頻率為2 560 Hz，通過計(jì)算得到外圈故障特征頻率為69.92 Hz，內(nèi)圈故障特征頻率為114.08 Hz，滾動體故障特征頻率為90.62 Hz。

表1 ER-10K軸承的結(jié)構(gòu)參數(shù)

圖14 軸承與轉(zhuǎn)子故障模擬實(shí)驗(yàn)臺

(1)滾動軸承外圈和內(nèi)圈復(fù)合故障

圖15為滾動軸承內(nèi)圈和外圈復(fù)合故障的原始時域波形圖和頻譜圖，時域波形圖具有一定的周期性沖擊成分，說明滾動軸承存在故障，而從頻譜圖中并未看到明顯的故障特征信息，同時伴隨著明顯的干擾成分。

對原始信號進(jìn)行最優(yōu)參數(shù)MCKD處理，首先利用包絡(luò)譜稀疏度確定較明顯的沖擊成分的周期值，并運(yùn)用排列熵篩選合適的濾波器長度，主故障確定的沖擊周期為37、濾波器長度為330，處理后的信號時域波形，如圖16所示。從圖16可知，噪聲成分明顯減少，周期成分變的明顯，對振動信號進(jìn)行主故障的分析和提取，其包絡(luò)功率譜如圖17所示?？芍?，頻率成分69.90 Hz及其倍頻，頻率成分與計(jì)算得到的外圈故障特征基頻及其倍頻成分相近，因此，可以判斷出滾動軸承復(fù)合故障中包含軸承外圈故障，故而利用本文方法可成功提取出反映軸承運(yùn)行狀態(tài)的信號特征頻率。

(a)

(b)

圖16 主故障時域波形

圖17 主故障包絡(luò)功率譜

對振動信號進(jìn)行ELMD處理，得到的模態(tài)分量信號如圖18(a)所示。其各分量的包絡(luò)譜圖如圖18(b)所示。根據(jù)包絡(luò)譜圖的頻率信息，將包含69.90 Hz頻率成分的PF2分量濾除，可以避免主故障沖擊成分妨礙次故障沖擊周期的確定，并且達(dá)到提取次故障的目的。

(a) 模態(tài)分量信號

(b) 各分量的包絡(luò)譜

其他分量進(jìn)行信號重構(gòu)，重構(gòu)后的信號再次運(yùn)用最優(yōu)參數(shù)MCKD處理，提取次故障信息，利用包絡(luò)譜稀疏度再次確定較明顯的沖擊成分的周期值為23，并運(yùn)用排列熵篩選合適的濾波器長度為380，次故障的時域波形如圖19所示。處理后的信號周期性沖擊成分突出，將其作為分析信號進(jìn)行后續(xù)處理，包絡(luò)功率譜如圖20所示。從圖19和圖20可知，頻率成分114 Hz，這與軸承內(nèi)圈故障特征頻率114.08 Hz相近，故可判斷復(fù)合故障中包含軸承內(nèi)圈故障。對比圖17和圖20，內(nèi)圈包絡(luò)功率譜圖的幅值遠(yuǎn)沒有外圈幅值高，說明了內(nèi)圈振動信號是被淹沒的弱沖擊成分，通過最優(yōu)參數(shù)MCKD參數(shù)和ELMD方法的結(jié)合，突出了內(nèi)圈的振動沖擊信號。經(jīng)過上述分析表明，本文方法可以有效診斷和分離滾動軸承復(fù)合故障。

圖19 次故障時域波形

圖20 次故障包絡(luò)功率譜

(2)滾動軸承外圈和滾動體復(fù)合故障

同樣觀察原始信號的時域波形圖和頻譜圖，從時域圖21(a)可知，基本看不到故障特征信息，而經(jīng)過傳統(tǒng)的傅里葉變換，圖21(b)中也沒有出現(xiàn)體現(xiàn)故障信息的故障特征頻率。

(a)

(b)

利用上文方法采用排列熵值和包絡(luò)譜稀疏度對原始信號進(jìn)行處理，首先確定的濾波器長度為330、沖擊周期值為37，與上文提到的外圈參數(shù)一致，故確定外圈故障為主故障，對原始信號利用上述確定的參數(shù)進(jìn)行MCKD處理并進(jìn)行故障特征提取，從圖22、23可知，時域波形較為干凈，瞬態(tài)沖擊成分突出，其包絡(luò)功率譜譜峰主要由69 Hz及其倍頻成分構(gòu)成，上述成分接近表征

圖22 主故障時域波形

外圈故障的特征頻率69.92 Hz及其倍頻成分，故可判斷出此復(fù)合故障包括外圈故障信息，為主故障。

圖23 主故障包絡(luò)功率譜

對振動信號利用ELMD分解方法進(jìn)行處理，各分量的包絡(luò)譜如圖24所示。從圖24可知，主故障基頻信息在PF1分量中，將PF1分量進(jìn)行濾除，防止其對次故障參數(shù)的確定造成干擾。

圖24 各分量的包絡(luò)譜

將其他分量重構(gòu)的信號利用上文方法進(jìn)行處理，確定次故障最佳濾波器參數(shù)和沖擊周期參數(shù)分別為320和28，再利用最優(yōu)參數(shù)MCKD方法進(jìn)行降噪處理，并突出次故障沖擊成分，(見圖25和圖26)。觀察圖26可知，信號的包絡(luò)功率譜主要頻率成分為91 Hz和其倍頻成分，接近滾動體故障特征頻率，故而可判斷復(fù)合故障中包括滾動體故障信息。上述實(shí)驗(yàn)數(shù)據(jù)分析進(jìn)一步表明，最優(yōu)參數(shù)MCKD和ELMD方法結(jié)合提取軸承復(fù)合故障效果是明顯的。

圖25 次故障時域波形

圖26 次故障包絡(luò)功率譜

5 結(jié) 論

本文針對MCKD中的沖擊周期T和濾波器長度L的選擇問題展開了研究，提出了最優(yōu)參數(shù)的選擇方法，然后將最優(yōu)參數(shù)MCKD和ELMD方法結(jié)合對滾動軸承的復(fù)合故障進(jìn)行了研究。通過對實(shí)驗(yàn)信號的分析，得出以下結(jié)論:

(1) 利用包絡(luò)譜稀疏度和排列熵可以很好的確定MCKD中的重要參數(shù)，從而使信號能夠有效降噪，提取期望得到的沖擊成分。

(2) 利用ELMD方法具有抑制模態(tài)混疊的特點(diǎn)，根據(jù)復(fù)合故障信號的特點(diǎn)，對其進(jìn)行分解，可以有效確定主故障頻率基頻成分所在的分量，方便信號做后續(xù)處理。

(3) 對MCKD算法設(shè)定不同的沖擊周期時，激發(fā)的頻率成分不同，根據(jù)這一原理將最優(yōu)參數(shù)MCKD方法引入到滾動軸承的復(fù)合故障診斷中去，提取不同的故障頻率成分。

(4) 沖擊成分明顯的分量會掩蓋信號弱的頻率成分，因此利用最優(yōu)參數(shù)MCKD和ELMD方法診斷復(fù)合故障時，應(yīng)首先分離主故障，再識別次故障。這樣可有效和準(zhǔn)確的將復(fù)合故障中的故障信息全部提取出來。