胡愛軍，南 冰

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

基于自適應(yīng)概率主成分分析的滾動(dòng)軸承故障特征增強(qiáng)方法

胡愛軍，南 冰

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

針對(duì)實(shí)際工程中滾動(dòng)軸承微弱故障信號(hào)特征難以提取的問(wèn)題，提出了一種新的自適應(yīng)概率主成分分析(Adaptive Probabilistic Principal Component Analysis, APPCA)的軸承故障特征增強(qiáng)方法。概率主成分分析(Probabilistic Principal Component Analysis, PPCA) 能夠提取信號(hào)的主要故障特征，去除背景噪聲干擾，但結(jié)果易受到主成分?jǐn)?shù)與原始變量維數(shù)選擇的影響。為了自適應(yīng)實(shí)現(xiàn)最佳分析結(jié)果，利用粒子群算法多參數(shù)尋優(yōu)特性，根據(jù)最大峭度準(zhǔn)則確定影響PPCA的最佳影響參數(shù)組合。原信號(hào)通過(guò)APPCA方法處理后，背景噪聲得到有效抑制，故障特征得到增強(qiáng)，最后通過(guò)包絡(luò)分析識(shí)別故障特征。仿真和實(shí)驗(yàn)結(jié)果證明了該方法的有效性。

滾動(dòng)軸承；概率主成分分析；故障診斷

滾動(dòng)軸承出現(xiàn)早期局部故障時(shí)，由于故障產(chǎn)生的沖擊成分通常比較微弱，并淹沒(méi)于強(qiáng)烈的背景噪聲中，導(dǎo)致難以提取到軸承故障特征頻率信息[1]。近些年，針對(duì)此類問(wèn)題，不少學(xué)者進(jìn)行了深入研究。唐貴基等[2]提出先對(duì)故障信號(hào)做最大相關(guān)峭度解卷積預(yù)處理增強(qiáng)信號(hào)沖擊特征，然后計(jì)算解卷積信號(hào)的包絡(luò)信號(hào)，最后對(duì)包絡(luò)信號(hào)做1.5維譜分析，通過(guò)分析譜圖中幅值突出的頻率成分準(zhǔn)確判斷了故障類型。文獻(xiàn)[3]首先利用小波包去除信號(hào)中的噪聲，后進(jìn)行LMD分解，選取有效PF集進(jìn)行功率譜分析，成功提取出了故障特征。文獻(xiàn)[4]提出了基于倒譜編輯預(yù)白化和形態(tài)學(xué)自互補(bǔ)Top-Hat變化的方法，先對(duì)信號(hào)預(yù)白化，后對(duì)白化信號(hào)形態(tài)學(xué)濾波消除背景噪聲的干擾，準(zhǔn)確的提取了軸承故障特征頻率。文獻(xiàn)[5]將EEMD與1.5維能量譜結(jié)合，對(duì)軸承內(nèi)圈故障信號(hào)分析取得了較為理想的效果。文獻(xiàn)[6]將EEMD、度量因子和快速峭度圖相結(jié)合，較好的提取了軸承故障特征頻率。上述方法在軸承故障診斷中均取得了一定的效果，但是這些方法均存在一定的局限性，由于缺少對(duì)原始數(shù)據(jù)的恰當(dāng)?shù)母怕誓Ｐ?，?dǎo)致一些高頻的噪聲不能夠被正確的分離，同時(shí)某些故障特征信息可能被當(dāng)作噪聲去除[7-8]。

概率主成分分析(Probabilistic Principal Component Analysis, PPCA)，是一種信號(hào)分析方法，其首先建立一個(gè)恰當(dāng)?shù)母怕誓Ｐ?，然后基于這個(gè)模型重新生成一個(gè)新的樣本數(shù)據(jù)，最后信號(hào)主成分可以通過(guò)正交投影的方法獲得。PPCA的本質(zhì)是將方差最大的方向作為主要特征，并且在各個(gè)正交方向上將數(shù)據(jù)“離相關(guān)”，也就是讓它們?cè)诓煌环较蛏蠜](méi)有相關(guān)性。因此PPCA不僅可以去除噪聲，還能增強(qiáng)對(duì)原始信號(hào)特征信息的保留，現(xiàn)已應(yīng)用于特征提取與模態(tài)識(shí)別等領(lǐng)域[9]。文獻(xiàn)[10]通過(guò)多次實(shí)驗(yàn)對(duì)比分析的方法，針對(duì)特定的故障信號(hào)選擇最優(yōu)的主成分?jǐn)?shù)k與原始變量維數(shù)n，較好的提取出了軸承故障頻率的邊頻帶。但由于其算法不具有自適應(yīng)性，針對(duì)不同的故障信號(hào)時(shí)，最優(yōu)參數(shù)取值需要重新進(jìn)行對(duì)比分析，因此在實(shí)際軸承故障診斷中受到了一定程度的限制。

鑒于上述問(wèn)題，本文提出了一種自適應(yīng)的概率主成分分析方法。為了自適應(yīng)實(shí)現(xiàn)最佳分析結(jié)果，利用粒子群多參數(shù)尋優(yōu)特性，根據(jù)最大峭度準(zhǔn)則確定影響PPCA的最佳影響參數(shù)組合，并應(yīng)用于滾動(dòng)軸承的故障特征增強(qiáng)。通過(guò)仿真和實(shí)驗(yàn)結(jié)果，證明了該方法能有效的增強(qiáng)軸承故障特征，適合用于軸承故障診斷。

1 PPCA基本原理和方法

PPCA作為一種信號(hào)分析方法，通過(guò)先將原始數(shù)據(jù)投影到其他的坐標(biāo)空間，后投影的方式來(lái)提取信號(hào)的主特征分量。其本質(zhì)是將方差最大的方向作為主要特征，并且在各個(gè)正交方向上將數(shù)據(jù)“離相關(guān)”，因此PPCA不僅可以去除噪聲，還能增強(qiáng)對(duì)原始信號(hào)特征信息的保留。

1.1PPCA模型

PPCA模型首先假設(shè)n維原始變量數(shù)據(jù)X滿足如下模型關(guān)系[11-12]

X=P·u+E

(1)

X～(0,PPT+σ2I)

(2)

其中u的先驗(yàn)分布為

(3)

且原始數(shù)據(jù)X在隱變量u條件下的先驗(yàn)概率分布為

(4)

根據(jù)式(3)和式(4)可得原始數(shù)據(jù)X的概率分布為

(5)

式中，C=PPT+σ2I為兩參數(shù)P與σ2決定的協(xié)方差矩陣。為了得到上述模型的P和σ2采用EM算法進(jìn)行估計(jì)，推導(dǎo)出其迭代公式

P=SP(σ2I+M-1PTSP)-1

(6)

(7)

式中：S為原始數(shù)據(jù)的協(xié)方差矩陣;M=PTP+σ2I；兩參數(shù)P與σ2可由式(6)和式(7)多次迭代至收斂求得，P與σ2計(jì)算出來(lái)后即可建立PPCA模型。

1.2PPCA降維

當(dāng)PPCA模型被建立后，可以用以下變換求取降維后的數(shù)據(jù)。

(8)

由式(8)可知，各主成分?jǐn)?shù)據(jù)(降維后的數(shù)據(jù))是原始變量數(shù)據(jù)X在相應(yīng)主成分向量pi的投影。PPCA去噪效果由主成分?jǐn)?shù)k與原始變量維數(shù)n確定，當(dāng)其中任何一個(gè)參數(shù)設(shè)置不合適，都難以達(dá)到理想的分析效果?，F(xiàn)階段，參數(shù)主要是根據(jù)人為經(jīng)驗(yàn)選取，缺乏自適應(yīng)性，因此在實(shí)際軸承故障診斷中受到了一定程度的限制。

2 APPCA方法原理與實(shí)現(xiàn)過(guò)程

粒子群算法[13]具有良好的全局尋優(yōu)能力，本文采用粒子群算法對(duì)PPCA算法的兩個(gè)參數(shù)進(jìn)行優(yōu)化，可實(shí)現(xiàn)主成分?jǐn)?shù)k與原始變量維數(shù)n的自適應(yīng)選取。假設(shè)D維空間中，種群X=(X1,X1,…,XM)包含M個(gè)粒子組成，其中第i個(gè)粒子表示一個(gè)D維向量Xi=(xi1,xi2,…,xiD)，代表第i個(gè)粒子在D維搜索空間中的位置。第i個(gè)粒子的速度為Vi=(vi1,vi2,…,viD)，Qi=(qi1,qi2,…,qiD)為個(gè)體局部均值，G=(g1,g2,…,gD)為種群全局極值，各粒子通過(guò)Pi和G迭代更新自身速度和位置，公式為

(9)

式中:ω為慣性權(quán)重；d=1,2,…,D；i=1,2,…,M；α為當(dāng)前迭代次數(shù)；c1和c2為加速度因子；η為介于[0,1]的隨機(jī)數(shù)。

粒子群算法尋優(yōu)時(shí)，需要確定一個(gè)適應(yīng)度函數(shù)，粒子每次更新位置時(shí)需計(jì)算當(dāng)前位置對(duì)應(yīng)的函數(shù)值，通過(guò)對(duì)比進(jìn)行更新。信號(hào)中沖擊成分的比重影響著峭度指標(biāo)的大小，信號(hào)所包含的沖擊成分越多，峭度值越大，因此本文將PPCA處理后信號(hào)的峭度值K作為適應(yīng)度函數(shù)，表達(dá)式為

(10)

式中:x為振動(dòng)信號(hào);μ為信號(hào)x的均值;δ為信號(hào)x的標(biāo)準(zhǔn)差。

利用粒子群參數(shù)尋優(yōu)的具體步驟如下：

步驟1初始化粒子群算法的各項(xiàng)參數(shù)。

步驟2初始化粒子種群，以影響參數(shù)組合[k,n]作為粒子的位置，隨機(jī)初始化各粒子的位置與移動(dòng)速度。

步驟3計(jì)算每個(gè)粒子的位置對(duì)應(yīng)的適應(yīng)度值。

步驟4對(duì)比適應(yīng)度值大小并更新個(gè)體局部極值和種群全局極值。

步驟5更新粒子的位置和速度。

步驟6循環(huán)迭代，轉(zhuǎn)置步驟3,直至迭代次數(shù)達(dá)到最大設(shè)定值，輸出最佳適應(yīng)度值及粒子的位置。

在設(shè)置粒子群尋優(yōu)參數(shù)時(shí)，參考了文獻(xiàn)[14-15]中的參數(shù)取值，如表1所示。

表1 粒子群算法各項(xiàng)參數(shù)Tab.1 Each parameter of particle swarm algorithm

自適應(yīng)概率主成分分析故障特征增強(qiáng)方法實(shí)現(xiàn)過(guò)程如圖1所示。利用粒子群多參數(shù)尋優(yōu)特性，根據(jù)最大峭度準(zhǔn)則確定影響PPCA的最佳影響參數(shù)組合，可有效避免參數(shù)設(shè)定時(shí)人為主觀因素帶來(lái)的弊端。原信號(hào)通過(guò)APPCA方法處理后，背景噪聲得到有效抑制，故障特征得到增強(qiáng)，最后通過(guò)包絡(luò)分析識(shí)別故障特征。

圖1 APPCA方法流程圖Fig.1 Flow chart of APPCA method

3 仿真分析及應(yīng)用

3.1仿真信號(hào)分析

采用文獻(xiàn)[16]中的滾動(dòng)軸承內(nèi)圈故障模型進(jìn)行模擬。故障數(shù)學(xué)模型如式(11)所示

(11)

式中，τi為第i次沖擊相對(duì)于平均周期T的微小波動(dòng)；Ai為以1/fr為周期的幅值調(diào)制；h(t)為指數(shù)衰減脈沖；B為系統(tǒng)的衰減系數(shù)；A0=2，CA=0；fr=20 Hz為軸承所在工作軸的轉(zhuǎn)頻；fi=150 Hz為內(nèi)圈故障通過(guò)頻率；fn=3 kHz為系統(tǒng)固有頻率；n(t)為信噪比-12 db的高斯白噪聲。設(shè)置采樣頻率為fs=12 800 Hz，取4 096點(diǎn)數(shù)據(jù)分析。加噪故障仿真信號(hào)的時(shí)域波形如圖2所示，圖3為直接對(duì)加噪軸承內(nèi)圈故障仿真信號(hào)做包絡(luò)譜的分析結(jié)果。圖3包絡(luò)譜中沒(méi)有找到幅值突出的頻率成分，說(shuō)明僅包絡(luò)分析難以提取到強(qiáng)背景噪聲下的軸承微弱故障信號(hào)。

圖2 加噪內(nèi)圈故障仿真信號(hào)時(shí)域波形Fig.2 Time domain waveform of inner ring fault simulation signal with noise

圖3 內(nèi)圈故障仿真信號(hào)包絡(luò)譜Fig.3 The envelope spectrum of the inner ring fault simulation signal

利用APPCA方法對(duì)故障信號(hào)進(jìn)行分析，首先將所選4 096點(diǎn)軸承故障數(shù)據(jù)xr去均值得到一維數(shù)據(jù)x，將x構(gòu)造n維原始變量數(shù)據(jù)X如式(12)所示

(12)

其次將參數(shù)P與σ2值初始化后按照式(6)式(7)經(jīng)過(guò)多次迭代求解參數(shù)值，當(dāng)兩參數(shù)取值確定后，按照式(1)建立PPCA模型。在模型中n與k的取值直接影響到主成分提取的效果，根據(jù)參考文獻(xiàn)[10]的取值，當(dāng)固定參數(shù)k的取值為2，參數(shù)n的取值大于20時(shí)，軸承故障信號(hào)的信噪比會(huì)降低，不利于軸承故障特征的增強(qiáng)，因此本文參數(shù)n的最大取值設(shè)置為20。另外由于PPCA算法的本質(zhì)是優(yōu)先將方差最大的方向作為主要特征，軸承故障信息主要集中在特征值較大的主成分中，背景噪聲主要分布在特征值較小的主成分中，k的取值較大時(shí)，特征值較小的噪聲成分會(huì)增加，同時(shí)計(jì)算負(fù)擔(dān)也會(huì)加重，所以參數(shù)k的最大值不宜過(guò)大。采用粒子群算法對(duì)主成分?jǐn)?shù)k與原始變量維數(shù)n進(jìn)行自適應(yīng)選取，圖4為峭度值隨進(jìn)化代數(shù)變化的關(guān)系曲線，峭度最大值3.65出現(xiàn)在了第11代進(jìn)化種群中，此時(shí)主成分?jǐn)?shù)k=2，原始變量維數(shù)n=13。

圖4 峭度值隨進(jìn)化代數(shù)變化的關(guān)系曲線Fig.4 Relation curve of kurtosis changing with evolutional generation

PPCA模型建立后，按照式(8)求解出主成分矩陣 (即降維后的數(shù)據(jù))，最后將主成分矩陣重構(gòu)軸承故障仿真信號(hào)時(shí)域波形如圖5所示。圖6為軸承故障仿真信號(hào)的包絡(luò)譜。

圖5 PPCA處理后仿真信號(hào)時(shí)域波形Fig.5 Time domain waveform of the simulation signal using PPCA

圖6 APPCA方法所得信號(hào)包絡(luò)譜Fig.6 The envelope spectrum of the signal by APPCA method

圖6包絡(luò)譜中可以較為清楚的找到150 Hz、300 Hz頻率成分，分別對(duì)應(yīng)軸承故障特征頻率及其二倍頻，故障特征頻率三、四倍頻譜線峰值也十分突出。說(shuō)明APPCA算法處理信號(hào)能夠去除大量背景噪聲，增強(qiáng)軸承的故障特征，效果較理想。

為驗(yàn)證所述方法的優(yōu)勢(shì)，與故障診斷領(lǐng)域常用的快速峭度圖方法結(jié)果作對(duì)比。圖7為故障信號(hào)快速峭度圖分析結(jié)果，可以看出軸承故障信號(hào)最大譜峭度處所對(duì)應(yīng)的分解層數(shù)為7，濾波器中心頻率和帶寬組合為[6 300,200]，即帶通濾波器的范圍為[6 200,6 400]，在此范圍內(nèi)信號(hào)的峭度值與信噪比均達(dá)到最大，按照上述的中心頻率和帶寬構(gòu)造濾波器對(duì)信號(hào)進(jìn)行濾波，所得包絡(luò)譜如圖8所示。通過(guò)對(duì)比分析發(fā)現(xiàn)：快速峭度圖濾波后包絡(luò)譜僅能勉強(qiáng)提取到軸承故障特征頻率(150 Hz)，倍頻成分均被背景噪聲淹沒(méi)且由于存在較多的干擾譜線(28.1 Hz、65.6 Hz)，無(wú)法判斷出是內(nèi)圈故障，與圖6相比，故障特征頻率幅值水平也存在較大差距。由此表明，APPCA方法在軸承故障特征增強(qiáng)方面具有一定的優(yōu)勢(shì)。

圖7 內(nèi)圈故障仿真信號(hào)快速峭度圖Fig.7 Fast kurtogram of the inner ring fault simulation signal

圖8 快速峭度圖濾波后信號(hào)包絡(luò)譜Fig.8 The envelope spectrum of the signal using fast kurtogram filter

3.2實(shí)驗(yàn)信號(hào)分析

為了進(jìn)一步驗(yàn)證該方法對(duì)實(shí)際軸承故障信號(hào)的處理效果，采用美國(guó)Case Western Reserve大學(xué)的滾動(dòng)軸承實(shí)驗(yàn)數(shù)據(jù)，軸承型號(hào)JEMSKF6023-2RS。故障源是滾動(dòng)體表面通過(guò)電火花加工的直徑分別為0.177 8 mm(0.007inch)、0.355 6 mm(0.014inch)、0.533 4 mm(0.021inch)的凹坑。采樣頻率12 kHz，軸的轉(zhuǎn)速為1 772 r/min。表2為軸承的各個(gè)故障特征頻率。

選用最輕微的滾動(dòng)體0.007inch故障數(shù)據(jù)進(jìn)行分析，分析點(diǎn)數(shù)取8 192點(diǎn)，圖9為軸承故障信號(hào)的時(shí)域波形。直接對(duì)軸承信號(hào)做包絡(luò)分析，結(jié)果如圖10所示。包絡(luò)譜中存在較多的干擾譜線，僅能提取到接近軸承轉(zhuǎn)頻的頻率成分(30.0 Hz)，115.7 Hz與滾動(dòng)體故障頻率118 Hz相差較大，且無(wú)法找到故障頻率的倍頻成分。因此對(duì)于此輕微故障，僅包絡(luò)分析效果欠佳。

表2 滾動(dòng)軸承故障特征頻率Tab.2 Fault feature frequency of rolling bearing

圖9 軸承故障信號(hào)時(shí)域波形Fig.9 Time domain waveform of the bearing fault signal

圖10 軸承故障信號(hào)包絡(luò)譜Fig.10 The envelope spectrum of the bearing fault signal

利用APPCA方法對(duì)故障信號(hào)進(jìn)行分析。圖11為峭度值隨進(jìn)化代數(shù)變化的關(guān)系曲線，峭度最大值3.812出現(xiàn)在了第10代進(jìn)化種群中，此時(shí)PPCA算法中主成分?jǐn)?shù)k=2，原始變量維數(shù)n=20。利用參數(shù)優(yōu)化后的PPCA算法對(duì)故障信號(hào)處理所得包絡(luò)譜如圖12所示。

圖11 峭度值隨進(jìn)化代數(shù)變化的關(guān)系曲線Fig.11 Relation curve of kurtosis changing with evolutional generation

圖12包絡(luò)譜中可以較為清楚的找到29.3 Hz、58.6 Hz、117.2 Hz、234.4 Hz等頻率成分。其中29.3 Hz、58.6 Hz分別對(duì)應(yīng)軸承轉(zhuǎn)頻及其倍頻。117.2 Hz、234.4 Hz與滾動(dòng)體故障特征頻率、二倍頻成分非常接近，因此可以判斷實(shí)際情況是軸承滾動(dòng)體存在故障。

圖12 本文方法所得信號(hào)包絡(luò)譜Fig.12 The envelope spectrum of the signal by the proposed method

同樣，將所述方法與快速峭度圖方法結(jié)果作對(duì)比，圖13為故障信號(hào)快速峭度圖分析結(jié)果?？梢钥闯鲚S承故障信號(hào)最大譜峭度處所對(duì)應(yīng)的分解層數(shù)為7，濾波器中心頻率和帶寬組合為[2 343.75,187.5]，即帶通濾波器的范圍為[2 250,2 437.5]，在此范圍內(nèi)信號(hào)的峭度值與信噪比均達(dá)到最大，按照上述的中心頻率和帶寬構(gòu)造濾波器對(duì)信號(hào)進(jìn)行濾波，所得包絡(luò)譜如圖14所示。

圖13 軸承故障信號(hào)快速峭度圖Fig.13 Fast kurtogram of the bearing fault signal

圖14 快速峭度圖濾波后信號(hào)包絡(luò)譜Fig.14 The envelope spectrum of the signal using fast kurtogram filter

由圖14可以看出包絡(luò)譜中存在較多的干擾譜線(8.8 Hz、90.8 Hz、114.3 Hz)，而軸承轉(zhuǎn)頻及其倍頻、滾動(dòng)體故障特征頻率及其倍頻被淹沒(méi)在了背景噪聲中，無(wú)法確定軸承存在滾動(dòng)體故障。而APPCA處理后的包絡(luò)譜圖相對(duì)干凈，基本不存在其他干擾成分，轉(zhuǎn)頻及其倍頻處譜線峰值十分突出，同時(shí)也較易提取到軸承的輕微故障特征。滾動(dòng)體輕微故障特征信號(hào)對(duì)比分析結(jié)果再次驗(yàn)證了本文所述方法在軸承故障特征增強(qiáng)方面的優(yōu)勢(shì)。

4 結(jié) 論

實(shí)際工程中滾動(dòng)軸承微弱故障信號(hào)易受到背景噪聲的干擾，導(dǎo)致僅包絡(luò)分析難以提取出軸承故障特征。PPCA能夠提取信號(hào)主要故障特征成分，去除背景噪聲干擾，但其算法中主成分?jǐn)?shù)k與原始變量維數(shù)n選擇起著十分關(guān)鍵的作用。利用粒子群多參數(shù)尋優(yōu)特性，根據(jù)最大峭度準(zhǔn)則確定影響PPCA的最佳影響參數(shù)組合的方法是有效的，可以用于軸承故障特征增強(qiáng)。仿真和實(shí)驗(yàn)對(duì)比結(jié)果證明了該方法的有效性。

[1] 胡愛軍, 馬萬(wàn)里, 唐貴基. 基于集成經(jīng)驗(yàn)?zāi)B(tài)分解和峭度準(zhǔn)則的滾動(dòng)軸承故障特征提取方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2012, 32(11):106-111.

HU Aijun, MA Wanli, TANG Guiji. Rolling bearing fault feature extraction method based on ensemble empirical mode decomposition and kurtosis criterion[J]. Proceedings of the CSEE, 2012, 32(11):106-111.

[2] 唐貴基, 王曉龍. 最大相關(guān)峭度解卷積結(jié)合1.5維譜的滾動(dòng)軸承早期故障特征提取方法[J]. 振動(dòng)與沖擊, 2015,34(12):79-84.

TANG Guiji, WANG Xiaolong. Feature extraction for rolling bearing incipient fault based on maximum correlated kurtosis deconvolution and 1.5 dimension spectrum[J]. Journal of Vibration and Shock, 2015,34(12):79-84.

[3] 孫偉, 熊邦書, 黃建萍,等. 小波包降噪與LMD相結(jié)合的滾動(dòng)軸承故障診斷方法[J]. 振動(dòng)與沖擊, 2012, 31(18):153-156.

SUN Wei, XIONG Bangshu, HUANG Jianping, et al. Fault diagnosis of a rolling bearing using Wavelet packet de-noising and LMD[J]. Journal of Vibration and Shock, 2012, 31(18):153-156.

[4] 鄧飛躍, 唐貴基, 何玉靈. 基于倒譜預(yù)白化和形態(tài)學(xué)自互補(bǔ)Top-Hat 變換的滾動(dòng)軸承故障特征提取[J]. 振動(dòng)與沖擊, 2015,34(15):77-81.

DENG Feiyue, TANG Guiji, HE Yuling. Fault feature extraction for rolling element bearings based on cepstrum pre-whitening and morphology self-complementary top-hat transformation[J]. Journal of Vibration and Shock, 2015,34(15):77-81.

[5] 唐貴基, 王曉龍. 基于EEMD降噪和1.5維能量譜的滾動(dòng)軸承故障診斷研究[J]. 振動(dòng)與沖擊, 2014,33(1):6-10.

TANG Guiji, WANG Xiaolong. Fault diagnosis for roller bearings based on EEMD de-noising and 1.5 dimensional energy spectrum[J]. Journal of Vibration and Shock, 2014,33(1):6-10.

[6] 彭暢, 柏林, 謝小亮. 基于EEMD、度量因子和快速峭度圖的滾動(dòng)軸承故障診斷方法[J]. 振動(dòng)與沖擊, 2012,31(20):143-146.

PENG Chang, BO Lin, XIE Xiaoliang. Fault diagnosis method of rolling element bearings based on EEMD, measure-factor and fast kurtogram[J]. Journal of Vibration and Shock, 2012,31(20):143-146.

[7] TIPPING M E, BISHOP C M. Probabilistic principal component analysis[J]. Journal of the Royal Statistical Society, 1999, 61(3):611-622.

[8] BISHOP C M, TIPPING M E. A hierarchical latent variable model for data visualization[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence, 1998, 20(3):281-293.

[9] 陸超, 陳捷, 洪榮晶. 采用概率主成分分析的回轉(zhuǎn)支承壽命狀態(tài)識(shí)別[J]. 西安交通大學(xué)學(xué)報(bào), 2015, 49(10):90-96.

LU Chao, CHEN Jie, HONG Rongjing. Recognition of life state for slewing bearings using probabilistic principal component analysis[J]. Journal of Xi’an Jiaotong University, 2015, 49(10):90-96.

[10] XIANG J, ZHONG Y, GAO H. Rolling element bearing fault detection using PPCA and spectral kurtosis[J]. Measurement, 2015, 75:180-191.

[11] BELLAS A, BOUVEYRON C, COTTRELL M, et al. Model-based clustering of high-dimensional data streams with online mixture of probabilistic PCA[J]. Advances in Data Analysis & Classification, 2013, 7(3):281-300.

[12] ZUCCOLOTTO P. Principal component analysis with interval imputed missing values[J]. Asta Advances in Statistical Analysis, 2012, 96(1):1-23.

[13] KENNEDY J, EBERHART R. Particle swarm optimization[C]∥Proceeding of IEEE International Conference on Neural Networks. Perth: IEEE, 1995:1942-1948.

[14] 沈伋, 韓麗川, 沈益斌. 基于粒子群算法的飛機(jī)總體參數(shù)優(yōu)化[J]. 航空學(xué)報(bào), 2008, 29(6):1538-1541.

SHEN Ji, HAN Lichuan, SHEN Yibin. Optimization of airplane primary parameters based on particle swarm algorithm[J]. Acta Aeronautica ET Astronautica Sinica, 2008, 29(6):1538-1541.

[15] 唐貴基, 王曉龍. 最大相關(guān)峭度解卷積結(jié)合稀疏編碼收縮的齒輪微弱故障特征提取[J]. 振動(dòng)工程學(xué)報(bào), 2015, 28(3):478-486.

TANG Guiji, WANG Xiaolong. Weak feature extraction of gear fault based on maximum correlated kurtosis deconvolution and sparse code shrinkage[J]. Journal of Vibration Engineering, 2015, 28(3):478-486.

[16] 王宏超, 陳進(jìn), 董廣明. 基于最小熵解卷積與稀疏分解的滾動(dòng)軸承微弱故障特征提取[J]. 機(jī)械工程學(xué)報(bào), 2013, 49(1):88-94.

WANG Hongchao, CHEN Jin, DONG Guangming. Fault diagnosis method for rolling bearing’s weak fault based on minimum entropy deconvolution and sparse decomposition[J]. Journal of Mechanical Engineering, 2013, 49(1):88-94.

Faultfeatureenhancementmethodforrollingbearingbasedonadaptiveprobabilisticprincipalcomponentanalysis

HU Aijun, NAN Bing

(School of Energy, Power and Mechanical Engineering, North China Electric Power University, Baoding 071003, China)

Aiming at the difficulty in extracting the features of weak fault signals of rolling element bearings in practical engineering, a new method named adaptive probabilistic principal component analysis (APPCA) was proposed to enhance the features of bearing faults. The method of PPCA is able to extract main fault features and remove background noise interferences, but is easily affected by the number of principal components and the dimension of original variables. In order to adaptively achieve the best analysis result, the particle swarm optimization algorithm with multi-parameter optimization characteristics was applied to search for the optimal combination of influencing parameters of PPCA based on the maximum kurtosis criterion. After the original signal was processed by the APPCA method, the background noise was effectively suppressed, and the fault features were enhanced. Finally, the signal envelope spectrum was analyzed to identify fault features. The simulation and experiment results show the effectiveness of the method.

rolling bearing; probabilistic principal component analysis; fault diagnosis

TH133.3；TH17

A

10.13465/j.cnki.jvs.2017.19.022

國(guó)家自然科學(xué)基金資助項(xiàng)目(51475164)

2016-06-23 修改稿收到日期：2016-08-20

胡愛軍 男，博士，副教授，1971年生