基于DSP+MCU的列車(chē)滾動(dòng)軸承故障診斷系統設計與應用
1.2.2分析處理程序
本系統利用DSP強大的數字信號處理功能,對采樣得到的數據進(jìn)行FFT運算和功率譜分析,更好地提取數據中的特征信息,加快系統的響應速度和提高準確度。
(1)FFT分析運算子程序
FFT分析運算子程序利用FFT計算相關(guān)函數。為防止發(fā)生頻疊現象,需要延長(cháng)線(xiàn)性相關(guān)中序列的長(cháng)度,即延長(cháng)到兩序列長(cháng)度之和2N。自相關(guān)函數的快速傅里葉變換計算過(guò)程如下:

(2)功率譜分析運算子程序
平均周期功率譜分析首先要把序列X(n)分成K段,每段長(cháng)為N,然后對每段進(jìn)行功率譜分析。平均周期法的每一段譜分析就是求該段的離散傅里葉變換,再除以分析點(diǎn)數。這樣的譜估計一共有K段,對K段譜估計求平均就得到平均周期功率譜分析。
(3)倒譜分析運算子程序
倒譜分析是對信號y(t)的功率譜的對數進(jìn)行傅里葉逆變換。倒譜分析的離散運算形式為:

2 故障診斷
2.1 提取軸承特征
小波分析利用時(shí)間平移和多分辨率的概念,可以同時(shí)處理時(shí)、頻分析,具有時(shí)頻局部化和多分辨功能。其基本思想是用一族函數去表示或逼近一信號或函數,通過(guò)滿(mǎn)足一定條件的基本小波函數的不同尺度的平移和展縮構成的。但在正交小波變換中,只對信號的低頻成分進(jìn)行了遞推分解,導致高頻成分的頻率分辨率較低,表現為時(shí)一頻分辨率在低頻處頻率分辨率高,在高頻處時(shí)間分辨率高,頻率分辨率卻降低。利用DSP強大的數字處理功能,本系統采用常見(jiàn)的Hilbert變換法來(lái)提取包絡(luò )信號,提高整個(gè)系統的可靠性和精確性。
實(shí)信號X(t)的Hilbert變換為:

2.2 小波奇異性檢測
函數f(x)的局部奇異性與其小波變換的漸進(jìn)衰減性之間的關(guān)系為:

式中:Wsf(x)為f(x)在尺度s上的小波變換。
本文根據小波變換各尺度上模極大值的傳遞性來(lái)判斷奇異點(diǎn)的位置以及作奇異性指數計算。奇異性指數的計算如下:
設s=2j,在尺度i上Xk處的極大值為Mi=| Wsf(x)|,則在各尺度相應位置處的模極大值可構成序列{Mi},在i較小時(shí),可以近似為:

由此可得:

根據上式計算幾個(gè)尺度上的α,然后求平均值,即可得到信號在此時(shí)刻的Lip指數。
3 實(shí)例分析
實(shí)驗用軸承參數如下:滾動(dòng)體直徑:O.84235英寸;支架直徑:7.5653英寸;輪子直徑:35.89英寸;接觸角α:10°;車(chē)速:30 km/h。
當軸承外圈滾道發(fā)生點(diǎn)蝕、裂紋及表面剝落等局部損傷故障后,滾動(dòng)軸承便產(chǎn)生沖擊振動(dòng)。利用加速度傳感器獲取軸承振動(dòng)信號,采樣頻率為261436SPS,滾動(dòng)軸承正常、滾子破裂、多處外圈剝落時(shí)振動(dòng)信號的時(shí)域波形如圖3所示。按照前述方法對外圈剝落振動(dòng)信號進(jìn)行包絡(luò )處理,并采用B樣條函數進(jìn)行7層小波變換,得到信號包絡(luò )在特征尺度

評論