李自剛　李明　江俊

摘要： 考慮轉(zhuǎn)子交角不對(duì)中和質(zhì)量不平衡等因素，研究了在滑動(dòng)軸承支承下柔性轉(zhuǎn)子-軸承耦合系統(tǒng)的非線性動(dòng)力學(xué)行為。首先，基于轉(zhuǎn)子間交角不對(duì)中的約束關(guān)系，利用第二類Lagrange方程推導(dǎo)了具有交角不對(duì)中故障的柔性多轉(zhuǎn)子系統(tǒng)運(yùn)動(dòng)微分方程。采用數(shù)值方法，分析系統(tǒng)的非線性振動(dòng)特性，例如，系統(tǒng)的軸心軌跡、響應(yīng)頻譜和最大Lyapunov 指數(shù)等。結(jié)果表明：在較低轉(zhuǎn)速時(shí)，系統(tǒng)主要呈現(xiàn)出與轉(zhuǎn)速同步的周期運(yùn)動(dòng)特性。隨著轉(zhuǎn)速的提高，穩(wěn)態(tài)響應(yīng)在某些參數(shù)下出現(xiàn)分叉、跳躍以及混沌等非線性現(xiàn)象。最后討論了交角不對(duì)中量以及質(zhì)量不平衡對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響。

關(guān)鍵詞： 非線性動(dòng)力學(xué); 轉(zhuǎn)子-軸承系統(tǒng); 交角不對(duì)中; 分叉與混沌; 質(zhì)量不平衡

中圖分類號(hào)： O322; O347.6; TH113.1 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)： 1004-4523（2019）03-0509-08

DOI：10.16385/j.cnki.issn.1004-4523.2019.03.016

引 言

轉(zhuǎn)子-軸承系統(tǒng)是旋轉(zhuǎn)機(jī)械的核心組成部分。在實(shí)際系統(tǒng)中，由于加工精度、安裝誤差等原因，相鄰轉(zhuǎn)子的軸線往往存在不對(duì)中現(xiàn)象，即使對(duì)中良好的轉(zhuǎn)子也會(huì)在其頻繁啟停過程中產(chǎn)生位置偏離。轉(zhuǎn)子不對(duì)中能夠引起機(jī)組的劇烈振動(dòng)和異常響聲，已成為機(jī)械故障最為重要的原因之一，如在大型發(fā)電機(jī)組中，不對(duì)中常常引起軸的撓曲變形、油膜失穩(wěn)以及轉(zhuǎn)子與定子間的碰磨等故障，危害很大[1-4]。文獻(xiàn)[2]通過對(duì)實(shí)際機(jī)組的故障診斷，指出引起系統(tǒng)振動(dòng)超標(biāo)的根本原因是發(fā)電機(jī)-勵(lì)磁機(jī)聯(lián)軸器不對(duì)中。在文獻(xiàn)[3-4]中，作者以實(shí)際發(fā)電機(jī)轉(zhuǎn)子系統(tǒng)為研究對(duì)象，分別研究了具有平行不對(duì)中和軸承不對(duì)中故障的轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)行為，揭示了該類系統(tǒng)的一些非線性振動(dòng)特征。然而，工程中往往也存在轉(zhuǎn)子間的幾何關(guān)系為交角不對(duì)中的故障形式。因此研究在滑動(dòng)軸承支承下具有交角不對(duì)中約束的多轉(zhuǎn)子-軸承耦合系統(tǒng)的非線性振動(dòng)特性及其參數(shù)影響對(duì)系統(tǒng)的動(dòng)力學(xué)設(shè)計(jì)、安裝以及運(yùn)行中的故障診斷具有重要意義。

關(guān)于交角不對(duì)中的研究最早可以追述到上個(gè)世紀(jì)70年代，但對(duì)此引起廣泛關(guān)注的大約始于80，90年代。近10年來，隨著國民經(jīng)濟(jì)和國防工業(yè)的快速發(fā)展，轉(zhuǎn)子系統(tǒng)趨于大型、高速、重載，而轉(zhuǎn)子與定子之間的間隙卻在相對(duì)縮小，轉(zhuǎn)子不對(duì)中所引起的振動(dòng)也就變得更加嚴(yán)重，這要求轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)設(shè)計(jì)更為精確合理。文獻(xiàn)[5]利用有限元法分析了一個(gè)簡化的不對(duì)中轉(zhuǎn)子模型，并指出二倍頻分量是這類轉(zhuǎn)子系統(tǒng)的特征頻率;文獻(xiàn)[6]中，作者通過分析聯(lián)軸節(jié)的變形和受力，研究了交角不對(duì)中軸系的扭振特征。文獻(xiàn)[7-8]利用數(shù)值方法分別研究了由球鉸聯(lián)軸器和柔性聯(lián)軸器連接的兩個(gè)剛性轉(zhuǎn)子系統(tǒng)的振動(dòng)特性，并分析了角度不對(duì)中量和軸承剛度等因素對(duì)系統(tǒng)的動(dòng)力學(xué)特性的影響，取得了一些重要的研究成果;文獻(xiàn)[9]則采用Lagrange方法分析了2個(gè)支承在線性化滑動(dòng)軸承上的剛性Jeffcott 轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性，其中不對(duì)中效應(yīng)主要體現(xiàn)在系統(tǒng)的彈性力上，而在系統(tǒng)的動(dòng)態(tài)力中未予以考慮。近年來，一些學(xué)者則更加關(guān)注于振動(dòng)信號(hào)中的諧波響應(yīng)。文獻(xiàn)[10]通過分析軸心軌跡和頻譜圖的特點(diǎn)，揭示了不對(duì)中轉(zhuǎn)子系統(tǒng)的一些非線性動(dòng)力學(xué)行為;文獻(xiàn)[11-12]則考慮了2個(gè)轉(zhuǎn)子間的連接剛度，建立了剛性多轉(zhuǎn)子系統(tǒng)的動(dòng)態(tài)不對(duì)中模型，研究表明當(dāng)轉(zhuǎn)角不對(duì)中時(shí)，系統(tǒng)的振動(dòng)主要以工頻為主，同時(shí)存在著組合頻率成分;文獻(xiàn)[13]從實(shí)驗(yàn)角度對(duì)雙跨轉(zhuǎn)子的不對(duì)中和不平衡耦合問題進(jìn)行了試驗(yàn)研究。文獻(xiàn)[14]采用隨機(jī)理論分析了在非線性支承下剛性轉(zhuǎn)子的角不對(duì)中問題。

以上關(guān)于交角不對(duì)中的研究中，所分析的機(jī)械結(jié)構(gòu)大多表現(xiàn)為轉(zhuǎn)子與軸承同時(shí)具有不對(duì)中，或允許較大角補(bǔ)償?shù)娜嵝月?lián)軸器-轉(zhuǎn)子系統(tǒng)。例如，文獻(xiàn)[15]建立了含萬向鉸偏斜旋轉(zhuǎn)軸的橫向振動(dòng)模型，該結(jié)構(gòu)允許較大的角不對(duì)中量。而對(duì)于軸承對(duì)中良好，但轉(zhuǎn)子交角不對(duì)中的情況研究還比較少（軸承不對(duì)中便于調(diào)節(jié)，轉(zhuǎn)子不對(duì)中則不便）。由于滑動(dòng)軸承間隙限制和長度效應(yīng)的影響，這種不對(duì)中情況不能允許轉(zhuǎn)子間有過大的角位移，但其對(duì)整個(gè)系統(tǒng)穩(wěn)定性產(chǎn)生的影響卻是不能忽視的。本文考慮了質(zhì)量不平衡等因素后，重點(diǎn)探討了轉(zhuǎn)子間具有交角不對(duì)中約束的柔性轉(zhuǎn)子-軸承耦合系統(tǒng)非線性動(dòng)力學(xué)行為。并深入分析了在不同轉(zhuǎn)速下交角不對(duì)中量以及質(zhì)量不平衡對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響，為該類轉(zhuǎn)子的故障檢測(cè)和在線監(jiān)測(cè)提供理論依據(jù)。

1 具有交角不對(duì)中故障的轉(zhuǎn)子-軸承系統(tǒng)動(dòng)力學(xué)模型

1.1 轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)模型

圖1為具有交角不對(duì)中故障的轉(zhuǎn)子-軸承系統(tǒng)示意圖。其中兩轉(zhuǎn)子間通過剛性聯(lián)軸器相連，并且在運(yùn)動(dòng)過程中保持大小不變的角不對(duì)中量α;m，M分別為軸頸和轉(zhuǎn)子本體的集中質(zhì)量;xi，yi為集中質(zhì)量點(diǎn)在坐標(biāo)系oxy下的坐標(biāo);k為轉(zhuǎn)軸剛度;a為圓盤的質(zhì)量偏心;Ω為轉(zhuǎn)子的轉(zhuǎn)速;t為時(shí)間。為方便問題的討論，現(xiàn)作如下基本假設(shè)：（1） 轉(zhuǎn)子間交角不對(duì)中量為微量;（2） 聯(lián)軸器兩端的相鄰軸承之間距離相對(duì)較小，不發(fā)生形變;（3） 2圓盤在轉(zhuǎn)軸的中間且近似地認(rèn)為作平面渦動(dòng)。

2.2 轉(zhuǎn)速對(duì)交角不對(duì)中轉(zhuǎn)子系統(tǒng)的影響

圖3為采用4-5階變步長Runge-Kutta法經(jīng)過充分衰減后求得的轉(zhuǎn)子系統(tǒng)在不對(duì)中軸承處的穩(wěn)態(tài)響應(yīng)隨轉(zhuǎn)速變化的分叉圖。從圖中能夠看出，對(duì)于對(duì)中情況（α=0），系統(tǒng)響應(yīng)在整個(gè)轉(zhuǎn)速區(qū)間內(nèi)始終保持周期運(yùn)動(dòng)。于此對(duì)比，微小不對(duì)中時(shí)（α=0.001 rad），系統(tǒng)在亞臨界和超臨界轉(zhuǎn)速區(qū)間內(nèi)呈現(xiàn)出一系列的分叉現(xiàn)象和混沌振蕩。在無量綱轉(zhuǎn)速ω≈2.45時(shí)系統(tǒng)通過倍周期分叉進(jìn)入混沌運(yùn)動(dòng)狀態(tài)，此時(shí)轉(zhuǎn)子的橫向振動(dòng)幅值較大;隨著轉(zhuǎn)速的提高，在線性共振區(qū)附近系統(tǒng)的混沌運(yùn)動(dòng)消失并重新回到周期軌道;而在ω≈4.40時(shí)轉(zhuǎn)子的穩(wěn)態(tài)響應(yīng)出現(xiàn)準(zhǔn)周期分叉，其中，當(dāng)ω≈4.80時(shí)，系統(tǒng)的振幅迅速增大，直至最后碰到軸承內(nèi)圈。圖4為其對(duì)應(yīng)的最大Lyapunov指數(shù)曲線。由圖4可見，轉(zhuǎn)子響應(yīng)分叉圖與最大Lyapunov指數(shù)曲線圖在周期運(yùn)動(dòng)、準(zhǔn)周期運(yùn)動(dòng)和混沌運(yùn)動(dòng)的描述上一一對(duì)應(yīng)，但在轉(zhuǎn)子響應(yīng)分叉圖難以區(qū)分開的準(zhǔn)周期運(yùn)動(dòng)和混沌運(yùn)動(dòng)，在最大Lyapunov指數(shù)曲線圖上則可以清楚區(qū)分開。

圖5為在不同轉(zhuǎn)速時(shí)交角不對(duì)中對(duì)轉(zhuǎn)子系統(tǒng)響應(yīng)頻率和軌跡的影響。由圖可見，在圖中的各轉(zhuǎn)速下對(duì)中轉(zhuǎn)子均以工頻振動(dòng)為主，呈現(xiàn)出周期運(yùn)動(dòng)特征，其軸心軌跡為一橢圓形，且幅值較小。對(duì)比來看，當(dāng)不對(duì)中量α=0.001 rad時(shí)，在亞臨界轉(zhuǎn)速ω=2.0下其響應(yīng)的Poincaré截面上呈現(xiàn)出1條封閉的曲線，表明系統(tǒng)在該轉(zhuǎn)速下作準(zhǔn)周期運(yùn)動(dòng)。從軌跡圖上可見，系統(tǒng)的穩(wěn)態(tài)響應(yīng)在1個(gè)環(huán)域內(nèi)作較為復(fù)雜的準(zhǔn)周期振蕩，頻域上主要表現(xiàn)為工頻及其倍頻分量;另外在0.423倍工頻等處有明顯的離散譜峰，這是由于不對(duì)中導(dǎo)致滑動(dòng)軸承的油膜壓力分配改變，從而可能誘發(fā)油膜渦動(dòng)[1]。當(dāng)ω=2.5時(shí)，轉(zhuǎn)子的軸心軌跡較為凌亂，頻域內(nèi)出現(xiàn)連續(xù)譜成分，Poincaré截面上呈現(xiàn)出無規(guī)則的點(diǎn)集，并且在該轉(zhuǎn)速下最大Lyapunov指數(shù)為0.0060，表明系統(tǒng)已進(jìn)入混沌運(yùn)動(dòng)。在臨界轉(zhuǎn)速附近ω=4.2時(shí)，轉(zhuǎn)子系統(tǒng)的響應(yīng)雖返回為周期1運(yùn)動(dòng)，但其軸心軌跡為明顯的“月牙形”，且在某些位置曲率半徑較大，頻譜圖中2倍頻成分突出，這些正是交角不對(duì)中故障具有的典型特征[5，10]。并且此時(shí)最大Lyapunov指數(shù)為-0.0059。當(dāng)轉(zhuǎn)子處于超臨界ω=4.7時(shí)，系統(tǒng)經(jīng)歷分叉再次進(jìn)入非周期運(yùn)動(dòng)。頻譜中相互不可公約的諧波分量出現(xiàn)，在該轉(zhuǎn)速下系統(tǒng)響應(yīng)的最大Lyapunov指數(shù)趨近于0，這些特征充分說明此時(shí)系統(tǒng)處于準(zhǔn)周期狀態(tài)。

2.3 不對(duì)中量對(duì)交角不對(duì)中轉(zhuǎn)子系統(tǒng)的影響

圖6為交角不對(duì)中轉(zhuǎn)子穩(wěn)態(tài)響應(yīng)在不同轉(zhuǎn)速下隨不對(duì)中量變化的分叉圖。從圖中可以看出，在較低轉(zhuǎn)速ω=2.0時(shí)，對(duì)于微小角不對(duì)中量，系統(tǒng)能夠保持在周期1運(yùn)動(dòng)。但隨著不對(duì)中量的逐漸增大，在α=0.0006 rad附近，轉(zhuǎn)子的穩(wěn)態(tài)響應(yīng)開始出現(xiàn)準(zhǔn)周期運(yùn)動(dòng)，其渦動(dòng)軌跡明顯增大。而在超臨界轉(zhuǎn)速ω=4.5下，隨著不對(duì)中量的增大，系統(tǒng)的動(dòng)力學(xué)行為表現(xiàn)的非常復(fù)雜，在α=0.0008 rad時(shí)直接從周期運(yùn)動(dòng)過渡到混沌振蕩。

為了進(jìn)一步研究不對(duì)中量對(duì)系統(tǒng)非線性振動(dòng)特性的影響，分析了在不同交角不對(duì)中量時(shí)轉(zhuǎn)子系統(tǒng)的穩(wěn)態(tài)響應(yīng)振幅隨轉(zhuǎn)速變化曲線。圖7中當(dāng)交角不對(duì)中量α=0時(shí)，即為對(duì)中情況下，系統(tǒng)的運(yùn)動(dòng)較為平穩(wěn)，其振動(dòng)幅度較小。然而不對(duì)中量的增大使轉(zhuǎn)子的穩(wěn)態(tài)響應(yīng)趨于復(fù)雜，出現(xiàn)明顯的跳躍現(xiàn)象。對(duì)比圖3可以看出，這種幅值跳躍現(xiàn)象是由于交角不對(duì)中故障而引起的轉(zhuǎn)子系統(tǒng)定性行為改變（解的數(shù)量變化）。而對(duì)于較大不對(duì)中量α=0.0010 rad時(shí)，隨著轉(zhuǎn)速的增加轉(zhuǎn)子的振幅迅速增大，最終碰到軸承內(nèi)圈。

2.4 質(zhì)量偏心對(duì)轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性的影響

下面主要討論質(zhì)量偏心變化對(duì)系統(tǒng)響應(yīng)的影響情況。圖8（a）和（b）分別為當(dāng)無量綱轉(zhuǎn)速ω=2.00和ω=4.50時(shí)，交角不對(duì)中轉(zhuǎn)子系統(tǒng)的穩(wěn)態(tài)響應(yīng)隨質(zhì)量偏心變化的分叉圖。從圖8（a）可以看出，隨著質(zhì)量偏心的增加，作用在轉(zhuǎn)子上的不平衡激勵(lì)增大，轉(zhuǎn)子的穩(wěn)態(tài)振動(dòng)中出現(xiàn)明顯的倍周期分叉現(xiàn)象。在不平衡量A≈0.10時(shí)，2倍周期中的每一個(gè)分支分別出現(xiàn)準(zhǔn)周期分叉，形成“肺”狀的兩片孤立區(qū)域，此時(shí)轉(zhuǎn)子的振動(dòng)幅度較大。從圖8（b）中可見，轉(zhuǎn)子的質(zhì)量不平衡能夠引起系統(tǒng)的周期分叉，在A=0.125時(shí)其穩(wěn)態(tài)響應(yīng)經(jīng)歷短暫的準(zhǔn)周期振蕩后，出現(xiàn)周期6運(yùn)動(dòng)，這預(yù)示著系統(tǒng)的運(yùn)動(dòng)將變得非常復(fù)雜。

圖9為在具有微小不對(duì)中量α=0.0005 rad時(shí)，不同質(zhì)量偏心下軸頸中心橫向振動(dòng)的幅頻特性曲線。經(jīng)過對(duì)比能夠看出，質(zhì)量偏心是影響系統(tǒng)振動(dòng)特性的重要因素。在低轉(zhuǎn)速時(shí)，較大的偏心量能夠使得轉(zhuǎn)子穩(wěn)態(tài)響應(yīng)出現(xiàn)明顯的幅值跳躍等非線性振動(dòng)現(xiàn)象。然而值得注意的是，在較高轉(zhuǎn)速時(shí)減小質(zhì)量不平衡量使得靠近聯(lián)軸器的軸承處振動(dòng)狀態(tài)并未改善，反而增大了轉(zhuǎn)子的振動(dòng)幅度。文獻(xiàn)[18]中通過實(shí)際機(jī)組的測(cè)試數(shù)據(jù)同樣反映了這種趨勢(shì)。

3 結(jié) 論

考慮了質(zhì)量不平衡等因素后，重點(diǎn)討論了轉(zhuǎn)子間具有交角不對(duì)中約束的柔性轉(zhuǎn)子-軸承耦合系統(tǒng)非線性動(dòng)力學(xué)行為。首先，根據(jù)描述交角不對(duì)中的約束關(guān)系，利用第二類Lagrange方程建立了具有交角不對(duì)中故障的柔性多轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)模型。理論分析表明，交角不對(duì)中轉(zhuǎn)子系統(tǒng)是1個(gè)具有強(qiáng)非線性特征的10自由度非自治系統(tǒng)。采用數(shù)值方法重點(diǎn)分析了不對(duì)中轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)行為，并深入討論了在不同轉(zhuǎn)速下交角不對(duì)中量和質(zhì)量不平衡對(duì)系統(tǒng)動(dòng)力學(xué)特性的影響，主要結(jié)論如下：

1）嵌入了非定常完整約束的轉(zhuǎn)子系統(tǒng)使得兩轉(zhuǎn)子間的位移通過角不對(duì)中量相互耦合，其耦合程度與轉(zhuǎn)子轉(zhuǎn)速密切相關(guān)。同時(shí)由于油膜間隙和長度效應(yīng)的影響，該約束影響軸頸內(nèi)油膜壓力重新分配，使得非線性油膜力與角不對(duì)中量α以及sin（Ωt）， cos（Ωt）相關(guān)，進(jìn)而引起系統(tǒng)的復(fù)雜運(yùn)動(dòng)。

2）在較低轉(zhuǎn)速時(shí)，系統(tǒng)主要呈現(xiàn)出與轉(zhuǎn)速同步的周期運(yùn)動(dòng)特性。隨著轉(zhuǎn)速的提高，穩(wěn)態(tài)響應(yīng)在某些參數(shù)下出現(xiàn)分叉、跳躍以及混沌等非線性現(xiàn)象。

3）系統(tǒng)的徑向振動(dòng)出現(xiàn)明顯的工頻和倍頻成分，并且在某些轉(zhuǎn)速和不對(duì)中量下可能出現(xiàn)高次和低次諧波分量，軸心軌跡呈現(xiàn)“香蕉”型。

4）交角不對(duì)中能夠引起系統(tǒng)失穩(wěn)，使轉(zhuǎn)子的橫向振幅增大，因此在實(shí)際系統(tǒng)中， 應(yīng)該嚴(yán)格控制轉(zhuǎn)子間的角不對(duì)中量。

5）在一定轉(zhuǎn)速下，隨著質(zhì)量偏心的增加，作用在轉(zhuǎn)子上的不平衡激勵(lì)增大，轉(zhuǎn)子的穩(wěn)態(tài)振動(dòng)中出現(xiàn)明顯的倍周期分叉以及振幅跳躍現(xiàn)象。

參考文獻(xiàn)：

[1] Lu Siliang， He Qingbo， Wang Jun.A review of stochastic resonance in rotating machine fault detection[J]. Mechanical Systems and Signal Processing，2019，116：230-260.

[2] 郭玉杰，李 克，石 峰.某臺(tái)汽輪發(fā)電機(jī)聯(lián)軸器不對(duì)中振動(dòng)故障的診斷處理[J].軸承，2013，8：56-58.

Guo Yujie，Li Ke，Shi Feng.Diagnosis and treatment for vibration fault of coupling misalignment on a turbine generator[J].Bearing，2013，8： 56-58.

[3] Li Zigang， Jiang Jun， Tian Zhui.Stochastic dynamics of a nonlinear misaligned rotor system subject to random fluid-induced forces[J].Journal of Computational and Nonlinear Dynamics，2017， 12： 011004.

[4] 李自剛，李 明.具有軸承不對(duì)中故障的柔性非圓截面多轉(zhuǎn)子系統(tǒng)非線性動(dòng)力學(xué)行為[J].振動(dòng)工程學(xué)報(bào)，2012，25（1）：68-73.

Li Zigang，Li Ming.Nonlinear dynamics of flexible multi-rotor system with a misaligned journal bearing and non-circular cross section of shaft[J].Journal of Vibration Engineering，2012，25（1）：68-73.

[5] Felipe Wenzel， Katia Lucchesi Cavalca， Vibration signatures of a rotor-coupling-bearing system under angular misalignment[J].Mechanism and Machine Theory， 2019， 133： 559-583.

[6] 陳予恕， 李 軍.汽輪發(fā)電機(jī)組軸系彎扭耦合振動(dòng)問題研究綜述 [J].汽輪機(jī)技術(shù)，2012，54（3）：161-164.

Chen Yushu，Li Jun.Research on bending-torsion coupling vibration of shaft system of turbine generator[J].Turbine Technology，2012，54（3）：161-164.

[7] Slim Bouaziz，Molka Attia Hili，Mohamed Mataar，et al.Dynamic behaviour of hydrodynamic journal bearings in presence of rotor spatial angular misalignment[J].Mechanism and Machine Theory，2009，44：1548-1559.

[8] Redmond I.Study of a misaligned flexibly coupled shaft system having nonlinear bearings and cyclic coupling stiffness—Theoretical model and analysis [J].Journal of Sound and Vibration，2010，329：700-720.

[9] Al-Hussain K M.Dynamic stability of two rigid rotors connected by a flexible coupling with angular misalignment[J].Journal of Sound and Vibration，2003，266（2）：217-234.

[10] Palolo Pennacchi， Andrea Vania， Steven Chatterton. Nonlinear effects caused by coupling misalignment in rotors equipped with journal bearings[J].Mechanical Systems and Signal Processing， 2012，30：306-322.

[11] 李 明.轉(zhuǎn)角不對(duì)中故障的轉(zhuǎn)子系統(tǒng)非線性動(dòng)力學(xué)特征[J].振動(dòng)、測(cè)試與診斷，2011，31（5）：552-556.

Li Ming.Nonlinear dynamics characteristics of rotor system with angular misalignment[J].Journal of Vibration， Measurement & Diagnosis，2011，31（5）：552-556.

[12] Li Ming.Nonlinear vibration of a multirotor system connected by a flexible coupling subjected to the holonomic constraint of dynamic angular misalignment [J].Mathematical Problems in Engineering，2012：243758.

[13] Wang Nanfei， Jiang Dongxiang.Vibration response characteristics of a dual-rotor with unbalance- misalignment coupling faults： theoretical analysis and experimental study[J].Mechanism and Machine Theory，2018，125： 207-219.

[14] Li Zigang，Jiang Jun，Tian Zhui.Nonlinear vibration of an angular-misaligned rotor system with uncertain parameters[J].Journal of Vibration and Control，2016，22（1）：129-144.

[15] 朱擁勇，王德石.含萬向鉸偏斜旋轉(zhuǎn)軸的組合共振及其穩(wěn)定性分析[J].應(yīng)用力學(xué)學(xué)報(bào)，2012，29（4）：361-367.

Zhu Yongyong，Wang Deshi.The stability analysis on combination resonance of misaligned rotary shafts driven through a universal joint[J].Chinese Journal of Applied Mechanics，2012，29（4）：361-367.

[16] Farzad Hemmati， Mohammad Miraskari， Mohamed S Gadala.Dynamic analysis of short and long journal bearings in laminar and turbulent regimes， application in critical shaft stiffness determination[J]. Applied Mathematical Modelling， 2017，48：451-475.

[17] Pawe Pietkiewicz，Sawomir Banaszek，Grzegorz ywica.The stability loss of the rotor-slide bearings system under random load variations[J].Journal of Vibroengineering，2017，19（7）：4921-4935.

[18] 王延博.大型汽輪發(fā)電機(jī)組軸系不對(duì)中振動(dòng)的研究[J].動(dòng)力工程，2004，24：768-774.

Wang Yanbo. Shaft mis-alignment induced vibration in high capacity steam-turbine generator sets[J]. Power Engineering，2004，24：768-774.

Abstract： The non-linear dynamical behaviors caused by the faults of angular misalignment and mass unbalance of the journal bearing supported flexible rotor-bearing system are investigated in this paper. Firstly， the motion equation of the rotor system is derived by the second class Lagrange method with the effects of the motion relation between the two rotors， in which a rheonomic holonomic constraint characterized by the angular misalignment is embedded. Then its nonlinear vibration features， such as the frequency spectra， rotor orbit and Poincaré section， and the largest Lyapunov exponent under different conditions are primarily analyzed based on numerical techniques. The results reveal that the vibration of the rotor system are synchronous at low speed， while as the rotating speed increase， some non-linear characteristics such as jumping， period-doubling and chaotic motions will occur. Furthermore， the parameter analysis for the effects of angular misalignment and mass unbalance is also discussed.

Key words： nonlinear dynamics; rotor-bearing systems; angular misalignment; bifurcation and chaos; mass unbalance

作者簡介： 李自剛 （1986-），男，講師。電話： （029）85583131; E-mail： lzghsfy@hotmail.com