国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

基于MKurt-MOMEDA的齒輪箱復(fù)合故障特征提取*

2017-09-12 07:07王志堅(jiān)王俊元趙志芳吳文軒張紀(jì)平寇彥飛
振動、測試與診斷 2017年4期
關(guān)鍵詞:峭度齒輪箱特征提取

王志堅(jiān),王俊元,趙志芳,吳文軒,張紀(jì)平,寇彥飛

(中北大學(xué)機(jī)械與動力工程學(xué)院 太原, 030051)

基于MKurt-MOMEDA的齒輪箱復(fù)合故障特征提取*

王志堅(jiān),王俊元,趙志芳,吳文軒,張紀(jì)平,寇彥飛

(中北大學(xué)機(jī)械與動力工程學(xué)院 太原, 030051)

針對齒輪箱中旋轉(zhuǎn)零部件的故障信號是周期性的沖擊信號這一特性,提出了一種基于多點(diǎn)峭度(multipoint kurtosis,簡稱MKurt)和多點(diǎn)最優(yōu)最小熵反褶積(multipoint optimal minimum entropy deconvolution adjusted,簡稱MOMEDA)的齒輪箱復(fù)合故障特征提取方法。利用MKurt可以有效提取齒輪箱中被噪聲淹沒的沖擊性振動信號的周期,實(shí)現(xiàn)對振動信號振動源的追蹤。根據(jù)故障的周期設(shè)置合理的周期區(qū)間,通過MOMEDA對原信號進(jìn)行降噪,進(jìn)一步提取原信號的周期性沖擊。通過仿真信號和實(shí)測數(shù)據(jù)的分析和驗(yàn)證,證明了MKurt-MOMEDA方法可以準(zhǔn)確有效地診斷齒輪箱復(fù)合故障故障特征。

多點(diǎn)峭度;最優(yōu)最小熵反褶積;復(fù)合故障;特征提取

引 言

齒輪箱是復(fù)雜機(jī)械系統(tǒng)使用最廣泛的機(jī)械,其中齒輪、軸承等零部件也最易損壞。當(dāng)齒輪箱中軸由于不對中、齒輪和軸承的內(nèi)外圈或滾動體發(fā)生故障時(shí),傳動系統(tǒng)將會受到影響,振動信號會出現(xiàn)周期性脈沖信號[1-3]。當(dāng)齒輪箱出現(xiàn)多種故障時(shí),傳感器測取的信號為各種源信號的疊加信號,這些源信號在傳遞過程中往往發(fā)生不同程度的畸變或混疊,在復(fù)合故障中表現(xiàn)的尤為明顯。因此,復(fù)合故障特征之間互相耦合、相互干擾,尤其在強(qiáng)背景噪聲下,微弱故障常被噪聲淹沒,難以提取特征信息??梢姡喙收咸卣鞯奶崛∪匀皇且粋€大的挑戰(zhàn)[4-6]。

峭度作為無綱量參數(shù),其值與軸承或齒輪的轉(zhuǎn)速、尺寸和載荷均無關(guān),但對沖擊信號特別敏感,適用于表面損傷、早期裂紋的故障診斷[7],可通過峭度反映信號的奇異程度,但其不能提供更多的關(guān)于故障位置等信息。文獻(xiàn)[8]提出了基于峭度圖的信號處理方法,并成功應(yīng)用于軸承和齒輪的故障診斷中。文獻(xiàn)[9]將最小熵反褶積(minimum entropy deconvolution,簡稱MED)用于滾動軸承與齒輪故障診斷,此方法的基本原理是利用最大峭度值作為迭代終止條件[10],通過解卷積最優(yōu)的濾波器,使信號突出少數(shù)大的尖脈沖。王志堅(jiān)等[11-12]通過MED和總體經(jīng)驗(yàn)?zāi)B(tài)分解(ensemble empirical mode decomposition,簡稱EEMD)二次降噪的方法提取軸承的微弱故障特征,但MED降噪只能突出少數(shù)大的尖脈沖,對復(fù)合故障的微弱成分的特征提取無能為力。為了彌補(bǔ)MED的局限性,McDonald等[13]提出最大相關(guān)峭度反褶積(maximum correlated kurtosis deconvolution,簡稱MCKD)算法。此方法引入相關(guān)峭度的概念,以相關(guān)峭度最大化作為目標(biāo)函數(shù),通過迭代算法優(yōu)化最佳濾波器,成功地運(yùn)用在齒輪的故障診斷中。利用MCKD方法可以有效提取故障信號中被噪聲淹沒的周期沖擊成分,抑制信號中的噪聲,實(shí)現(xiàn)信號降噪,提升原信號的峭度[14]。MCKD受兩個參數(shù)的限制,分別為濾波器的長度和故障周期,使其不具有自適應(yīng)性。通過大量仿真信號分析發(fā)現(xiàn),濾波器長度增大將增加計(jì)算量,故障周期必須取整,若為非整數(shù),需要對其重采樣并取整。由于MCKD只能提取出有限個沖擊成分,因此該方法在很大程度上受到約束。McDonald等[15]提出了MOMEDA的齒輪箱故障診斷方法,此方法適用于周期性故障特征提取,它使用一個時(shí)間目標(biāo)向量定義解反褶積確定所得的脈沖位置。該算法不需要迭代即可獲取最佳濾波器,也不需先驗(yàn)確定故障的周期。筆者通過多點(diǎn)峭度譜圖確定故障周期的大小,設(shè)置合理的故障區(qū)間,通過MOMEDA提取故障特征。該方法有效識別封閉式功率流試驗(yàn)臺的復(fù)合故障特征。

1 多點(diǎn)優(yōu)化最小熵反褶積

Cabrelli[16]對脈沖反褶積提出了一種新的準(zhǔn)則,稱為D范數(shù),并用幾何學(xué)證明了反褶積問題。D范數(shù)反褶積問題可以通過一個準(zhǔn)確的非迭代求解方法對濾波器系數(shù)進(jìn)行求解。文獻(xiàn)[15]提出一種基于多重D范數(shù)(multi D-norm,簡稱MDN)的多脈沖目標(biāo)識別的反褶積算法主要應(yīng)用于周期性脈沖的故障特征提取,并引入了MOMEDA。

假設(shè)采集到的響應(yīng)信號為

(1)

其中:e(n)為噪聲;x(n)為沖擊序列;h(n)為傳遞函數(shù);y(n)為采集到的振動信號。

MOMEDA算法本質(zhì)就是尋找一個FIR濾波器,通過輸出信號x(n)盡可能恢復(fù)輸入信號y(n)。

MOMEDA是一種針對位置已知的多脈沖目標(biāo)識別反褶積算法,能夠識別連續(xù)的沖擊脈沖。引入的最大值問題的解即為MOMEDA。

(2)

(3)

其中:t為確定脈沖位置和權(quán)重的常數(shù)矢量。

利用歸一化水平表示達(dá)到的最佳目標(biāo)解,可提取不同采樣率的故障周期,也可以識別同一采樣頻率下不同故障特征的周期。因此,通過目標(biāo)矢量t可實(shí)現(xiàn)脈沖信號的分離和位置的確定。

通過對濾波器系數(shù)(f=f1,f2,…,fL)求導(dǎo)數(shù),得到式(3)的極值為

(4)

(5)

進(jìn)一步簡化為

t1M1+t2M2+…+tN-LMN-L=X0t

通過導(dǎo)數(shù)等于0解出極值,式(5)變?yōu)?/p>

(6)

MOMEDA濾波器和輸出解可概括為

(7)

(8)

(9)

目標(biāo)矢量t和輸出結(jié)果具有相同的長度(N-L+1),目標(biāo)矢量表示在輸出中解反卷積脈沖的位置和權(quán)重,控制了脈沖的位置,且MOMEDA方法通過非迭代取得最優(yōu)解,為旋轉(zhuǎn)機(jī)械周期性故障診斷提供了基礎(chǔ)。OMEDA可對M個連續(xù)的目標(biāo)矢量進(jìn)行如下計(jì)算,此時(shí)式(7),(9)變成式(10),(11)。

(10)

(11)

通過MOMEDA降噪,求解反褶積的位置是唯一的,對單一故障而言,不需要考慮周期是否為整數(shù)或?yàn)V波器的長度對降噪的影響。在整個采樣區(qū)間內(nèi),每個周期都會出現(xiàn)一個周期性沖擊脈沖。多級傳動的齒輪箱內(nèi)部特征頻率分布較寬,故障周期為多個時(shí),對復(fù)合故障的追蹤效果不佳。為了準(zhǔn)確提取復(fù)合故障特征,引入了多點(diǎn)峭度作為特征提取的度量標(biāo)準(zhǔn)。

(12)

當(dāng)輸出結(jié)果y與目標(biāo)矢量t相同時(shí),MKurt歸一化得到

(13)

得到標(biāo)準(zhǔn)化因子

(14)

最終,標(biāo)準(zhǔn)化的多點(diǎn)峭度被定義為

(15)

此定義是基于峭度提出的,但其目標(biāo)矢量在受控位置處被擴(kuò)展為多個脈沖,進(jìn)一步被歸一化。實(shí)際上旋轉(zhuǎn)機(jī)械每轉(zhuǎn)動一周,可能存在兩個或多個故障沖擊成分,多點(diǎn)峭度達(dá)到峰值時(shí)可能對應(yīng)的故障周期并非只有一個。當(dāng)多點(diǎn)峭度達(dá)到峰值時(shí),對應(yīng)的采樣點(diǎn)數(shù)(周期)即為故障周期,也可能為周期的整數(shù)倍或半倍等。因此,多點(diǎn)峭度可以對故障周期和周圍的非故障進(jìn)行區(qū)分,其中多點(diǎn)峭度和采樣點(diǎn)數(shù)均為無量綱量。

為了提取復(fù)合故障特征,首先,計(jì)算原信號的多點(diǎn)峭度,確定復(fù)合故障的周期;其次,設(shè)置合理的故障周期區(qū)間,假設(shè)某復(fù)合故障存在兩個故障周期,分別為40和120,如果設(shè)置一個故障周期區(qū)間為[10~200],通過MOMEDA降噪時(shí)只能提取較強(qiáng)的沖擊成分,另一故障將被誤認(rèn)為噪聲造成誤診斷現(xiàn)象,此時(shí)可設(shè)置兩個周期區(qū)間,分別為[10~80]和[100~150],每個周期區(qū)間含有唯一的故障周期;最后,分別對原信號進(jìn)行MOMEDA降噪,即可分別獲取這兩個周期區(qū)間的沖擊脈沖。

仿真信號如圖1所示。復(fù)合信號e由噪聲信號a、正弦信號b和沖擊信號1組成。其中,沖擊性周期為90,信號e可知沖擊信號1已被噪聲淹沒。為了確定沖擊性信號1的周期,對其求解多點(diǎn)峭度,結(jié)果如圖2所示。在周期為30,45,90,180和270處均出現(xiàn)峰值,分別對應(yīng)于沖擊性信號1周期的半倍、1/3和整數(shù)倍關(guān)系。如圖3所示,為了進(jìn)一步提取沖擊脈沖信號1,對組合信號e進(jìn)行MOMEDA降噪,成功提取周期為90的沖擊信號。

文獻(xiàn)[15]和上述仿真信號分別對單一故障特征提取,為了驗(yàn)證此方法提取多故障特征的能力,對含有兩個沖擊成分如圖1的仿真信號f進(jìn)行分析,它包含了兩個沖擊成分,周期分別為90和140。圖4對應(yīng)其多點(diǎn)峭度譜圖,顯然在譜峰處周期為30,45,90,180,270,360分別對應(yīng)沖擊信號1的1/3倍、0.5倍和整數(shù)倍;70,140,280分別對應(yīng)沖擊信號2的0.5倍和整數(shù)倍。多點(diǎn)峭度可以確定多故障對應(yīng)的周期,為了進(jìn)一步提取故障信號,需要設(shè)定恰當(dāng)?shù)闹芷诜秶?。因?yàn)閮蓚€沖擊周期分別為90和140,因此分別設(shè)定周期為[50~100],[100,150],設(shè)定步長ΔT=0.1,經(jīng)過MOMEDA降噪成功提取出沖擊信號1和2,如圖5,6所示。

圖1 仿真信號Fig.1 Simulation signal

圖2 仿真信號e的多點(diǎn)峭度譜圖Fig.2 Multipoint kurtosis spectrum of simulation signal e

圖3 仿真信號e的MOMEDA降噪Fig.3 MOMEDA noise reduction of Simulation signal e

圖4 仿真信號f的多點(diǎn)峭度譜圖Fig.4 Multipoint kurtosis spectrum of simulation signal f

圖5 仿真信號f的MOMEDA降噪提取沖擊信號1Fig.5 MOMEDA noise reduction of Simulation signal f to extract the impact signal 1

圖6 仿真信號f的MOMEDA降噪提取沖擊信號2Fig.6 MOMEDA noise reduction of Simulation signal f to extract the impact signal 2

2 齒輪箱多故障特征識別

試驗(yàn)裝置如圖7所示,試驗(yàn)齒輪為18個,試驗(yàn)軸承型號為32212,試驗(yàn)齒輪的傳動比為1∶1,采取半齒嚙合,轉(zhuǎn)速為1 200r/min,采樣頻率為8kHz,正常和故障的齒輪和軸承試驗(yàn)負(fù)載均加載到1kNm,采集信號的傳感器型號為YD77SA三向加速度傳感器(靈敏度為0.01V/ms2),采樣數(shù)為2 048個,軸的轉(zhuǎn)動周期為400,軸承滾珠故障周期為111.1,齒輪的嚙合頻率為360Hz,齒輪的嚙合周期為22.4。齒輪箱故障類型為兩個:第1個為齒輪點(diǎn)蝕故障,如圖8所示;第2個為軸對不中,齒輪軸每轉(zhuǎn)動一周可產(chǎn)生一個周期性沖擊[12]。

圖7 齒輪傳動試驗(yàn)臺Fig.7 Rig for gear transmission testing

圖8 軸承齒輪故障示意圖Fig.8 Gear fault diagram

圖9 健康齒輪時(shí)域波形Fig.9 Time-domain waveform of healthy gear

圖10 點(diǎn)蝕齒輪時(shí)域波形Fig.10 Time-domain waveform of pitting gear

圖11 健康齒輪箱多點(diǎn)峭度譜圖Fig.11 Multipoint kurtosis spectrum of healthy gear box

圖12 故障齒輪箱多點(diǎn)峭度譜圖Fig.12 Multipoint kurtosis spectrum of faulty gear box

圖13 故障齒輪箱MOMEDA降噪Fig.13 MOMEDA noise reduction of faulty gear box

同一個加速度傳感器同一個方向測得齒輪箱健康和故障的振動信號時(shí)域波形如圖9,10所示。正常齒輪振動波形比較平滑且幅值較小,當(dāng)齒輪箱出現(xiàn)故障時(shí),會出現(xiàn)明顯沖擊振動,幅值有所增加,基本上是正常齒輪幅值的4倍左右,且有明顯周期成分出現(xiàn),周期為400,對應(yīng)軸的轉(zhuǎn)動周期,另一沖擊成分在時(shí)域波形中并不能確定。為了確定故障位置,分別對健康和故障下的振動數(shù)據(jù)進(jìn)行多點(diǎn)峭度譜圖分析,如圖11,12所示。顯然,周期為22.4,44.8,67.2和134.4代表齒輪的嚙合周期及其倍數(shù)的關(guān)系,100.8,201.6和403.2代表軸的轉(zhuǎn)動周期及其因子。通過對比可知,齒輪故障多點(diǎn)峭度的峰值明顯強(qiáng)于健康齒輪箱的峰值,峰值保持在兩倍左右,說明齒輪有缺陷。如圖13所示,為了進(jìn)一步提取沖擊成分,設(shè)置周期區(qū)間為[10~50]和[350~450],通過MOMEDA降噪,成功提取齒輪的點(diǎn)蝕故障和軸的不對中故障。

3 結(jié) 論

1) MOMEDA克服了MED和MCKD的缺點(diǎn),可以將每個周期的故障沖擊提取出來,但只能對單一故障特征提取。

2) MKurt能夠求解復(fù)合故障沖擊性周期,在周期的倍數(shù)或半倍處均有譜峰出現(xiàn)。在噪聲環(huán)境中此方法具有免疫性,為了進(jìn)一步提取沖擊性脈沖,可通過設(shè)置降噪周期的區(qū)間來控制解卷積的脈沖位置。

3) 該方法為復(fù)合故障特征提取提供了一種新思路,具有一定的工程應(yīng)用價(jià)值。

[1] Georgoulas G,Loutas T,Stylios C D,et al.Bearing fault detection based on hybrid ensemble detector and empirical mode decomposition[J].Mechanical Systems &Signal Processing,2013,41(1-2):510-525.

[2] Randall R B,Antoni J.Rolling element bearing diagnostics-a tutorial[J].Mechanical Systems &Signal Processing,2011,25(2):485-520.

[3] Lei Yaguo,He Zhengjia,Zi Yangyang.EEMD method and WNN for fault diagnosis of locomotive roller bearings[J].Expert Systems with Applications an International Journal,2011,38(6):7334-7341.

[4] Bachschmid N,Pennacchi P,Vania A.Identification of multiple faults in rotor systems[J].Journal of Sound &Vibration,2002,254(2):327-366.

[5] Patel T H,Darpe A K.Coupled bending-torsional vibration analysis of rotor with rub and crack[J].Journal of Sound &Vibration,2009,326(3):740-752.

[6] Wang Yanxue,Liang Ming.Identification of multiple transient faults based on the adaptive spectral kurtosis method[J].Journal of Sound &Vibration,2012,331(2):470-486.

[7] 張超,陳建軍.基于LMD和Lempel-Ziv指標(biāo)的滾動軸承故障損傷程度研究[J].振動與沖擊,2012(16):77-82.

Zhang Chao,Chen Jianjun.Fault severity assessment for rolling element bearings based on LMD and Lempel-Ziv index[J].Journal of Vibration and Shock,2012(16):77-82.(in Chinese)

[8] 彭暢,柏林,劉小峰.基于魯棒性小波包峭度圖的滾動軸承故障診斷[J].振動、測試與診斷,2016(1):11-16.

Peng Chang,Bo Lin,Liu Xiaofeng.Robust wallet transform-based kurtogram for the fault diagnostics of rolling element bearing[J].Journal of Vibration,Measurement &Diagnosis,2016(1):11-16.(in Chinese)

[9] Sawalhi N,Randall R B,Endo H.The enhancement of fault detection and diagnosis in rolling element bearings using minimum entropy deconvolution combined with spectral kurtosis[J].Mechanical Systems &Signal Processing,2007,21(6):2616-2633.

[10]張丹,隋文濤.MED方法及在滾動軸承故障診斷中的應(yīng)用[J].機(jī)床與液壓,2014(23):185-187.

Zhang Dan,Sui Wentao.MED method and application in fault diagnosis of rolling element bearings[J].Machine Tool &Hydraulics,2014(23):185-187.(in Chinese)

[11]王志堅(jiān),韓振南,劉邱祖,等.基于MED-EEMD的滾動軸承微弱故障特征提取[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(23):70-78.

Wang Zhijian,Han Zhennan,Liu Qiuzu,et al.Weak fault feature extraction of rolling bearing based on MED-EEMD[J].Transactions of the Chinese Society of Agricultural Engineering,2014,30 (23):70-78.(in Chinese)

[12]Wang Zhijian,Han Zhennan.A novel procedure for diagnosing multiple faults in rotating machinery[J].ISA Transactions,2015 (55):208-218.

[13]McDonald G L,Zhao Qing,Zuo Ming.Maximum correlated kurtosis deconvolution and application on gear tooth chip fault detection[J].Mechanical Systems and Signal Processing,2012,33:237-255.

[14]武兵,賈峰,熊曉燕.基于最大相關(guān)峭度反褶積的軸承故障診斷方法[J].振動、測試與診斷,2014(3):570-575.

Wu Bing,Jia Feng,Xiong Xiaoyan.Bearing fault diagnosis based on maximum correlated kurtosis deconvolution[J].Journal of Vibration,Measurement &Diagnosis,2014(3):570-575.(in Chinese)

[15]Mcdonald G L,Zhao Qing.Multipoint optimal minimum entropy deconvolution and convolution fix:application to vibration fault detection[J].Mechanical Systems &Signal Processing,2016,82:461-477.

[16]Cabrelli C A.Minimum entropy deconvolution and simplicity:a noniterative algorithm[J].Geophysics,1985,50(3):394-413.

10.16450/j.cnki.issn.1004-6801.2017.04.030

* 山西省自然科學(xué)基金資助項(xiàng)目(2015011063)

2017-02-27;

2017-04-26

TH113.1

王志堅(jiān),男,1985年1月生,博士、講師。主要研究方向?yàn)樾D(zhuǎn)機(jī)械復(fù)合故障診斷、動態(tài)測試。曾發(fā)表《基于MED-EEMD滾動軸承微弱故障特征提取》(《農(nóng)業(yè)工程學(xué)報(bào)》2014年第30卷第23期)等論文。 E-mail:wangzhijian1013@163.com

猜你喜歡
峭度齒輪箱特征提取
基于重加權(quán)譜峭度方法的航空發(fā)動機(jī)故障診斷
風(fēng)電齒輪箱軸承用鋼100CrMnSi6-4的開發(fā)
CJ-1型齒輪箱箱體強(qiáng)度分析
聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對軸承復(fù)合故障檢測研究
空間目標(biāo)的ISAR成像及輪廓特征提取
基于Gazebo仿真環(huán)境的ORB特征提取與比對的研究
基于特征提取的繪本閱讀機(jī)器人設(shè)計(jì)方案
基于Daubechies(dbN)的飛行器音頻特征提取
基于加權(quán)峭度的滾動軸承故障特征提取
譜峭度在軸承故障振動信號共振頻帶優(yōu)選中的應(yīng)用