（江蘇省特種設(shè)備安全監(jiān)督檢驗研究院江陰分院 江陰 214400）

隨著國民經(jīng)濟和工業(yè)的發(fā)展，起重設(shè)備朝著大型化和高速化發(fā)展，滾動軸承作為起升機構(gòu)電機中的重要部件，承擔(dān)著傳遞轉(zhuǎn)矩和支承整個轉(zhuǎn)子系統(tǒng)的重要作用，由于經(jīng)常受到頻繁起制動和高速運轉(zhuǎn)的影響，經(jīng)常超負(fù)荷運行導(dǎo)致?lián)p壞。調(diào)查發(fā)現(xiàn)大約有30%的旋轉(zhuǎn)機械安全事故是由軸承故障引發(fā)的[1]，所以建立一套有關(guān)軸承故障檢測系統(tǒng)顯得十分重要。軸承的破壞一般從表面損傷開始，由破壞表面與其他接觸體之間的摩擦形成的脈沖振蕩信號判斷軸承故障的形式與位置，是一項十分有效的故障檢測方法[2]。特種設(shè)備用電機由于其工作特點的復(fù)雜性，環(huán)境噪聲將對振動信號的采集和傳輸產(chǎn)生影響，因此采用何種降噪技術(shù)對檢驗精度的提高有重要意義[3]。目前常用的兩種降噪技術(shù)為FFT降噪和小波去噪，F(xiàn)FT非平信號濾波降噪會導(dǎo)致一部分尖端信號被當(dāng)作噪聲去除，且無法去除特定頻率信號，而基于小波去噪的軸承故障檢測雖然避免了此誤差，但由于其按指數(shù)間隔進(jìn)行時頻分解，對高頻段的分辨率低于檢驗所需標(biāo)準(zhǔn)。因此本文基于小波包能量法對電機軸承實施信號處理和故障診斷，以獲得較高的精度和檢驗效率。

1 小波包與Hilbert 分析

小波變換由于其在時頻分解方法上具有很強的局部優(yōu)化特點，因此在非平穩(wěn)故障信號的提取上越來越多地被使用，只是在高頻信號區(qū)域頻率和低頻時域的分辨率較差[4]。故本文采用的小波包分析法通過對頻帶進(jìn)行多層次分解，并由不同頻帶信號特征展開信號頻譜的匹配，以解決分辨率不足的問題[5-7]。小波包能量法的原理如圖1所示，其中：U0為采集的滾動軸承源信號；A，D分別表示信號的低頻與高頻區(qū)域。在一個完整的小波包分解中，將小波子空間Wj與尺度空間Vj相結(jié)合起用一個新的子空間來表達(dá)：

圖1 三層小波包分解樹

小波包能量譜即采用此方法最終得到的能量分布，根據(jù)Parseval原理，故障軸承的原始采集信號以及經(jīng)過小波包變換之后在能量層具有如下關(guān)系：

式中：

f（x)——原始采集信號；

C（j,k)——經(jīng)小波包能量分解后第j層第k個節(jié)點的幅值。

由上式可知，將各個頻信號作如下的平方加運算，即可求得指定頻帶中的小波包能量譜能量表達(dá)式Ej，k：

式中：

N——被采集信號長度值大??；

所有的Ej，k之和所構(gòu)成的小波能量譜E為：

軸承表面的損傷特征不同，其與其他元件之間在旋轉(zhuǎn)時產(chǎn)生的振動沖擊信號也就不同，據(jù)此可判斷故障所屬構(gòu)件以及故障破壞類型。高速旋轉(zhuǎn)下故障引起的高頻間歇性沖擊帶來軸承系統(tǒng)的固有振動，因此信號面臨解調(diào)問題。為了提高故障檢測的精度，需要對最終重構(gòu)的故障振動特征信號進(jìn)行包絡(luò)解調(diào)。Hilbert基于包絡(luò)解調(diào)的信號變化方法，可快速得到檢測所需的瞬態(tài)時域、頻域、相位圖等特征。若軸承檢測中發(fā)現(xiàn)某一信號x（t)，則針對其的Hilbert變換可表示如下：

則x（t)的解析信號r（t)可按下式求出：

A（t)為r（t)的幅值，也即x（t)的包絡(luò)：

最終根據(jù)被求解析信號x（t)可得到調(diào)制信號即相應(yīng)故障的特征信號詳細(xì)信息，此時的頻率則可視為起重機減速器滾動軸承的發(fā)生故障時的特征頻率。

2 試驗方案

美國凱斯西儲大學(xué)對軸承故障進(jìn)行試驗研究并提供各種詳盡的試驗數(shù)據(jù)庫，本文采用橋式起重機常用滾動軸承型號6205故障數(shù)據(jù)。信號采集過程中將振動傳感器置于所用電機的驅(qū)動端，采集時的電機實際轉(zhuǎn)速n為1772 rpm，采樣頻率設(shè)為12kHz。軸承故障信號的不同代表著所首先破壞的構(gòu)件不同，依據(jù)內(nèi)圈、滾動體和外圈的位置可提前計算出不同構(gòu)件發(fā)生故障時的固有特征頻率，計算式如下：

式中：

f0——外圈故障的特征頻率；

f1——內(nèi)圈故障的特征頻率；

f——轉(zhuǎn)軸的轉(zhuǎn)動頻率（f=s/60，s為轉(zhuǎn)軸轉(zhuǎn)速)；

Z——滾動體個數(shù)；

d——滾動體直徑；

B——軸承節(jié)徑。

根據(jù)上式帶入6205型軸承的相關(guān)參數(shù)計算得內(nèi)圈的固有特征頻率為159.96 Hz，外圈則為105.89 Hz，滾動體特征故障頻率為70.55Hz，保持架的特征故障頻率最低為11.92Hz，在MATLAB中針對此軸承的正常工作、內(nèi)外圈及滾動體發(fā)生故障時的非正常信號作數(shù)據(jù)輸出，其時域波形結(jié)果如圖2～圖5所示。橫坐標(biāo)為時域信號的時間，縱坐標(biāo)則為壓電式采集傳感器輸出的振動信號值。

圖2 正常運行時的時域

圖3 內(nèi)圈故障時的時域信號

圖4 外圈故障時的時域信號

由圖2～圖5可以看出：相比正常工作的軸承信號，起重機故障滾動軸承的輸出信號時域波形明顯較為密集，且波形峰值規(guī)律有所不同。具體表現(xiàn)為：當(dāng)軸承內(nèi)圈發(fā)生故障時，波形峰值分布較為分散且與其他振動壓電信號差別不大；軸承外圈故障時，波形分布更加密集，峰值較為稀疏且易于分辨；當(dāng)滾動體發(fā)生故障時，時域波形表現(xiàn)為只有一個明顯的波峰，且最大峰值較小，只有正常壓電信號的兩倍。其在頻域上，故障振幅的表現(xiàn)也有相同特征且高頻特征明顯，如圖5所示?；谛〔ò芰糠纸夥ǖ钠鹬貦C滾動軸承檢測流程圖如圖6所示。

圖5 滾動體故障時的時域信號

圖6 檢驗方案流程圖

3 檢測結(jié)果與故障分析

為更加清晰地對軸承故障類別進(jìn)行判斷，利用小波包分解法進(jìn)行能量譜圖的繪制，采用db2小波對以上分析的起重機滾動軸承正常工作、內(nèi)圈故障、外圈故障及滾動體故障幾種采集信號進(jìn)行能量的集中頻段提取表達(dá)，第三層的分解結(jié)果如圖7～圖10所示。

圖7 正常振動的能量分布

圖8 內(nèi)圈故障振動的能量分布

圖9 外圈故障振動的能量分布

圖10 軸承滾動體故障信號的能量分布

由圖7～圖10可知，不同構(gòu)件故障時表現(xiàn)出的能量分布有很大差異，軸承無故障正常運行時前四個低頻區(qū)域大約占據(jù)了總能量的90%以上，這主要歸因于周期性振動引發(fā)的規(guī)律性動態(tài)響應(yīng)。而當(dāng)各構(gòu)件發(fā)生故障時，其中若干單個頻帶的能量明顯高于其他值，此例中頻帶7最為顯著，主要是由于當(dāng)各構(gòu)件發(fā)生點蝕、磨損等損傷時帶來的周期性沖擊將破壞原來的低頻振動，以致引起支承轉(zhuǎn)子對沖擊的高頻響應(yīng)和固有頻率的改變。對比圖8和圖9、圖10還可發(fā)現(xiàn)，于頻帶7能量而言，軸承內(nèi)圈和滾動體發(fā)生故障時較外圈故障時其值稍低，這是因為當(dāng)內(nèi)圈和滾動體之間的面面接觸摩擦造成故障脈沖力信號能量損失，因此能量分布也會顯現(xiàn)出幅度降低趨勢。由此可知，利用小波包能量法可以有效地對起重機用軸承的運行狀態(tài)及故障原因進(jìn)行檢測和判斷。

接下來采用能量熵標(biāo)準(zhǔn)對滾動軸承的低頻到高頻的振動故障信號進(jìn)行采集和系數(shù)重構(gòu)，歸一化處理后的結(jié)果如圖11所示，重構(gòu)了各分解節(jié)點（3，0），（3，1），（3，2），（3，3），（3，4），（3，5），（3，6），（3，7）所對應(yīng)的第三層對應(yīng)節(jié)點各個頻段能量值總和。由圖11可知，軸承故障節(jié)點主要分布在節(jié)點（3，6）處，故障軸承的能量分布較為集中且高頻，對重構(gòu)信號作EMD分解，求得9個IMF分量和1個殘余量，將各IMF分量與采集特征頻率相結(jié)合后，結(jié)果見表1。

圖11 故障頻段能量圖

由表1可知IMF1的相關(guān)系數(shù)值最大，故該分量應(yīng)該最能準(zhǔn)確地反映出故障信號特征，因此選擇IMF1做Hilbert包絡(luò)譜分析，圖12為IMF1分量的Hilbert頻譜圖，最大峰值頻率分布在154.35Hz附近，最接近6205型軸承之前計算得到的內(nèi)圈的固有特征頻率為f=159.96Hz，故可以判斷此時滾動軸承為內(nèi)圈故障。相比理論計算值，利用小波包檢測的故障特征頻率的相對誤差為3.5%，誤差較小，并不會影響對故障類型的判別，因此結(jié)果有一定的可信度。

表1 IMF分量的相關(guān)系數(shù)

雖然本文利用的是實驗室數(shù)據(jù)，但將此方法對實際起重機起升機構(gòu)的減速器滾動軸承進(jìn)行檢測，可以對支承電機轉(zhuǎn)子的滾動軸承的運行狀態(tài)預(yù)測和故障甄別提高準(zhǔn)確度和效率，勢必也將提高電機及起重機整體安全性能和降低事故發(fā)生率。

圖12 IMF1分量的 Hilbert頻譜圖

4 結(jié)論

利用小波包能量法對滾動軸承展開運行狀態(tài)監(jiān)測和故障原因判斷，通過Hilbert頻譜分析對比采集信號的能量分布形式和特征頻率，可以有效檢驗軸承的正常運行、內(nèi)圈、外圈及滾動體三類構(gòu)件發(fā)生故障的表現(xiàn)，進(jìn)而找出故障位置，可為起重機減速器用滾動軸承的安全監(jiān)測和故障判斷提供一種行之有效的方法，也可用于電機的壽命預(yù)測和整機安全評估。