基于微分算子的局部特征尺度分解方法及其應(yīng)用

吳占濤　程軍圣　李寶慶　鄭近德

1. 湖南大學(xué)汽車車身先進(jìn)設(shè)計制造國家重點實驗室，長沙，4100822.安徽工業(yè)大學(xué)，馬鞍山，243032

基于微分算子的局部特征尺度分解方法及其應(yīng)用

吳占濤1程軍圣1李寶慶1鄭近德2

1. 湖南大學(xué)汽車車身先進(jìn)設(shè)計制造國家重點實驗室，長沙，4100822.安徽工業(yè)大學(xué)，馬鞍山，243032

針對局部特征尺度分解(LCD)方法在目標(biāo)信號所含頻率分量較接近時分解能力降低，易出現(xiàn)模態(tài)混淆現(xiàn)象，從而導(dǎo)致內(nèi)稟尺度分量失去原有物理意義的問題，提出了基于微分算子的局部特征尺度分解(DOLCD)。DOLCD在對目標(biāo)信號進(jìn)行分解前，先將目標(biāo)信號進(jìn)行一階微分，則在頻率比一定的情況下，可提升分解能力及抑制模態(tài)混淆能力。研究了DOLCD方法的原理，通過仿真信號模型將DOLCD與LCD的分解能力進(jìn)行對比分析，結(jié)果表明，DOLCD方法在提高分解能力，抑制模態(tài)混淆等方面具有一定的優(yōu)越性，并將DOLCD方法應(yīng)用于轉(zhuǎn)子不對中故障的診斷，結(jié)果表明該方法有效。

局部特征尺度分解；微分算子；故障診斷；內(nèi)稟尺度分量；轉(zhuǎn)子

0　引言

經(jīng)驗?zāi)B(tài)分解(empirical mode decomposition,EMD)是一種自適應(yīng)的非平穩(wěn)信號處理方法[1-2]，EMD方法自提出后在很多領(lǐng)域得到了廣泛的應(yīng)用[3]。然而，EMD方法存在包絡(luò)過沖和欠包絡(luò)、迭代計算量大、端點效應(yīng)以及模態(tài)混淆等問題[4]。作為對EMD方法的改進(jìn)，程軍圣等[5]提出了另一種自適應(yīng)時頻分析方法——局部特征尺度分解(local characteristic-scale decomposition, LCD)。LCD能夠自適應(yīng)地將一個復(fù)雜信號分解為若干個瞬時頻率具有物理意義的內(nèi)稟尺度分量(intrinsic scale component，ISC)之和，從而得到目標(biāo)信號完整的時頻分布。與EMD相比，LCD在減小擬合誤差、提高計算速度和抑制模態(tài)混淆等方面有一定優(yōu)勢，已被應(yīng)用于信號分析和機(jī)械故障診斷等領(lǐng)域，取得了較好的效果[6-7]。LCD方法理論基礎(chǔ)尚不完善，實際應(yīng)用中難免出現(xiàn)模態(tài)混淆現(xiàn)象，即同一尺度或頻率的信號被分解到多個不同的ISC分量當(dāng)中，或者一個ISC分量當(dāng)中出現(xiàn)了尺度或頻率差異較大的信號，降低了LCD方法的分解能力，導(dǎo)致分解得到的ISC分量失去原有的物理意義。

為解決LCD的模態(tài)混淆問題，鄭近德等[8]基于向目標(biāo)信號中成對地添加符號相反的白噪聲信號，對加噪信號依據(jù)頻率高低逐層進(jìn)行集成平均分解的思路，提出了部分集成局部特征分解(partly ensemble local characteristic-scale decomposition，PELCD)方法。PELCD對LCD分解中的模態(tài)混淆有很好的抑制作用，但由于PELCD不能自適應(yīng)地選擇添加白噪聲的幅值和數(shù)目，若參數(shù)選擇不合適，可能會產(chǎn)生偽分量。本文從另一種思路，解決LCD模態(tài)混淆問題。

改變目標(biāo)信號中低頻分量與高頻分量的幅值比，可提升LCD方法的分解能力，進(jìn)而有效抑制模態(tài)混淆的產(chǎn)生。借鑒劉海波等[9]在對EMD改進(jìn)時提出的利用微分算子增強(qiáng)EMD分解能力的思路，在對目標(biāo)信號進(jìn)行LCD分解前，先將目標(biāo)信號進(jìn)行一階微分，則在頻率比一定的情況下，可降低低頻分量信號與高頻分量信號的幅值比，提升LCD方法的分解能力及抑制模態(tài)混淆能力?；诖耍P者提出了基于微分算子的局部特征尺度分解(differential operator based local characteristic-scale decomposition, DOLCD)。

本文研究了DOLCD方法的原理，通過仿真信號將DOLCD與LCD的分解能力進(jìn)行對比分析，結(jié)果表明，DOLCD方法分解能力范圍要明顯大于LCD方法的分解能力范圍，增強(qiáng)了抑制模態(tài)混淆能力。采用DOLCD方法對具有不對中故障的轉(zhuǎn)子振動位移信號進(jìn)行了分析，結(jié)果表明DOLCD能夠有效地將高頻不對中成分與轉(zhuǎn)頻成分等進(jìn)行分離，從而實現(xiàn)轉(zhuǎn)子故障診斷。

1　DOLCD方法

楊宇等[10]對LCD方法的分解能力從理論層面上進(jìn)行了深入研究。研究結(jié)果表明LCD方法的分解能力主要由各分量信號的頻率比、幅值比決定，在頻率比一定的情況下，降低低頻分量信號與高頻分量信號幅值比，可提高LCD方法的分解能力，抑制模態(tài)混淆；而各分量信號的初始相位以及單個分量迭代次數(shù)對LCD方法的分解能力影響均不大。為提高LCD方法的分解能力，抑制模態(tài)混淆，本文提出了基于微分算子的局部特征尺度分解方法。對目標(biāo)信號x(t)(t>0)，DOLCD方法的主要分解步驟如下：

(3)對y1(t)進(jìn)行積分，得到z1(t)。

(4)對z1(t)進(jìn)行LCD分解，得到目標(biāo)信號x(t)的第一個有效ISC分量I1(t)與剩余信號r1(t)。

(5)將分量I1(t)與目標(biāo)信號x(t)分離，得到剩余信號s1(t)，即

s1(t)=x(t)-I1(t)

(1)

再將s1(t)視為目標(biāo)信號，重復(fù)上述步驟(1)～(5)，直到sn(t)為單調(diào)或常函數(shù)，依次得到ISC分量I1(t)，I2(t)，…，In(t)和趨勢項r1(t)，r2(t)，…，rn(t)。令

R(t)=r1(t)+r2(t)+…+rn(t)

(2)

(6)原目標(biāo)信號x(t)被分解為n個ISC分量和一個趨勢項R(t)之和，即

(3)

由于信號的初始相位對LCD的分解能力影響不大，為了討論方便、不失一般性，設(shè)輸入目標(biāo)信號x(t)由兩種模態(tài)構(gòu)成，即

x(t)=x1(t)+x2(t)=A1sin(2πf1t)+A2sin(2πf2t)

(4)

其中，f1>f2。低頻分量信號x2(t)與高頻分量信號x1(t)的頻率比F為

F=f2/f1

(5)

低頻分量信號x2(t)與高頻分量信號x1(t)的幅值比a為

a=A2/A1

(6)

對式(4)兩邊進(jìn)行一階微分得

2πf2A2cos(2πf2t)

(7)

b=(f2/f1)(A2/A1)

(8)

由于f1>f2，則b

2　DOLCD分解能力

2.1分解能力評估模型與評價指標(biāo)

為了評價DOLCD方法的分解能力，需確定一個分解能力的評價指標(biāo)。參考文獻(xiàn)[10]定義式(9)為分解能力評價指標(biāo)(evaluationindexofdecompositioncapacity)，記為E：

(9)

式中，xi(t)為單分量信號的真實值；Ii(t)為分解得到的單分量信號；T為信號長度。

為了保證分解得到的分量全部準(zhǔn)確，分解能力評價指標(biāo)由下式確定：

E=max[Ei]

(10)

由于式(10)可能出現(xiàn)E>1的情況，為了便于研究，對E>1的情況，則令E=1，即E∈[0,1]。E值越大說明分解得到的Ii(t)越偏離真實分量xi(t)，分解能力越弱；E值越小說明分解得到的分量Ii(t)越接近于真實分量xi(t)，分解能力越強(qiáng)。文獻(xiàn)[10-11]研究表明，當(dāng)E≤0.05時，可認(rèn)為混合信號被準(zhǔn)確分解。

2.2分解能力對比

為了不失一般性，仍采用式(4)所示的信號模型為分析對象，并設(shè)A1=4，f1=80，頻率比F∈[0.1,1]，幅值比a∈[0.1,∞]，信號采樣頻率fs=8192Hz。則在不同的幅值比a、頻率比F條件下得到DOLCD方法的分解能力評價指標(biāo)E值如圖1所示。

圖1　幅值比a和頻率比F對E值的影響

由圖1可以看出，當(dāng)a大致滿足a<1.1時，E值曲面形成一個坡面，該坡面在0.1≤F<0.7的部分較為平坦，此時E≤0.05；隨著F增大，坡面逐漸變陡，在F接近1的部分坡面又變得相對平坦，但E值接近1。當(dāng)a>1.1，F(xiàn)≤0.37時，E值曲面也較平坦，此時E≤0.05。當(dāng)a>1.1，F(xiàn)>0.37時，隨著F增大，會出現(xiàn)非常陡峭的曲面。

當(dāng)E≤0.05時，表示DOLCD方法具有較好的分解效果，此時可認(rèn)為兩個分量信號能夠被有效分離。當(dāng)0.05

圖2　DOLCD方法分解能力

由圖2可以看出，E=1和E=0.05兩條等高線將整個a-F平面劃分為3個區(qū)域：Ⅰ區(qū)為E=1區(qū)域，表示DOLCD方法完全不能將兩個單分量分離；Ⅱ區(qū)為0.05

采用LCD方法對式(4)所示的信號模型進(jìn)行分解，運行參數(shù)與圖2相同，其分解能力區(qū)域分布圖如圖3所示。參照圖2，繪制E值分別為0.05和1的兩條等高線，將a-F平面劃分為3個區(qū)域。與圖2對比可以發(fā)現(xiàn)，圖3中的區(qū)域Ⅲ明顯小于圖2中的區(qū)域Ⅲ。表明在相同的采樣頻率fs、最大頻率fmax和相同最大迭代次數(shù)等運行條件下，DOLCD方法分解能力范圍要明顯大于LCD方法的分解能力范圍，增強(qiáng)了抑制模態(tài)混淆的能力。

圖3　LCD方法分解能力

2.3分解能力的影響因素

需要說明的是，圖3所示的LCD分解能力圖與文獻(xiàn)[10]有所差異，這主要是因為信號模型的采樣頻率fs，以及fs與最大頻率fmax的關(guān)系對LCD分解效果的影響。圖2和圖3中，采樣信號模型的采樣頻率fs=8192Hz，高頻分量信號最大頻率fmax=80Hz。仍采用式(4)所示的信號模型，分別更改信號模型的采樣頻率fs和最大頻率fmax，分析DOLCD方法和LCD方法的分解能力。當(dāng)fs=4096Hz，高頻分量信號最大頻率fmax=80Hz時，以及當(dāng)fs=8192Hz，高頻分量信號最大頻率fmax=200 Hz時，DOLCD方法和LCD方法的分解能力分別如圖4a、圖4c和圖4b、圖4d所示。

對圖4分析可得，盡管圖4a和圖4c的區(qū)域Ⅲ與圖2中的區(qū)域Ⅲ，以及圖4b和圖4d的區(qū)域Ⅲ和圖3中的區(qū)域Ⅲ均有所不同，但圖4a和圖4c的區(qū)域Ⅲ均明顯大于圖4b和圖4d的區(qū)域Ⅲ。這也說明了DOLCD方法的分解能力要強(qiáng)于LCD方法。

3　仿真分析與應(yīng)用實例

3.1仿真分析

為了說明DOLCD方法的優(yōu)越性、不失一般性，仍采用式(4)所示的信號模型為分析對象，其中A1=2，A2=18，f1=16，f2=4，即x(t)=2sin(32πt)+18sin(8πt)，t∈[0,1]，其時域波形如圖5所示。

(a)fs=4096 Hz、fmax=80 Hz時，DOLCD方法的分解能力

(b)fs=4096 Hz、fmax=80 Hz時，LCD方法的分解能力

(c)fs=8192 Hz、fmax=200 Hz時，DOLCD方法的分解能力

(d)fs=8192 Hz、fmax=200 Hz時，LCD方法的分解能力圖4　兩組不同運行參數(shù)對應(yīng)的兩種方法的分解能力

(a)仿真信號x1(t)

(b)仿真信號x2(t)

(c)仿真信號x(t)圖5　仿真信號x(t)及其各成分的時域波形

信號x(t)是一個大幅值的低頻信號與小幅值的頻率略高信號的疊加信號，其幅值比a=A2/A1=9，頻率比F=f2/f1=0.25。根據(jù)對圖3的分析，信號x(t)處于圖3的區(qū)域Ⅰ，LCD無法將信號x(t)的兩分量有效分離。對圖6所示信號x(t)的LCD分解結(jié)果分析可得，LCD分解得到第一個分量I1(t)與目標(biāo)信號x(t)的幅值、頻率均區(qū)別不大；第二個分量I2(t)的幅值很小，從其時域圖上無法看出有物理意義的頻率成分。對其第一個分量I1(t)的幅值譜(圖7)進(jìn)行分析，可明顯看出分量I1(t)包含了目標(biāo)信號x(t)的兩種頻率成分。可見采用LCD對信號x(t)進(jìn)行分解時，第一個ISC分量I1(t)中包含有頻率差異較大的不同信號成分，出現(xiàn)了嚴(yán)重的模態(tài)混淆，LCD未能實現(xiàn)對仿真信號x(t)的有效分解，這與前面的分析是一致的。

(a)第一個ISC分量I1(t)

(b)第二個ISC分量I2(t)

(c)剩余項R2(t)圖6　仿真信號x(t)的LCD分解結(jié)果

圖7　LCD分解信號x(t)得到的分量I1(t)的幅值譜

對圖2進(jìn)行分析，信號x(t)處于圖2中區(qū)域Ⅲ，采用DOLCD可有效實現(xiàn)信號x(t)的兩分量分離，信號x(t)的DOLCD的分解結(jié)果如圖8所示。

(a)第一個ISC分量I1(t)

(b)第二個ISC分量I2(t)

為了進(jìn)一步比較兩種方法的分解效果，論文還對兩種方法得到的兩分量與對應(yīng)真實分量的相關(guān)系數(shù)[11]進(jìn)行了對比，相關(guān)系數(shù)指標(biāo)值越大，表示分解的準(zhǔn)確性越高。DOLCD得到的兩分量的相關(guān)系數(shù)指標(biāo)值為0.9999和1，LCD得到的兩分量的相關(guān)系數(shù)指標(biāo)值為0.1104和0.1713。這說明與LCD方法相比，DOLCD方法在分解精確度方面表現(xiàn)出一定的優(yōu)越性。

3.2應(yīng)用實例

為了進(jìn)一步說明DOLCD方法的有效性與實用性，論文將其應(yīng)用于不對中故障[12]的轉(zhuǎn)子振動位移實驗信號分析。實驗裝置如圖9所示，轉(zhuǎn)子轉(zhuǎn)速為3000 r/min，轉(zhuǎn)頻fr=50 Hz，實驗采樣頻率fs=2048 Hz，采樣時長為0.5 s。實驗數(shù)據(jù)的時域波形如圖10所示。其幅值譜如圖11所示。采用剛性聯(lián)軸器連接的轉(zhuǎn)子不對中時，強(qiáng)制連接所產(chǎn)生的力矩，使轉(zhuǎn)子發(fā)生彎曲變形；轉(zhuǎn)軸內(nèi)阻以及轉(zhuǎn)軸表面與旋轉(zhuǎn)體內(nèi)表面之間的摩擦而產(chǎn)生的相對滑動，使轉(zhuǎn)子產(chǎn)生自激旋轉(zhuǎn)周期性振動，其振動頻率為轉(zhuǎn)子轉(zhuǎn)頻fr的2倍，常伴頻率為轉(zhuǎn)頻fr的1倍頻，轉(zhuǎn)子振動位移信號會產(chǎn)生一個以轉(zhuǎn)頻為幅值調(diào)制頻率的高頻分量[13]。由圖10可以看出，高頻的分量信號被淹沒在強(qiáng)大的背景信號中。從圖11中主要看到的是與轉(zhuǎn)頻fr相關(guān)的基頻分量，看不出不對中故障信息。為提取高頻不對中故障信息，分別采用DOLCD和LCD對實驗數(shù)據(jù)進(jìn)行分解，分解結(jié)果如圖12和圖13所示，限于篇幅只畫出前三個分量。

圖9　轉(zhuǎn)子實驗裝置示意圖

圖10　實驗數(shù)據(jù)的時域波形

圖11　實驗數(shù)據(jù)的幅值譜

(a)第一個ISC分量I1(t)

(b)第二個ISC分量I2(t)

(c)第三個ISC分量I3(t)圖12　實驗數(shù)據(jù)的DOLCD分解結(jié)果

(a)第一個ISC分量I1(t)

(b)第二個ISC分量I2(t)

(c)第三個ISC分量I3(t)圖13　實驗數(shù)據(jù)的LCD分解結(jié)果

由圖12可以看出，DOLCD方法對實驗數(shù)據(jù)分解得到的第一個分量I1(t)是高頻背景噪聲信號，第二個分量I2(t)具有明顯的調(diào)制特征，包含了主要故障信息，分量I3(t)是與轉(zhuǎn)頻有關(guān)的背景信號。從圖13中可以看出，LCD方法對實驗數(shù)據(jù)分解得到的第一個分量I1(t)也是高頻背景噪聲信號，第二個分量I2(t)也能看出調(diào)制特征，第三個分量I3(t)是與轉(zhuǎn)頻有關(guān)的背景信號；但分量I2(t)和I3(t)均出現(xiàn)了嚴(yán)重的波形失真和模態(tài)混淆(如圖中實線所圈部分)。采用希爾伯特變換(Hilbert transform, HT)分別估計兩種方法得到的第二個分量的包絡(luò)譜，結(jié)果分別如圖14和圖15所示。

圖14　DOLCD的分量I2(t)的包絡(luò)譜

圖15　LCD的分量I2(t)的包絡(luò)譜

由圖14所示的DOLCD分量I2(t)的包絡(luò)譜圖可以看到明顯的轉(zhuǎn)頻fr的2倍頻和轉(zhuǎn)頻fr的1倍頻，這與轉(zhuǎn)子發(fā)生不對中故障時的頻率特征相吻合[14]。圖15所示的LCD分量I2(t)的包絡(luò)譜圖中未出現(xiàn)明顯的轉(zhuǎn)頻fr的2倍頻，無法有效識別實驗數(shù)據(jù)所包含的故障類型。

對上述實驗數(shù)據(jù)分析說明，相對LCD方法，DOLCD能更為有效地將高頻不對中信號從強(qiáng)大的背景信號中提取出來，實現(xiàn)不對中故障信號、背景信號和噪聲信號的分離，從而實現(xiàn)轉(zhuǎn)子故障診斷，同時也有效抑制了模態(tài)混淆。

4　結(jié)論

(1)采用仿真信號模型將DOLCD與LCD分解能力進(jìn)行對比分析，發(fā)現(xiàn)DOLCD方法分解能力范圍要明顯大于LCD方法的分解能力范圍，分解結(jié)果準(zhǔn)確性更高，也有效地抑制了模態(tài)混淆。

(2)將DOLCD方法應(yīng)用于轉(zhuǎn)子不對中故障信號的分析，分解得到了高頻調(diào)制特征分量和低頻分量，有效地實現(xiàn)了轉(zhuǎn)子不對中故障的診斷。

論文提出的DOLCD方法不需要事先選擇基函數(shù)，是一種自適應(yīng)的信號分析方法，可有效地應(yīng)用于機(jī)械設(shè)備故障診斷。本文僅是簡單的從信號模型兩組不同的采樣頻率和最大頻率對兩種方法分解能力的影響進(jìn)行了對比分析。筆者后續(xù)將進(jìn)一步對信號模型不同的運行參數(shù)對LCD分解能力的影響進(jìn)行系統(tǒng)性的研究。

[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 of London. Series A: Mathematical, Physical and Engineering Sciences, 1998, 454(1971): 903-995.

[2]Huang N E, Wu Z.A Review on Hilbert Huang Transform: Method and Its Applications to Geophysical Studies[J]. Reviews of Geophysics, 2008, 47(2): 1-23.

[3]Lei Yaguo, Lin Jing, He Zhengjia, et al. A Review on Empirical Mode Decomposition in Fault Diagnosis of Rotating Machinery[J]. Mechanical Systems and Signal Processing, 2013, 35(1): 108-126.

[4]時培明, 蔣金水, 劉彬, 等. 基于邊界特征尺度匹配延拓的EMD改進(jìn)方法及應(yīng)用[J]. 中國機(jī)械工程, 2014，25(12): 1616-1623.

Shi Peiming, Jiang Jinshui, Liu Bin, et al. Improved Method of EMD and Its Applications Based on Boundary Characteristic Scale Matching Extension Method[J]. China Mechanical Engineering, 2014，25(12): 1616-1623.

[5]程軍圣, 鄭近德, 楊宇. 一種新的非平穩(wěn)信號處理方法——局部特征尺度分解法[J]. 振動工程學(xué)報, 2012, 25(2): 215-220.

Cheng Junsheng, Zheng Jinde, Yang Yu. A Nonstationary Signal Analysis Approach——the Local Characteristic-scale Decomposition Method[J]. Journal of Vibration Engineering, 2012, 25(2): 215-220.

[6]Zheng Jinde, Cheng Junsheng, Yang Yu. Generalized Empirical Mode Decomposition and Its Applications to Rolling Element Bearing Fault Diagnosis[J]. Mechanical Systems and Signal Processing, 2013, 40(1): 136-153.

[7]鄭超, 郭奇, 郭麗杰, 等. 基于局部特征尺度分解的旋轉(zhuǎn)機(jī)械故障欠定盲源分離方法研究[J]. 燕山大學(xué)學(xué)報, 2014, 38(2): 168-174.

Zheng Chao, Guo Qi, Guo Lijie, et al. Underdetermined Blind Source Separation Method of Rotating Machinery Faults Based on Local Characteristic-scale Decomposition[J]. Journal of Yanshan University, 2014, 38(2):168-174.

[8]鄭近德, 程軍圣, 楊宇. 部分集成局部特征尺度分解：一種新的基于噪聲輔助數(shù)據(jù)分析方法[J]. 電子學(xué)報, 2013, 41(5): 1130-1135.

Zheng Jinde, Cheng Junsheng, Yang Yu. Partly Ensemble Local Characteristic-scale Decomposition: A New Noise Assisted Data Analysis Method[J]. Acta Electronica Sinica, 2013, 41(5): 1130-1135.

[9]劉海波, 玄志武. 利用微分算子增強(qiáng)EMD算法頻帶分解能力[J]. 振動與沖擊, 2013, 32(18): 133-135.

Liu Haibo, Xuan Zhiwu. Improving Frequency-band Separating Ability of EMD with a Differential Operator[J]. Journal of Vibration and Shock, 2013, 32(18): 133-135.

[10]楊宇, 曾鳴, 程軍圣. 局部特征尺度分解方法及其分解能力研究[J]. 振動工程學(xué)報, 2012, 25(5): 602-609.

Yang Yu, Zeng Ming, Cheng Junsheng. Research on Local Characteristic-scale Decomposition and Its Capacities[J]. Journal of Vibration Engineering, 2012, 25(5):602-609.

[11]Rilling G, Flandrin P. One or Two Frequencies? The Empirical Mode Decomposition Answers[J]. IEEE Transactions on Signal Processing, 2008, 56(1): 85-95.

[12]孫斌, 劉立遠(yuǎn), 雷偉. 基于正交局部保持映射的轉(zhuǎn)子故障診斷方法[J]. 中國機(jī)械工程, 2014, 25(16): 2219-2224.

Sun Bin, Liu Liyuan, Lei Wei. Rotor Fault Diagnosis Method on Orthogonal Locality Preserving Projection[J]. China Mechanical Engineering, 2014, 25(16): 2219-2224.

[13]盛兆順, 尹琦嶺. 設(shè)備狀態(tài)監(jiān)測與故障診斷技術(shù)及應(yīng)用[M]. 北京: 化學(xué)化工出版社, 2003.

[14]劉楊, 太興宇, 姚紅良, 等. 雙盤轉(zhuǎn)子軸承系統(tǒng)不對中-碰摩耦合故障分析[J]. 振動、測試與診斷, 2013, 33(5): 819-823.

Liu Yang, TaiXingyu, Yao Hongliang, et al. Study on Misalignment-rubbing Coupling Fault of Dual-disk Rotor-bearing System Including the Impact of Motor Coupling[J]. Journal of Vibration, Measurement & Diagnosis, 2013, 33(5): 819-823.

(編輯王旻玥)

Local Characteristic-scale Decomposition Method Based on Differential Operator and Its Applications

Wu Zhantao1Cheng Junsheng1Li Baoqing1Zheng Jinde2

1.State key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University,Changsha,410082 2.Anhui University of Technology,Ma’anshan,Anhui,243002

A novel non-stationary signal method differential operator based LCD(DOLCD) was proposed for improving the LCD method, of which, when the target signal contained components with frequencies close to each other, the decomposition ability was decreased, mode mixing phenomenon would easily occur and the physical meaning of intrinsic scale components decomposed were lost. A differential operator was used to differentiate the target signal in time domain, then the differentiated target signal was decomposed with DOLCD, then in the certain frequency ratio, the decomposition ability was enhanced, and the mode mixing was inhibited. The paper firstly studied the theory of DOLCD, then simulation experiments were used to compare the decomposition ability of DOLCD with LCD. The results indicate that DOLCD is more efficient in improving the decomposition ability, and inhibiting the mode mixing than LCD. Finally, the proposed method was applied to diagnose the rotor with misalignment faults successfully which indicated the effectiveness of DOLCD.

local characteristic-scale decomposition(LCD); differential operator; fault diagnosis; intrinsic scale component(ISC); rotor

2015-03-15

國家自然科學(xué)基金資助項目(51375152)；湖南省科技計劃資助項目(2014WK3005)

TH165.3；TH911.7DOI：10.3969/j.issn.1004-132X.2015.17.003

吳占濤，男，1982年生。湖南大學(xué)機(jī)械與運載工程學(xué)院博士研究生。主要研究方向為動態(tài)信號處理及機(jī)械設(shè)備故障診斷。程軍圣(通信作者)，男，1968年生。湖南大學(xué)機(jī)械與運載工程學(xué)院教授、博士、博士研究生導(dǎo)師。李寶慶，男，1984年生。湖南大學(xué)機(jī)械與運載工程學(xué)院博士研究生。鄭近德，男，1986年生。安徽工業(yè)大學(xué)機(jī)械工程學(xué)院講師、博士。