呂明珠， 劉世勛， 蘇曉明， 陳長(zhǎng)征

(1.沈陽(yáng)工業(yè)大學(xué) 機(jī)械工程學(xué)院，沈陽(yáng) 110870；2.遼寧裝備制造職業(yè)技術(shù)學(xué)院 自控學(xué)院，沈陽(yáng) 110161；3.中認(rèn)(沈陽(yáng))北方實(shí)驗(yàn)室有限公司，沈陽(yáng) 110164)

滾動(dòng)軸承是旋轉(zhuǎn)機(jī)械的主要部件之一，一旦發(fā)生故障會(huì)導(dǎo)致整機(jī)運(yùn)行失效，因此對(duì)軸承進(jìn)行狀態(tài)監(jiān)測(cè)和故障預(yù)測(cè)具有重要的工程意義[1]。狀態(tài)監(jiān)測(cè)的主要目的是判斷軸承的健康狀況，以便于及時(shí)處理已經(jīng)發(fā)生或即將發(fā)生的故障，保證運(yùn)行的可靠性和降低維修成本。

振動(dòng)信號(hào)分析在滾動(dòng)軸承狀態(tài)監(jiān)測(cè)中應(yīng)用最為普遍[2]。當(dāng)發(fā)生局部故障時(shí)，振動(dòng)信號(hào)會(huì)產(chǎn)生周期性的脈沖信號(hào)。然而，在早期故障階段，脈沖信號(hào)非常微弱且極易受到噪聲干擾，難以確定退化的起始點(diǎn)[3]。因此，振動(dòng)信號(hào)處理技術(shù)和故障特征提取技術(shù)成為了早期故障診斷領(lǐng)域的研究熱點(diǎn)。

變分模態(tài)分解(variational mode decompositon，VMD)是一種于2014年提出的信號(hào)分析方法[4]，它能將信號(hào)分解為若干個(gè)具有一定帶寬和中心頻率的模態(tài)分量，而且能夠克服經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decompositon，EMD)的模態(tài)混疊及集成經(jīng)驗(yàn)?zāi)B(tài)分解(ensemble empirical mode decompositon，EEMD)的時(shí)頻聚集性不高等問題。因此，近年來(lái)VMD方法在機(jī)械設(shè)備振動(dòng)信號(hào)處理方面得到了廣泛的應(yīng)用[5-7]。然而，在使用VMD進(jìn)行信號(hào)分解之前，需要預(yù)先設(shè)置分解層數(shù)K和二次懲罰因子α等參數(shù)。Zheng等[8]采用觀察中心頻率法確定VMD的分解層數(shù)，這種方法需要有一定的經(jīng)驗(yàn)且很難快速準(zhǔn)確地確定參數(shù)取值。任學(xué)平等[9]提出了用能量差曲線確定分解層數(shù)的方法，但需要對(duì)所有分解結(jié)果進(jìn)行比較，計(jì)算量大且沒有考慮到懲罰因子對(duì)分解結(jié)果的影響。近年來(lái)一些傳統(tǒng)的優(yōu)化算法如粒子群、遺傳算法、蟻群算法等開始逐漸應(yīng)用于VMD的參數(shù)優(yōu)化中，但普遍存在算法復(fù)雜、搜索速度慢、易于陷入局部最優(yōu)等問題?；依莾?yōu)化算法(grey wolf optimization,GWO)是一種新型智能算法，具有較強(qiáng)的搜索能力和回避局部最優(yōu)的特點(diǎn)[10]。因此，本文將GWO算法引入到VMD的參數(shù)選擇中，從而自適應(yīng)地獲得最優(yōu)參數(shù)組合，稱為自適應(yīng)變分模態(tài)分解(adaptive variational mode decompositon，AVMD)。在早期故障特征提取方面，傳統(tǒng)的峭度指標(biāo)對(duì)非周期干擾敏感，但隨著沖擊的密集這種敏感會(huì)逐漸失效。Chegini等[11]提出了一種新的特征提取方法，叫做包絡(luò)諧噪比(envelope harmonic-to-noise ratio,EHNR)，雖然能夠發(fā)現(xiàn)軸承早期故障異常點(diǎn)，但由于沒有很好地濾除干擾信號(hào)，造成在軸承運(yùn)行初期出現(xiàn)誤報(bào)警(虛警)問題。

針對(duì)以上問題，本文提出了一種基于AVMD和EHNR的滾動(dòng)軸承早期退化檢測(cè)方法。該方法包括信號(hào)處理和特征提取兩個(gè)步驟。在信號(hào)處理部分，先引入GWO算法優(yōu)化VMD的參數(shù)選擇，再?gòu)腣MD分解后的若干子模態(tài)中分理處有效模態(tài)分量并重構(gòu)信號(hào)；在特征提取部分，利用EHNR劃分軸承全壽命運(yùn)行的不同階段，通過設(shè)置報(bào)警閾值判斷是否出現(xiàn)異常點(diǎn)，實(shí)現(xiàn)了滾動(dòng)軸承早期退化點(diǎn)的檢測(cè)且沒有出現(xiàn)虛警問題。實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了本文所提方法的有效性。

1 理論基礎(chǔ)

1.1 變分模態(tài)分解(VMD)

VMD是一種自適應(yīng)、準(zhǔn)正交且完全非遞歸的信號(hào)分解和時(shí)頻分布估計(jì)方法，它可以將輸入信號(hào)分解為若干個(gè)具有帶寬限制的本征模態(tài)函數(shù)(intrinsic mode function,IMF),并使每一組IMF集中在不同的中心頻率附近，具有良好的噪聲魯棒性。VMD方法的基本思路是將信號(hào)的分解過程轉(zhuǎn)化為變分問題求解，可以劃分為變分問題的構(gòu)造和求解兩部分。

先對(duì)每個(gè)IMF分量做希爾伯特變換，具體如下

(1)

式中：δ(t)為脈沖函數(shù)；uk(t)是第k個(gè)IMF分量。

再將每個(gè)模態(tài)調(diào)制到相應(yīng)的基頻帶，即將預(yù)估的中心頻率添加到每個(gè)IMF分量的分析信號(hào)中，如下所示

(2)

式中,ωk為每個(gè)分量估計(jì)的中心頻率。

通過求解調(diào)制信號(hào)梯度L2范數(shù)的平方構(gòu)造如下約束變分問題

(3)

式中,f為原始輸入信號(hào)。

為了把約束變分問題轉(zhuǎn)化為無(wú)約束問題，引入了二次懲罰因子α和拉格朗日乘子λ(t)，得到擴(kuò)展拉格朗日表達(dá)式

(4)

利用乘法算子交替方向法搜尋上述變分問題的最優(yōu)解，得到uk、ωk及λ的更新公式為

(5)

(6)

(7)

式中：n為迭代次數(shù)；∧表示傅里葉變換；τ為噪聲容限系數(shù)。

迭代上述式(2)～(4)，直到滿足收斂條件

(8)

式中,ε為判別精度。

通常，在對(duì)信號(hào)進(jìn)行VMD分解前需要人為地設(shè)定4個(gè)參數(shù)，分別為分解層數(shù)K，二次懲罰因子α，噪聲容限系數(shù)τ和判別精度ε，其中τ和ε對(duì)VMD分解結(jié)果影響不大，一般可采用默認(rèn)值[12]。而K和α是決定VMD分解能否達(dá)到良好效果的兩個(gè)重要參數(shù)，若K值選擇不當(dāng)會(huì)造成欠分解或過分解現(xiàn)象，若α值選擇不當(dāng)會(huì)影響模態(tài)分量帶寬的大小，因此有必要尋找最優(yōu)的參數(shù)組合。

1.2 灰狼優(yōu)化算法(GWO)

GWO是一種新型智能搜索算法，具有原理簡(jiǎn)單、搜索能力強(qiáng)、回避局部最優(yōu)等優(yōu)點(diǎn)，目前已被證明該算法與傳統(tǒng)的遺傳算法、粒子群算法和蟻群算法相比，在收斂精度和收斂速度方面更具優(yōu)勢(shì)[13]。在GWO算法中，按照評(píng)估候選解的適應(yīng)度優(yōu)劣程度從高到低分為α,β,δ,ω4個(gè)等級(jí)，狼群的位置通過計(jì)算與獵物之間的距離不斷更新和進(jìn)化，最終實(shí)現(xiàn)狩獵，具體操作如下

(1) 圍捕策略

D=|C·Xp(t)-X(t)|

(9)

X(t+1)=Xp(t)-A·D

(10)

C=2r1

(11)

A=2a·r2-a

(12)

(2) 位置更新

Dk=|Ci·Xk(t)-X(t)|

(13)

Xi(t)=Xk(t)-Ai·Dk

(14)

(15)

式中：k=α,β,δ代表候選解評(píng)估排名前三的優(yōu)質(zhì)狼；i=1,2,3代表使用α,β,δ更新的其它個(gè)體狼；X(t+1)為迭代后狼群的位置。

2 所提方法

2.1 AVMD算法

由于分解層數(shù)K和二次懲罰因子α對(duì)VMD分解結(jié)果有顯著影響，若參數(shù)選擇不當(dāng)會(huì)導(dǎo)致模態(tài)分量丟失或模態(tài)混疊現(xiàn)象[14]，因此本文將GWO算法引入VMD的參數(shù)優(yōu)化過程，自適應(yīng)地獲得最佳參數(shù)組合，解決了依賴經(jīng)驗(yàn)或先驗(yàn)知識(shí)人為設(shè)定的問題。

適應(yīng)度的選擇是優(yōu)化算法的核心問題，由于包絡(luò)熵(envelope entropy,EE)可以反映信號(hào)的稀疏性，故障脈沖越規(guī)律則包絡(luò)熵越小[15]。因此，本文將最小平均包絡(luò)熵(minimum average envelope entropy,MAEE )為AVMD算法的目標(biāo)函數(shù)，定義如下

(16)

(17)

(18)

根據(jù)前面的理論分析，AVMD算法的基本流程如圖1所示。

圖1 AVMD算法流程圖

2.2 有效加權(quán)稀疏峭度指標(biāo)

利用AVMD算法可將原信號(hào)分解成K個(gè)IMF模態(tài)分量，這些分量有些包含了豐富的故障信息，有些則是噪聲干擾，為了便于分離，常采用峭度準(zhǔn)則或相關(guān)系數(shù)準(zhǔn)則進(jìn)行識(shí)別[16]。然而，峭度雖然能反映沖擊信號(hào)的分布密度但忽略了振幅大、分布分散的分量；而相關(guān)系數(shù)雖然能表達(dá)分量對(duì)原信號(hào)的復(fù)現(xiàn)程度卻易受噪聲干擾，這說(shuō)明單一的指標(biāo)并不能完全詮釋信號(hào)的振動(dòng)特性。因此，本文提出了一個(gè)新的評(píng)價(jià)指標(biāo)——有效加權(quán)稀疏峭度指標(biāo)(efficient weighted sparseness kurtosis index,EWSK),定義如下

WSK=Spa·Kur·Cor

(19)

(20)

(21)

(22)

(23)

按式(23)計(jì)算所得的EWSK大于0，則認(rèn)為是有效模態(tài)分量，否則為噪聲模態(tài)分量。EWSK結(jié)合了稀疏度、峭度、相關(guān)系數(shù)的優(yōu)點(diǎn)，使重構(gòu)信號(hào)不僅保留了原始信號(hào)的故障信息而且有效地濾除了噪聲。

2.3 包絡(luò)諧噪比(EHNR)

EHNR是指對(duì)數(shù)據(jù)樣本做包絡(luò)處理后諧波能量與噪聲能量的比值，它的大小不僅可以表征故障的強(qiáng)弱，還可以充分反映故障引起的周期性沖擊，是一種用于滾動(dòng)軸承故障診斷和預(yù)測(cè)的新指標(biāo)。

(24)

(25)

(26)

(27)

式中,τ為自相關(guān)函數(shù)的時(shí)間滯后量。

確定原始信號(hào)的自相關(guān)函數(shù)在時(shí)間滯后域中的最大位置，即EHNR定義如下

(28)

式中：τmax為rx(τ)自相關(guān)函數(shù)的最大位置；rx(τmax)為最大位置對(duì)應(yīng)的幅值；rx(0)為包絡(luò)總能量。

2.4 基于AVMD-EHNR的滾動(dòng)軸承早期退化檢測(cè)方法

滾動(dòng)軸承的全壽命運(yùn)行周期通?？煞譃檎Ｆ凇⑤p微退化、中度退化、嚴(yán)重退化和失效期5個(gè)階段，其中退化期是故障從產(chǎn)生到擴(kuò)展的劣化過程，若能檢測(cè)到退化的起點(diǎn)，不僅有助于早期故障診斷防止缺陷進(jìn)一步擴(kuò)大，而且為剩余壽命預(yù)測(cè)機(jī)制的及時(shí)啟動(dòng)提供依據(jù)，具有重要的工程意義和經(jīng)濟(jì)價(jià)值。

AVMD-EHNR方法的實(shí)現(xiàn)過程如圖2所示。首先，將原始信號(hào)進(jìn)行AVMD分解，采用EWSK指數(shù)從K個(gè)模態(tài)分量中選擇有效IMF分量并重構(gòu)信號(hào)，然后計(jì)算重構(gòu)信號(hào)的EHNR，由高斯分布的4σ準(zhǔn)則[17]，設(shè)置報(bào)警閾值為正常期時(shí)相對(duì)穩(wěn)定的EHNR均值μ加上4倍標(biāo)準(zhǔn)差σ，超出這個(gè)區(qū)間的數(shù)據(jù)則為異常點(diǎn)，根據(jù)報(bào)警閾值判斷軸承是否進(jìn)入退化期，并找出退化的起始點(diǎn)，實(shí)現(xiàn)軸承早期退化檢測(cè)。

圖2 AVMD-EHNR實(shí)現(xiàn)過程

3 仿真驗(yàn)證

通過構(gòu)造一組仿真信號(hào)來(lái)驗(yàn)證所提AVMD方法及EWSK指標(biāo)的有效性，如式(29)所示

(29)

式中:x1(t)為低頻信號(hào)；x2(t)為間歇性信號(hào)，用來(lái)模擬故障脈沖；x3(t)為高頻信號(hào)；x4(t)為均值為0方差為0.1的隨機(jī)噪聲；仿真信號(hào)的采樣頻率為1 000 Hz，采樣點(diǎn)數(shù)為1 200，時(shí)間為1 s，如圖3所示。

圖3 仿真信號(hào)及其分量

首先，利用傳統(tǒng)的EMD和EEMD方法對(duì)原信號(hào)進(jìn)行分解，各分量的時(shí)頻信息分別如圖4和圖5所示。從圖中可以看出，EMD和EEMD將原信號(hào)分別分解成7個(gè)分量和6個(gè)分量，都存在不同程度的模態(tài)混疊現(xiàn)象而且沒有分離出高頻分量并都有模態(tài)重復(fù)問題。

(a) 時(shí)域結(jié)果

(a) 時(shí)域結(jié)果

下面采用所提AVMD方法對(duì)原信號(hào)進(jìn)行分解，設(shè)置GWO狼群數(shù)量為20，最大迭代次數(shù)為10。根據(jù)MAEE原則得到的K和α的最佳組合為〈3,1 556〉，再對(duì)原信號(hào)進(jìn)行分解，如圖7所示。原信號(hào)被分解成中心頻率分別為42 Hz，113 Hz和250 Hz的3個(gè)分量，與原信號(hào)對(duì)比可以看出，不僅避免了模態(tài)混疊問題還濾除了噪聲分量，充分體現(xiàn)了AVMD方法的優(yōu)越性。

接下來(lái)，利用所提的EWSK指標(biāo)選擇有效IMF分量，如表1所示。為了驗(yàn)證EWSK指標(biāo)的有效性，表1中的前2行分別列出了各分量的峰值信噪比(peak signal to noise ratio,PSNR)和均方誤差(mean square error,MSE)，這是用來(lái)評(píng)價(jià)重構(gòu)信號(hào)質(zhì)量的常用指標(biāo)，見文獻(xiàn)[15]。從表中可以看出，IMF2的PSNR最高(基準(zhǔn)為30 dB，越高越好)且MSE最小，這說(shuō)明該分量包含了豐富的故障信息和原信號(hào)的波形和頻譜中尖峰特征。同時(shí)，表1中還列出了各分量的EWSK計(jì)算結(jié)果，選擇其值大于0的IMF2為有效特征分量，而其值小于0的IMF1和IMF3為偽分量。

表1 有效IMF分量的選擇

然后，采用固定參數(shù)的VMD方法對(duì)原信號(hào)進(jìn)行分解，選擇K和α分別為4和2 000，分解結(jié)果如圖6所示。從圖6中可以看出，由于沒有合理選擇VMD的參數(shù)，不僅出現(xiàn)了過分解問題，且原信號(hào)中的3個(gè)中心頻率也沒有準(zhǔn)確提取出來(lái)。

(a) 時(shí)域結(jié)果

(a) 時(shí)域結(jié)果

由于所提EWSK準(zhǔn)則憑借故障呈現(xiàn)的統(tǒng)計(jì)特性及其值的正負(fù)便能準(zhǔn)確分離出有效的IMF分量，無(wú)需故障特征的先驗(yàn)知識(shí)，所以即使在缺乏故障精確信息的不利條件下也能很好地完成分離任務(wù)。

4 案例分析

4.1 滾動(dòng)軸承全壽命實(shí)驗(yàn)數(shù)據(jù)

本文使用的實(shí)驗(yàn)數(shù)據(jù)來(lái)自辛辛那提大學(xué)智能維護(hù)中心(IMS)提供的滾動(dòng)軸承全壽命數(shù)據(jù)集[18]，實(shí)驗(yàn)裝置如圖8所示。

(a) 測(cè)試系統(tǒng)

電機(jī)轉(zhuǎn)速為2 000 r/min通過皮帶驅(qū)動(dòng)軸轉(zhuǎn)動(dòng)，軸上安裝了4個(gè)相同的ZA-2115型滾動(dòng)軸承，施加26 690 N的徑向載荷，由PCB353B33型加速度傳感器收集振動(dòng)信號(hào)，采樣頻率為20 kHz，每10 min采集1次。4個(gè)軸承均采用強(qiáng)制潤(rùn)滑，通過磁塞從油循環(huán)中收集碎屑，當(dāng)碎屑超過一定水平時(shí)，開關(guān)關(guān)閉實(shí)驗(yàn)停止。數(shù)據(jù)集共有3個(gè)測(cè)試案例，表2記錄了失效軸承的編號(hào)及故障類型，全壽命振動(dòng)數(shù)據(jù)如圖9所示。

表2 全壽命實(shí)驗(yàn)案例簡(jiǎn)要信息

(a) 案例1

4.2 振動(dòng)信號(hào)處理

以案例2中的100個(gè)樣本數(shù)據(jù)為例，設(shè)置GWO的狼群數(shù)量為20，最大迭代次數(shù)為10，模態(tài)分量K的取值范圍為[2,12]，二次懲罰因子α的范圍為[1 000,3 000]。搜索結(jié)果如圖10所示，迭代了4次就已收斂，此時(shí)對(duì)應(yīng)的K和α分別為6和2 048，并與K和α分別取4和2 000的固定參數(shù)VMD進(jìn)行對(duì)比，如圖11所示。在圖11(a)中丟失了0～1 000 Hz和5 000～6 000 Hz范圍的模態(tài)，屬于欠分解。在圖11(b)中，模態(tài)分量在各頻帶分布均勻且沒有出現(xiàn)模態(tài)混疊現(xiàn)象，說(shuō)明振動(dòng)信號(hào)沒有被欠分解和過分解。再根據(jù)EWSK指標(biāo)選擇IMF3和IMF4為有效IMF分量，如表3所示。重構(gòu)信號(hào)的PSNR和MSE分別為73.139 2 dB和0.003 2，說(shuō)明保留了大量原信號(hào)故障信息。

圖10 案例2的GWO收斂曲線

(a) 固定參數(shù)VMD分解后的頻率譜

表3 案例2每個(gè)模態(tài)分量的EWSK值

按照上述方法，分別對(duì)案例1和案例3中的數(shù)據(jù)進(jìn)行處理，經(jīng)GWO優(yōu)化后得到的最優(yōu)參數(shù)組合分別為〈6,2 078〉和〈6,2 114〉，選取有效分量分別為IMF3、4、5和IMF1、4，對(duì)有效分量進(jìn)行信號(hào)重構(gòu)，這樣既去除了噪聲干擾又保留了豐富的故障信息。

4.3 早期故障點(diǎn)檢測(cè)

為了驗(yàn)證EHNR指標(biāo)對(duì)軸承早期故障的敏感性，先對(duì)案例2中的軸承1的全壽命數(shù)據(jù)采用有效值和峭度指標(biāo)分析，如圖12所示。從圖12中可以看出，有效值雖然能表征軸承退化過程的整體趨勢(shì)，但對(duì)早期故障并不敏感。而峭度指標(biāo)對(duì)較大的干擾沖擊敏感，并隨著沖擊的密集這種敏感會(huì)逐漸失效，兩種指標(biāo)都很難確定退化的起始點(diǎn)。

(a) 有效值指標(biāo)

然后再根據(jù)所提EHNR指標(biāo)進(jìn)行分析，如圖13所示。從圖13(a)中可以明顯識(shí)別出軸承1從運(yùn)行到失效的5個(gè)階段，這說(shuō)明該指標(biāo)趨勢(shì)性好。1為正常期，歷經(jīng)的時(shí)間最長(zhǎng)(約為89 h),這一階段軸承處于穩(wěn)定運(yùn)行狀態(tài)，EHNR的值僅有微小的波動(dòng)，選擇這一時(shí)期EHNR的均值加上4倍標(biāo)準(zhǔn)差確定報(bào)警閾值，如圖13中紅線所示；2為輕度退化期，從EHNR超越報(bào)警閾值開始累計(jì)，此時(shí)的故障特征往往被強(qiáng)噪聲掩蓋，而首個(gè)認(rèn)定的異常點(diǎn)即為故障的起點(diǎn)，如圖13(b)所示；3為中度退化期，是故障擴(kuò)展階段，這一階段故障的尺寸會(huì)加深加大，故障特征也變得更加明顯，但也會(huì)由于連續(xù)的接觸磨損使故障部位變平滑，出現(xiàn)“自愈”的假象，EHNR的幅值出現(xiàn)先升高后降低的現(xiàn)象；4為嚴(yán)重退化期，此時(shí)軸承故障缺陷增長(zhǎng)加快，EHNR幅值劇烈振蕩，工作的可靠度顯著降低，軸承瀕臨失效；5為失效期，此時(shí)軸承已完全不能可靠工作，甚至?xí)?duì)整個(gè)機(jī)械系統(tǒng)造成嚴(yán)重危害。

(a) 軸承運(yùn)行狀態(tài)劃分

采用本文所提方法對(duì)原始信號(hào)進(jìn)行分解和重構(gòu)，再計(jì)算重構(gòu)信號(hào)的EHNR，如圖14所示。

圖14 案例2的AVMD-EHNR方法檢測(cè)結(jié)果

從圖14中可以看出，AVMD-EHNR方法識(shí)別的異常點(diǎn)為第536個(gè)樣本點(diǎn)，比原始信號(hào)識(shí)別的547點(diǎn)早了11個(gè)樣本(約為1.8 h),雖然兩個(gè)結(jié)果中均沒有出現(xiàn)虛警點(diǎn)，但是明顯看出所提方法得到的EHNR幅值在運(yùn)行前期更加平穩(wěn)，整體振蕩不大。

采用同樣方法對(duì)案例1和案例3中的軸承3的全壽命數(shù)據(jù)進(jìn)行分析對(duì)比，得到的結(jié)果如圖15和圖16所示，其中圖15(a)和圖16(a)分別表示原信號(hào)的檢測(cè)結(jié)果，圖15(b)和圖16(b)分別表示采用本文所提方法得到的檢測(cè)結(jié)果。再將本文所提AVMD-EHNR方法與文獻(xiàn)[11]所提的快速集成經(jīng)驗(yàn)?zāi)Ｊ椒纸?fast ensemble empirical mode decomposition,FEEMD)方法的檢測(cè)結(jié)果進(jìn)行比較，如表4所示。

(a) 案例1原信號(hào)的檢測(cè)結(jié)果

(a) 案例3原信號(hào)的檢測(cè)結(jié)果

表4 三種方法的早期退化點(diǎn)檢測(cè)結(jié)果對(duì)比

結(jié)合圖15、圖16及表4可以看出，案例1的原始信號(hào)EHNR值在運(yùn)行初期極不穩(wěn)定，造成虛警率很高，F(xiàn)EEMD-EHNR方法雖然解決了原始信號(hào)的虛警問題，但是檢測(cè)結(jié)果卻比本文方法晚了105個(gè)樣本點(diǎn)(約為17.5 h)；案例3的原始信號(hào)EHNR值在運(yùn)行初期有3處虛警點(diǎn)，而FEEMD-EHNR方法雖然將虛警點(diǎn)降為1個(gè)但檢測(cè)結(jié)果仍比本文方法晚(約為0.83 h)。經(jīng)過對(duì)比發(fā)現(xiàn)，本文所提方法在3個(gè)案例的檢測(cè)中均沒有虛警問題，且檢測(cè)的早期故障點(diǎn)比原始信號(hào)和FEEMD-EHNR方法都要早，這說(shuō)明AVMD-EHNR方法不僅可以有效濾除噪聲分量，而且可以更好地識(shí)別早期故障點(diǎn)，兼具魯棒性和敏感性，適合滾動(dòng)軸承在線運(yùn)行檢測(cè)和故障診斷。

5 結(jié) 論

為了檢測(cè)滾動(dòng)軸承的早期故障點(diǎn)，提出了一種基于自適應(yīng)變分模態(tài)分解(AVMD)和包絡(luò)諧噪比(EHNR)的方法。

(1) 采用GWO算法得到了VMD的最佳參數(shù)組合，解決了人為選擇參數(shù)困難及參數(shù)設(shè)置不當(dāng)帶來(lái)的模態(tài)分量丟失和模態(tài)混疊問題。

(2) 提出了一種有效加權(quán)稀疏峭度指數(shù)(EWSK)，該指標(biāo)可對(duì)故障模態(tài)分量和噪聲模態(tài)分量進(jìn)行區(qū)分，從而獲得含噪少且故障信息豐富的重構(gòu)信號(hào)。

(3) EHNR是一種新的故障特征提取技術(shù)，與傳統(tǒng)的有效值和峭度指標(biāo)相比，它對(duì)周期性故障脈沖具有更強(qiáng)的敏感性、魯棒性和趨勢(shì)性。

(4) 結(jié)合AVMD和EHNR的優(yōu)點(diǎn)對(duì)重構(gòu)信號(hào)進(jìn)行分析，可以剔除運(yùn)行初期的虛假異常點(diǎn)，有效避免“虛警”現(xiàn)象，有利于進(jìn)一步識(shí)別出真實(shí)的異常點(diǎn)，從而解決了軸承退化過程的起始點(diǎn)難以確定的問題，為滾動(dòng)軸承早期故障診斷和剩余壽命預(yù)測(cè)提供依據(jù)，有重要的工程意義。