李勇發(fā)，左小清，楊芳，徐晶

(1.昆明理工大學(xué) 國(guó)土資源工程學(xué)院，昆明 650093;2.黑龍江科技大學(xué) 理學(xué)院，哈爾濱 150027)

利用信號(hào)處理技術(shù)對(duì)滾動(dòng)軸承進(jìn)行故障檢測(cè)主要涉及2個(gè)問題：1)提取故障信號(hào)的特征值，實(shí)現(xiàn)降維和去噪的目的；2)利用訓(xùn)練得到的模型進(jìn)行故障檢測(cè)[1-2]。故障數(shù)據(jù)缺乏是故障檢測(cè)的一大難題，因難以獲取樣本數(shù)據(jù)而導(dǎo)致傳統(tǒng)的檢測(cè)方法效果不理想。

支持向量數(shù)據(jù)描述(Support Vector Data Des-cription,SVDD)具有訓(xùn)練速度快、普適性強(qiáng)等優(yōu)點(diǎn)，而且在分類過(guò)程中僅需一類樣本，能夠有效處理小樣本數(shù)據(jù)[2-3]。但SVDD算法需要對(duì)二次規(guī)劃問題進(jìn)行求解，求解過(guò)程中的計(jì)算量與數(shù)據(jù)維數(shù)成正比，數(shù)據(jù)維數(shù)越高，就需要越多的樣本數(shù)去準(zhǔn)確描述目標(biāo)類[3]。因此，嘗試?yán)眯〔ǚ治龇椒ǚ纸獬跏夹盘?hào)，計(jì)算小波奇異譜值并進(jìn)行優(yōu)選，實(shí)現(xiàn)信號(hào)的降維和去噪，然后再使用SVDD進(jìn)行故障診斷，以取得更好的分類效果。

1 故障診斷流程

1.1 小波奇異值優(yōu)選

在限定時(shí)間序列長(zhǎng)度為N=2J的條件下，由信號(hào)J個(gè)尺度的分解構(gòu)成的(J+1)×N矩陣為D(J+1)×N，其中：DJ表示信號(hào)在j(j=1,2,…,J)尺度下的小波分解細(xì)節(jié)，SJ表示信號(hào)J層小波分解的尺度。

利用小波變換得到下列矩陣

(1)

首先利用小波變換對(duì)信號(hào)進(jìn)行分解并計(jì)算出奇異譜值，然后以該奇異譜值作為特征進(jìn)行故障檢測(cè)[2]。為了提高SVDD算法的分類效果，需對(duì)奇異值進(jìn)行優(yōu)選。獲得的數(shù)據(jù)不同，奇異值優(yōu)選的方法也不同，如果擁有部分故障數(shù)據(jù)，應(yīng)當(dāng)選擇最適合故障檢測(cè)的奇異值；如果沒有故障數(shù)據(jù)，則首先假定一個(gè)閾值K，將序號(hào)小于這個(gè)閾值的奇異值作為特征值。

1.2 支持向量數(shù)據(jù)描述

SVDD是一種單值分類方法，其理論基礎(chǔ)源于支持向量機(jī)[6]，但在構(gòu)建最優(yōu)超平面上與有監(jiān)督的支持向量機(jī)不同。SVDD算法構(gòu)建一個(gè)最優(yōu)超球體將目標(biāo)類數(shù)據(jù)(球內(nèi))與非目標(biāo)類數(shù)據(jù)(球外)區(qū)分開，對(duì)于輸入空間不可分時(shí)，可引入核函數(shù)將其轉(zhuǎn)化為線性可分，然后再構(gòu)造最優(yōu)超球體，并使其半徑最小且包含大多數(shù)據(jù)點(diǎn)。該超球體就是區(qū)分目標(biāo)類與非目標(biāo)類的描述模型，如圖1所示。

圖1 超球體示意圖

對(duì)于給定的含有N個(gè)數(shù)據(jù)的訓(xùn)練樣本集{xi,i=1,2,…,N}，構(gòu)造中心為a和半徑為R的超球體，定義結(jié)構(gòu)風(fēng)險(xiǎn)ε=(R,a)=R2，限制條件為xi-a≤R2，引入松弛因子ξi和懲罰參數(shù)C，對(duì)初始函數(shù)進(jìn)行最小化，即

(2)

s.txi-a2≤R2+ξi,

ξi≥0(i=1,2,…,N)。

引入Lagrange乘子進(jìn)行化簡(jiǎn)后得

(3)

對(duì)Lagrange函數(shù)求關(guān)于R，a，ξi的偏導(dǎo)數(shù)，并令其等于0得

。(4)

求出當(dāng)L達(dá)到最小值時(shí)的αi。當(dāng)αi=0時(shí)對(duì)應(yīng)的樣本點(diǎn)xi為內(nèi)部點(diǎn)，即目標(biāo)樣本點(diǎn)；αi=C時(shí)對(duì)應(yīng)的xi為外部點(diǎn)，即非目標(biāo)點(diǎn);0≤αi

。(5)

對(duì)于某一待測(cè)樣本Z，可用決策函數(shù)劃分其所屬類別，即

fSVDD(z,a,R)=(z-a)(z-a)T=

(6)

(7)

式中：Q為指示函數(shù)。若Q為真，則標(biāo)記為-1，說(shuō)明是目標(biāo)類樣本；反之標(biāo)記為1，說(shuō)明是非目標(biāo)類樣本。

檢驗(yàn)一組新樣本數(shù)據(jù)的決策函數(shù)為

(8)

從上述分析可知，SVDD模型只需通過(guò)訓(xùn)練得到球體中心a和球體半徑R，即可構(gòu)成故障檢測(cè)模型，用于新樣本的識(shí)別。

核函數(shù)是將輸入空間通過(guò)非線性函數(shù)映射到高維特征空間，可以用原始輸入空間中核函數(shù)的計(jì)算替代高維特征空間內(nèi)積的計(jì)算，從而大大減少計(jì)算量[9-10]。選用高斯核函數(shù)代替內(nèi)積，即

(x.y)→K(x,y)=exp(-‖x-y‖2/s2)，

(9)

則檢驗(yàn)一個(gè)新數(shù)據(jù)對(duì)象的決策函數(shù)變?yōu)?/p>

1.3 故障檢測(cè)流程

首先，利用db5小波對(duì)正常振動(dòng)信號(hào)進(jìn)行5層小波分解，計(jì)算出奇異值進(jìn)行優(yōu)選；然后，將優(yōu)選出來(lái)的奇異值輸入SVDD中進(jìn)行訓(xùn)練，并利用(5)式計(jì)算超球體中心和半徑；最后，利用判別函數(shù)進(jìn)行分類(正?；蚬收?，故障檢測(cè)流程如圖2所示。

圖2 基于SVDD的故障檢測(cè)流程

2 試驗(yàn)分析

某型號(hào)軸承正常狀態(tài)及內(nèi)圈、外圈、滾動(dòng)體發(fā)生磨損的故障狀態(tài)下原始信號(hào)的時(shí)域波形如圖3所示[11]。從圖中可以看出，正常狀態(tài)下振動(dòng)信號(hào)比較平穩(wěn)，幅值較低；故障信號(hào)的幅值較大，其中滾動(dòng)體振動(dòng)信號(hào)中包含周期性的瞬時(shí)沖擊成分。

圖3 原始信號(hào)的時(shí)域波形圖

2.1 線性區(qū)分能力分析

每種狀態(tài)取20組樣本，利用db5小波進(jìn)行分解并計(jì)算小波奇異譜值，結(jié)果如圖4所示。從圖中可以看出，利用小波變換得到的奇異值按遞減順序排列，在前10個(gè)奇異值特征上，正常樣本與內(nèi)、外圈故障數(shù)據(jù)的線性區(qū)分能力較為明顯，但與滾動(dòng)體之間的線性區(qū)分能力不明顯。

圖4 不同狀態(tài)下振動(dòng)信號(hào)的奇異譜均值

利用SVDD對(duì)前20個(gè)奇異值進(jìn)行線性區(qū)分能力因子測(cè)試，正常狀態(tài)與不同故障狀態(tài)間的線性區(qū)分能力如圖5所示。結(jié)果表明：正常數(shù)據(jù)與內(nèi)、外圈故障數(shù)據(jù)間的線性區(qū)分能力隨著奇異值序號(hào)的增加而降低，與滾動(dòng)體故障間的線性區(qū)分能力則隨奇異值序號(hào)的增加而逐漸上升。

圖5 不同奇異值的線性區(qū)分能力

定義線性檢測(cè)能力為

(16)

式中：p為故障數(shù)據(jù)；k為正常樣本數(shù)據(jù)；ri為第i個(gè)奇異值。不同奇異值維度的線性區(qū)分能力如圖6所示。

圖6 不同奇異值維度的線性區(qū)分能力

從圖中可以看出，第3個(gè)奇異值處線性區(qū)分能力較低，第8和第9維奇異值處噪聲對(duì)線性區(qū)分能力影響最強(qiáng)。為使正常信號(hào)與滾動(dòng)體信號(hào)之間的區(qū)分能力提高，將[8，9，7，6，10，5，11，4，2]9個(gè)奇異值維度作為特征向量輸入SVDD進(jìn)行故障檢測(cè)。

2.2 基于SVDD的軸承故障檢測(cè)方法

不同窗口大小對(duì)檢測(cè)正確率的影響如圖7所示。結(jié)果表明：數(shù)據(jù)窗口越大，檢測(cè)正確率越高；反之，檢測(cè)正確率越低。但數(shù)據(jù)量的加大會(huì)加倍的增加計(jì)算量，并造成系統(tǒng)延遲，因此應(yīng)根據(jù)實(shí)際情況選擇窗口大小。從圖中可以發(fā)現(xiàn)，當(dāng)窗口大小為1 024時(shí)，就已經(jīng)能夠達(dá)到較好的效果。

圖7 窗口大小對(duì)檢測(cè)性能的影響

高斯核參數(shù)對(duì)SVDD檢測(cè)正確率的影響如圖8所示。從圖中可以看出：高斯核寬度與正常樣本的檢測(cè)正確率成正比，與故障樣本的檢測(cè)正確率成反比。當(dāng)高斯核寬度為0.3時(shí)，正常及故障樣本的檢測(cè)正確率均達(dá)到最大。

圖8 不同核函數(shù)參數(shù)對(duì)檢測(cè)率的影響

利用優(yōu)選得到的9個(gè)奇異值作為特征向量輸入SVDD分類器進(jìn)行訓(xùn)練。選擇高斯核函數(shù)寬度為0.3，構(gòu)建好SVDD模型后，按 (7) 式求出超球體中心和半徑，即可對(duì)新樣本進(jìn)行分類。經(jīng)過(guò)訓(xùn)練得R=3.091 912，a=0.418 414。利用測(cè)試樣本(測(cè)試樣本數(shù)為2 936，其中內(nèi)圈、外圈、滾動(dòng)體故障樣本數(shù)分別為582，586，588)進(jìn)行模型檢驗(yàn)，分類效果如圖9所示，檢測(cè)正確率達(dá)到了98.6%。

3 結(jié)束語(yǔ)

采用小波分析提取奇異值并進(jìn)行優(yōu)選能夠有效提取可充分反映故障本質(zhì)的敏感特征，從而使訓(xùn)練和分類效果得到很大的提高。而通過(guò)選擇合理的窗口大小及高斯核寬度，SVDD算法的檢測(cè)正確率達(dá)到了98.6%，完全能夠滿足軸承故障診斷的要求。但上述方法主要檢測(cè)了有無(wú)故障發(fā)生，還不能準(zhǔn)確判斷故障發(fā)生部位，需做進(jìn)一步研究。

圖9 SVDD分類效果圖