陳芳

摘要：滾動(dòng)軸承發(fā)生故障時(shí)，其信號(hào)具有非平穩(wěn)性特點(diǎn)，因此本文介紹了基于Hilbert-Huang變換的信號(hào)處理方法。該方法以EMD為基礎(chǔ)，得到由高頻到低頻的IMF分量，并對(duì)每個(gè)IMF分量進(jìn)行時(shí)頻譜分析，最終確定故障類型。由實(shí)驗(yàn)可知，該方法分析外圈故障得到的故障頻率與理論計(jì)算值相符，在滾動(dòng)軸承故障診斷方面具有很強(qiáng)的實(shí)踐性。

Abstract： When the rolling bearing fails， its signal has non-stationary characteristics. Therefore， this paper introduces the signal processing method based on Hilbert-Huang transform. Based on EMD， the IMF components from high frequency to low frequency are obtained， The time-frequency spectrum of each IMF component is analyzed to determine the fault type. The experimental results show that the fault frequency from the outer ring obtained by this method is consistent with the theoretical calculation value， and it has strong practicality in fault diagnosis of rolling bearing.

關(guān)鍵詞：滾動(dòng)軸承;故障分析;Hilbert-Huang變換

Key words： rolling bearing;fault analysis;Hilbert-Huang transform

中圖分類號(hào)：TH165.3? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-957X（2020）24-0178-02

0? 引言

就滾動(dòng)軸承這一零件的應(yīng)用特性來講，有易損壞的特征。根據(jù)相關(guān)調(diào)查數(shù)據(jù)資料信息統(tǒng)計(jì)結(jié)果來看，機(jī)械設(shè)備在運(yùn)行期間發(fā)生的故障問題，有70%的幾率在于振動(dòng)故障，而從這70%的振動(dòng)故障中，大約有30%的故障原因在于軸承，就此，本文對(duì)滾動(dòng)軸承的狀態(tài)進(jìn)行檢測(cè)和故障診斷具有現(xiàn)實(shí)意義。軸承故障診斷最關(guān)鍵的是從滾動(dòng)軸承故障振動(dòng)信號(hào)中提取故障特征和狀態(tài)識(shí)別。滾動(dòng)軸承發(fā)生故障時(shí)，故障信號(hào)為非平穩(wěn)信號(hào)，那么對(duì)其進(jìn)行故障診斷的關(guān)鍵就在于從非平穩(wěn)信號(hào)中提取故障頻率。Hilbert-Huang變換可以對(duì)非線性非平穩(wěn)信號(hào)進(jìn)行處理，并且在工程實(shí)踐中診斷準(zhǔn)確率很高。

1? Hilbert-Huang變換原理

Hilbert-Huang變換主要是由美國(guó)國(guó)家宇航局Norden E. Huang提出的新型信號(hào)分析方法，該信號(hào)分析方法的整個(gè)過程包括：①將采集到的信號(hào)，用經(jīng)驗(yàn)?zāi)B(tài)分解法（EMD）分解成若干個(gè)固有模態(tài)函數(shù)（簡(jiǎn)稱IMF），這樣任何一個(gè)信號(hào)就可以被分解為有限個(gè)IMF分量之和;②對(duì)每個(gè)IMF分量進(jìn)行Hilbert變換，求出每個(gè)IMF分量隨時(shí)間變化的瞬時(shí)頻率與瞬時(shí)幅值;③對(duì)每個(gè)IMF的頻率和幅值構(gòu)造Hilbert譜，并進(jìn)行分析。HHT是一種新的能夠有效處理非平穩(wěn)信號(hào)的時(shí)頻分析方法，這種方法對(duì)信號(hào)能進(jìn)行自適應(yīng)的時(shí)頻分解，并具有很高的時(shí)域和頻域分辨率。

1.1 經(jīng)驗(yàn)?zāi)B(tài)分解法

對(duì)信號(hào)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解是HHT變換的核心，EMD是借助時(shí)間序列上下包絡(luò)平均值來對(duì)瞬時(shí)平衡位置進(jìn)行確定，然后再將非平穩(wěn)信號(hào)進(jìn)行分解，使其形成一組穩(wěn)態(tài)和線性數(shù)據(jù)序列集，也就是固有模態(tài)函數(shù)（IMF）。運(yùn)用EMD方法把一個(gè)復(fù)雜的非平穩(wěn)信號(hào)分解為有限個(gè)IMF分量之和，其具體的步驟如下：

①確定信號(hào)中所有局部極大值和極小值，形成上下包絡(luò)線，上下包絡(luò)線的均值為m1，原始信號(hào)與m1的差值定義為分量h1，h1=x（t）-m1，若分量h1為一個(gè)固有模態(tài)函數(shù)，則h1就是這個(gè)信號(hào)中x（t）的第一個(gè)分量。

②若分量h1沒有滿足IMF所提出的相應(yīng)條件，可將h1作為初始數(shù)據(jù)，然后對(duì)上述步驟重復(fù)，以此來獲取上下包絡(luò)線平均值m11，然后對(duì)h11=h1-m11是否符合上述固有模態(tài)函數(shù)所提出的要求和條件進(jìn)行判斷，若還是不滿足的情況下，則需要循環(huán)重復(fù)k次，最終獲取到h1k=h1 （k-1）-m1k，使得h1k滿足IMF的條件，記c1（t）=h1k，則c1（t）為信號(hào)x（t）第1個(gè)滿足IMF條件的分量，代表原數(shù)據(jù)序列中的高頻部分。

③將c1（t）從信號(hào)x（t）中分離出來，得到r1=x （t）-c1（t）將r1作為初始數(shù)據(jù)，然后對(duì)步驟①～③循環(huán)重復(fù)進(jìn)行，從而獲取到信號(hào)x（t）的第2個(gè)滿足固有模態(tài)函數(shù)條件分量c2（t），然后繼續(xù)進(jìn)行n次循環(huán)，直到獲取的第n個(gè)信號(hào)滿足固有模態(tài)函數(shù)所提出的條件分量，則：

當(dāng)rn成為一個(gè)單調(diào)函數(shù)不能再?gòu)闹刑崛M足IMF條件的分量時(shí)，循環(huán)結(jié)束。

其中rn殘余分量，表示信號(hào)的平均趨勢(shì);ci（t）表示信號(hào)的各個(gè)IMF分量。因此可以把任何一個(gè)信號(hào)分解為n個(gè)滿足IMF條件的分量和一個(gè)殘量之和，從c1（t），c2（t）直到cn（t），信號(hào)由高至低不同頻率段組成，同時(shí)帶寬也不等[1]。

1.2 Hilbert譜

EMD方法，簡(jiǎn)單來講就是將信號(hào)進(jìn)行自適應(yīng)分解，然后將各個(gè)IMF分量相加求和，將其中各個(gè)IMF分量進(jìn)行相應(yīng)的變換，促使瞬時(shí)頻率概念有實(shí)際存在意義，最終獲取Hilbert譜[2]。因此對(duì)上面公式中的ci? （t）進(jìn)行Hilbert變換后，得到的公式結(jié)果為：

2? 滾動(dòng)軸承故障分析

滾動(dòng)軸承由內(nèi)圈、外圈、滾動(dòng)體、保持架組成。大部分時(shí)候外圈處于靜止不動(dòng)，內(nèi)圈跟隨軸承狀態(tài)進(jìn)行轉(zhuǎn)動(dòng)。在滾動(dòng)軸承中，滾動(dòng)體作為其中核心部分的元件，能夠?qū)⑤S承原本的滑動(dòng)摩擦狀態(tài)轉(zhuǎn)變?yōu)闈L動(dòng)摩擦，摩擦質(zhì)量對(duì)軸承使用壽命有直接性影響。基于以上，若滾動(dòng)軸承在運(yùn)行期間出現(xiàn)故障情況下，損傷部位與之接觸面相互碰撞會(huì)產(chǎn)生低頻成分特征頻率，這類信號(hào)具有周期性，其中頻率可運(yùn)用相應(yīng)的公式進(jìn)行計(jì)算，另外對(duì)于故障發(fā)生位置不同，其特征頻率也會(huì)不同?；诖?，可通過頻譜特征頻率定位故障類型并進(jìn)行相應(yīng)診斷。假設(shè)外圈處于固定狀態(tài)，內(nèi)圈轉(zhuǎn)動(dòng)速度、軸承節(jié)圓直徑、滾動(dòng)體直徑、接觸角以及滾動(dòng)體個(gè)數(shù)分別為n、D、d、α、Z，通過以上，最終獲取到的各類故障的特征頻率為[3]：

當(dāng)滾動(dòng)軸承在運(yùn)行期間發(fā)生損傷性的故障事故情況下，通過分析振動(dòng)信號(hào)頻譜，尋找其中比較大的峰值頻率，即為故障特征頻率。

本次實(shí)驗(yàn)采用加速度傳感器對(duì)6205-JEMSKF 深溝型軸承外圈故障進(jìn)行分析，其中徑為1.537英寸，接觸角為90，滾子直徑為0.3126英寸，滾子個(gè)數(shù)為9，轉(zhuǎn)速為1730，采樣頻率為12KHz，計(jì)算該軸承外圈故障特征頻率為103.4Hz。當(dāng)外圈發(fā)生故障時(shí)，其振動(dòng)信號(hào)的時(shí)域波形圖如圖1所示。

經(jīng)過降噪后，振動(dòng)信號(hào)進(jìn)行EMD分解，得到8個(gè)IMF， 并對(duì)所有IMF進(jìn)相關(guān)系數(shù)、峭度值進(jìn)行相應(yīng)的計(jì)算，通過對(duì)計(jì)算結(jié)果篩選后，將其中存在故障的IMF1找出來，然后通過對(duì)其進(jìn)行希爾伯特變換，即可得到圖2的故障包絡(luò)普。

就圖2所示：通過Hilbert變換后，IMF1包絡(luò)譜中存在非常明顯的故障發(fā)生特征頻率，故障發(fā)生的頻率為104Hz，忽略軸承本身和安裝誤差，與理論計(jì)算數(shù)值相符，即可將滾動(dòng)軸承外圈故障問題檢測(cè)出來。

3? 結(jié)論

HHT將采集到的信號(hào)進(jìn)行EMD分解，把原始信號(hào)分解成幾個(gè)IFM之和，然后對(duì)其中所有IMF分量進(jìn)行Hilbert 變換，然后獲取Hilbert譜和Hilbert邊際譜，最終確定滾動(dòng)軸承的故障類型。從整體上來看，HHT擺脫以往傳統(tǒng)模式下Fourier在變化時(shí)的局限性約束，將非平穩(wěn)信號(hào)中頻率成分準(zhǔn)確分離出來，在實(shí)際應(yīng)用中具有一定的實(shí)用性優(yōu)勢(shì)。根據(jù)實(shí)驗(yàn)研究，將這種方法應(yīng)用于軸承故障檢測(cè)以及故障診斷作業(yè)中，具有效率高、準(zhǔn)確率高的特點(diǎn)。

參考文獻(xiàn)：

[1]程軍圣，于德介，楊宇.基于EMD和奇異值分解技術(shù)的滾動(dòng)軸承故障診斷方法[J].數(shù)據(jù)采集與處理，2004（5）：90-95.

[2]任玥.基于Hilbert Huang變換的滾動(dòng)軸承智能診斷方法的研究[D].成都：西南交通大學(xué)，2007：11-17.

[3]張灝繽.基于振動(dòng)分析的滾動(dòng)軸承故障診斷研究[D].遼寧科技大學(xué)，2015：11-12.