陳 娜， 楊紹普， 潘存治

(1.石家莊鐵道大學(xué) 信息科學(xué)與技術(shù)學(xué)院，石家莊 050043；2.石家莊鐵道大學(xué) 機(jī)械工程學(xué)院，石家莊 050043； 3.北京交通大學(xué) 機(jī)械與電子控制工程學(xué)院，北京 100044)

基于匹配追蹤預(yù)處理的鐵道車(chē)輛滾動(dòng)軸承故障診斷方法

陳 娜1，3， 楊紹普2， 潘存治2

(1.石家莊鐵道大學(xué) 信息科學(xué)與技術(shù)學(xué)院，石家莊 050043；2.石家莊鐵道大學(xué) 機(jī)械工程學(xué)院，石家莊 050043； 3.北京交通大學(xué) 機(jī)械與電子控制工程學(xué)院，北京 100044)

提出一種基于匹配追蹤預(yù)處理的軸承故障診斷算法，避免了傳統(tǒng)共振解調(diào)方法對(duì)帶通濾波的經(jīng)驗(yàn)依賴，消除了傳統(tǒng)方法在處理環(huán)節(jié)產(chǎn)生的互調(diào)干擾頻率成分，能夠精確提取損傷軸承的故障頻率。通過(guò)在輪對(duì)跑合實(shí)驗(yàn)臺(tái)上進(jìn)行試驗(yàn)驗(yàn)證，結(jié)果顯示比傳統(tǒng)共振解調(diào)方法有一定優(yōu)勢(shì)，可以更有效提取輪對(duì)軸承的損傷故障。

匹配追蹤；故障診斷；互調(diào)；輪對(duì)軸承

輪對(duì)軸承的工作狀況對(duì)鐵道車(chē)輛的安全運(yùn)行影響很大。而對(duì)于軸承故障診斷，最直觀的信息是沖擊脈沖發(fā)生頻率。但是由于實(shí)際系統(tǒng)的復(fù)雜性，直接提取沖擊脈沖發(fā)生頻率幾乎是不可能的。共振解調(diào)方法可以通過(guò)采集軸承元件的固有振動(dòng)信號(hào)，提取到?jīng)_擊脈沖的發(fā)生頻率，以此進(jìn)行軸承故障診斷。

傳統(tǒng)共振解調(diào)方法是用帶通濾波濾除噪聲和干擾，保留有意義的固有振動(dòng)信號(hào)。實(shí)際中各元件的固有振動(dòng)頻率往往變化很大，而且當(dāng)脈沖寬度不是足夠窄時(shí)，沖擊脈沖將難以引起頻率較高的固有振動(dòng)[1]。因此，振動(dòng)系統(tǒng)傳遞規(guī)律的復(fù)雜性，使得對(duì)帶通濾波的頻率設(shè)置很大程度上依賴于工程經(jīng)驗(yàn)，這種依賴性限制了共振解調(diào)技術(shù)的普及[2]。

對(duì)于固有振動(dòng)信號(hào)的提取，頻域上要求嚴(yán)格對(duì)應(yīng)各元件固有振動(dòng)頻率，時(shí)域上要求嚴(yán)格對(duì)應(yīng)沖擊脈沖發(fā)生時(shí)刻。因此對(duì)于軸承故障振動(dòng)信號(hào)的處理而言，在時(shí)域和頻域均有較高的分辨率要求。小波變換、Wigner分布等方法雖然在時(shí)頻域作了一些分辨率上的權(quán)衡[3-4]，但仍不能滿足實(shí)際輪對(duì)系統(tǒng)軸承故障診斷對(duì)特征提取的要求。

匹配追蹤算法采用過(guò)完備字典替代正交分解對(duì)信號(hào)進(jìn)行表示，尤其適合于描述沖擊波形等非穩(wěn)態(tài)信號(hào)[5]。文獻(xiàn)[6]采用非參數(shù)波形估計(jì)方法，將得到的一組波形對(duì)采集信號(hào)進(jìn)行表示，可以對(duì)信號(hào)和噪聲進(jìn)行有效分離。文獻(xiàn)[7]依據(jù)軸承轉(zhuǎn)速、尺寸等參數(shù)建立字典模型，并且提出了逐次改變特性參數(shù)的方法，在保證原子與實(shí)際信號(hào)能夠較好匹配的前提下，大幅度提升了字典的使用效率。文獻(xiàn)[8]利用改進(jìn)型K-SVD算法構(gòu)造自適應(yīng)字典，在稀疏分解的階段，選用峭度值最大的原子對(duì)信號(hào)進(jìn)行逼近，最終基于重構(gòu)信號(hào)的包絡(luò)譜分析獲取故障類型。

如果提取的沖擊振動(dòng)信號(hào)中包含多個(gè)固有振動(dòng)頻率，在解調(diào)處理中會(huì)產(chǎn)生互調(diào)成分，影響故障頻率判別，而以往的相關(guān)文獻(xiàn)并未對(duì)此進(jìn)行分析。本文首先對(duì)互調(diào)問(wèn)題進(jìn)行說(shuō)明，并提出相應(yīng)的解決辦法。然后簡(jiǎn)單介紹匹配追蹤算法，并描述該算法實(shí)現(xiàn)固有振動(dòng)信號(hào)提取的處理過(guò)程。最后將本文算法應(yīng)用于輪對(duì)跑合實(shí)驗(yàn)臺(tái)軸承故障信號(hào)的處理，并進(jìn)行對(duì)比分析。結(jié)果表明，本文方法可以有效提取輪對(duì)軸承的損傷故障。

1 互調(diào)成分的產(chǎn)生原理

假定傳感器采集到的是兩種元件(比如軸承外圈和軸承座)的振動(dòng)信號(hào)，根據(jù)滾動(dòng)軸承損傷類故障振動(dòng)機(jī)理，可以認(rèn)為兩者受到?jīng)_擊脈沖的時(shí)間間隔相同，振動(dòng)信號(hào)傳遞到傳感器的路徑不同，可以認(rèn)為引起振動(dòng)的沖擊幅度不同。此外，阻尼衰減特性的不同只會(huì)影響解調(diào)后故障頻率諧波的相對(duì)幅度，不會(huì)產(chǎn)生新的頻率分量。因此為了分析的方便，可以認(rèn)為兩者在沖擊脈沖的激勵(lì)下產(chǎn)生的脈沖調(diào)制信號(hào)波形類似，相差僅為幅度，分別為A(t)和αA(t)。兩者的固有振動(dòng)頻率分別為f1和f2，由沖擊脈沖激勵(lì)產(chǎn)生的固有振動(dòng)信號(hào)表示成A(t)ej2πf1t和αA(t)ej2πf2t。對(duì)混合振動(dòng)信號(hào)進(jìn)行包絡(luò)解調(diào)的過(guò)程如式(1)所示[9]，其中A*(t)代表A(t)的復(fù)共軛?？梢钥闯?，結(jié)果中前兩項(xiàng)為脈沖調(diào)制幅度信號(hào)的疊加，后兩項(xiàng)為兩種固有振動(dòng)信號(hào)的互調(diào)產(chǎn)物。

[A(t)ej2πf1t+αA(t)ej2πf2t][A*(t)e-j2πf1t+

(1)

軸承系統(tǒng)各元件的固有振動(dòng)雖然是高頻信號(hào)，但并不能保證不同元件固有振動(dòng)頻率的差值仍為一個(gè)高頻。如果軸承系統(tǒng)包含多個(gè)相同元件，由于實(shí)際工況或加工等原因，這些元件更可能有不一樣但相近的固有振動(dòng)頻率。如此一來(lái)，互調(diào)成分就會(huì)在低頻段產(chǎn)生干擾，從而影響故障頻率判別。

因此，更為合理的包絡(luò)解調(diào)過(guò)程應(yīng)該先把不同頻率的共振信號(hào)分離開(kāi)來(lái)，然后分別解調(diào)，最后將解調(diào)結(jié)果求和。這樣既保證了多個(gè)元件的脈沖調(diào)制幅度信號(hào)的疊加，提升了故障信號(hào)的信噪比，又避免了異頻共振信號(hào)互調(diào)對(duì)低頻故障信號(hào)的干擾。

2 匹配追蹤原理

2.1正交基與超完備原子庫(kù)

人們可以針對(duì)類似的概念準(zhǔn)確而方便地表達(dá)出細(xì)微的差別，得益于語(yǔ)言系統(tǒng)中包含大量意義相近的詞匯。若非如此，單調(diào)的詞匯只能表達(dá)大概的含義。想要表達(dá)精確的含義，則必須附加大量的解釋。也就是說(shuō)，冗余度的減小，使得利用詞匯來(lái)傳遞信息的能力被“稀釋”[10-11]。

信號(hào)處理中，傳統(tǒng)的正交分解會(huì)由于這種“冗余度不足”而造成信號(hào)分解不具有簡(jiǎn)潔性，進(jìn)而導(dǎo)致后續(xù)處理的復(fù)雜甚至得到錯(cuò)誤的信號(hào)解釋。Mallat和Zhang提出了應(yīng)用超完備冗余字典對(duì)信號(hào)進(jìn)行稀疏分解的思想，并引入了匹配追蹤算法。作為一種新的信號(hào)表示理論，匹配追蹤算法利用超完備的冗余函數(shù)系統(tǒng)代替?zhèn)鹘y(tǒng)的正交基函數(shù)，采用自適應(yīng)方法將信號(hào)分解為一系列基本信號(hào)(即原子)的線性組合，利用原子庫(kù)的高度冗余性保證信號(hào)分解能夠更準(zhǔn)確地反映信號(hào)的內(nèi)在屬性。

2.2匹配追蹤算法

匹配追蹤算法是一種逐步迭代取得信號(hào)稀疏表達(dá)的貪婪算法，具體流程描述如下：

H表示Hilbert空間，待處理信號(hào)f∈H。D={gγ}γ∈Г為過(guò)完備原子庫(kù)，gγ=1。為逼近f，首先在原字庫(kù)D中選擇一個(gè)原子gγ0，使其能最大限度匹配信號(hào)f結(jié)構(gòu)，即：

〈f,gγ0〉=sup〈f,gγ〉

(2)

信號(hào)分解成如下形式：

f=〈f,gγ0〉gγ0+R1f

(3)

〈f,gγ0〉gγ0是f對(duì)gγ0的投影，R1f是用gγ0表示信號(hào)f后產(chǎn)生的殘余。顯然gγ0與R1f是正交的，可得：

(4)

對(duì)殘余R1f進(jìn)行類似的分解，得到：

R1f=〈R1f,gγ1〉gγ1+R2f

(5)

其中

〈f,gγ1〉=sup〈f,gγ〉

(6)

依次類推，匹配追蹤就是不斷將殘余信號(hào)在最匹配的原子上進(jìn)行投影，然后得到新的殘余信號(hào)。假設(shè)已完成k次投影，殘余信號(hào)Rkf與原子庫(kù)中g(shù)γk匹配度最好，則：

Rkf=〈Rkf,gγk〉gγk+Rk+1f

(7)

gγk與Rkf是正交的，即：

(8)

當(dāng)信號(hào)殘余能量小于門(mén)限時(shí)，迭代結(jié)束，此時(shí)已經(jīng)對(duì)信號(hào)進(jìn)行了n次分解：

(9)

式(9)為信號(hào)f在原字庫(kù)D中n階逼近的表達(dá)式，Rnf為信號(hào)進(jìn)行n階逼近后的殘余誤差，誤差能量會(huì)隨分解的進(jìn)行迅速減小。對(duì)于信號(hào)f∈H滿足有限長(zhǎng)度條件時(shí)，Mallat已證明當(dāng)?shù)螖?shù)n趨于無(wú)窮時(shí)，Rnf將按指數(shù)收斂于0。

3 基于匹配追蹤預(yù)處理的軸承故障診斷

3.1沖擊字典的構(gòu)造

當(dāng)滾動(dòng)軸承的表面有損傷故障時(shí)，損傷部分與其它元件表面接觸會(huì)產(chǎn)生突變的沖擊脈沖力，進(jìn)而引起軸承系統(tǒng)高頻固有振動(dòng)頻率的諧振。這種由于損傷故障而產(chǎn)生的沖擊振動(dòng)信號(hào)主要涉及兩部分信息：一是沖擊產(chǎn)生的時(shí)刻，如果等間隔連續(xù)沖擊則表現(xiàn)為故障通過(guò)頻率；二是沖擊作用引發(fā)的系統(tǒng)高頻固有振動(dòng)頻率，如軸承外圈的徑向彎曲固有振動(dòng)或是軸承座的固有振動(dòng)。此外，沖擊振動(dòng)信號(hào)還具有單邊性和指數(shù)衰減的特點(diǎn)。因此可以選定以下函數(shù)作為沖擊字典的基元函數(shù)[12]。

(10)

式中：pn表現(xiàn)了沖擊響應(yīng)的阻尼衰減特性，un給出了沖擊響應(yīng)事件發(fā)生的初始時(shí)刻，fn是系統(tǒng)中不同元件的固有頻率，Kn保證了基元函數(shù)具有單位能量。

3.2故障信息提取

(11)

根據(jù)式(10)，gn(t)不僅包含沖擊脈沖發(fā)生時(shí)刻的信息，還包含元件固有振動(dòng)頻率的信息。直接對(duì)式(11)進(jìn)行包絡(luò)解調(diào)，雖然可以保留下沖擊脈沖發(fā)生時(shí)刻的信息，但由于第1節(jié)分析結(jié)果，會(huì)產(chǎn)生各元件固有振動(dòng)頻率的互調(diào)成分，在包絡(luò)信號(hào)頻譜低頻段產(chǎn)生干擾，影響故障頻率判別。本文對(duì)沖擊振動(dòng)成分的處理過(guò)程如下所示。

首先分別求取各固有振動(dòng)成分gn(t)的包絡(luò)信號(hào)：

xn(t)={[angn(t)]×[angn(t)]*}1/2

(12)

然后再將各包絡(luò)信號(hào)xn(t)求和：

(13)

最后根據(jù)和信號(hào)y(t)的頻譜判斷故障頻率。

4 輪對(duì)跑合實(shí)驗(yàn)臺(tái)采集信號(hào)處理與分析

4.1實(shí)驗(yàn)設(shè)備

為驗(yàn)證本文方法的有效性，基于輪對(duì)跑合實(shí)驗(yàn)臺(tái)進(jìn)行了故障軸承振動(dòng)信號(hào)采集和處理。實(shí)驗(yàn)臺(tái)采用液壓控制系統(tǒng)，利用橡膠輪帶動(dòng)被測(cè)輪對(duì)的輪緣進(jìn)行跑合，如圖1所示。振動(dòng)信號(hào)采集存儲(chǔ)設(shè)備包括CA-YD-189型壓電式加速度傳感器和INV36DF型信號(hào)采集儀，其中加速度傳感器磁座吸附在機(jī)架上，如圖2所示。通過(guò)掃頻試驗(yàn)，得到傳感器的諧振頻率為5.9 kHz。采集設(shè)備采樣頻率為51.2 kHz，采樣時(shí)長(zhǎng)為2 s。

圖1 輪對(duì)跑合實(shí)驗(yàn)臺(tái)

圖2 傳感器布置

實(shí)驗(yàn)軸承為鐵路貨運(yùn)列車(chē)較常用的197726型雙列圓錐滾子軸承，主要參數(shù)如表1所示。故障類型為外圈剝落，故障軸承實(shí)物圖如圖3所示。實(shí)驗(yàn)中輪對(duì)轉(zhuǎn)速為465 r/min，理論計(jì)算可得到外圈故障特征頻率為fo=67.4 Hz。

表1 197726型軸承的主要參數(shù)

圖3 故障軸承實(shí)物圖

4.2傳統(tǒng)共振解調(diào)處理與分析

采集到的振動(dòng)加速度信號(hào)時(shí)域波形如圖4所示，圖5為10 kHz以下能量按頻率的分布圖。觀察時(shí)域波形或頻域分布，均不能獲得故障頻率的信息。

將同一工況下，使用無(wú)故障軸承采集到的振動(dòng)信號(hào)功率譜(本文未給出)與圖5作比較，均可以在1 000～2 000 Hz之間頻段發(fā)現(xiàn)較強(qiáng)的振動(dòng)信號(hào)。對(duì)故障軸承運(yùn)轉(zhuǎn)情況下采集信號(hào)進(jìn)行1 000～2 000 Hz的帶通濾波，然后進(jìn)行包絡(luò)解調(diào)處理，得到500 Hz以下的包絡(luò)譜如圖6所示?？梢钥闯?，該頻段信號(hào)的包絡(luò)能量主要集中在40 Hz以下，與軸承外圈故障頻率無(wú)關(guān)。結(jié)合其他分析，本文認(rèn)為該頻段的振動(dòng)主要是由于輪對(duì)轉(zhuǎn)動(dòng)引起的機(jī)架等設(shè)備的振動(dòng)，對(duì)于故障提取來(lái)說(shuō)是干擾信號(hào)。

圖4 采集信號(hào)時(shí)域波形

圖5 采集信號(hào)功率譜

輪對(duì)軸承轉(zhuǎn)動(dòng)時(shí)，滾子每通過(guò)一次外圈損傷點(diǎn)，就會(huì)產(chǎn)生一次沖擊振動(dòng)。因此，外圈和滾子的固有振動(dòng)信號(hào)中應(yīng)當(dāng)承載最豐富的故障信息。實(shí)踐證明，工作時(shí)軸承各元件的固有振動(dòng)頻率很難準(zhǔn)確估計(jì)。根據(jù)載荷、材料和物理尺寸等情況進(jìn)行理論估算[13-14]，197726型軸承外圈和滾子的固有振動(dòng)頻率大致范圍分別是3 000～7 000 Hz和2 000～4 000 Hz。

圖6 主要干擾成分包絡(luò)譜

將采集信號(hào)進(jìn)行2 000～7 000 Hz的帶通濾波，然后進(jìn)行包絡(luò)解調(diào)，500 Hz以下的包絡(luò)譜如圖7所示。圖7的包絡(luò)譜中雖然可以找到故障頻率fo=67.4 Hz，但附近存在大量干擾，最高譜線位置與轉(zhuǎn)頻fr(7.75 Hz)非常接近，而且故障頻率二次諧波分量過(guò)大，不利于故障分析。

圖7 帶通濾波后包絡(luò)譜(5 000 Hz帶寬)

傳感器在安裝諧振頻率附近有較高的靈敏度，因此對(duì)采集信號(hào)按照5 900 Hz為中心，分別做2 000 Hz和500 Hz帶寬的帶通濾波，然后進(jìn)行包絡(luò)解調(diào)，分別得到包絡(luò)譜如圖8和圖9所示。對(duì)比外圈故障頻率處譜線能量的絕對(duì)值，兩種濾波策略的處理結(jié)果變化不大，但較窄濾波后的包絡(luò)譜中干擾成分稍顯降低。即便如此，圖9所示的包絡(luò)解調(diào)結(jié)果仍不利于故障分析。

4.3本文方法處理過(guò)程

下面使用本文方法進(jìn)行數(shù)據(jù)處理。首先仍舊對(duì)采集信號(hào)進(jìn)行2 000～7 000 Hz的帶通濾波，然后使用匹配追蹤算法提取各個(gè)脈沖沖擊成分，最后按照式(12)和式(13)所述過(guò)程進(jìn)行包絡(luò)解調(diào)。

圖8 帶通濾波后包絡(luò)譜(2 000 Hz帶寬)

圖9 帶通濾波后包絡(luò)譜(500 Hz帶寬)

本次試驗(yàn)中，un搜索步進(jìn)為10個(gè)采樣點(diǎn)，fn的搜索范圍為[2 000 Hz，7 000 Hz]，步進(jìn)為20 Hz。固有振動(dòng)阻尼通常比較大，各個(gè)衰減振動(dòng)基本上是獨(dú)立的，因此把原子的支撐區(qū)間長(zhǎng)度設(shè)置為沖擊發(fā)生間隔的三分之一。由于帶通濾波后信號(hào)中仍存在大量非沖擊振動(dòng)成分，因此迭代終止的條件不能設(shè)置為殘余信號(hào)能量小于某門(mén)限。本文的迭代終止條件如(14)式所示，即如果本次搜索結(jié)果的能量遠(yuǎn)小于之前每次搜索結(jié)果的平均能量時(shí)，迭代結(jié)束。

(14)

將提取出的沖擊成分按照式(11)的形式進(jìn)行求和，直接解調(diào)得到包絡(luò)譜如圖10所示。將提取出的沖擊成分分別按照式(12)求包絡(luò)，再按照式(13)將包絡(luò)求和，得到和信號(hào)的頻譜如圖11所示。

4.4本文方法處理結(jié)果分析

4.4.1 濾波處理的影響分析

相對(duì)于圖10，圖7的包絡(luò)譜中包含大量干擾信號(hào)。這是因?yàn)榻庹{(diào)前的帶通濾波器較寬，原本不是因?yàn)閾p傷故障沖擊引起的振動(dòng)信號(hào)在解調(diào)后落入了低頻段。使用匹配追蹤處理較好的屏蔽了非沖擊振動(dòng)信號(hào)，因此圖10包絡(luò)譜中的干擾成分大大減少。

圖10 直接解調(diào)得到的包絡(luò)譜

圖11 本文方法得到的包絡(luò)譜

觀察圖8和圖9中信號(hào)能量的絕對(duì)幅度，看到故障頻率處信號(hào)能量比圖7和圖10的少了大概10倍。這是因?yàn)檩^窄的帶通濾波使得原本承載沖擊振動(dòng)的信號(hào)被濾掉了。而匹配追蹤的處理過(guò)程使得解調(diào)前數(shù)據(jù)可以具有較寬的帶寬，即保證了盡量豐富的固有振動(dòng)信號(hào)承載了更多的沖擊脈沖成分，進(jìn)而保證了故障信息的完整性。

4.4.2 信號(hào)互調(diào)的影響分析

圖11是對(duì)提取出來(lái)的每個(gè)沖擊成分分別求包絡(luò)得到的總包絡(luò)譜，圖10是將全部沖擊成分混合信號(hào)直接求包絡(luò)得到的包絡(luò)譜。對(duì)比圖11，圖10的頻譜中出現(xiàn)了一些干擾成分，根據(jù)第1節(jié)的分析，某些干擾源自于不同頻率沖擊信號(hào)的互調(diào)成分。本文在匹配追蹤處理時(shí)使用了20 Hz的頻率搜索步進(jìn)，導(dǎo)致了提取出的各個(gè)沖擊信號(hào)的頻率差集中在20 Hz附近及其倍頻上。圖10中20 Hz處譜線較為明顯，是由故障包絡(luò)信號(hào)的直流成分調(diào)制20 Hz的互調(diào)成分產(chǎn)生的，107 Hz處的信號(hào)是由故障頻率調(diào)制40 Hz的互調(diào)成分產(chǎn)生的。圖11的結(jié)果驗(yàn)證了對(duì)每個(gè)沖擊成分分別求包絡(luò)可以避免互調(diào)產(chǎn)生，同時(shí)不丟失沖擊信號(hào)的包絡(luò)信息，使得故障頻率明顯，干擾減少，利于故障診斷。

5 結(jié) 論

本文通過(guò)對(duì)轉(zhuǎn)動(dòng)軸承損傷故障作用機(jī)理的簡(jiǎn)單介紹，說(shuō)明了共振解調(diào)方法對(duì)提取固有振動(dòng)信號(hào)的處理要求。在此基礎(chǔ)上，本文提出了基于匹配追蹤預(yù)處理的軸承故障診斷算法，并詳細(xì)說(shuō)明了如何準(zhǔn)確提取固有振動(dòng)信號(hào)以及避免解調(diào)過(guò)程產(chǎn)生互調(diào)干擾的處理過(guò)程。通過(guò)在輪對(duì)跑合實(shí)驗(yàn)臺(tái)上進(jìn)行實(shí)采信號(hào)驗(yàn)證，本文方法在對(duì)帶通濾波設(shè)置規(guī)律經(jīng)驗(yàn)依賴的規(guī)避、沖擊脈沖信號(hào)提取的完整性以及故障頻率附近干擾剔出等方面，均比傳統(tǒng)共振解調(diào)方法具有一定優(yōu)勢(shì)。

[1] 梅宏斌. 滾動(dòng)軸承振動(dòng)監(jiān)測(cè)與診斷 理論·方法·系統(tǒng)[M]. 北京：機(jī)械工業(yè)出版社，1995.

[2] 劉金朝，丁夏完，王成國(guó). 自適應(yīng)共振解調(diào)法及其在滾動(dòng)軸承故障診斷中的應(yīng)用[J]. 振動(dòng)與沖擊，2007, 26(1): 38-41.

LIU Jinzhao, DING Xiawan, WANG Chengguo. Adaptive resonance demolation method and its application to fault diagnosis of freight carrolling bearings[J]. Journal of Vibration and Shock, 2007, 26(1): 38-41.

[3] 張進(jìn)，馮志鵬，褚福磊. 滾動(dòng)軸承故障特征的時(shí)間—小波能量譜提取方法[J]. 機(jī)械工程學(xué)報(bào)，2011，47(17): 44-49.

ZHANG Jin, FENG Zhipeng, CHU Fulei. Extraction of rolling bearing fault feature based on time-wavelet energy spectrum[J]. Journal of Mechanical Engineering, 2011, 47(17): 44-49.

[4] 王詩(shī)彬，許佳，朱忠奎. 瞬態(tài)成分參數(shù)的最小二乘法辨識(shí)及其軸承故障特征提取應(yīng)用[J]. 機(jī)械工程學(xué)報(bào)，2012, 48(7):68-76.

WANG Shibin, XU Jia, ZHU Zhongkui. LSM-based transient parameter identification and its application in feature extraction of bearing fault[J]. Journal of Mechanical Engineering,2012,48(7):68-76.

[5] MALLAT S G, ZHANG Z. Matching pursuits with time-frequency dictionaries[J]. IEEE Transactions on Signal Processing, 1993, 41(12): 3377-3415.

[6] 范虹，孟慶豐，張優(yōu)云,等. 基于改進(jìn)匹配追蹤算法的特征提取及其應(yīng)用[J]. 機(jī)械工程學(xué)報(bào)，2007,43(7):115-119.

FAN Hong, MENG Qingfeng, ZHANG Youyun,et al. Feature extraction based on modified matching pursuit and its application[J]. Journal of Mechanical Engineering, 2007, 43(7): 115-119.

[7] 崔玲麗，王婧，鄔娜，等. 基于軸承故障信號(hào)特征的自適應(yīng)沖擊字典匹配追蹤方法及應(yīng)用[J].振動(dòng)與沖擊，2014，33(11):54-60.

CUI Lingli, WANG Jing, WU Na, et al. Bearing fault diagnosis based on self-adaptive impulse dictionary matching pursuit[J]. Journal of Vibration and Shock, 2014, 33(11):54-60.

[8] 余發(fā)軍，周鳳星，嚴(yán)?？? 基于字典學(xué)習(xí)的軸承早期故障稀疏特征提取[J].振動(dòng)與沖擊，2016,35(6):181-186.

YU Fajun, ZHOU Fengxing, YAN Baokang. Bearing initial fault feature extraction via sparse representation based on dictionary learning[J]. Journal of Vibration and Shock, 2016,35(6):181-186.

[9] DIMAGGIO S J, SAKO B H. Basic system identification for condition monitoring of turbopumps[C]. 2001 IEEE Aerospace Conference, 73189-73200.

[10] QIAN S,CHEN D. Signal representation using adaptive normalized Gaussian functions[J]. Signal Processing, 1994, 36(1): 329-355.

[11] FAROOK S, GORAN S. Nonparametric waveform estimation using filter bank[J]. IEEE Transactions on Signal Processing, 1996,44(2): 240-247.

[12] 費(fèi)曉琪，孟慶豐，何正嘉. 基于沖擊時(shí)頻原子的匹配追蹤信號(hào)分解及機(jī)械故障特征提取技術(shù)[J].振動(dòng)與沖擊，2003, 22(2): 26-29.

FEI Xiaoqi, MENG Qingfeng, HE Zhengjia.Signal decomposition with matching pursuits and technology of extraction machinery fault feature based on impulse time-frequency atom[J]. Journal of Vibration and Shock, 2003, 22(2): 26-29.

[13] 王江萍. 機(jī)械設(shè)備故障診斷技術(shù)及應(yīng)用[M]. 西安:西北工業(yè)大學(xué)出版社，2004.

[14] 梁瑜，賈利民，蔡國(guó)強(qiáng)，等. 滾動(dòng)軸承的非線性動(dòng)力學(xué)故障模型研究[J]. 中國(guó)鐵道科學(xué)，2014,35(1):98-104.

LIANG Yu, JIA Limin, CAI Guoqiang, et al. Research on nonlinear dynamics fault model of rolling bearing[J]. China Railway Science, 2014, 35(1): 98-104.

RailwayvehiclerollingbearingfaultdiagnosismethodbasedonthematchingpursuitpretreatmentCHENNa1,3,YANG

Shaopu2,PANCunzhi2

(1.School of Computing and Informatics, Shijiazhuang Tiedao University, Shijiazhuang 050043, China； 2. School of Mechanical Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China； 3. School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China)

A bearing fault diagnosis algorithm was proposed based on the matching pursuit pretreatment. The method can avoid the dependence on the bandpass filtering in traditional envelope spectrum methods. The method also can remove the additional interference frequency components from the processing chain of traditional methods, and it can precisely extract the failure frequencies of damaged bearings. Through the experimental verification by using the real signals collected on a running-wheel bench, it is shown the method has some advantages over the traditional method of envelope spectrum, and it can effectively extract the failure features of wheel bearing damages.

matching pursuit; fault diagnosis; intermodulation; wheel bearing

TH165+.3, TN911.6

A

10.13465/j.cnki.jvs.2017.21.020

國(guó)家自然科學(xué)基金(11372197)；國(guó)家自然科學(xué)基金儀器專項(xiàng)(11227201)；教育部新世紀(jì)優(yōu)秀人才計(jì)劃(NCET-13-0913)

2016-02-23 修改稿收到日期：2016-09-02

陳娜 女，博士生，副教授，1979年生

楊紹普 男，教授，博士生導(dǎo)師，1962年生