楊杰，鄭海起，劉曉平，王彥剛

(石家莊軍械工程學(xué)院，石家莊 050003)

滾動(dòng)軸承存在局部缺陷時(shí)，其振動(dòng)信號(hào)中的

脈沖信號(hào)含有豐富的缺陷信息。若能夠有效地將缺陷引起的脈沖信號(hào)提取出來(lái)，便可診斷出缺陷存在的部位[1]。

形態(tài)學(xué)濾波器作為一種非線性濾波器，可以有效地提取信號(hào)的邊緣輪廓以及形狀特征，克服了線性濾波的不足。目前，形態(tài)濾波在旋轉(zhuǎn)機(jī)械振動(dòng)信號(hào)濾波方面已經(jīng)有一定的應(yīng)用[2]，但是經(jīng)典的形態(tài)學(xué)濾波在取單一結(jié)構(gòu)元素的情況下存在缺陷[3]，因此各種針對(duì)結(jié)構(gòu)元素的改進(jìn)算法相繼出現(xiàn)，如基于權(quán)重的形態(tài)學(xué)運(yùn)算[4]、多尺度結(jié)構(gòu)元[5]和結(jié)構(gòu)元的自適應(yīng)優(yōu)化算法[6]等。

文獻(xiàn)[7]提出了柔性形態(tài)學(xué)(Soft Mathematical Morphology)的概念，研究表明柔性形態(tài)學(xué)對(duì)加性噪聲和信號(hào)的微小變化不敏感，能同時(shí)去除信號(hào)中的正負(fù)噪聲，并保留細(xì)節(jié)。文獻(xiàn)[8]比較了柔性形態(tài)濾波(Soft Morphological Filter，SMF)和小波包的濾波效果，表明SMF在診斷外圈故障和滾動(dòng)體故障時(shí)的優(yōu)越性。

在此提出了一種多尺度柔性形態(tài)濾波(Multi-scale Soft Morphological Filter，MSMF)方法，并將其應(yīng)用于軸承故障診斷。

1 數(shù)學(xué)形態(tài)濾波

1.1 基本形態(tài)運(yùn)算

數(shù)學(xué)形態(tài)學(xué)的基本運(yùn)算包括腐蝕、膨脹、開(kāi)運(yùn)算和閉運(yùn)算。設(shè)一維信號(hào)f(n)為定義在F=(0,1,…,N-1)上的離散函數(shù)，結(jié)構(gòu)元素g(n)為定義在G=(0,1,…,M-1)上的離散函數(shù)，且N≥M，則g(n)對(duì)f(n)的腐蝕和膨脹運(yùn)算分別定義為

(fΘg)(n)=min[f(n+m)-g(m)],

(1)

(f⊕g)(n)=max[f(n-m)+g(m)],

(2)

式中：m∈(0,1,…,M-1)。

g(n)對(duì)f(n)的開(kāi)運(yùn)算和閉運(yùn)算分別定義為

(f°g)(n)=(fΘg⊕g)(n),

(3)

(f·g)(n)=(f⊕gΘg)(n)。

(4)

通常開(kāi)、閉運(yùn)算用于構(gòu)成各種形態(tài)濾波器，而它們本身就是最基本的形態(tài)濾波器，且以不同的方式平滑信號(hào)。開(kāi)運(yùn)算可抑制信號(hào)中的正噪聲；相反，閉運(yùn)算可抑制信號(hào)中的負(fù)噪聲。

1.2 柔性形態(tài)濾波

在柔性數(shù)學(xué)形態(tài)學(xué)中，結(jié)構(gòu)元素被分割成硬核和柔性邊緣兩部分，經(jīng)典形態(tài)算子中的最大最小運(yùn)算在柔性形態(tài)算子中被排序統(tǒng)計(jì)所代替，且最主要的思想是引入了定義一種重復(fù)集。

柔性形態(tài)學(xué)的結(jié)構(gòu)系統(tǒng)由[B，A，k]3個(gè)參數(shù)組成，有限集A和B為定義在Z上的凸集，k滿足于1

{k<>f(a)}={f(a),f(a),…,f(a)}，

式中：k為正整數(shù)，且1≤k≤min{|B|/2,|BA|}，a∈B。

[B，A，k]對(duì)f的腐蝕和膨脹分別定義為

fΘ[B,A,k](x)=min(k){k<>f(a):a∈Ax}∪

{f(b):b∈(BA)x},

(5)

f⊕[B，A，k](x)=max(k){k<>f(a):a∈Ax}∪

{f(b):b∈(BA)x},

(6)

其中，min(k)表示取第k小的值，同理max(k)表示取第k大的值。當(dāng)k=1時(shí)，柔性形態(tài)變換與經(jīng)典形態(tài)變換相同；當(dāng)|A|=0時(shí)，柔性形態(tài)變換與順序形態(tài)變換相同。

設(shè)A，B的映像集分別為AS，BS，則結(jié)構(gòu)元素[B，A，k]對(duì)f進(jìn)行形態(tài)開(kāi)和形態(tài)閉運(yùn)算定義如下

f°[B,A,k]=(fΘ[B,A,k])⊕[BS,AS,k],

(7)

f·[B,A,k]=(f⊕[B,A,k])Θ[BS,AS,k]。

(8)

由以上定義可知，在柔性形態(tài)學(xué)運(yùn)算中，落入到硬核區(qū)域中的信號(hào)數(shù)據(jù)權(quán)重較大，而落入到軟核區(qū)域中的信號(hào)數(shù)據(jù)權(quán)重較小。因此，和標(biāo)準(zhǔn)形態(tài)學(xué)運(yùn)算比較，柔性形態(tài)學(xué)運(yùn)算具有更好的保持細(xì)節(jié)和抗噪聲功能[9]。

2 多尺度柔性形態(tài)濾波

2.1 多尺度形態(tài)學(xué)

多尺度運(yùn)算是一個(gè)由粗到細(xì)的多層次描述和處理方法，如同人的感知過(guò)程，先通過(guò)對(duì)大范圍特征進(jìn)行粗分析，再通過(guò)捕捉細(xì)節(jié)使分析越來(lái)越精細(xì)，最后獲得對(duì)感知對(duì)象的準(zhǔn)確認(rèn)識(shí)。在具體實(shí)施過(guò)程中，其采用不同尺寸的結(jié)構(gòu)元素對(duì)信號(hào)進(jìn)行變換。

假設(shè)給定了形態(tài)運(yùn)算T和信號(hào)X，則基于T的多尺度形態(tài)學(xué)運(yùn)算可定義為一族形態(tài)學(xué)變換{Tλ|λ>0}，其中Tλ定義為

Tλ(X)=λT(X/λ),λ>0,

(9)

式中：λ>0為刻度，即可推導(dǎo)出多尺度膨脹、腐蝕、開(kāi)和閉運(yùn)算分別為

(fΘg)λ(n)=fΘλg1/λ(n),

(10)

(f⊕g)λ=f⊕λg1/λ(n)，

(11)

(f°g)λ=f°λg1/λ(n)，

(12)

(f·g)λ=f·λg1/λ(n)，

(13)

式中：g1/λ(n)=g(n/λ)。若結(jié)構(gòu)函數(shù)取凸函數(shù)，隨著刻度的增大，多刻度運(yùn)算會(huì)濾去信號(hào)更大的變化，使信號(hào)越來(lái)越簡(jiǎn)單。

2.2 多尺度柔性形態(tài)濾波

將柔性形態(tài)學(xué)中的硬核A進(jìn)行多尺度化，即實(shí)現(xiàn)了多尺度柔性形態(tài)運(yùn)算，其膨脹和腐蝕分別定義為

{fΘ[B,A,k]}λ(x)=min(k){k<>λf(a/λ):a∈Ax}∪{f(b):b∈(BA)x},

(14)

f⊕[B,A,k](x)=max(k){k<>λf(a/λ):a∈Ax}∪{f(b):b∈(BA)x}。

(15)

由此可推導(dǎo)出開(kāi)運(yùn)算和閉運(yùn)算的表達(dá)式，這里不再贅述。MSMF對(duì)SMF中的核結(jié)構(gòu)元素多尺度化使得形態(tài)運(yùn)算更加靈活，從而達(dá)到進(jìn)一步改善形態(tài)濾波器性能的目的。

3 在軸承故障診斷中的應(yīng)用

在某型單級(jí)齒輪箱上進(jìn)行試驗(yàn)驗(yàn)證。在該系統(tǒng)中，由電動(dòng)機(jī)帶動(dòng)輸入軸，輸出軸帶動(dòng)負(fù)載。主動(dòng)齒輪齒數(shù)Z1為30，被動(dòng)齒輪齒數(shù)Z2為50，輸入端滾動(dòng)軸承型號(hào)為6206，局部損傷故障分別在鋼球、內(nèi)圈和外圈表面采用電火花機(jī)加工制作，損傷直徑為0.178 mm。試驗(yàn)時(shí)采樣頻率fs為12 800 Hz，采樣時(shí)長(zhǎng)為1 s，轉(zhuǎn)速為1 136 r/min，計(jì)算可得內(nèi)圈故障的特征頻率為102.5 Hz，外圈故障的特征頻率為67.9 Hz，鋼球故障的特征頻率為89.2 Hz。

測(cè)得軸承故障信號(hào)的時(shí)域波形如圖 1所示。圖 2是對(duì)原始信號(hào)直接進(jìn)行Fourier變換得到的頻譜圖，由圖 2可知，軸承故障特征頻率完全被噪聲所淹沒(méi)，無(wú)法診斷出任何故障。

圖1 軸承3種故障的原始信號(hào)

圖2 原始信號(hào)的頻譜圖

選擇結(jié)構(gòu)元素B的長(zhǎng)度為20，A的長(zhǎng)度為7，k為2，λ為2，取MSMF開(kāi)和閉運(yùn)算的平均值，濾波后的時(shí)域如圖3所示，頻譜如圖4所示。由圖4可知，經(jīng)過(guò)MSMF濾波后，轉(zhuǎn)頻及內(nèi)圈、外圈和鋼球的故障特征頻率及其倍頻清晰可見(jiàn)，從而能診斷出軸承的3種故障。

圖3 多尺度柔性形態(tài)濾波時(shí)域圖

圖4 多尺度柔性形態(tài)濾波頻譜圖

4 與柔性形態(tài)濾波的比較

文獻(xiàn)[1]給出了形態(tài)濾波解調(diào)和Hilbert變換解調(diào)比較結(jié)果，表明形態(tài)濾波解調(diào)和Hilbert變換解調(diào)相比能獲得更好的信噪比，在此不再贅述。下文以外圈損傷為例對(duì)多尺度柔性形態(tài)濾波和柔性形態(tài)濾波的效果進(jìn)行比較，對(duì)圖 1b所示軸承外圈損傷信號(hào)進(jìn)行柔性形態(tài)濾波，A，B，k的取值同上，取開(kāi)運(yùn)算和閉運(yùn)算的平均值，濾波后的時(shí)域和頻域如圖5所示。

圖5 柔性形態(tài)濾波效果圖

由圖5可知，柔性形態(tài)濾波同樣能達(dá)到濾除噪聲的目的，但是濾波后信號(hào)的幅值明顯降低，說(shuō)明柔性形態(tài)濾波在濾除噪聲的同時(shí)也濾除了部分的沖擊脈沖，且濾波后仍有較大噪聲，這為進(jìn)一步的故障診斷帶來(lái)不便，相比之下，多尺度柔性形態(tài)濾波在濾除噪聲的同時(shí)能更好地保留原始信號(hào)中的沖擊成分，因此其濾波效果優(yōu)于柔性形態(tài)濾波的濾波效果。

5 結(jié)束語(yǔ)

數(shù)學(xué)形態(tài)濾波作為一種很好的非線性濾波方法在機(jī)械故障診斷領(lǐng)域獲得了越來(lái)越廣泛的應(yīng)用，但是經(jīng)典的形態(tài)學(xué)濾波由于取單一結(jié)構(gòu)元素，故在濾波效果上存在很大不足，常常會(huì)濾除對(duì)診斷具有很大價(jià)值的脈沖信號(hào)。柔性形態(tài)學(xué)在魯棒性和保留細(xì)節(jié)方面對(duì)經(jīng)典形態(tài)學(xué)做了改進(jìn)，在此基礎(chǔ)上，提出了一種基于多尺度的柔性形態(tài)濾波算法，該算法進(jìn)一步優(yōu)化了柔性形態(tài)運(yùn)算中的結(jié)構(gòu)元素，目的是更好地提取軸承故障信號(hào)的沖擊脈沖。試驗(yàn)表明，多尺度柔性形態(tài)濾波在獲得較高信噪比的同時(shí)，比經(jīng)典的柔性形態(tài)濾波能夠更好地保留故障軸承原始信號(hào)中的沖擊成分。

除了形態(tài)算子的設(shè)計(jì)影響濾波效果外，結(jié)構(gòu)元素、k和λ值的選擇對(duì)結(jié)果的影響也較大，針對(duì)軸承振動(dòng)信號(hào)如何選取適當(dāng)?shù)慕Y(jié)構(gòu)元素及其內(nèi)核的大小以及k，λ值，以獲得更好的處理效果，還有待進(jìn)一步研究。