基于ANSYS的齒輪根部裂痕故障診斷法*

王 龍，田懷谷，李 森，劉守法，董 鋒

(1. 西京學(xué)院機(jī)械工程學(xué)院，陜西西安710123; 2. 西安熱工研究院有限公司，陜西西安710032)

?

基于ANSYS的齒輪根部裂痕故障診斷法*

王 龍1，田懷谷1，李 森1，劉守法1，董 鋒2

(1. 西京學(xué)院機(jī)械工程學(xué)院，陜西西安710123; 2. 西安熱工研究院有限公司，陜西西安710032)

利用集總參數(shù)法建立了齒輪動力學(xué)模型，并以龍格-庫塔數(shù)值法求解齒輪的振動響應(yīng)。使用ANSYS 軟件建立了無缺陷齒輪組與帶齒根裂紋齒輪組模型，分析了兩種齒輪組在多種轉(zhuǎn)速下的振動特性。研究表明，齒根處裂紋對齒輪嚙合剛度等效彈簧的振動頻率影響不明顯，且隨著齒輪轉(zhuǎn)速增大，帶裂紋齒輪組振幅的均方根值略高于無缺陷齒輪組；帶裂紋齒輪組在最大振幅附近產(chǎn)生邊頻現(xiàn)象，且齒輪轉(zhuǎn)速越高邊頻越明顯，此特征可用來判斷齒輪組是否存在故障。

集總參數(shù)法 有限元分析 動態(tài)響應(yīng) 等效彈簧

0 引言

齒輪廣泛應(yīng)用于工業(yè)機(jī)械、農(nóng)業(yè)機(jī)械、航空航天和車輛的動力傳輸系統(tǒng)中，產(chǎn)生故障的形式有輪齒折斷、疲勞點蝕和磨損等，這些故障會使齒輪產(chǎn)生周期性脈動力，在振動信號上體現(xiàn)為調(diào)制現(xiàn)象，在頻譜圖上會出現(xiàn)一系列的調(diào)制邊頻帶[1-4]。振動的產(chǎn)生會使設(shè)備磨損嚴(yán)重、產(chǎn)生振動和噪音，造成嚴(yán)重經(jīng)濟(jì)損失，因此近年來齒輪動態(tài)特性一直都是工程領(lǐng)域的研究重點。

本文利用集總參數(shù)法和有限元法，對時域和頻域內(nèi)齒輪嚙合振動特點進(jìn)行了深入分析，討論了齒根處有裂紋齒輪的振動特性，以期利用振動信號來判斷齒根處是否出現(xiàn)裂紋。

1 動力學(xué)模型和運(yùn)動方程

圖1 齒輪系統(tǒng)數(shù)學(xué)模型

假設(shè)傳動軸與軸承為剛性，利用集總參數(shù)法，使用等效質(zhì)量、等效阻尼與嚙合剛度等效彈簧的形式表達(dá)系統(tǒng)，忽略嚙合點處的摩擦力與扭矩T1、T2的波動，所建模型如圖1所示[5]，根據(jù)牛頓運(yùn)動定律，由文獻(xiàn)[6]得知齒輪1和2的系統(tǒng)運(yùn)動方程式分別為：

(1)

(2)

其中I1、I2為齒輪的質(zhì)量慣性矩，θ1、θ2為齒輪角位移，c為嚙合阻尼系數(shù)，rb1、rb2為齒輪基圓半徑，Kc(θ1)為嚙合剛度，T1、T2為齒輪扭矩。

等效彈簧的位移為x=rb1θ1-rb2θ2，即為傳遞誤差，將(1)與(2)組合可得等效彈簧的系統(tǒng)運(yùn)動方程為：

(3)

表1 動力學(xué)模型各參數(shù)值

其中y1、y2、y3和y4的初始值由系統(tǒng)初始狀態(tài)決定。

2 齒輪動態(tài)響應(yīng)數(shù)值分析

齒輪為漸開線標(biāo)準(zhǔn)直齒圓柱齒輪，模數(shù)m=6 mm，齒數(shù)z1=23，z2=23，齒寬b=15 mm。在無缺陷齒輪的齒根處切出一條裂紋，裂紋的深度與齒寬比值為η=a/b，裂紋的位置系數(shù)β =vd/h(vd為裂紋與齒根圓距離，h為全齒高)，式中η = 0.1、β=0.11，生成的齒根處有裂紋齒輪模型如圖2所示。齒輪使用的材料為鋼，密度為7 850 kg/m3、楊氏模量為210 GPa、泊松比為0.3、材料屈服強(qiáng)度為250 MPa。

在ANSYS中設(shè)定齒面接觸條件為無摩擦接觸，接觸力算法選擇增廣拉格朗日法，接觸力Fc=knx+λ，其中kn為接觸剛度，x為接觸位移響應(yīng)，λ為一外加的接觸力變量，理想的接觸是沒有變形的，所以kn越大結(jié)果越精確，但收斂越慢。設(shè)定一齒輪軸孔為固定支撐，另一齒輪的軸孔為無摩擦支撐加扭矩100 N·m。

將幾何模型劃分成六面形狀單元體，并在接觸點附近劃分加密網(wǎng)格如圖3所示，加密區(qū)域網(wǎng)格尺寸設(shè)為0.2 mm。

圖2 齒根處帶裂紋的 圖3 齒輪模型劃分齒輪模型

2.1 嚙合剛度等效彈簧的位移響應(yīng)

給定齒輪一個初始轉(zhuǎn)速，初始角位移為零，設(shè)定相對誤差為10-4，時間步階為10-4s，計算系統(tǒng)0～5 s的振動信號，截取系統(tǒng)穩(wěn)態(tài)的振動信號。

取ζ =0.06，在不同轉(zhuǎn)速下得到θ1和θ2的數(shù)值解，計算嚙合剛度等效彈簧的位移響應(yīng)x。齒輪旋轉(zhuǎn)一周，則每對齒均嚙合一次，則剛度變化呈周期性，所以在齒輪的振動成分中，存在嚙合頻率的成分。嚙合頻率的公式為fmesh=zω/60，圖4為不同轉(zhuǎn)速下位移響應(yīng)x的時域圖，可以發(fā)現(xiàn)嚙合頻率為45 Hz和90 Hz時振幅較大，嚙合頻率90 Hz時振幅最大。

圖4 無缺陷齒輪位移響應(yīng)時域圖

同樣選取ζ =0.06，得到齒根帶裂紋齒輪不同轉(zhuǎn)速下嚙合剛度等效彈簧的位移響應(yīng)x的時域圖，如圖5所示，嚙合頻率為45 Hz和90 Hz時振幅較大，嚙合頻率90 Hz時振幅最大，與無缺陷齒輪相似。由于只有一個輪齒根部有裂紋，所以振動信號每經(jīng)過23周期就出現(xiàn)一次較大的振幅。由此可見，無論齒輪有無缺陷，齒輪嚙合頻率接近嚙合剛度等效彈簧固有頻率或者固有頻率的1/2時，等效彈簧的位移響應(yīng)會很大，隨著齒輪轉(zhuǎn)速增大等效彈簧的位移響應(yīng)迅速減小。因此齒輪的轉(zhuǎn)速應(yīng)該避免使嚙合頻率接近等效彈簧的固有頻率或其半頻，才能避免系統(tǒng)產(chǎn)生很大的振幅，齒輪工作時轉(zhuǎn)速應(yīng)遠(yuǎn)離這個區(qū)域。

科胡特曾說：“個體如果想要發(fā)展出健康、統(tǒng)合、連續(xù)的個體，關(guān)系背景是十分重要的；只有在這個關(guān)系中嬰兒的自體客體需要通過養(yǎng)育者的情感協(xié)調(diào)得以滿足，個體才有可能發(fā)展出健康的自我。也正是因為如此，如果養(yǎng)育者給予嬰兒的反應(yīng)并不協(xié)調(diào)，個體就有可能出現(xiàn)心理問題。”

圖5 有缺陷齒輪位移響應(yīng)時域圖

2.2 振動特性對比分析

圖6 等效彈簧位移的均方根曲線

獲取振動達(dá)到穩(wěn)態(tài)時的信號，截取3 s～5 s的時域信號，各轉(zhuǎn)速下等效彈簧位移的均方根曲線如圖6所示，可見轉(zhuǎn)速為114.8 rpm(45 Hz)和227 rpm(90 Hz)對應(yīng)的均方根處出現(xiàn)峰值，兩種齒輪出現(xiàn)峰值的位置相同，但有裂紋齒輪組的均方根在很大范圍內(nèi)略大于無缺陷齒輪組，可見齒根裂紋會使齒輪嚙合出現(xiàn)較大的振動。

2.3 有齒根缺陷齒輪振動頻域分析

圖7為帶裂紋齒輪組在不同轉(zhuǎn)速下的嚙合頻率頻譜圖，可見每個較大振幅旁均有邊頻出現(xiàn)，且齒輪轉(zhuǎn)速越高邊頻越明顯，邊頻的出現(xiàn)是由振幅調(diào)制現(xiàn)象引起的，振幅調(diào)制由數(shù)學(xué)原理解釋為兩種信號在時域上相乘相當(dāng)于信號在頻域上折積，若數(shù)值模擬的結(jié)果在頻域嚙合頻率兩側(cè)出現(xiàn)邊頻，說明嚙合頻率為載波成分，齒輪轉(zhuǎn)速頻率為調(diào)制波成分。

隨著齒輪轉(zhuǎn)速增加，最大振幅先增大后迅速減小，齒輪轉(zhuǎn)速227 rpm對應(yīng)的最大振幅達(dá)到了14.69 μm，為所有情況下振幅最大值。

圖7 有裂紋齒輪組的嚙合頻率頻譜圖

3 結(jié)論

1)無論齒輪有無齒根裂紋，齒輪嚙合頻率接近嚙合剛度等效彈簧固有頻率或者固有頻率的1/2時，等效彈簧的位移響應(yīng)均會增大，但隨著齒輪轉(zhuǎn)速提高，等效彈簧的位移響應(yīng)迅速減小，可知使齒輪的嚙合頻率遠(yuǎn)離固有頻率可避免共振產(chǎn)生。

3) 帶齒根裂紋的齒輪組在最大振幅對應(yīng)的頻率附近出現(xiàn)邊頻現(xiàn)象，且齒輪轉(zhuǎn)速越高邊頻越明顯，此為帶齒根裂紋的齒輪組的故障特征，這些都是無缺陷齒輪不具備的特征， 可依此來判斷齒輪的齒根

處是否出現(xiàn)裂紋。

[1] 康傳章，陳兵奎，張 靖，等. 齒輪系統(tǒng)振動信號的邊頻帶研究及其應(yīng)用[J]. 機(jī)械傳動，2013，37(4): 41-43+76.

[2] Gosselin C M, Pereault L, Vailancourt C. Simulation and computeraided kinematics design of three-degreed-of freedom spherical parallel manipulators[J]. Journal of Robotic Systems, 1995, 12(12) : 857- 869.

[3] 劉辛軍, 張立杰, 高 峰. 基于AutoCAD 平臺的六自由度并聯(lián)機(jī)器人位置工作空間的機(jī)械求解方法[J].機(jī)器人，2000，22(6)：457-464.

[4] 陳修龍, 趙永生. 并聯(lián)機(jī)床運(yùn)動仿真與運(yùn)動學(xué)分析的新方法[J]. 中國機(jī)械工程, 2006, 17(5) : 467- 471.

[5] Kahraman A， Singh R．Non-linear dynamics of a geared rotor-bearing system with multiple clearances[J]．Journal of Sound and Vibration, 1991, 144(3)：469-506.

[6] Howard I, Jia S. The dynamic modeling of a spur gear in mesh including friction and crack[J]. Mechanical Systems and Signal Processing, 2001, 15(5):831-853.

Fault diagnosis method of gear root with rifts based on ANSYS

WANG Long, TIAN Huaigu, LI Sen, LIU Shoufa, DONG Feng

The dynamic model of the spur gear pair is built by lumped parameter method and the motion equations are solved by using Runge-Kutta method to acquire vibration response of the gear pair. The finite element model of defect-free gear pair and gear pair with root crack are established by ANSYS to analyze the vibration characteristics of the two gear pairs under conditions of different rotating speed. Results show that effects of the root crack on the vibration frequency of the meshing stiffness equivalent spring is not obvious. The RMS of the amplitude of gear pair with root crack is slightly higher than the defect-free gear pair as the speed increasing. There is side frequency in the vicinity of the maximum amplitude during meshing transmission of the gear pair with root crack, the higher the speed of the gear, the more obvious the side frequency, which can determine whether the gear pair has root crack.

lumped parameter method, finite element analysis, dynamic response, equivalent spring

TH132.413

A

1002-6886(2016)06-0048-04

陜西教育廳科學(xué)研究計劃項目 (12JK1075)；西京學(xué)院教改項目(JGYB1616)。

王龍(1995- )，男，漢族，河南南陽人，校企合作試點培養(yǎng)學(xué)生。 劉守法(1980- )，男，漢族，河北滄州人，講師，工學(xué)碩士，研究方向: 航空用有色合金FSW和 FSP技術(shù)及機(jī)床加工精度預(yù)測研究。

2016-06-20