張 偉，吳偉亮

(上海交通大學(xué) 機(jī)械與動力工程學(xué)院，上海200240)

0 引 言

軸系回轉(zhuǎn)振動的實(shí)質(zhì)是轉(zhuǎn)軸的進(jìn)動，軸系發(fā)生回轉(zhuǎn)運(yùn)動時一方面圍繞其幾何軸線旋轉(zhuǎn)，另一方面彎曲的幾何軸線圍繞其支撐中心線旋轉(zhuǎn)。現(xiàn)代船舶朝著越來越大型化的方向發(fā)展，使得船體剛度下降，導(dǎo)致橫向振動頻率下降［1］。對于氣墊船軸系，作為支撐體的船體浮箱剛度較弱，使得軸系的橫向振動響應(yīng)更為強(qiáng)烈，所以對氣墊船軸系回旋振動的研究意義更大。

國內(nèi)外學(xué)者對軸系回旋振動展開了一些研究。陳之炎等［2－3］對回旋振動的機(jī)理進(jìn)行了深入地探討，并對回旋振動計算方法進(jìn)行了研究。王磊等［4］研究分析了慣性力矩、應(yīng)力剛化效應(yīng)及旋轉(zhuǎn)軟化效應(yīng)等對軸系回旋振動的影響。此外，王傳溥、Y Hori和Al-Bedoor B O 等［5－7］在回旋振動研究方面也做了大量工作。

本文以某氣墊船墊升軸系為對象，考慮到其與傳統(tǒng)船舶軸系的不同，需要對其回旋特性進(jìn)行計算分析，通過計算了解其振動特性，為軸系安全運(yùn)行提供指導(dǎo)。

1 力學(xué)模型

為便于分析，本文以某氣墊船墊升軸系為例，建立有限元計算模型，進(jìn)行有關(guān)共性問題的研究，軸系主要參數(shù)如表1所示。在有限元計算軟件Ansys中建立模型，計算模型如圖1所示。

圖1 墊升軸系有限元模型Fig.1 The finite element model of the shafting

考慮到墊升軸系結(jié)構(gòu)的復(fù)雜性，在進(jìn)行船舶軸系振動特性計算時一般都要將其簡化，以便于求解，針對這一軸系，作以下適當(dāng)簡化，計算結(jié)果顯示可以反映實(shí)際振動情況:

1)選取梁單元模擬軸系軸段，彈簧單元模擬彈性支承系統(tǒng)的軸承。彈簧單元的一端與軸系的節(jié)點(diǎn)相連接，另一端設(shè)置為固定端，這樣可以沿某一方向傳遞力或力矩。

2)風(fēng)機(jī)、齒輪在有限元計算中都簡化為質(zhì)量單元，并考慮其轉(zhuǎn)動慣量。軸段端部聯(lián)軸節(jié)帶有附加質(zhì)量也由質(zhì)量單元模擬。

3)由于法蘭連接的彎曲剛度比軸的剛度大得多，直接將連接法蘭作為軸段元件，忽略其彎曲變形。

4)齒輪耦合單元采用矩陣單元模擬，其平均嚙合剛度取為2 ×106N/m［8］。

2 計算分析與討論

2.1 墊升軸系回旋振動計算

軸承的支撐剛度是軸系振動特性的重要影響因素，它與很多因素有關(guān)，如軸承結(jié)構(gòu)、軸承材料、軸承間隙、油膜等，所以對支撐剛度的確定是一件相當(dāng)復(fù)雜的工作，一般很難單純由計算獲得。根據(jù)國內(nèi)外實(shí)船數(shù)據(jù)，軸承支撐剛度在108～5 ×109N/m范圍內(nèi)。因此在本文中，采用考慮不同軸承支撐剛度的方法進(jìn)行計算，給出轉(zhuǎn)速范圍的參考值。

由于墊升軸系減速器速比為1.53，墊升軸系水平分軸與豎直分軸轉(zhuǎn)速不同，故分別列出其回旋振動計算結(jié)果。

當(dāng)軸承剛度為1 ×108N/m，其坎貝爾圖如圖2所示，可以得到一次回旋轉(zhuǎn)速及葉片次回旋轉(zhuǎn)速。通過改變軸承剛度，計算得到不同的回旋轉(zhuǎn)速，如表2和表3所示。

考慮到風(fēng)機(jī)安裝區(qū)域內(nèi)船體剛性較弱，船體剛

圖2 軸承剛度1 ×108 N/m 下墊升風(fēng)機(jī)豎直分軸坎貝爾圖Fig.2 Campbell diagram of the vertical shaft under the bearing stiffness of 1 ×108 N/m

表2 未考慮船體剛度時墊升軸系豎直分軸回旋振動計算Tab.2 Lateral vibration of the vertical shaft taking no consideration of the hull stiffness

表3 墊升軸系水平分軸回旋振動計算Tab.3 Lateral vibration of the horizontal shaft

度在2 ×108～2 ×109N/m 內(nèi)取值，將該船體剛度串聯(lián)至風(fēng)機(jī)端的上述3個剛度上，得到總的安裝區(qū)域軸承剛度分布在6.7 ×107～6.7 ×108N/m 內(nèi)，在此范圍內(nèi)取一系列剛度值，計算結(jié)果如表4所示。

表4 考慮船體剛度時墊升軸系豎直分軸回旋振動計算Tab.4 Lateral vibration of the vertical shaft taking consideration of the hull stiffness

通過表2～表4 可看出，隨著軸承支撐剛度的增大，風(fēng)機(jī)軸系的共振轉(zhuǎn)速逐漸升高。分析表2 可看出，在未考慮船體剛度時，墊升軸系豎直分軸一次共振轉(zhuǎn)速在1 736.86～3 154.62 r/min 范圍內(nèi);一次正回旋共振轉(zhuǎn)速在1 809.66～3 503.13 r/min 范圍內(nèi);葉片次共振轉(zhuǎn)速在146.01～275.82 r/min 范圍內(nèi)。當(dāng)考慮船體剛度時，由于串聯(lián)船體剛度，總軸承剛度的下降，相應(yīng)共振轉(zhuǎn)速也隨之下降。表4 顯示墊升風(fēng)機(jī)豎直分軸一次共振轉(zhuǎn)速在1 471.37～3 026.24 r/min 范圍內(nèi)，葉片次共振轉(zhuǎn)速在123.72～261.27 r/min 范圍內(nèi)。墊升軸系水平分軸的一次共振轉(zhuǎn)速為3 200.70～5 252.05 r/min，一次正回旋轉(zhuǎn)速也在3 211.30 r/min 以上。

墊升軸系豎直分軸的額定轉(zhuǎn)速為850 r/min，水平分軸的額定轉(zhuǎn)速在1 300 r/min，根據(jù)規(guī)定，軸系在1.15 倍額定轉(zhuǎn)速下無一次回旋振動臨界轉(zhuǎn)速，一次葉片次臨界轉(zhuǎn)速不應(yīng)在0.8～1.2 倍額定轉(zhuǎn)速范圍內(nèi)，故墊升軸系回旋振動安全。

2.2 對比陀螺效應(yīng)影響

對于一個帶有質(zhì)量的圓盤，不考慮陀螺力矩時自由振動方程

當(dāng)考慮陀螺效應(yīng)時，系統(tǒng)的運(yùn)動方程中會出現(xiàn)一個反對稱的陀螺矩陣。一般情況下，靜止坐標(biāo)系下考慮陀螺效應(yīng)的轉(zhuǎn)子系統(tǒng)運(yùn)動方程式可寫為

式中:M，C，K 分別為質(zhì)量、阻尼、剛度矩陣;G為陀螺力矩。

對上述實(shí)例水平分軸進(jìn)行回旋振動有限元計算，針對考慮陀螺效應(yīng)與不考慮陀螺效應(yīng)在不同轉(zhuǎn)速下進(jìn)行計算，結(jié)果如表5～表7所示。

表5 轉(zhuǎn)速944 r/min 下陀螺效應(yīng)對固有頻率的影響Tab.5 The gyroscopic effect on natural frequency under 944 r/min

表6 轉(zhuǎn)速2 360 r/min 下陀螺效應(yīng)對固有頻率的影響Tab.6 The gyroscopic effect on natural frequency under 2 360 r/min

表7 轉(zhuǎn)速3 776 r/min 下陀螺效應(yīng)對固有頻率的影響Tab.7 The gyroscopic effect on natural frequency under 3 776 r/min

當(dāng)軸系發(fā)生正進(jìn)動時，此時陀螺力矩為正值，它使轉(zhuǎn)軸的變形變小，因而提高了軸系的固有頻率。反之，當(dāng)發(fā)生反進(jìn)動時，軸系的固有頻率下降。此外陀螺效應(yīng)與軸系的轉(zhuǎn)速有關(guān)，轉(zhuǎn)速越高，陀螺效應(yīng)表現(xiàn)的越顯著，隨著軸系轉(zhuǎn)速從944 r/min 到3 776 r/min，正進(jìn)動差值也從0.15% 上升到0.31%。這是因?yàn)橥勇萘匕缡綉T性力矩及牽連慣性力矩，二者都與軸系旋轉(zhuǎn)角速度有關(guān)，軸系旋轉(zhuǎn)角速度越大，陀螺效應(yīng)也就愈顯著［1］。

3 模型驗(yàn)證

文獻(xiàn)［9］通過建立氣墊船軸系縮比臺架，進(jìn)行回旋自由振動實(shí)驗(yàn)。針對氣墊船軸系試驗(yàn)臺架，通過Ansys 軟件，采用上述簡化方法對其建模計算，對比實(shí)驗(yàn)值與計算值，如表8所示，證明了使用有限元方法對軸系的簡化計算能較好地反映振動特性，同時也說明了上述分析合理。

表8 實(shí)驗(yàn)與計算結(jié)果對比Tab.8 Comparison between the experimental data and the numerical results

4 結(jié) 語

本文以氣墊船墊升軸系為對象，建立有限元計算模型，計算驗(yàn)證了陀螺效應(yīng)對軸系回旋振動的影響，當(dāng)陀螺力矩為正值，軸系的固有頻率提高，反之下降;并分別討論了在不考慮船體剛度及考慮船體剛度下軸系的共振轉(zhuǎn)速及葉片次共振轉(zhuǎn)速，對照工作轉(zhuǎn)速，給出其回旋振動安全的結(jié)論，通過與文獻(xiàn)對比，也驗(yàn)證了計算模型的準(zhǔn)確性，為氣墊船墊升軸系的安全運(yùn)行提供了指導(dǎo)。

［1］陳之炎.船舶推進(jìn)軸系振動［M］.上海:上海交通大學(xué)出版社，1987:172－233.CHEN Zhi-yan.The vibration of ship's propulsion shafting［M］.Shanghai:Shanghai Jiaotong University Press，1987:172－233.

［2］陳之炎，趙玫，駱振黃.船舶推進(jìn)軸系的回旋振動［J］.上海交通大學(xué)學(xué)報，1984(2):91－99.CHEN Zhi-yan，ZHAO Mei，LUO Zhen-huang.Whirling of ship's propulsion shafting［J］.Journal of Shanghai Jiaotong University，1984(2):91－99.

［3］陳之炎，駱振黃，熊四昌.船舶推進(jìn)軸系回旋振動計算方法研究［J］.中國造船，1988(3):33－43.CHEN Zhi-yan，LUO Zhen-huang，XIONG Si-chang.A study of the calculation method of transverse vibration for ship's shafting system［J］.Shipbuilding of China，1988(3):33－43.

［4］王磊，謝俊超，周瑞平.大型船舶推進(jìn)軸系回旋振動特性分析研究［J］.江蘇船舶，2010，27(1):14－17.WANG Lei，XIE Jun-chao，ZHOU Rui-ping.Research of the lateral vibration of super ship propulsion shafting［J］.Jiangsu Ship，2010，27(1):14－17.

［5］王傳溥，劉志剛，張洪田，等.船舶軸系橫向振動共振轉(zhuǎn)速的實(shí)驗(yàn)［J］.船舶工程，1995(4):23－24，59.WANG Chuan-pu，LIU Zhi-gang，ZHANG Hong-tian，et al.Experimental investigation of resonance speed for ship propeller shafting in lateral vibration［J］.Ship Engineering，1995(4):23－24，59.

［6］HORI Y，et al.Lateral vibrations of propeller shaft systems［M］.Bulletin of the Marine Engineering Society in Japan，1978.

［7］AL-BEDOOR B O.Modeling the coupled torsional and lateral vibrations of unbalanced rotors［J］.Computer Methods in Applied Mechanics and Engineering，2001，190(45):5999－6008.

［8］宋雪萍，劉樹英，聞邦椿.齒輪－轉(zhuǎn)子系統(tǒng)的振動特性分析［J］.機(jī)械科學(xué)與技術(shù)，2006，25(2):153－157.SONG Xue-ping，LIU Shu-ying，WEN Bang-chun.Vibration characteristics analysis of a gear rotor system ［J］.Mechanical Science and Technology，2006，25(2):153－157.

［9］李寧.具有彈性基礎(chǔ)的氣墊船細(xì)長軸系橫向振動特性研究［D］.哈爾濱:哈爾濱工程大學(xué)，2007.LI Ning.Research of lateral vibration for slightness shafting with elastic base on the hovercraft［D］.Harbin:Harbin Engineering University，2007.