AR雙譜的電梯機(jī)械故障診斷

陳建燦，劉曉梅

（華僑大學(xué) 機(jī)電及自動(dòng)化學(xué)院，福建 廈門361021）

目前，廣泛使用的升降電梯主要有曳引式和液壓式兩種，高層樓宇多使用曳引式電梯.曳引式電梯通過安全裝置，防止轎箱因鋼纜斷裂、制動(dòng)失靈等原因造成的墜落.電梯主要是由機(jī)械系統(tǒng)、拖動(dòng)系統(tǒng)和電氣控制部分組成.拖動(dòng)系統(tǒng)也可以歸入電氣系統(tǒng)，所以電梯的故障就可以分為機(jī)械故障和電氣故障兩大類.高階譜分析能有效地消除信號(hào)中的高斯噪聲，在反映信號(hào)能量的同時(shí)，還保留信號(hào)的相位信息，因此能定量描述非線性相位耦合.本文通過分析正常和故障情況下的振動(dòng)信號(hào)，研究其非線性特征及區(qū)別，從而實(shí)現(xiàn)故障模式識(shí)別和故障診斷.

1 AR模型和AR雙譜

1.1 AR模型［1-4］

假設(shè)a（t）為電梯輸出振動(dòng)信號(hào)中的隨機(jī)信號(hào)受到均值為零的非高斯白噪聲的干擾，x（t）為零均值有色非高斯噪聲，則其AR模型為

式（1）中：?i（i＝1，2，…，p）為自回歸系數(shù)；p為自回歸模型的階數(shù).

用滯后量m和n表示變量x（t）的三階累積量函數(shù)c（m，n），有

令m＝n，可得｛x（t）｝三階自相關(guān)c（m，n）的對(duì)角切片值.當(dāng)m＝n＝1，2，…，p，可得矩陣方程［5-6］為

1.2 AR雙譜

隨機(jī)量x（t）的雙譜為三階矩的二維Fourier變換［1，7-8］，即

式（4）中：ω＝2πf；｜ω1｜≤π；｜ω2｜≤π；｜ω1+ω2｜≤π；F（ω）為x（t）的Fourier變換；F＊（ω）為F（ω）的共軛復(fù)數(shù).

雙譜B（ω1，ω2）為復(fù)函數(shù)［9］，可表示為

式（5）中：H（ω）為X（n）在三階統(tǒng)計(jì)意義下的 AR模型；H＊（ω1+ω2）為H（ω1+ω2）的共軛復(fù)數(shù).

當(dāng)ω1＝ω2時(shí)，可得雙譜的一維對(duì)角切片譜，其定義［10］為

2 實(shí)驗(yàn)原理

實(shí)驗(yàn)采用福建泉州中僑富士電梯有限公司提供的實(shí)驗(yàn)電梯.在轎廂運(yùn)動(dòng)的垂直方向、前后方向分別安裝的傳感器，采集電梯在正常工作和故障狀態(tài)下的加速度振動(dòng)信號(hào)，并對(duì)其進(jìn)行分析、處理.實(shí)驗(yàn)電梯的工作狀態(tài)有3種：正常運(yùn)行（無故障）；故障1（繩頭張緊）；故障2（導(dǎo)靴左右夾緊）.在實(shí)驗(yàn)中設(shè)置采樣頻率為1 024Hz，讀取頻率為512Hz.利用NI軟件［11］構(gòu)建檢測(cè)平臺(tái)，使用PCI-604數(shù)據(jù)采集卡、加速度傳感器、信號(hào)放大器等設(shè)備進(jìn)行數(shù)據(jù)采集.

3 實(shí)驗(yàn)數(shù)據(jù)分析

對(duì)采集到的振動(dòng)加速度信號(hào)進(jìn)行預(yù)處理后分析，數(shù)據(jù)長(zhǎng)度為1 024點(diǎn)，繪制出電梯正常工作和故障情況下的振動(dòng)加速度信號(hào)的雙譜圖、等高線圖及其對(duì)角切片圖.垂直方向加速度振動(dòng)信號(hào)雙譜圖、雙譜等高線圖、及雙譜對(duì)角切片圖，如圖1～3所示.前后方向加速度振動(dòng)信號(hào)雙譜圖、雙譜等高線圖、及雙譜對(duì)角切片，如圖4～6所示.

圖1 垂直方向加速度振動(dòng)信號(hào)的雙譜圖Fig.1 Double spectrum of the acceleration vibration signal in the vertical direction

圖2 垂直方向的雙譜等高線圖Fig.2 Dual spectrum contour in the vertical direction

圖3 垂直方向的雙譜對(duì)角切片圖Fig.3 Bispectrum diagonal slice in the vertical direction

圖4 前后方向加速度振動(dòng)信號(hào)的雙譜圖Fig.4 Double spectrum of the acceleration vibration signal in the fore-and-aft direction

圖5 前后方向的雙譜等高線圖Fig.5 Dual spectrum contour in the fore-and-aft direction

圖6 前后方向的雙譜對(duì)角切片圖Fig.6 Bispectrum diagonal slice in the fore-and-aft direction

從圖1，4中可以看出：電梯振動(dòng)信號(hào)的雙譜圖都存在明顯的譜峰；在正常情況下，垂直方向和前后方向加速度信號(hào)的雙譜譜峰的數(shù)目都比較少，位置集中在兩對(duì)角處，譜峰明顯，小譜峰少；而故障1的情況，垂直方向和前后方向加速度信號(hào)的雙譜圖譜的譜峰分散在四角處，譜峰的數(shù)目較多，能量比較分散，小譜峰多；故障2的情況，垂直方向和前后方向加速度信號(hào)的譜峰數(shù)目較正常情況的多且分散，垂直方向信號(hào)的雙譜譜峰在兩對(duì)角處分散開來，有3個(gè)主峰，但譜峰和能量較故障1的集中，而前后方向信號(hào)的雙譜譜峰分散在四角處，較故障1的峰陡且多.

從圖2，5中可以看出：雙譜的能量都集中（-π，π）之間，在正常情況下，能量較為集中，說明系統(tǒng)運(yùn)行較平穩(wěn)；在故障的情況下，由于系統(tǒng)運(yùn)行的不平穩(wěn)，能量較為分散；雙譜的對(duì)角切片更直觀的表達(dá)了峰值所對(duì)應(yīng)的頻率位置.

從圖3，6中可以看出：譜峰所對(duì)應(yīng)的頻率位置都不同，故障1峰值所占的范圍最窄，正常情況次之，而故障2峰值所占的范圍最寬；正常情況下峰值所占的范圍最窄，故障2次之，而故障1峰值所占的范圍最寬.

4 結(jié)論

1）高階譜能夠抑制高斯信號(hào)，提取系統(tǒng)的特征信號(hào)，是處理非高斯，非線性的有效方法.AR雙譜譜峰的形狀和分布反映了能量分布的差異，可以整體的觀察電梯系統(tǒng)的動(dòng)態(tài)特性.電梯在正常和故障工作時(shí)的AR雙譜譜峰形狀和分布差異明顯.

2）AR雙譜的等高線和對(duì)角切片圖更能直觀地表達(dá)雙譜的能量分布情況、譜峰的形狀和譜峰所對(duì)應(yīng)的頻率位置，更能清晰分辨出正常工作和各故障之間的區(qū)別，可以很好的判斷故障的存在.

3）通過分析雙譜、等高線及雙譜對(duì)角切片差異，可以對(duì)電梯的故障進(jìn)行有效地診斷.

［1］ 張賢達(dá).時(shí)間序列分析：高階統(tǒng)計(jì)量方法［M］.北京：清華大學(xué)出版社，1996：6-13，20-24，52-55.

［2］ 高宇，黃宜堅(jiān).AR 雙譜在減壓閥故障診斷中的應(yīng)用［J］.機(jī)械設(shè)計(jì)與制造，2011（11）：70-72

［3］ ZHANG Yan，WANG Shu-xun.Slice spectra approach to nonlinear phase coupling estimation［J］.Electronics Letters of IEEE，1998，34（4）：340-341.

［4］ 楊叔子，吳雅，軒建平，等.時(shí)間序列分析的工程應(yīng)用［M］.武漢：華中科技大學(xué)出版社，2007：24-52，184-218，222-226.

［5］ 柴茂，黃宜堅(jiān).基于 AR雙譜的調(diào)速閥故障診斷［J］.福州大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，36（6）：825-830.

［6］ 吳福森，黃宜堅(jiān).砌塊成形機(jī)液壓振動(dòng)信號(hào)的AR雙譜分析［J］.中國(guó)機(jī)械工程，2011，22（16）：1954-1959.

［7］ RAGH U M，NIKIAS C.Bispectrum estimation：A parametric approach［J］.IEEE Transactions on Acoustics，Speech and Signal Processing，1985，33（5）：1213-1230.

［8］ NIKIAS C L，RAGHUVEER M R.Bispectrum estimation：A digital signal processing framework［J］.Proceedings of the IEEE，1987，75（7）：869-891.

［9］ 呂金虎，陸君安，陳士華.混沌時(shí)間序列分析及其應(yīng)用［M］.武漢大學(xué)出版社，2001：114-183.

［10］ 湯琴，黃宜堅(jiān).采用AR模型雙譜估計(jì)的概率篩篩分效率［J］.華僑大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，32（3）：253-257.

［11］ 侯國(guó)屏，王珅，葉齊鑫，等.LabVIEW 7.1編程與虛擬儀器設(shè)計(jì)［M］.北京：清華大學(xué)出版社，2005：178-189.