間隙對齒輪-轉(zhuǎn)子-軸承系統(tǒng)彎扭耦合振動的影響分析

向 玲, 賈 軼, 高雪媛, 邸薇薇, 李媛媛

(華北電力大學(xué) 機(jī)械工程系，保定 071003)

?

間隙對齒輪-轉(zhuǎn)子-軸承系統(tǒng)彎扭耦合振動的影響分析

向 玲, 賈 軼, 高雪媛, 邸薇薇, 李媛媛

(華北電力大學(xué) 機(jī)械工程系，保定 071003)

綜合了動態(tài)側(cè)隙、齒面摩擦、齒輪偏心及時(shí)變嚙合剛度等因素，建立了齒輪-轉(zhuǎn)子-滾動軸承系統(tǒng)的彎扭耦合非線性動力學(xué)模型，通過分析動態(tài)側(cè)隙及軸承間隙對系統(tǒng)的影響來探究軸承端與齒輪端振動之間的耦合作用關(guān)系。結(jié)果表明：相對于無間隙系統(tǒng)，動態(tài)側(cè)隙下軸承端的徑向振動在高速區(qū)較為強(qiáng)烈，而齒輪的扭轉(zhuǎn)振動在整個(gè)轉(zhuǎn)速區(qū)幅度較大，隨轉(zhuǎn)速變化時(shí)系統(tǒng)提前通過非線性跳躍進(jìn)入主共振區(qū)。動態(tài)側(cè)隙的改變對軸承端的振動影響不大，但對扭振作用明顯；軸承間隙的大小對系統(tǒng)徑向和扭轉(zhuǎn)振動有著顯著影響，隨著軸承間隙的變化，兩者的時(shí)域特征及頻譜存在著規(guī)律性的變化；另外，軸承間隙直接影響著動態(tài)側(cè)隙的大小。分析結(jié)果對含間隙齒輪轉(zhuǎn)子系統(tǒng)的研究具有重要的理論與工程價(jià)值。

動態(tài)側(cè)隙；彎扭耦合；軸承間隙；偏心；嚙合剛度

齒輪系統(tǒng)是機(jī)械系統(tǒng)的重要組成部分，長期以來，人們對齒輪傳動系統(tǒng)的振動特性進(jìn)行了各種各樣的理論和實(shí)驗(yàn)研究，也不斷揭示了齒輪系統(tǒng)的運(yùn)行機(jī)理及存在的一些復(fù)雜的非線性現(xiàn)象。其中KAHRAMAN[1-2]通過建立多樣化的動力學(xué)模型對含間隙齒輪系統(tǒng)的非線性動力學(xué)特性進(jìn)行了一系列深入的研究，為后續(xù)學(xué)者進(jìn)行相關(guān)研究奠定了理論基礎(chǔ)。MORADI[3]等利用多尺度法研究分析了設(shè)計(jì)參數(shù)對齒輪非線性系統(tǒng)主共振、超諧和次諧共振的影響。LI等[4]研究了輪齒表面粗糙度，潤滑溫度及其他工況參數(shù)對齒輪系統(tǒng)振動特性的影響，結(jié)果表明這些因素對系統(tǒng)的彎扭耦合振動及動態(tài)嚙合起著重要作用。FARSHIDIANFAR等[5]建立了含齒側(cè)間隙和徑向間隙的動力學(xué)模型，著重分析研究了齒輪系統(tǒng)的混沌和分岔現(xiàn)象。高洪波等[6]綜合考慮了動態(tài)側(cè)隙、偏心及摩擦等因素著重研究了齒輪全齒磨損和偏心磨損下系統(tǒng)的振動特性。李朝峰等[7]建立了含齒輪偏心及軸承非線性的齒輪-轉(zhuǎn)子-軸承系統(tǒng)彎扭耦合動力學(xué)模型，并分析了各設(shè)計(jì)參量對系統(tǒng)振動響應(yīng)的影響，但該模型中并沒有考慮嚙合剛度的時(shí)變性、齒側(cè)間隙以及齒面摩擦的作用。

本文基于上述研究，綜合考慮了動態(tài)側(cè)隙、齒輪偏心、軸及軸承等非線性因素，建立了齒輪-轉(zhuǎn)子-軸承傳動系統(tǒng)的非線性動力學(xué)模型，并對其進(jìn)行了數(shù)值仿真分析。著重研究了齒側(cè)間隙、軸承間隙對系統(tǒng)彎扭耦合振動特性的影響及兩者的耦合作用關(guān)系。

1 動力學(xué)模型

圖1所示為動態(tài)條件下齒輪-轉(zhuǎn)子-軸承傳動系統(tǒng)的集中質(zhì)量模型，各軸承部分幾何坐標(biāo)原點(diǎn)為A、B、C、D。主/從動齒輪的理想中心O和O’均位于軸承孔的理想中心線上，按圖示主坐標(biāo)系(x,y,z)的方向可得到各軸承及齒輪的固定坐標(biāo)系。Oi,Gi(i= 1,2)分別為主從動齒輪的形心和質(zhì)心。兩軸與各相應(yīng)軸承的集中質(zhì)量部分對應(yīng)的當(dāng)量質(zhì)量為mbi(i= 1,2,3,4)。主/從動齒輪的質(zhì)量為mi(i= 1,2)，對應(yīng)質(zhì)心的轉(zhuǎn)動慣量分別為Ji(i= 1,2)。系統(tǒng)輸入端電動機(jī)轉(zhuǎn)子及輸出端負(fù)載轉(zhuǎn)子對應(yīng)的轉(zhuǎn)動慣量為Jm、Jp。

假設(shè)不考慮齒輪的軸向竄動，在xoy平面上，現(xiàn)設(shè)主、從動齒輪的形心坐標(biāo)為O1(x1,y1),O2(x2,y2)，對應(yīng)的質(zhì)心坐標(biāo)為G1(xo1,yo1)),G2(xo2,yo2)，偏心距O1G1=ρ1,O2G2=ρ2。主/從動輪的轉(zhuǎn)角為φ1,φ2，對應(yīng)的振動角位移分別為θ1,θ2。設(shè)主從動輪的初始旋轉(zhuǎn)角均為0。以主動輪轉(zhuǎn)動方向?yàn)檎较?，則由幾何關(guān)系可得

(1)

圖1 齒輪傳動系統(tǒng)的動力學(xué)模型Fig.1 Dynamic model of gear transmission system

另設(shè)軸承旋轉(zhuǎn)中心的位移坐標(biāo) (xbi,ybi) (i= 1,2,3,4)，根據(jù)軸的變形協(xié)調(diào)條件可知，主被動軸在兩齒輪旋轉(zhuǎn)中心位置x,y方向的彈性變形大小可表示為

(2)

式中:ξ1=lbi/lj(i= 1,2,j= 2;i= 3,4,j= 2)，如圖1中所示，l1,l2為主從動軸在兩軸承間的長度；lbi(i= 1,2,3,4)為主從動軸上齒輪旋轉(zhuǎn)中心到相應(yīng)軸承質(zhì)心的距離。由系統(tǒng)的幾何關(guān)系可得，在某個(gè)嚙合時(shí)刻輪齒的嚙合力如式(3)；f(…)為間隙性非線性函數(shù)。

DTE=rb1φ1-rb2φ2+y2-y1+ρ2sinφ2-ρ1sinφ1

(3)

1.1 摩擦力及摩擦力矩

齒面摩擦力可由齒輪嚙合的動態(tài)嚙合力得到，如式(4)，其中，λi,μi分別為摩擦力方向系數(shù)和時(shí)變摩擦系數(shù)，對于λi，嚙合點(diǎn)在節(jié)點(diǎn)以上取“+1”，反之則取“-1”。

(4)

對于μi，其大小隨嚙合位置作周期變化，為了接近現(xiàn)實(shí)，本文采用經(jīng)典的Buckingham半經(jīng)驗(yàn)公式[8]：

(5)

其中Vsi為齒廓接觸點(diǎn)的相對滑動速度，由文獻(xiàn)[8]可知

Vsi=w1×S1i-w2×S2i

(6)

式中:S1i,S2i為齒對i之間的摩擦力對齒輪1、齒輪2的摩擦力臂，由圖1(c)可知，w1,w2分別為兩齒輪的角速度；對于S1i/S2i，當(dāng)i=1時(shí)為K1N1/K1N2，當(dāng)i=2時(shí)為K2N1/K2N2；當(dāng)0<ε<1, 力臂可表示為

(7)

式中:ra1,ra2為主從動輪的齒頂高半徑，t1為一個(gè)雙齒嚙合時(shí)間，Pb為基圓節(jié)距，α為壓力角。

1.2 動態(tài)側(cè)隙

齒輪安裝時(shí)要留有一定的齒側(cè)間隙，其隨著齒輪中心距或齒厚的變化而變化，當(dāng)齒輪系統(tǒng)由于軸承振動或偏心引起中心距變化時(shí)，齒側(cè)間隙可表示為:

b=2b0+2a0cosα0(tanα′-tanα0+α0-α′)

α′=arccos(a0cosα0/a′)

a′=

(8)

式中:a0,α0為初始中心距和壓力角，a’,α’為實(shí)際嚙合過程中的中心距和壓力角，b0為初始單邊齒側(cè)間隙，這樣動態(tài)嚙合力式(3)中的非線性函數(shù)表示為

(9)

1.3 時(shí)變嚙合剛度

在直齒輪傳動時(shí)，輪齒的變形及剛度均作周期變化且存在突變；有不少研究將嚙合剛度近似為矩形波或傅里葉級數(shù)的形式，但與實(shí)際結(jié)果存在一定差異，本文采用文獻(xiàn)[9]的能量法并結(jié)合文獻(xiàn)[10]的修正法進(jìn)行計(jì)算。這里僅給出一對齒綜合嚙合剛度的計(jì)算公式：

k=

(10)

各參數(shù)的定義及計(jì)算公式詳見文獻(xiàn)[9-10]。當(dāng)兩對齒同時(shí)參與嚙合時(shí)，另一對齒的嚙合剛度計(jì)算同上式，將兩對齒嚙合剛度相加即為齒輪副的綜合嚙合剛度。為了后續(xù)分析，這里給出經(jīng)計(jì)算得到的一個(gè)嚙合周期內(nèi)的時(shí)變剛度曲線(圖2)，齒輪參數(shù)如表1所示。

圖2 時(shí)變嚙合剛度Fig.2 Time-varying mesh stiffness

1.4 非線性軸承力模型

圖3為滾動軸承模型示意圖，假定軸承的內(nèi)圈與旋轉(zhuǎn)軸、外圈與軸承座均為剛性連接，滾動體等間距排列且與滾道之間為純滾動。軸承變形主要是滾珠和滾道的接觸變形。

圖3 滾動軸承模型Fig.3 Rolling bearing model

圖3中第j個(gè)滾珠與滾道的法向接觸變形可表示為

δj=xbcosφj+ybsinφj-r0

(11)

式中:xb,yb分別為內(nèi)圈中心在x,y方向的振動位移，r0為滾動軸承間隙，φj為第j個(gè)滾珠的位置角。軸承的非線性力可由軸承運(yùn)動學(xué)和赫茲基理論推得[11]。

(12)

式中:nb為滾珠個(gè)數(shù)；Kc為赫茲接觸剛度系數(shù)；H(x)為Heaviside函數(shù)，當(dāng)x>0，該函數(shù)返回值為1，反之，返回值為0，其反映著接觸變形與接觸力的關(guān)系。

根據(jù)以上分析，進(jìn)一步引入系統(tǒng)的廣義坐標(biāo)列陣[φmx1y1φ1x2y2φ2xb1yb1xb2yb2xb3yb3xb4yb4φb]T，求得系統(tǒng)的動能，勢能和耗散函數(shù)，根據(jù)拉格朗日方程可得系統(tǒng)的非線性動力學(xué)方程，如式(13)。

式中xi,yi為主/從動輪沿x,y方向的振動位移；xbi,ybi則對應(yīng)主/從動軸軸承在x,y方向的振動位移；ksi,kti(i=1,2)分別為輸入軸和輸出軸的彎曲剛度和和扭轉(zhuǎn)剛度；csi,cti(i=1,2)則對應(yīng)軸的彎曲阻尼和扭轉(zhuǎn)阻尼；cbxi,cbyi(i=1,2,3,4)為軸承的徑向阻尼；Fbxi,Fbyi為軸承的非線性力；M1,M2為輸入端和輸出端扭矩。

2 間隙對系統(tǒng)彎扭耦合振動的影響

間隙對系統(tǒng)振動特性及嚙合狀態(tài)有著重要影響，有必要對系統(tǒng)的動力學(xué)響應(yīng)及各間隙的作用進(jìn)行詳細(xì)分析。因?yàn)?，?13)綜合考慮了時(shí)變剛度和動態(tài)側(cè)隙等強(qiáng)非線性因素，方程參數(shù)量級差別較大，求解時(shí)易出現(xiàn)病態(tài)和不收斂問題，因此，采用求解穩(wěn)定性及精度較高的ode15s對其進(jìn)行求解，首先將方程寫成狀態(tài)方程的形式(降階)，然后設(shè)置求解精度(相對和絕對精度均為1.0×10-9)進(jìn)行計(jì)算。系統(tǒng)參數(shù)如表1所示，軸承均選擇SKF6200，查閱軸承手冊可得軸承的詳細(xì)參數(shù)。

2.1 齒側(cè)間隙的影響分析

齒側(cè)間隙直接影響著輪齒的嚙合狀態(tài)，為了分析齒側(cè)間隙大小對系統(tǒng)動力學(xué)響應(yīng)的影響規(guī)律，分別取初始齒側(cè)間隙b0為0 μm，0.05 μm，0.1 μm，并以無間隙作為對比對系統(tǒng)進(jìn)行數(shù)值仿真，此時(shí)主從動軸軸承的間隙分別為1 μm，2 μm。

首先在整個(gè)轉(zhuǎn)速內(nèi)分析系統(tǒng)的動力學(xué)響應(yīng)，圖4是主從動軸上軸承(編號1和3)沿x,y方向振動位移、動態(tài)傳動誤差(DTE)隨主動軸轉(zhuǎn)頻(f)變化的RMS(均方根值)曲線。從圖中看出，無間隙時(shí)，系統(tǒng)在f=100(嚙頻fm=1 900) 和f=192(fm=3 648)發(fā)生扭轉(zhuǎn)方向的共振，軸承端在f=100和f=202(fm=3 838)發(fā)生x方向的共振在y方向未發(fā)現(xiàn)共振峰的存在，注意到扭振和軸承端徑向振動主共振頻率略有差別。而對于動態(tài)間隙系統(tǒng)，同樣在f=100和f=166處出現(xiàn)了共振峰，且彎振和扭振的峰值頻率基本接近；當(dāng)增大b0時(shí)，峰值及固頻變化微弱。另外，扭振的RMS隨b0增大明顯增大(圖4(a))，同時(shí)動態(tài)側(cè)隙也隨之增大(如圖4(f))，但軸承端的徑向振動變化較小，說明了初始齒側(cè)間隙對于軸承端振動影響較弱。其次，由圖可知，對于扭振，動態(tài)側(cè)隙下的RMS比無間隙時(shí)大很多，這是因?yàn)閯討B(tài)側(cè)隙綜合考慮了橫向和徑向振動，其促使了扭振的加強(qiáng)，由式(8)也可判斷；而對于軸承端振動，在低速區(qū)，動態(tài)側(cè)隙下的均方根值和無間隙時(shí)相差不大，而在高速區(qū)，前者明顯比后者大，如圖4(b)～圖4(e)所示?？傊?，相對于無間隙模型，動態(tài)側(cè)隙在整個(gè)轉(zhuǎn)速區(qū)內(nèi)對系統(tǒng)的扭振影響顯著，而在高速區(qū)則對軸承端徑向振動影響較大。另外，系統(tǒng)在進(jìn)入主共振區(qū)(第二個(gè)峰值區(qū)域)時(shí)產(chǎn)生了跳躍，這種跳躍并沒有隨b0的增大而變化，但對于無間隙系統(tǒng)，系統(tǒng)在進(jìn)入此區(qū)域時(shí)是連續(xù)的，這側(cè)面反映了間隙的非線性作用。

對于輪齒的嚙合狀態(tài)，對比圖5曲線可知，在主共振區(qū)域，無間隙系統(tǒng)嚙合力的各項(xiàng)指標(biāo)明顯偏大，但共振區(qū)出現(xiàn)滯后；而在其他轉(zhuǎn)速區(qū)域，嚙合力各指標(biāo)相差較小。并且，由圖5(c)可知，在共振區(qū)，有/無間隙系統(tǒng)均處于雙邊沖擊狀態(tài)，而動態(tài)間隙的嚙合力相對較小，其存在減弱了共振區(qū)的輪齒沖擊；在其他區(qū)域，輪齒處于無沖擊狀態(tài)，動態(tài)側(cè)隙的改變并未影響該狀態(tài)。

(13)

圖4 各參量的RMS曲線Fig.4 RMS curve of all vibration variables

圖5 嚙合力的均方根值、最大值及最小值曲線Fig.5 RMS、Maximum and Minimum curve of dynamic mesh force

進(jìn)一步分析動態(tài)側(cè)隙對系統(tǒng)振動特性的影響，設(shè)置f=30 Hz，b0=0.1 μm，這里從頻譜的角度進(jìn)行分析。首先確定該工況下各振動量的頻率成分，圖6為軸承端徑向振動的時(shí)間歷程圖及頻譜圖，以下主/被動軸軸承端分別用B1和B2表示。對于橫向(x)振動：B1處以2f1(f1=30 Hz，軸頻)頻率為主，并出現(xiàn)了高倍軸頻、嚙頻(fm=N1*f1=570 Hz)及其二倍頻成分，如圖6(a2)；與之類似，B2處以2f2(f1*N1/N2=11.875 Hz)頻率為主，并出現(xiàn)了fm及其二倍頻成分，如圖6(c2)。對于縱向(y)振動：B1處以2f1頻為主，且存在fm及其倍頻成分，如圖(b2)；B2處以2f2頻為主，且存在fm及其倍頻成分，如圖6(d2)。另外注意到，B1/B2徑向(x,y)振動在低頻處的頻帶較寬且連續(xù)，同時(shí)存在其他頻率成分，這應(yīng)該與軸承的非線性有關(guān)。對于扭轉(zhuǎn)振動：時(shí)域圖(圖7)中可明顯看到?jīng)_擊的存在，且DTE及動態(tài)嚙合力(DMF) 的頻帶較寬。

DTE頻譜以fm及其倍頻為主，且存在2f1、2f2及其他隨機(jī)頻率成分(表現(xiàn)為低頻)。DMF頻譜以fm及其倍頻為主，且出現(xiàn)了f1、f2及以其倍頻為頻率間隔的頻率成分，另外也存在一些隨機(jī)頻率成分，總之DMF頻帶的頻率成分較復(fù)雜。

圖6 軸承端徑向振動的時(shí)域波形及頻譜Fig.6 Transverse vibration waveform and spectrum diagram of bearing end

圖7 DTE和動態(tài)嚙合力的時(shí)域波形及頻譜Fig.7 Waveform and spectrum diagram of DTE and DMF

為了全面了解動態(tài)側(cè)隙的作用，圖8給出了各振動量隨初始齒側(cè)間隙b0變化的三維頻譜圖。對比各圖可知，當(dāng)b0變化時(shí)，無論軸承端的徑向振動或是齒輪的扭振，各頻譜中主要頻譜成分沒有明顯的增減。結(jié)合圖4可知，b0的變化對于軸承端的振動影響并不大。另外，b0對于扭振 (DTE)幅度有顯著的影響，但對DMF影響較小，同時(shí)結(jié)合式(3)及式(9)可知，初始側(cè)隙對動態(tài)側(cè)隙系統(tǒng)扭振的影響類似于一個(gè)線性增加的常量，其并未對DMF產(chǎn)生較大影響，同時(shí)該常量在頻譜中相當(dāng)于直流分量，因此各頻率成分并未隨b0的變化明顯改變。當(dāng)然，上述結(jié)果僅限于初始側(cè)隙在允許范圍內(nèi)變化時(shí)得到。

圖8 各振動參量的三維頻譜圖Fig.8 Three-dimensional spectrum of all vibration variables

2.2 軸承間隙的影響分析

首先主/從動軸軸承用C1/C2表示，主/從動齒輪用G1/G2表示。為了分析軸承間隙的影響規(guī)律，分別取C1間隙r1為0.5 μm，1 μm，2 μm及4 μm，為了保持C1/C2間隙的定量關(guān)系，C2間隙設(shè)為2r1，b0設(shè)為0.1 μm。圖9為G1/G2徑向振動位移、DTE及動態(tài)側(cè)隙隨主動軸軸頻f變化的RMS曲線。

由圖可見，對于G1/G2的橫向(x)振動(如圖9(b)、(d) 所示)：系統(tǒng)在f=100 Hz(fm=1 900 Hz) 和f=166 Hz(fm=3 154 Hz)處出現(xiàn)共振峰，隨著C1/C2間隙的增加，在中高轉(zhuǎn)速區(qū)，曲線均小幅上移，對應(yīng)振動的加強(qiáng)；而在低速區(qū)(多為工作轉(zhuǎn)速區(qū))，G1/G2的各振動曲線相差較小；僅在r1=4 μm時(shí)，主動輪振動明顯較弱，而被動輪振動較強(qiáng)，這是由于間隙相對過大而引起軸承非線性力劇烈變化而造成的。對于G1/G2的縱向(y)振動(圖9(c)、(e))：振動在該方向上隨轉(zhuǎn)速變化比較平穩(wěn)且未出現(xiàn)共振峰；隨著C1/C2間隙的增加，曲線整體上移，振動明顯加強(qiáng)，說明了間隙的增加使得該向的軸承力或軸的彎曲力明顯增加。對于扭振 (圖9(a),(f))：當(dāng)軸承間隙增加時(shí)，DTE及動態(tài)側(cè)隙的RMS曲線均整體下移，扭振幅度逐步減弱。

關(guān)于軸承間隙對嚙合狀態(tài)的影響，圖10給出了DMF均方根值，最大值及最小值隨軸頻f的變化曲線。圖中可看到共振峰及跳躍現(xiàn)象的存在；隨著C1/C2間隙的增加，各曲線并未明顯變化，嚙合沖擊狀態(tài)變化很小，說明軸承間隙雖引起徑向振動的改變，但這種改變所引起力或力矩的變化并不足以使輪齒嚙合狀態(tài)發(fā)生大的改變。

同樣，從頻譜的角度進(jìn)一步分析軸承間隙的影響，這里設(shè)f=30 Hz，b0=0.1 μm；圖11為各振動量隨C1間隙r1(r2=2r1)變化的三維頻譜圖。由于G1/G2徑向振動頻譜成分基本與軸承端(圖6,7)相同，在此不再贅述。

圖9 各振動參量的RMS曲線Fig.9 RMS curve of all vibration variables

圖10 嚙合力曲線Fig.10 RMS、Maximum and Minimum curve of dynamic mesh force

對于G1/G2的橫向振動(圖11(c), (e))，各倍頻成分變化不大，以主頻和嚙頻變化為主。主動輪：主頻2f1幅值隨r1呈現(xiàn)先增后減再增的趨勢，在r1=3.5 μm取最小，如圖12(a)；嚙頻fm幅值變化如圖12(c)，呈現(xiàn)先減后增的趨勢，在r1=1.5 μm取最小，但幅值變化不大；被動輪：主頻2f2幅值變化同主動輪主頻2f1，在r1=1.5 μm取最小(圖12(b))，fm幅值在r1=0.4 μm取最小(圖12(d))。對比可知G2主頻平均幅值比G1大，而G1嚙頻平均幅值較G2大，G2振動相對較強(qiáng)。

對于G1/G2的縱向振動(圖11(d), (f))，仍以主頻和嚙頻為主。主動輪：主頻2f1幅值隨r1呈現(xiàn)先增后減的趨勢，在r1=2.8 μm取最大，如圖13(a)所示；嚙頻fm幅值變化規(guī)律性不強(qiáng)，但波動較小(圖13(c))。被動輪：主頻2f2幅值先增后減再增，在r1=2.8 μm最小(圖13(b))；而嚙頻fm幅值則隨r1不斷跳變，但變化幅度較小(圖13(d))。

圖11 各振動量的三維FFT譜圖Fig.11 Three-dimensional spectrum of all vibration variables

對于G1/G2的扭振(圖11(a), (b))，振動頻帶較寬，以嚙頻及其倍頻為主。當(dāng)r1變化時(shí)，扭振頻帶以G1/G2徑向振動的主頻(2f1,2f2)變化為主。主動輪：2f1幅值呈現(xiàn)先增后減再增的趨勢，同圖12(a)中x1向振動，在r1=3.2 μm取最大，如圖14(a)；2f2幅值也呈現(xiàn)先增后減再增的趨勢，同圖12(b) 中x2振動，在r1=1.5 μm取最小 (圖14(b))；嚙頻fm幅值總體上呈遞增趨勢，幅值較大(圖14(c))。對于嚙合力，嚙頻fm幅值則隨r1不斷跳變，無明顯規(guī)律，變化幅度較小(圖14(d))。

圖12 主/被動輪沿x方向振動的主頻幅值隨軸承間隙的變化Fig.12 The main frequency amplitude curve via bearing clearance along the X direction

圖13 主/被動輪沿y方向振動的主頻幅值隨軸承間隙的變化Fig.13 The main frequency amplitude curve via bearing clearance along the Y direction

圖14 扭振的主頻幅值隨軸承間隙的變化Fig.14 The main frequency amplitude curve of torsional vibration via bearing clearance

綜上所述，合理范圍內(nèi)增大軸承間隙，齒輪縱向(y向)振動明顯增強(qiáng)；橫向(x向)振動在中高速區(qū)內(nèi)增強(qiáng)，而在低速區(qū)內(nèi)變化較小，而當(dāng)間隙相對較大時(shí)，才會有較大變化；對于扭振，其幅度隨軸承間隙的增大則逐漸減小，動態(tài)側(cè)隙隨之也逐漸減小，但輪齒的嚙合狀態(tài)并未明顯變化；另外，軸承間隙的改變并未對系統(tǒng)扭振或彎振的固有頻率有所影響；而在頻域內(nèi)，齒輪徑向振動的主頻(2f1,2f2)幅值受軸承間隙影響較大，其隨間隙變化時(shí)出現(xiàn)了最值，嚙頻(fm)幅值同樣受間隙的影響，但其變化幅度較小。扭振頻譜中包含齒輪徑向振動的主頻成分(2f1,2f2)，軸承間隙對其影響主要表現(xiàn)在對這兩種頻率成分的影響，而對于譜中嚙頻及其倍頻成分的影響較弱。

3 結(jié) 論

本文建立了直齒輪傳動系統(tǒng)的彎扭耦合非線性動力學(xué)模型，主要研究了初始側(cè)隙對齒輪扭振和軸承端振動的影響、軸承間隙對齒輪扭振和徑向振動影響，以此來探究軸承端振動與齒輪端振動之間的作用關(guān)系。所得結(jié)論如下：

(1) 相對于無間隙系統(tǒng)，動態(tài)側(cè)隙下齒輪扭振幅度明顯較強(qiáng)；對于軸承端的徑向振動，在高速區(qū)，動態(tài)側(cè)隙系統(tǒng)的振動幅度較無間隙系統(tǒng)大，而在低速區(qū)(大部分系統(tǒng)的工作區(qū)域)，兩者相差較??；動態(tài)側(cè)隙的影響作用與齒輪轉(zhuǎn)速有著密切的關(guān)系。另外，動態(tài)側(cè)隙下系統(tǒng)的主共振較無間隙時(shí)有所提前，且在進(jìn)入主共振區(qū)域時(shí)出現(xiàn)常見的非線性跳躍現(xiàn)象。再者，隨著初始側(cè)隙的增加，系統(tǒng)的扭振幅度隨之增大，但軸承端的徑向振動變化較??；而初始側(cè)隙的增加對系統(tǒng)振動的各主要頻率成分影響并不大。

(2) 軸承間隙對齒輪的徑向振動和扭轉(zhuǎn)振動影響較大，隨著軸承間隙的增加，齒輪的徑向振動整體上均有所加強(qiáng)，但扭振的幅度有所降低；軸承間隙對齒輪的嚙合力、嚙合狀態(tài)及系統(tǒng)的彎、扭固有頻率影響很弱。而在頻域內(nèi)，隨著軸承間隙的增加，齒輪徑向振動的主頻成分均呈現(xiàn)規(guī)律性的變化，且在此過程中有最值的出現(xiàn)；扭振頻譜的變化表現(xiàn)在以徑向振動主頻為頻率成分的規(guī)律性變化，同樣有最值的存在。而在各頻譜中，嚙頻及其倍頻成分變化較小。這對軸承間隙的選擇具有一定的理論指導(dǎo)意義。

(3) 總體上，初始側(cè)隙的變化(動態(tài)側(cè)隙隨之變化)對軸承端的徑向振動影響不大，而軸承端的間隙變化則對齒輪端的徑向和扭轉(zhuǎn)振動的影響顯著，且軸承間隙直接影響著動態(tài)側(cè)隙的大小。

[1] KAHRAMAN A, SINGH R. Interactions between time-varying mesh stiffness and clearance non-linearities in a geared system[J]. Journal of Sound & Vibration, 1991, 146(1):135-156.

[2] KAHRAMAN A, SINGH R. Non-linear dynamics of a geared rotor-bearing system with multiple clearances[J]. Journal of Sound & Vibration, 1991, 144(3):469-506.

[3] MORADI H, SALARIEH H. Analysis of nonlinear oscillations in spur gear pairs with approximated modelling of backlash nonlinearity[J]. Mechanism and Machine Theory, 2012, 51: 14-31.

[4] LI S, KAHRAMAN A. A tribo-dynamic modelof a spur gear pair[J]. Journal of Sound and Vibration, 2013, 332(20): 4963-4978.

[5] FARSHIDIANFAR A, SAGHAFI A. Bifurcation and chaos prediction in nonlinear gear systems[J]. Shock and Vibration, 2014(2):1-8.

[6] 高洪波, 李允公, 劉杰. 基于動態(tài)側(cè)隙的齒輪系統(tǒng)齒面磨損故障動力學(xué)分析[J]. 振動與沖擊, 2014, 33(18): 221-226. GAO Hongbo, LI Yungong, LIU Jie. Dynamic analysis of a spur gear system with tooth-wear faults based on dynamic backlash[J]. Journal of Vibration and Shock,2014,33(18): 221-226.

[7] 李朝峰, 周世華, 劉文明, 等. 齒輪-轉(zhuǎn)子-滾動軸承傳動系統(tǒng)的彎扭耦合振動[J]. 航空動力學(xué)報(bào), 2014, 29(7): 1543-1555. LI Chaofeng, ZHOU Shihua, LIU Wenming, et al. Vibrationof bending-torsion coupling gear-rotor-rolling bearing transmission system[J]. Journal of Aerospace Power, 2014, 29(7): 1543-1555.

[8] BUCKINGHAM E. Analytical mechanics of gears[M]. Courier Corporation, 1988.

[9] TIAN X. Dynamic simulation for system response of gearbox including localized gear faults[M]. Library and Archives Canada= Bibliothèque et Archives Canada, 2005.

[10] CHAARI F, FAKHFAKH T, HADDAR M. Analytical modelling of spur gear tooth crack and influence on gearmesh stiffness[J]. European Journal of Mechanics-A/Solids, 2009, 28(3): 461-468.

[11] SAWALHI N, RANDALL R B. Simulating gear and bearing interactions in the presence of faults: Part I. The combined gear bearing dynamic model and the simulation of localised bearing faults[J]. Mechanical Systems and Signal Processing, 2008, 22(8): 1924-1951.

Effects of backlash on bending-torsion coupled vibration of a gear-rotor bearing system

XIANG Ling, JIA Yi, GAO Xueyuan, DI Weiwei, LI Yuanyuan

(Mechanical Engineering Department, North China Electric Power University, Baoding 071003, China)

A dynamic model for bending-torsion coupled vibration of a gear-rotor-bearing system was developed considering dynamic backlash, tooth face friction, gear eccentricity and time-varying mesh stiffness. The coupled vibration relationship between bearing end and gear end was investigated by analyzing effects of dynamic backlash and bearing clearance on the system. The results showed that compared with a system without backlash, the radial vibration of the bearing end of the system with dynamic backlash is stronger within a high rotating speed zone, while the amplitude of the gear torsional vibration is larger within the whole rotating speed range; the system with dynamic backlash enters the main resonance region in advance through nonlinear jump; the change of dynamic backlash has little effect on the vibration of bearing end, but it has an obvious effect on the gear torsional vibration; the magnitude of bearing clearance has a significant influence on the radial and torsional vibrations of the system; with the variation of bearing clearance, the time-domain characteristics and frequency spectra of two vibrations vary regularly; moreover, the bearing clearance directly affects the magnitude of dynamic backlash. The results had important theoretical and engineering values for studying gear-rotor systems with clearances.

dynamic backlash; bending-torsion coupling; bearing clearance; eccentricity; mesh stiffness

國家自然科學(xué)基金資助項(xiàng)目(51475164)

2015-09-07 修改稿收到日期：2015-11-07

向玲 女，博士,教授，1971年生

TH113

A

10.13465/j.cnki.jvs.2016.21.001