周志朋 謝冬梅

摘 要：針對風電齒輪箱的復合故障診斷問題，本文提出了基于改進變分模態(tài)的風電齒輪箱復合故障的特征提取方法，并采用改進鯨魚算法優(yōu)化的最小二乘支持向量機（DEWOA-LSSVM）故障診斷模型進行故障診斷。通過變分模態(tài)方法將收集的信號分解為K個模態(tài)分量，通過加權(quán)排序熵對其進行量化處理，得到特征向量，將特征向量輸入到經(jīng)過改進鯨魚算法優(yōu)化后的最小二乘支持向量機模型完成故障診斷，并通過試驗驗證該方法的有效性。

關鍵詞：故障診斷;變分模態(tài)分解;加權(quán)排序熵;鯨魚算法;最小二乘支持向量機

中圖分類號：TH165.3文獻標識碼：A文章編號：1003-5168（2020）17-0072-04

Research on Compound Fault Diagnosis Method of Wind Power Gearbox

ZHOU Zhipeng1 XIE Dongmei2

（1. Graduate Department of Shenyang Institute of Engineering，Shenyang Liaoning 110136;2. School of Electric Power， Shenyang Institute of Engineering，Shenyang Liaoning 110136）

Abstract： Aiming at the problem of compound fault diagnosis of wind power gearboxes， this paper proposed a feature extraction method for compound faults of wind power gearboxes based on improved variational modes， and used DEWOA-LSSVM fault diagnosis model for fault diagnosis. The variational modal method was used to decompose the collected signal into K modal components， the weighted ranking entropy was used to quantize it to obtain the feature vector， the feature vector was input to the least square support vector machine model optimized by the improved whale algorithm to complete the fault diagnosis， and the effectiveness of the method was verified through experiments.

Keywords： fault diagnosis;variational modal decomposition;weighted ranking entropy;whale algorithm;least squares support vector machine

齒輪箱是風力發(fā)電機組的核心部件，如果齒輪箱發(fā)生故障，容易造成內(nèi)部結(jié)構(gòu)的損害，人們需要停機進行維修，而且故障維修時間較長，工作量較大。這對風力發(fā)電系統(tǒng)工作的穩(wěn)定和安全性來說也是一種很大的威脅，嚴重時會造成事故和財產(chǎn)損失。

風電機組齒輪箱在實際應用過程中存在多種故障狀態(tài)，在實際工作中也可能會存在多種工況，而在實際工況中對風電齒輪箱的故障進行檢測時，風電齒輪箱噪聲和外界環(huán)境噪聲會影響齒輪箱的正常運行[1]。所以，風電齒輪箱的故障應該是非線性的，這些噪聲都會對風電齒輪箱的故障診斷產(chǎn)生誤判，因此人們需要對故障信號進行處理。目前，風電齒輪箱故障信號的處理主要集中在時頻域分析，主要有小波包分解、傅里葉變化、經(jīng)驗模態(tài)分解等時頻分析算法[2]。風電齒輪箱的故障信號存在多個頻域，人們可以將故障信號分解為多個不同頻域的信號進行分析。通過加權(quán)排序熵對故障信號進行處理，使得相同排序模式可以表示其特征的差異性，保留多重數(shù)據(jù)的微弱特征，同時沒有增加整體計算方法的計算量[3]。所以，加權(quán)排序熵可以更好地解決多種信號特征下的狀態(tài)分析問題。因此，本文采用變分模態(tài)對復雜故障進行分解，通過加權(quán)排序熵對數(shù)據(jù)進行分析，最后結(jié)合改進鯨魚算法優(yōu)化的最小二乘支持向量機，獲得故障診斷指標，完成故障診斷。

1 基于VMMWPE的故障特征提取

1.1 變分模態(tài)算法

變分模態(tài)分解算法（VMD）對模態(tài)的定義進行了更新，主要將模態(tài)函數(shù)變?yōu)榫哂兄行念l率的信號[4]。具體的表達式如式（1）所示。

式中，[Ak（t）]為[uk（t）]的瞬時值;[wk（t）]為[uk（t）]的瞬時頻率，[dφkk/dt=wk（t）]。如果確定信號分量[uk（t）]中以中心頻率為核心的帶寬，保證信號帶寬最小，則模態(tài)函數(shù)即可代替信號。

VMD算法的具體流程為：一是參數(shù)初始化，對信號[u1k]和[ω1k]進行初始化;按照上述的計算步驟開始循環(huán)，并對參數(shù)進行更新;判斷是否達到迭代條件[kun+1k-unk22unk22<ε]，當未達到時，則轉(zhuǎn)到第二步驟，條件成立時則終止循環(huán)并輸出IMF分量。

1.2 加權(quán)排序熵的基本原理

排序熵（PE）可以處理多種數(shù)據(jù)在時間序列上的排序分析，具有優(yōu)異的處理速度，能夠分析多重數(shù)據(jù)的特征變化。因此，排序熵可以很好地分析各種故障信息組合，并從中提煉出各個特征的狀態(tài)指標值[5]。

對一組時間序列數(shù)據(jù)[x=xi，1，2，…，N]進行計算時，首先需要確定一組時間序列的長度[N]，并得到[m]重數(shù)據(jù)的狀態(tài)向量[X=xi，xi+1，…，xi+m-1τ]。然后，對多維向量[X]進行升序排列并記錄各種數(shù)據(jù)的時間值，將數(shù)值相同的數(shù)據(jù)按照時間序列順序進行排列。此時，狀態(tài)向量的排序模式共有[m]次，而不同的模式對應不同的映射向量，可以計算不同模式下的概率，如式（2）所示。

式中，C為不同的排序模式次數(shù)。

下面可以進一步定義排序熵，如式（3）所示。

為使熵值取值介于0～1，對式（3）進行歸一化計算，具體計算公式如式（4）所示。

式中，[PE]為一組時間序列數(shù)據(jù)下的熵值，熵值可以表示該組數(shù)據(jù)在時間序列上的波動性。

當[PE]=0時，數(shù)據(jù)較為平穩(wěn)，表示隨著時間變化，數(shù)據(jù)整體沒有上下變動性，即這組數(shù)據(jù)中的各種數(shù)據(jù)值在這一時間序列下具有唯一性，沒有間隔差異時間下的相同數(shù)據(jù)。當[PE]=1時，該組數(shù)據(jù)在時間序列上具有劇烈的波動性。

1.3 基于變分模態(tài)和加權(quán)排序熵的信號特征提取

對于風電機組齒輪箱中的復雜故障變化，首先通過VMD算法對故障信號進行分解，得到具有不同特征信號的IMF分量，然后利用加權(quán)排序熵（WPE）計算方法對不同的IMF分量進行分析，使得不同IMF分量的故障信息實現(xiàn)過量化，最后通過對比各種故障信息的指標，確定齒輪箱的具體故障狀態(tài)。

基于VMMWPE的齒輪箱故障分析步驟如下：一是通過傳感器對風電機組齒輪箱的振動信息進行采集;二是利用VMD算法對采集的振動信息進行分解，在得到的IMF分量中，選取含有故障信息的前m個IMF分量，為了區(qū)分不同的IMF分量，設定IMF分量的尺度[s]，其中第[i]個IMF分量對應的尺度為[s-i+1];三是基于VMMWPE的非線性故障特征進行量化分析。按照WPE方法，對步驟二中選取的IMF分量進行計算，通過觀察熵值-尺度曲線，人們就可以了解不同故障信號的特征信息。

2 DEWOA-LSSVM的故障診斷模型

LSSVM是在SVM算法的基礎上針對SVM在優(yōu)化求解上存在的問題做了改進，核心算法原理是相同的。在構(gòu)建目標函數(shù)的約束條件時，利用等式約束取代傳統(tǒng)SVM中的不等式約束，這種等式約束條件的變化能夠求解復雜的一次非線性規(guī)劃問題[6]。通過構(gòu)建等式約束條件求解線性方程組，不僅降低了算法求解的難度和條件的復雜度，還提高了在求解時算法的收斂精度、速度、運行效率[7]。最小二乘支持向量機中的核心工作之一就是確定回歸方程的核函數(shù)，針對最小二乘支持向量機的核函數(shù)和懲罰項的優(yōu)化，本文提出了差分進化鯨魚算法（DEWOA）。其遵循如下原理。

2.1 基于差分進化鯨魚算法的變異操作過程

針對每一個目標個體[Xi，Gi=1，2，…，N]，從當前代[G]中隨機選擇3個沒有相同點的個體進行變異，相關函數(shù)表達式為：

式中，[r1，r2，r3?1，2，…，N]，為隨機產(chǎn)生的整數(shù)，并且滿足[r1≠r≠r3≠i];[F]為縮放比例因子，取值范圍為[0，1]，這里取0.5。

2.2 基于差分進化鯨魚算法的交叉操作

為實現(xiàn)目標個體[Xi，Gi=1，2，…，N]的進化，通過交叉操作，隨機選擇可使試驗向量[Ui，G+1U1，i，G+1，U2，i，G+1，…，Un，i，G+1]中至少含有一個貢獻的目標個體[Xi，G]，通過貢獻的目標個體改變種群。具體公式為：

式中，[randj?0，1]為均勻分布的隨機數(shù);[rnbi?1，2，…，n]，為隨機整數(shù);[CR?0，1]，為交叉概率，這里取0.1。

2.3 適應度函數(shù)的變異交叉

完成交叉操作后，對比交叉操作后的個體與當前個體的值，哪個適應度函數(shù)小，哪個將進入一下代的種群中進行變異交叉。

基于改進變分模態(tài)和加權(quán)排序熵的最小二乘支持向量機的雙饋風電機組復合故障診斷流程如圖1所示。

3 試驗分析

3.1 基于VMMWPE的故障特征提取

本文在不同試驗條件下輸入不同的原始信號，具體原始信號的輸入如表1所示。其間分別在三種不同的齒輪箱障礙和正常狀態(tài)下進行試驗，其中，齒輪箱總共分成五個機械部件。筆者通過傳感器采集試驗數(shù)據(jù)，采集頻率設置為10 kHz。

為了使得加權(quán)排序熵可以匹配到適宜的數(shù)據(jù)長度和維度，本研究開展了數(shù)據(jù)匹配性試驗，設置的數(shù)據(jù)長度從10 000降低到500。數(shù)據(jù)長度會影響計算結(jié)果的穩(wěn)定性，因此本文取數(shù)據(jù)長度為2 048。接下來分析延遲信號時間[τ]對計算結(jié)果的影響，[τ]從1開始取值，當[τ]取1時，不同狀態(tài)的加權(quán)排序熵值可以得到良好的特征值，因此本文中[τ]取1。最后確定數(shù)據(jù)維度和加權(quán)排序熵值的關系，在維度3～7區(qū)間內(nèi)取值?？梢钥闯?，區(qū)分度最高的數(shù)據(jù)維度為5，故[m]值取5。

筆者通過試驗平臺提取齒輪箱四種故障狀態(tài)數(shù)據(jù)各60組，并利用VMMWPE處理原始信號來求取分量熵值，熵值分布情況如圖2所示。

3.2 基于DEWOA-LSSVM的故障診斷

將所得到的60組熵值作為數(shù)據(jù)樣本，以40組作為訓練樣本，其余20組作為測試樣本?；趫D1所示的改進變分模態(tài)和加權(quán)排序熵的支持向量機的雙饋風電機組復合故障診斷流程，利用訓練樣本得到最優(yōu)的LSSVM模型，再通過測試樣本的輸入驗證，得到診斷結(jié)果為100%，如圖3所示。

4 結(jié)論

本文利用變分模態(tài)和加權(quán)排序熵對風電齒輪箱復合特征進行提取分析，得到歸一化熵值。其間通過最小二乘支持向量機對所得到的熵值進行分類，完成故障診斷，通過差分進化鯨魚算法來優(yōu)化最小二乘支持向量機的核函數(shù)和懲罰項，通過試驗驗證了基于VMMWPE-DEWOA-LSSVM算法的風電齒輪箱復合故障診斷的有效性。

參考文獻：

[1]王世棟，蘇欣，劉義虎，等.齒輪箱振動噪聲仿真分析與結(jié)構(gòu)優(yōu)化[J].熱能動力工程，2019（1）：92-97.

[2]丁文潔，趙武云，張征，等.變分模態(tài)分解方法濾波特性及其在齒輪箱故障檢測中的應用[J].機電工程，2019（8）：32-35.

[3]任國春，趙永東，馮輔周，等.基于自適應EEMD樣本熵的行星齒輪箱特征提取方法[J].裝甲兵工程學院學報，2017（6）：49-55.

[4]胡蔦慶，陳徽鵬，程哲，等.基于經(jīng)驗模態(tài)分解和深度卷積神經(jīng)網(wǎng)絡的行星齒輪箱故障診斷方法[J].機械工程學報，2019（7）：9-18.

[5]Sun C，Wang Y，Yan Z.Vibration diagnostic method based on improved intrinsic time-scale decomposition and energy operator demodulation for fault diagnosis of helicopter planetary gearboxes[C]//2016 Prognostics and System Health Management Conference.2017.

[6]徐可，陳宗海，張陳斌，等.基于經(jīng)驗模態(tài)分解和支持向量機的滾動軸承故障診斷[J].控制理論與應用，2019（6）：915-922.

[7]侯志華.基于最小二乘復頻域法的某齒輪箱測試研究[J].汽車維修，2017（2）：8-11.

[8]徐可，陳宗海，張陳斌，等.基于經(jīng)驗模態(tài)分解和支持向量機的滾動軸承故障診斷[J].控制理論與應用，2019（6）：915-922.