劉麗群，劉 挺，聞邦椿

（1.東北石油大學(xué) 電子科學(xué)學(xué)院，黑龍江 大慶 163318；2.東北大學(xué) 機械工程與自動化學(xué)院，遼寧 沈陽 110819）

轉(zhuǎn)子－軸承系統(tǒng)的松動是轉(zhuǎn)子系統(tǒng)中常見的故障形式，通?？煞譃樾D(zhuǎn)部件松動和支承松動兩種形式.旋轉(zhuǎn)部件松動是指支承系統(tǒng)配合間隙過大或配合過盈量不足引起的振動，而支承松動是指配合面的連接螺栓緊固不牢發(fā)生的異常振動.在發(fā)生支承松動的轉(zhuǎn)子－軸承系統(tǒng)中，支承松動會隨著轉(zhuǎn)子－軸承系統(tǒng)的運轉(zhuǎn)而不斷發(fā)展，系統(tǒng)的運動狀態(tài)會隨著松動的不斷發(fā)展而發(fā)生顯著的改變.當(dāng)系統(tǒng)運轉(zhuǎn)時所產(chǎn)生的不平衡力大于支座所受到的重力時，支座就會跳動并產(chǎn)生碰摩，導(dǎo)致系統(tǒng)的剛度發(fā)生變化，從而使得整個轉(zhuǎn)子系統(tǒng)的運動變得更加復(fù)雜.

現(xiàn)有的松動故障模擬，多采用分段線性剛度和阻尼模型，主要研究內(nèi)容集中在松動－碰摩耦合故障非線性動力學(xué)特性［1］.羅躍綱等［2］研究了松動與碰摩耦合故障對轉(zhuǎn)子動力學(xué)特性的影響，昂雪野等［3］考慮了定子本身的運動，王宗勇等［4］系統(tǒng)地研究了質(zhì)量慢變對于松動－碰摩故障的影響.另外考慮到軸承的影響，很多學(xué)者采用了非線性油膜力模型來模擬軸承，劉長利等［5～8］研究了短軸承非線性油膜力的影響，采用分段線性光滑碰摩模型，運用數(shù)值仿真方法模擬了松動－碰摩耦合故障.陳果［9］研究了滾動軸承的間隙、非線性赫茲接觸力以及變?nèi)嵝哉駝樱\用數(shù)值積分方法分析了轉(zhuǎn)子旋轉(zhuǎn)速度、滾動軸承間隙、碰摩剛度、轉(zhuǎn)子偏心量及軸承座質(zhì)量對系統(tǒng)動力學(xué)響應(yīng)的影響.馬輝［10］等采用有限元模型基于接觸理論研究了轉(zhuǎn)定子系統(tǒng)松動－碰摩故障.

從以上文獻分析可以看出，有關(guān)松動—碰摩故障的研究多限于固定的松動間隙，不同松動間隙的研究尚不多見，而轉(zhuǎn)子系統(tǒng)松動故障在發(fā)展過程中，其松動間隙會不斷的增大，由于松動間隙的增大導(dǎo)致轉(zhuǎn)子系統(tǒng)的松動—碰摩耦合故障，與固定松動間隙的碰摩耦合故障具有很大的不同，因此研究轉(zhuǎn)子系統(tǒng)松動發(fā)展過程中的碰摩問題，具有很重要的實際應(yīng)用價值.

1 動力學(xué)模型和微分方程

圖1為帶有支承松動－碰摩的轉(zhuǎn)子軸承系統(tǒng)，兩端由相同的油膜支承.假定右端支座出現(xiàn)了松動，O1為圓盤的幾何中心，Oe為轉(zhuǎn)盤的質(zhì)心，ma為轉(zhuǎn)盤處的等效集中質(zhì)量，O2和O3分別為兩段軸頸的幾何中心，兩軸頸處的等效集中質(zhì)量相等，均為mb，支座質(zhì)量為m，而松動的最大間隙值為δ態(tài)，假設(shè)軸頸與圓盤之間為無質(zhì)量彈性軸，其剛度系數(shù)用k來表示，c1為轉(zhuǎn)盤處阻尼系數(shù)，c2為軸頸處阻尼系數(shù)，基礎(chǔ)作用于支座O4的阻尼和剛度系數(shù)分別為c4和k4.

設(shè)盤心O1和軸頸中心O2，O3的位移分別為（X1，Y1），（X2，Y3），（X3，Y3），由于松動端支座在水平方向的位移很小，在此僅考慮其在鉛垂方向的位移Y4.設(shè)ω為軸的轉(zhuǎn)動角速度，e為轉(zhuǎn)子的偏心量，μ為油黏度，R為軸頸半徑，L為軸頸長，F(xiàn)X，F(xiàn)Y為滑動軸承作用在轉(zhuǎn)軸上的非線性油膜力，PX（x，y）和PY（x，y）分別為轉(zhuǎn)子系統(tǒng)碰摩力在x，y方向的分量，s為轉(zhuǎn)子的徑向位移s2＝x2＋y2.當(dāng)徑向位移大于轉(zhuǎn)子與定子之間的最大間隙時，轉(zhuǎn)子和定子發(fā)生碰摩.而此時的碰摩力可以表示為

圖1 支承松動－碰摩轉(zhuǎn)子系統(tǒng)示意圖Fig.1 Pedestal looseness and rub-impact rotor system schematic

此時系統(tǒng)的運動微分方程可以表示為

2 松動不同發(fā)展過程中的實驗分析

在實驗室模擬松動發(fā)生發(fā)展的不同的階段，根據(jù)其松動間隙，將其分為4個階段進行研究，即松動開始（松動幾乎不起作用），松動初期（松動特征不明顯），松動中期（松動變大以后），松動后期（松動繼續(xù)增大）.從圖2松動各個階段的單蹤時域分析中，我們可以看到，松動初期系統(tǒng)的振幅增大，而當(dāng)系統(tǒng)發(fā)生碰摩時，由于系統(tǒng)剛度和阻尼的增大，此時系統(tǒng)的最大振幅減小，而一旦系統(tǒng)運行平穩(wěn)，在松動的影響下，又會很快的出現(xiàn)碰摩，松動開始時的振幅曲線與系統(tǒng)發(fā)生碰摩時的振幅曲線明顯不同.

從圖3松動的各個階段的三維譜陣分析中可以看出，隨著松動的不斷發(fā)展，由于松動的最大間隙值增大，發(fā)生碰摩的轉(zhuǎn)速范圍越來越大，而支座發(fā)生較大幅值振動的轉(zhuǎn)速范圍越來越小.隨著轉(zhuǎn)子系統(tǒng)碰摩的發(fā)生，系統(tǒng)會出現(xiàn)高頻振動.

3 結(jié)論

（1）對于松動故障初期而言，主要影響松動端的振動，相對于整個轉(zhuǎn)子系統(tǒng)的振動影響不大.

（2）由于松動的不斷發(fā)展，松動的最大間隙值增大，發(fā)生碰摩的轉(zhuǎn)速范圍越來越大，而支座發(fā)生較大幅值振動的轉(zhuǎn)速范圍越來越小.

圖2 轉(zhuǎn)子松動各個階段的單蹤時域分析Fig.2 Analysis of single track time at all stages of the loose rotor

圖3 轉(zhuǎn)子松動的各個階段的三維譜陣分析（最大熵法）Fig.3 Three dimensional spectrum matrix analysis（MEM）with each stage of the loose rotor

（3）在相同參數(shù)條件下，對支承松動慢變引起碰摩的耦合故障轉(zhuǎn)子系統(tǒng)，與只存在支承松動或碰摩的單一慢變故障轉(zhuǎn)子系統(tǒng)相比，在支承松動慢變引起的碰摩轉(zhuǎn)子系統(tǒng)中，碰摩發(fā)生的轉(zhuǎn)速范圍和振動形式會有所不同，但總體運動趨勢基本一致，轉(zhuǎn)子碰摩能夠減小松動引起的低頻振動，主要激發(fā)高頻振動.

［1］聞邦椿，武新華，丁千，等.故障旋轉(zhuǎn)機械非線性動力學(xué)的理論與試驗［M］.北京：科學(xué)出版社，2004.

WEN Bangchun，WU Xinhua，DING Qian，et al.Theory and experiment of nonlinear dynamics for rotating machinery with faults［M］.Beijing：Science Press，2004.

［2］羅躍綱，曾海泉，李振平，等.基礎(chǔ)松動－碰摩轉(zhuǎn)子系統(tǒng)的混沌特性研究［J］.振動工程學(xué)報，2003，16（2）：184－188.

LUO Yuegang，ZENG Haiquan，LI Zhenping，et al.Study on chaos behaviors of rotor systems with local rubbing and pedestal looseness faults［J］.Journal of Vibration Engineering，2003，16（2）：184－188.

［3］昂雪野，張松鶴.轉(zhuǎn)子－定子系統(tǒng)松動－碰摩耦合故障的動態(tài)響應(yīng)分析［J］.汽輪機技術(shù)，2004，46（4）：275－277.

ANG Xueye，ZHANG Songhe.Analysis on dynamic responses of rotor-stator system with coupling faults of pedestal looseness and rub-impact［J］.Turbine Technology，2004，46（4）：275－277.

［4］王宗勇，吳敬東，聞邦椿.質(zhì)量慢變轉(zhuǎn)子系統(tǒng)的松動與碰摩故障研究［J］.振動工程學(xué)報，2005，18（2）：167－171.

WANG Zongyong，WU Jingdong ，WEN Bangchun.Research on pedestal looseness and rub-impact faults of rotor system with slowly varying mass［J］.Journal of Vibration Engineering，2005，18（2）：167－171.

［5］劉長利，姚紅良，羅躍綱，等.松動碰摩轉(zhuǎn)子軸承系統(tǒng)周期運動穩(wěn)定性研究［J］.振動工程學(xué)報，2004，17（3）：336－340.

LIU Changli，YAO Hongliang，LUO Yuegang，et al.Dynamics of rotor-bearing system with coupling faults of pedestal looseness and rub-impact［J］.Journal of Vibration Engineering，2004，17（3）：336－340.

［6］羅躍綱，何韶君，聞邦椿.松動對碰摩轉(zhuǎn)子－軸承系統(tǒng)非線性特性的影響研究［J］.機械科學(xué)與技術(shù)，2005，24（1）：75－78.

LUO Yuegang，HE Shaojun，WEN Bangchun.Study on effects of pedestal looseness on nonlinear characteristics of rotor-bearing system with rub-impact fault［J］.Mechanical Science and Technology，2005，24（1）：75－78.

［7］孫保蒼，葉啟海，邱飛宇.含松動與碰摩的轉(zhuǎn)子－軸承系統(tǒng)非線性行為分析［J］.潤滑與密封，2007，32（5）：30－33.

SUN Baocang，YE Qihai，QIU Feiyu.Nonlinear dynamical behaviors of the rotor-bearing system with rub and pedestal looseness faults［J］.Lubrication Engineering，2007，32（5）：30－33.

［8］張靖，聞邦椿.松動－碰摩耦合故障轉(zhuǎn)子系統(tǒng)振動特性分析［J］.機械科學(xué)與技術(shù)，2005，24（2）：182－185.

ZHANG Jing，WEN Bangchun.Vibration characteristic of a rotor-bearing system with pedestal looseness and rub-impact［J］.Mechanical Science and Technology，2005，24（2）：182－185.

［9］陳果.含不平衡－碰摩－基礎(chǔ)松動耦合故障的轉(zhuǎn)子—滾動軸承系統(tǒng)非線性動力響應(yīng)分析［J］.振動與沖擊，2008，27（9）：100－104.

CHEN Guo.Nonlinear dynamic response analysis of rotorball bearing system including unbalance-rubbing-looseness coupled faults［J］.Journal of Vibration and Shock，2008，27（9）：100－104.

［10］馬輝，太興宇，汪博，等.松動—碰摩耦合故障轉(zhuǎn)子系統(tǒng)動力學(xué)特性分析［J］.機械工程學(xué)報，2012，48（19）：80－86.

MA Hui，TAI Xingyu，WANG Bo，et al.Dynamic characteristic analysis of a rotor system with pedestal looseness coupled rub-impact fault［J］.Journal of Mechanical Engineering，2012，48（9）：80－86.