劉繼華，陳立，羅載奇，李奕新，錢興良

（中國(guó)燃?xì)鉁u輪研究院，成都610500）

非線性彈性轉(zhuǎn)子-軸承系統(tǒng)穩(wěn)定裕度數(shù)值分析

劉繼華，陳立，羅載奇，李奕新，錢興良

（中國(guó)燃?xì)鉁u輪研究院，成都610500）

基于互補(bǔ)群群際能量壁壘準(zhǔn)則量化理論，提出了軌跡加速度-位移擴(kuò)展相平面穩(wěn)定裕度分析法，從理論分析、數(shù)值計(jì)算兩方面驗(yàn)證了該方法的正確性。建立了采用Capone圓軸承非線性油膜力模型的彈性轉(zhuǎn)子-軸承系統(tǒng)模型，并用數(shù)值積分和龐加萊映射方法，得到系統(tǒng)在某些參數(shù)域中的分岔圖、軸心軌跡圖、龐加萊映射圖、時(shí)間歷程圖和頻譜圖，得出分岔失穩(wěn)速度隨質(zhì)量偏心的變化規(guī)律：分岔速度隨質(zhì)量偏心的趨勢(shì)是先減小后增大，擬合曲線符合三次高斯公式，擬合精度達(dá)0.998 2。數(shù)值分析結(jié)果為該類轉(zhuǎn)子-軸承系統(tǒng)的設(shè)計(jì)和安全運(yùn)行提供了理論參考。

旋轉(zhuǎn)機(jī)械；轉(zhuǎn)子-軸承系統(tǒng)；加速度-位移擴(kuò)展相平面；非線性；穩(wěn)定裕度；分岔失穩(wěn)轉(zhuǎn)速；質(zhì)量偏心

1　引言

高轉(zhuǎn)速輕結(jié)構(gòu)是高速旋轉(zhuǎn)機(jī)械的設(shè)計(jì)趨勢(shì)，它提高了旋轉(zhuǎn)機(jī)械的性能，但也引發(fā)出很多嚴(yán)重的問題［1］。比如原來能用線性理論得出的結(jié)果不僅誤差大，而且無法解釋實(shí)際中出現(xiàn)的自激振動(dòng)、分岔、跳躍等現(xiàn)象［2］。動(dòng)力機(jī)械工程中，具有強(qiáng)烈非線性的重要組成部分是油膜軸承轉(zhuǎn)子系統(tǒng)［3-6］。

定性的穩(wěn)定性分析只能確定某些特殊點(diǎn)系統(tǒng)是否穩(wěn)定，而無法判斷一個(gè)原來穩(wěn)定的系統(tǒng)在參數(shù)小變化下是否還能保持穩(wěn)定，也無法估計(jì)一個(gè)不穩(wěn)定系統(tǒng)的參數(shù)如何變化、變化多少才能使其穩(wěn)定。這時(shí)需采用穩(wěn)定定量分析方法［7］。

薛禹勝院士在電力系統(tǒng)穩(wěn)定性分析中提出的互補(bǔ)群群際能量壁壘準(zhǔn)則量化理論，在國(guó)內(nèi)外電力系統(tǒng)中得到成功應(yīng)用［8］，后來又被引入到非線性轉(zhuǎn)子系統(tǒng)中［9-10］。文獻(xiàn)［11］采用動(dòng)能差序列衰減指數(shù)法來衡量工頻周期運(yùn)動(dòng)的穩(wěn)定裕度。由于分析的是系統(tǒng)瞬態(tài)響應(yīng)的數(shù)據(jù)，每次實(shí)驗(yàn)需用小重物敲擊轉(zhuǎn)子圓盤幾次。這樣不僅對(duì)轉(zhuǎn)子系統(tǒng)有損壞，而且每次敲擊的力度不同，所測(cè)的數(shù)據(jù)也會(huì)不同，從而造成一定的誤差。

本文提出了基于軌跡的加速度-位移擴(kuò)展相平面穩(wěn)定裕度分析法，并從數(shù)學(xué)理論和數(shù)值計(jì)算兩個(gè)方面求證其正確性。這種方法的好處是直接從穩(wěn)態(tài)數(shù)據(jù)著手，避免敲擊轉(zhuǎn)子-軸承系統(tǒng)而產(chǎn)生的系統(tǒng)損壞?；贑apone圓軸承非線性油膜力模型，對(duì)彈性轉(zhuǎn)子-軸承系統(tǒng)進(jìn)行數(shù)值計(jì)算，得出系統(tǒng)隨參數(shù)變化的穩(wěn)定裕度變化。采用分岔圖、軸心軌跡圖、龐加萊映射圖、時(shí)間歷程圖和頻譜圖分析轉(zhuǎn)子盤中心振動(dòng)響應(yīng)，得出分岔速度隨質(zhì)量偏心的變化規(guī)律。

2　數(shù)學(xué)模型

分析對(duì)象為單盤對(duì)稱轉(zhuǎn)子，見圖1，兩端支撐為滑動(dòng)短軸承。采用簡(jiǎn)單離散轉(zhuǎn)子模型［12］，將圖中的轉(zhuǎn)子模型離散為三個(gè)質(zhì)點(diǎn)，軸段處的質(zhì)量分別離散到軸頸和輪盤處。單盤轉(zhuǎn)子系統(tǒng)的離散模型如圖2所示，兩軸段等效為剛度k，輪盤處的阻尼為c1，非線性油膜力分別作用在兩個(gè)軸頸處的質(zhì)點(diǎn)上。

圖1　單盤對(duì)稱轉(zhuǎn)子模型Fig.1 Mode of single-disk rotor

圖2　單盤轉(zhuǎn)子系統(tǒng)的離散模型Fig.2 Discrete mode of single-disk rotor bearing

參數(shù)采用無量綱處理，無量綱時(shí)間τ=wt；無量綱位移x=X/c，y=Y/c，c為軸承半徑間隙（mm）。則無量綱速度，其中“?”表示d/dt，“′”表示d/dτ；無量綱偏心量ρ=e/c，e為偏心量（mm）；無量綱非線性油膜力分量［7］，；轉(zhuǎn)子受到的無量綱重力；無量綱質(zhì)量，其中δ為Sommerfeld修正數(shù)，反映潤(rùn)滑油粘度、軸承徑向間隙及長(zhǎng)徑比等多因素的影響，，R為軸承半徑（mm），L為軸承有效寬度（mm），η為潤(rùn)滑動(dòng)力粘度（Pa·s）?？傻贸鱿铝袩o量綱方程：

式中：fx、fy為無量綱非線性油膜力，其詳細(xì)表達(dá)式參見文獻(xiàn)［13］。

3　加速度-位移擴(kuò)展相平面穩(wěn)定裕度計(jì)算法

單自由度剛體系統(tǒng)在一個(gè)自由度上的運(yùn)動(dòng)方程為：

式（2）可變?yōu)椋?/p>

將式（3）兩邊分別對(duì)δ積分，得到：

式中：M為廣義慣量，Pm、Pe分別為廣義驅(qū)動(dòng)力和廣義制動(dòng)力，δ為剛體位移，w為剛體運(yùn)動(dòng)速度，為單剛體運(yùn)動(dòng)的等值加速度。式（5）的物理意義為：不平衡力對(duì)剛體所做的功轉(zhuǎn)換為剛體的動(dòng)能。剛體動(dòng)能的變化可以用外力-位移擴(kuò)展相平面（P-δ平面）上的面積來表達(dá)。

當(dāng)M為常數(shù)時(shí)：

將式（6）代入式（5）可得出：

對(duì)于M為常數(shù)的系統(tǒng)，剛體動(dòng)能變化也可用加速度-位移擴(kuò)展相平面（γ-δ平面）上的面積來表達(dá)。，由于x存在二階導(dǎo)數(shù)，一階導(dǎo)數(shù)必連續(xù)，，所以。這就從理論分析上說明了。

圖3為轉(zhuǎn)速w=300 rad/s的周期解的加速度-位移擴(kuò)展相平面，用MATLAB計(jì)算其閉合圖形面積，得出S正=0.001 2，S負(fù)=0.001 2。這就從數(shù)值方面說明加速度-位移擴(kuò)展相平面所圍成的面積代數(shù)和為零。

圖3　轉(zhuǎn)速300 rad/s時(shí)的加速度-位移擴(kuò)展相平面Fig.3 Acceleration-displacement expansion phase plane dynamics（w=300 rad/s）

從數(shù)學(xué)理論和數(shù)值計(jì)算得出，對(duì)于M為常數(shù)的系統(tǒng)，對(duì)于穩(wěn)態(tài)周期1解，系統(tǒng)每個(gè)周期T內(nèi)，在γ-δ平面上所圍成面積的代數(shù)和為零，即動(dòng)能增加面積與動(dòng)能減小面積相等。實(shí)際上，正面積的大小反應(yīng)了系統(tǒng)貯存能量的多少，所以可用加速度-位移擴(kuò)展相平面的正面積表示系統(tǒng)的穩(wěn)定性。

4　數(shù)值計(jì)算分析

式（1）描述了一個(gè)非線性動(dòng)力學(xué)系統(tǒng)的微分耦合方程，可以由數(shù)值計(jì)算得出近似解。本文使用MATLAB編程，由于油膜力和基座支撐力的強(qiáng)非線性作用，采用Rung-Kutta法求解。選用的時(shí)間步長(zhǎng)為2π/100，誤差小于10-6，計(jì)算1 500周期，取系統(tǒng)后500周期的響應(yīng)，以消除其瞬態(tài)響應(yīng)。系統(tǒng)穩(wěn)定性分析時(shí)，選取如下轉(zhuǎn)子模型參數(shù)［13］：m1=420 kg，m2= 50 kg，L=28.5 mm，c=0.2 mm，R=57 mm，η=0.018 Pa·s，k=2.105×108，c1=3×10 N·s/m。

4.1隨轉(zhuǎn)速增大的動(dòng)態(tài)響應(yīng)的分析

轉(zhuǎn)軸和葉輪的加工及安裝過程中，總會(huì)出現(xiàn)殘余偏心量，使得轉(zhuǎn)子系統(tǒng)在運(yùn)轉(zhuǎn)過程中產(chǎn)生振動(dòng)。本節(jié)重點(diǎn)分析油膜支撐轉(zhuǎn)子-軸承系統(tǒng)在偏心量-轉(zhuǎn)速參數(shù)域內(nèi)的穩(wěn)定性及其分岔規(guī)律。ρ=0.1、w= 300 rad/s時(shí)系統(tǒng)的動(dòng)態(tài)響應(yīng)如圖4所示。此時(shí)軸心軌跡為一個(gè)較規(guī)則的橢圓；Poincare圖只出現(xiàn)一個(gè)離散的點(diǎn)；頻譜圖上只出現(xiàn)了一個(gè)工頻；時(shí)間歷程圖為較規(guī)則的周期運(yùn)動(dòng)。由此可以判斷，系統(tǒng)此時(shí)的運(yùn)動(dòng)表現(xiàn)為工頻周期運(yùn)動(dòng)。

圖4　無量綱偏心量為0.1轉(zhuǎn)速為300 rad/s時(shí)的動(dòng)態(tài)響應(yīng)圖Fig.4 Dynamic response diagram（ρ=0.1，w=300 rad/s）

進(jìn)一步增大系統(tǒng)轉(zhuǎn)速，系統(tǒng)同頻運(yùn)動(dòng)周期發(fā)生分岔而失穩(wěn)。ρ=0.1、w=500 rad/s時(shí)系統(tǒng)的動(dòng)態(tài)響應(yīng)如圖5所示。軸心軌跡為不相交的內(nèi)凹封閉曲線；Poincare圖上表現(xiàn)為兩個(gè)離散的點(diǎn)；頻譜圖上出現(xiàn)一個(gè)基頻率和一個(gè)半頻率，且半頻率幅值比基頻率幅值大；時(shí)間歷程圖為較規(guī)則的倍周期運(yùn)動(dòng)。由此可以得出，系統(tǒng)此時(shí)的運(yùn)動(dòng)表現(xiàn)為倍周期運(yùn)動(dòng)。

當(dāng)轉(zhuǎn)速升高到w=650 rad/s時(shí)，系統(tǒng)動(dòng)態(tài)響應(yīng)如圖6所示。此時(shí)軸心軌跡為在一個(gè)圓環(huán)內(nèi)相交的曲線；Poincare圖為一些離散的點(diǎn)，但這些點(diǎn)形成一個(gè)封閉的形狀；頻譜圖上出現(xiàn)了一些不可約分的頻率。這些特點(diǎn)都說明，系統(tǒng)的運(yùn)動(dòng)表現(xiàn)為概周期運(yùn)動(dòng)。

4.2分岔圖對(duì)比分析

進(jìn)一步用數(shù)值計(jì)算方法得出系統(tǒng)在ρ=0.1時(shí)的全局分岔圖，見圖7?？梢姡瑆≤423 rad/s時(shí)都表現(xiàn)為周期運(yùn)動(dòng)，w＞423 rad/s時(shí)系統(tǒng)的同頻周期運(yùn)動(dòng)發(fā)生分岔而失穩(wěn)，424 rad/s＜w＜560 rad/s時(shí)為周期運(yùn)動(dòng)，570 rad/s＜w＜720 rad/s時(shí)系統(tǒng)由周期運(yùn)動(dòng)逐漸分岔為概周期運(yùn)動(dòng)。

圖5　無量綱偏心量為0.1轉(zhuǎn)速為500 rad/s時(shí)的動(dòng)態(tài)響應(yīng)圖Fig.5 Dynamic response diagram（ρ=0.1，w=500 rad/s）

圖6　無量綱偏心量為0.1轉(zhuǎn)速為650 rad/s時(shí)的動(dòng)態(tài)響應(yīng)圖Fig.6 Dynamic response diagram（ρ=0.1，w=650 rad/s）

圖7　無量綱偏心量為0.1時(shí)的全局分岔圖Fig.7 Global bifurcation diagram（ρ=0.1）

圖8示出了ρ=0.4時(shí)的分岔細(xì)化圖?？梢?，轉(zhuǎn)速358 rad/s與360 rad/s之間只相差2 rad/s，但幅值卻相差0.5個(gè)無量綱值，此時(shí)系統(tǒng)也沒出現(xiàn)分岔現(xiàn)象。這種由外力強(qiáng)迫頻率發(fā)生微小變化而振幅出現(xiàn)劇烈變化的現(xiàn)象，稱為跳躍現(xiàn)象，為非線性所特有。

圖8　無量綱偏心量為0.4時(shí)的分岔細(xì)化圖Fig.8 Local bifurcation diagram（ρ=0.4）

由于不同的偏心量系統(tǒng)將出現(xiàn)不同的失穩(wěn)特性，因此有必要研究偏心量變化對(duì)系統(tǒng)穩(wěn)定的影響。從圖9可知，隨著偏心量的增大，開始時(shí)失穩(wěn)轉(zhuǎn)速不斷降低；當(dāng)無量綱偏心量增至0.25時(shí)，系統(tǒng)失穩(wěn)轉(zhuǎn)速又逐漸升高。較小和較大的不平衡量都會(huì)增大系統(tǒng)的穩(wěn)定性，但較大的不平衡量會(huì)引起轉(zhuǎn)子軸承受力增大，且有可能出現(xiàn)跳躍現(xiàn)象。因此，為確保系統(tǒng)安全運(yùn)行，必須對(duì)轉(zhuǎn)子進(jìn)行較高精度的動(dòng)平衡。

圖9　轉(zhuǎn)子軸承系統(tǒng)偏心量-轉(zhuǎn)速參數(shù)域分岔集Fig.9 Bifurcation set of rotor bearing system with the parameter region of eccentricity and speed

對(duì)轉(zhuǎn)子軸承系統(tǒng)偏心量-轉(zhuǎn)速參數(shù)域分岔集進(jìn)行擬合得出，三次高斯公式有很高的擬合精度，擬合度達(dá)0.998 2。失穩(wěn)轉(zhuǎn)速與偏心率、潤(rùn)滑粘度、軸承間隙、軸承寬徑比等參數(shù)，形成一個(gè)三次高斯公式的表達(dá)式。

5　系統(tǒng)穩(wěn)定性裕度的計(jì)算分析

5.1周期運(yùn)動(dòng)穩(wěn)定裕度

周期運(yùn)動(dòng)的穩(wěn)定裕度用加速度-位移擴(kuò)展相平面上相比正面積準(zhǔn)則，來衡量周期運(yùn)動(dòng)的穩(wěn)定情況。首先畫出某個(gè)周期運(yùn)動(dòng)的加速度-位移相平面圖，計(jì)算出正面積；然后找出第一個(gè)周期分岔點(diǎn)，計(jì)算出臨界正面積Acr；最后得出周期運(yùn)動(dòng)穩(wěn)定裕度β。β值越小，系統(tǒng)的穩(wěn)定性程度越差；臨界狀態(tài)時(shí)，β=0。

5.2倍周期運(yùn)動(dòng)穩(wěn)定裕度

隨著參數(shù)的變化，系統(tǒng)經(jīng)倍周期分岔，其運(yùn)動(dòng)狀態(tài)由周期運(yùn)動(dòng)變?yōu)楸吨芷谶\(yùn)動(dòng)。倍周期運(yùn)動(dòng)的明顯特征是穩(wěn)態(tài)時(shí)動(dòng)能差序列為常數(shù)C1，該數(shù)值為一個(gè)周期內(nèi)最大擺次中的最大動(dòng)能與同向擺動(dòng)的另一個(gè)擺次中最大動(dòng)能的差值的絕對(duì)值，反映了倍周期運(yùn)動(dòng)中低頻分量的大小及倍周期運(yùn)動(dòng)的劇烈程度。故把穩(wěn)態(tài)時(shí)動(dòng)能差序列之值加上負(fù)號(hào)，定義為倍周期運(yùn)動(dòng)的穩(wěn)定裕度，即有：

C1越大，表示倍周期運(yùn)動(dòng)中的低頻分量越大，倍周期運(yùn)動(dòng)越劇烈，系統(tǒng)的失穩(wěn)程度越嚴(yán)重。臨界狀態(tài)時(shí)，β=0。

5.3概周期運(yùn)動(dòng)穩(wěn)定裕度

鑒于概周期運(yùn)動(dòng)動(dòng)能差序列的非周期特征，研究時(shí)必須選擇合適的樣本長(zhǎng)度。取n個(gè)近似周期（本研究取10個(gè)近似周期），在動(dòng)能差序列上依次截取形成一個(gè)樣本序列。計(jì)算對(duì)應(yīng)于每一個(gè)長(zhǎng)度內(nèi)動(dòng)能差序列的數(shù)學(xué)期望（Ex），從而形成一個(gè)均值序列，該值反映了概周期運(yùn)動(dòng)動(dòng)態(tài)中心點(diǎn)動(dòng)能差的平均變化程度。因此，概周期運(yùn)動(dòng)的穩(wěn)定裕度可以表示為：

β值越小，說明系統(tǒng)失穩(wěn)程度越嚴(yán)重，概周期運(yùn)動(dòng)越劇烈。臨界狀態(tài)時(shí)，β=0。

5.4加速度-位移擴(kuò)展相平面穩(wěn)定裕度計(jì)算法求解

流程

通過上述分析，得出了不同運(yùn)動(dòng)的穩(wěn)定裕度算法，按照以下流程，可以得出加速度-位移擴(kuò)展相平面穩(wěn)定裕度，從而確定系統(tǒng)的穩(wěn)定裕度。

（1）做出系統(tǒng)的位移-加速擴(kuò)展平面圖，通過圖上動(dòng)態(tài)中心點(diǎn)的個(gè)數(shù)來判斷系統(tǒng)所處運(yùn)動(dòng)狀態(tài)。

（2）找出周期分岔點(diǎn)，求出其加速度-位移擴(kuò)展正面積A+，并令其為臨界面積。

（4）系統(tǒng)處于倍周期運(yùn)動(dòng)時(shí)，做出其速度序列圖，求出速度峰值。

（5）系統(tǒng)處于概周期運(yùn)動(dòng)時(shí)，做出其速度序列圖，取n個(gè)近似周期，求出速度峰值v1、v2…，動(dòng)能差序列；然后求出Si的數(shù)學(xué)期望E（Si），β=-E（Si）。

5.5穩(wěn)定裕度計(jì)算算例

采用加速度-位移擴(kuò)展相平面穩(wěn)定裕度計(jì)算法，求解式（1）中ρ=0.1時(shí)系統(tǒng)的穩(wěn)定裕度變化規(guī)律。周期運(yùn)動(dòng)的穩(wěn)定裕度變化規(guī)律如圖10所示，可見系統(tǒng)周期運(yùn)動(dòng)時(shí)，隨著轉(zhuǎn)速的增大穩(wěn)定裕度不斷降低，具有很好的單調(diào)性，符合穩(wěn)定性的概念。

圖10　無量綱偏心量為0.1時(shí)周期運(yùn)動(dòng)的穩(wěn)定裕度規(guī)律圖Fig.10 Periodic motion stability margin（ρ=0.1）

圖11　無量綱偏心量為0.1時(shí)倍周期運(yùn)動(dòng)的穩(wěn)定裕度規(guī)律圖Fig.11 Double-periodic motion stability margin（ρ=0.1）

系統(tǒng)倍周期穩(wěn)定裕度變化規(guī)律如圖11所示?？梢?，轉(zhuǎn)速為423 rad/s時(shí)，穩(wěn)定裕度為-0.000 196 7，可近似為零。分岔圖上的分岔失穩(wěn)轉(zhuǎn)速為423 rad/s，分岔點(diǎn)的穩(wěn)定裕度為零。隨著轉(zhuǎn)速的增大，穩(wěn)定裕度逐漸減小，當(dāng)轉(zhuǎn)速增大到480 rad/s時(shí)，穩(wěn)定裕度的單調(diào)性開始出現(xiàn)轉(zhuǎn)折；隨著轉(zhuǎn)速的繼續(xù)增大，穩(wěn)定裕度開始逐漸增大，當(dāng)轉(zhuǎn)速增大到560 rad/s時(shí)，穩(wěn)定裕度逐漸趨于零。由前文分岔圖可知，此時(shí)處在分岔點(diǎn)，重新回到工頻周期運(yùn)動(dòng)。

系統(tǒng)概周期的穩(wěn)定裕度變化規(guī)律如圖12所示。可見，當(dāng)轉(zhuǎn)速低于580 rad/s時(shí)，系統(tǒng)的穩(wěn)定裕度為零，結(jié)合分岔圖可得出此時(shí)系統(tǒng)表現(xiàn)為工頻周期運(yùn)動(dòng)。當(dāng)轉(zhuǎn)速超過580 rad/s時(shí)，系統(tǒng)開始逐漸變?yōu)楦胖芷谶\(yùn)動(dòng)，穩(wěn)定裕度逐漸下降。當(dāng)轉(zhuǎn)速增大到660 rad/s時(shí)，穩(wěn)定裕度的單調(diào)性開始出現(xiàn)轉(zhuǎn)折，隨著轉(zhuǎn)速的增大穩(wěn)定裕度逐漸增大。

圖12　無量綱偏心量為0.1時(shí)概周期運(yùn)動(dòng)的穩(wěn)定裕度規(guī)律圖Fig.12 Almost-periodic motion stability margin（ρ=0.1）

從以上分析可得出，基于軌跡的加速度-位移擴(kuò)展相平面得出的穩(wěn)定裕度為零的點(diǎn)，與通過Poincare映射法做出的分岔圖上的分岔點(diǎn)，其對(duì)應(yīng)轉(zhuǎn)速基本相同，說明了這種方法的正確性。基于軌跡的加速度-位移擴(kuò)展相平面方法不僅可以判斷系統(tǒng)的分岔點(diǎn)，還具有其獨(dú)特的優(yōu)點(diǎn)。通過這種方法得出的穩(wěn)定裕度可用來判斷系統(tǒng)的穩(wěn)定度和失穩(wěn)度，適當(dāng)改變系統(tǒng)的質(zhì)量偏心、潤(rùn)滑油粘度、軸承間隙等參數(shù)可增加系統(tǒng)的穩(wěn)定性；在保證系統(tǒng)穩(wěn)定的條件下，增加系統(tǒng)的穩(wěn)定裕度，使轉(zhuǎn)子-軸承系統(tǒng)處于最佳的穩(wěn)定性和安全性運(yùn)行狀態(tài)。

6　結(jié)論

（1）提出了軌跡加速度-位移擴(kuò)展相平面穩(wěn)定裕度分析法，從理論分析、數(shù)值計(jì)算兩方面驗(yàn)證了該方法的正確性。

（2）偏心量較大的系統(tǒng)的失穩(wěn)分岔速度較大，但偏心量較大其轉(zhuǎn)子系統(tǒng)會(huì)存在跳躍現(xiàn)象的可能性，軸承受力增加；需對(duì)轉(zhuǎn)子系統(tǒng)進(jìn)行較高精度的動(dòng)平衡，以保證系統(tǒng)的穩(wěn)定性和安全性。

（3）對(duì)于穩(wěn)態(tài)周期運(yùn)動(dòng)，系統(tǒng)每個(gè)周期內(nèi)在加速度-位移擴(kuò)展相平面上所圍成面積的代數(shù)和為零。

（4）通過基于軌跡的加速度-位移擴(kuò)展相平面得出的穩(wěn)定裕度為零的點(diǎn)，與通過Poincare映射法做出的分岔圖上的分岔點(diǎn)，其對(duì)應(yīng)轉(zhuǎn)速基本相同，說明了基于軌跡的加速度-位移擴(kuò)展相平面方法的正確性。

（5）基于軌跡的加速度-位移擴(kuò)展相平面方法，不僅可以判斷系統(tǒng)的分岔點(diǎn)，還可以得出穩(wěn)定裕度，從而判斷系統(tǒng)的穩(wěn)定程度。

［1］聞邦椿，顧佳柳，夏松波，等.高等轉(zhuǎn)子動(dòng)力學(xué)——理論、技術(shù)與應(yīng)用［M］.北京：機(jī)械工業(yè)出版社，2000.

［2］劉延柱，陳立群.非線性振動(dòng)［M］.北京：高等教育出版社，2001.

［3］Muszynska A.Partial lateral rotor to stator rubs［C］//.Proceedings of the 3rdInternational Conference on Vibration in Rotating Machinery.York UK：1984.

［4］Helio Fiori de C，Katia Lucchesi C，Rainer N.Whirl and whip instabilities in rotor-bearing system considering a nonlinear force model［J］.Journal of Sound and Vibration，2008，371（1-2）：273—293.

［5］Yau H T，Chen C K，Chen C L.Chaos and bifurcation analysis of a flexible rotor supported by short journal bearings with non-linear suspension［J］.Journal of Mechanical Engineering Science，2000，214（7）：931—947.

［6］張新江，武新華.彈性轉(zhuǎn)子-軸承-基礎(chǔ)系統(tǒng)的非線性振動(dòng)研究［J］.振動(dòng)工程學(xué)報(bào)，2001，14（2）：228—232.

［7］黃文虎，夏松波.旋轉(zhuǎn)機(jī)械非線性動(dòng)力學(xué)設(shè)計(jì)基礎(chǔ)理論與方法［M］.北京：科學(xué)出版社，2006.

［8］薛禹勝.運(yùn)動(dòng)穩(wěn)定性量化理論［M］.南京：江蘇科學(xué)技術(shù)出版社，1999.

［9］Xue Y S.The stability-preserving trajectory-reduction methodology for analyzing nonlinear dynamics［C］//.KeynoteinInternationalConferenceonInfo-tech&Info-net（IEEE Catalog Number 01EX479）.Beijing：2001.

［10］李銀山，陳予恕，薛禹勝.非線性油膜力軸承上不平衡彈性轉(zhuǎn)子的穩(wěn)定裕度［J］.機(jī)械工程學(xué)報(bào)，2002，38（9）：27—32.

［11］陸永杰.滑動(dòng)軸承轉(zhuǎn)子系統(tǒng)穩(wěn)定性量化分析方法的試驗(yàn)研究［D］.石家莊：河北科技大學(xué)，2009.

［12］聞邦椿，武新華.故障旋轉(zhuǎn)機(jī)械非線性動(dòng)力學(xué)的理論與實(shí)驗(yàn)［M］.北京：科學(xué)出版社，2004.

［13］Adiletta G，Guido A.Chaotic motions of a rigid rotor in short journal bearings［J］.Journal of Nonlinear Dynamics，1996，（10）：251—269.

Numerical analysis of stability margin for nonlinear vibration elastic rotor-bearing system

LIU Ji-hua，CHEN Li，LUO Zai-qi，LI Yi-xin，QIAN Xing-liang
（China Gas Turbine Establishment，Chengdu 610500，China）

A stability margin analysis method for the acceleration displacement of trajectory extension plane was proposed based on complementary-cluster energy-barrier criterion（CCEBC）of the stability of quantification theory.The correctness of this method is verified through theoretical analysis and numerical simulation.An elastic rotor bearing system was investigated based on the assumption of Capone bearing mode.The bifurcation diagrams，the shaft centerline orbit，Poincare maps，time histories and frequency spectrums of the system in some typical parameter regions were acquired by using the Poincare maps theory and numerical integration method.The results show that bifurcation speed varies as mass eccentricity changes.The speed decreased first and then increased，curve fitting with three times the Gauss formula and the fitting accuracy reached 0.998 2.The analysis results provide theoretical reference for designing and safely operating of this kind of systems.

rotating machine；rotor-bearing system；acceleration-displacement expansion phase plane；nonlinear；stability margin；bifurcation speed；mass eccentricity

O322；V231.96

A

1672-2620（2015）06-0039-06

2015-03-25；

2015-09-12

劉繼華（1988-），男，貴州榕江縣人，碩士，主要從事轉(zhuǎn)子動(dòng)力學(xué)、航空發(fā)動(dòng)機(jī)起動(dòng)研究。