馮志鵬， 秦嗣峰(.北京科技大學(xué) 機(jī)械工程學(xué)院，北京　00083； .南車(chē)青島四方機(jī)車(chē)車(chē)輛股份有限公司，山東 青島　66)

?

基于Hilbert振動(dòng)分解和高階能量算子的行星齒輪箱故障診斷研究

馮志鵬1， 秦嗣峰2(1.北京科技大學(xué) 機(jī)械工程學(xué)院，北京100083； 2.南車(chē)青島四方機(jī)車(chē)車(chē)輛股份有限公司，山東 青島266111)

摘要：行星齒輪箱廣泛應(yīng)用于各種機(jī)械設(shè)備中，其故障診斷問(wèn)題是近年來(lái)的研究熱點(diǎn)之一。提出了基于Hilbert振動(dòng)分解和高階微分能量算子的故障診斷方法。Hilbert振動(dòng)分解計(jì)算復(fù)雜性低，能夠?qū)?fù)雜信號(hào)分解為單分量，應(yīng)用該方法對(duì)信號(hào)進(jìn)行分解，滿(mǎn)足高階微分能量算子的要求。高階微分能量算子的時(shí)間分辨率高，對(duì)信號(hào)的瞬態(tài)變化具有良好的自適應(yīng)性，應(yīng)用該方法檢測(cè)故障引起的瞬態(tài)沖擊，估計(jì)信號(hào)的幅值包絡(luò)和瞬時(shí)頻率。對(duì)高階微分能量算子輸出以及幅值包絡(luò)和瞬時(shí)頻率進(jìn)行Fourier變換，通過(guò)頻譜識(shí)別特征頻率，從而診斷行星齒輪箱故障。分析了行星齒輪箱的仿真信號(hào)和實(shí)驗(yàn)信號(hào)，準(zhǔn)確地診斷了太陽(yáng)輪、行星輪和齒圈的故障，驗(yàn)證了該方法的有效性。

關(guān)鍵詞：行星齒輪箱；故障診斷；高階能量算子；Hilbert振動(dòng)分解

行星齒輪箱結(jié)構(gòu)緊湊，傳動(dòng)比大，廣泛應(yīng)用于風(fēng)力發(fā)電機(jī)、直升機(jī)、船舶和車(chē)輛等裝備中[1]。實(shí)際運(yùn)行中，行星齒輪箱工作環(huán)境惡劣，容易出現(xiàn)故障。因此，開(kāi)展行星齒輪箱故障診斷研究具有重要意義。

目前，國(guó)內(nèi)外學(xué)者已經(jīng)提出了一些行星齒輪箱故障診斷方法。McFadden[2-3]推廣了時(shí)域平均方法用于分析行星齒輪傳動(dòng)系統(tǒng)振動(dòng)信號(hào)。Mosher[4]研究了行星齒輪箱振動(dòng)信號(hào)的頻譜結(jié)構(gòu)，發(fā)現(xiàn)頻譜具有不對(duì)稱(chēng)特點(diǎn)。Inalpolat等[5-6]建立了行星齒輪箱非線(xiàn)性時(shí)變動(dòng)力學(xué)模型，研究了行星架旋轉(zhuǎn)和齒輪制造誤差對(duì)振動(dòng)信號(hào)的調(diào)制作用及其邊帶特征。Mark[7-8]研究行星輪載荷的不均勻分配和行星架旋轉(zhuǎn)引起的調(diào)制作用以及振動(dòng)信號(hào)的邊帶特征。文獻(xiàn)[9-14]考慮行星齒輪傳動(dòng)機(jī)構(gòu)的動(dòng)力學(xué)特點(diǎn)以及齒輪的常見(jiàn)故障形式，建立了正常、分布式故障和局部故障狀態(tài)下的振動(dòng)信號(hào)模型，推導(dǎo)了太陽(yáng)輪、行星輪和齒圈分布式故障和局部故障特征頻率計(jì)算公式，推導(dǎo)了振動(dòng)信號(hào)頻譜的解析表達(dá)式，提出了幅值解調(diào)和頻率解調(diào)分析方法，總結(jié)歸納了太陽(yáng)輪、行星輪和齒圈等全部三種齒輪分布式故障和局部故障的征兆規(guī)律，針對(duì)工程實(shí)踐中常見(jiàn)的時(shí)變工況，提出了時(shí)頻分析診斷方法。

行星齒輪箱故障診斷的實(shí)質(zhì)是根據(jù)故障特征頻率及其幅值的變化識(shí)別故障原因。雖然通過(guò)Fourier頻譜分析能夠診斷各個(gè)齒輪的故障，但是其中涉及復(fù)雜的邊帶分析。齒輪故障引起的周期性沖擊的重復(fù)頻率、以及沖擊對(duì)嚙合振動(dòng)幅值和頻率的調(diào)制頻率與齒輪故障的特征頻率密切相關(guān)[9-13]。分別獨(dú)立分析周期性沖擊或調(diào)幅和調(diào)頻成分，其Fourier頻譜不涉及卷積，結(jié)構(gòu)簡(jiǎn)單，有望避免行星齒輪箱振動(dòng)信號(hào)固有的調(diào)幅調(diào)頻現(xiàn)象在Fourier頻譜中引起的復(fù)雜邊帶結(jié)構(gòu)，有效地診斷行星齒輪箱故障。

近年來(lái)，Maragos等[15]在Teager能量算子的基礎(chǔ)上提出了高階微分能量算子，通過(guò)信號(hào)的瞬時(shí)幅值及其微分的非線(xiàn)性組合可以估計(jì)信號(hào)的‘能量’。另外，基于高階微分能量算子的能量分離算法還可以估計(jì)任意調(diào)制信號(hào)的幅值包絡(luò)和瞬時(shí)頻率等參數(shù)。高階微分能量算子在信號(hào)及其導(dǎo)數(shù)的基礎(chǔ)上進(jìn)行計(jì)算，算法簡(jiǎn)單，時(shí)間分辨率高，對(duì)信號(hào)的瞬態(tài)變化具有良好的適應(yīng)性，因此為檢測(cè)信號(hào)瞬態(tài)特征以及實(shí)現(xiàn)幅值頻率解調(diào)提供了一種有效途徑。但是，高階微分能量算子僅適用于單分量信號(hào)。對(duì)于成分復(fù)雜的行星齒輪箱振動(dòng)信號(hào)，需要將其分解為單分量再進(jìn)行分析。

Hilbert振動(dòng)分解可以將多分量信號(hào)分解為多個(gè)單分量疊加的形式。相對(duì)于經(jīng)驗(yàn)?zāi)Ｊ椒纸夥椒ǘ?，Hilbert振動(dòng)分解得到的單分量不會(huì)出現(xiàn)模態(tài)混淆現(xiàn)象，而且計(jì)算簡(jiǎn)單，工作效率高[16]。

考慮高階微分能量算子和Hilbert振動(dòng)分解方法的獨(dú)特優(yōu)勢(shì)，本文將二者結(jié)合，發(fā)揮它們各自的優(yōu)勢(shì)，進(jìn)行行星齒輪箱故障診斷。

1高階微分能量算子

高階微分能量算子只需要信號(hào)幅值及其導(dǎo)數(shù)便可計(jì)算信號(hào)任意時(shí)刻的‘能量’，計(jì)算復(fù)雜性低，算法效率高，對(duì)信號(hào)的瞬態(tài)變化具有良好的時(shí)間分辨率和自適應(yīng)性，能夠準(zhǔn)確估計(jì)任意調(diào)制信號(hào)的幅值包絡(luò)和瞬時(shí)頻率。另外，高階能量算子的輸出為信號(hào)瞬時(shí)幅值的平方與瞬時(shí)頻率的k次方之積，相對(duì)于Teager能量算子，其輸出增加了和頻率的k-2次方的乘積，更能夠突出沖擊的瞬態(tài)特征。

對(duì)于任意連續(xù)信號(hào)x(t)，高階能量算子定義為[15]

k=0,±1,±2,…

(1)

式中，k為高階能量算子的階次，為任意整數(shù)，x(k)為x(t)的第k階導(dǎo)數(shù)

(2)

本文關(guān)注k為正數(shù)時(shí)的高階能量算子。通常，高階能量算子可由低階能量算子通過(guò)遞歸推導(dǎo)得到

(3)

(4)

其中

(5)

Ek=-ω2Ek-2,Ek=γk[Acos(ωt+θ)]

(6)

初始條件為E0=-A2,E1=0。由式(6)遞推，可得

(7)

由式(7)可見(jiàn)，k階能量算子的輸出相當(dāng)于在二階Teager能量算子輸出的基礎(chǔ)上又乘以信號(hào)的頻率的k-2次方ωk-2，

對(duì)于簡(jiǎn)諧信號(hào)x(t)=Acos(ωt+φ)，由四階和二階能量算子可以估計(jì)信號(hào)的頻率和幅值

(8)

式(1)針對(duì)連續(xù)時(shí)間信號(hào)定義，對(duì)于離散時(shí)間信號(hào)x(n)，用差分代替微分，可得離散高階能量算子

γk(x[n])=x[n]x[n+k-2]-

x[n-1]x[n+k-1]

(9)

其中，k=0,1,2,…。

對(duì)于離散時(shí)間簡(jiǎn)諧信號(hào)x(n)=Acos(Ωn+θ)，高階能量算子輸出

γk[Acos(Ωn+θ)]=A2sin(Ω)sin[(k-1)Ω]

(10)

由三階和二階離散能量算子可以估計(jì)簡(jiǎn)諧信號(hào)的頻率和幅值

(11)

2Hilbert振動(dòng)分解

Hilbert振動(dòng)分解應(yīng)用Hilbert變換和同步檢測(cè)解調(diào)估計(jì)信號(hào)分量的瞬時(shí)頻率和包絡(luò)幅值，按照信號(hào)分量幅值包絡(luò)由大到小的順序，將復(fù)雜信號(hào)分解為單分量成分。該方法通過(guò)Hilbert變換實(shí)現(xiàn)，不涉及其他復(fù)雜的信號(hào)處理技術(shù)，算法簡(jiǎn)單[16]。

Hilbert振動(dòng)分解假設(shè)信號(hào)滿(mǎn)足以下條件：① 信號(hào)由若干個(gè)準(zhǔn)諧波疊加而成；② 每個(gè)信號(hào)分量的幅值包絡(luò)大小不同；③ 信號(hào)長(zhǎng)度涵蓋變化最緩慢分量的多個(gè)周期。

對(duì)于多分量信號(hào)，在基于Hilbert變換和解析信號(hào)估計(jì)的瞬時(shí)頻率中，以及通過(guò)同步檢測(cè)解調(diào)得到的同相位信號(hào)和相位移動(dòng)正交信號(hào)中，都包含與幅值最大分量對(duì)應(yīng)的慢變部分以及其它分量綜合作用引起的快變部分，而且快變部分對(duì)時(shí)間積分等于零，因此可以通過(guò)低通濾波估計(jì)幅值最大分量的慢變瞬時(shí)頻率、幅值包絡(luò)和瞬時(shí)相位。在此基礎(chǔ)上，Hilbert振動(dòng)分解的基本過(guò)程如下：

(1) 對(duì)于信號(hào)x(t)，基于Hilbert變換構(gòu)造其解析信號(hào)，并估計(jì)瞬時(shí)頻率ω(t)和瞬時(shí)幅值a(t)。

(2) 對(duì)瞬時(shí)頻率ω(t)進(jìn)行低通濾波得到幅值最大分量對(duì)應(yīng)的瞬時(shí)頻率ω1(t)。

(3) 以ω1(t)作為參考信號(hào)的頻率，對(duì)信號(hào)x(t)進(jìn)行同步檢測(cè)解調(diào)和低通濾波，得到幅值最大分量的幅值包絡(luò)a1(t)和瞬時(shí)相位φ1(t)，重構(gòu)幅值最大分量x1(t)=a1(t)cos[∫(t)dt]。

(4) 從原始信號(hào)減去幅值最大分量，得到新信號(hào)xl-1(t)=x(t)-x1(t)，令x(t)=xl-1(t)，重復(fù)步驟(1)～(4)，直到先后兩次迭代結(jié)果之間的標(biāo)準(zhǔn)差小于設(shè)定的閾值(本文中為0.001)，停止迭代。

每次迭代中，瞬時(shí)頻率和包絡(luò)幅值通過(guò)低通濾波得到，濾波器的截止頻率決定了Hilbert振動(dòng)分解方法的頻率分辨率。因此，在滿(mǎn)足濾波器精度和穩(wěn)定性要求的前提下，截止頻率越小越好。

3分析過(guò)程

基于Hilbert振動(dòng)分解和高階微分能量算子的分析方法基本步驟如下：

(1) 應(yīng)用Hilbert振動(dòng)分解方法將信號(hào)分解為若干單分量成分。

(2) 應(yīng)用高階能量算子估計(jì)各分量的瞬時(shí)頻率、幅值包絡(luò)和‘能量’。

(3) 選擇敏感的單分量進(jìn)行后續(xù)分析。齒輪故障振動(dòng)信號(hào)具有調(diào)幅-調(diào)頻特征，幅值和頻率調(diào)制頻率為故障齒輪的特征頻率及其倍頻，載波頻率為嚙合頻率及其倍頻，瞬時(shí)頻率圍繞嚙合頻率或其倍頻上下波動(dòng)的分量包含齒輪故障信息；另外，Hilbert振動(dòng)分解按照幅值包絡(luò)由大到小的順序分解出信號(hào)的組成分量。因此選擇最先分解得到的且瞬時(shí)頻率圍繞嚙合頻率或其倍頻上下波動(dòng)的單分量作為敏感分量進(jìn)行分析。

(4) 對(duì)敏感單分量的瞬時(shí)頻率、幅值包絡(luò)及‘能量’進(jìn)行Fourier變換，根據(jù)頻譜中的峰值頻率以及齒輪的故障特征頻率診斷故障。

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

行星齒輪箱振動(dòng)信號(hào)可以用調(diào)幅-調(diào)頻信號(hào)模型描述[9-14,17]。不失一般性，假設(shè)太陽(yáng)輪出現(xiàn)故障，只考慮載波頻率和調(diào)制頻率的基頻，即齒輪嚙合頻率和故障齒輪的特征頻率，則振動(dòng)信號(hào)模型簡(jiǎn)化為

x(t)=[1-cos(2πfsrt)][1+Acos(2πfst)]×

cos[2πfmt+Bsin(2πfst)+φ]+n(t)

(12)

式中，fsr=12 Hz為太陽(yáng)輪絕對(duì)旋轉(zhuǎn)頻率，fs=30 Hz為太陽(yáng)輪故障特征頻率，fm=500 Hz為嚙合頻率，A=B=1分別為調(diào)幅和調(diào)頻的調(diào)制系數(shù)，φ=0位初相位。為了模擬實(shí)際測(cè)試中的背景噪聲干擾，在仿真信號(hào)中加入信噪比為6 dB的Gauss白噪聲n(t)。信號(hào)采樣頻率為4 096 Hz。

信號(hào)分析結(jié)果如圖1所示，圖1(a)為信號(hào)的時(shí)域波形，圖1(b)為Hilbert振動(dòng)分解得到的前7個(gè)分量。計(jì)算各分量的瞬時(shí)頻率，發(fā)現(xiàn)第一個(gè)分量的瞬時(shí)頻率(圖1(c))圍繞齒輪嚙合頻率波動(dòng)，因此選擇第一個(gè)分量進(jìn)一步分析?；谒碾A和二階能量算子計(jì)算得到所選分量的幅值包絡(luò)、瞬時(shí)頻率及四階能量，并進(jìn)行Fourier變換，結(jié)果如圖1(d)～(f)所示。在圖1(d)包絡(luò)譜中，太陽(yáng)輪故障特征頻率fs、太陽(yáng)輪旋轉(zhuǎn)頻率fsr以及它們的組合頻率fs±fsr處的峰值明顯。在圖1(e)瞬時(shí)頻率Fourier頻譜中，太陽(yáng)輪故障特征頻率fs及其2倍頻2fs峰值明顯。在圖1(f)四階能量Fourier頻譜中，太陽(yáng)輪故障特征頻率fs、太陽(yáng)輪旋轉(zhuǎn)頻率fsr以及它們的組合頻率fs±fsr、還有太陽(yáng)輪故障特征頻率的2倍頻2fs以及它和太陽(yáng)輪旋轉(zhuǎn)頻率fsr的組合頻率2fs±fsr處的峰值明顯。上述分析結(jié)果符合太陽(yáng)輪故障的振動(dòng)信號(hào)特征。

圖1　仿真信號(hào)分析Fig.1 Analysis result of a simulated signal

5實(shí)驗(yàn)信號(hào)分析

5.1實(shí)驗(yàn)說(shuō)明

某行星齒輪箱實(shí)驗(yàn)系統(tǒng)如圖2示，其中行星齒輪箱結(jié)構(gòu)參數(shù)見(jiàn)表1。為測(cè)試振動(dòng)信號(hào)，在齒輪箱底座和箱體表面布置了多個(gè)加速度傳感器，本文分析與齒圈相連的箱體頂部的傳感器信號(hào)，信號(hào)采樣頻率為16 384 Hz。實(shí)驗(yàn)過(guò)程中，和太陽(yáng)輪相連的輸入軸轉(zhuǎn)速為15.95 Hz。根據(jù)齒輪箱結(jié)構(gòu)參數(shù)和轉(zhuǎn)速，計(jì)算得到齒輪局部故障特征頻率，見(jiàn)表2。

表1　行星齒輪箱參數(shù)(單位：Hz)

表2　行星齒輪箱特征頻率(單位：Hz)

為模擬行星齒輪箱中的齒輪局部故障，在太陽(yáng)輪、行星輪和齒圈的一個(gè)輪齒上加工了局部磨損剝落，見(jiàn)圖3。

圖2　實(shí)驗(yàn)系統(tǒng)Fig.2 Experimental system

圖3　齒輪局部損傷Fig.3 Gear localized damage

5.2實(shí)驗(yàn)信號(hào)分析

5.2.1正常信號(hào)

圖4為正常信號(hào)分析結(jié)果，其中圖4 (a)～(c)分別為信號(hào)時(shí)域波形、Hilbert振動(dòng)分解得到的前7個(gè)單分量成分、第一個(gè)分量的瞬時(shí)頻率。由于第一個(gè)分量的瞬時(shí)頻率圍繞嚙合頻率的9倍頻波動(dòng)，因此選擇第一個(gè)分量進(jìn)一步分析，在此基礎(chǔ)上應(yīng)用四階和二階能量算子和Fourier變換得到第一個(gè)分量的包絡(luò)譜、瞬時(shí)頻率Fourier頻譜、四階能量Fourier頻譜，如圖4(d)～(f)示。其中的峰值主要出現(xiàn)在行星架旋轉(zhuǎn)頻率fc及其倍頻nfc處，這是由于行星架旋轉(zhuǎn)對(duì)齒輪箱嚙合振動(dòng)的調(diào)制效應(yīng)及齒輪制造安裝誤差所造成的，這些現(xiàn)象符合正常信號(hào)的理論特征。

5.2.2太陽(yáng)輪故障信號(hào)

圖5為太陽(yáng)輪故障信號(hào)的分析結(jié)果。經(jīng)過(guò)Hilbert振動(dòng)分解得到的第一個(gè)分量的瞬時(shí)頻率圍繞嚙合頻率的13倍頻波動(dòng)，因此仍然選擇第一個(gè)分量進(jìn)一步分析。由圖5(d)包絡(luò)譜可見(jiàn)，太陽(yáng)輪旋轉(zhuǎn)頻率的2倍頻2fsr及3倍頻3fsr占主導(dǎo)地位，其它峰值主要出現(xiàn)在太陽(yáng)輪故障特征頻率fs及其倍頻nfs、太陽(yáng)輪旋轉(zhuǎn)頻率fsr及其倍頻mfsr、以及它們的頻率組合nfs±m(xù)fsr處。由圖5(e)瞬時(shí)頻率Fourier頻譜可見(jiàn)，頻譜特征更為簡(jiǎn)單，峰值主要出現(xiàn)在太陽(yáng)輪故障特征頻率fs及其倍頻nfs、行星架旋轉(zhuǎn)頻率的2倍頻2fc以及它們的組合nfs±m(xù)fc處。由圖5(f)四階能量Fourier頻譜可見(jiàn)，太陽(yáng)輪故障特征頻率fs及其倍頻nfs、太陽(yáng)輪旋轉(zhuǎn)頻率fsr及其倍頻mfsr、以及它們的頻率組合nfs±m(xù)fsr處的峰值明顯，而且太陽(yáng)輪故障特征頻率更突出。上述特征說(shuō)明太陽(yáng)輪出現(xiàn)了故障，符合實(shí)驗(yàn)中的實(shí)際情況。

5.2.3行星輪故障信號(hào)

圖6為行星輪故障信號(hào)的分析結(jié)果。經(jīng)過(guò)Hilbert振動(dòng)分解得到的第一個(gè)分量的瞬時(shí)頻率圍繞嚙合頻率的7倍頻波動(dòng)，因此選擇該分量進(jìn)一步分析。由圖6(d)包絡(luò)譜和圖6(e)瞬時(shí)頻率Fourier頻譜可見(jiàn)，行星輪故障特征頻率的3倍頻與行星架旋轉(zhuǎn)頻率的組合3fp+fc占主導(dǎo)地位，其他峰值主要出現(xiàn)在行星輪故障特征頻率fp及其倍頻nfp以及行星架旋轉(zhuǎn)頻率fc及其倍頻nfc處。一些峰值和行星架旋轉(zhuǎn)頻率及其倍頻有關(guān)，是因?yàn)樾行禽喌木植抗收峡赡茉斐尚行羌茌d荷分配不均，使得行星架旋轉(zhuǎn)運(yùn)動(dòng)對(duì)嚙合振動(dòng)的調(diào)制作用增強(qiáng)，導(dǎo)致行星架旋轉(zhuǎn)頻率及其倍頻處峰值增大。由圖6(f)四階能量Fourier頻譜可見(jiàn)，行星輪故障特征頻率fp及其倍頻nfp以及行星架旋轉(zhuǎn)頻率fc及其倍頻nfc處峰值明顯，而且行星輪故障特征頻率及其倍頻突出，表明行星輪出現(xiàn)了損傷，符合實(shí)驗(yàn)中的實(shí)際情況。

圖4　正常齒輪信號(hào)Fig. 4 Normal gear signal

圖5　太陽(yáng)輪故障信號(hào)Fig.5 Faulty sun gear signal

圖6　行星輪故障信號(hào)Fig.6 Faulty planet gear signal

5.2.4齒圈故障信號(hào)

圖7為齒圈故障信號(hào)的分析結(jié)果。經(jīng)過(guò)Hilbert振動(dòng)分解得到的第一個(gè)分量的瞬時(shí)頻率圍繞嚙合頻率的14倍頻波動(dòng)，因此選擇該分量進(jìn)一步分析。由于制造誤差在所難免，和齒圈嚙合的三個(gè)行星輪不可能完全相同，因此考慮齒圈故障特征頻率的1/3及其倍頻nfr/3。由圖7(d)包絡(luò)譜和圖7(e)瞬時(shí)頻率Fourier頻譜可見(jiàn)，齒圈故障特征頻率的2/3占主導(dǎo)地位，其余峰值主要出現(xiàn)在齒圈故障特征頻率的1/3及其倍頻nfr/3處。由圖7(f)四階能量Fourier頻譜可見(jiàn)，齒圈故障特征頻率的1/3及其倍頻nfr/3處的峰值更為突出，說(shuō)明齒圈出現(xiàn)了故障，符合實(shí)驗(yàn)中的實(shí)際情況。

圖7　齒圈故障信號(hào)Fig.7 Faulty ring gear signal

6結(jié)論

行星齒輪箱的安裝制造誤差、振動(dòng)傳遞路徑變化和齒輪損傷會(huì)對(duì)嚙合振動(dòng)信號(hào)產(chǎn)生調(diào)幅調(diào)頻效應(yīng)，導(dǎo)致信號(hào)的Fourier頻譜十分復(fù)雜。提出了基于Hilbert振動(dòng)分解與高階微分能量算子的分析方法，可以有效避免復(fù)雜的邊帶分析。提出了敏感信號(hào)分量的選取原則。分析了行星齒輪箱故障仿真信號(hào)和實(shí)驗(yàn)信號(hào)，準(zhǔn)確診斷出了太陽(yáng)輪、行星輪和齒圈的局部損傷故障，驗(yàn)證了該方法的有效性。

參 考 文 獻(xiàn)

[ 1 ] 雷亞國(guó)，何正嘉，林京,等．行星齒輪箱故障診斷技術(shù)的研究進(jìn)展[J] ．機(jī)械工程學(xué)報(bào)，2011，47(19) ：59-67．

LEI Ya-guo，HE Zheng-jia，LIN Jing，et al．Research advances of fault diagnosis technique for planetary gearboxes[J]．Journal of Mechanical Engineering，2011，47(19)：59-67.

[ 2 ] McFadden P D． A technique for calculating the time domain averages of the vibration of the individual planet gears and the sun gear in an epicyclic gearbox[J]．Journal of Sound and Vibration，1991，144：163-172．

[ 3 ] McFadden P D．Window functions for the calculation of the time domain averages of the vibration of the individual planet gears and sun gear in an epicyclic gearbox[J]．Journal of Sound and Vibration，1994，144：179-187．

[ 4 ] Mosher M．Understanding vibration spectra of planetary gear systems for fault detection[C]//Proceedings of ASME Design Engineering Technical Conferences，2003：1-8．

[ 5 ] Inalpolat M，Kahraman A．A dynamic model to predict modulation sidebands of a planetary gear set having manufacturing errors[J]．Journal of Sound and Vibration，2010，329：371-393．

[ 6 ] Inalpolat M，Kahraman A．A theoretical and experimental investigation of modulation sidebands of planetary gear sets[J]．Journal of Sound and Vibration ，2009，329：371-393．

[ 7 ] Mark W D．Stationary transducer response to planetary-gear vibration excitation with non-uniform planet loading[J]．Mechanical Systems and Signal Processing，2009，23：1366-1381．

[ 8 ] Mark W D．Stationary transducer response to planetary-gear vibration excitation II: Effects of torque modulations[J]．Mechanical Systems and Signal Processing，2009，23：2253-2259．

[ 9 ] 馮志鵬，褚福磊．行星齒輪箱齒輪分布式故障振動(dòng)頻譜特征[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(2)：118-125．

FENG Zhi-peng，CHU Fu-lei．Vibration spectral characteristics of distributed gear fault of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(2)：118-125.

[10] 馮志鵬，趙鐳鐳，褚福磊．行星齒輪箱齒輪局部故障振動(dòng)頻譜特征[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(5)：119-127．

FENG Zhi-peng，ZHAO Lei-lei，CHU Fu-lei．Vibration spectral characteristics of localized gear fault of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(5)：119-127．

[11] 馮志鵬，趙鐳鐳，褚福磊．行星齒輪箱故障診斷的幅值解調(diào)分析方法[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(8)：107-111．

FENG Zhi-peng，ZHAO Lei-lei，CHU Fu-lei．Amplitude demodulation analysis method for fault diagnosis of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(8)：107-111.

[12] 馮志鵬，褚福磊．行星齒輪箱故障診斷的頻率解調(diào)分析方法[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(11)：112-117．

FENG Zhi-peng，CHU Fu-lei．Frequency demodulation analysis method for fault diagnosis of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(11)：112-117.

[13] 馮志鵬，褚福磊．行星齒輪箱故障診斷的扭轉(zhuǎn)振動(dòng)信號(hào)分析方法[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(14)：101-106．

FENG Zhi-peng，CHU Fu-lei．Torisonal vibration signal analysis method for fault diagnosis of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(14)：101-106.

[14] 馮志鵬，范寅夕，Liang Ming，等．行星齒輪箱故障診斷的非平穩(wěn)振動(dòng)信號(hào)分析方法[J]．中國(guó)電機(jī)工程學(xué)報(bào)，2013，33(17)：105-110．

FENG Zhi-peng，F(xiàn)AN Yin-xi，LIANG Ming，et al．A nonstationary vibration signal analysis method for fault diagnosis of planetary gearboxes[J]．Proceedings of the CSEE，2013，33(17)：105-110.

[15] Maragos P，Potamianos A．Higher order differential energy operators[J]．IEEE Signal Processing，1995，2(8)：152-154．

[16] Feldman M．Time-varying vibration decomposition and analysis based on the Hilbert transform[J]．Journal of Sound and Vibration，2006，295：518-530．

[17] McFadden P D．Detecting fatigue cracks in gears by amplitude and phase demodulation of the meshing vibration[J]．Journal of Vibration Acoustics Stress and Reliability in Design-Transactions of the ASME，1986， 108：165-170．

Planetary gearbox fault diagnosis based on Hilbert vibration decomposition and higher order differential energy operator

FENGZhi-peng1,QINSi-feng2(1. School of Mechanical Engineering, University of Science and Technology Beijing, Beijing 100083, China; 2. CSR Qingdao Sifang Co., Ltd, Qingdao 266111, China)

Abstract:Planetary gearboxes are widely used in various types of machinery. Fault diagnosis of planetary gearboxes is very important. A method based on Hilbert vibration decomposition and higher order differential energy operator was proposed for planetary gearbox fault diagnosis. By adopting its advantages of low computational complexity and mono-component decomposition ability, Hilbert vibration decomposition was used to decompose a signal and meet the mono-component requirement for using the higher order differential energy operator. Considering its merits of fine time resolution and good adaptability to instantaneous changes of signals, the higher order differential energy operator was applied to detect the transient impulses induced by gear faults and to estimate the corresponding amplitude envelope and instantaneous frequency. Fourier transformation was applied in the higher order differential energy operator output and the obtained amplitude envelope and instantaneous frequency to identify the gear fault characteristic frequencies and thereby diagnose planetary gearbox faults. The effectiveness of the proposed method was validated with both simulated signals and test ones of a planetary gearbox.

Key words:planetary gearbox; fault diagnosis; higher order differential energy operator; Hilbert vibration decomposition

中圖分類(lèi)號(hào):TH165.3

文獻(xiàn)標(biāo)志碼：A

DOI:10.13465/j.cnki.jvs.2016.05.007

收稿日期：2015-01-27修改稿收到日期：2015-03-20

基金項(xiàng)目：國(guó)家自然科學(xué)基金(11272047);教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃(NCET-12-0775)

第一作者 馮志鵬 男，教授, 博士生導(dǎo)師，1973年12月生