丁建明， 王 晗 ， 林建輝 ， 黃晨光

(1.西南交通大學(xué) 牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室，成都 610031；2. 青島四方機(jī)車車輛股份有限公司 高速列車系統(tǒng)集成國家工程實(shí)驗(yàn)室，青島 266111)

基于EMD-Hankel-SVD的高速列車萬向軸動(dòng)不平衡檢測

丁建明1， 王 晗2， 林建輝1， 黃晨光1

(1.西南交通大學(xué) 牽引動(dòng)力國家重點(diǎn)實(shí)驗(yàn)室，成都 610031；2. 青島四方機(jī)車車輛股份有限公司 高速列車系統(tǒng)集成國家工程實(shí)驗(yàn)室，青島 266111)

針對(duì)EMD(Empirical Model Decomposition)存在模式頻率混疊帶來的頻譜雜亂的根本缺陷，提出一種高速列車萬向軸動(dòng)不平衡動(dòng)態(tài)檢測的新方法。該方法的核心是：對(duì)萬向節(jié)安裝機(jī)座的振動(dòng)信號(hào)進(jìn)行EMD分解得到基本模式分量，應(yīng)用基本模式分量信號(hào)來構(gòu)造Hankel矩陣，對(duì)該矩陣進(jìn)行奇異值正交化分解，以奇異值關(guān)鍵疊層作為奇異值的選擇準(zhǔn)則對(duì)信號(hào)進(jìn)行重構(gòu)，應(yīng)用重構(gòu)信號(hào)的傅里葉譜來檢測高速列車萬向軸的動(dòng)不平衡，消除EMD分解模式頻率混疊帶來頻譜雜亂，提高了譜的清晰度，凸顯了故障特征。應(yīng)用萬向軸動(dòng)不平衡試驗(yàn)數(shù)據(jù)對(duì)該方法進(jìn)行試驗(yàn)驗(yàn)證，結(jié)果表明：該方法能夠有效檢測萬向軸動(dòng)不平衡引起故障特征和萬向軸的固有振動(dòng)特性，與純EMD方法相比，該方法在譜的清晰度和故障表征力上得到了顯著提高。

高速列車；萬向軸動(dòng)不平衡；EMD；Hankel矩陣；SVD；動(dòng)態(tài)檢測

轉(zhuǎn)向架是高速列車的核心部件，起到承載、導(dǎo)向、支撐、動(dòng)力傳遞的關(guān)鍵作用，直接決定了高速列車的運(yùn)行品質(zhì)和運(yùn)行安全。萬向傳動(dòng)軸是動(dòng)力傳遞的核心元件，在CRH5型動(dòng)車組的動(dòng)力傳遞結(jié)構(gòu)中，牽引電機(jī)采用體懸結(jié)構(gòu)，齒輪箱采用抱軸式結(jié)構(gòu)[1]，萬向軸既要傳遞牽引力矩，又要適應(yīng)復(fù)雜的運(yùn)動(dòng)關(guān)系[2]，其工作環(huán)境惡劣，同時(shí)萬向軸為細(xì)長結(jié)構(gòu)，其彎曲剛度、扭轉(zhuǎn)剛度都很小[3]，這些原因促使車輛在高速運(yùn)行過程中萬向軸極易產(chǎn)生偏心。加之，車輛在長期的運(yùn)行中，萬向節(jié)軸的磨損間隙、萬向傳動(dòng)軸平衡滑塊的松動(dòng)等也會(huì)導(dǎo)致萬向傳動(dòng)軸的偏心。偏心使萬向傳動(dòng)軸產(chǎn)生動(dòng)不平衡附加力矩，加劇傳動(dòng)系統(tǒng)的振動(dòng)、極易快速破壞傳動(dòng)系統(tǒng)的軸承、萬向節(jié)等動(dòng)力傳遞部件，嚴(yán)重者導(dǎo)致動(dòng)力中斷、機(jī)破、車毀人亡等重特大事故的發(fā)生。為了保障動(dòng)力的安全傳遞，開展萬向軸動(dòng)不平衡車載動(dòng)態(tài)檢測技術(shù)的研究顯得十分迫切和必要。

EMD分解十分適合處理非線性、非平穩(wěn)信號(hào)，突破了基于傅里葉分解、小波分解受測不準(zhǔn)原理的約束，能夠無限提高頻率、時(shí)間的分辨率，克服了傅里葉變換分析非平穩(wěn)信號(hào)失去物理意義的根本缺陷，克服了小波分解對(duì)小波的基函數(shù)依賴，是一種自適應(yīng)的數(shù)據(jù)分解方法[4-6]。這些優(yōu)良特性決定了EMD具有廣闊的應(yīng)用前景，已在齒輪箱、軸承、發(fā)電機(jī)組碰擦、橋梁健康等故障診斷中得到了成功應(yīng)用[7-9]。

盡管如此，作者應(yīng)用EMD來提取高速列車萬向軸不平衡故障的研究中，發(fā)現(xiàn)頻譜十分雜亂，很難檢測萬向軸不平衡引起的故障特征和萬向軸的固有振動(dòng)特征，極易導(dǎo)致誤警、虛警。探究其原因，EMD分解同小波分解一樣，同樣存在相鄰分解模式間存在頻率重疊混淆的現(xiàn)象。造成同一內(nèi)稟模態(tài)函數(shù)IMF的細(xì)節(jié)信號(hào)的頻帶變寬，相鄰內(nèi)稟模態(tài)函數(shù)的頻率重疊較大，分解得到的內(nèi)稟模態(tài)函數(shù)的噪聲變大，頻率成分復(fù)雜，淹沒和混淆故障特征和固有振動(dòng)特征。為了充分發(fā)揮EMD在處理非線性、非平穩(wěn)信號(hào)中的獨(dú)特優(yōu)勢，同時(shí)克服EMD分解與經(jīng)典小波變換一樣存在IMF分量之間頻帶重疊的共有缺陷，將Hankel矩陣、奇異值分解理論引入到EMD分解中，利用不同特性信號(hào)的奇異值的分布規(guī)律的不同來減少頻率重疊效應(yīng)，突顯故障特征和固有振動(dòng)特征。

為此，為了克服EMD存在模式頻率混疊帶來的頻譜雜亂的根本缺陷，提出一種高速列車萬向軸動(dòng)不平衡動(dòng)態(tài)檢測的新方法。該方法的核心是：對(duì)萬向節(jié)安裝機(jī)座的振動(dòng)信號(hào)進(jìn)行EMD分解得到基本模式分量，應(yīng)用基本模式分量信號(hào)來構(gòu)造Hankel矩陣，對(duì)該矩陣進(jìn)行奇異值正交化分解，以奇異值關(guān)鍵疊層作為奇異值的選擇準(zhǔn)則對(duì)信號(hào)進(jìn)行重構(gòu)，應(yīng)用重構(gòu)信號(hào)的傅里葉譜來檢測高速列車萬向軸的動(dòng)不平衡，消除EMD分解模式頻率混疊帶來頻譜雜亂，提高了譜的清晰度，凸顯了故障特征。應(yīng)用萬向軸動(dòng)不平衡試驗(yàn)數(shù)據(jù)對(duì)該方法進(jìn)行試驗(yàn)驗(yàn)證。

1 EMD分解及其等效濾波特性分析

1.1 振動(dòng)信號(hào)的EMD分解

經(jīng)驗(yàn)?zāi)Ｊ椒纸馐且环N新的非平穩(wěn)信號(hào)分析方法，將萬向節(jié)安裝機(jī)座振動(dòng)信號(hào)的這種非線性非平穩(wěn)信號(hào)分解為若干個(gè)基本模式分量(IMF)和一個(gè)余項(xiàng)?；灸Ｊ椒至勘仨殱M足兩個(gè)基本條件：條件1，在數(shù)據(jù)列中，數(shù)據(jù)的極值點(diǎn)數(shù)量和過零點(diǎn)數(shù)量必須相等或最多相差一個(gè)；條件2，任何一點(diǎn)，極大值和極小值的包絡(luò)的平均值為零。經(jīng)驗(yàn)?zāi)Ｊ椒纸獾脑敿?xì)過程如下：

第一步：給定被分解的振動(dòng)信號(hào)x(k)，k=1,2,…N，N表示信號(hào)的長度。

第二步：計(jì)算信號(hào)的極大值點(diǎn)xmax(k)和極大值坐標(biāo)kmax(k)，k=1,2,…Nm，Nm表示出去端點(diǎn)外的極大值點(diǎn)的個(gè)數(shù)。

第三步：采用線性插值外推糾正信號(hào)的極大值點(diǎn)的端點(diǎn)效應(yīng)，修正原理具體如圖1的(a)和(b)所示。

對(duì)圖1中所示的接近端點(diǎn)的極值點(diǎn)進(jìn)行線性插值，則有過點(diǎn)1和點(diǎn)2的直線方程為

(1)

則圖1中點(diǎn)3的縱坐標(biāo)為

(2)

圖1 線性插值外推控制端點(diǎn)效應(yīng)的基本原理Fig.1The basic principle of correcting end effect applying linear interpolation

如果x3≥x4，x3為端點(diǎn)坐標(biāo)，否則，x4為端點(diǎn)坐標(biāo)，因此圖中紅色的點(diǎn)為被選端點(diǎn)，與第二步得到的最大點(diǎn)整合得到包含端點(diǎn)的所有最大極值點(diǎn)。

第四步：計(jì)算信號(hào)的極小值點(diǎn)xmin(k)和極小值坐標(biāo)kmin(k)，k=1,2,…,Ni，Ni表示出去端點(diǎn)外的極小值點(diǎn)的個(gè)數(shù)。

第五步：采用線性插值外推糾正信號(hào)的極小值點(diǎn)的端點(diǎn)效應(yīng)，修正原理具體如圖1的(c)和(d)所示。同理應(yīng)用式(1)和(2)來確定x3的值的大小。如果x3≤x4，x3為端點(diǎn)坐標(biāo)，否則，x4為端點(diǎn)坐標(biāo)，因此圖中紅色的點(diǎn)為被選端點(diǎn)，與第四步得到的最小點(diǎn)整合得到包含端點(diǎn)的所有最小極值點(diǎn)。

第六步：采用三次樣條對(duì)所有極大值點(diǎn)進(jìn)行插值得到上包絡(luò)，采用三次樣條對(duì)所有極小值點(diǎn)進(jìn)行插值得到下包絡(luò)，得到上下包絡(luò)的平均值m1。則有

x(t)-m1=h1

第七步：判斷h1是否滿足IMF的條件，如果不滿足，對(duì)h1重復(fù)第二步到第六步得到包絡(luò)的平均值m11，在判斷h11=h1-m11是否滿足IMF的條件，如不滿足則循環(huán)s次，使得h1s=h1(s-1)-m1s滿足IMF的條件，記c1=h1s，則c1為信號(hào)x(t)的第一個(gè)滿足IMF條件的分量。

第八步：把c1從x(t)中分離出來，得到r1=x(t)-c1，將r1作為原始數(shù)據(jù)重復(fù)第二步到第七步得到信號(hào)x(t)的第二個(gè)滿足IMF條件的分量c2。重復(fù)循環(huán)n次使得rn為單調(diào)函數(shù)，得信號(hào)x(t)分解為n個(gè)滿足IMF的分量，則有

(3)

因此，信號(hào)分解為n個(gè)基本模式分量IMF和一個(gè)殘余項(xiàng)，則有

1.2 EMD分解的等效濾波特性分析

設(shè)定10 000個(gè)獨(dú)立的高斯白噪聲數(shù)據(jù)(其中均值為零，方差為1)，對(duì)每個(gè)高斯白噪聲數(shù)據(jù)應(yīng)用式(3)進(jìn)行EMD分解得到不同的IMF分量，計(jì)算每個(gè)分量的傅里葉頻譜，最好計(jì)算10 000個(gè)信號(hào)對(duì)應(yīng)傅里葉譜的平均值并進(jìn)行幅值歸一化，計(jì)算結(jié)果為EMD的等效濾波特性，具體如圖2所示。

圖2 EMD分解的等效濾波特性Fig.2 The equivalent filter characteristics of EMD

從圖2中的EMD分解的等效濾波特性可知，EMD分解與小波變換一樣，不同IMF之間依然存在頻率重疊現(xiàn)象，導(dǎo)致不同的IMF分量信號(hào)的噪聲變大，分解信號(hào)的頻率成分更為復(fù)雜，頻率重疊導(dǎo)致頻譜雜亂，混淆和淹沒故障特征和固有振動(dòng)特性。這對(duì)萬向軸動(dòng)不平衡的檢測、故障判定和故障溯源帶來了困難。為了消除這一缺陷，將Hankel矩陣、奇異值變換理論引入到EMD分解中，提出一種融合EMD、Hankel矩陣、SVD的故障檢測方法。

2 Hankel矩陣的構(gòu)造

為了方便應(yīng)用信號(hào)奇異值分解方法，先對(duì)EMD分解得到的IMF來構(gòu)造一個(gè)Hankel矩陣，設(shè)定近似信號(hào)c(i),i=1,3,…Nc，Nc表示信號(hào)的長度，利用此信號(hào)構(gòu)造的Hankel矩陣為公式(4)所示[10-11]

(4)

式中:1

3 奇異值分解理論及其對(duì)信號(hào)篩選特性分析

3.1 奇異值分解理論

奇異值分解是一種正交化分解方法，對(duì)于任意實(shí)矩陣A∈Rm×n，無論矩陣行向量或者列向量是否相關(guān)，必定存在正交矩陣U∈Rm×m和正交矩陣V∈Rn×n使得

A=UDVT

(5)

式中:D=(diag(σ1,σ2,…,σp),0)或者其轉(zhuǎn)置，0表示零矩陣，p=min(n,m)，并滿足σi≥σi+1

式(5)叫做奇異值分解[10]，σi就是分解得到的奇異值。

3.2 奇異值分解對(duì)信號(hào)的篩選特性分析

假設(shè)分解信號(hào)IMF由兩個(gè)部分組成，一是有用信號(hào)，二是測量、頻率重疊帶來的信號(hào)噪聲

c(i)=s(i)+ξ(i)

(6)

則有由c(i)構(gòu)造的矩陣可以分解為分別由s(i)和ξ(i)構(gòu)造的Hankel矩陣As和Aξ，即為

Ac=As+Aξ

(7)

式中:ξ(i)表示均值為0，方差為1的白噪聲

計(jì)算式(7)所描述信號(hào)的奇異值如圖3所示。

圖3 含噪聲不同頻率特性信號(hào)的奇異值分布規(guī)律Fig.3 Singular value distribution regularity of signal with noise

(8)

計(jì)算式 (8)所描述信號(hào)的奇異值如圖4所示。

圖4 不含噪聲不同頻率特性信號(hào)的奇異值分布規(guī)律Fig.4 Singular value distribution regularity of signal without noise

從圖4中可以看出，信號(hào)的行相關(guān)性大，矩陣呈病態(tài)，只有前面幾個(gè)關(guān)鍵特征值，其余特征值為零，與前面理論分析一致，同時(shí)針對(duì)不同頻率的合成信號(hào)的奇異值正好是單一信號(hào)的移位組合，這說明不同頻率信號(hào)具有的奇異值錯(cuò)位特性，奇異值分布規(guī)律上呈現(xiàn)疊層結(jié)構(gòu)，這正是消除頻率重疊噪聲的關(guān)鍵理論基礎(chǔ)；從圖3中可以看出，噪聲的奇異值分布均勻，矩陣呈現(xiàn)良性特性，同時(shí)不同頻率信號(hào)也表現(xiàn)出奇異值的錯(cuò)位特性。為此，本文將EMD與、Hankel矩陣、奇異值理論相結(jié)合，利用奇異值對(duì)不同頻率信號(hào)的錯(cuò)位特性來消除EMD分解的頻率重疊的缺陷，提高譜的清晰度和對(duì)故障的表現(xiàn)力。

4 萬向軸動(dòng)不平衡EMD-Hankel-SVD檢測模型

依據(jù)不平衡故障頻率可能存在的范圍，對(duì)信號(hào)進(jìn)行EMD分解，對(duì)單個(gè)IMF信構(gòu)造Hankel矩陣，應(yīng)用奇異值理論對(duì)Hankel矩陣進(jìn)行奇異值分解，在奇異值的階梯結(jié)構(gòu)中選擇貢獻(xiàn)量較大的前3階奇異值重構(gòu)信號(hào)，用重構(gòu)信號(hào)的傅里葉變換來檢測萬向軸的動(dòng)不平衡。檢測模型如圖5所示。

圖5 萬向軸動(dòng)不平衡EMD-Hankel-SVD檢測模型Fig.5 The detection model of dynamic imbalance with cardan shaft applying the EMD-Hankel-SVD

從圖中5可以看出，萬向軸動(dòng)不平衡檢測模型主要包括三個(gè)關(guān)鍵步驟：第一步，對(duì)振動(dòng)加速度信號(hào)EMD分解得到5個(gè)IMF分量；第二步，構(gòu)造每個(gè)分量的Hankel矩陣；第三步，對(duì)每個(gè)Hankel矩陣進(jìn)行奇異值分解與重構(gòu)；第四步，對(duì)重構(gòu)信號(hào)做FFT變換，并合并變換信號(hào)的頻率軸。合并信號(hào)的傅里葉譜來檢測萬向軸的動(dòng)不平衡。

5 萬向傳動(dòng)軸動(dòng)不平衡檢測模型的試驗(yàn)驗(yàn)證

為了驗(yàn)證本文提出的萬向軸動(dòng)不平衡檢測的EMD-Hankel-SVD方法和模型的有效性，搭建了如圖6所示的萬向軸動(dòng)不平衡試驗(yàn)臺(tái)。

圖6 萬向傳動(dòng)軸動(dòng)不平衡試驗(yàn)平臺(tái)Fig.6 Test platform of cardan shaft dynamic unbalance

該試驗(yàn)平臺(tái)的動(dòng)力傳遞途徑為：電機(jī)、齒輪箱、萬向節(jié)、萬向傳動(dòng)軸、萬向節(jié)。不平衡軸選用專項(xiàng)修軸(不平衡量超過車輛的使用標(biāo)準(zhǔn))和新軸(不存在不平衡量)兩種。采集萬向軸靠近電機(jī)端的萬向節(jié)軸的垂向振動(dòng)加速度信號(hào)作為萬向軸不平衡檢測的信號(hào)源。萬向軸的試驗(yàn)轉(zhuǎn)速為2 750 r/min，對(duì)存在動(dòng)不平衡的專修軸做故障試驗(yàn)，采集數(shù)據(jù)如圖7所示。

采用EMD對(duì)信號(hào)分解得到IMF分量信號(hào)分別如圖8～圖12所示。

“我什么意思你不知道？好吧，讓我告訴你，憑我對(duì)張仲平的了解，他絕對(duì)不會(huì)站在你現(xiàn)在站的這個(gè)位置上。年輕人，我給你的忠告是這樣，人在高處，別兩邊沒有依靠。得防著有人從你背后下手?！?/p>

圖7 專修軸的不平衡試驗(yàn)數(shù)據(jù)Fig．7Unbalanceexperimentdataofspecialrevisedshaft圖8 IMF1分量信號(hào)Fig．8ComponentsignalofIMF1圖9 IMF2分量信號(hào)Fig．9ComponentsignalofIMF2

圖10 IMF3分量信號(hào)Fig．10ComponentsignalofIMF3圖11 IMF4分量信號(hào)Fig．11ComponentsignalofIMF4圖12 IMF5分量信號(hào)Fig．12ComponentsignalofIMF5

圖13 IMF1的SVD分解重構(gòu)信號(hào)Fig．13IMF1reconstructionsignalthroughSVD圖14 IMF2的SVD分解重構(gòu)信號(hào)Fig．14IMF2reconstructionsignalthroughSVD圖15 IMF3的SVD分解重構(gòu)信號(hào)Fig．15IMF3reconstructionsignalthroughSVD

圖16 IMF4的SVD分解重構(gòu)信號(hào)Fig.16 IMF4reconstruction signal through SVD

圖17 IMF5的SVD分解重構(gòu)信號(hào)Fig.17 IMF5 reconstruction signal through SVD

分別對(duì)信號(hào)分量IMF1到信號(hào)分量IMF5構(gòu)造Hankel矩陣，對(duì)矩陣進(jìn)行奇異值分解并進(jìn)行重構(gòu)，得到各個(gè)分量的重構(gòu)信號(hào)分別如圖13～圖17所示。

分別對(duì)信號(hào)IMF1～I(xiàn)MF5的SVD分解重構(gòu)信號(hào)進(jìn)行傅里葉變換,并合并其頻率軸,具體如圖18所示。

圖18中的關(guān)鍵固有頻率結(jié)構(gòu)如表1所示。

表1 專修軸的主要故障頻率

結(jié)合萬向軸的試驗(yàn)轉(zhuǎn)速可知，第1頻率45.90 Hz是萬向軸的試驗(yàn)轉(zhuǎn)頻(2750/60=45.83), 第3頻率137．70正好是轉(zhuǎn)頻的3倍頻。因此通過EMD分解、Hankel矩陣和奇異值分解的萬向軸動(dòng)不平衡檢測方法能夠有效檢測不平衡故障引起的轉(zhuǎn)速基頻、倍頻等故障特征。為萬向軸的不平衡故障檢測提供了有效的檢測方法。

空心軸的一階固有頻率通用的計(jì)算公式得到萬向軸的一階固有頻率

(10)

式中：da為萬向軸的外徑，大小為142 mm；db為萬向軸的內(nèi)徑，大小為130 mm；Lc為萬向軸長度，大小為1830 mm；η為空心軸系數(shù)，取0.8。

則有

f=93.5 Hz

圖18中和表1中所示出的第2個(gè)頻率與萬向軸的一階固有頻率十分接近，第2個(gè)頻率恰好為一階固有頻率的2倍頻。因此該方法有效的提取了萬向軸的固有振動(dòng)特性。

在同樣的試驗(yàn)條件下，對(duì)不存在動(dòng)不平衡的萬向軸新軸做故障試驗(yàn)，采集的試驗(yàn)數(shù)據(jù)如圖19所示。對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行EMD-Hankel-SVD處理后得到的傅里葉如圖20所示(鑒于篇幅，不詳細(xì)列舉EMD、SVD分解的圖形)。

圖18 專修軸振動(dòng)數(shù)據(jù)EMD?Hankel?SVD處理后的傅里葉譜Fig．18FourierspectrumofEMD?Hankel?SVDinthespecialrevisedshaft圖19 新軸的不平衡試驗(yàn)數(shù)據(jù)Fig．19Unbalanceexperimentdataofnewcardanshaft圖20 新軸振動(dòng)數(shù)據(jù)EMD?Hankel?SVD處理后的傅里葉譜Fig．20FourierspectrumofEMD?Hankel?SVDinthenewshaft

圖20中的關(guān)鍵固有頻率結(jié)構(gòu)如表2所示。

表2 新軸的主要故障頻率

從表2中可見，沒有出現(xiàn)反映萬向軸的不平衡帶來的故障頻率，僅在一階固有頻率附近、0.5倍固有頻率附近、1.5倍固有頻率附近處出現(xiàn)振動(dòng)，說明該檢測方法虛警幾率低，具有較高的可靠性。

6 與純EMD方法在萬向軸動(dòng)不平衡檢測中的應(yīng)用比較

直接應(yīng)用EMD分解得到5個(gè)IMF信號(hào)分量，對(duì)每個(gè)分量進(jìn)行傅里葉變換并合并其頻率軸，專修軸數(shù)據(jù)處理結(jié)果如圖21所示。

圖21 專修軸試驗(yàn)數(shù)據(jù)純EMD處理結(jié)果Fig.21 Processing of special revised shaft experiment data by EMD

從圖21中可以看出，單純EMD分解模式能夠提取萬向軸的不平衡振動(dòng)的基頻、一階固有頻率、一階固有頻率的2倍頻，原因是EMD本身是一種優(yōu)良的非平穩(wěn)信號(hào)處理方法。但很難提取出不平衡振動(dòng)倍頻的微弱振動(dòng)特征，頻率點(diǎn)137.7 Hz～139.6 Hz處無能量峰值，是一條直線，原因是微弱特征被頻率重疊噪聲淹沒了。新軸數(shù)據(jù)處理結(jié)果如圖22所示。

圖22 新軸試驗(yàn)數(shù)據(jù)純EMD處理結(jié)果Fig.22 The processing of the new shaft experiment data by EMD

從圖22中可以看出，單純EMD分解模式能夠提取萬向軸的一階固有頻率、0.5倍固有頻率、1.5倍固有頻率，但是頻率雜亂，在軸的旋轉(zhuǎn)頻率附近依然存在一定能量的振動(dòng),容易誤判萬向軸的動(dòng)不平衡故障，造成虛警。

因此與EMD方法相比，EMD-Hankel-SVD檢測方法能夠提取微弱的不平衡特征，同時(shí)譜的清晰度顯著提高，除了能夠完整提取萬向軸的動(dòng)不平衡的固有特征及其倍頻外，也能夠提取萬向軸承一階固有頻率及其分頻、倍頻的振動(dòng)特征，提高診斷的準(zhǔn)確性，減少了錯(cuò)誤和虛警。

7 結(jié) 論

本文將EMD分解、Hankel矩陣、SVD分解相結(jié)合，提出一種高速列車萬向軸故障檢測方法和模型。概括起來具有以下特點(diǎn)：

(1) EMD 分解與小波變換相比，是一種自適應(yīng)的信號(hào)分解方法，十分適合處理非平穩(wěn)信號(hào)，但是EMD分解與小波分解一樣，同樣存在相鄰IMF分量的頻率重疊混淆的缺陷，易致頻率結(jié)構(gòu)混雜，淹沒故障特征。

(2) 奇異值分解針對(duì)不同頻率結(jié)構(gòu)的信號(hào)在奇異值分布規(guī)律上呈現(xiàn)錯(cuò)位特性，為關(guān)鍵振動(dòng)模式提取的奇異值選擇提供了理論依據(jù)。

(3) EMD、Hankel、SVD的有機(jī)結(jié)合，能夠有效提取萬向軸動(dòng)不平衡的基頻、倍頻的故障特性，同時(shí)也能夠提取萬向軸的一階固有頻率、0.5倍一階固有頻率、1.5倍一階固有頻率的振動(dòng)特征，基本消去相鄰IMF頻率重疊對(duì)信號(hào)特征的混淆和淹沒效應(yīng)。

(4) 該方法與純EMD相比，其譜的分辨率、清晰度、故障表征力得到了顯著提高，對(duì)減少虛警漏檢、提高故障檢測品質(zhì)具有重要的工程意義。

[1] 冷揚(yáng)立，李秋澤，李慶國,等. CRH5 型動(dòng)車組萬向軸結(jié)構(gòu)及臨界轉(zhuǎn)速分析[J]. 鐵道車輛，2010，48(12)：6-11. LENG Yang-li，LI Qiu-ze，LI Qing-guo, et al. Analysis of the structure and critical rotating speed of the cardan Shafts for CRH5 multiple units[J]. Railway Vehicles， 2010，48(12)：6-11.

[2] 張紅軍，姚遠(yuǎn)，羅赟,等. CRH5 型動(dòng)車萬向軸傳動(dòng)系統(tǒng)技術(shù)特征分析[J].鐵道學(xué)報(bào)，20009，31(2)：115-119. ZHANG Hong-jun，YAO Yuan，LUO Yun, et al. Analysis on technical characteristics of CRH5 cardan drive system[J]. Journal of the China Railway Society， 20009，31(2)：115-119.

[3] 姚遠(yuǎn)，張紅軍，羅赟. CRH5型動(dòng)車萬向軸扭轉(zhuǎn)振動(dòng)分析[J].中國鐵道科學(xué)，2009，30(2)：82-86. YAO Yuan，ZhANG Hong-jun，LUO Yun. Analysis on the torsional vibration of cardan shaft on CRH5 motor Car[J]. China Railway Science，2009，30(2)：82-86.

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

[5] 張超, 陳建軍, 徐亞蘭. 基于EMD 分解和奇異值差分譜理論的軸承故障診斷方法[J].振動(dòng)工程學(xué)報(bào)，2011,24(5):539-545. ZHANG Chao, CHEN Jian-jun, XU Ya-lan. A bearing fault diagnosis method based on EMD and difference spectrum theory of singular value[J].Journal of Vibration Engineering，2011,24(5):539-545.

[6] Huang N E, WU Zhao-hua，Long S R, et al. On instantaneous frequency[J]. Advance in Adaptive Data Analysis， 2009，1(2)：177-229.

[7] 時(shí)培明,丁雪娟,李庚,等. 一種EMD 改進(jìn)方法及其在旋轉(zhuǎn)機(jī)械故障診斷中的應(yīng)用[J].振動(dòng)與沖擊，2013,32(4):185-190. SHI Pei-ming, DING Xue-juan, LI Geng, et al. An improved method of EMD and its applications in rotating machinery fault diagnosis[J].Journal of Vibration and Shock，2013,32(4):185-190.

[8] Cheng Jun-sheng, Yu De-jie, Yang Yu. The application of energy operator demodulation approach based on EMD in machinery fault diagnosis[J], Mechanical Systems and Signal Processing,2007,21:668-677.

[9] Hea X H, Huab X G, Chenb Z Q. EMD-based random decrement technique for modal parameter identification of an existing railway bridge[J]. Engineering Structs, 2011,33:1348-1356.

[10] 戈盧布G H，范洛恩 C F. 矩陣計(jì)算[M]. 袁亞湘，譯.北京：科學(xué)出版社，2001.

[11] 趙學(xué)智，葉邦彥，陳統(tǒng)堅(jiān).奇異值差分譜理論及其在車床主軸箱故障診斷中的應(yīng)用[J]. 機(jī)械工程報(bào)，2010，46(1)：100-108. ZHAO Xue-zhi，YE Bang-yan， CHEN Tong-jian. Difference spectrum theory of singular value and its application to the fault diagnosis of headstock of lathe[J]. Journal of Mechanical Engineering， 2010，46(1)：100-108.

Detection of dynamic imbalance due to cardan shaft in high-speed train based on EMD-Hankel-SVD method

DING Jian-ming1, WANG Han2, LIN Jian-hui1, HUANG Chen-guang1

(1. Southwest Jiaotong University, State Key Laboratory of Traction Power, Chengdu 610031, China;2. CSR Qingdao Sifang Co. Ltd., High-Speed Train System Integration State Engineering Laboratory，Qingdao 266111, China)

A new method of detecting dynamic imbalance due to cardan shaft in high-speed train was proposed applying a synthesized method of EMD (Empirical Modal Decomposition), Hankel matrix and SVD (singular value decomposition) against the aliasing defect between adjacent intrinsic model functions existing in EMD. The vibration signals of gimbal installed base were decomposed through EMD to get the different intrinsic model components. The Hankel matrix, constructed of the single decomposition model components, was orthogonally decomposed through SVD and the key singular values were selected to reconstruct the vibration signal on the base of the key stack of singular values. The Fourier spectrum of the reconstructed signal was applied to detect dynamic imbalance of shaft and the clutter spectrum caused by the aliasing defect between adjacent intrinsic model functions was eliminated so as to highlight the failure characteristics. The method was verified by test data in the condition of dynamic imbalance. The results show the method can effectively detect the fault vibration characteristics caused by cardan shaft dynamic imbalance and extract the natural vibration features. By comparison to the simple EMD, the clarity and failure characterization ability are significantly improved.

high-speed train; dynamic imbalance due to cardan shaft; EMD; Hankel matrix; SVD; dynamic detection

國家自然科學(xué)基金(61134002;51305358);863計(jì)劃(2011AA110501);精密測試技術(shù)及儀器國家重點(diǎn)實(shí)驗(yàn)室開放課題(PIL1303)

2014-01-08 修改稿收到日期：2014-05-13

丁建明 男，博士，助理研究員，1981年生

U211.5；U270.331.1

A

10.13465/j.cnki.jvs.2015.09.030