變分模態(tài)分解在轉(zhuǎn)子故障診斷中的應(yīng)用*

向 玲,張力佳

(華北電力大學(xué)機(jī)械工程系 保定，071003)

變分模態(tài)分解在轉(zhuǎn)子故障診斷中的應(yīng)用*

向 玲,張力佳

(華北電力大學(xué)機(jī)械工程系 保定，071003)

針對(duì)轉(zhuǎn)子故障診斷問(wèn)題，提出一種基于變分模態(tài)分解(variational mode decomposition，簡(jiǎn)稱VMD)的信號(hào)處理方法。該方法在獲取分解分量的過(guò)程中通過(guò)迭代搜尋變分模型最優(yōu)解來(lái)確定每個(gè)分量的頻率中心及帶寬，從而能夠自適應(yīng)地實(shí)現(xiàn)信號(hào)的頻域剖分及各分量的有效分離，對(duì)各單分量信號(hào)進(jìn)行希爾伯特變換，即可得到瞬時(shí)的頻率和幅值信息。對(duì)仿真信號(hào)和典型轉(zhuǎn)子故障信號(hào)進(jìn)行VMD方法和經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition,簡(jiǎn)稱EMD)方法的分析比較，以驗(yàn)證所提方法的有效性。仿真信號(hào)的分解結(jié)果表明，變分模態(tài)能夠準(zhǔn)確分離出信號(hào)中的固有模態(tài)分量且不存在模態(tài)混疊；轉(zhuǎn)子故障實(shí)驗(yàn)信號(hào)的分析結(jié)果表明，所提方法能夠有效提取出明顯的故障特征，從而準(zhǔn)確診斷出轉(zhuǎn)子存在的故障。

變分模態(tài)分解； 轉(zhuǎn)子故障； 故障特征； 時(shí)頻分析

引 言

轉(zhuǎn)子故障診斷的方法有多種，將時(shí)域和頻域相結(jié)合反映信號(hào)全局特征的時(shí)頻分析是一種重要的非線性、非平穩(wěn)信號(hào)的處理方法。常用的時(shí)頻分析方法如希爾伯特黃變換[1-2](Hilbert-Huang transform,簡(jiǎn)稱HHT)能夠在時(shí)頻域中清晰展示出信號(hào)中的故障特征，從而有效進(jìn)行故障診斷。HHT由EMD和希爾伯特變換(Hilbert transform,簡(jiǎn)稱HT)組成，即對(duì)一個(gè)采集到的信號(hào)采用EMD進(jìn)行分解，得到有限固有模態(tài)函數(shù)(intrinsic mode function,簡(jiǎn)稱IMF)，然后使用HT對(duì)各分量進(jìn)行解調(diào)，得到原始信號(hào)的時(shí)頻信息。HHT時(shí)頻分析方法在近幾年的機(jī)械故障診斷中得到了有效應(yīng)用[3-7]，但是EMD分解過(guò)程中的模態(tài)混疊現(xiàn)象影響了后續(xù)的時(shí)頻分析[8-9]。VMD[10]是一種新的自適應(yīng)信號(hào)預(yù)處理方法，分解過(guò)程實(shí)質(zhì)是變分問(wèn)題的迭代求解，可分為變分模型的建立和求解，是建立在Wiener濾波、Hilbert變換和外差解調(diào)等理論基礎(chǔ)上的自適應(yīng)信號(hào)分解方法[11-13]，對(duì)其分解得到的單分量進(jìn)行Hilbert變換即可得到時(shí)頻圖。

筆者提出基于變分模態(tài)分解的轉(zhuǎn)子故障信號(hào)時(shí)頻分析方法，對(duì)仿真信號(hào)和轉(zhuǎn)子故障實(shí)驗(yàn)信號(hào)進(jìn)行了基于EMD方法與VMD方法的對(duì)比。結(jié)果表明，VMD方法能更準(zhǔn)確地分離出轉(zhuǎn)子故障信號(hào)分量，且故障特征值明顯，可有效診斷出轉(zhuǎn)子的故障。

1 變分模態(tài)分解原理

1.1 變分模型的建立

假設(shè)輸入信號(hào)x(t)是由有限個(gè)中心頻率不同、帶寬有限的分量組成，將信號(hào)分解問(wèn)題轉(zhuǎn)換到變分模型中進(jìn)行分解。在各分量之和等于原始輸入信號(hào)x(t)約束下，尋求各分量的聚集帶寬之和最小。構(gòu)造步驟如下[10]。

3) 計(jì)算解調(diào)信號(hào)梯度的平方L2范數(shù)，估計(jì)出模態(tài)信號(hào)帶寬。受約束的最優(yōu)化變分模型為

(1)

1.2 變分模型的求解

使用二次懲罰因子α和拉格朗日乘法算子λ(t)，將約束性問(wèn)題變?yōu)榉羌s束性問(wèn)題。二次懲罰因子可以保證含高斯噪聲的情況下信號(hào)的重構(gòu)精度，拉格朗日乘子可以保證約束條件的嚴(yán)格性。擴(kuò)展的拉格朗日函數(shù)表達(dá)式為

(2)

(3)

利用Parseval/Plancherel傅里葉等距離變換，將式(3)變換到頻域，得到各模態(tài)的頻域更新

（3）另取0.1 mol/L的氫氧化鈉溶液2mL，滴加2滴1%的酚酞溶液，用水稀釋至20mL，將0.1 mol/L的鹽酸改為0.01 mol/L的稀鹽酸重復(fù)上述實(shí)驗(yàn)操作，直至溶液的pH為2.00為止，用蒸餾水沖洗pH電極。

(4)

根據(jù)同樣計(jì)算過(guò)程，將中心頻率取值問(wèn)題轉(zhuǎn)換到頻域，得到中心頻率的更新

(5)

1.3 VMD算法具體過(guò)程

VMD算法具體步驟[13]如下：

2) 根據(jù)式(4)，(5)，更新uk和ωk；

3) 更新λ為

2 時(shí)頻譜與邊際譜

(6)

構(gòu)造解析信號(hào)，得到幅值和相位函數(shù)為

瞬時(shí)頻率為

(9)

(10)

其中：Re表示取實(shí)部并忽略殘余項(xiàng)；ak(t)為解析信號(hào)的瞬時(shí)幅值；ωk(t)為解析信號(hào)的瞬時(shí)頻率。

式(10)描述了信號(hào)中每個(gè)時(shí)刻包含的頻率成分，據(jù)此可作出Hilbert時(shí)頻譜[14]。

進(jìn)一步定義邊際譜為

(11)

邊際譜的幅值能準(zhǔn)確反映原始信號(hào)中某頻率分量是否真實(shí)存在，表示信號(hào)在某頻率下的時(shí)間段內(nèi)的幅值之和[15]。

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

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

利用EMD方法對(duì)仿真信號(hào)進(jìn)行分解，得到的分量時(shí)域圖如圖2所示,從上到下依次為IMF1，IMF2，IMF3，IMF4，IMF5。EMD分解受仿真信號(hào)中隨機(jī)噪聲信號(hào)的影響，分解得到的信號(hào)分量存在明顯的模態(tài)混疊現(xiàn)象。

圖2 仿真信號(hào)的EMD分解結(jié)果Fig.2 EMD result generated from the simulation signal

VMD分解與EMD分解的最大不同在于VMD需要提前設(shè)置所需處理信號(hào)的分解個(gè)數(shù)k。若k值偏大，會(huì)導(dǎo)致分解結(jié)果出現(xiàn)過(guò)分解；若k值偏小，則會(huì)出現(xiàn)欠分解，因此合理設(shè)置k值是使用VMD的關(guān)鍵。筆者采用中心頻率相近原則[16]來(lái)選擇k值，同時(shí)，懲罰因子α采用默認(rèn)值2 000。不同k值時(shí)對(duì)仿真信號(hào)進(jìn)行分解，并計(jì)算各模態(tài)分量的中心頻率，如表1所示。由表1可知，當(dāng)k=4時(shí)出現(xiàn)了頻率中心相近的模態(tài)分量，即出現(xiàn)了過(guò)分解。因此，使用VMD對(duì)仿真信號(hào)進(jìn)行3層分解，得到的各分量時(shí)域圖如圖3所示，從上到下依次為F1，F(xiàn)2，F(xiàn)3。圖3中，F(xiàn)1，F(xiàn)2分量分別對(duì)應(yīng)仿真信號(hào)中的正弦信號(hào)與調(diào)幅調(diào)頻信號(hào)，信號(hào)分量時(shí)域圖規(guī)律清晰，不存在模態(tài)混疊現(xiàn)象；F3分量為仿真信號(hào)中的隨機(jī)噪聲信號(hào)。通過(guò)對(duì)比EMD，VMD兩種方法對(duì)仿真信號(hào)的分解，表明VMD能夠準(zhǔn)確將仿真信號(hào)中的不同模態(tài)分量分離。

表1 不同k值各模態(tài)分量的中心頻率Tab.1 Center frequency corresponding to different k

圖3 仿真信號(hào)的VMD分解結(jié)果Fig.3 VMD result generated from the simulation signal

4 實(shí)驗(yàn)故障信號(hào)分析

使用Bebntly RK-4轉(zhuǎn)子實(shí)驗(yàn)臺(tái)模擬轉(zhuǎn)子不平衡、碰摩、油膜渦動(dòng)等典型故障。故障信號(hào)的記錄設(shè)備采用美國(guó)Iotech公司的ZonicBook/618E。

4.1 轉(zhuǎn)子不平衡故障分析

加載不平衡質(zhì)量1 g在轉(zhuǎn)子圓盤45°角處進(jìn)行轉(zhuǎn)子不平衡實(shí)驗(yàn)。調(diào)整轉(zhuǎn)子的轉(zhuǎn)速為2 480 r/min，模擬轉(zhuǎn)子不平衡故障。采樣頻率為1 280 Hz，選取記錄數(shù)據(jù)中的1 024點(diǎn)進(jìn)行分析。

圖4為轉(zhuǎn)子不平衡故障信號(hào)的時(shí)域圖和頻譜圖。由圖4(a)可知，轉(zhuǎn)子不平衡故障信號(hào)的時(shí)域波形呈現(xiàn)正弦波形。在頻譜圖4(b)中，41.33 Hz基頻幅值明顯突出，有二倍頻等高頻分量出現(xiàn),但幅值較小，不能通過(guò)頻譜圖判斷出是否有轉(zhuǎn)子故障發(fā)生。

圖4 轉(zhuǎn)子不平衡故障信號(hào)的時(shí)域圖和頻譜圖Fig.4 Time domain and spectral diagram of rotor imbalance signal

圖5 轉(zhuǎn)子不平衡信號(hào)的時(shí)頻圖Fig.5 Time-frequency spectrum of rotor imbalance signal

圖6 轉(zhuǎn)子不平衡信號(hào)的邊際譜Fig.6 Marginal spectrum of rotor imbalance signal

圖5(a)為基于EMD方法的轉(zhuǎn)子不平衡故障信號(hào)時(shí)頻分析圖。圖中可辨識(shí)出轉(zhuǎn)子不平衡故障特征頻率為41.33 Hz，但故障特征頻率的調(diào)頻現(xiàn)象明顯，且在0.6～0.8 s存在模態(tài)混疊。根據(jù)中心頻率相近原則，轉(zhuǎn)子不平衡故障信號(hào)選擇k值為4進(jìn)行VMD分解分析，得到的時(shí)頻圖如圖5(b)所示。圖中不平衡故障特征頻率41.33 Hz清晰可見(jiàn)，且不存在模態(tài)混疊現(xiàn)象。圖6為兩種方法得到的邊際譜。結(jié)果表明，VMD能夠準(zhǔn)確分離出原始信號(hào)中的不平衡故障信號(hào)分量，且故障特征的幅值明顯突出。在轉(zhuǎn)子不平衡故障信號(hào)的時(shí)頻分析中，基于VMD方法的時(shí)頻圖中故障特征頻率清晰明顯，可有效診斷出轉(zhuǎn)子不平衡故障。

4.2 轉(zhuǎn)子碰摩故障分析

轉(zhuǎn)子局部碰摩實(shí)驗(yàn)時(shí)轉(zhuǎn)子轉(zhuǎn)速為1 260 r/min，用碰摩針與轉(zhuǎn)軸進(jìn)行局部碰摩故障實(shí)驗(yàn)。采樣頻率設(shè)置為1 280 Hz，選取記錄數(shù)據(jù)中的1 024點(diǎn)進(jìn)行分析。圖7為轉(zhuǎn)子局部碰摩故障信號(hào)的時(shí)域圖和頻譜圖。由圖7(a)可看出，轉(zhuǎn)子局部碰摩時(shí)域圖為規(guī)則正弦，在7(b)中，21 Hz基頻分量及2倍頻42 Hz幅值明顯，其他高倍頻分量幅值較小且不易辨識(shí)，簡(jiǎn)單的頻譜分析不能有效判斷出轉(zhuǎn)子故障。

圖7 轉(zhuǎn)子局部碰摩信號(hào)的時(shí)域和頻譜圖Fig.7 Time domain and spectral diagram of local rotor rubbing signal

圖8 轉(zhuǎn)子局部碰摩信號(hào)的時(shí)頻圖Fig.8 Time-frequency spectrum of rotor rubbing fault signal

圖9 轉(zhuǎn)子碰摩故障信號(hào)的邊際譜Fig.9 Marginal spectrum of local rotor rubbing fault signal

圖8(a)為基于EMD方法的轉(zhuǎn)子碰摩故障信號(hào)時(shí)頻分析圖，可辨識(shí)出轉(zhuǎn)子碰摩故障特征頻率21 Hz及不穩(wěn)定的倍頻成分。根據(jù)中心頻率相近原則，轉(zhuǎn)子碰摩故障信號(hào)選擇k值為5進(jìn)行VMD分解分析，得到的時(shí)頻圖如圖8(b)所示。圖中碰摩故障特征頻率21 Hz及2倍頻42 Hz清晰明顯，不存在模態(tài)的混疊。進(jìn)一步對(duì)比兩種方法得到的邊際譜圖9(a),(b)可知，VMD準(zhǔn)確分離出原始信號(hào)中的轉(zhuǎn)子碰摩故障信號(hào)分量，且故障特征頻率的幅值明顯突出?；赩MD分解的時(shí)頻分析方法提取出的碰摩故障特征明顯，可有效判斷出轉(zhuǎn)子局部碰摩故障。

4.3 轉(zhuǎn)子油膜渦動(dòng)故障分析

轉(zhuǎn)子產(chǎn)生油膜渦動(dòng)故障時(shí)，理論上會(huì)產(chǎn)生頻率值為轉(zhuǎn)子轉(zhuǎn)頻0.43～0.48倍的故障特征頻率。通過(guò)轉(zhuǎn)子實(shí)驗(yàn)平臺(tái)的調(diào)速電機(jī)控制器，使轉(zhuǎn)子轉(zhuǎn)速達(dá)到2 800 r/min，轉(zhuǎn)子產(chǎn)生動(dòng)力失穩(wěn)從而模擬油膜渦動(dòng)故障。采樣頻率設(shè)為1 280 Hz，選取數(shù)據(jù)中的1 024點(diǎn)進(jìn)行分析。

圖10為轉(zhuǎn)子油膜渦動(dòng)故障信號(hào)的時(shí)域圖和頻譜圖。由圖10(a)可知，時(shí)域圖為規(guī)則信號(hào)的波動(dòng)，在頻譜圖10(b)中，幅值最大的頻率分量為半頻24 Hz，基頻46.67 Hz幅值也有明顯突出，同時(shí)還出現(xiàn)了幅值不明顯的高倍頻分量，依靠單一頻譜圖分析無(wú)法準(zhǔn)確判斷出轉(zhuǎn)子是否發(fā)生油膜渦動(dòng)故障。

圖10 轉(zhuǎn)子油膜渦動(dòng)故障信號(hào)的時(shí)域圖和頻譜圖Fig.10 Time domain and spectral diagram of rotor oil whirl signal

圖11 轉(zhuǎn)子油膜渦動(dòng)信號(hào)的時(shí)頻圖Fig.11 Time-frequency spectrum of rotor oil whirl signal

圖11(a)為基于EMD方法得到的轉(zhuǎn)子油膜渦動(dòng)故障信號(hào)時(shí)頻分析圖。圖中可以辨識(shí)出轉(zhuǎn)子油膜渦動(dòng)的故障特征頻率24 Hz，但是0.4～0.6 s有明顯的模態(tài)混疊。根據(jù)中心頻率相近原則，轉(zhuǎn)子油膜渦動(dòng)故障信號(hào)選擇k值為4進(jìn)行VMD分解分析，得到的時(shí)頻圖如圖11(b)所示。圖中可以有效辨識(shí)出轉(zhuǎn)子基頻46.67 Hz與油膜渦動(dòng)故障的特征頻率24 Hz，且不存在模態(tài)混疊。12(a)，(b)為通過(guò)兩種方法得到的邊際譜圖。結(jié)果表明，VMD可準(zhǔn)確分離出原始信號(hào)中的故障分量信號(hào)，且故障特征的幅值明顯突出，不存在模態(tài)混疊?；赩MD方法的時(shí)頻分析得到的故障特征清晰明顯，油膜渦動(dòng)故障特征頻率提取效果優(yōu)于基于EMD方法的時(shí)頻分析，故障診斷效果可得到提高。

圖12 轉(zhuǎn)子不平衡信號(hào)的邊際譜Fig.12 Marginal spectrum of rotor oil whirl signal

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

提出一種基于變分模態(tài)分解的轉(zhuǎn)子故障信號(hào)時(shí)頻分析方法，針對(duì)仿真信號(hào)與轉(zhuǎn)子實(shí)驗(yàn)故障信號(hào)對(duì)比分析了EMD和VMD方法的有效性。仿真信號(hào)的分解結(jié)果表明，EMD分解得到的分量信號(hào)存在模態(tài)混疊現(xiàn)象，VMD方法分解得到的分量信號(hào)與仿真信號(hào)中存在的分量信號(hào)一致。轉(zhuǎn)子典型故障信號(hào)的分析結(jié)果表明，基于VMD的時(shí)頻分析方法能夠有效避免出現(xiàn)模態(tài)混疊，提取出的故障特征清晰明顯，可準(zhǔn)確診斷出轉(zhuǎn)子存在的故障。這為轉(zhuǎn)子故障的有效診斷提供了一種新的時(shí)頻分析方法。

[1] Huang N E,Shen Z,Long S R,et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceedings of the Royal Society A，1998,454:903-995.

[2] Huang N E,Shen Z,Long S R,et al.A new view of nonlinear water waves:the Hilbert spectrum[J].Annual Review of Fluid Mechanics,1999,31:417-457.

[3] Yang Zhensheng,Yu Zhonghua,Chao Xie,et al.Application of Hilbert-Huang transform to acoustic emission signal for burn feature extraction in surface grinding process[J].Measurement,2014,47:14-21.

[4] 鄭近德，程軍圣.改進(jìn)的希爾伯特-黃變換及其在滾動(dòng)軸承故障診斷中的應(yīng)用[J].機(jī)械工程學(xué)報(bào),2015,51(1):138-145.

Zheng Jinde,Cheng Junsheng.Improved Hilbert-Huang transform and its applications to rolling bearing fault diagnosis[J].Journal of Mechanical Engineering,2015,51(1):138-145.(in Chinese)

[5] 馬文朋，張俊紅，馬梁，等.改進(jìn)的經(jīng)驗(yàn)?zāi)Ｊ椒纸庠跈C(jī)械故障診斷中的應(yīng)用[J].振動(dòng)、測(cè)試與診斷,2015,35(4):637-644.

Ma Wenpeng,Zhang Junhong,Ma Liang,et al.Applications of improved empirical mode decomposition in machinery diagnosis[J].Journal of Vibration,Measurement &Diagnosis,2015,35(4):637-644.(in Chinese)

[6] 熊炘，楊世錫，甘春標(biāo).運(yùn)用BF-HHT的轉(zhuǎn)子徑向碰摩故障特征提取方法[J].振動(dòng)、測(cè)試與診斷,2012,32(2):192-196.

Xiong Xin,Yang Shixi,Gan Chunbiao.Fault feature extraction method for rotor system under radial impact based on BE-HHT[J].Journal of Vibration,Measurement &Diagnosis,2012,32(2):192-196.(in Chinese)

[7] 孟宗，李?yuàn)檴?基于小波改進(jìn)閾值去噪和HHT的滾動(dòng)軸承故障診斷[J].振動(dòng)與沖擊,2013,32(14):204-208.

Meng Zong,Li Shanshan.Rolling bearing fault diagnosis based on improved wavelet threshold de-noising method and HHT[J].Journal of Vibration and Shock,2013,32(14):204-208.(in Chinese)

[8] 曹瑩，段玉波，劉繼承.Hilbert-Huang變換中的模態(tài)混疊問(wèn)題[J].振動(dòng)、測(cè)試與診斷,2016,36(3):518-523.

Cao Ying,Duan Yubo,Liu Jicheng.Research and application of mode-mixing in Hilbert-Huang transform[J].Journal of Vibration,Measurement &Diagnosis,2016,36(3):518-523.(in Chinese)

[9] 向玲，鄢小安.汽輪機(jī)轉(zhuǎn)子故障診斷中LMD法和EMD法的性能對(duì)比研究[J].動(dòng)力工程學(xué)報(bào),2014,34(12):945-951.

Xiang Ling,Yan Xiaoan.Performance contrast between LMD and EMD in fault diagnosis of turbine rotors[J].Journal of Chinese Society of Power Engineering,2014,34(12):945-951.(in Chinese)

[10]Dragomiretskiy K,Zosso D． Variational mode decomposition[J].Transactions on Signal Processing,2014,62(3):531-544．

[11]Aneesh C,Sachin K,Hisham P M,et al.Performance comparison of variational mode decomposition over empirical wavelet transform for the classification of power quality disturbances using support vector machine[J].Procedia Computer Science,2015(46):372-380.

[12]Salim L.Comparative study of signal denoising by wavelet threshold in empirical and variational mode decomposition domains[J].Healthcare Technology Letters,2014,1(3):104-109.

[13]劉長(zhǎng)良，武英杰，甄成剛.基于變分模態(tài)分解和模糊C聚類均值聚類的滾動(dòng)軸承故障診斷[J].中國(guó)電機(jī)工程學(xué)報(bào),2015,35(13):3358-3365.

Liu Changliang,Wu Yingjie,Zhen Chenggang.Rolling bearing fault diagnosis based on variational mode decomposition and fuzzy C means clustering[J].Proceedings of the CSEE,2015,35(13):3358-3365.(in Chinese)

[14]張帥，楊勇，韓清凱，等.基于HHT的轉(zhuǎn)子系統(tǒng)定點(diǎn)碰摩實(shí)驗(yàn)研究[J].振動(dòng)與沖擊,2010,29(7):121-125.

Zhang Shuai,Yang Yong,Han Qingkai,et al.Experimental study on a rotor system with rub-impact at fixed limiter based on HHT[J].Journal of Vibration and Shock,2010,29(7):121-125.(in Chinese)

[15]馬百雪，潘宏俠，楊素梅.基于改進(jìn)的HHT邊際譜齒輪箱故障診斷[J].煤礦機(jī)械,2013,34(10):251-254.

Ma Baixue,Pan Hongxia,Yang Sumei.Gearbox fault diagnosis based on improved Hilbert-Huang transform marginal spectrum[J].Coal Mine Machinery,2013,34(10):251-254.(in Chinese)

[16]馬增強(qiáng)，李亞超，劉政，等.基于變分模態(tài)分解和Teager能量算子的滾動(dòng)軸承故障特征提取[J].振動(dòng)與沖擊,2016,35(13):134-139.

Ma Zengqiang,Li Yachao,Liu Zheng,et al.Rolling bearings′ fault feature extraction based on variational mode decomposition and Teager energy operator[J].Journal of Vibration and Shock,2016,35(13):134-139.(in Chinese)

10.16450/j.cnki.issn.1004-6801.2017.04.025

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

2016-06-24；

2016-08-08

TH165.3； TK288

向玲，女，1971年4月生，博士、教授。主要研究方向?yàn)榉蔷€性動(dòng)力學(xué)和故障診斷。曾發(fā)表《汽輪發(fā)電機(jī)組軸系扭振的時(shí)頻特征分析》(《動(dòng)力工程學(xué)報(bào)》2011年第31卷第9期)等論文。E-mail:xiangling5042@126.com