徐亞南， 徐玉秀

（1.天津工業(yè)大學(xué) 機械工程學(xué)院， 天津 300387； 2.天津工業(yè)大學(xué) 現(xiàn)代機電裝備技術(shù)重點實驗室，天津 300387）

引言

由于風(fēng)機齒輪箱工況復(fù)雜，易發(fā)生兩種及以上的耦合故障，其實驗信號中表現(xiàn)出了區(qū)別于兩種單獨故障時的故障特征。常用的時頻分析方法難以滿足對耦合故障中行星齒輪故障特征的識別。

在進(jìn)行齒輪箱故障識別與診斷工作時，國內(nèi)外學(xué)者大多從齒輪箱傳動系統(tǒng)動力學(xué)建模、箱體表面振動響應(yīng)的信號處理和齒輪箱系統(tǒng)傳遞特性研究三方面入手，通過建立剛?cè)狁詈蟿恿W(xué)模型，求解得到了內(nèi)部動態(tài)激勵信號[1]，推導(dǎo)出行星齒輪箱振動信號方程式，并進(jìn)行相應(yīng)的頻譜分析[2-3]。在基于行星齒輪箱傳遞機理上，建立了單級行星齒輪齒輪傳動系統(tǒng)的振動信號模型。對于行星齒輪公轉(zhuǎn)調(diào)制作用而產(chǎn)生的動態(tài)信號[4]，建立信號傳播模型，應(yīng)用最小二乘法類似方法，對信號傳播模型進(jìn)行參數(shù)估計，得到單一故障、耦合故障狀態(tài)下多信號傳播路徑的傳遞函數(shù)。

1 風(fēng)機齒輪箱模擬實驗

由于風(fēng)電齒輪箱信號采集存在一定困難，所以采用圖1所示的模擬試驗臺作為研究對象，1為電機，2為轉(zhuǎn)頻調(diào)節(jié)器，3為定軸齒輪箱，4為測點，5為行星齒輪箱，6為負(fù)載加載器。通過對試驗臺的研究，能模擬出正常、故障狀態(tài)下的運行狀態(tài)。

圖1 齒輪故障模擬試驗臺

通過以故障齒輪替代正常齒輪的方法，分別實驗并測取正常狀態(tài)、行星齒輪缺齒狀態(tài)、行星齒輪均勻磨損狀態(tài)、行星齒輪缺齒-行星齒輪磨損的耦合狀態(tài)下齒輪箱箱體表面加速度響應(yīng)信號。試驗時，輸入軸轉(zhuǎn)頻fd=40 Hz，輸出軸加載電壓為10 V，加速度傳感器靈敏度為100.2 mV/（m·s-2）。通過計算可知該轉(zhuǎn)頻下各級嚙合頻率分別為高速級嚙合頻率f1=1160 Hz、中速級嚙合頻率f2=417 Hz、低速級嚙合頻率f3=101.5 Hz。

通過觀察圖2中正常信號2-1與下頁故障信號2-2、2-3、2-4可知，當(dāng)行星齒輪發(fā)生缺齒故障時，引起行星齒輪嚙合頻率101.5 Hz幅值的升高，同時280 Hz幅值也發(fā)生了升高。當(dāng)發(fā)生行星齒輪磨損故障時，除406 Hz被激起外，1200 Hz也升高明顯。當(dāng)發(fā)生兩種故障同時存在的耦合故障時，除40 Hz、101.5 Hz、406 Hz以及1200 Hz的增高外，還存在異于兩種單獨故障時的故障頻率460 Hz的發(fā)生。由于故障并不是線性疊加，該故障頻率的產(chǎn)生無法通過正常頻譜分析解釋，因此引入傳遞函數(shù)進(jìn)行分析。

2 齒輪故障模擬試驗臺仿真模型

由于故障齒輪模擬實驗只能獲取箱體表面的振動信號，對于箱體內(nèi)部齒輪相互嚙合時的接觸力信號無法直接通過傳感器測量獲取，因此需要通過動力學(xué)仿真的方式對風(fēng)機傳動系統(tǒng)進(jìn)行建模。為滿足建模精度要求，采用建模軟件UG進(jìn)行試驗臺建模，通過ADAMS/Exchange接口與專業(yè)的動力學(xué)仿真軟件ADAMS進(jìn)行數(shù)據(jù)轉(zhuǎn)換。在建模時，為了提高精度，將定軸齒輪箱及行星齒輪箱與軸的接觸點設(shè)為剛性區(qū)域，其他部分為柔性區(qū)域。

對于建立好的傳動系統(tǒng)剛?cè)狁詈夏Ｐ停瑢ζ溥M(jìn)行Dynamic仿真，提取正常狀態(tài)、兩種單故障狀態(tài)以及行星齒輪缺齒-行星齒輪均勻磨損耦合故障狀態(tài)下齒輪箱中高速級嚙合、中速級嚙合、低速級內(nèi)外嚙合的接觸力模擬信號。

圖2 各工況下實驗信號

3 傳遞函數(shù)及ARX模型選擇求解

3.1 傳遞函數(shù)的定義

對于實際的振動系統(tǒng)來說，建立系統(tǒng)的多自由度振動微分方程，得：

式中，[Kg]為復(fù)剛度矩陣。并對公式（1）進(jìn)行拉氏變換，得：

式中，[Z（s）]為系統(tǒng)的動剛度。在一定的激勵作用下，動剛度與系統(tǒng)的響應(yīng){X（s）}成反比，即動剛度會阻礙的系統(tǒng)振動，故又稱[Z（s）]為系統(tǒng)的機械阻抗，其倒數(shù)稱為傳遞函數(shù)，用[H（s）]表示，即：

由公式（1）得：

對于實際的振動系統(tǒng)來說，用ω代替s，不會失去有用的信息，故對公式（2）兩邊進(jìn)行傅氏變換，得：

3.2 傳遞函數(shù)ARX模型的選擇及求解

對于給定的系統(tǒng)，可以通過選擇時域響應(yīng)或頻域響應(yīng)，對系統(tǒng)進(jìn)行辨識。由于連續(xù)的系統(tǒng)傳遞模型在擬合時存在解不唯一或擬合結(jié)果不穩(wěn)定等缺點，通常不對連續(xù)模型直接辨識，更多的是將連續(xù)的系統(tǒng)模型進(jìn)行離散化處理，通過系統(tǒng)辨識獲得離散的傳遞函數(shù)后，再對離散模型進(jìn)行連續(xù)化處理[5]。

離散系統(tǒng)傳遞函數(shù)為：

它對應(yīng)的差分方程為：

式中：ε（t）為殘差信號，這里為方便起見，輸出信號簡記為 y（t），且用 y（t-1）表示輸出信號 y（t）在前一個采樣周期處的函數(shù)值，這種模型又稱為自回歸歷遍模型（ARX模型）。

對于多輸入單輸出ARX模型，其對應(yīng)的差分方程為：

式中：d 為各個延遲構(gòu)成的矩陣；A（z-1）和 B（z-1）均為p×q多項式矩陣，且

即拉氏域內(nèi)的傳遞函數(shù)變?yōu)楦凳嫌騼?nèi)的頻響函數(shù)。傳遞函數(shù)由系統(tǒng)本身所決定的，與輸入量沒有關(guān)系。當(dāng)獲取系統(tǒng)的傳遞函數(shù)之后，可以根據(jù)輸入量獲取系統(tǒng)的輸出量，或者可以根據(jù)所需的輸出量來求取輸入量。通過頻響函數(shù)[H（ω）]可充分描述系統(tǒng)的動態(tài)特性。

系統(tǒng)辨識常用的軟件為MATLAB，其優(yōu)點在于內(nèi)置辨識工具箱中包含各類系統(tǒng)辨識模型，并且擁有強大的計算能力。本文所運用的ARX模型辨識系統(tǒng)可以通過函數(shù)arx（）實現(xiàn)，在確定了分子多項式階次m-1，分母多項式階次n以及系統(tǒng)純滯后d之后，使用arx（）函數(shù)可以直接辨識出系統(tǒng)的Ai和Bi矩陣。

4 故障特征傳遞特性分析

4.1 傳遞函數(shù)分析

以正常狀態(tài)、行星齒輪缺齒故障狀態(tài)、行星齒輪均勻磨損故障狀態(tài)、行星齒輪缺齒-行星齒輪均勻磨損耦合故障狀態(tài)的嚙合接觸力為輸入信號，各種狀態(tài)的實驗信號為輸出信號，分別建立以上四種狀態(tài)下的傳遞函數(shù)，如圖3所示。

比較圖3-2和3-1，當(dāng)發(fā)生行星齒輪缺齒故障時，高速軸轉(zhuǎn)頻fd（40 Hz）幅值增加，同時未知頻率280 Hz幅值也發(fā)生升高。

圖3 各狀態(tài)下的傳遞函數(shù)

比較圖3-3和3-1，當(dāng)發(fā)生行星齒輪單磨損故障時，未知頻率406 Hz以及1200 Hz幅值增大。通過簡單計算406 Hz恰好為低速級嚙合頻率f3（101.5 Hz）的四倍頻。

比較圖3-4和3-1，當(dāng)發(fā)生行星齒輪缺齒-行星齒輪均勻磨損耦合故障時，主要的故障特征頻率有 40 Hz、406 Hz、460 Hz和 1200 Hz。

傳遞函數(shù)中表現(xiàn)出了與頻譜相一致的特征頻率，驗證了傳遞函數(shù)的正確性。通過以上分析，280 Hz為行星齒輪缺齒所特有的故障頻率，406 Hz為均勻磨損所特有的故障特征頻率，460 Hz為耦合故障所特有的故障頻率，因此針對這幾個頻率，對其貢獻(xiàn)量進(jìn)行詳細(xì)分析。

4.2 貢獻(xiàn)量分析

為了更清楚地描述傳遞函數(shù)中各頻率的產(chǎn)生原因，分別對各故嚙合接觸力在故障頻率處的貢獻(xiàn)量進(jìn)行分析，如圖4所示。

圖4 不同狀態(tài)下各嚙合接觸力對特定頻率貢獻(xiàn)量統(tǒng)計

通過圖4-1可知，在280 Hz處，高速級嚙合接觸力貢獻(xiàn)量由正常狀態(tài)下的33%增加到66.3%，缺齒故障引起了高速級嚙合接觸力傳遞路徑的改變。

通過圖4-2，在406 Hz處，低速級內(nèi)嚙合接觸力和外嚙合接觸力分別由正常狀態(tài)下的26.1%和22.4%升高到37.1%和31.6%，即由于行星齒輪均勻磨損故障，引起了低速級內(nèi)、外嚙合接觸力傳遞路徑的改變。

通過圖4-3可知，當(dāng)發(fā)生耦合故障時，在460 Hz處高速級嚙合接觸力由正常狀態(tài)下的28.2%升高到45.9%，即由于耦合故障的發(fā)生，導(dǎo)致高速級嚙合接觸力傳遞路徑的改變。對于耦合故障，由傳遞函數(shù)不難發(fā)現(xiàn)，在400 Hz到500 Hz左右的頻率范圍內(nèi)，高速級嚙合所對應(yīng)的接觸力對傳遞路徑的貢獻(xiàn)量最大，說明行星齒輪耦合故障發(fā)生時對高速級嚙合振動造成了一定影響，引發(fā)了460 Hz幅值的增大。

5 結(jié)語

通過ARX模型，建立了正常狀態(tài)、行星齒輪缺齒、行星齒輪均勻磨損和行星齒輪缺齒-行星齒輪均勻磨損耦合故障狀態(tài)下的傳遞函數(shù)。通過對兩種單故障的分析，引出了耦合故障中的連帶故障特征頻率460 Hz。通過對各類故障傳遞函數(shù)中各嚙合接觸力幅值及貢獻(xiàn)量的變化情況，總結(jié)了耦合故障的傳遞特性，解決了460 Hz的來源問題，實現(xiàn)了對行星齒輪缺齒-行星齒輪均勻磨損該耦合故障的診斷。

[1]徐玉秀，趙曉菲，熊一奇.基于傳遞路徑的多級齒輪箱齒輪裂紋故障識別[J].儀器儀表學(xué)報，2016（5）：1018-1024.

[2]Inalpolat M,Kahraman A.A theoretical and experimental investigation of modulation sidebands of planetary gear sets[J].Journal ofSound and Vibration,2009,323(4):677-696.

[3]Inalpolat M,Kahraman A.A dynamic model to predict modulation sidebands of a planetary gear set having manufacturing errors[J].Journal ofSound and Vibration,2010,329(4):371-393.

[4]熊一奇,徐玉秀.基于信號傳播模型的行星齒輪缺齒故障識別[J].儀器儀表學(xué)報，2016（2）：249-255.

[5]齊曉慧，田慶民，董海瑞.基于Matlab系統(tǒng)辨識工具箱的系統(tǒng)建模[J].兵工自動化，2006，25（10）：88-90.