復累量切片在故障診斷中的應用

吳文兵，黃宜堅

1.福州外語外貿學院，福州 350018

2.廈門大學 信息科技學院，福建 廈門 361000

3.華僑大學 機電工程學院，福建 泉州 362021

復累量切片在故障診斷中的應用

吳文兵1，2，黃宜堅3

1.福州外語外貿學院，福州 350018

2.廈門大學 信息科技學院，福建 廈門 361000

3.華僑大學 機電工程學院，福建 泉州 362021

1 引言

利用高階累積量可以自動地抑制高斯背景噪聲（有色或白色）的影響，正因為這樣，高階累積量這一統計量已日益受到人們的重視并已成為信號處理中一種非常有用的工具。一般而言，機械振動過程中工件的振動信號是非線性、非高斯的，環境中的噪聲可以近似當做高斯噪聲處理，因而用高階累量分析振動信號更容易提取與工件有關的特征信息。因此，高階累量已經在機械振動、故障診斷等領域得到了廣泛應用[1-4]。M.R.Raughveer和國內學者王樹勛等人對高階累量的耦合性質進行了深入研究，明確指出了各種復數高階累量的不同定義形式的耦合特征[5-8]。本文在此基礎上，利用復數三階累量切片進行故障診斷，實驗結果表明，由于定義形式的不同而導致的復數三階累量切片中包含了不同的耦合信息，這些不同的耦合信息使得故障診斷的準確率不一樣，說明了不同耦合方式的復數三階累量在保持同類信號特征這一點上具有不同的特點。

2 高階累積量

設{x(n)}為零均值k階平穩隨機過程，則該過程的k階累積量ckx(τ1，τ2,…，τk-1)定義為隨機變量{x(n)，x(n+τ1)，…，x(n+τk-1)}的k階聯合累積量，即

而該過程的k階矩mkx(τ1，τ2，…，τk-1)則定義為隨機變量{x(n)，x(n+τ1)，…，x(n+τk-1)}的k階聯合矩，即

這里，mom()表示聯合矩，三階累積量為[9]：

3 復數三階累量切片

根據文獻[7]，在式（3）中，令x(n)為復數信號，并進行如下定義：

其中x*(n)是x(n)的共軛復數，本文的復數信號由采集的原始信號進行Hilbert變換得到。文獻[6]指出，按定義1的方式，可以用復數諧波形式表示的發生了二次相位耦合的信號（如式（7），其中?3=?2+?1，諧波分量ω3是由諧波分量ω1和ω2通過二次相位耦合而成的，且ω3=ω1+ω2）其三階累量為0，以定義2的方式其三階累量如式（8）所示，以定義3的方式其三階累量如式（9）所示。

從式（3）可以看出，當τ2=0時所求出的即為三階累量切片，式（10）和式（11）分別表明，當信號可以由復數諧波形式表示時，按定義2和定義3方式得出的三階累量切片，分別只包含參加耦合的信號和耦合后的信號。

4 數據采集

本實驗研究的減壓閥為先導式減壓閥。當減壓閥進出油口有異物，壓力偏高或偏低都會影響到減壓閥的正常運行，為了獲取減壓閥在故障狀態下的運行信號，本文進行的實驗人為設置了如下故障：

故障一：在減壓閥出油口加?3 mm的鐵芯；

故障二：在減壓閥進，出油口都加了?3 mm鐵芯。

通過實驗可以近似模擬減壓閥工作故障狀況。

本文利用LabVIEW軟件及PCI-6014的數據采集卡和一個加速度傳感器，依次采集減壓閥在正常和故障狀態下的振動信號。在每種測量中，將油路壓力從1～5 MPa分五個壓力等級。采樣頻率250 Hz，讀取頻率125 Hz，采樣過程時間約2 min。本次實驗分別在減壓閥正常工作狀態和兩種故障狀態各自采集了18組共54組數據。本文實驗使用的數據個數為1 536個。

由于測試過程中系統外部和內部各種因素的影響必然在輸出過程中夾雜著不需要的成分，本文采用中值法對采集的振動信號進行預處理剔除混雜在信號中的干擾噪聲。濾波后的信號如圖1所示。圖中橫坐標表示信號個數，單位為個，縱坐標表示信號振幅，無量綱。

圖1 濾波后信號

5 三階累量切片計算

本文為了進行故障識別，將實驗中測得的正常狀態和故障一狀態的36組數據分別按式（4）～（6），并令τ2=0，計算出復數三階累量切片并取其模值，經過對比，在計算時令 -8≤τ1≤8，并以此計算結果作為特征向量輸入支持向量機進行故障判別，可以取得較為滿意的結果。其中按定義1計算出的累量切片示于表1。

表1 正常狀態和故障一狀態三階累量切片值

6 實驗結果

本文根據文獻[10]，建立了相同的LSSVM。在訓練LSSVM時，分別將正常狀態和故障狀態編碼為1和 -1，由于正常數據和故障數據分別是在五種不同的油壓下獲得的，訓練時分別在表1的正常數據和故障一數據中取出不同油壓狀態下相對應的5組數據共10組，對所建立的最小二乘支持向量機進行訓練，訓練結束后再利用該支持向量機進行故障識別，利用訓練后的網絡對所剩下的26組數據進行模擬識別，如圖2所示，圖中顯示26組數據的模擬后的狀態，其中空心圓圈表示訓練數據預先設定的狀態，比如前13組數據是測量的正常數據，應該處于圖中上一條橫線所表示的1狀態，后13組故障數據應該處于下一條橫線所表示的 -1狀態，紅色實心圓表示經過LSSVM模擬后對數據所進行的分類結果，如果某組數據的空心圓圈和紅色實心圓重合，則表示這組數據的模擬結果正確，否則錯誤。模擬后的結果如圖2所示，圖中的橫坐標表示信號的組數，單位為組，縱坐標表示信號模擬后的數值，無量綱。結果顯示，按定義1的方式進行識別，只有2組識別錯誤，正確率達到了90%以上。為了便于對比，將按定義2和定義3的方式進行相同的實驗，其識別錯誤數分別為2組和4組。同時將故障一的數據更換為故障二，也在三種不同的定義方式下進行了上述完全相同的實驗，其識別錯誤數分別為3組、2組和5組。

圖2 按定義1方式正常狀態和故障一狀態模擬結果圖

7 實驗結果分析

如前所述，按定義1的方式，可以用復數諧波形式表示的發生了二次相位耦合的信號其三階累量為0，按照定義2，三階累量中只包含參加二次耦合的諧波分量，而按定義3，則三階累量中既包含參加耦合的分量，也包含耦合后的分量。由上述不同的定義方式而計算出的三階累量切片，則同樣具有不同的耦合性質，這種不同的耦合性質就必然包含了不一樣的信號信息，在進行故障診斷時，由其得出的結果就可能會呈現出不同的正確率。

8 譜圖

為了對上述三種定義方式的三階累量切片在故障診斷中的性質有一個直觀上的認識，本文定義如下形式的三階累量切片譜：

圖3 定義1方式正常狀態三階累量切片譜圖

圖4 定義1方式故障狀態三階累量切片譜圖

圖5 定義2方式正常狀態三階累量切片譜圖

由于僅僅進行傅里葉變換，不會改變原信號中的頻率信息，因此三階累量切片譜有利于在直觀上對各種定義方式下的三階累量切片進行頻率信息的分析。本實驗首先將正常狀態和故障一狀態的36組數據，分別按前述的三種定義方式得出每組數據的三階累量切片共108組數據，再按式（12）求出其復數三階累量切片譜。本文在正常狀態和故障一狀態下，當油壓分別為1 MPa，3 MPa，5 MPa時，將兩種狀態的三階累量切片譜各選取1組示于圖3至圖8。

圖中橫軸表示圓頻率，單位為π，縱軸表示歸一化后的幅值大小，無量綱。由圖3至圖8可以看出，總體上來說，在三種定義方式下，故障狀態的三階累量切片譜與正常狀態相比，其譜峰的分布要更復雜一些，同時無論在正常狀態還是故障狀態，定義2和定義3方式下的三階累量切片譜譜峰的分布都比定義1復雜些，這應該是由于前述的不同的定義方式包含了不同的耦合信息的原因。

圖6 定義2方式故障狀態三階累量切片譜圖

圖7 定義3方式正常狀態三階累量切片譜圖

圖8 定義3方式故障狀態三階累量切片譜圖

9 結論與展望

復數三階累量的不同定義形式決定了其包含了不同的耦合信息，這種不同的耦合信息必然反應到由其得出的復數三階累量切片中。本文通過故障診斷實驗對這種不同的耦合特征進行分析，得出了具有不同耦合性質的復數三階累量切片其保持同類信號特征的能力是不一樣的這一結論，也為故障診斷提供了一種可以選擇的方法。本文提供的故障診斷方法可以通過單片機系統予以實現，在系統中可以先存儲一定數量的數據，再根據實時采集到的數據進行故障診斷。

[1]Kocur D，Stanko R.Order bispectrum：a new tool for reciprocated machinecondition monitoring[J].MechanicalSystems and Signal Processing，2000，14（6）：87l-890.

[2]曹精明，邵忍平，胡文濤.HOC與 EMD結合的齒輪損傷檢測研究[J].儀器儀表學報，2011，32（4）：729-735.

[3]邵忍平，黃欣娜.基于高階累積量的齒輪系統故障檢測與診斷[J].機械工程學報，2008，44（6）.

[4]蔡奇志，黃宜堅.自回歸三譜切片研究[J].儀器儀表學報，2009，30（2）.

[5]Hinich M J.Higher-order cumulants and cumulant spectra[J]. Circuits Systems and Signal Processing，1994，13（4）：391-420.

[6]Raghuveer M R.Τime-domain approaches to quadratic phase coupling estimation[J].IEEE Τransactions on Automatic Control，1990，35（1）.

[7]張嚴，王樹勛.非線性相位耦合的切片譜分析方法[J].電子學報，1998，26（10）：104-109.

[9]吳正國，夏立，尹為民.現代信號處理技術：高階譜、時頻分析與小波變換[M].武漢：武漢大學出版社，2002：21-22.

[10]吳文兵，黃宜堅.基于雙相干譜的減壓閥故障診斷[J].計算機測量與控制，2011（10）.

WU Wenbing1，2,HUANG Yijian3

1.Fuzhou Τechnical College of Foreign Studies,Fuzhou 350018,China
2.School of Information Science and Engineering,Xiamen University,Xiamen,Fujian 361000,China
3.College of Mechanical Engineering,Huaqiao University,Quanzhou,Fujian 362021,China

Complex third-order cumulant has different definition forms.Different forms have different coupling properties,and the generated complex cumulants slices contain different coupling information of signals.In experiments,using the different definitions,the same coupling method is applied to both specific fault signals and normal signals.Furthermore,complex third-order cumulant slices spectrum is defined,and it is used to analyse the coupling features of normal signals and fault signals.Experiments indicate,the detection accuracy rate on the same fault are not the same with the different coupling methods,thus,it provides an alternative method to diagnose the specific fault.

complex third-order cumulant slices;coupling;fault diagnosis;spectrum of slices

復數三階累量具有不同的定義形式，定義形式不同，其耦合性質也不相同，則由其產生的復數累量切片也包含了信號間不同的耦合信息，實驗中利用這些不同的定義方式，對特定的故障信號和正常信號同時采用相同的耦合方式進行故障判別。同時定義了復數三階累量切片譜，通過切片譜對正常信號和故障信號的耦合特征進行了分析。實驗結果表明，隨著耦合方式的不同，同一故障的診斷正確率也不相同，從而為特定的故障判別提供了一種可供選擇的方法。

復數三階累量切片；耦合；故障診斷；切片譜

A

ΤG156

10.3778/j.issn.1002-8331.1112-0267

WU Wenbing,HUANG Yijian.Application of complex cumulants slices in mechanical fault diagnosis.Computer Engineering and Applications,2013,49（15）：238-241.

國家自然科學基金（No.50975098）；福建省教育廳A類科技項目（No.JA11342）。

吳文兵，男，在讀博士，講師，主要研究方向為信號檢測、故障診斷。

2011-12-14

2012-04-16

1002-8331（2013）15-0238-04

CNKI出版日期：2012-05-31 http://www.cnki.net/kcms/detail/11.2127.ΤP.20120531.1542.003.html