基于EWT-OPRCMDE-ELM的風(fēng)電機(jī)組齒輪箱故障診斷研究

李 輝,李 宣,賈 嶸,羅興琦,白 亮

(1.西安理工大學(xué)電氣工程學(xué)院，陜西 西安 710048；2.西安理工大學(xué)水利水電學(xué)院，陜西 西安 710048)

0 引言

惡劣的工作環(huán)境導(dǎo)致風(fēng)機(jī)故障頻發(fā)[1]。齒輪箱作為風(fēng)電機(jī)組傳動(dòng)系統(tǒng)的主要設(shè)備，更是故障的高發(fā)部位。據(jù)統(tǒng)計(jì)，齒輪箱故障占機(jī)組總故障60%以上[2]，而一臺(tái)風(fēng)力機(jī)故障停機(jī)時(shí)間的20%是由齒輪箱故障引起的。一旦發(fā)生故障，機(jī)組將面臨長(zhǎng)時(shí)間停機(jī)和昂貴的維修費(fèi)用，損失巨大[3]。因此，準(zhǔn)確、高效地對(duì)風(fēng)機(jī)齒輪箱進(jìn)行故障診斷，對(duì)于保障機(jī)組的安全、穩(wěn)定運(yùn)行具有重要意義。

由于齒輪箱故障振動(dòng)信號(hào)受背景噪聲、交變載荷和傳輸路徑等影響，呈現(xiàn)非平穩(wěn)、非線性的特點(diǎn)，通常用時(shí)頻分析方法對(duì)其處理，如小波變換、經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition,EMD)、變分模態(tài)分解(variational mode decomposition,VMD)等。但是小波變換中小波基的選擇不當(dāng)會(huì)對(duì)處理效果產(chǎn)生不良影響。EMD方法缺乏理論基礎(chǔ)且存在模態(tài)混疊、端點(diǎn)效應(yīng)等問(wèn)題，所以無(wú)法有效分析故障信號(hào)。VMD方法依賴(lài)于分解層數(shù)和懲罰因子的設(shè)定，如參數(shù)設(shè)置不當(dāng)，將極大地影響地處理效果。針對(duì)傳統(tǒng)信號(hào)處理方法的不足，Gilles[4]結(jié)合經(jīng)驗(yàn)?zāi)B(tài)分解的自適應(yīng)性和小波分析的框架理論，提出了經(jīng)驗(yàn)小波變換(empirical wavelet transform,EWT)方法,有效地改善了傳統(tǒng)時(shí)頻分析方法的弊端。文獻(xiàn)[5]將EWT應(yīng)用于旋轉(zhuǎn)機(jī)械故障診斷中，并證明EWT處理效果優(yōu)于傳統(tǒng)EMD方法。

高效地提取區(qū)分度大的穩(wěn)定特征是風(fēng)電機(jī)組齒輪箱故障診斷的重要環(huán)節(jié)。學(xué)者們通常提取故障信號(hào)的時(shí)域、頻域或時(shí)頻域特征，進(jìn)行故障診斷。文獻(xiàn)[6]提取軸承振動(dòng)信號(hào)的均值、峭度等時(shí)域特征，進(jìn)行故障診斷，取得較好的診斷效果。文獻(xiàn)[7]提取軸承振動(dòng)信號(hào)的時(shí)域、頻域和時(shí)頻域的多個(gè)混合特征參量，利用XGBoost進(jìn)行風(fēng)機(jī)軸承故障診斷。文獻(xiàn)[8]融合風(fēng)機(jī)齒輪箱振動(dòng)信號(hào)的時(shí)域、頻域和時(shí)頻域特征，特征降維后利用隨機(jī)森林算法進(jìn)行故障模式識(shí)別。但是，由于振動(dòng)信號(hào)受強(qiáng)烈噪聲干擾、多振動(dòng)源耦合和傳輸路徑等影響，導(dǎo)致基于傳統(tǒng)時(shí)域、頻域和時(shí)頻域特征提取的故障診斷方法效果不佳，故障診斷正確率和診斷時(shí)間受特征參量的維數(shù)影響較大。精細(xì)復(fù)合多尺度散布熵[9](refined composite multi-scale dispersion entropy,RCMDE)作為一種衡量時(shí)間序列復(fù)雜度的非線性動(dòng)力學(xué)參數(shù)，已被成功應(yīng)用于心電信號(hào)和腦電信號(hào)的分析中，并被證明在計(jì)算速度和穩(wěn)定性方面優(yōu)于樣本熵和排列熵。其特征提取效果優(yōu)于傳統(tǒng)時(shí)域、頻域特征。

極限學(xué)習(xí)機(jī)(extreme learning machine,ELM)[10]是一種泛化的單隱層前饋神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)算法。其輸入層和隱含層之間的連接權(quán)值和隱含層閾值隨機(jī)產(chǎn)生，訓(xùn)練過(guò)程無(wú)需迭代調(diào)整，與傳統(tǒng)反向傳播(back propagation,BP)神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)等智能分類(lèi)算法相比，輸入?yún)?shù)少、運(yùn)算速度更快。

為提高風(fēng)電機(jī)組齒輪箱故障診斷效果，本文提出一種基于EWT重構(gòu)、最優(yōu)參數(shù)精細(xì)復(fù)合多尺度散布熵和ELM結(jié)合的故障診斷方法，并將其應(yīng)用于風(fēng)電機(jī)組齒輪箱。

1 經(jīng)驗(yàn)小波變換

EWT是通過(guò)對(duì)信號(hào)的傅里葉頻譜進(jìn)行自適應(yīng)分割，根據(jù)分割邊界，構(gòu)建正交小波濾波器組以提取包含故障信息的調(diào)幅調(diào)頻分量。步驟如下。

①對(duì)故障信號(hào)進(jìn)行傅里葉變換，將信號(hào)的傅里葉頻譜定義在[0，π]范圍內(nèi)，并將[0，π]分成N個(gè)帶寬不等的頻帶，則每個(gè)頻帶可以表示為Λn=[ωn-1,ωn]。以ωn為中心，定義Tn=2τn為過(guò)渡帶。傅里葉軸分割如圖1所示。

圖1 傅里葉軸分割Fig.1 Partitioning of the Fourier axis

②經(jīng)驗(yàn)小波的尺度函數(shù)φn(ω)和小波函數(shù)ψn(ω)在頻域里的計(jì)算式如下：

(1)

(2)

③根據(jù)經(jīng)典小波變換的方法構(gòu)造經(jīng)驗(yàn)小波變換，細(xì)節(jié)系數(shù)由經(jīng)驗(yàn)小波函數(shù)內(nèi)積計(jì)算得出，近似系數(shù)由經(jīng)驗(yàn)尺度函數(shù)內(nèi)積計(jì)算得出：

(3)

(4)

原信號(hào)重構(gòu)表達(dá)式如下：

(5)

式中：∧表示傅里葉變換;∨表示傅里葉逆變換;—表示復(fù)共軛;?表示卷積。

經(jīng)驗(yàn)?zāi)B(tài)fk(t)定義如下：

(6)

2 精細(xì)復(fù)合多尺度散布熵

散布熵計(jì)算步驟如下。

①x={xj|j=1,2,…,N}為一時(shí)間序列，利用正態(tài)分布函數(shù)將時(shí)間序列x映射到y(tǒng)={yj|j=1,2,…,N}，yj∈(0,1)，即：

(7)

式中：μ、σ分別為期望和標(biāo)準(zhǔn)差。

②通過(guò)線性變換將y映射到{1,2，…，c}范圍內(nèi)：

(8)

(9)

式中：m、d分別為嵌入維數(shù)和時(shí)延。

④計(jì)算散布模式πv0v1…vm-1，v=1,2,…,c。

(10)

⑤計(jì)算每種散布模式的概率πv0v1…vm-1：

(11)

⑥根據(jù)香農(nóng)熵定義,有：

(12)

散步熵(dispersion entropy,DE)的值越大，信號(hào)時(shí)間序列不規(guī)則程度越高；反之，不規(guī)則程度越低。

RCMDE計(jì)算步驟如下。

RCMDE是對(duì)原始數(shù)據(jù)精細(xì)化和多尺度化，其作用是減小初始點(diǎn)位置對(duì)計(jì)算結(jié)果的影響和充分利用數(shù)據(jù)序列的信息。計(jì)算過(guò)程如下。當(dāng)尺度因子為S時(shí)：首先，將原始信號(hào)按初始點(diǎn)連續(xù)分割成長(zhǎng)度為S的小段，求取每一段的平均值，再將這些平均值按順序排列成S個(gè)粗?；蛄?；其次，計(jì)算每個(gè)粗?；蛄猩⒉寄Ｊ溅械母怕?，并計(jì)算這些散布模式概率的平均值；最后，根據(jù)香農(nóng)熵的定義計(jì)算RCMDE的值。

(13)

②對(duì)于每個(gè)尺度S，RCMDE值定義如下：

(14)

3 Relief-F與ELM

提取齒輪箱故障特征后，由于存在特征冗余，會(huì)使診斷時(shí)長(zhǎng)和精度受影響。故采用Relief-F算法[11]計(jì)算分類(lèi)權(quán)重，剔除冗余。Relief-F算法是從Relief算法擴(kuò)展而來(lái)的，可以處理多類(lèi)別分類(lèi)權(quán)重問(wèn)題。該算法計(jì)算過(guò)程如下。

①?gòu)乃袠颖局?，隨機(jī)取出一個(gè)樣本a。

②在與樣本a相同分類(lèi)的樣本組內(nèi)，取出k個(gè)最近鄰樣本。

③在所有其他與樣本a不同分類(lèi)的樣本組內(nèi)， 也分別取出k個(gè)最近鄰樣本。

④計(jì)算每個(gè)特征的權(quán)重。

極限學(xué)習(xí)機(jī)是一種單隱層前饋神經(jīng)網(wǎng)絡(luò)，由輸入層、隱含層和輸出層構(gòu)成。與傳統(tǒng)神經(jīng)網(wǎng)絡(luò)不同的是，極限學(xué)習(xí)機(jī)的輸入層與隱含層之間的連接權(quán)值和隱含層神經(jīng)元的閾值為隨機(jī)設(shè)定，在開(kāi)始創(chuàng)建網(wǎng)絡(luò)時(shí)隨機(jī)產(chǎn)生，且設(shè)置完成后無(wú)需調(diào)整。隱含層與輸出層之間的連接權(quán)值由方程組依次解出，無(wú)需迭代計(jì)算。故極限學(xué)習(xí)機(jī)可在保證學(xué)習(xí)精度的前提下極大地加快學(xué)習(xí)速度。

極限學(xué)習(xí)機(jī)訓(xùn)練過(guò)程如下。假設(shè)有N個(gè)任意樣本(Xi,ti)。其中，Xi=[xi1xi2…xin]T∈Rn，ti=[ti1ti2…tin]∈Rm，n為輸入層的維度，m為輸出層的維度。

當(dāng)隱層節(jié)點(diǎn)的個(gè)數(shù)為L(zhǎng)時(shí)，單隱層前饋神經(jīng)網(wǎng)絡(luò)輸出可以表示為：

(15)

式中：i=1,2,…,N;Wj=[wj1wj2…wjn]T為輸入權(quán)值;βj=[βj1βj2…βjm]T為輸出權(quán)重;bj是第j個(gè)隱層單元的偏置;WjXj為Wj與Xj的內(nèi)積;g()為激活函數(shù)。

單隱層前饋神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)的目標(biāo)是在最小的誤差下逼近N個(gè)樣本，可以表示為：

(16)

即存在Wj、βj、bj，使得：

(17)

式(17)可簡(jiǎn)化為Hβ=T。其中：H為隱層節(jié)點(diǎn)輸出矩陣；β為輸出權(quán)重；T為期望輸出。

(18)

1,2,…,L

(19)

當(dāng)輸入權(quán)重Wj和隱含層偏置bj隨機(jī)確定后，隱層的輸出矩陣H就會(huì)被確定。此時(shí)，單隱層神經(jīng)網(wǎng)絡(luò)求解的問(wèn)題被轉(zhuǎn)換成線性系統(tǒng)Hβ=T的求解，則輸出權(quán)重β可由式(20)確定。

(20)

極限學(xué)習(xí)機(jī)診斷過(guò)程如下。

①設(shè)定激活函數(shù)g()和隱層神經(jīng)元個(gè)數(shù)N。

②輸入訓(xùn)練集進(jìn)行網(wǎng)絡(luò)訓(xùn)練。

③利用訓(xùn)練好的ELM，輸入測(cè)試集進(jìn)行故障診斷。

4 故障診斷流程

故障診斷流程如下。

①利用振動(dòng)傳感器采集風(fēng)電機(jī)組齒輪箱軸承故障振動(dòng)信號(hào)。

②利用EWT分解原始振動(dòng)信號(hào)，將得到若干子模態(tài)分量EWF。

③計(jì)算各子模態(tài)分量與原始信號(hào)的相關(guān)系數(shù)Rn。設(shè)定相關(guān)系數(shù)閾值為0.3，選取大于閾值的子模態(tài)分量EWF進(jìn)行信號(hào)重構(gòu)，提取故障主要信息。

相關(guān)系數(shù)是一種度量2個(gè)反映之間相關(guān)程度的系數(shù)，其值介于1和-1之間。其中：1表示變量完全正相關(guān)；0表示無(wú)關(guān)；-1表示完全負(fù)相關(guān)。皮爾遜相關(guān)系數(shù)定義為兩個(gè)變量之間的協(xié)方差和標(biāo)準(zhǔn)差的商。函數(shù)表達(dá)式為：

(21)

④以重構(gòu)信號(hào)RCMDE偏度值的平方函數(shù)為適應(yīng)度函數(shù)，計(jì)算重構(gòu)信號(hào)最優(yōu)RCMDE參數(shù)，提取重構(gòu)信號(hào)最優(yōu)參數(shù)RCMDE(optimal parameter RCMDE,OPRCMDE)構(gòu)成故障特征向量。

⑤利用Relief-F計(jì)算OPRCMDE各尺度特征分類(lèi)權(quán)重，并進(jìn)行特征參量篩選，剔除冗余特征。

⑥將篩選得到的特征向量作為輸入，利用ELM進(jìn)行故障診斷。

5 信號(hào)仿真驗(yàn)證

為驗(yàn)證EWT分解在處理復(fù)雜信號(hào)時(shí)的優(yōu)越性。采用仿真信號(hào)進(jìn)行驗(yàn)證。仿真信號(hào)f由指數(shù)信號(hào)f1、調(diào)頻信號(hào)f2、噪聲信號(hào)f3(強(qiáng)度為0.2 dB·W)構(gòu)成。仿真信號(hào)時(shí)域如圖2所示。

圖2 仿真信號(hào)時(shí)域圖Fig.2 Time domain diagram of simulation signal

利用EWT分解信號(hào)f，同時(shí)與EMD方法對(duì)比。EWT分解結(jié)果時(shí)域如圖3所示。EWT分解結(jié)果的Hilbert譜如圖4所示。

圖3 EWT分解結(jié)果時(shí)域圖Fig.3 Time domain diagram of EWT decomposition results

圖4 EWT分解結(jié)果的Hilbert譜圖Fig.4 Hilbert spectrum of EWT decomposition results

(22)

根據(jù)圖3和圖4，由 EWT分解結(jié)果可得，EWT將模擬信號(hào)f分解為4個(gè)子模態(tài)分量。 其中，EWF1在[0,1]上單調(diào)增加，最大值為5，對(duì)應(yīng)于f1的變化趨勢(shì)。 EWF2的周期約為0.1 s，其波形和周期對(duì)應(yīng)于f2中的sin(20πt)。而EWF3的周期約為0.05 s，對(duì)應(yīng)于f2中cos(40πt)。 EWF4沒(méi)有明顯的特征，并且在時(shí)域中很雜亂，因此被判斷為噪聲信號(hào)。 同時(shí)，在Hilbert頻譜中可以觀察到清晰、穩(wěn)定的10 Hz和20 Hz對(duì)應(yīng)于f2的特征頻率。EMD分解結(jié)果時(shí)域如圖5所示。

圖5 EMD分解結(jié)果時(shí)域圖Fig.5 Time domain diagram of EMD decomposition results

EMD分解結(jié)果的Hilbert譜圖如圖6所示。觀察圖5和圖6，EMD將模擬信號(hào)分解為8個(gè)子模態(tài)分量，但只有Res分量具有與其相同的變化趨勢(shì)，并且很難清楚地觀察到其余IMF分量中與模擬信號(hào)分量相對(duì)應(yīng)的分量。 而且，由于噪聲干擾，Hilbert頻譜中的特征頻率不規(guī)律且雜亂混疊，并且難以觀察到穩(wěn)定的模擬信號(hào)中包含的特征頻率。

圖6 EMD分解結(jié)果的Hilbert譜圖Fig.6 Hilbert spectrum of EMD decomposition results

上述仿真試驗(yàn)中，通過(guò)比較EWT和EMD的分解結(jié)果，可以發(fā)現(xiàn)EMD方法會(huì)過(guò)度分解信號(hào)，并且存在模態(tài)混疊現(xiàn)象，不能有效地提取故障特征分量。 EWT可以在嘈雜的環(huán)境中有效地提取信號(hào)的主要成分，以獲得更好的特征提取結(jié)果。

6 試驗(yàn)分析

本節(jié)結(jié)合了EWT信號(hào)重構(gòu)和OPRCMDE，并使用極限學(xué)習(xí)機(jī)進(jìn)行故障診斷。采用基于千鵬 QPZZ-Ⅱ 旋轉(zhuǎn)機(jī)械振動(dòng)及故障模擬試驗(yàn)平臺(tái)的齒輪箱振動(dòng)數(shù)據(jù)進(jìn)行分析。試驗(yàn)平臺(tái)由變速驅(qū)動(dòng)電機(jī)(0.75 kW)、軸承、齒輪箱、軸、偏重轉(zhuǎn)盤(pán)、調(diào)速器等組成。通過(guò)調(diào)節(jié)配重、調(diào)節(jié)部分的安裝位置以及組件的有機(jī)組合，快速模擬各種齒輪箱故障。

齒輪箱參數(shù)如下。

①齒輪。大齒輪：模數(shù)2、齒數(shù)75、材質(zhì)S45C(3只，其中2只備件齒輪)；小齒輪：模數(shù)2、齒數(shù)55、材質(zhì)S45C(2只，其中1只備件齒輪)。

②潤(rùn)滑：浸油式。

③軸承：滾動(dòng)軸承(N205)。

設(shè)定轉(zhuǎn)速為880 r/min，采樣頻率為5 120 Hz，通過(guò)安裝在齒輪箱輸出軸負(fù)載側(cè)軸承X方向的加速度傳感器，采集齒輪箱正常、點(diǎn)蝕、磨損和斷齒這4類(lèi)故障狀態(tài)的振動(dòng)信號(hào)各25組樣本。每個(gè)樣本長(zhǎng)度為2 048。

6.1 EWT分解與重構(gòu)

以齒輪箱齒輪磨損故障為例。EWT將磨損故障信號(hào)分解為10個(gè)EWF。根據(jù)本文所提到的相關(guān)系數(shù)閾值，選取EWF2、EWF6為最優(yōu)分量重構(gòu)故障信號(hào)并進(jìn)行下一步的故障特征提取及模式識(shí)別。

EWF相關(guān)系數(shù)如圖7所示。

圖7 EWF相關(guān)系數(shù)Fig.7 Correlation coefficient of EWF

6.2 OPRCMDE參數(shù)計(jì)算及特征提取

選取合適的RCMDE參數(shù)，能夠有效提取區(qū)分度較大的故障特征向量，同時(shí)減少計(jì)算時(shí)間。針對(duì)不同的故障類(lèi)型選取對(duì)應(yīng)的最佳參數(shù)，可獲取最佳的故障特征區(qū)分度。

在精細(xì)復(fù)合多尺度散布熵的計(jì)算中，涉及4個(gè)參數(shù)的選取，即尺度因子S、嵌入維數(shù)m、類(lèi)別個(gè)數(shù)C和時(shí)延d。為充分利用軸承故障數(shù)據(jù)，設(shè)置尺度因子S=15，時(shí)延d對(duì)算法影響較小，故設(shè)置為1。嵌入維數(shù)m、類(lèi)別個(gè)數(shù)C的設(shè)定對(duì)特征提取效果影響重大。m過(guò)大，將導(dǎo)致RCMDE運(yùn)算時(shí)間加長(zhǎng)；m過(guò)小，則算法檢測(cè)信號(hào)突變性減弱。C過(guò)大，將導(dǎo)致RCMDE對(duì)噪聲十分敏感；C過(guò)小，則彼此相距較遠(yuǎn)的兩個(gè)振幅可能被分到相同的類(lèi)別內(nèi)。考慮到m和C之間的交互作用，使用網(wǎng)格搜索算法對(duì)范圍內(nèi)的最優(yōu)m值和C值進(jìn)行同步搜索。

為分析多尺度熵值的集中趨勢(shì)，可計(jì)算其均值。但僅憑均值不足以表征多尺度熵值的總體趨勢(shì)。本文引入精細(xì)復(fù)合多尺度散布熵的偏度。其絕對(duì)值越大，表明均值的效能越有問(wèn)題；其絕對(duì)值越小，表明均值越可信賴(lài)。本文以EWT重構(gòu)信號(hào)精細(xì)復(fù)合多尺度散布熵偏度的平方函數(shù)[12]作為適應(yīng)度函數(shù)，求取其最小值。偏度計(jì)算過(guò)程如下：

信號(hào)序列X={xi},i=1,2,…,N。15個(gè)尺度下的精細(xì)復(fù)合多尺度散布熵值組成序列：

BP(X)={BP(1),BP(2),…,BP(15)}

(23)

偏度值為：

(24)

適應(yīng)度函數(shù)取為：

(25)

設(shè)定m搜索范圍為[2,6]、C搜索范圍為[3,8]、設(shè)定步長(zhǎng)均為1。通過(guò)網(wǎng)格搜索算法，得到4類(lèi)故障OPRCMDE參數(shù)，如表1所示。

表1 4類(lèi)故障OPRCMDE參數(shù)

提取4類(lèi)故障狀態(tài)EWT重構(gòu)信號(hào)的OPRCMDE值，得到4組25×15故障特征矩陣。原始信號(hào)推薦參數(shù)RCMDE如圖8所示(4類(lèi)故障狀態(tài)均設(shè)定m=3、C=6)。EWT重構(gòu)信號(hào)OPRCMDE如圖9所示。

圖8 原始信號(hào)推薦參數(shù)RCMDEFig.8 RCMDE of recommended parameters of original signal

圖9 EWT重構(gòu)信號(hào)OPRCMDEFig.9 OPRCMDE of the EWT reconstructed signal

對(duì)比圖8、圖9，可觀察到本文方法所提取的4種故障狀態(tài)特征區(qū)別明顯，有利于后期的故障特征提取和故障診斷。

由于高維OPRCMDE特征向量存在冗余，利用Relief-F算法(參數(shù)設(shè)置迭代次數(shù)M=20,近鄰值KL=5)計(jì)算每個(gè)尺度因子下OPRCMDE值的特征權(quán)重，并進(jìn)行特征向量二次篩選和約簡(jiǎn)。各尺度因子分類(lèi)特征權(quán)重如圖10所示。

圖10 各尺度因子分類(lèi)特征權(quán)重Fig.10 Feature weights at each scale

根據(jù)權(quán)重大小，選取第1、第2、第6共3個(gè)尺度因子下的OPRCMDE值作為故障特征向量，約簡(jiǎn)后得到4組25×3的故障特征矩陣。Relief-F約簡(jiǎn)特征散點(diǎn)圖如圖11所示。觀察圖11可發(fā)現(xiàn)：各類(lèi)故障間聚集緊密從而可將4種不同故障狀態(tài)明顯地區(qū)分開(kāi)。

圖11 Relief-F約簡(jiǎn)特征散點(diǎn)圖Fig.11 Relief-F reduced feature scatter

6.3 極限學(xué)習(xí)機(jī)故障診斷

為4故障狀態(tài)貼上標(biāo)簽，正常、點(diǎn)蝕、磨損、斷齒的標(biāo)簽分別為1、2、3、4。選取60%樣本作為訓(xùn)練樣本，剩余40%樣本作為測(cè)試樣本。測(cè)試集40個(gè)樣本分類(lèi)結(jié)果與實(shí)際類(lèi)別一致，故障診斷正確率100%，對(duì)應(yīng)類(lèi)別未出現(xiàn)錯(cuò)誤分類(lèi)。ELM診斷結(jié)果如圖12所示。

圖12 ELM診斷結(jié)果Fig.12 Results of ELM diagnosis

根據(jù)試驗(yàn)結(jié)果，證明本文提出的方法能夠提取風(fēng)機(jī)齒輪箱區(qū)分度較大的故障特征，實(shí)現(xiàn)了對(duì)齒輪箱不同故障的準(zhǔn)確診斷。為進(jìn)一步說(shuō)明本文方法的優(yōu)越性，將本文方法與多尺度排列熵(multi-scale permutation entropy,MPE)-ELM(MPE參數(shù)設(shè)置為m=5、S=12、T=1)、MPE-ReliefF-ELM、RCMDE-ELM(RCMDE參數(shù)推薦設(shè)置為m=3、C=6)、時(shí)頻域融合特征-ReliefF-ELM(時(shí)頻域融合特征共11個(gè)，其中8個(gè)時(shí)域特征分別為平均值、有效值、峰值、峰值因子、峭度、波形因子、脈沖因子、裕度因子，3個(gè)頻域特征分別為重心頻率、均方頻率、頻率方差)進(jìn)行對(duì)比。每種方法運(yùn)行10次，取診斷正確率平均值。各診斷模型診斷效果對(duì)比如表2所示。

表2 各診斷模型診斷效果對(duì)比

由表2中結(jié)果可知，當(dāng)本文方法所提取的齒輪箱故障特征應(yīng)用于齒輪箱故障診斷時(shí)，診斷準(zhǔn)確率更高且更加穩(wěn)定。

7 結(jié)論

本文針對(duì)傳統(tǒng)風(fēng)電機(jī)組齒輪箱故障診斷方法的不足，提出一種基于EWT重構(gòu)、最優(yōu)參數(shù)精細(xì)復(fù)合多尺度散布熵和極限學(xué)習(xí)機(jī)結(jié)合的故障診斷方法。具體結(jié)論如下。

①EWT解決了傳統(tǒng)時(shí)頻處理方法存在的模態(tài)混疊、端點(diǎn)效應(yīng)、參數(shù)設(shè)置等問(wèn)題。通過(guò)相關(guān)系數(shù)選取子模態(tài)分量進(jìn)行信號(hào)重構(gòu)，提取了故障的主要信息。

②為提高精細(xì)多尺度散布熵算法的特征提取性能，以其偏度值的平方函數(shù)作為適應(yīng)度函數(shù)，通過(guò)網(wǎng)格搜索算法同步搜索計(jì)算2個(gè)關(guān)鍵參數(shù)(即m和C)。通過(guò)試驗(yàn)對(duì)比，證明EWT重構(gòu)信號(hào)最優(yōu)參數(shù)精細(xì)復(fù)合多尺度散布熵在提取各類(lèi)故障特征時(shí)區(qū)分度更好，診斷結(jié)果更穩(wěn)定準(zhǔn)確。

③針對(duì)特征向量冗余和一般分類(lèi)算法參數(shù)多且參數(shù)設(shè)定影響分類(lèi)準(zhǔn)確率的問(wèn)題，采用Relief-F算法計(jì)算特征向量的分類(lèi)權(quán)重，選擇權(quán)重大者構(gòu)成最終的特征向量，剔除了冗余特征。最后，再利用運(yùn)算速度快、參數(shù)設(shè)置少的ELM進(jìn)行故障診斷。試驗(yàn)結(jié)果證明，該方法能夠有效運(yùn)用于風(fēng)電機(jī)組齒輪箱故障診斷中。