花純利， 塔 娜， 饒柱石

(上海交通大學(xué) 振動(dòng)、沖擊噪聲研究所 機(jī)械系統(tǒng)與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，上海 200240)

轉(zhuǎn)軸系統(tǒng)廣泛用于實(shí)際工程。為提高轉(zhuǎn)軸系統(tǒng)機(jī)械效率，將轉(zhuǎn)軸與軸承間隙設(shè)計(jì)的越來(lái)越小，但會(huì)使轉(zhuǎn)軸系統(tǒng)發(fā)生碰摩的可能性提高。轉(zhuǎn)軸與軸承間發(fā)生碰摩將致轉(zhuǎn)軸系統(tǒng)局部發(fā)熱甚至嚴(yán)重磨損，誘發(fā)機(jī)械劇烈振動(dòng)，嚴(yán)重時(shí)會(huì)出現(xiàn)反向渦動(dòng)失穩(wěn)造成整個(gè)機(jī)械系統(tǒng)破壞。因此，研究轉(zhuǎn)軸系統(tǒng)的碰摩復(fù)雜非線性動(dòng)力學(xué)行為、確定其與系統(tǒng)參數(shù)間關(guān)系、揭示系統(tǒng)穩(wěn)定性邊界條件，對(duì)優(yōu)化轉(zhuǎn)軸系統(tǒng)設(shè)計(jì)與故障診斷意義甚為重要。

對(duì)轉(zhuǎn)軸系統(tǒng)碰摩研究已有大量成果。自上世紀(jì)80年代以來(lái)，在實(shí)驗(yàn)、數(shù)值模擬、理論分析等方面對(duì)轉(zhuǎn)子系統(tǒng)碰摩響應(yīng)特性進(jìn)行廣泛深入研究。如同頻全周碰摩運(yùn)動(dòng)、“跳躍”現(xiàn)象[1-6]，準(zhǔn)周期局部碰摩運(yùn)動(dòng)[7]及混沌行為[8-10]。文獻(xiàn)[11-13]中所用模型雖與本文相似，但其接觸力為線性模型且未深入討論偏心率對(duì)系統(tǒng)動(dòng)態(tài)特性影響，亦未解釋、闡明跳躍現(xiàn)象產(chǎn)生的條件。文獻(xiàn)[2]亦用兩自由度轉(zhuǎn)子模型，但其定子具有耦合剛度且接觸力為線性模型。止今，大部分研究主要集中在將軸承簡(jiǎn)化為無(wú)質(zhì)量線性彈簧的轉(zhuǎn)軸碰摩模型上，但該簡(jiǎn)化不能充分反映具有超彈性材料特性的橡膠軸承動(dòng)態(tài)特性。國(guó)內(nèi)外船舶及深井泵等設(shè)備中因大量使用橡膠軸承，致使橡膠體磨損嚴(yán)重。此與螺旋槳軸承支承的整個(gè)轉(zhuǎn)子系統(tǒng)靜態(tài)特性與動(dòng)力響應(yīng)密切相關(guān)。由于對(duì)橡膠軸承支撐的轉(zhuǎn)軸/軸承系統(tǒng)研究較少，且對(duì)全面考慮碰摩主要參數(shù)的非線性彈性支撐的轉(zhuǎn)軸/軸承系統(tǒng)模型動(dòng)力學(xué)特性認(rèn)識(shí)不足，尤其對(duì)轉(zhuǎn)軸/軸承偏心率及阻尼比對(duì)系統(tǒng)動(dòng)力響應(yīng)所致影響了解更少。因此，本文以橡膠軸承支撐的轉(zhuǎn)軸/軸承系統(tǒng)為對(duì)象，將其簡(jiǎn)化為非線性彈性支撐的轉(zhuǎn)軸/軸承系統(tǒng)進(jìn)行研究，并分析該模型動(dòng)力學(xué)特性。即分析轉(zhuǎn)軸系統(tǒng)發(fā)生碰摩、鞍結(jié)分叉及Hopf分叉的邊界條件，討論阻尼比與偏心率等系統(tǒng)參數(shù)對(duì)系統(tǒng)穩(wěn)定性影響。

1 模型與運(yùn)動(dòng)控制方程

橡膠軸承/轉(zhuǎn)軸系統(tǒng)示意圖見(jiàn)圖1，橡膠軸承由軸承襯套與橡膠襯套兩部分組成。將其簡(jiǎn)化成Jeffcott轉(zhuǎn)子系統(tǒng)，見(jiàn)圖2，轉(zhuǎn)軸簡(jiǎn)化為一支撐在無(wú)質(zhì)量剛度k，阻尼c，轉(zhuǎn)軸中間質(zhì)量m的剛性轉(zhuǎn)軸上，轉(zhuǎn)軸與軸承間間隙δ，轉(zhuǎn)軸質(zhì)心與幾何中心偏心距e。

圖1 轉(zhuǎn)軸/軸承系統(tǒng)示意圖

圖2 Jeffcott轉(zhuǎn)軸/軸承系統(tǒng)模型

考慮橡膠軸承具有的非線性特征，計(jì)入支承剛度的非線性特性。由于軸系的橡膠艉軸承均為新設(shè)計(jì)產(chǎn)品，未經(jīng)實(shí)驗(yàn)測(cè)得準(zhǔn)確剛度曲線。因此，其非線性載荷-變形關(guān)系參照德國(guó)國(guó)防軍艦艇建造規(guī)范(BV043)進(jìn)行估算：

F=krr(1+100r)=krr+αr2

(1)

其中：F為載荷(N)；r為變形量(m)；kr為軸承剛度值(N/m)。故軸承內(nèi)環(huán)面上轉(zhuǎn)軸/軸承接觸剛度為kr+αr。則轉(zhuǎn)軸與軸承間的摩擦力及接觸力為：

(2)

(3)

橡膠軸承/轉(zhuǎn)軸系統(tǒng)動(dòng)力學(xué)方程可表達(dá)為：

(4)

式中：Θ為heaveside函數(shù)，即：

將橡膠軸承/轉(zhuǎn)軸系統(tǒng)動(dòng)力學(xué)方程無(wú)量綱化：

(5)

轉(zhuǎn)軸旋轉(zhuǎn)過(guò)程中，轉(zhuǎn)軸與軸承接觸、不接觸狀態(tài)均有穩(wěn)定的周期解。設(shè)其解的形式為：

(6)

在軸承與轉(zhuǎn)軸未發(fā)生接觸狀態(tài)下Θ=0，將式(6)代入控制方程并求解得幅值、相位角分別為：

(7)

由于軸承與轉(zhuǎn)軸間隙有限，故非接觸狀態(tài)下所求幅值A(chǔ)須滿足A≤1條件。即：

(ρ2-1)Ω4+2(β-2ξ2)Ω2-β2≤0

(8)

通過(guò)求解式(8)可兩實(shí)根Ωl，Ωu，且記為線性轉(zhuǎn)軸開(kāi)始發(fā)生碰摩時(shí)低、高轉(zhuǎn)速。當(dāng)Ω<Ωl或Ω>Ωu時(shí)，轉(zhuǎn)軸系統(tǒng)存在無(wú)碰摩同頻周期響應(yīng)。

軸承與轉(zhuǎn)軸發(fā)生接觸狀態(tài)下Θ=1，將式(6)代入式(5)得關(guān)于振動(dòng)幅值A(chǔ)的方程式：

a4A4+a3A3+a2A2+a1A+a0=0

(9)

其中：

a4=g2+g2f2

a3=2g[(1-2g+β-Ω2)+f(2ξΩ+f-2gf)]

a2=(1-2g+β-Ω2)2+(2ξΩ+f-2gf)2+

2g(g-1)(1+f2)

a1=2(g-1)[(1-2g+β-Ω2)+f(2ξΩ+f-2gf)]

a0=(g-1)2(1+f2)-ρ2Ω4。

一元四次方程最多有兩個(gè)互異的正實(shí)數(shù)根，復(fù)數(shù)與負(fù)數(shù)實(shí)根均無(wú)實(shí)際意義。在式(5)Θ=1條件下，方程有二重正實(shí)根時(shí)，即為系統(tǒng)鞍結(jié)分叉的邊界條件。為保證軸承與轉(zhuǎn)軸處于接觸狀態(tài)，幅值A(chǔ)不僅需滿足正實(shí)數(shù)條件且需滿足幅值A(chǔ)>1。

2 周期解穩(wěn)定性分析

(10)

(11)

式中：

為簡(jiǎn)便，引入穩(wěn)態(tài)周期解形式：

(12)

其中：A為非接觸狀態(tài)下式(7)與接觸狀態(tài)下式(9)的穩(wěn)定周期解。

通過(guò)解式(11)得：

(13)

當(dāng)Θ=0時(shí)轉(zhuǎn)軸與軸承處于非接觸狀態(tài)，雅可比方程J恰好是矩陣B，對(duì)應(yīng)的特征方程為|B-λI|=0，將其展開(kāi)：

λ4+4ξλ3+(2β+4ξ2)λ2+4ξβλ+β2=0

(14)

根據(jù)Routh-Hurwitz(勞斯-霍爾維茨)穩(wěn)定性判據(jù)，式(5)(Θ=0時(shí))非接觸狀態(tài)解的穩(wěn)定性條件為：

ξ>0

(15)

由式(15)知，當(dāng)系統(tǒng)阻尼為正阻尼時(shí)，對(duì)應(yīng)的周期解是穩(wěn)定的。當(dāng)系統(tǒng)振動(dòng)幅值大于間隙δ時(shí)，軸承與轉(zhuǎn)軸會(huì)發(fā)生碰摩。此時(shí)轉(zhuǎn)軸系統(tǒng)控制方程中Θ=1。式(5)的解將為非線性的周期解。其雅可比矩陣J可表達(dá)為：

(16)

式中：

由式(16)知，雅可比矩陣J為周期性的時(shí)間依賴(lài)矩陣，不能直接推導(dǎo)、分析其解的穩(wěn)定性，需作變換：

(17)

式中：轉(zhuǎn)換矩陣T為：

將式(17)代入式(13)，得：

(18)

(19)

式(19)表明雅可比矩陣Jc與時(shí)間參數(shù)無(wú)關(guān)。δU解與式(11)解的穩(wěn)定性取決于矩陣Jc特征值實(shí)部符號(hào)。對(duì)應(yīng)的特征方程滿足|Jc-λI|=0，將其展開(kāi)為：

b4λ4+b3λ3+b2λ2+b1λ+b0=0

(20)

其中：

b4=1

b3=4ξ

諸多系數(shù)均為幅值A(chǔ)的函數(shù)，故Routh-Hurwitz(勞斯-霍爾維茨)穩(wěn)定性判據(jù)可用于判斷式(11)非線性穩(wěn)態(tài)周期解的穩(wěn)定性。

基于分叉理論,分析周期解的分叉邊界，給出參數(shù)空間穩(wěn)定區(qū)域。若雅可比矩陣Jc有一零特征值，系統(tǒng)將出現(xiàn)鞍結(jié)分叉，此時(shí)式(20)中b0=0，即：

(21)

通過(guò)消除幅值A(chǔ)的符號(hào)計(jì)算，同時(shí)求解式(9)、(21)獲得關(guān)于Ω的12次多項(xiàng)式。求解參數(shù)方程，可得式(9)發(fā)生鞍結(jié)分叉條件的參數(shù)空間。全周碰摩解的鞍結(jié)分叉點(diǎn)處幅值A(chǔ)為大于1的正實(shí)數(shù)。

基于Hopf分叉理論，系統(tǒng)會(huì)有一對(duì)共軛純虛數(shù)特征值。將λ=+iωυ代入式(20)得：

(22)

消去參數(shù)ωυ得(代入λ=-iωυ得同樣結(jié)果)：

(23)

需滿足不等式：

b1/b3>0

(24)

用式(9)、(23)消去幅值A(chǔ)，求解參數(shù)方程，可得系統(tǒng)發(fā)生Hopf分叉的邊界線HF。同頻全周碰摩解Hopf分叉點(diǎn)處幅值A(chǔ)為大于1的正實(shí)數(shù)。

3 系統(tǒng)參數(shù)對(duì)轉(zhuǎn)軸響應(yīng)影響

由式(15)知，轉(zhuǎn)軸系統(tǒng)在無(wú)碰摩狀態(tài)下周期解總是穩(wěn)定的，但轉(zhuǎn)軸與軸承間間隙有限，因此周期解不可能在系統(tǒng)的所有參數(shù)下一直存在。由式(8)及分析知，線性橡膠軸承/轉(zhuǎn)軸系統(tǒng)激振頻率分別由高、低轉(zhuǎn)速趨向系統(tǒng)固有頻率時(shí)，系統(tǒng)振動(dòng)幅值A(chǔ)會(huì)在Ω=Ωl與Ω=Ωu處趨近1。當(dāng)轉(zhuǎn)軸/軸承系統(tǒng)處于接觸狀態(tài)、激振頻率分別向低、高轉(zhuǎn)速變化時(shí)，系統(tǒng)振動(dòng)幅值A(chǔ)會(huì)在Ω=SNl與Ω=SNu處趨近1。由式(9)、(21)知，碰摩狀態(tài)存在于轉(zhuǎn)速SNl與SNu之間，SNl與Ωl基本相等，SNu與Ωu或相等或不等，取決于系統(tǒng)參數(shù)，兩者不等時(shí)系統(tǒng)振動(dòng)幅值將隨轉(zhuǎn)速的變化出現(xiàn)跳躍現(xiàn)象。如將曲線Ω=Ωl，Ω=Ωu，Ω=SNl，Ω=SNu，Hopf分叉邊界HF繪制于(Ω,ξ)平面，則該平面將被分隔成若干部分，不同部分代表系統(tǒng)不同運(yùn)動(dòng)狀態(tài)，即(A)，(B)為無(wú)碰摩周期運(yùn)動(dòng)；(C)為同頻全周碰摩運(yùn)動(dòng)；(D)為局部碰摩運(yùn)動(dòng)及反向渦動(dòng)失穩(wěn)運(yùn)動(dòng)；(E)為無(wú)碰摩周期運(yùn)動(dòng)及同頻全周碰摩運(yùn)動(dòng)共存，即處于接觸狀態(tài)條件下系統(tǒng)做同頻全周碰摩運(yùn)動(dòng)，而處于非接觸狀態(tài)條件下系統(tǒng)做無(wú)碰摩周期運(yùn)動(dòng)；(F)為無(wú)碰摩周期運(yùn)動(dòng)、局部碰摩運(yùn)動(dòng)或反向渦動(dòng)失穩(wěn)運(yùn)動(dòng)，即處于接觸狀態(tài)條件下系統(tǒng)做局部碰摩運(yùn)動(dòng)或反向渦動(dòng)失穩(wěn)運(yùn)動(dòng)，而處于非接觸狀態(tài)條件下系統(tǒng)做無(wú)碰摩周期運(yùn)動(dòng)。

對(duì)不同參數(shù)取值下系統(tǒng)運(yùn)動(dòng)狀態(tài)演變過(guò)程分別進(jìn)行討論。取偏心率ρ=0.8，系統(tǒng)各參數(shù)分別為：f=0.2,β=0.5，g=0.01。此時(shí)曲線Ω=Ωl，Ω=Ωu，Ω=SNl，Ω=SNu，分叉邊界HF將(Ω,ξ)平面分成四部分，轉(zhuǎn)軸系統(tǒng)在轉(zhuǎn)速Ω，阻尼比ξ參數(shù)平面內(nèi)系統(tǒng)響應(yīng)特性區(qū)域見(jiàn)圖3。因Ωl=SNl，Ωu=SNu，系統(tǒng)不會(huì)發(fā)生跳躍現(xiàn)象。由圖3知，在區(qū)域(C)內(nèi)轉(zhuǎn)軸系統(tǒng)做同頻全周碰摩運(yùn)動(dòng)，在區(qū)域(D)內(nèi)系統(tǒng)做非穩(wěn)定周期碰摩運(yùn)動(dòng)。

將系統(tǒng)偏心率減小為ρ=0.5時(shí)，轉(zhuǎn)軸系統(tǒng)響應(yīng)穩(wěn)定區(qū)域發(fā)生明顯變化，見(jiàn)圖4。

與偏心率ρ=0.8情況相比，系統(tǒng)響應(yīng)特性在(Ω,ξ)平面內(nèi)發(fā)生顯著變化，邊界Ω=Ωl=SNl向右移動(dòng)，邊界Ω=Ωu與Ω=SNu向左移動(dòng)，且曲線Ωu，SNu不重合。Hopf分叉邊界線HF無(wú)明顯變化。因此，非穩(wěn)定碰摩響應(yīng)區(qū)域(D)及同頻全周碰摩運(yùn)動(dòng)區(qū)域(C)的兩側(cè)邊界向中間靠攏而縮小，表明系統(tǒng)出現(xiàn)全周碰摩轉(zhuǎn)速范圍變?。怀霈F(xiàn)同頻全周碰摩運(yùn)動(dòng)與無(wú)碰摩周期運(yùn)動(dòng)共存區(qū)域(E)及無(wú)碰摩周期運(yùn)動(dòng)、局部碰摩運(yùn)動(dòng)或反向渦動(dòng)失穩(wěn)運(yùn)動(dòng)區(qū)域(F)。此時(shí)轉(zhuǎn)軸系統(tǒng)同頻全周碰摩運(yùn)動(dòng)會(huì)存在于Ω=Ωl與Ω=SNu區(qū)間內(nèi)；區(qū)間Ω>Ωu存在無(wú)碰摩周期運(yùn)動(dòng)，故在區(qū)間Ωu，SNu間同頻全周碰摩運(yùn)動(dòng)與無(wú)碰摩周期運(yùn)動(dòng)均可能存在。當(dāng)轉(zhuǎn)軸系統(tǒng)轉(zhuǎn)速Ω發(fā)生變化時(shí)會(huì)有跳躍現(xiàn)象發(fā)生。

由討論可知，偏心率降低會(huì)使碰摩區(qū)域發(fā)生明顯變化并出現(xiàn)跳躍現(xiàn)象，因此，需進(jìn)一步討論其對(duì)轉(zhuǎn)軸系統(tǒng)響應(yīng)影響。轉(zhuǎn)軸系統(tǒng)在(Ω,ρ)平面的響應(yīng)特性見(jiàn)圖5。

圖3 轉(zhuǎn)軸碰摩響應(yīng)在(Ω, ξ)平面穩(wěn)定區(qū)域圖(ρ=0.8)

由式(8)所得Ωu，Ωl值可給出轉(zhuǎn)軸與軸承由無(wú)接觸向接觸過(guò)渡時(shí)臨界轉(zhuǎn)速上、下邊界。由圖5看出，①系統(tǒng)共振振幅小于轉(zhuǎn)軸與軸承間間隙即ρ<0.17時(shí)，轉(zhuǎn)軸與軸承間不發(fā)生碰摩；隨偏心率的增加Ωl逐漸減小而Ωu卻逐漸增大，因此Ωl與Ωu間區(qū)域增大，即發(fā)生碰摩的轉(zhuǎn)速區(qū)間增大。②同理可得轉(zhuǎn)軸與軸承由接觸狀態(tài)向脫離接觸狀態(tài)過(guò)渡時(shí)轉(zhuǎn)速的上、下邊界值SNu及SNl，而Ωl與SNl總相等；Ωu與SNu在0.17<ρ<0.65范圍內(nèi)不相等，該區(qū)域內(nèi)轉(zhuǎn)軸系統(tǒng)振幅隨轉(zhuǎn)速的變化會(huì)發(fā)生跳躍現(xiàn)象，且隨偏心率的增加發(fā)生跳躍現(xiàn)象的轉(zhuǎn)速區(qū)間呈現(xiàn)先增大后減小直至消失的趨勢(shì)；③由于轉(zhuǎn)軸系統(tǒng)Hopf分叉邊界HF與鞍結(jié)分叉邊界SNu在ρ=0.3處相交一點(diǎn)，ρ>0.3時(shí)，隨偏心率的增加，開(kāi)始發(fā)生Hopf分叉的轉(zhuǎn)速Ω逐漸減小，導(dǎo)致落在區(qū)域(E)內(nèi)的轉(zhuǎn)速區(qū)間不斷縮小并最終消失；0.17<ρ<0.3時(shí)，系統(tǒng)不會(huì)出現(xiàn)Hopf分叉及失穩(wěn)現(xiàn)象，只會(huì)出現(xiàn)無(wú)碰摩周期運(yùn)動(dòng)、同頻全周碰摩運(yùn)動(dòng)或兩運(yùn)動(dòng)狀態(tài)跳躍轉(zhuǎn)變現(xiàn)象。

4 結(jié) 論

本文對(duì)軸承簡(jiǎn)化為非線性彈性支撐轉(zhuǎn)軸系統(tǒng)發(fā)生周期無(wú)碰摩運(yùn)動(dòng)與同頻全周碰摩運(yùn)動(dòng)存在的區(qū)域及穩(wěn)定性進(jìn)行研究，并分析阻尼比及偏心率等系統(tǒng)參數(shù)對(duì)轉(zhuǎn)軸系統(tǒng)動(dòng)態(tài)特性影響、對(duì)不同轉(zhuǎn)軸系統(tǒng)響應(yīng)特性邊界條件影響。結(jié)論如下：

(1)轉(zhuǎn)軸系統(tǒng)Hopf分叉邊界HF為轉(zhuǎn)速、摩擦系數(shù)、阻尼比及偏心率的函數(shù)。Hopf分叉邊界將碰摩響應(yīng)區(qū)域劃分成同頻全周碰摩響應(yīng)區(qū)域與局部碰摩運(yùn)動(dòng)、反向渦動(dòng)失穩(wěn)運(yùn)動(dòng)區(qū)域兩部分。

(2)阻尼比是影響轉(zhuǎn)軸碰摩響應(yīng)的主要參數(shù)之一。大阻尼比系數(shù)下，轉(zhuǎn)軸系統(tǒng)響應(yīng)主要為無(wú)碰摩周期運(yùn)動(dòng)、同頻全周碰摩運(yùn)動(dòng)及兩種響應(yīng)共存；小阻尼比系數(shù)會(huì)出現(xiàn)不穩(wěn)定碰摩運(yùn)動(dòng)。因此，大阻尼比轉(zhuǎn)軸系統(tǒng)更穩(wěn)定。

(3)偏心率是轉(zhuǎn)軸系統(tǒng)重要參數(shù)之一，偏心率小轉(zhuǎn)軸系統(tǒng)穩(wěn)定。偏心率小至轉(zhuǎn)軸系統(tǒng)共振幅值小于間隙時(shí)，轉(zhuǎn)軸與軸承間不發(fā)生碰摩。隨偏心率的增加轉(zhuǎn)軸與軸承會(huì)在某狀態(tài)下發(fā)生碰摩；小偏心率時(shí)，系統(tǒng)有周期無(wú)碰摩運(yùn)動(dòng)及同頻全周碰摩運(yùn)動(dòng)且會(huì)發(fā)生在兩運(yùn)動(dòng)狀態(tài)間跳躍現(xiàn)象，較大偏心時(shí)系統(tǒng)會(huì)在某轉(zhuǎn)速下發(fā)生Hopf分叉,且隨偏心率的增大,發(fā)生Hopf分叉所需轉(zhuǎn)速會(huì)降低。

參 考 文 獻(xiàn)

[1]Ehrich F F. High order sub-harmonic response of high speed rotors in bearing clearance[J]. Journal of Vibration, Acoustics, Stress, and Reliability in Design, 1988, 110: 9-16.

[2]Jiang J,Ulbrich H. Stability analysis of sliding whirl in a nonlinear jeffcott rotor with cross-coupling stiffness coefficients[J]. Nonlinear Dynamics,2001,24(3):269-283.

[3]Zhang H B, Chen Y S, Li J. Bifurcation on synchronous full annular rub of rigid-rotor elastic-support system[J]. Applied Mathematics and Mechanics,2012,33(7):865-880.

[4]Zhang H B, Chen Y S. Bifurcation analysis on full annular rub of a nonlinear rotor system [J]. Science China Technological Sciences, 2011, 54(8): 1977-1985.

[5]Shang Z Y, Jiang J, Hong L. The global responses characteristics of a rotor/stator rubbing system with dry friction effects [J]. Journal of Sound and Vibration, 2011, 330(10) :2150-2160.

[6]李朝峰，戴繼雙，聞邦椿. 油膜支撐雙盤(pán)轉(zhuǎn)子-軸承系統(tǒng)周期運(yùn)動(dòng)穩(wěn)定性與分岔[J]. 力學(xué)學(xué)報(bào),2011,43(1):208-216.

LI Chao-feng, DAI Ji-shuang, WEN Bang-chun. Stability and bifurcation of the rotor-bearing system with double disks supported by cylindrical bearings[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(1):208-216.

[7]Kim Y B, Noah S T. Quasi-periodic response and stability analysis for a nonlinear jeffcott rotor[J]. Journal of Sound and Vibration, 1996, 190(2): 239-253.

[8]Goldman P, Muszynska A. Chaotic behavior of rotor/stator system with rubs[J]. Journal of Engineering for Gas Turbines and Power, 1994, 116: 693-701.

[9]Evgueni V K, Wiercigroch M, Cartmell M P. Regular and chaotic dynamics of a discontinuously nonlinear rotor system[J]. Chaos,Solitons and Fractals,2002,13(6): 1231-1242.

[10]褚福磊, 張正松, 馮冠平. 碰摩轉(zhuǎn)子系統(tǒng)的混沌特[J]. 清華大學(xué)學(xué)報(bào)( 自然科學(xué)版),1996,36(7):52-57.

CHU Fu-lei, ZHANG Zheng-song, FENG Guan-ping. Chaotic behavior of a rub rotor model[J]. Journal of Tsinghua University ( Sci & Tech), 1996,36(7):52-57.

[11]劉 林, 江 俊. 轉(zhuǎn)子/定子碰摩響應(yīng)的全局動(dòng)力學(xué)特性研究[J]. 應(yīng)用力學(xué)學(xué)報(bào),2006, 23(3):351-356.

LIU Lin, JIANG Jun. Global dynamical characteristics of rotor/ stator rubbing [J]. Chinese Journal of Applied Mechanics, 2006, 23(3):351-356.

[12]Zhang W M, Meng G. Nonlinear dynamics of a rub-impact micro-rotor system with scale-dependent friction model [J]. Journal of Sound and Vibration, 2008, 309(3-5): 756-777.

[13]張義民, 劉巧伶, 聞邦椿. 旋轉(zhuǎn)機(jī)械系統(tǒng)碰摩故障的失效靈敏度研究[J]. 振動(dòng)工程學(xué)報(bào),2007,20(2):189-193.

ZHANG Yi-min, LIU Qiao-ling , WEN Bang-chun. Failure sensitivity investigation in impact-rub malfunction for rotating machinery[J].Journal of Vibration Engineering,2007, 20(2):189-193.