董 偉,黃 文,邢宗義,龍 靜
(1.南京理工大學(xué)機(jī)械工程學(xué)院,南京 210094; 2.南京理工大學(xué)自動(dòng)化學(xué)院,南京 210094; 3.廣州地鐵集團(tuán)有限公司運(yùn)營(yíng)總部,廣州 510380)
基于小波包能量熵的鋼軌波磨檢測(cè)方法
董 偉1,黃 文1,邢宗義2,龍 靜3
(1.南京理工大學(xué)機(jī)械工程學(xué)院,南京 210094; 2.南京理工大學(xué)自動(dòng)化學(xué)院,南京 210094; 3.廣州地鐵集團(tuán)有限公司運(yùn)營(yíng)總部,廣州 510380)
針對(duì)目前軌道波浪型磨損檢測(cè)效率慢、成功率低等問(wèn)題,提出一種基于小波包能量熵的軌道波磨檢測(cè)技術(shù)。首先建立軌道波磨簡(jiǎn)化模型及車(chē)輛集總化參數(shù)模型并搭建SIMPACK仿真動(dòng)力學(xué)模型,其次創(chuàng)新性地將小波包和能量熵相結(jié)合應(yīng)用在軌道波磨檢測(cè)領(lǐng)域,再通過(guò)仿真模型得到不同波磨分組下的軸箱振動(dòng)響應(yīng),然后對(duì)振動(dòng)加速度信號(hào)進(jìn)行4層小波包分解并計(jì)算各節(jié)點(diǎn)的小波包能量和小波包能量熵,最后通過(guò)對(duì)比分析仿真結(jié)果確定軌道波磨故障類(lèi)型。
鋼軌;波浪型磨損;故障診斷;軸箱振動(dòng)加速度;小波包;能量熵
軌道作為承載列車(chē)運(yùn)行的關(guān)鍵部件,其質(zhì)量情況影響著乘客乘坐的舒適性以及列車(chē)運(yùn)營(yíng)的安全性。而在日常的運(yùn)營(yíng)過(guò)程中,軌道不斷接受來(lái)自輪對(duì)的沖擊和壓迫,這使得軌道會(huì)出現(xiàn)橫向和垂向的形變,嚴(yán)重的形變會(huì)對(duì)列車(chē)運(yùn)營(yíng)造成巨大的危害。因此對(duì)軌道質(zhì)量狀態(tài)進(jìn)行監(jiān)測(cè),保障軌道安全可靠是地鐵公司的檢修項(xiàng)目之一[1]。
目前各國(guó)均將基于運(yùn)營(yíng)車(chē)輛的軌道檢測(cè)技術(shù)和裝備研制作為研究的重點(diǎn)。Grassie[2]通過(guò)對(duì)運(yùn)營(yíng)車(chē)輛軸箱振動(dòng)以及其他部件振動(dòng)信號(hào)進(jìn)行分析,并在1996年最早提出基于運(yùn)營(yíng)車(chē)輛的軌道動(dòng)態(tài)檢測(cè)構(gòu)想;HitoshiTsunashima[3]等通過(guò)采集車(chē)體振動(dòng)信號(hào)并進(jìn)行小波包分析,發(fā)現(xiàn)車(chē)體振動(dòng)信號(hào)在某個(gè)頻域上能夠完整表征出軌道波磨故障;曹西寧[4]等通過(guò)對(duì)采集的軸箱振動(dòng)信號(hào)進(jìn)行希爾伯特黃變換分析,發(fā)現(xiàn)該方法能夠?qū)壍啦黄巾樳M(jìn)行大致的定位;李宏峰[5]等采用SIMPACK軟件進(jìn)行列車(chē)仿真建模,對(duì)仿真振動(dòng)信號(hào)進(jìn)行二次積分發(fā)現(xiàn)軸箱振動(dòng)信號(hào)能夠很好地表征軌道波磨故障;史紅梅[6]等研究了基于運(yùn)營(yíng)列車(chē)振動(dòng)信號(hào)的軌道剛度不平順和幾何不平順的感知算法,實(shí)驗(yàn)結(jié)果表明該方法能有效識(shí)別和區(qū)分軌道故障。
本文提出了基于小波包能量熵的鋼軌波磨檢測(cè),通過(guò)對(duì)軸箱振動(dòng)加速度信號(hào)進(jìn)行小波包能量熵分析,利用SIMPACK對(duì)搭建好的模型進(jìn)行仿真,仿真結(jié)果表明小波包能量熵能夠很好地確定鋼軌波磨故障。
軌道波浪形磨損是一種周期性的垂向不平順,因此可以用諧波函數(shù)來(lái)近似描述軌道波磨??紤]到波磨現(xiàn)象的多樣性,若只考慮一種波長(zhǎng)的波磨,即可用單諧波函數(shù)來(lái)描述,如圖1(a)所示;若考慮多種波長(zhǎng)符合的軌道波磨,則可用多諧波函數(shù)描述,如圖1(b)所示。
圖1 諧波激擾簡(jiǎn)化模型
對(duì)于單諧波函數(shù),其表達(dá)式為
對(duì)于多諧波函數(shù),其表達(dá)式為
式中,L為波磨波長(zhǎng);a為波磨波深;ω=2πv/L;n為激擾波數(shù)。
為了便于分析軌道波磨對(duì)軸箱引起的振動(dòng)變化,需要建立列車(chē)集總化參數(shù)模型??紤]到簧上部分引起的軸箱振動(dòng)變化遠(yuǎn)沒(méi)有輪軌部分引起的振動(dòng)沖擊明顯[7],因此著重對(duì)簧下質(zhì)量部分進(jìn)行建模。建立的模型如圖2所示。
圖2 列車(chē)集總化模型
圖2中,將輪對(duì)上部統(tǒng)一稱為車(chē)體部分,并考慮用有剛度有阻尼的二系簧相連接。輪軌作用部分轉(zhuǎn)化為具有一定質(zhì)量的輪對(duì)通過(guò)具有剛度的一系簧和軌道相連接,其中軌道是一個(gè)有剛度有阻尼的單自由度系統(tǒng)。圖中v為車(chē)速,m1為簧下質(zhì)量,z1為靜態(tài)下的輪對(duì)垂向位移,K1為一系簧的垂向剛度,m2為軌道質(zhì)量,K2為軌道垂向剛度,c2為軌道阻尼,z2為靜態(tài)時(shí)軌道垂向位移,所有參數(shù)向下為正方向。
輪軌的垂向動(dòng)力學(xué)方程在車(chē)速v下為
假設(shè)軌道波磨為η,則軌道的靜平衡位移變?yōu)閆2+η,輪軌的垂向動(dòng)力學(xué)方程在車(chē)速v下為
將式(4)進(jìn)行常微分方程求解,可得
以單諧波波形為例,對(duì)式(1)進(jìn)行拉氏變換有
式中,z1(s)和z2(s)為z1、z2的拉氏變換;HZ1(s)和HZ2(s)分別為z1、z2的傳遞函數(shù)。
將式(5)進(jìn)行拉普拉斯變換,并將式(7)代入計(jì)算得
聯(lián)立以上各式有
HZ1(s)=
式中,Ai,B,C,pi為常數(shù),其大小由ω1、ω2、K1、K2確定。
對(duì)式(11)做拉氏反變換得
特別地,當(dāng)列車(chē)以車(chē)速v行駛,軌道存在波長(zhǎng)為λ的波磨故障,則有波磨的頻率為ω=2πv/λ,而軸箱的響應(yīng)頻率為f=v/λ,因此只要知道車(chē)速和軸箱響應(yīng)頻率即可求出軌道波磨波長(zhǎng)。
SIMPACK軟件是由德國(guó)航空宇航局于1985年開(kāi)發(fā),由于其使用了基于相對(duì)坐標(biāo)系遞歸的優(yōu)秀算法,使得其相對(duì)于其他動(dòng)力學(xué)軟件計(jì)算方程更少,速度更快[8]。本文使用的是SIMPACK Wheel/Rail模塊,即鐵路行業(yè)專業(yè)模塊。
SIMPACK建立軌道車(chē)輛模型首先要做的就是構(gòu)建模型拓?fù)鋱D,拓?fù)鋱D構(gòu)建的好壞直接影響建模的好壞、后期計(jì)算速度和仿真結(jié)果的正確性。SIMPACK中涉及幾個(gè)建模的重要概念,即:體(Bodies)、約束(Cons-traints)、鉸(Joints)、力元(Force Elements)、標(biāo)記(Makers)等,通常的軌道車(chē)輛建模就是先利用以上一些參數(shù)進(jìn)行拓?fù)鋱D的構(gòu)建,再創(chuàng)建具體模型。考慮到列車(chē)部件的重復(fù)性和復(fù)雜性,采用自底向上的結(jié)構(gòu)化建模方式,將車(chē)體模型按照輪對(duì)-轉(zhuǎn)向架-車(chē)輛三步進(jìn)行建模,軌道模型采用剛體模型[9]。其拓?fù)浣Y(jié)構(gòu)如圖3所示。
圖3 軌道車(chē)輛拓?fù)鋱D
由圖3可看出,將整節(jié)車(chē)廂分解成車(chē)體和子結(jié)構(gòu)兩大塊,其中子結(jié)構(gòu)中封裝了虛車(chē)體、轉(zhuǎn)向架構(gòu)架、輪對(duì)。車(chē)體和子結(jié)構(gòu)之間采用二系懸掛力元連接;子結(jié)構(gòu)內(nèi)部轉(zhuǎn)向架構(gòu)架和輪對(duì)之間采用一系懸掛力元連接,并且構(gòu)架上輔以相應(yīng)鉸接(07號(hào));輪對(duì)和軌道通過(guò)09號(hào)約束、07號(hào)鉸接以及89力元相連接。建立拓?fù)鋱D后,依據(jù)軌道車(chē)輛參數(shù)進(jìn)行三維實(shí)體建模。
SIMPACK建模的過(guò)程是參數(shù)化的,必須依據(jù)實(shí)際軌道車(chē)輛參數(shù)進(jìn)行相應(yīng)部件參數(shù)設(shè)置。選用某型軌道車(chē)輛動(dòng)力學(xué)參數(shù)建立如圖4所示單節(jié)車(chē)輛模型。
圖4 單節(jié)車(chē)輛模型
多分辨率分析(Multi-Resolution Analysis)[10]是由S.Mallat于1988年提出的概念,它的核心思想就是不斷對(duì)信號(hào)進(jìn)行降采樣處理從而得到細(xì)化的信號(hào)特征。關(guān)于多分辨率分析可以結(jié)合圖5用三層結(jié)構(gòu)樹(shù)進(jìn)行說(shuō)明。圖中S為原始信號(hào),Ai和Di(i=1,2,3)分別為相應(yīng)信號(hào)段分解后的低頻信號(hào)和高頻信號(hào)。由圖5可看出一個(gè)原始信號(hào)在多分辨率分析中只對(duì)低頻部分做了持續(xù)分解,而信號(hào)的高頻部分在第一層分解后就不再進(jìn)行任何處理。一個(gè)信號(hào)經(jīng)過(guò)三層多分辨率分析后的關(guān)系可表示為:S=A3+D3+D2+D1。
圖5 三層多分辨率分析結(jié)構(gòu)樹(shù)
雖然多分辨率分析的出現(xiàn)是時(shí)頻信號(hào)處理的一次革命,但由于它自身算法對(duì)高頻信號(hào)處理的忽視性,使得它先天就存在對(duì)高頻信號(hào)處理能力的不足。因此針對(duì)這一缺陷,M.V.Wickerhauser和R.Rcoifman在小波分析和多分辨率分析的基礎(chǔ)上提出了小波包分析。
小波包分析[11](Wavelet Packet Analysis)在原有多尺度分析的基礎(chǔ)上進(jìn)一步優(yōu)化,將多分辨率分析中忽略的高頻信號(hào)也做同低頻信號(hào)一樣的逐層細(xì)化分解,即小波包分解能夠?qū)π盘?hào)低頻和高頻部分同時(shí)進(jìn)行多尺度分析,提供多尺度分解[12,13]。原始信號(hào)經(jīng)過(guò)三層小波包分解后的結(jié)構(gòu)樹(shù)[14]如圖6所示。
圖6 三層小波包分解結(jié)構(gòu)樹(shù)
在信息論中,熵被用來(lái)衡量一個(gè)隨機(jī)變量出現(xiàn)的期望值[15]。信息熵理論是香農(nóng)于1948年提出,他認(rèn)為一個(gè)信號(hào)的信息熵越大,則它出現(xiàn)的各種可能性就越高[16]。根據(jù)信息熵理論:假設(shè)有一隨機(jī)變量X={x1,x2,x3…xn}組成的信號(hào),信息熵的計(jì)算公式如下
式中,H(X)為信息熵測(cè)度;pi為取值概率。
由于信息熵能夠反映信號(hào)出現(xiàn)隨機(jī)變化的程度,對(duì)于某個(gè)信號(hào)源,其信息熵越大即表示該信號(hào)在采樣范圍內(nèi)的幅值波動(dòng)越大,相反信息熵越小該信號(hào)越趨于平穩(wěn)。因此一段信號(hào)的平均統(tǒng)計(jì)特征大小通??梢杂眠@段信號(hào)的信息熵大小來(lái)表示。因此其在故障診斷領(lǐng)域有著重要的作用。
假設(shè)有一待分析信號(hào),其信號(hào)長(zhǎng)度為N,該段信號(hào)的小波包能量熵可由以下步驟求得。
(1)對(duì)待測(cè)信號(hào)進(jìn)行j層小波包分解,經(jīng)過(guò)分解后得到一系列不同頻段的子信號(hào)序列。
(2)將得到的高低頻序列的小波包重構(gòu)系數(shù)表示為Sjk,(k=0,1,…2j-1)。
(3)計(jì)算每個(gè)子序列的小波包能量值Ejk,(k=0,1,…2j-1)。
式中,Ai(t)為該節(jié)點(diǎn)處的幅值最大值;ti-1和ti分別是該節(jié)點(diǎn)信號(hào)的起止時(shí)間。
(4)求取總能量值E
(5)計(jì)算每個(gè)子節(jié)點(diǎn)能量相對(duì)于總能量的概率pjk
(6)求取小波包能量熵值測(cè)度
式中,Hjk為信號(hào)j層分解后的第k個(gè)小波包能量熵測(cè)度。
基于小波包能量熵的軌道波磨故障檢測(cè)算法流程如圖7所示,其主要步驟為:
(1)對(duì)采集的軸箱垂向振動(dòng)信號(hào)進(jìn)行濾波處理,濾除噪聲和其他干擾項(xiàng),只保留和軌道波磨相關(guān)的振動(dòng)信息;
(2)對(duì)濾波后的軸箱垂向振動(dòng)信號(hào)進(jìn)行j層小波包分解,使信號(hào)分解成2j個(gè)包含不同頻率帶的子節(jié)點(diǎn)信號(hào);
(3)計(jì)算分解后信號(hào)的小波包能量,這其中包括各節(jié)點(diǎn)的小波包能量以及信號(hào)總的能量值;
(4)依據(jù)式(16)和式(17)計(jì)算各節(jié)點(diǎn)的小波包能量熵Hjk;
(5)計(jì)算分解信號(hào)的總熵H并和閾值進(jìn)行比較,若總熵超出閾值,則判定該段軸箱振動(dòng)信號(hào)包含軌道波磨故障,相反則不存在軌道波磨故障;
(6)若信號(hào)存在軌道波磨故障,對(duì)小波包分解信號(hào)進(jìn)行頻率分析,找出時(shí)域上的頻率高集中區(qū)域,計(jì)算出軌道波磨在軌道上的具體故障位置以及根據(jù)頻率集中帶計(jì)算出波磨故障的波長(zhǎng)范圍。
圖7 小波包能量熵檢測(cè)軌道波磨算法流程
在實(shí)際工程應(yīng)用中,由于背景噪聲的存在,必須對(duì)軸箱振動(dòng)信號(hào)進(jìn)行濾波處理,基于SIMPACK軟件仿真出的軸箱振動(dòng)信號(hào)并沒(méi)有受到噪聲信號(hào)干擾,所以算法流程中的濾波處理可以暫不考慮。
國(guó)內(nèi)外鐵路行業(yè)習(xí)慣使用功率譜密度函數(shù)(power spectral density)來(lái)表示正常軌道的不平順程度變化[17]。它是以頻率為橫坐標(biāo),PSD為縱坐標(biāo)來(lái)表示的一條曲線。目前我國(guó)學(xué)者普遍認(rèn)為美國(guó)五級(jí)軌道譜較為符合我國(guó)現(xiàn)有軌道情況。本文采用在軌道譜上疊加波磨數(shù)據(jù)的方法進(jìn)行軌道波磨設(shè)置。
為了更為全面地驗(yàn)證檢測(cè)算法的可行性,選取14組波長(zhǎng)不同的波磨疊加五級(jí)軌道譜垂向信號(hào)進(jìn)行仿真實(shí)驗(yàn)。SIMPACK仿真實(shí)驗(yàn)時(shí),設(shè)定軌道長(zhǎng)度為95 m,采樣頻率為1 kHz,在0~50 m處左右軌均為不做任何處理的軌道,50~95 m處右側(cè)軌道分別設(shè)定14組不同的波磨故障信號(hào),左側(cè)軌道只添加美國(guó)五級(jí)軌道譜垂向不平順,具體設(shè)置參數(shù)如表1所示。
表1 軌道波磨分組設(shè)置參數(shù)
通過(guò)對(duì)上面分組線路得到的仿真軸箱振動(dòng)信號(hào)進(jìn)行數(shù)據(jù)分析,從兩個(gè)方向進(jìn)行軌道波磨故障診斷。首先是對(duì)軸箱振動(dòng)信號(hào)進(jìn)行小波包能量熵計(jì)算,通過(guò)閾值判斷軌道是否存在波磨故障,然后對(duì)故障軌道進(jìn)行小波包時(shí)頻圖分析,定位軌道波磨故障時(shí)域位置及確定波磨波長(zhǎng)范圍。
(1)軌道波磨檢測(cè)閾值的設(shè)定與故障判斷
根據(jù)圖7軸箱振動(dòng)信號(hào)進(jìn)行小波包分解,選用db4小波作為小波母函數(shù),分解層數(shù)暫定為4層,計(jì)算得到各節(jié)點(diǎn)的小波包能量熵后,根據(jù)式(17)計(jì)算4層小波包分解后的總熵值,計(jì)算結(jié)果見(jiàn)表2。
表2 各分組實(shí)驗(yàn)總熵計(jì)算結(jié)果
對(duì)表2進(jìn)行分析可看出,第一組軌道僅包含美國(guó)五級(jí)軌道譜,其小波包能量熵為1.835 6×10-5,而其余13組包含波磨的軌道小波包能量熵值范圍在0.3~2.5,由此可看出,包含軌道波磨段的軌道其小波包能量熵遠(yuǎn)大于無(wú)故障軌道;進(jìn)一步分析波磨波長(zhǎng)相同但波深不同的軌道,例如當(dāng)波長(zhǎng)為150 mm,波深分別為0.05、0.1、0.15 mm時(shí),小波包能量分別為1.996 4、2.310 5和2.345 0,小波包能量熵基本和波深呈正相關(guān)關(guān)系,其余波長(zhǎng)下也有相同情況;交叉分析不同波長(zhǎng),相同波深情況下的小波包能量熵,例如當(dāng)波長(zhǎng)分別為150、200、250 mm,波深為0.10 mm時(shí),小波包能量熵值分別為:2.310 5、1.833 4和1.497 7,可發(fā)現(xiàn)在波深相同情況下波長(zhǎng)越小,小波包能量熵越大,且其余波深下也會(huì)有類(lèi)似情況;對(duì)比第9組、第11組和第12組數(shù)據(jù),相同波磨故障狀況,不同列車(chē)運(yùn)行速度,在速度分別為10、15,20 m/s情況下,小波包能量熵分別為1.497 7、1.773 6和1.794 8,可看出小波包能量熵隨著速度的增加而增加。對(duì)比第13組和第14組實(shí)驗(yàn),當(dāng)有兩組不同波長(zhǎng)波磨分別疊加軌道不同位置時(shí),疊加波磨波長(zhǎng)越大,小波包能量越小。
綜合以上分析結(jié)果可以得出,小波包能量熵在軌道波磨故障檢測(cè)上有著明顯的度量性,其值的大小受波磨故障的波長(zhǎng)、波深以及行車(chē)速度影響,且與這三者均有著明顯的線性關(guān)系,因此可以將軸箱振動(dòng)信號(hào)的小波包能量熵值作為軌道波磨故障檢測(cè)的閾值參考。在本次仿真中綜合以上數(shù)據(jù)結(jié)果,建議將閾值設(shè)定為0.3,若計(jì)算得小波包能量熵值大于0.3,則判斷發(fā)生波磨故障,反之則無(wú)波磨故障。由于閾值的選取與車(chē)輛運(yùn)行狀態(tài)、車(chē)輛速度等相關(guān),因此需要通過(guò)后期根據(jù)現(xiàn)場(chǎng)實(shí)際情況進(jìn)行標(biāo)定調(diào)整。
(2)軌道波磨故障波長(zhǎng)的檢測(cè)
通過(guò)對(duì)信號(hào)的小波包能量熵計(jì)算可以判斷軌道波磨故障是否存在,本節(jié)通過(guò)對(duì)信號(hào)進(jìn)一步分析來(lái)確定故障段的時(shí)域位置及波磨波長(zhǎng)。
選取表1中的第3組、第6組和第9組數(shù)據(jù)進(jìn)行分析。對(duì)分解后的子節(jié)點(diǎn)信號(hào)進(jìn)行小波包時(shí)頻分析,3組數(shù)據(jù)的二維時(shí)頻圖如圖8~圖10所示。
圖8 第3組實(shí)驗(yàn)時(shí)頻圖分析
圖9 第6組實(shí)驗(yàn)時(shí)頻圖分析
圖10 第9組實(shí)驗(yàn)時(shí)頻圖分析
圖11 第13組實(shí)驗(yàn)時(shí)頻圖分析
根據(jù)小波包分解可知4層小波包分解的16個(gè)節(jié)點(diǎn)分別對(duì)應(yīng)一個(gè)頻率范圍,由于采樣頻率為1 kHz,因此每個(gè)頻帶的范圍為500/16=31.25 Hz。根據(jù)第一節(jié)中的推導(dǎo)可知,波磨長(zhǎng)度和軸箱振動(dòng)頻率有如下關(guān)系
式中,λ為波磨波長(zhǎng);ω為振動(dòng)頻率,v為行駛車(chē)速。
由圖8~圖10分別得出波磨故障引起的振動(dòng)中心頻率為:66.67、50.3、39.91 Hz,得出對(duì)應(yīng)波磨波長(zhǎng)為149.99、198.81、250.56 mm。觀察圖11,a處故障中心頻率為65.22 Hz,b處故障中心頻率為49.78 Hz,計(jì)算對(duì)應(yīng)波磨波長(zhǎng)分別為153.32、200.88 mm。可見(jiàn)診斷結(jié)果和實(shí)際波磨設(shè)置參數(shù)基本吻合,最大誤差為2.21%。
觀察圖8~圖10可看出,3幅圖中均存在不止一處的高頻率集中區(qū)域,如圖8~圖10中方框a、b、c所示。而在仿真中只設(shè)定了1處波磨故障,b、c兩處的高頻率集中區(qū)域呈現(xiàn)很明顯的橫向集中性,即頻率跨度很大而里程跨度很小,并且兩個(gè)高頻率集中區(qū)域分別是在設(shè)定波磨段的首尾段。因此可以認(rèn)為這兩處的頻率集中區(qū)域是由于車(chē)輛經(jīng)過(guò)波磨和非波磨段過(guò)度區(qū)域時(shí)引起的沖擊變化,而真正的波磨段引起的時(shí)頻圖變化區(qū)域應(yīng)該是a區(qū)域所示。
本文基于小波包分解和信息熵理論研究了一種基于小波包能量熵的鋼軌波磨故障檢測(cè)方法。首先建立了軌道波磨簡(jiǎn)化模型和列車(chē)集總參數(shù)簡(jiǎn)化模型,驗(yàn)證了基于軸箱振動(dòng)信號(hào)檢測(cè)軌道波磨的可行性。然后利用SIMPACK多體動(dòng)力學(xué)軟件對(duì)模型進(jìn)行仿真。并創(chuàng)新性地將小波包能量熵運(yùn)用于軌道波磨檢測(cè)上;隨后通過(guò)設(shè)定14組不同類(lèi)型的軌道波磨故障形式進(jìn)行仿真實(shí)驗(yàn)。仿真實(shí)驗(yàn)結(jié)果表明,小波包能量熵能夠很好地判斷鋼軌波磨故障的情況。
[1] 高雨峰,王福田,劉仍奎.鐵路軌道狀態(tài)分析信息系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J].中國(guó)安全科學(xué)學(xué)報(bào),2006,16(7):128-132.
[2] Grassie S L. Measurement of railhead longitudinal profiles: a comparison of different techniques[J]. Wear, 1996, 191(1):245-251.
[3] Tsunashima H, Mori H, Yanagisawa K, et al. Condition monitoring of railway tracks using compact size on-board monitoring device[C]∥Railway Condition Monitoring (RCM 2014), 6th IET Conference on. IET,2014:1-5.
[4] 曹西寧,柴曉冬,鄭樹(shù)彬.基于Hilbert-Huang變換的軌道車(chē)輛軸箱加速度信號(hào)分析[J].儀表技術(shù)與傳感器,2015(3):92-95.
[5] 李宏鋒.地鐵軌道波浪形磨耗檢測(cè)系統(tǒng)研究[D].成都:西南交通大學(xué),2014.
[6] 史紅梅.基于車(chē)輛動(dòng)態(tài)響應(yīng)的軌道不平順智能感知算法研究[D].北京:北京交通大學(xué),2013.
[7] 張健,金學(xué)松,肖新標(biāo),等.車(chē)輛-軌道耦合動(dòng)力學(xué)鋼軌模型求解方法[J].交通運(yùn)輸工程學(xué)報(bào),2011,11(2):32-38.
[8] 繆炳榮,方向華,傅秀通.動(dòng)力學(xué)分析基礎(chǔ)教程[M].成都:西南交通大學(xué)出版社,2008.
[9] Kurzeck B. Combined friction induced oscillations of wheelset and track during the curving of metros and their influence on corrugation[J]. Wear, 2011,271(1):299-310.
[10] 劉珊,曹海泉,于海,洪剛.基于小波包能量熵的變壓器振動(dòng)信號(hào)特征研究[J].電網(wǎng)與清潔能源,2010(5):35-38.
[11] 楊江天,趙明元,張志強(qiáng),李平康.基于定子電流小波包分析的牽引電機(jī)軸承故障診斷[J].鐵道學(xué)報(bào),2013(2):32-36.
[12] 肖玉飛,劉祖潤(rùn),李目.基于小波包能量熵與SVM的模擬電路故障診斷[J].電子測(cè)量技術(shù),2011(6):110-113.
[13] 武兵,曹勇敢.基于小波包能量熵的輸電線路兩相接地故障甄別[J].河南電力,2016(1):15-19.
[14] 趙志宏,楊紹普.基于小波包變換與樣本熵的滾動(dòng)軸承故障診斷[J].振動(dòng).測(cè)試與診斷,2012(4):640-644,692.
[15] 安春蘭,甘方成,羅微,秦發(fā)園.提速道岔小波包能量熵故障診斷方法[J].鐵道科學(xué)與工程學(xué)報(bào),2015(2):269-274.
[16] 李彥蒼,彭?yè)P(yáng).基于信息熵的改進(jìn)人工蜂群算法[J].控制與決策,2015(6):1121-1125.
[17] 趙永超.無(wú)砟軌道結(jié)構(gòu)振動(dòng)響應(yīng)的概率分析[D].長(zhǎng)沙:中南大學(xué),2009.
A Rail Corrugation Detection Method Based on Wavelet Packet-energy entropy
DONG Wei1, HUANG Wen1, XING Zong-yi2, LONG Jing3
(1. School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China; 2. School of Automation, Nanjing University of Science and Technology, Nanjing 210094, China; 3. Guangzhou City Underground Railway Corporation, Guangzhou 510380, China)
Aiming at the problem of low efficiency and low success rate in rail corrugation detection, a method based on wavelet packet-energy entropy is proposed. Firstly, physical fitting models of rail corrugation and the train are established and a coupling dynamics model is set up based on SIMAPCK. Then wavelet packet-energy entropy is introduced creatively in the area of rail corrugation detection. Furthermore, four-layer wavelet packet decomposition of the axle box vibration acceleration signal is performed, and wavelet pocket-energy entropy and wavelet energy of each nod are calculated. Finally the simulation results are compared and analyzed to determine the type of rail corrugation.
Rail; Corrugation; Fault diagnosis; Axis box vibration acceleration; Wavelet packet; Energy entropy
1004-2954(2018)01-0052-07
2017-02-27;
2017-04-07
國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFB1200402);國(guó)家工程實(shí)驗(yàn)室建設(shè)項(xiàng)目(2016582)
董 偉(1992—),男,碩士研究生,主要從事運(yùn)輸信息工程與安全保障研究,E-mail:njustdongwei@163.com。
U213.4+2
A
10.13238/j.issn.1004-2954.201702270001