王冬云，張文志，張建剛

（1.燕山大學(xué) 軋機(jī)研究所，河北 秦皇島 066004；2.秦皇島職業(yè)技術(shù)學(xué)院，河北 秦皇島 066100；3.秦皇島首鋼板材有限公司，河北 秦皇島 066003）

滾動軸承是機(jī)械中最常用的部件之一，其運(yùn)行情況直接影響到主機(jī)的機(jī)能。振動作為一種信息因子直接預(yù)示著運(yùn)行狀態(tài)的正常與否。因而對軸承進(jìn)行振動監(jiān)測和診斷是可行的，且必將獲得較大的經(jīng)濟(jì)效益［1］。

當(dāng)軸承系統(tǒng)發(fā)生故障時(shí)，振動信號能量的空間分布與正常系統(tǒng)的振動信號相比會發(fā)生相應(yīng)的變化，即振動信號的能量改變包含著豐富的故障特征信息。小波變換的時(shí)域和頻域的局部化和可變分辨率的特點(diǎn)使得用它分析瞬變信號具有比傳統(tǒng)Fourier分析更為顯著的優(yōu)點(diǎn)［2］。目前，小波變換法被廣泛應(yīng)用于軸承故障診斷系統(tǒng)。該方法的主要缺點(diǎn)是對高頻部分的頻率分辨率相對較低。而軸承出現(xiàn)故障時(shí)，高頻帶分解卻又是至關(guān)重要的，因而往往難以得到滿意的識別效果。

小波包分析能為信號提供一種更精細(xì)的分析方法。其對信號的低頻和高頻部分都進(jìn)行分解，可對非平穩(wěn)和突變信號進(jìn)行精確的特征提取，更有效地反映信號的時(shí)頻特征。文獻(xiàn)［3-4］將小波包分析方法應(yīng)用于軸承故障診斷，并建立了人工神經(jīng)網(wǎng)絡(luò)診斷模型，展示了小波包分析在軸承故障特征提取的有效性。但軸承運(yùn)行狀態(tài)從正常到異常是一個(gè)漸變過程，所表現(xiàn)出來的征兆在很多情況下具有模糊性，神經(jīng)網(wǎng)絡(luò)模型不能夠?qū)ζ涮卣餍盘柕哪：匦猿浞煮w現(xiàn)。而以模糊數(shù)學(xué)為理論基礎(chǔ)的聚類分析方法為解決這類模糊性問題提供了新的途徑［5］。鑒于此，采用基于小波包能量譜的滾動軸承故障特征提取方法，并提出應(yīng)用模糊聚類法進(jìn)行診斷。試驗(yàn)證明，此方法能有效地對軸承進(jìn)行狀態(tài)監(jiān)測與故障診斷。

1 小波包能量譜特征量提取

1.1 小波包分解

3層小波包分解樹如圖1所示，其中，S為原始信號，H和G分別為正交共軛低通濾波器和高通濾波器，通過小波包分解后可以把信號無泄漏、不重疊地分解到相鄰的獨(dú)立頻段上。分解后得到第3層8個(gè)頻帶的小波包分解系數(shù)，每個(gè)頻帶的頻率范圍為：［（n-1）2-j-1fs，n2-j-1fs］，n＝1，2，3，...8，其中fs為采樣頻率。

圖1 小波包分解結(jié)構(gòu)圖

可見，小波包分解［6］可以以更高的分辨率展示信號在不同頻段的時(shí)域信息，為準(zhǔn)確提取信號特征提供了強(qiáng)有力的工具。

1.2 能量譜特征量提取

由于軸承出現(xiàn)故障時(shí)會對各頻帶內(nèi)信號能量有較大的影響，根據(jù)不同頻段內(nèi)能量的分布情況可以診斷出故障的類型。從小波包分解結(jié)果中提取頻段能量作為信號特征可有效反映軸承運(yùn)行狀態(tài)的變化。這種以能量方式表示的小波包分解結(jié)果稱為小波包能量譜。因此，在小波包能量譜中，可以選取各頻帶內(nèi)信號的平方和作為能量的標(biāo)志，為進(jìn)一步的故障識別提供了診斷的前提條件。

以3層小波包分解為例，說明故障信號小波包能量譜特征量提取方法［7］：

（1）首先對振動能量信號進(jìn)行3層小波包分解，從而得到第3層從低頻到高頻8個(gè)子頻帶的小波包分解系數(shù)

（4）構(gòu)造特征向量。由于系統(tǒng)出現(xiàn)故障時(shí)會對各頻帶內(nèi)信號能量有較大的影響，因此，以能量為元素可以構(gòu)造一個(gè)特征向量。特征向量T構(gòu)造如下：

根據(jù)各頻帶能量所占比重的變化即可實(shí)現(xiàn)對軸承運(yùn)行狀態(tài)的監(jiān)測。

2 試驗(yàn)及數(shù)據(jù)處理分析

2.1 試驗(yàn)裝置

為了驗(yàn)證提出的診斷方法的有效性，進(jìn)行了軸承故障診斷試驗(yàn)研究。根據(jù)故障診斷理論模型，建立了由直流電動機(jī)、齒輪箱、采集系統(tǒng)、信號分析系統(tǒng)組成的試驗(yàn)系統(tǒng)見圖2。通過該試驗(yàn)平臺，可以模擬齒輪箱軸承的各種故障狀態(tài)，包括軸承的內(nèi)圈、外圈故障，并采集故障信號。用于試驗(yàn)的軸承型號為6406，通過安裝在齒輪減速器中軸承座外側(cè)的加速度傳感器，測取軸承處于不同故障時(shí)的加速度信號。

圖2 故障診斷試驗(yàn)裝置

2.2 試驗(yàn)步驟及結(jié)果

首先，用電火花蝕刻技術(shù)在一套軸承的內(nèi)滾道和另一套軸承的外滾道上分別加工直徑為4 mm的孔，產(chǎn)生類似凹痕的表面損傷來模擬軸承內(nèi)圈故障與外圈故障。啟動試驗(yàn)平臺，軸承內(nèi)圈轉(zhuǎn)速為900 r／min，而外圈保持靜止。采樣頻率為fs＝5 000 Hz，采樣時(shí)間為0.248 s。正常工況與故障軸承的時(shí)域波形信號如圖3所示。

圖3 軸承振動信號時(shí)域波形

2.3 信號處理

根據(jù)小波包能量譜特征量提取方法，計(jì)算出各個(gè)頻段的能量，并以此作為特征參數(shù)，依次排列做為特征向量。在此分別提取正常工況、外圈故障和內(nèi)圈故障3種運(yùn)行狀態(tài)的特征向量各3組，作為診斷的標(biāo)準(zhǔn)樣本，如表2所示。

表1 3層小波包分解后節(jié)點(diǎn)對應(yīng)的頻率范圍

表2 不同運(yùn)行狀態(tài)的特征向量

利用小波包分解提取的能量在不同頻帶內(nèi)數(shù)值相差也很大，必須對原始數(shù)據(jù)進(jìn)行無量綱化處理，使每一指標(biāo)值統(tǒng)一于某種共同的數(shù)據(jù)特性范圍。本案采用極值歸一化公式把數(shù)據(jù)壓縮到［0，1］之間。極值歸一化公式為

式中：X為原始數(shù)據(jù)；Xmax為原始數(shù)據(jù)最大值；Xmin為原始數(shù)據(jù)最小值。

歸一化后軸承正常和故障時(shí)特征向量對比如圖4所示。從圖4可以看出，軸承正常和故障運(yùn)行狀態(tài)的特征向量有很大的差別。正常工況時(shí)，振動加速度信號的能量主要分布在低頻段，如頻帶2和頻帶4中，這是由周期性振源引起的響應(yīng)。當(dāng)軸承的內(nèi)、外圈出現(xiàn)故障時(shí)，振動信號中包含了相應(yīng)的沖擊成分，產(chǎn)生脈沖沖擊響應(yīng)，由于脈沖沖擊是一種瞬態(tài)激振，將引起軸承系統(tǒng)固有頻率的共振［8］。軸承故障產(chǎn)生的沖擊響應(yīng)特征在高頻段表現(xiàn)十分明顯，這時(shí)能量主要分布在頻段5、頻段6和頻段8。故障運(yùn)行時(shí)，軸承內(nèi)圈故障和外圈故障能量分布也有所區(qū)別：軸承外圈出現(xiàn)故障時(shí)頻段8的高頻能量激增值大；而當(dāng)軸承內(nèi)圈出現(xiàn)故障時(shí)頻段8的高頻能量激增值相對較小。這是由于當(dāng)滾動體與內(nèi)圈損傷點(diǎn)接觸時(shí)產(chǎn)生的脈沖力要通過滾動體與外圈的界面?zhèn)鞑ズ蟛抛饔糜谕馊Γa(chǎn)生了能量損失，所以脈沖力幅度肯定要比前者小得多。軸承不同的故障沖擊將會導(dǎo)致不同的振動響應(yīng)信號，也就是說，不同的故障類型對應(yīng)有不同的頻帶能量的分布特征，這表明應(yīng)用小波包分解提取信號的能量分布信息作為信號的特征向量能很好地表征軸承的運(yùn)行狀態(tài)。

圖4 軸承正常和故障時(shí)特征向量對比圖

同理，將待診斷的軸承運(yùn)行狀態(tài)信號進(jìn)行特征量提取，作為待診斷樣本，最后將9個(gè)標(biāo)準(zhǔn)樣本和1個(gè)待診斷樣本共同組成診斷樣本集：

提取的不同運(yùn)行狀態(tài)的特征向量矩陣E為：

式中：E1，E2，E3為正常運(yùn)行狀態(tài)的標(biāo)準(zhǔn)樣本；E4，E5，E6為內(nèi)圈故障運(yùn)行狀態(tài)的標(biāo)準(zhǔn)樣本；E7，E8，E9為外圈故障運(yùn)行狀態(tài)的標(biāo)準(zhǔn)樣本；Ed為待診斷樣本。

3 故障診斷

采用模糊聚類分析方法［9-10］完成故障診斷。首先采用絕對值減數(shù)法對矩陣E進(jìn)行標(biāo)定，其定義為：

式中：rst為相似系數(shù)，表示樣本Es與Et之間的親疏關(guān)系；C為系數(shù)，應(yīng)適當(dāng)選取，為了使0≤rst≤1，取C＝0.308。求出診斷樣本集E中任意兩個(gè)樣本之間的相似系數(shù)，并組成模糊相似矩陣R，再根據(jù)傳遞閉包法得到診斷樣本集的模糊等價(jià)矩陣R＝R4：

采用λ截矩陣法對模糊等價(jià)矩陣R進(jìn)行聚類分析。對于任一λ∈［0，1］，分類時(shí)Rλ中各元素改為這樣取值：大于或等于λ的元素都改取為1，小于λ的元素都改取為0，根據(jù)Rλ中1與0的排列情況進(jìn)行分類。結(jié)果表明：λ∈（0.850 65，0.858 83）時(shí)分類和故障診斷的結(jié)果最為準(zhǔn)確。如λ＝0.858 8，得到λ截矩陣R0.8588為：

可見，截矩陣R0.8588將診斷樣本集分為3類：｛E1，E2，E3｝，｛E4，E5，E6｝，和｛E7，E8，E9，Ed｝，其中待診斷樣本Ed屬于軸承外圈故障運(yùn)行狀態(tài)。

4 結(jié)束語

針對滾動軸承故障診斷的高頻帶分解至關(guān)重要的特點(diǎn)，提出一種把信號無泄漏、不重疊地分解到相鄰的獨(dú)立頻段上的小波包能量譜技術(shù)，以提取軸承運(yùn)行狀態(tài)特征。同時(shí)針對軸承運(yùn)行狀態(tài)所表現(xiàn)出來的征兆具有模糊性的特點(diǎn)，提出一種基于模糊聚類的故障診斷方法。試驗(yàn)表明，該方法獲得的結(jié)果與軸承運(yùn)行狀態(tài)有很好的一致性，利用此方法能有效地監(jiān)測軸承運(yùn)行狀態(tài)，對故障類型做出準(zhǔn)確判斷，效果理想。