李同杰，靳廣虎，朱如鵬，安魯陵

?

滑動(dòng)軸承支撐下齒輪耦合轉(zhuǎn)子系統(tǒng)彎扭耦合振動(dòng)特性分析

李同杰1, 2，靳廣虎1，朱如鵬1，安魯陵1

(1. 南京航空航天大學(xué) 機(jī)電學(xué)院，江蘇 南京，210016； 2. 安徽科技學(xué)院 機(jī)械工程學(xué)院，安徽 滁州，233100)

在綜合考慮齒輪副齒側(cè)間隙以及滑動(dòng)軸承的油膜力等非線(xiàn)性因素的基礎(chǔ)上，建立滑動(dòng)軸承?雙轉(zhuǎn)子?齒輪耦合系統(tǒng)的非線(xiàn)性動(dòng)力學(xué)模型。通過(guò)數(shù)值仿真的方法研究耦合系統(tǒng)參數(shù)對(duì)齒輪副嚙合沖擊特性的影響規(guī)律，以及系統(tǒng)穩(wěn)定性隨轉(zhuǎn)速的分岔規(guī)律。研究結(jié)果表明，滑動(dòng)軸承的油膜對(duì)齒輪耦合轉(zhuǎn)子系統(tǒng)的混沌運(yùn)動(dòng)具有顯著的鎮(zhèn)定作用；滑動(dòng)軸承間隙以及轉(zhuǎn)子質(zhì)量偏心設(shè)計(jì)不當(dāng)將會(huì)導(dǎo)致系統(tǒng)齒輪副產(chǎn)生單邊沖擊現(xiàn)象。

齒輪耦合轉(zhuǎn)子；滑動(dòng)軸承；耦合非線(xiàn)性動(dòng)力學(xué)模型；嚙合沖擊；分岔

齒輪耦合轉(zhuǎn)子系統(tǒng)是各種傳動(dòng)系統(tǒng)的核心部分，其動(dòng)特性的優(yōu)劣直接影響著減速器的運(yùn)行品質(zhì)，所以齒輪轉(zhuǎn)子的研究一直都是一個(gè)學(xué)術(shù)熱點(diǎn)[1?8]。由于大型減速器基本上都是采用滑動(dòng)軸承支撐，所以把滑動(dòng)軸承集成到齒輪耦合轉(zhuǎn)子系統(tǒng)模型中，研究滑動(dòng)軸承?齒輪?轉(zhuǎn)子耦合系統(tǒng)的動(dòng)力學(xué)就具有重要的工程應(yīng)用價(jià)值。該耦合系統(tǒng)的早期研究多以線(xiàn)性模型為主，或者僅僅考慮某單一方面的非線(xiàn)性因素，例如在對(duì)滑動(dòng)軸承的處理上，眾多學(xué)者把軸承假設(shè)為支撐力隨轉(zhuǎn)子橫向位移線(xiàn)性變化的線(xiàn)彈性軸承[9?10]；在對(duì)齒輪嚙合力的處理上，忽略了齒側(cè)間隙可能導(dǎo)致的齒輪副的嚙合沖擊非線(xiàn)性[11?12]。同時(shí)綜合考慮了滑動(dòng)軸承油膜力非線(xiàn)性以及齒輪副嚙合沖擊非線(xiàn)性的研究較少。BAGUET等[13]利用有限元研究了單對(duì)齒輪?轉(zhuǎn)子?軸承耦合模型，考慮了非線(xiàn)性時(shí)變剛度和非線(xiàn)性支撐，對(duì)比了非線(xiàn)性支撐和線(xiàn)性支撐下系統(tǒng)動(dòng)態(tài)響應(yīng)之間的區(qū)別，并研究了軸承和齒輪載荷的動(dòng)態(tài)特性。THEODOSSIADES等[14]研究了油膜力和嚙合力共同作用下齒輪耦合轉(zhuǎn)子系統(tǒng)的動(dòng)態(tài)響應(yīng)，分析了油膜穩(wěn)定性，研究了不同參數(shù)組合下系統(tǒng)的幅頻響應(yīng)。GAO等[15]在綜合考慮了齒側(cè)間隙、傳動(dòng)誤差和時(shí)變嚙合剛度等非線(xiàn)性因素的基礎(chǔ)上，研究了滑動(dòng)軸承與滾動(dòng)軸承混合支撐下匯流傳動(dòng)齒輪轉(zhuǎn)子系統(tǒng)的運(yùn)動(dòng)分岔特性。然而，目前針對(duì)滑動(dòng)軸承支撐下齒輪耦合轉(zhuǎn)子系統(tǒng)振動(dòng)特性的研究還不完善，尤其是滑動(dòng)軸承非線(xiàn)性油膜力對(duì)齒輪副嚙合沖擊特性的影響規(guī)律方面的研究尚未見(jiàn)報(bào)道。本文作者將綜合考慮滑動(dòng)軸承的非線(xiàn)性油膜力以及齒輪副的嚙合沖擊力等非線(xiàn)性因素，建立滑動(dòng)軸承支撐下齒輪耦合雙軸轉(zhuǎn)子的集中參數(shù)形式的非線(xiàn)性動(dòng)力學(xué)模型，推導(dǎo)系統(tǒng)關(guān)于齒輪副嚙合點(diǎn)相對(duì)位移的運(yùn)動(dòng)微分方程，通過(guò)數(shù)值仿真的手段初步研究滑動(dòng)軸承非線(xiàn)性油膜力對(duì)齒輪副嚙合沖擊特性的影響，并對(duì)系統(tǒng)運(yùn)動(dòng)狀態(tài)的穩(wěn)定性隨轉(zhuǎn)速的分岔規(guī)律進(jìn)行研究。

1 動(dòng)力學(xué)模型

滑動(dòng)軸承支撐下齒輪耦合雙軸轉(zhuǎn)子系統(tǒng)的非線(xiàn)性動(dòng)力學(xué)模型如圖1所示。

圖1(a)中，主動(dòng)轉(zhuǎn)子(下標(biāo)1標(biāo)注)以及從動(dòng)轉(zhuǎn)子(下標(biāo)2標(biāo)注)均受滑動(dòng)軸承支撐，主、從動(dòng)齒輪均固定于轉(zhuǎn)軸跨中，2個(gè)轉(zhuǎn)子均為Jeffcott轉(zhuǎn)子，即轉(zhuǎn)子的2個(gè)軸頸運(yùn)動(dòng)規(guī)律完全對(duì)稱(chēng)?；瑒?dòng)軸承作用給轉(zhuǎn)子的非線(xiàn)性油膜力分解到固定坐標(biāo)x和y(=1, 2)上的投影分別為F，F；表示相應(yīng)齒輪的齒數(shù)，表示相應(yīng)轉(zhuǎn)子所承受的扭矩。每個(gè)轉(zhuǎn)子的平面運(yùn)動(dòng)均存在2個(gè)橫向振動(dòng)位移(,)和1個(gè)扭轉(zhuǎn)振動(dòng)位移，這樣整個(gè)系統(tǒng)的動(dòng)力學(xué)描述就需要6個(gè)廣義坐標(biāo)，即(1,1,1,2,2,2)。為齒輪副的嚙合力，齒輪副的時(shí)變嚙合剛度以及齒側(cè)間隙決定了分段非線(xiàn)性力的本質(zhì)；為轉(zhuǎn)子質(zhì)量偏心；為轉(zhuǎn)子的偏心變化相位；為轉(zhuǎn)子的轉(zhuǎn)動(dòng)角速度；g為轉(zhuǎn)子重力；b為齒輪基圓半徑；為齒輪壓力角。

運(yùn)用剛體平面運(yùn)動(dòng)微分方程可以建立系統(tǒng)廣義坐標(biāo)(1,1,1,2,2,2)下的運(yùn)動(dòng)微分方程，如式(1)所示。

圖1 滑動(dòng)軸承支撐下齒輪耦合雙軸轉(zhuǎn)子系統(tǒng)的非線(xiàn)性動(dòng)力學(xué)模型

不考慮動(dòng)力學(xué)模型中轉(zhuǎn)子所承受外扭矩的波動(dòng)，即扭轉(zhuǎn)振動(dòng)方程中的外扭矩1和2為定值；相較于嚙合力而言，重力對(duì)各自齒輪產(chǎn)生力矩的轉(zhuǎn)動(dòng)效應(yīng)非常微弱，所以扭轉(zhuǎn)振動(dòng)方程中忽略了重力對(duì)齒輪回轉(zhuǎn)中心的力矩。按照齒輪傳動(dòng)的物理關(guān)系，有傳動(dòng)比=2/1，21，1=101，2=202，其中，0為初始相位。

1.1 齒輪副嚙合力的非線(xiàn)性表達(dá)

齒輪副動(dòng)態(tài)嚙合力與嚙合點(diǎn)的相對(duì)位移有關(guān)，為了方便表達(dá)，引入齒輪副嚙合點(diǎn)的相對(duì)位移

式(2)關(guān)于時(shí)間2次求導(dǎo)后，結(jié)合方程(1)可得齒輪副嚙合點(diǎn)相對(duì)位移的振動(dòng)微分方程

這樣，齒輪副的動(dòng)態(tài)嚙合力可以表達(dá)為

式(4)中，近似按矩形波周期時(shí)變的嚙合剛度()可以展開(kāi)成為以齒輪副嚙頻(=11)為基頻的Fourier級(jí)數(shù)，取一次諧波項(xiàng)為

其中：m為嚙合剛度的均值，a為嚙合剛度波動(dòng)幅值，為剛度波動(dòng)初相位。

式(4)中的間隙非線(xiàn)性函數(shù)(,)充分考慮了齒輪副可能會(huì)出現(xiàn)的3種嚙合狀態(tài)，即正常嚙合、空嚙合以及齒背嚙合[16]

1.2 滑動(dòng)軸承油膜力的非線(xiàn)性表達(dá)

作用在轉(zhuǎn)子軸頸上的非線(xiàn)性油膜力F，F經(jīng)Sommerfeld 系數(shù)可實(shí)現(xiàn)量綱一化，即

國(guó)內(nèi)外專(zhuān)家學(xué)者的研究發(fā)現(xiàn)，滑動(dòng)軸承的量綱一油膜力與轉(zhuǎn)子徑向位移以及徑向速度成非線(xiàn)性關(guān) 系[17?18]，其表達(dá)如式(8)所示：

式中：

1.3 相對(duì)坐標(biāo)下的量綱一化方程

2 運(yùn)動(dòng)狀態(tài)穩(wěn)定性的分岔特性

仿真案例中的基本參數(shù)如下：齒輪模數(shù)=3 mm，齒輪1=40，2=80，齒輪壓力角=20°，齒側(cè)間隙=10 μm，主動(dòng)輪轉(zhuǎn)子質(zhì)量1=10.6 kg，從動(dòng)輪轉(zhuǎn)子質(zhì)量1=14.6 kg，主動(dòng)轉(zhuǎn)子半徑1=20 mm，從動(dòng)輪轉(zhuǎn)子直徑2=30 mm，根據(jù)減速器滑動(dòng)軸承參數(shù)的工程取值范圍[19]，軸承半徑相對(duì)間隙取=0.2‰~1.5‰，軸承寬徑比/=2，潤(rùn)滑油動(dòng)力黏度=0.018 N·s/m2。

忽略系統(tǒng)質(zhì)量偏心的影響(1=2=0)，方程(9)就退化為僅受齒輪副嚙合頻率激勵(lì)的非光滑非線(xiàn)性系統(tǒng)。取主動(dòng)轉(zhuǎn)子軸承半徑間隙為1=20 μm，從動(dòng)轉(zhuǎn)子軸承半徑間隙為1=30 μm，考查主動(dòng)轉(zhuǎn)子轉(zhuǎn)速1在500~900 rad/s之間變化時(shí)系統(tǒng)運(yùn)動(dòng)狀態(tài)穩(wěn)定性的分岔規(guī)律。

圖2表明：總體上，在500~810 rad/s之間的低速區(qū)間系統(tǒng)振動(dòng)呈現(xiàn)周期狀態(tài)，在810~900 rad/s之間的高速區(qū)間內(nèi)系統(tǒng)振動(dòng)呈現(xiàn)復(fù)雜運(yùn)動(dòng)狀態(tài)。

圖2 非線(xiàn)性油膜力作用下系統(tǒng)全局分岔圖

圖3 線(xiàn)彈性軸承支撐下系統(tǒng)的全局分岔圖

在500~810 rad/s的低速區(qū)間內(nèi)，隨著轉(zhuǎn)速的變化系統(tǒng)的周期運(yùn)動(dòng)的形態(tài)也不是一成不變的。在轉(zhuǎn)速區(qū)間500~610 rad/s內(nèi)，系統(tǒng)振動(dòng)表現(xiàn)為周期一狀態(tài)，其中當(dāng)1=550 rad/s時(shí)，周期一運(yùn)動(dòng)的時(shí)域圖、相圖以及Poincare?截面圖如圖4所示；在轉(zhuǎn)速區(qū)間610~ 690 rad/s內(nèi)，系統(tǒng)振動(dòng)以周期二狀態(tài)為主，期間出現(xiàn)了短暫的周期三運(yùn)動(dòng)狀態(tài)，當(dāng)1=620 rad/s時(shí)，周期二運(yùn)動(dòng)的時(shí)域圖、相圖以及Poincare?截面圖如圖5所示；在轉(zhuǎn)速區(qū)間690~740 rad/s內(nèi)，系統(tǒng)振動(dòng)又回到周期一狀態(tài)，當(dāng)轉(zhuǎn)速增大到740~810 rad/s區(qū)間，周期一運(yùn)動(dòng)又分岔為周期三運(yùn)動(dòng)，其中當(dāng)1=770 rad/s時(shí)，周期一運(yùn)動(dòng)的時(shí)域圖、相圖以及Poincare?截面圖如圖6所示。圖7所示為1=900 rad/s時(shí)，系統(tǒng)復(fù)雜運(yùn)動(dòng)的時(shí)域圖、相圖以及Poincare?截面圖，其中Poincare?截面圖所呈現(xiàn)的封閉圈說(shuō)明發(fā)生在810~900 rad/s高速區(qū)間的復(fù)雜運(yùn)動(dòng)為擬周期運(yùn)動(dòng)。

圖8所示為1=680 rad/s時(shí)，線(xiàn)彈性軸承支撐下齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)的時(shí)域圖、相圖以及Poincare?截面圖，其中Poincare?截面圖所呈現(xiàn)的分形現(xiàn)象說(shuō)明線(xiàn)彈性軸承支撐下齒輪耦合轉(zhuǎn)子系統(tǒng)的復(fù)雜運(yùn)動(dòng)為混沌。

圖4 ω1=550 rad/s時(shí)滑動(dòng)軸承齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)吸引子的形態(tài)

圖5 ω1=620 rad/s時(shí)滑動(dòng)軸承齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)吸引子的形態(tài)

圖6 ω1=770 rad/s時(shí)滑動(dòng)軸承齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)吸引子的形態(tài)

圖7 ω1=900 rad/s時(shí)滑動(dòng)軸承齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)吸引子的形態(tài)

圖8 ω1=680 rad/s時(shí)線(xiàn)彈性軸承支撐下齒輪耦合轉(zhuǎn)子系統(tǒng)穩(wěn)定運(yùn)動(dòng)吸引子的形態(tài)

綜合以上仿真圖形分析可知：在較低轉(zhuǎn)速的工況下，滑動(dòng)軸承支撐下齒輪耦合轉(zhuǎn)子的振動(dòng)表現(xiàn)為形態(tài)各異的周期狀態(tài)，在較高轉(zhuǎn)速工況下系統(tǒng)振動(dòng)則會(huì)呈現(xiàn)出擬周期運(yùn)動(dòng)形態(tài)；滑動(dòng)軸承油膜力對(duì)齒輪耦合轉(zhuǎn)子的混沌運(yùn)動(dòng)具有顯著的鎮(zhèn)定作用。

3 系統(tǒng)參數(shù)對(duì)齒輪副嚙合沖擊特性的影響

3.1 滑動(dòng)軸承間隙的影響規(guī)律

固定轉(zhuǎn)速1=720 rad/s，轉(zhuǎn)子質(zhì)量偏心=0，為了研究方便，主從動(dòng)轉(zhuǎn)子滑動(dòng)軸承間隙取同一值，即1=2=。圖9所示為不同滑動(dòng)軸承間隙下系統(tǒng)齒輪副動(dòng)態(tài)嚙合力的波動(dòng)規(guī)律。

從圖9可以看到：當(dāng)軸承間隙=100 μm時(shí)系統(tǒng)齒輪副能夠正常嚙合；當(dāng)軸承間隙減小到=10 μm時(shí)系統(tǒng)齒輪副發(fā)生了單邊沖擊現(xiàn)象，即齒輪副出現(xiàn)了空嚙合到正常嚙合的循環(huán)沖擊現(xiàn)象；當(dāng)軸承間隙繼續(xù)減小到=8 μm時(shí)，系統(tǒng)齒輪副的單邊沖擊現(xiàn)象進(jìn)一步加劇。以上仿真結(jié)果說(shuō)明：滑動(dòng)軸承間隙對(duì)齒輪的嚙合狀態(tài)影響非常顯著，間隙選取不當(dāng)可能導(dǎo)致齒輪單邊沖擊現(xiàn)象的發(fā)生。

在軸承間隙區(qū)間4~100 μm范圍內(nèi)離散選取有限多個(gè)間隙點(diǎn)，重復(fù)數(shù)值計(jì)算得到每個(gè)軸承間隙值對(duì)應(yīng)的齒輪副嚙合力波動(dòng)的最小值，把上述數(shù)據(jù)點(diǎn)描成曲線(xiàn)即可獲得滑動(dòng)軸承間隙對(duì)齒輪嚙合狀態(tài)具體的影響規(guī)律，如圖10所示。顯然，在4~12 μm之間的較小間隙區(qū)間內(nèi)系統(tǒng)齒輪副動(dòng)態(tài)嚙合力波動(dòng)的最小值為0，意味著該軸承間隙區(qū)間內(nèi)齒輪副出現(xiàn)了單邊沖擊現(xiàn)象；在12~100 μm之間的較大軸承間隙區(qū)間內(nèi)齒輪副正常嚙合，并且隨著間隙的增大嚙合力波動(dòng)的最小值按照一種上凸單調(diào)曲線(xiàn)遞增。

圖9 不同滑動(dòng)軸承間隙時(shí)系統(tǒng)齒輪副的動(dòng)態(tài)嚙合力

圖10 系統(tǒng)齒輪副最小嚙合力隨軸承間隙的變化規(guī)律

3.2 轉(zhuǎn)子偏心的影響規(guī)律

以下固定轉(zhuǎn)速1=720 rad/s，主從動(dòng)轉(zhuǎn)子的軸承間隙均為50 μm，為了研究方便，主從動(dòng)轉(zhuǎn)子的質(zhì)量偏心假定為同一值，即1=2=。圖11所示為不同質(zhì)量偏心下系統(tǒng)齒輪副的動(dòng)態(tài)嚙合力的波動(dòng)規(guī)律，

從圖11可以看到，當(dāng)轉(zhuǎn)子質(zhì)量偏心取=3.0 μm時(shí)系統(tǒng)齒輪副能夠正常嚙合；當(dāng)質(zhì)量偏心增大為=3.1 μm時(shí)系統(tǒng)齒輪副開(kāi)始出現(xiàn)單邊沖擊現(xiàn)象；當(dāng)質(zhì)量偏心增大到=4.0 μm時(shí)，系統(tǒng)齒輪副的單邊沖擊現(xiàn)象進(jìn)一步加劇。以上仿真結(jié)果說(shuō)明，轉(zhuǎn)子質(zhì)量偏心對(duì)齒輪的嚙合狀態(tài)也有著非常顯著的影響，間隙值選取不當(dāng)可能導(dǎo)致齒輪副單邊沖擊現(xiàn)象的發(fā)生。

在轉(zhuǎn)子質(zhì)量偏心區(qū)間0~4 μm范圍內(nèi)離散選取有限多個(gè)偏心點(diǎn)，重復(fù)數(shù)值計(jì)算得到每個(gè)轉(zhuǎn)子質(zhì)量偏心對(duì)應(yīng)的齒輪副嚙合力波動(dòng)的最小值，把上述數(shù)據(jù)點(diǎn)描成曲線(xiàn)即可獲得轉(zhuǎn)子質(zhì)量偏心對(duì)齒輪嚙合狀態(tài)具體的影響規(guī)律，如圖12所示。顯然，在3.1~4.0 μm之間的較大質(zhì)量偏心區(qū)間內(nèi)系統(tǒng)齒輪副動(dòng)態(tài)嚙合力波動(dòng)的最小值為0，意味著該偏心區(qū)間內(nèi)齒輪副出現(xiàn)了單邊沖擊現(xiàn)象；在0~3.0 μm之間較小的偏心區(qū)間內(nèi)齒輪副正常嚙合。需要特別說(shuō)明的是，隨著轉(zhuǎn)子偏心量的逐漸增大，齒輪副動(dòng)態(tài)嚙合力的最小值并非逐漸減小而達(dá)到零值，而是在偏心增大到3.1 μm附近時(shí)，經(jīng)過(guò)一個(gè)“階躍”突然歸零，所以找出“階躍”發(fā)生時(shí)的偏心閥值對(duì)避免齒輪副單邊沖擊現(xiàn)象的發(fā)生有著重要的價(jià)值。

圖11 不同質(zhì)量偏心下系統(tǒng)齒輪副的動(dòng)態(tài)嚙合力

圖12 系統(tǒng)齒輪副最小嚙合力隨質(zhì)量偏心的變化

4 結(jié)論

1) 在綜合考慮滑動(dòng)軸承的非線(xiàn)性油膜力以及齒輪副的嚙合沖擊力等非線(xiàn)性因素的基礎(chǔ)上，建立了滑動(dòng)軸承支撐下齒輪耦合雙軸轉(zhuǎn)子系統(tǒng)的非線(xiàn)性動(dòng)力學(xué)模型，推導(dǎo)了系統(tǒng)關(guān)于齒輪副嚙合點(diǎn)相對(duì)位移的運(yùn)動(dòng)微分方程。

2) 在較低轉(zhuǎn)速的工況下，滑動(dòng)軸承支撐下齒輪耦合轉(zhuǎn)子的振動(dòng)表現(xiàn)為形態(tài)各異的周期狀態(tài)，在較高轉(zhuǎn)速工況下系統(tǒng)振動(dòng)則會(huì)呈現(xiàn)出擬周期運(yùn)動(dòng)形態(tài)；滑動(dòng)軸承油膜力對(duì)齒輪耦合轉(zhuǎn)子的混沌運(yùn)動(dòng)具有顯著的鎮(zhèn)定作用。

3) 滑動(dòng)軸承間隙以及轉(zhuǎn)子質(zhì)量偏心對(duì)齒輪的嚙合狀態(tài)影響非常顯著，軸承間隙以及轉(zhuǎn)子偏心量選取不當(dāng)會(huì)導(dǎo)致齒輪單邊沖擊現(xiàn)象的發(fā)生；轉(zhuǎn)子質(zhì)量偏心所導(dǎo)致的單邊沖擊具有突然“階躍”發(fā)生的特征，所以找出“階躍”發(fā)生時(shí)的偏心閥值對(duì)避免齒輪單邊沖擊現(xiàn)象的發(fā)生有著重要的指導(dǎo)意義。

[1] 唐進(jìn)元, 陳思雨, 鐘掘. 一種改進(jìn)的齒輪非線(xiàn)性動(dòng)力學(xué)模型[J]. 工程力學(xué), 2008, 25(1): 217?223.TANG Jinyuan, CHEN Siyu, ZHONG Jue. An improved nonlinear model for a spur gear pair system[J]. Engineering Mechanics, 2008, 25(1): 217?223.

[2] LEE A S, HA J W, CHOI D H. Coupled lateral and torsional vibration characteristics of a speed increasing geared rotor bearing system[J]. Journal of Sound and Vibration, 2003, 263: 725?742.

[3] 陳兵奎, 易文翠, 鐘暉, 等. 線(xiàn)面共軛嚙合原理及齒面構(gòu)建方法[J]. 機(jī)械工程學(xué)報(bào), 2012, 48(19): 17?22. CHEN Bingkui, YI Wencui, ZHONG Hui, et al. Theory of curve-surface Conjugated and method of teeth surfaces constructed[J]. Journal of Mechanical Engineer, 2012, 48(19): 17?22.

[4] HSIEH S C, CHEN J H, LEE A. A modified transfer matrix method for the coupled lateral and torsional vibrations of symmetric rotor-bearing systems[J]. Journal of Sound and Vibration, 2006, 289(1/2): 294?333

[5] CHEN Z, SHAO Y. Dynamic simulation of spur gear with tooth root crack propagating along tooth width and crack depth[J]. Engineering Failure Analysis, 2011, 18(8): 2149?2164.

[6] CHEN Z, SHAO Y. Mesh stiffness calculation of a spur gear pair with tooth profile modification and tooth root crack[J]. Mechanism and Machine Theory, 2013, 62: 63?74.

[7] ABBES M S, TRIGUI M, CHAARI F, et al. Dynamic behaviour modeling of a flexible gear system by the elastic foundation theory in presence of defects[J]. European Journal of Mechanics-A/Solids, 2010, 29(5): 887?896.

[8] ZHANG Y, WANG Q, MA H, et al. Dynamic analysis of three-dimensional helical geared rotor system with geometric eccentricity[J]. Journal of Mechanical Science and Technology, 2013, 27(11): 3231?3242.

[9] 李立, 鄭鐵生, 許慶余. 齒輪–轉(zhuǎn)子–滑動(dòng)軸承系統(tǒng)時(shí)變非線(xiàn)性動(dòng)力特性研究[J]. 應(yīng)用力學(xué)學(xué)報(bào), 1995, 12(1): 15?23.LI Li, ZHENG Tiesheng, XU Qingyu. The study of nonlinear varying characteristics in the geared rotor-bearing system[J]. Chinese Journal of Applied Mechanics, 1995, 12(1): 15?23.

[10] 張鎖懷, 謝振宇, 丘大謀. 軸承參數(shù)對(duì)平行轉(zhuǎn)子?軸承系統(tǒng)動(dòng)力特性的影響[J]. 機(jī)械工程學(xué)報(bào), 2000, 36(9): 16?19.ZHANG Suohuai, XIE Zhenyu, QIU Damou. Influence of the bearing parameters on dynamic characteristic of a parallel rotor-bearing system[J]. Chinese Journal of Mechanical Engineering, 2000, 36(9): 16?19.

[11] KUBUR M, KAHRAMA A, ZINI D M. Dynamic analysis of a multi-shaft helical gear transmission by finite elements: model and experiment[J]. Transactions of the American Society Mechanical Engineering, 2004, 126: 398?406.

[12] MAHARATHI B B. Dynamic behaviour analysis of linear rotor-bearing systems using the complex transfer matrix technique[J]. International Journal of Acoustics and Vibrations, 2005, 10(3): 413?417.

[13] BAGUET S, JACQUENOT G. Nonlinear couplings in a gear-shaft bearing system[J]. Mechanism and Machine Theory, 2010, 45(12): 1777?1796.

[14] THEODOSSIADEA S, NATSIAVAS S. On geared rotor dynamic systems with oil journal bearings[J]. Journal of Sound and Vibration, 2001, 243(4): 721?745.

[15] GAO H, ZHANG Yidu. Nonlinear behavior analysis of geared rotor bearing system featuring confluence transmission[J]. Nonlinear Dynamics, 2014, 76: 2025?2039.

[16] 李同杰, 朱如鵬, 鮑和云, 等. 行星齒輪系扭轉(zhuǎn)非線(xiàn)性振動(dòng)建模與運(yùn)動(dòng)分岔特性研究[J]. 機(jī)械工程學(xué)報(bào), 2011, 47(21): 76?83. LI Tongjie, ZHU Rupeng, BAO Heyun, et al. Nonlinear torsional vibration modeling and bifurcation characteristic study of a planetary gear train[J]. Journal of Mechanical Engineering, 2011, 47(21): 76?83.

[17] ADILETTA G., GUIDO A R, ROSSI C. Chaotic motions of a rigid rotor in short journal bearings[J]. Nonlinear Dynamics, 1996, 10(3): 251?269.

[18] 焦映厚, 陳照波, 劉福利, 等. Jeffcott轉(zhuǎn)子?可傾瓦滑動(dòng)軸承系統(tǒng)不平衡響應(yīng)的非線(xiàn)性分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2004, 24(12): 227?232.JIAO Yinghou, CHEN Zhaobo, LIU Fuli, et al. Nonlinear analysis of unbalance response for Jeffcott rotor-tilting pad bearing system[J]. Proceedings of the CSEE, 2004, 24(12): 227?232.

[19] 成大先. 機(jī)械設(shè)計(jì)手冊(cè)[M]. 5版. 北京: 化學(xué)工業(yè)出版社, 2008.CHENG Daxian. Hand book of mechanical design[M]. 5th ed. Beijing: Chemical Industry Press, 2008.

(編輯 趙俊)

Nonlinear bending-torsional dynamics of geared rotor system supported by sliding bearing

LI Tongjie1, 2, JIN Guanghu1, ZHU Rupeng1, AN Luling1

(1. College of Mechanical and Electrical Engineering, Nanjing University of Aeronautics and Astronautics,Nanjing 210016, China; 2. College of Mechanical Engineering, Anhui Science & Technology University, Chuzhou 233100, China)

A nonlinear dynamic model of geared rotor system coupling with sliding bearing was established taking nonlinear oil film force of the bearing and backlash of the system into account. By using the method of numerical simulation, some influence laws of the system parameters on meshing impact properties of the gear pair and system bifurcation characteristics about motion stability were studied based on the nonlinear dynamic model. The results reveal that slide bearing nonlinear oil film force can calm the chaos of the system, and improper sliding bearing clearance value and rotor eccentricity will lead to a unilateral impact phenomenon.

geared rotor system; sliding bearing; coupled nonlinear dynamic model; meshing impact; bifurcation

TH132

A

10.11817/j.issn.1672-7207.2018.03.008

1672?7207(2018)03?0566?08

2017?03?09；

2017?05?18

國(guó)家自然科學(xué)基金資助項(xiàng)目(51475226)；安徽省教育廳自然科學(xué)重點(diǎn)項(xiàng)目(KJ2015A179)；安徽省科技公關(guān)計(jì)劃項(xiàng)目(1604a0902134，1501031095) (Project(51475226) supported by the National Natural Science Foundation of China; Project(KJ2015A179) supported by Key Natural Science Fund of Education Department of Anhui Province; Projects(1604a0902134, 1501031095) supported by Science and Technology Agency of Anhui Province)

李同杰，博士，副教授，從事齒輪傳動(dòng)系統(tǒng)動(dòng)力學(xué)研究；E-mail: litongjie2000@163.com