轉(zhuǎn)子-滑動軸承系統(tǒng)動力學相似性研究

王永亮， 崔 穎， 韓 聿， 曾之祿

(1.大連海事大學 輪機工程學院，大連 116026； 2.中核集團中核核電運行管理有限公司，海鹽 314300)

轉(zhuǎn)子-滑動軸承系統(tǒng)動力學相似性研究

王永亮1， 崔 穎1， 韓 聿1， 曾之祿2

(1.大連海事大學 輪機工程學院，大連 116026； 2.中核集團中核核電運行管理有限公司，海鹽 314300)

針對轉(zhuǎn)子-滑動軸承系統(tǒng)縮比模型與原型是否滿足動力學相似的問題，采用量綱分析法建立了考慮陀螺力矩和滑動軸承非線性油膜力的轉(zhuǎn)子-軸承系統(tǒng)相似準則，確立了模型與原型各物理量相似比。理論研究表明，通過采用模化轉(zhuǎn)子滑動軸承靜載荷補償措施，可使轉(zhuǎn)子-軸承系統(tǒng)滿足動力學相似要求。補償處理后的模型和原型轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速、失穩(wěn)轉(zhuǎn)速、不平衡響應均具有相似性。并通過算例對比分析轉(zhuǎn)子幾何比、材料密度?；群蛷椥阅Ａ磕；葘S系不平衡響應特性相似性的影響規(guī)律，驗證了所推導的轉(zhuǎn)子動力學相似準則的正確性。

動力學相似；轉(zhuǎn)子動力學；陀螺效應；滑動軸承；非線性油膜力

轉(zhuǎn)子是大型汽輪發(fā)電機組、給水泵、風機、重型燃氣輪機等旋轉(zhuǎn)機械的核心部件，其動力學行為關(guān)系到設備運行穩(wěn)定性和安全性。轉(zhuǎn)子動力學的理論分析和數(shù)值仿真方面的研究已經(jīng)取得了大量有意義的成果，但由于某些聯(lián)接件動力學特性的復雜性和非線性，模型實驗成為揭示轉(zhuǎn)子動力學現(xiàn)象的科學規(guī)律和定性分析動力學特性的必不可少的手段[1-3]。對于實際大型旋轉(zhuǎn)機械轉(zhuǎn)子系統(tǒng)，如果能夠通過縮比模型實驗定量反映原型的動力學特性，將對其動力學設計和故障分析的實踐具有重要意義。

對于如何建立縮比轉(zhuǎn)子-軸承系統(tǒng)和原型動力學特性之間定量關(guān)系的問題，國內(nèi)外公開發(fā)表的相關(guān)理論方法和分析方面的資料均不多見。胡培民[4]理論分析了轉(zhuǎn)子彎曲振動的試驗相似律，研究表明，轉(zhuǎn)子相似模型固有頻率與其幾何相似比成反比關(guān)系。文獻[5]通過對燃氣輪機模塊式氦冷反應堆(GT-MHR)軸系縮比模型進行試驗，分析其支撐電磁軸承在整個工作轉(zhuǎn)速范圍內(nèi)的冗余特性和在線維護功能。WU[6]根據(jù)轉(zhuǎn)子-軸承系統(tǒng)橫向運動方程，推導了縮比模型和原型轉(zhuǎn)子之間相似參數(shù)的比例，以及動力學相似比。羅忠等[7]基于相似理論，結(jié)合量綱分析和方程分析法，建立了轉(zhuǎn)子-軸承系統(tǒng)的動力學相似關(guān)系，求解了相似因子，通過有限元法數(shù)值驗證了相似關(guān)系，考察了軸承參數(shù)對固有頻率和振幅的影響。陳廣凱等[8]采用方程分析法，建立了考慮轉(zhuǎn)子非線性剛度和局部碰摩的Jeffcott轉(zhuǎn)子系統(tǒng)碰摩故障特性相似關(guān)系，對比分析了模型和原型轉(zhuǎn)子系統(tǒng)隨轉(zhuǎn)速變化的分岔特性，研究結(jié)果表明，兩相似轉(zhuǎn)子系統(tǒng)分岔圖趨勢一致。

然而以上轉(zhuǎn)子動力學相似研究中均未考慮陀螺力矩對臨界轉(zhuǎn)速和響應特性等的影響。從動力學方程角度看，轉(zhuǎn)子動力學與結(jié)構(gòu)動力學方程的主要區(qū)別在于：①轉(zhuǎn)子動力學方程中包含隨轉(zhuǎn)速變化的陀螺力矩陣；②轉(zhuǎn)子系統(tǒng)不平衡質(zhì)量力隨轉(zhuǎn)速變化；③轉(zhuǎn)子支撐剛度和阻尼系數(shù)可能隨轉(zhuǎn)子轉(zhuǎn)速變化(對于滑動軸承來說)。由此引出兩個問題：一是陀螺效應和隨轉(zhuǎn)速變化的不平衡質(zhì)量力對轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速和響應特性相似性的影響;二是相似理論是否適用于滑動軸承支撐的轉(zhuǎn)子系統(tǒng)動力學特性分析。

本文針對上述問題進行探討，分析考慮陀螺力矩和滑動軸承非線性油膜力的轉(zhuǎn)子動力學相似規(guī)律。利用Timoshenko梁理論和短圓瓦滑動軸承非線性力模型，建立轉(zhuǎn)子-滑動軸承系統(tǒng)的動力學方程。采用量綱分析法導出考慮陀螺力矩和非線性油膜力的轉(zhuǎn)子-軸承系統(tǒng)動力學相似準則，給出轉(zhuǎn)子系統(tǒng)模型與原型各個物理量相似比。通過對一實際實驗轉(zhuǎn)子-圓柱瓦滑動軸承系統(tǒng)進行原型和不同?；认履Ｐ偷膭恿W特性數(shù)值仿真，對比分析幅頻特性和分岔特性，驗證所給出的考慮陀螺效應和滑動軸承非線性油膜力的轉(zhuǎn)子動力學相似準則的正確性。

1 轉(zhuǎn)子-滑動軸承系統(tǒng)動力學模型

1.1 轉(zhuǎn)子模型

目前，有限元法已廣泛被應用于轉(zhuǎn)子動力學建模中，其中的二維軸對稱準三維單元和三維實體單元已被應用于轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速和振型分析中[9-10]，而一維梁單元具有力學概念清晰、建模方便，總自由度相對較少，在以往轉(zhuǎn)子-系統(tǒng)臨界轉(zhuǎn)速和振型分析、線性或非線性響應特性研究中占據(jù)著主導地位。本文在對轉(zhuǎn)子進行有限元離散后，采用考慮了轉(zhuǎn)動慣量效應和橫截面剪切變形影響的Timoshenko梁單元對軸段進行動力學建模[11]，并將輪盤作為剛體進行動力學分析，得到整個軸系的動力學模型如下：

(1)

式中：[M]、[C]、[J]、[K]分別為轉(zhuǎn)子的質(zhì)量矩陣、阻尼矩陣、陀螺力矩陣和剛度矩陣，均為4N*4N矩陣，N為節(jié)點數(shù)。{G}、{Q(t)}、{Fb}分別為各節(jié)點重力分量、質(zhì)量偏心產(chǎn)生的不平衡力和滑動軸承支撐油膜力。Ω為轉(zhuǎn)軸旋轉(zhuǎn)速度。

廣義位移向量為：

xi,yi,-θxi,θyi表示軸段第i節(jié)點橫截面水平和豎直方向的位移和轉(zhuǎn)角。

重力為：

不平衡質(zhì)量力為：

{Q(t)}=[{m}eΩ2sin(Ωt),{m}eΩ2cos(Ωt)]T，其中{m}={m1,0,m2,0,…,mN,0}。

質(zhì)量矩陣、陀螺矩陣和剛度矩陣形式如下[12]：

式中:mi、Jdi和Jpi分別為第i節(jié)點附加質(zhì)量、附加直徑轉(zhuǎn)動慣量和極轉(zhuǎn)動慣量。[M1]為整體質(zhì)量矩陣，其具體形式如圖1所示，為2N*2N階對稱帶狀矩陣。矩陣[J1]和[K1]的形式和[M1]的形式類似，對于[J1]來說，用[J11]s(i)、[J12]s(i)、[J21]s(i)和[J22]s(i)替換圖中[m11]s(i)、[m12]s(i)、[m21]s(i)和[m22]s(i)，并將[Md](i)改為[Gd](i)即可；而對于[K1]來說，用[k11]s(i)、[k12]s(i)、[k21]s(i)和[k22]s(i)替換圖中[m11]s(i)、[m12]s(i)、[m21]s(i)和[m22]s(i)，并令[Md](i)=[0]即可。

圖1 轉(zhuǎn)子系統(tǒng)質(zhì)量矩陣[M1]Fig.1 Mass matrix of rotor

轉(zhuǎn)子第i個軸段的質(zhì)量矩陣、剛度矩陣、陀螺力矩陣表達式為[13]：

ρ、A、l、E、I、G分別為密度、橫截面面積、軸段長度、彈性模量、截面慣性矩、剪切模量。

對于橫截面為圓或圓環(huán)的轉(zhuǎn)軸，剪切校正因子為[14]：

而轉(zhuǎn)軸阻尼矩陣[C]采用經(jīng)典的Rayleigh阻尼模型。其表達式如下：

[C]=α[M]+β[K]

式中：ω1、ω2分別為轉(zhuǎn)子的第1、2階固有頻率；ξ1、ξ2分別為第1、2階模態(tài)的阻尼系數(shù)。

1.2 滑動軸承非線性油膜力模型

圓瓦徑向滑動軸承非線性油膜力求解方法有數(shù)值求解雷諾方程，通過數(shù)值積分油膜壓力得到油膜力、油膜力數(shù)據(jù)庫法、變分法模型、長軸承模型、短軸承模型、有限長軸承模型等[15]?？紤]到本文研究的重點在于轉(zhuǎn)子動力學相似性理論及其驗證，因此，采用計算快速簡便、已被廣泛應用于非線性動力學現(xiàn)象分析的短軸承模型[16]，其模型表達式為：

式中：

2 轉(zhuǎn)子動力學相似關(guān)系

針對轉(zhuǎn)子-滑動軸承系統(tǒng)動力學方程(1)，本節(jié)通過量綱理論建立各物理量間相似關(guān)系，分析其相似特性。

根據(jù)相似性理論[17]和轉(zhuǎn)子-軸承系統(tǒng)動力學理論，式(1)轉(zhuǎn)子系統(tǒng)中各物理量共有12個，分別為彈性模量E、材料密度ρ、幾何尺寸l、泊松比ν，重力G、潤滑油粘度μ，時間t、阻尼比ξ、固有頻率ωn、位移響應u、轉(zhuǎn)速Ω和特征轉(zhuǎn)速(臨界轉(zhuǎn)速、油膜失穩(wěn)轉(zhuǎn)速)Ωcha，對于線彈性范圍轉(zhuǎn)子系統(tǒng)振動，各物理量間的關(guān)系可表達為：

D(E,ρ,l,ν,G,μ,t,ξ,ωn,u,Ω,Ωcha)=0

(2)

選擇力[F]、長度[L]和時間[T]作為基本量綱，根據(jù)相似理論，獨立π準則的數(shù)量為9個。由相似第二定理可將式(2)轉(zhuǎn)化為無因次的相似準則π1、π2、…、π9之間的函數(shù)：

f(π1,π2,π3,π4,π5,π6,π7,π8,π9)=0

(3)

用a1、a2、…、a12分別代表物理量E、ρ、l、ν、G、μ、t、ξ、ωn、u、Ω和Ωcha中的指數(shù)，則πi的形式為：

(4)

轉(zhuǎn)子系統(tǒng)各物理量間量綱矩陣可表示為：

表1 各物理量間量綱矩陣

根據(jù)量綱矩陣，可得以下三個線性齊次代數(shù)方程：

量綱[F]：

a1+a2+a5+a6=0

(5)

量綱[L]：

-2a1-4a2+a3-2a6+a10=0

(6)

量綱[T]：

2a2+a6+a7-a9-a11-a12=0

(7)

3個方程無法解出12個未知數(shù)，將未知數(shù)中的a1、a2、a3表達為a4、a5、…、a12的函數(shù)，則有：

(8)

(9)

a3=-2a5-a6-a7+a9-a10+a11+a12

(10)

轉(zhuǎn)子各物理量指數(shù)可以由π矩陣來表示，根據(jù)式(8)～(10)將各個物理量指數(shù)列成表的形式，如表2所示。

由表2的π矩陣，每一行代表無量綱乘積的一組指數(shù)，因此可建立各獨立的9個π準則，如下所示：

根據(jù)相似定理，縮比模型與原型相似，則兩者的各相似指標為1。

表2 轉(zhuǎn)子系統(tǒng)的π矩陣

Tab.2 π matrix of rotor system

a1a2a3a4a5a6a7a8a9a10a11a12EρlνGμtξωnuΩΩchaπ1000100000000π2-10-2010000000π3-1/2-1/2-1001000000π41/2-1/2-1000100000π5000000010000π6-1/21/21000001000π700-1000000100π8-1/21/21000000010π9-1/21/21000000001

設模型與原型的E、ρ、l、ν、G、μ、t、ξ、ωn、u、Ω和Ωcha的相似比分別為：CE、Cρ、Cl、Cν、CG、Cμ、Ct、Cξ、Cu、CΩ、CΩcha。根據(jù)相似定理，縮比模型與原型相似，則兩者的各相似指標為1，可得：

(11)

對于轉(zhuǎn)子系統(tǒng)模型設計，如果已知幾何相似比Cl、密度相似比Cρ、彈性模量相似比CE，則其它物理量相似比可表示為：

由于重力加速度g無法進行?；?，由此導致實際?；D(zhuǎn)子動力學方程中重力項不滿足表3中重力相似比要求。對于轉(zhuǎn)子-滑動軸承系統(tǒng)，重力與滑動軸承靜載荷相關(guān)，直接影響著軸承的剛度阻尼特性，重力不滿足相似比要求意味著該軸系模型和原型不具有動力學相似性。

表3 轉(zhuǎn)子系統(tǒng)模型與原型各個物理量相似比

然而可通過在模化轉(zhuǎn)子豎直方向上施加額外的補償力(如通過?；D(zhuǎn)子軸線傾斜使轉(zhuǎn)軸在軸承上分力匹配，或通過電磁力等實現(xiàn)補償)，使得滑動軸承的靜載荷滿足相似比要求。

設原型的重力向量為：

{Gp}={m}g

根據(jù)表3中的重力相似比，?；D(zhuǎn)子重力向量應該為：

然而實際?；D(zhuǎn)子重力值為：

兩者之間差即為補償力：

載荷補償處理后的?；D(zhuǎn)子動力學方程可表達為如下形式：

通過補償處理，原型和模型參數(shù)滿足了相似比要求。由表3可知，在幾何、密度和彈性模量相似比確定的情況下，固有頻率和轉(zhuǎn)速相似比相等，即模型在ΩCΩ轉(zhuǎn)速下固有頻率等于原型在Ω轉(zhuǎn)速下固有頻率的Cωn(即CΩ)倍，模型的特征轉(zhuǎn)速(臨界轉(zhuǎn)速、失穩(wěn)轉(zhuǎn)速等)必然為原型的Cωn(或CΩ)倍。

3 實例驗證

以文獻[18]中實際實驗臺轉(zhuǎn)子為例，對比分析模化參數(shù)比對轉(zhuǎn)子-軸承系統(tǒng)非線性動力學特性影響規(guī)律，以驗證所建立的轉(zhuǎn)子動力學相似律。

3.1 轉(zhuǎn)子-軸承系統(tǒng)參數(shù)

圖2為轉(zhuǎn)子系統(tǒng)結(jié)構(gòu)圖，該轉(zhuǎn)子四個輪盤通過錐套套裝在轉(zhuǎn)軸上，總長0.993 m，跨距0.738 m，總質(zhì)量約168 kg。

圖2 實驗臺轉(zhuǎn)子結(jié)構(gòu)參數(shù)Fig.2 Configuration of the rotor

對圖2轉(zhuǎn)子進行有限元離散，圖3為轉(zhuǎn)子有限元模型示意圖，表4為軸端劃分數(shù)據(jù)列表，共分為11個軸段，12個節(jié)點，四個剛性輪盤的質(zhì)量和轉(zhuǎn)動慣量作為附加質(zhì)量加入動力學方程中，3號節(jié)點和10號節(jié)點分別為#1和#2圓瓦滑動軸承支撐位置。表5為圓瓦滑動軸承參數(shù)列表。

圖3 模化轉(zhuǎn)子節(jié)點劃分示意圖Fig.3 Structure sketch of finite element model of the rotor bearing system

Tab.4 Finite element model of rotor modeling data list

段號段長/m外徑/m段重/kg節(jié)點編號附加質(zhì)量/kg節(jié)點總質(zhì)量/kg附加直徑轉(zhuǎn)動慣量/(kg·m2)支撐10.0630.03800.6054100.3025020.0620.04000.6116200.6085030.0650.04000.6412300.62650*40.1000.04461.2264400.8905050.1400.04551.7869539.140.4800.3960.1280.04571.6482639.140.6390.3970.1400.04551.7869739.140.6390.3980.1000.04461.2263839.140.4800.3990.0650.04000.6412900.89050100.0620.04000.61161000.62650*110.0680.03800.60541100.608501200.30250sum0.993168.257

表5 圓瓦滑動軸承參數(shù)

3.2 不平衡響應仿真結(jié)果

采用1轉(zhuǎn)子-滑動軸承系統(tǒng)動力學模型相關(guān)理論對上述實驗轉(zhuǎn)子系統(tǒng)進行動力學建模。其中前兩階模態(tài)阻尼系數(shù)取ξ1=0.02，ξ2=0.05，設定#3輪盤處存在大小為203.195(g·cm)的不平衡質(zhì)量。

采用Newmark法[19]在整個時間域內(nèi)對轉(zhuǎn)子-軸承系統(tǒng)動力學方程(1)進行逐步積分求解，獲取系統(tǒng)各個時刻的響應數(shù)據(jù)。數(shù)值求解過程中，計算轉(zhuǎn)速為1 000～6 000 r/min, 各轉(zhuǎn)速下計算200個周期，計算步長為每轉(zhuǎn)速周期均分160 000個計算點，每個周期輸出1 000個計算點，取后50個周期作為分析數(shù)據(jù)，Newmark程序迭代過程中采用前后兩次迭代計算方程右端力相對偏差作為收斂判據(jù)，收斂偏差標準為1E-006。

對獲得的各時刻動力學響應數(shù)據(jù)進行處理分析，利用幅頻特性曲線、分岔圖、三維譜圖、頻譜圖和軸心軌跡圖來分析軸系的非線性動力學響應特性和穩(wěn)定性。

圖4為#1軸承處(節(jié)點3)軸頸x向振動圖譜。

圖4 #1軸承處軸頸x方向振動圖譜Fig.4 Variation of the vibration of the rotor with the journal angular velocity

由圖4(a)可知，圓瓦滑動軸承支撐下，軸系臨界轉(zhuǎn)速為2 070 r/min(同時在1 800 r/min時出現(xiàn)小幅度峰值，可能是由于軸承x向和y向主剛度的差異及交叉剛度引起的)。當轉(zhuǎn)速超過5 150 r/min后，振幅突增。從分岔圖可明顯看出，分岔點約為5 050 r/min，該轉(zhuǎn)速以前軸系處于周期1運動，當轉(zhuǎn)速大于5 050后，系統(tǒng)失去了原有動力穩(wěn)定性，軸系運動變的復雜。從圖4(c)三維譜圖可以看出，5 050 r/min以前軸系振動只有基頻，而轉(zhuǎn)速升至5 050后，開始出現(xiàn)明顯低頻分量，且振動幅值急劇增大，之后低頻振動頻率基本鎖定在系統(tǒng)一階彎曲固有頻率值，即發(fā)生鎖頻現(xiàn)象。

圖5和6分別為不同轉(zhuǎn)速時#1軸承處軸頸振動頻譜圖和軸心軌跡圖。

由圖可知，轉(zhuǎn)速為5 000 r/min時，軸系以基頻振動為主，軸心軌跡為橢圓，軸心處于第四象限。當轉(zhuǎn)速升至5 150 r/min時，出現(xiàn)明顯低頻分量，軸心軌跡為概周期運動。轉(zhuǎn)速超過5 250 r/min后(發(fā)生油膜振蕩的轉(zhuǎn)速)，低頻振動幅值超過基頻振動振幅，開始占具主導地位，轉(zhuǎn)速為5 950 r/min時，軸心軌跡近似為圓形，此時軸系運動為周期2運動。

圖5 不同轉(zhuǎn)速下#1軸承處軸頸處x方向頻譜圖(橫坐標為倍頻數(shù)；縱坐標為振動幅值，單位μm)Fig.5Frequency spectrum of bearing #1 in x direction for different journal angular velocity (The horizontal direction is normalized frequency, and vertical is amplitude/μm)

圖6 不同轉(zhuǎn)速下#1軸承處軸頸軸心軌跡圖(單位μm)Fig.6 Orbit for different journal angular velocity of bearing #1 (unit: μm)

3.3 ?；瘏?shù)對動力學特性影響

表6給出了6個不同轉(zhuǎn)子材料密度比、彈性模量比和幾何比的數(shù)值仿真方案，獲取各方案下轉(zhuǎn)子-軸承系統(tǒng)的響應特性，分析?；葘S系幅頻特性和分岔特性的影響規(guī)律，驗證所推導的考慮陀螺力矩和非線性油膜力的轉(zhuǎn)子動力學相似性準則的正確性。

數(shù)值仿真中初始位移、滑動軸承參數(shù)、轉(zhuǎn)子幾何參數(shù)等均按照表3中?；瘻蕜t計算，數(shù)值積分迭代參數(shù)與3.2一致。

表6 不同參數(shù)?；鹊姆桨噶斜?/p>

圖7為不同物理參數(shù)?；认罗D(zhuǎn)子-軸承系統(tǒng)的幅頻特性曲線，圖8為對應的振動分岔特性隨轉(zhuǎn)速變化圖。

圖7 不同方案下#1軸承處軸頸x向幅頻特性(橫坐標為轉(zhuǎn)速，單位r/min；縱坐標為通頻振動幅值，單位μm)Fig.7 Amplitude frequency characteristics of bearing #1 in x direction with different scale ratio (The horizontal direction represents speed/r/min, and vertical is amplitude/μm)

由圖可知，各方案下幅頻特性和分岔特性隨轉(zhuǎn)速變化趨勢相似，臨界轉(zhuǎn)速和失穩(wěn)轉(zhuǎn)速相似，各轉(zhuǎn)速下振動幅值相似。表7為各方案下轉(zhuǎn)子臨界轉(zhuǎn)速、過臨界時振動幅值、發(fā)生油膜振蕩轉(zhuǎn)速和振幅對比，可以看出，模型與原型的特征參數(shù)比與密度、彈性模量和幾何比例相關(guān)，各參數(shù)模化比下轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速與通過原型推算出的臨界轉(zhuǎn)速誤差最大為4.8‰，過臨界時振動幅值誤差最大為0.037‰，發(fā)生油膜振動轉(zhuǎn)速和振幅最大誤差分別為0.014‰和0.251‰。

因此，考慮陀螺力矩和滑動軸承非線性油膜力的轉(zhuǎn)子系統(tǒng)在對重力進行相似性處理后，其動力學特性滿足相似律，具有表3所示的相似關(guān)系。

圖8 不同方案下#1軸承處軸頸x向分岔圖(橫坐標為轉(zhuǎn)速，單位r/min；縱坐標為振幅，單位μm)Fig.8 Bifurcation diagram of bearing #1 in x direction with different scale ratio (The horizontal direction represents speed(r/min, and vertical is amplitude μm)

方案各方案下特征參數(shù)(與通過原型推算出的特征參數(shù)誤差千分比/‰)臨界轉(zhuǎn)速過臨界振幅油膜振蕩轉(zhuǎn)速油膜振蕩振幅12575(4.8)19.14(0.002)6562.5(0.000)36.95(0.016)21030(1.9)47.85(0.006)2625(0.000)92.36(0.046)3927(1.7)47.85(0.006)2362.5(0.000)92.36(0.046)41133(2.1)47.85(0.006)2887.5(0.000)92.36(0.046)51144(2.3)47.77(0.037)2916.7(0.013)91.83(0.251)6936(1.9)47.77(0.037)2386.4(0.014)91.84(0.247)

4 結(jié) 論

(1)采用量綱分析法推導了考慮陀螺效應和非線性油膜力的轉(zhuǎn)子-軸承系統(tǒng)動力學相似準則，確立了轉(zhuǎn)子系統(tǒng)模型與原型各個參數(shù)相似比。

(2)對于轉(zhuǎn)子-軸承系統(tǒng)，轉(zhuǎn)軸的重力對滑動軸承動特性有著重要的影響，由于重力加速度不滿足相似性?；?，致使轉(zhuǎn)子-軸承系統(tǒng)方程理論上不滿足相似性要求。

(3)可通過采取對?；D(zhuǎn)子在軸承處施加外力，或改變軸系傾角等措施，使得模化滑動軸承靜載荷滿足相似性要求，進而使得轉(zhuǎn)子-滑動軸承系統(tǒng)和處理后的?；D(zhuǎn)子與原型具有動力學相似性，其臨界轉(zhuǎn)速、失穩(wěn)轉(zhuǎn)速、不平衡響應等特性參數(shù)之比均具有相似性關(guān)系。

[ 1 ] 馬輝, 太興宇, 汪博, 等. 柔性轉(zhuǎn)子系統(tǒng)輪盤外緣定點碰摩動力學特性分析[J]. 中國電機工程學報,2012,32(17): 89-96. MA Hui, TAI Xingyu, WANG bo, et al. Dynamic characteristic analysis of a flexible rotor system with fixed-point rubbing fault at a wheel edge[J]. Proceedings of the CSEE,2012, 32(17): 89-96.

[ 2 ] 黃梓嫄, 韓邦成, 周銀鋒. 非線性接觸下磁懸浮電機柔性轉(zhuǎn)子系統(tǒng)模態(tài)分析[J]. 中國電機工程學報,2014, 34(15): 2438-2444. HUANG Ziyuan, HAN Bangcheng, ZHOU Yinfeng. Modal analysis of the flexible rotor system of magnetic levitation motors under nonlinear contact[J]. Proceedings of the CSEE,2014, 34(15): 2438-2444.

[ 3 ] WANG Yongliang, LIU Zhansheng, ZHANG Guanghui, et al. Theoretical and experimental study on stability of rotor system with cylindrical bearing and elliptical bearing[R]. ASME Paper, GT2010-23120, 2010.

[ 4 ] 胡培民. 大型轉(zhuǎn)子彎曲振動的比例模型試驗分析[J]. 機械研究與應用,1998, 11(3): 36-38. HU Peimin. Scale model test of bending vibration for large rotor[J]. Mechanical Research & Application,1998,11(3): 36-38.

[ 5 ] BAXI C B, TELENGATOR A, RAZVI J. Rotor scale model tests for power conversion unit of GT-MHR[J]. Nuclear Engineering and Design,2012, 251: 344-348.

[ 6 ] WU J J. Prediction of lateral vibration characteristics of a full-size rotor-bearing system by using those of its scale models[J]. Finite Elements in Analysis and Design, 2007, 43: 803-816.

[ 7 ] 羅忠, 閆宇龍, 韓清凱. 軸承-轉(zhuǎn)子系統(tǒng)的動力學相似及其軸承參數(shù)的影響分析[J]. 振動與沖擊,2012,31(增刊1): 12-16. LUO Zhong, YAN Yulong, HAN Qingkai. Analysis of dynamic similarity of the rotor-bearing system and it’s effects of bearing parameters[J]. Journal of Vibration and Shock, 2012,31(Sup1): 12-16.

[ 8 ] 陳廣凱, 郭健, 時海麗, 等. 考慮局部碰摩的轉(zhuǎn)子系統(tǒng)動力學相似研究[C]//第11屆全國轉(zhuǎn)子動力學學術(shù)討論會, 大連, 2014年8月.

[ 9 ] 龔建政, 鐘芳明, 賀星, 等. 利用ANSYS和坎貝爾圖對燃氣輪機壓氣機轉(zhuǎn)子模態(tài)及臨界轉(zhuǎn)速的分析計算[J]. 現(xiàn)代制造工程,2012(4): 73-77. GONG Jianzheng, ZHONG Fangming, HE Xing, et al. Mode analysis and critical speed calculation on a low-pressure compressor rotor of marine gas turbine based on ANSYS[J]. Modern Manufacturing Engineering,2012(4): 73-77.

[10] 金宏偉, 孫立權(quán). 某型風機轉(zhuǎn)子臨界轉(zhuǎn)速計算和對中研究[J]. 汽輪機技術(shù),2014, 56(5): 343-344. JIN Hongwei, SUN Liquan. Critical speed and alignment research of a certain type blower rotor[J]. Turbine Technology,2014, 56(5): 343-344.

[11] HAN S M, BENAROYA H,WEI T. Dynamics of transversely vibrating beams using four engineering theories[J]. Journal of Sound and Vibration,1999, 225(5): 935-988.

[12] 鐘一諤, 何衍宗, 王正, 李方澤. 轉(zhuǎn)子動力學[M]. 北京：清華大學出版社,1984.

[13] GIANCARLO G. Dynamics of rotating systems[M]. S SpringerVerlag,2005.

[14] COWPER G R. The shear coefficient in timoshenko’s beam theory[J].Journal of Applied Mechanics,1966, 33(2): 335-340.

[15] 王永亮, 劉占生. 圓瓦滑動軸承油膜力近似解析模型[J]. 中國電機工程學報, 2011(29): 110-117. WANG Yongliang, LIU Zhansheng. Approximate analytical model of oil-film force for cylindrical journal bearing[J]. Proceedings of the CSEE,2011(29): 110-117.

[16] CAPONE G.Descizione analitica del campo di forze fluidodinamico nei cuscinetti cillindrici lubrificati[J]. L’Energia Elettrica,1991, 3: 105-110.

[17] 楊俊杰. 相似理論與結(jié)構(gòu)模型試驗[M]. 武漢:武漢理工大學出版社,2005.

[18] 王永亮. 滑動軸承油膜力近似解析模型研究及軸系動力穩(wěn)定性分析[D]. 哈爾濱:哈爾濱工業(yè)大學, 2011.

[19] NEWMARK N M. A method of computation for structural dynamics[J]. Journal of Engineering Mechanics Division,1959, 85(3): 67-94.

Dynamic similarity of rotor-sliding bearing systems

WANG Yongliang1, CUI Ying1, HAN Yu1, ZENG Zhilu2

(1. Marine Engineering College, Dalian Maritime University, Dalian 116026, China;2.CNNC Nuclear Power Operation Management Co, Ltd. Haiyan 314300, China)

Are dynamic characteristics of a rotor-sliding bearing system scaled model similar to those of the original system? To answer this problem, dynamic similarity criteria for rotor-sliding bearing systems considering gyroscopic moment and nonlinear oil film force were derived by using the dimensional analysis method, and the similarity ratios of physical variables of the scaled model to those of the original system were obtained. The theoretical study showed that the dynamic similarity requirements of rotor-sliding bearing systems can be met by using the static load compensation measures of sliding journal bearings of the scaled rotor system; the critical speed, stability threshold speed and imbalance response of the scaled system after compensation and those of the original system have a similarity; the correctness of dynamic similarity criteria of rotor systems proposed here is verified using comparative analysis for imbalance response characteristics of the scaled system model and those of the original system with different rotor geometric ratios, material density ratios and elastic modulus ratios.

dynamic similarity; rotor dynamics; gyroscopic effect; sliding bearing; nonlinear oil film force

國家自然科學基金(51606023;51436002)；遼寧省自然科學基金(2015020130)；中央高?；究蒲袠I(yè)務費專項資金資助(3132016016)

2015-11-05 修改稿收到日期：2016-02-02

王永亮 男，博士，講師，1983年10月生

TH133.3

A

10.13465/j.cnki.jvs.2017.01.023