史永峰，李 明，劉 剛，何 琳

（1.西安科技大學(xué) 理學(xué)院力學(xué)系，西安 710054；2.海軍工程大學(xué) 振動(dòng)與噪聲研究所，武漢 430033）

水潤滑橡膠軸承支承轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性研究

史永峰1，李 明1，劉 剛1，何 琳2

（1.西安科技大學(xué) 理學(xué)院力學(xué)系，西安 710054；2.海軍工程大學(xué) 振動(dòng)與噪聲研究所，武漢 430033）

文章研究了水潤滑橡膠軸承的動(dòng)態(tài)特性及其在橡膠軸承支承下轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)行為。首先，考慮了動(dòng)壓潤滑和橡膠體的彈性變形等因素后，建立了水潤滑橡膠軸承的水彈性動(dòng)力潤滑數(shù)學(xué)模型，根據(jù)載荷增量法與差分法，分析了軸承的剛度和阻尼。然后，選擇集總參數(shù)法離散化轉(zhuǎn)子系統(tǒng)，根據(jù)分析力學(xué)方法推導(dǎo)了水潤滑橡膠軸承耦合轉(zhuǎn)子系統(tǒng)的運(yùn)動(dòng)微分方程。最后，通過廣義逆迭代法，研究了水潤滑橡膠軸承支承的轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性和穩(wěn)定性，得到了轉(zhuǎn)子-軸承系統(tǒng)的坎貝爾曲線，臨界轉(zhuǎn)速與振型以及失穩(wěn)轉(zhuǎn)速等。結(jié)果表明偏心率、間隙比、轉(zhuǎn)速及長徑比對(duì)軸承的動(dòng)態(tài)特性的影響比較顯著；軸承的水彈性潤滑特性對(duì)系統(tǒng)的臨界轉(zhuǎn)速影響較大，橡膠軸承的彈性對(duì)系統(tǒng)的穩(wěn)定性影響次之。

水潤滑橡膠軸承；轉(zhuǎn)子系統(tǒng)；動(dòng)態(tài)特性；臨界轉(zhuǎn)速；穩(wěn)定性

0 引 言

水潤滑與油潤滑相比具有來源方便、使用安全、無污染等優(yōu)點(diǎn)，并且由于水具有較大的比熱容，其冷卻效果要比油好得多[1]，因而水潤滑橡膠軸承被廣泛應(yīng)用于大型船舶的動(dòng)力推進(jìn)系統(tǒng)中[2-3]。近十余年來，隨著我國國防工業(yè)的需求和國民經(jīng)濟(jì)的快速發(fā)展，船舶特別是軍用水面、水下艦艇日趨大型、重載和高速，這就對(duì)水潤滑橡膠軸承在推進(jìn)軸系的設(shè)計(jì)、制造和安裝方面提出了更高的要求。

在水潤滑橡膠軸承方面，考慮了軸承彈性變形的流體潤滑理論研究也有一些進(jìn)展[4]。文獻(xiàn)[5-6]分析了水潤滑橡膠軸承的結(jié)構(gòu)和工作特點(diǎn)，針對(duì)橡膠軸承存在的鳴音問題，提出了軸承在設(shè)計(jì)時(shí)需注意的一些問題；文獻(xiàn)[7]通過數(shù)值方法討論了軸頸傾斜對(duì)徑向滑動(dòng)軸承潤滑性能的影響；文獻(xiàn)[8]基于有限元法提出一種混合水彈潤滑的有限長彈性滑動(dòng)軸承模型，并采用改進(jìn)的松弛迭代方法求解水潤滑軸承的潤滑特性問題；文獻(xiàn)[9-12]分別采用實(shí)驗(yàn)和數(shù)值分析方法研究了不同工況不同類型的水潤滑橡膠軸承的水彈潤滑特性；文獻(xiàn)[13]建立了考慮多溝槽潤滑結(jié)構(gòu)和實(shí)際工況邊界條件的水潤滑橡膠合金軸承的數(shù)學(xué)模型，研究了溝槽半徑對(duì)潤滑性能的影響；文獻(xiàn)[14]以水潤滑橡膠板條結(jié)構(gòu)參數(shù)為切入點(diǎn)，著重研究了橡膠板條厚度對(duì)軸承潤滑性能的影響；文獻(xiàn)[15]針對(duì)水潤滑橡膠軸承動(dòng)態(tài)特性系數(shù)剛度與阻尼測試過程中信號(hào)混雜等問題，采用了頻譜分析法來測定其動(dòng)力參數(shù)，計(jì)算出水潤滑橡膠軸承的動(dòng)力特性；文獻(xiàn)[16]在考慮橡膠軸瓦彈性變形的條件下，建立水潤滑軸承的水彈潤滑模型并進(jìn)行了數(shù)值計(jì)算，從理論上分析了水潤滑橡膠層厚度對(duì)軸承彈流潤滑性能的影響。

在軸承—轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)研究領(lǐng)域，近年來對(duì)各種類型軸承支承轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性研究取得了許多成果，例如，文獻(xiàn)[17]研究了滑動(dòng)軸承—轉(zhuǎn)子系統(tǒng)的動(dòng)力特性，利用Riccati傳遞矩陣法對(duì)轉(zhuǎn)子系統(tǒng)的動(dòng)態(tài)特性進(jìn)行了求解，開發(fā)了相應(yīng)的計(jì)算程序；文獻(xiàn)[18]建立了水潤滑軸承轉(zhuǎn)子系統(tǒng)有限元模型，計(jì)算了水軸承支承下的軸系動(dòng)力學(xué)響應(yīng)；文獻(xiàn)[19]研究了高速動(dòng)靜混合氣體軸承的靜動(dòng)態(tài)特性，同時(shí)結(jié)合非線性動(dòng)力學(xué)的理論和方法，利用數(shù)值仿真和試驗(yàn)手段，對(duì)系統(tǒng)的耦合動(dòng)力學(xué)特性進(jìn)行了研究；文獻(xiàn)[20]以某高壓多級(jí)泵為例，創(chuàng)建了分析該泵軸承—轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性的模型；利用ANSYS計(jì)算生成了剛性支承的坎貝爾圖，獲得了臨界轉(zhuǎn)速。文獻(xiàn)[21]提出了一種考慮空間軸承分布力的新模型，并給出其動(dòng)態(tài)系數(shù)。

綜上所述，國內(nèi)外學(xué)者對(duì)油潤滑軸承支承的轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)研究比較多，而水潤滑橡膠軸承支承轉(zhuǎn)子系統(tǒng)作為一種十分重要的傳動(dòng)系統(tǒng)，有關(guān)其動(dòng)力學(xué)特性方面的研究還相對(duì)薄弱，隨著工程應(yīng)用的不斷深入以及轉(zhuǎn)子轉(zhuǎn)速的不斷提高，水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速、穩(wěn)定性及各種載荷作用下的動(dòng)態(tài)響應(yīng)將有待進(jìn)一步深入研究。本文將主要討論水潤滑橡膠軸承的靜動(dòng)態(tài)特性及其對(duì)水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性的影響。

1 水潤滑橡膠軸承的水彈潤滑特性

水潤滑橡膠軸承結(jié)構(gòu)如圖1所示。其中轉(zhuǎn)子以角速度ω繞軸承中心O轉(zhuǎn)動(dòng)，O′為軸頸中心，e為偏心距，φ為偏位角，hmax為最大水膜間隙，hmin為最小水膜間隙。轉(zhuǎn)子在旋轉(zhuǎn)時(shí)，所形成的水膜將轉(zhuǎn)軸與橡膠襯層隔開，流體動(dòng)壓潤滑在橡膠軸承間隙中形成。

水潤滑橡膠軸承水膜壓力的無量綱形式可表示為

圖1 水潤滑橡膠軸承示意圖Fig.1 Schematic diagram of the water lubrication rubber bearing

式中：θ為軸承在圓柱坐標(biāo)下周向坐標(biāo)，λ為無量綱軸向坐標(biāo)，H為無量綱水膜厚度，P為無量綱水膜力，L為軸承寬度或長度，d為軸頸直徑，τ為時(shí)間。

無量綱所用特征量

與油潤滑軸承的金屬軸套相比，水潤滑橡膠軸承的橡膠襯層變形要大得多。Winkler假設(shè)[22]在連續(xù)分布載荷作用下，材料各點(diǎn)的位移與所受載荷成正比。Conway和Lee[23]根據(jù)Winkler假定，將橡膠內(nèi)襯的承壓作用等效為無窮多個(gè)彈簧的并聯(lián)，固定在剛性外殼與水膜之間，各節(jié)點(diǎn)的無量綱彈性變形可以表示為

式中：ΔH為內(nèi)襯變形量，u為軸頸轉(zhuǎn)速，μ為水的動(dòng)力粘度，R為軸頸半徑，l為橡膠內(nèi)部襯層的厚度，E為橡膠材料的彈性模量，c為軸承絕對(duì)間隙，v為泊松比。

橡膠軸承的水膜厚度包括軸承與轉(zhuǎn)子間的間隙和橡膠材料的變形。以軸承間隙c為特征長度進(jìn)行無量綱化，可得到全圓周式水潤滑軸承的無量綱水膜厚度為

式中：H為橡膠軸承的無量綱水膜厚度；ε為偏心率。

下面采用有限差分法求解（1）式。圖2為求解域網(wǎng)格圖，其中沿軸向和周向?qū)ⅵ撕挺葎澐譃閚和m個(gè)等距區(qū)間。

圖2 網(wǎng)格劃分示意圖Fig.2 Schematic diagram of meshing

采用中心差分格式，則（1）式可表示為

考慮自然破裂邊界條件，由（5）式可得各節(jié)點(diǎn)的Pi，j，再通過積分可求得水膜力FRx、FRy。

2 水潤滑橡膠軸承的動(dòng)態(tài)特性

2.1 水潤滑橡膠軸承的阻尼和剛度

當(dāng)軸承在工作中受到微小擾動(dòng)時(shí)，水膜通常被視為無窮多個(gè)并聯(lián)的彈簧與阻尼，其線性化的動(dòng)力學(xué)特性即為水膜的剛度和阻尼。圖3為水潤滑橡膠軸承的受力圖，其中FRx、FRy分別表示水膜力在x、y坐標(biāo)軸上的分力。

水膜剛度與阻尼為[24]

2.2 差分法求解動(dòng)態(tài)特性系數(shù)

同理可得Kxy、Kyy、Cxx、Cxy、Cyx和Cyy，其中無量綱位移擾動(dòng)量和速度擾動(dòng)量均取0.01。

圖3 水潤滑橡膠軸承彈性變形示意圖Fig.3 Schematic diagram of water lubrication rubber bearing’s elastic deformation

2.3 數(shù)值分析

下面主要研究偏心率、長徑比、間隙比及轉(zhuǎn)速對(duì)動(dòng)特性系數(shù)的影響，其中軸承的參數(shù)如表1所示。

表1 軸承參數(shù)Tab.1 Bearing parameters

（1）偏心率對(duì)動(dòng)力特性系數(shù)旳影響

圖4 無量綱動(dòng)態(tài)特性系數(shù)隨偏心率的變化：（a）阻尼系數(shù)Cxx，Cxy，Cyx，Cyy；（b）剛度系數(shù)Kxx，Kxy，Kyx，KyyFig.4 A curve on dimensionless coefficient with the change of eccentricity:(a)C;(b)K

由上圖可知，偏心率對(duì)軸承的動(dòng)態(tài)特性系數(shù)Kyy和Cyy的影響較大。隨著偏心率的增大，Kxx和Cxx剛度和阻尼系數(shù)隨之增大。在小偏心率下，交叉阻尼數(shù)值幾乎相等，而在大偏心率下，交叉阻尼數(shù)值不等。

（2）間隙比對(duì)動(dòng)力特性系數(shù)旳影響

圖5為水潤滑橡膠軸承的動(dòng)特性系數(shù)隨間隙比的變化情況。由圖可知，改變軸承與軸頸間的間隙比，無量綱動(dòng)態(tài)特性系數(shù)變化顯著。在小偏心率下，由于軸頸與軸承間的不完全接觸，間隙比對(duì)無量綱系數(shù)的影響較小；然而在大偏心率下，由于軸頸與軸承橡膠襯層的接觸面積較大，無量綱系數(shù)變化顯著，隨著偏心率的增加而顯著增大。

圖5 無量綱動(dòng)態(tài)特性系數(shù)隨間隙比的變化：（a）阻尼系數(shù)Cyy；（b）剛度系數(shù)KyyFig.5 A curve on dimensionless coefficient with the change of relative clearance ratio:(a)Cyy;(b)Kyy

圖6 無量綱動(dòng)態(tài)特性系數(shù)隨轉(zhuǎn)速的變化：（a）Cyy；（b）KyyFig.6 A curve on dimensionless coefficient with the change of speed:(a)Cyy;(b)Kyy

圖7 無量綱動(dòng)態(tài)特性系數(shù)隨長徑比的變化：（a）Cyy；（b）KyyFig.7 A curve on dimensionless coefficient with the change of length-diameter ratio:(a)Cyy;(b)Kyy

（3）轉(zhuǎn)速對(duì)動(dòng)態(tài)特性系數(shù)旳影響

圖6為水潤滑橡膠軸承的動(dòng)特性系數(shù)隨轉(zhuǎn)速的變化情況。從中可見，轉(zhuǎn)速變化對(duì)軸承的阻尼和剛度系數(shù)影響較大。在小偏心率下，轉(zhuǎn)速對(duì)無量綱的剛度系數(shù)與阻尼系數(shù)的影響較??；然而在大偏心率下，無量綱剛度系數(shù)與阻尼系數(shù)隨著轉(zhuǎn)速增加而顯著增大。

（4）長徑比對(duì)動(dòng)力特性系數(shù)旳影響

圖7為水潤滑橡膠軸承的動(dòng)特性系數(shù)隨長徑比的變化情況。從中可知，隨著長徑比的增大，無量綱動(dòng)態(tài)特性系數(shù)增大。在小偏心率下，阻尼和剛度系數(shù)增加的幅度較小，而在大偏心率下，則幅度增加較大。

3 水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)的動(dòng)態(tài)特性

下面分析轉(zhuǎn)子—軸承系統(tǒng)的動(dòng)力學(xué)特性，圖8為某大型船用橡膠軸承支承的螺旋槳軸—艉軸耦合系統(tǒng)，該系統(tǒng)由一個(gè)螺旋槳、三個(gè)軸承、二個(gè)轉(zhuǎn)子和一個(gè)聯(lián)軸器組成。

圖8 水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)模型Fig.8 A model of rotor system guided by water lubricated rubber bearings

轉(zhuǎn)子總長度超過14 m；螺旋槳質(zhì)量超過14 t，與轉(zhuǎn)子的自重合計(jì)近29 t。由于軸承約束反力作用點(diǎn)位置對(duì)計(jì)算結(jié)果影響較大，按照船舶水潤滑尾管橡膠軸承的設(shè)計(jì)規(guī)范[25]，螺旋槳軸承的支反力作用位置可以取在l1/5～l1/2軸承長度之間，本文取l1/4符合要求，亦符合美國軍用推進(jìn)軸系設(shè)計(jì)規(guī)范[26]。艉軸后軸承和艉軸前軸承約束反力的作用點(diǎn)分別為各自軸承中點(diǎn)處。

3.1 運(yùn)動(dòng)方程

下面討論水潤滑橡膠軸承支承的轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)模型。其中橡膠軸承利用8參數(shù)線性化的阻尼與剛度?；?，轉(zhuǎn)子采用集總參數(shù)法[27]離散。

（1）系統(tǒng)動(dòng)能

將轉(zhuǎn)子和圓盤結(jié)構(gòu)劃分為n個(gè)單元，系統(tǒng)的動(dòng)能[2]

忽略高階小量后，系統(tǒng)的動(dòng)能可以寫成為

式中：mi（i=1，2，…，n）是第i個(gè)單元的集中質(zhì)量，Vic是質(zhì)心速度，ωi是繞質(zhì)心的角速度，Gi是繞質(zhì)心的動(dòng)量矩，Jij（i=1，2，…，n；j=d，z）則分別是第i個(gè)單元的剖面轉(zhuǎn)動(dòng)慣量和繞z軸的轉(zhuǎn)動(dòng)慣量，xi、yi、φi、ψi是第i個(gè)結(jié)點(diǎn)的廣義位移，Ω為轉(zhuǎn)子的轉(zhuǎn)動(dòng)角速度。

（2）系統(tǒng)勢能

第i個(gè)梁單元的變形如圖9所示，其中Si-1，Mi-1，Qi-1，Ni-1，Si，Mi，Qi，Ni分別是作用在第i個(gè)單元上的剪力和彎矩。根據(jù)平衡方程，則在x方向可以得到

式中：E是轉(zhuǎn)子的彈性模量，I是截面的慣性矩，l是單元的長度。

圖9 轉(zhuǎn)子的第i個(gè)單元Fig.9 The ith element of the rotor

第i個(gè)單元在x方向的橫向勢能為

（3）運(yùn)動(dòng)微分方程

根據(jù)Lagrange方程

式中：廣義坐標(biāo)q= ｛…，xi，yi，zi，φi，ψi，…｝T，則系統(tǒng)的運(yùn)動(dòng)方程為

式中：系統(tǒng)的質(zhì)量矩陣M和剛度矩陣K是對(duì)稱的，而陀螺矩陣G則是反對(duì)稱的，C為阻尼矩陣，Q為廣義力矢量。一般而言，在廣義力矩陣Q中包含了各種外激勵(lì)力，例如不平衡力、扭轉(zhuǎn)激振力矩等，而質(zhì)量矩陣M和剛度矩陣K則主要與系統(tǒng)的物理參數(shù)有關(guān)，軸承的阻尼矩陣C和剛度矩陣K為

式中：cxx、cxy、cyx、cyy、kxx、kxy、kyx和kyy分別是水潤滑橡膠軸承的阻尼和剛度系數(shù)，這些系數(shù)主要與軸承的靜載荷和轉(zhuǎn)子的轉(zhuǎn)速有關(guān)。

3.2 動(dòng)態(tài)特性

下面以圖8所示的水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)為例進(jìn)行分析，其中的相關(guān)參數(shù)如表2所示。在下面的計(jì)算中，共劃分了28個(gè)梁單元，其中如圖8所示在軸承相對(duì)應(yīng)的位置設(shè)置了三個(gè)水潤滑橡膠軸承單元。

表2 轉(zhuǎn)子系統(tǒng)的結(jié)構(gòu)參數(shù)Tab.2 The structural parameters of rotor-bearing system

（1）臨界轉(zhuǎn)速及其振型

當(dāng)（14）式中的廣義力Q=0時(shí)，系統(tǒng)運(yùn)動(dòng)為自由振動(dòng)，相應(yīng)的特征值問題可以采用廣義逆迭代法[27]求解。通過對(duì)方程（15）進(jìn)行求解獲得其特征值，進(jìn)而得到轉(zhuǎn)子系統(tǒng)在剛性支承和彈性支承下的坎貝爾曲線，分別見圖10和圖11，對(duì)應(yīng)臨界轉(zhuǎn)速如表3所示。

圖10 剛性支承下的坎貝爾曲線Fig.10 Campbell curve under rigid support

圖11 彈性支承下的坎貝爾曲線Fig.11 Campbell curve under elastic support

表3 臨界轉(zhuǎn)速Tab.3 The critical speed

從表3可知，對(duì)于同一轉(zhuǎn)子而言，彈性支承相對(duì)于剛性支承來說，其支承方式變得更加柔性，因而其總剛度K下降，同時(shí)導(dǎo)致其固有頻率下降，從而導(dǎo)致其臨界轉(zhuǎn)速降低。圖12為彈性支承下各階臨界轉(zhuǎn)速對(duì)應(yīng)的振型。

圖12 彈性支承下系統(tǒng)的各階振型：（a）第一階；（b）第二階；（c）第三階；（d）第四階；（e）第五階；（f）第六階Fig.12 The modes of system under the elastic:(a)The first mode;(b)The second mode;(c)The third mode; (d)The fourth mode;(e)The five mode;(f)The sixth mode

由圖12可知，彈性支承下一階與二階臨界轉(zhuǎn)速對(duì)應(yīng)的模態(tài)主要以螺旋槳的振動(dòng)為主，三階—六階模態(tài)主要以軸段振動(dòng)為主。

（2）失穩(wěn)轉(zhuǎn)速

設(shè)系統(tǒng)方程的解具有一般形式x=Aeλt，系統(tǒng)穩(wěn)定性由其特性值λ=-u+iv決定，取決于其在復(fù)平面上的實(shí)際分布狀況。當(dāng)u＜0時(shí)，此時(shí)λ具有正實(shí)部，系統(tǒng)處于不穩(wěn)定狀態(tài)，即失穩(wěn)，即把u中出現(xiàn)u=0時(shí)所對(duì)應(yīng)轉(zhuǎn)子的轉(zhuǎn)速定義為失穩(wěn)轉(zhuǎn)速。

通過計(jì)算，發(fā)現(xiàn)轉(zhuǎn)子系統(tǒng)的失穩(wěn)轉(zhuǎn)速為1 380 r/min左右，經(jīng)分析可知系統(tǒng)失穩(wěn)轉(zhuǎn)速較低主要是由于螺旋槳集中質(zhì)量及轉(zhuǎn)子的型號(hào)相對(duì)較大而造成的，然而對(duì)于軸承直徑較大的大型軸承來說，這樣的速度已經(jīng)是相當(dāng)可觀的。

（3）橡膠襯層厚度的影響

橡膠軸承的支承剛度是一個(gè)很重要的參數(shù)，下面重點(diǎn)研究橡膠襯層的彈性對(duì)轉(zhuǎn)子系統(tǒng)臨界轉(zhuǎn)速與失穩(wěn)轉(zhuǎn)速的影響，結(jié)果如表4-5所示。

表4 橡膠襯層對(duì)臨界轉(zhuǎn)速的影響Tab.4 The impact of the support stiffness on the critical speed(r/min)

表5 橡膠襯層對(duì)失穩(wěn)轉(zhuǎn)速的影響Tab.5 The impact of support stiffness on threshold speed of stability(r/min)

由表4可知，水潤滑橡膠軸承的襯層厚度對(duì)水潤滑橡膠軸承支承的轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速影響較小。隨著襯層厚度的增加，系統(tǒng)的臨界轉(zhuǎn)速降低。其中第一階臨界轉(zhuǎn)速與第四階臨界轉(zhuǎn)速變化比較明顯；其它的各臨界轉(zhuǎn)速變化不大。

由表5可知，水潤滑橡膠軸承的襯層厚度對(duì)水潤滑橡膠軸承支承的轉(zhuǎn)子系統(tǒng)的失穩(wěn)轉(zhuǎn)速影響較小，但隨著襯層厚度的增加，系統(tǒng)的失穩(wěn)轉(zhuǎn)速有增大的趨勢。

4 結(jié) 論

本文主要研究了水潤滑橡膠軸承的動(dòng)態(tài)特性及水潤滑橡膠軸承支承下轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速及穩(wěn)定性。首先基于Winkler假定采用有限差分法分析了水潤滑橡膠軸承的剛度與阻尼；然后根據(jù)拉格朗日方程建立了水潤滑橡膠軸承—轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)模型，并針對(duì)某大型船舶推進(jìn)軸系研究了系統(tǒng)的動(dòng)態(tài)特性。結(jié)果表明水潤滑橡膠軸承的偏心率、間隙比、轉(zhuǎn)速及長徑比對(duì)軸承的阻尼和剛度影響較大；軸承的水彈潤滑特性對(duì)系統(tǒng)的動(dòng)態(tài)特性影響顯著，而水潤滑橡膠軸承的彈性變化對(duì)系統(tǒng)的失穩(wěn)轉(zhuǎn)速影響不大。上述研究可以為水潤滑橡膠軸承支承轉(zhuǎn)子系統(tǒng)在各種激勵(lì)作用下的動(dòng)力學(xué)分析提供理論準(zhǔn)備。

[1]彭晉民,朱志宏,王家序.水潤滑軸承的研究現(xiàn)狀及展望[J].潤滑與密封,2004(5):124-126. Peng Jinming,Zhu Zhihong,Wang Jiaxu.Water lubricated bearings present and future[J].Lubrication Engineering,2004 (5):124-126.

[2]李 方,帥長庚,何 琳,楊 雪,李 明.橡膠軸承耦合轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)研究[J].噪聲與振動(dòng)控制,2011,31(7):37-41. Li Fang,Shuai Changgeng,He Lin,et al.Dynamic analysis of rotor system coupled by rubber bearings[J].Noise and Vibration Control,2011,31(7):37-41.

[3]岑少起,陳 瀚,唐照民.長重載徑向軸承潤滑力學(xué)問題數(shù)值解[J].西安交通大學(xué)學(xué)報(bào),1982,16(2):119-126. Cen Shaoqi,Chen Han,Tang Zhaomin.The numerical solution of lubrication-mechanical equations for heavy-loaded jounal bearings[J].Jounal of Xi’an Jiaotong University,1982,16(2):119-126.

[4]陳伯賢,裘祖干.動(dòng)載滑動(dòng)軸承潤滑狀況的預(yù)測—彈性流體動(dòng)力潤滑理論的應(yīng)用[J].內(nèi)燃機(jī)學(xué)報(bào),1983,1(4):47-58. Chen Boxian,Qiuzugan.A calculation for dynamically loaded journal bearings-An application of elastohydrodynamic lubrication theory[J].Transaction of CSICE,1983,1(4):47-58.

[5]王優(yōu)強(qiáng),宋 玲,李鴻琦,等.水潤滑橡膠軸承設(shè)計(jì)研究[J].潤滑與密封,2003(2):21-22. Wang Youqiang,Song Ling,Li Hongqi,et al.Design of the water-lubricated rubber bearings[J].Lubrication Engineering, 2003(2):21-22.

[6]姚世衛(wèi),胡宗成,馬 斌,等.橡膠軸承研究進(jìn)展及在艦艇上的應(yīng)用分析[J].艦船科學(xué)技術(shù),2005,(27):27-30. Yao Shiwei,Hu Zongcheng,Ma Bin,et al.The new development of rubber bearing and its application in warships[J]. Ship Science and Technology,2005,(27):27-30.

[7]賈小俊,范世東.考慮軸頸傾斜的徑向滑動(dòng)軸承動(dòng)態(tài)特性研究[J].船海工程,2008,37(5):54-57. Jia Xiaojun,Fan Shidong.Dynamic lubrication characteristics analysis of journal bearing considering misalignment[J]. Ship and Ocean Engineering,2008,37(5):54-57.

[8]de Kraker A,van Ostayen R A J,Rixen D J.Calculation of stribeck curves for(water)lubricated journal bearings[J].Tribology International,2007,40:459-469.

[9]Wang Youqiang,Wang Meishu.Experimental study on the tribological properties of water-lubricated rubber bearings[J]. Advanced Materials Research,2011,291-294:1454-1457.

[10]Wang Youqiang.Numerical simulation of water-lubricated rubber/chromium coating pair test[J].Advanced Materials Research,2012,502:351-354.

[11]Zhao Hui,Wang Youqiang.The elastohydrodynamic lubrication analysis of starved water-lubricated rubber bearing[J]. Applied Mechanics and Materials,2012,182-183:1529-1532.

[12]Wang Jiaxu,Hua Xijin,Xiao Ke,Zhu Juanjuan.Two-dimensional study of lubrication mechanism of water-lubricated rubber alloy bearing[J].International Journal of Computer Applications in Technology,2011,40(1):98-106.

[13]周廣武,王家序,王戰(zhàn)江,等.多溝槽水潤滑橡膠合金軸承潤滑特性研究[J].摩擦學(xué)學(xué)報(bào),2013,33(6):631-637. Zhou Guangwu,Wang Jiaxu,Wang Zhanjiang,et al.Analysis of multi-grooves water lubricated rubber alloy bearing considering the elastohydrodynamic lubrication[J].Tribology,2013,33(6):631-637.

[14]王 娟,王 雋,楊 俊.水潤滑橡膠軸承板條設(shè)計(jì)參數(shù)分析[J].船舶科學(xué)技術(shù),2013,35(3):126-128. Wang Juan,Wang Jun,Yang Jun.Design parameter research of water-lubricated rubber bearings[J].Ship Science and Technology,2013,35(3):126-128.

[15]周廣武,王家序,周 憶,等.水潤滑橡膠軸承動(dòng)態(tài)剛度和阻尼測試方法研究[J].四川大學(xué)學(xué)報(bào)(工程科學(xué)版),2014, 46(3):194-198. Zhou Guangwu,Wang Jiaxu,Zhou Yi,et al.Experimental study on the dynamic stiffness and damping of water lubricated rubber bearings[J].Jounal of Sichuan University(Engineering Science Edition),2014,46(3):194-198.

[16]李文鋒,王家序,丁行武,程美玲.橡膠層厚度對(duì)水潤滑軸承潤滑性能的影響[J].潤滑與密封,2014,39(1):65-79. Li Wenfeng,Wang Jiaxu,Ding Xingwu,et al.Effect of thickness of rubber layer on lubrication performance of water-lubricated bearings[J].Lubrication Engineering,2014,39(1):65-79.

[17]趙 娟.高速滑動(dòng)軸承—轉(zhuǎn)子系統(tǒng)的動(dòng)態(tài)特性[D].哈爾濱:哈爾濱工業(yè)大學(xué),2010. Zhao Juan.Dynamic characteristics analysis of high speed journal bearing and rotor system[D].Harbin:Harbin Institute of Technology,2010.

[18]劉鎮(zhèn)星.轉(zhuǎn)子—水潤滑橡膠軸承系統(tǒng)動(dòng)力學(xué)特性研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2013. Liu Zhenxing.Research on dynamics characteristic of a rotor water lubricated rubber bearing system[D].Harbin:Harbin Institute of Technology,2013.

[19]張廣輝.高速動(dòng)靜壓混合氣體軸承轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特性研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2010. Zhang Guanghui.Research of dynamic characteristics for high speed hybrid gas bearing rotor system[D].Harbin:Harbin Institute of Technology,2010.

[20]蔣小平,施衛(wèi)東,李 偉,等.泵水潤滑軸承轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)特性研究[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(1):63-66. Jiang Xiaoping,Shi Weidong,Li Wei,et al.Rotor dynamic characteristics of water lubricated bearing for pumps[J].Transactions of the Chinese Society for Agricultural Machinery,2013,44(1):63-66.

[21]Liu Shibing,Yang Bingen.A new model of water-lubricated rubber bearings for vibration analysis of flexible multistage rotor systems[J].Journal of Sound and Vibration,2015,349:230-258.

[22]Winkler S T.Strength of materials:Advanced theory and problems[M].D.Van Nostrand Company,1956.

[23]Conway H D,Lee H C.The analysis of the lubrication of a flexible journal bearing[J].Journal of Lubrication Technology, 1975(97):599-604.

[24]鐘一諤,何衍宗,等.轉(zhuǎn)子動(dòng)力學(xué)[M].北京:清華大學(xué)出版社,1987:43-44.

[25]楊和廷,唐育民.船舶水潤滑尾管橡膠軸承的設(shè)計(jì)[J].武漢造船,2000,2:19-22. Yang Heting,Tang Yumin.Design of marine stern rubber bearings[J].Wuhan Shipbuilding,2000,2:19-22.

[26]U.S.Department of defense military handbook[M].Design methods for naval shipboard systems,Section 243-1,Part 1, Propulsion shafting,1994.

[27]虞 烈,劉 恒.軸承—轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)[M].西安:西安交通大學(xué)出版社,2001:129-154.

Study on the dynamic characteristics of rotor system supported on water-lubricated rubber bearings

SHI Yong-feng1,LI Ming1,LIU Gang1,HE Lin2
(1.Department of Mechanics,Xi’an University of Science and Technology,Xi’an 710054,China; 2.Institute of Noise and Vibration,Naval University of Engineering,Wuhan 430033,China)

The dynamic characteristics of water-lubricated rubber bearing and dynamic characteristics of rotor system supported on water-lubricated rubber bearings are studied.The mathematical model of waterlubricated rubber bearing is built based on the theory of elastohydrodynamic lubrication.Its dynamic characteristics are analyzed according to the load increment method and finite difference method.Then,the differential equation of rotor system coupled with the water-lubricated rubber bearing is deduced based on Lumped parameter method and the theory of Analytical Mechanics.Finally,the stability and dynamic characteristics of rotor system supported on water-lubricated rubber bearings are analyzed by using the generalized inverse iteration method.The critical speeds and its corresponding modes,and the threshold speed of stability,and so on are acquired.The results show that the parameters,such as the eccentricity,the length-diameter ratio,the relative clearance ratio and the speed,have a greater influence on dynamic characteristics of water-lubricated rubber bearings;the elastohydrodynamic lubrication characteristics of rubber bearing has a large effect on the critical speed of system,but the elastic deformation of rubber bearings has a little influence on the stability of system.

water-lubricated rubber bearing;rotor system;dynamic characteristics;the critical speed;the stability

U664.81+3

：Adoi:10.3969/j.issn.1007-7294.2017.05.008

1007-7294（2017）05-0584-11

2016-12-16

國家自然科學(xué)基金資助項(xiàng)目（11372245）；陜西省自然科學(xué)基金資助項(xiàng)目（2014JM1015）

史永峰（1988-），男，碩士研究生，E-mail：S634821083@126.com；李 明（1963-），男，教授，博士生導(dǎo)師。