王 偉

（山西省煤炭職業(yè)中等專業(yè)學(xué)校，山西太原 030032）

0 引言

齒輪系統(tǒng)是機(jī)械設(shè)備的主要零部件，應(yīng)用十分廣泛[1]，但復(fù)雜的齒輪傳動系統(tǒng)往往包含著行星齒輪和定軸齒輪，由于工作環(huán)境的不確定性和各種振動的影響，故障信號特征調(diào)制現(xiàn)象嚴(yán)重，增加了信號分離和故障診斷的難度[2]，而且定軸齒輪故障特征往往會被淹沒在行星齒輪的信號特征中。目前對復(fù)雜齒輪系統(tǒng)的多故障診斷研究尚少。

目前，分解故障信號的常用辦法有局部均值分解法、EMD（Empirical Mode Decomposition，經(jīng)驗(yàn)?zāi)B(tài)分解法）和小波變換法等[3]。徐先峰等人通過對軸承故障信號分解，實(shí)現(xiàn)對軸承故障的有效診斷以及故障程度的基本判斷[4]。郭向玲等人利用EMD信號分解法對汽輪機(jī)轉(zhuǎn)子的故障進(jìn)行診斷，結(jié)果顯示，該方法存在故障信號重疊、信號包絡(luò)不充分等缺陷，不利于故障信號的提取[5]，而且小波變換需要提前選定小波基值和分解層數(shù)，自適應(yīng)力弱；局部均值分解法運(yùn)算速度慢、信號沖突等問題無法避免[6]。

因此本文針對復(fù)雜行星齒輪系統(tǒng)的多故障診斷問題，基于不同齒輪副嚙合頻率的差異，引入VMD（Variational Mode Decomposition，變分模態(tài)信號分解法）實(shí)現(xiàn)了行星齒輪和定軸齒輪信號特征的分離，并分別對其進(jìn)行包絡(luò)解調(diào)，以實(shí)現(xiàn)對不同齒輪特征信號的提取。

VDM[7-8]是一種自適應(yīng)能力很強(qiáng)的信號分解方法，相對于其他信號分解方法，該方法分解精度高、收斂速度快等優(yōu)點(diǎn)[9]，適用于復(fù)雜齒輪系統(tǒng)故障的信號故障處理。

1 變分模態(tài)分解原理

設(shè)VMD 分解得到的分量為uk，即：

式中

Ak（t）——uk的幅值

φk（t）——uk的瞬時頻率

VMD 具體迭代運(yùn)算過程如下：

（2）n=n+1，開始計(jì)算。

（6）重復(fù)步驟（2）～（5），直到達(dá)到迭代停止條件Σk，循環(huán)結(jié)束，得到和ωk；通過傅里葉逆變換計(jì)算uk。

2 仿真分析

2.1 仿真信號構(gòu)建

為了驗(yàn)證所提方法的有效性，本文模擬了三級齒輪的振動信號，模擬參數(shù)見表1。在其中設(shè)有局部故障，并添加了一個高斯信號N（t）作為噪聲干擾。

表1 仿真參數(shù)Hz

模擬信號模型見式（2）。

xs（t）——正常齒輪嚙合頻率

N（t）——噪聲信號

齒輪信號的傳遞路徑如圖1 所示。

圖1 行星齒輪箱振動信號的傳輸路徑

2.2 故障齒輪的特征頻率

2.2.1 行星齒輪

齒輪系統(tǒng)中行星齒輪存在故障時，故障的特征頻率與齒輪之間的嚙合次數(shù)有關(guān)。

在行星齒輪系統(tǒng)中，行星架轉(zhuǎn)動一周太陽輪時和所有的行星齒輪嚙合一次，故太陽輪的故障頻率×N。行星架轉(zhuǎn)動一周，行星輪與太陽輪和齒圈各發(fā)生一次嚙合，所以行星輪的故障特征頻率為。

2.2.2 固定軸齒輪故障

定軸齒輪故障的特征頻率相對行星齒輪要簡單，當(dāng)局部產(chǎn)生故障時定軸齒輪特征頻率和齒輪的轉(zhuǎn)動頻率相同，所以故障特征頻率為fl=fr。

2.3 信號處理流程

VMD 和包絡(luò)解調(diào)的多故障診斷法，根據(jù)多級傳動時嚙合頻率不同的特點(diǎn)，先將振動信號按不同頻率段分解，完成行星齒輪和直齒齒輪信號特征的分離，其次對不同的信號分量進(jìn)行包絡(luò)解調(diào)提取故障特征。具體流程如圖2 所示。

圖2 故障診斷流程

（1）首先采集行星齒輪箱上的振動信號，對原始信號濾波預(yù)處理，減小高頻諧波的影響。

（2）對濾波后的信號VMD 處理提取，獲得特征頻率的不同分量。

（3）對各頻率分量進(jìn)行包絡(luò)處理，得到各自的包絡(luò)譜。

（4）對包絡(luò)譜進(jìn)行特征分析，得到齒輪系統(tǒng)的故障信息。

2.4 仿真分析結(jié)果

圖3～圖5 為模擬信號的處理結(jié)果圖、VMD 分解信號的時域圖和時頻圖。

圖3 模擬信號的處理結(jié)果

齒輪系統(tǒng)的復(fù)雜性和齒輪轉(zhuǎn)動過程中信號重疊等問題，導(dǎo)致當(dāng)多個齒輪同時出現(xiàn)問題時信號互相干擾，無法有效提取齒輪故障信息。通過對齒輪系統(tǒng)的信號行星包絡(luò)處理，可以在圖3c）中得到行星架的轉(zhuǎn)動頻率fc和太陽輪故障的特征頻率fc等。但是定軸齒輪的故障信號無法在包絡(luò)圖中得到，故對模擬信號通過VMD 分解得到圖4。由圖4 可以看出，直齒輪和太陽輪的故障特征被淹沒在VIMF 中，所以分別對VIMF1 和VIMF3 做包絡(luò)譜分析（圖5）。由圖5 可以清楚地看到齒輪所對應(yīng)的故障特征，由此可證該方法對于提取多級行星齒輪箱多齒輪故障的有效性。

圖4 通過VMD 分解模擬信號的時域圖和時頻圖

圖5 VIMF1 和VIMF3 的包絡(luò)譜分析

3 實(shí)驗(yàn)分析

3.1 實(shí)驗(yàn)設(shè)置

實(shí)驗(yàn)所使用的行星齒輪箱試驗(yàn)臺如圖6 所示，模擬太陽輪和直齒輪同時斷齒時情況。

圖6 多級行星齒輪箱試驗(yàn)臺

在該實(shí)驗(yàn)中，人為設(shè)置斷齒的故障齒輪。圖7 為實(shí)驗(yàn)用不同太陽輪和定軸齒輪照片。

圖7 太陽輪和定軸齒輪

3.2 實(shí)驗(yàn)結(jié)果分析

圖8～圖10 為當(dāng)行星齒輪的太陽輪和二級定軸直齒輪的主動輪同時出現(xiàn)斷齒故障時的振動信號分析結(jié)果。

由圖8b）可知，定軸齒輪斷齒會在直齒輪1 倍到4 倍嚙合頻率間產(chǎn)生連續(xù)的峰值，但是信號特征并不明顯，而且出現(xiàn)了太陽輪故障頻率和嚙合頻率疊加的復(fù)合頻率（fm1±4fs），這是齒輪在實(shí)際制造過程中由于制造誤差導(dǎo)致的。由圖8c）可以看到太陽輪故障特征頻率4fs及其相關(guān)頻率，但不影響定軸齒輪的故障特征信號。

圖8 模擬信號的處理結(jié)果

分析VMD 分解后的特征頻率信號，由圖9 可看出不同齒輪的振動信號被分解開。從圖10a）可以明顯看到定軸直齒輪的故障特征頻率fr1，便于定軸直齒輪的故障分析。在圖10b）中可以清楚地觀察到太陽輪的故障特征頻率4fs以及一些相關(guān)頻率信息，表示太陽輪的斷齒故障。

圖9 通過VMD 分解模擬信號的時域圖和時頻圖

圖10 VIMF1 和VIMF3 的包絡(luò)譜分析

利用VMD 分解振動信號，使得振動信號按照頻段分解。便于對不同齒輪的故障特性進(jìn)行提取和分析，有助于齒輪故障的判斷。通過實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證所提出的方法的有效性。

4 結(jié)論

針對復(fù)雜齒輪系統(tǒng)發(fā)生多故障時特征信號的提取和故障診斷困難等問題，提出了一種基于VMD 的信號分解方法和特征頻率的解調(diào)方法來識別齒輪系統(tǒng)的故障。該方法首先用VMD 將信號分解得到不同齒輪的振動信號，然后對所需要的信號包絡(luò)計(jì)算得到包絡(luò)譜圖，通過不同特征頻率的峰值識別不同的齒輪故障。本研究結(jié)合仿真分析和實(shí)驗(yàn)驗(yàn)證，證明了該方法的可行性。

（1）結(jié)合VMD 與包絡(luò)解調(diào)法，通過VMD 方法分離不同的故障特征，再用包絡(luò)解調(diào)法可以實(shí)現(xiàn)故障特征的有效提取。

（2）建立多級齒輪系統(tǒng)存在故障時的振動信號模型，并對其信號進(jìn)行仿真分析和實(shí)驗(yàn)驗(yàn)證，證明了該方法的有效性。