基于流固耦合的船舶軸-槳耦合振動特性分析

李小軍，朱漢華，熊 維，吳繼東

(1. 武漢理工大學 能源與動力工程學院，湖北 武漢 430063；2. 武昌船舶重工集團有限公司，湖北 武漢 430060)

基于流固耦合的船舶軸-槳耦合振動特性分析

李小軍1，朱漢華1，熊 維2，吳繼東2

(1. 武漢理工大學 能源與動力工程學院，湖北 武漢 430063；2. 武昌船舶重工集團有限公司，湖北 武漢 430060)

以研究螺旋槳水動力和離心力對船舶軸-槳組合振動特性的影響為本文的研究目的，基于W o r kbench平臺，采用流固耦合有限元分析方法，進行船舶軸-槳組合模態分析。在CFX中計算螺旋槳敞水性能，并在Ansys中將螺旋槳葉面水壓力和離心力作為預應力分析軸槳組合振動的固有頻率和振型，比較軸系、螺旋槳單獨模型和軸-槳組合模型在固有頻率上的區別。計算結果表明，軸槳組合的固有頻率遠遠低于軸系和螺旋槳獨立模型的固有頻率；軸-槳旋轉產生的離心力對其固有頻率影響不大；螺旋槳在流場中產生的水壓力略微提高縱向振動固有頻率，但影響很小，在實際應用中可以忽略。

流固耦合；螺旋槳；軸系；Ansys；模態分析

0 引 言

船舶螺旋槳和軸系在實際工作時，流場中不均勻載荷會影響槳葉表面壓力分布和槳葉振動特性；同時流場中的載荷會通過螺旋槳葉面傳遞給軸系，引起軸系應力、應變以及振動特性的變化，因此研究軸-槳組合的流固耦合振動特性對探索船舶尾部振動和噪聲的產生原因非常有意義。

邊界元和有限元理論的出現以及計算機的發展，使得研究軸-槳組合的流固耦合振動數值方法得到了較快發展。W.Detter等[1]利用有限元軟件和流體計算軟件來建立固體域和流體域的計算模型，在2個模型間進行位移和壓力邊界條件的傳遞從而達到求解流固耦問題的目的。流固耦合模態分析在水輪機、軸流泵等旋轉機械中應用較廣，施衛東等[2]運用Ansys WORKBENCH和APDL，命令流耦合的方法，對潛水軸流泵的軸系轉動部件進行了考慮預應力情況下的濕模態數值模擬。針對螺旋槳和軸系的振動特性，國內外一些學者進行了相關研究，文獻[3]采用CFD方法對螺旋槳敞水特征進行模擬，提取螺旋槳表面的壓力波動作為外部載荷，再利用有限元方法對螺旋槳進行振動響應分析；文獻[4]采用Ansys軟件將螺旋槳簡化為圓盤，計算軸系的振動特性。

然而，將軸系和螺旋槳作為組合結構，考慮其相互耦合作用的流固耦合模態分析的文獻卻不多見。本文流固耦合的有限元分析方法，建立螺旋槳和軸系的組合模型，利用CFX對螺旋槳進行螺旋槳的敞水性能分析，將模擬結果壓力項數據通過流固耦合面傳遞給Ansys求解器，進行軸-槳結構的結構分析和模態分析，研究軸-槳結構在離心力和螺旋槳濕表面壓力作為預應力下的固有振動特性。

1 力學模型及計算方法

1.1 螺旋槳的不可壓縮粘性流體運動方程

螺旋槳在流體中旋轉，不涉及能量方程，其連續性方程和動量方程可表示為：

連續方程

本文中流體是水，不可壓縮，則

1.2 軸-槳流固耦合運動方程

考慮到流體作用，軸系和螺旋槳以恒定轉速繞軸線運轉的運動方程如下[5]：

式中：M為結構質量矩陣；C為結構阻尼矩陣；K為結構剛度矩陣；u為結構廣義位移矢量；f0為流體作用在結構濕表面上的壓力矢量，f為結構力矢量。

1.3 求解方法與過程

在Pro/E中建立螺旋槳幾何模型，并在CFX中采用RNGk～ε湍流模型和SIMPLE耦合算法，進行敞水性能分析。隨后在Ansys中建立船舶軸-槳組合模型，其中離心力通過設置旋轉速度進行模擬，螺旋槳葉面上濕表面壓力通過流固耦合面（Fluid Solid Interface）載入到Ansys Multiphysics中的螺旋槳葉面上，采用Lanczos法，進行預應力的模態分析。

2 應用分析

本文采用西江干線某船的軸系及其螺旋槳作為對象，應用流固耦合方法進行分析。螺旋槳為MAU型等螺距槳，其尺寸如表1所示；軸系包含一個尾軸和3個中間軸，包括前后尾軸承各1個和3個中間軸軸承，軸系尺寸如表2所示。螺旋槳和軸系的材料性能參數如下：密度7 850×103kg/m3，楊氏模量E=2.06 GPa，泊松比μ=0.3。軸-槳模型如圖1所示。

表 1 螺旋槳主要參數Tab. 1 Main parameters of the propeller

表 2 軸系主要尺寸參數Tab. 2 Main dimension parameters of the shaft

通過已知的螺旋槳葉二維型值點坐標，計算出螺旋槳的三維型值點，將其轉化為pts格式導入Pro/E軟件中，連成曲線，再生成空間曲面，通過邊界混合命令曲面混合生成螺旋槳葉實體，通過圓周陣列將其與軸轂合成一體，建立螺旋槳的三維模型[6]，最后與軸系進行合并。

2.1 流體域的建立及網格劃分

流體域采用分區網格劃分，計算域由固定域和旋轉域組成[7]。固定域為外部大圓柱體（直徑5D，長10D），采用結構化六面體網格劃分，旋轉域為內部小圓柱體（直徑1.2D，長1.5D），2個區采用CFX的GGI方式連接。對螺旋槳葉梢和葉根處進行加密，在螺旋槳槳葉外設置3層邊界層。進口邊界為速度進口，給定進速系數；出口采用壓力出口，壓力為常數，內部旋轉域和外部固定域之間采用interface面傳遞數據，網格關聯方式為GGI；螺旋槳和輪轂部分表面設置為無滑移、不可穿透壁面，外部邊界設置為開放邊界。旋轉域的運動形式為MRF，繞Z軸旋轉；固定域和螺旋槳靜止。螺旋槳周圍流場計算域如圖2所示，旋轉域非結構化網格如圖3所示。

2.2 軸-槳結構場的模型建立及網格劃分

螺旋槳與軸系的組合模型采用四面體非結構化劃分，網格劃分后的有限元模型如圖4所示。

3 計算結果分析

3.1 水動力計算

進速系數分別取0.3，0.4，0.5，0.6，0.7，由CFX數值模擬得到螺旋槳所受到的力和力矩值，通過換算[8]得到推力系數KT、轉矩系數KQ以及敞水效率η，如圖5所示。

由圖5可知隨著J的增大，推力和扭矩均降低，敞水效率η增大。當J=0.3時，推力與扭矩較大，表3列出了此時螺旋槳所受力與力矩，圖6給出了螺旋槳葉背和葉面的壓力分布圖。

表 3 螺旋槳所受的力和力矩Tab. 3 Force and torque of propeller

由表3和圖6可知，螺旋槳在軸向受到了較大的推力和扭矩，其中槳葉的最大壓力主要出現在旋轉正方向的槳葉邊緣；X和Y方向由于葉片的均勻分布而受力動平衡。

3.2 結構振動分析

螺旋槳在進行流固耦合時，以槳葉片處所受載荷為主，通過流固耦合面，將流體域葉片的表面壓力施加到對應固體域葉片上，如圖7所示。軸承采用BEARING進行Body-Ground連接，考慮彈性支撐，水平剛度和垂直剛度均為4.6×109N/m，不考慮交叉剛度。

3.2.1 旋轉預應力和水壓力預應力對軸槳振動影響

當J=0.3時，旋轉速度分別為0 r/min和327.3 r/min（額定轉速），進行有預應力和無預應力的模態分析，其中轉速為327.3 r/min時，軸-槳結構的前6階振型如圖8所示；軸-槳組合結構前10階固有頻率，如表4和表5所示。

從圖8可以看出，軸-槳的前幾階振型均出現在靠近螺旋槳處，這是由于近推力軸承端被設置為全約束，而螺旋槳重量和慣量較大，因此其1階振動時螺旋槳處振幅較大，這與實際相符。由于葉片柔度較大，因此多表現為葉片劇烈振動。

圖8中，1階振型為扭轉1階振動，表現為螺旋槳繞軸線扭轉；2階振型為橫向1階（水平）振動，表現為螺旋槳處水平彎曲；3階振型為橫向1階（垂直）振動，表現為螺旋槳處垂直彎曲；4階陣型為縱向1階振動，表現為1，3槳葉與2，4槳葉沿著軸向反向擺動；5階振型為2階縱向振動，主要表現為螺旋槳處軸向運動；6階振型為橫向2階（水平）振動，主要表現為螺旋槳和2#軸承處水平彎曲。

由表4和表5可以看出，當考慮螺旋槳對軸系的預應力時，軸-槳的縱向固有頻率會略微提高，增長幅度很小；橫向和扭轉固有頻率不僅沒有提高，反而略微下降。總的來說，螺旋槳在均勻來流中產生的水動力對于軸-槳結構的振動模態影響較小。

3.2.2 獨立模型與組合模型對比分析

為探究螺旋槳與軸系之間耦合振動的關系，對約束軸轂近軸端面的螺旋槳進行有預應力的模態分析，對不含螺旋槳的軸系進行模態分析，并將兩者與軸槳組合模型預應力的模態分析進行對比，其固有頻率關系如圖9所示。

分析圖9可知，由于是四葉槳，螺旋槳的振動多以4階為一個周期；由于螺旋槳的存在相當于增大了軸系的懸臂長度，降低了其相對剛度，因此軸-槳組合固有頻率遠遠低于軸、槳獨立模型的固有頻率，且差距隨著階次的增加而擴大。

表 4 預應力對軸-槳組合結構固有頻率的影響對比（轉速為0 r/min）Tab. 4 Effect of pre-stress on natural frequency of shaftpropeller（rotational speed=0 r/min）

表 5 預應力對軸-槳組合結構固有頻率的影響對比（轉速為327.3 r/min）Tab. 5 Effect of pre-stress on natural frequency of shaftpropeller（rotational speed =327.273 r/min）

4 結 語

本文在Pro/E中建立了某船舶軸系與螺旋槳的組合模型，在Workbench中完成了基于流固耦合的軸-槳組合模型振動特性分析。通過CFX對螺旋槳進行單項流固耦合計算，通過流固耦合面將槳葉上的壓力載荷施加Ansys中的軸槳組合模型上，進行了軸-槳組合在預應力下的模態分析，并比較了軸系、螺旋槳單獨模型和軸-槳組合模型在固有頻率上的區別，主要結論如下：

1）離心力對軸槳組合的固有頻率影響不大；螺旋槳在流場中所受的水壓力略微提高了軸槳組合的縱向剛度，進而略提高縱向固有頻率，但影響不大，在實際應用中可以忽略。

2）由于螺旋槳的懸臂效應，軸-槳組合的固有頻率遠遠低于軸系和螺旋槳獨立模型的固有頻率。

3）螺旋槳在水中工作，其振動時會帶動周圍流體一起振動，同時流體會對螺旋槳產生阻尼效應，因此其濕模態可能與干模態有所差距，需進一步分析。

[1]W DETTER, D PERIC. A computational framework for fluid–structure interaction: Finite element formulation and applications[J]. Comput Methods Appl. Mech. Engrg, 2006, 195(41–43): 5754–5779.

[2]施衛東, 郭艷磊, 張德勝, 等. 大型潛水軸流泵轉子部件濕模態數值模擬[J]. 農業工程學報, 2013, 29(24): 72–78. SHI Wei-dong, GUO Yan-lei, ZHANG De-sheng, et al. Numericalsimulation on modal of large submersible axial-flow pumprotor[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(24): 72–78.

[3]吳思遠, 黎勝. 船用螺旋槳葉片振動輻射噪聲數值分析[J]. 振動與沖擊, 2014, 33(12): 207–210. WU Si-yuan, LI Sheng. Numerical analysis for noise induced by vibration of propeller blades [J]. Journal of Vibration and Shock, 2014, 33(12): 207–210.

[4]朱軍超. 大型船舶支撐狀態對軸系振動的影響研究[D]. 武漢:武漢理工大學, 2012. ZHU Jun-chao. Reasearch on shafting vibration of large vessel based on the bearing support state[D]. Wu Han: Wuhan University of Technology, 2012.

[5]楊光, 熊鷹, 黃政. 復合材料螺旋槳雙向流固耦合計算[J]. 艦船科學技術, 2015, 37(10): 16–20. YANG Guang, XIONG Ying, HUANG Zheng. Computation of composite propeller’s two-way fluid-structure coupling[J]. Ship Science and Technology, 2015, 37(10): 16–20.

[6]陳艷鋒, 吳新躍. 螺旋槳槳葉計算機實體造型方法研究[J]. 海軍工程大學學報, 2005, 17(4): 104–107. CHEN Yan-feng, WU Xin-yue. On way of 3D solid modeling of propeller leaf[J]. Journal of Naval University of Engineering, 2005, 17(4): 104–107.

[7]蔡榮泉, 陳鳳明, 馮學梅. 使用Fluent軟件的螺旋槳敞水性能計算分析[J]. 船舶力學, 2006, 10(5): 41–48. CAI Rong-quan, CHEN Feng-ming, FENG Xue-mei. Calculation and analysis of the open water performance of propeller by CFD software Fluent[J]. Journal of Ship Mechanics, 2006, 10(5): 41–48.

[8]顧鋮璋, 鄭百林. 船用螺旋槳敞水性能與槳葉應力的數值分析[J]. 力學季刊, 2011, 32(3): 440–443. GU Cheng-zhang, ZHENG Bai-lin. Numerical analysis of propeller open—water performance and stress distribution in the blade of ship propeller[J]. Journal of Chinese Quarterly of Mechanics, 2011, 32(3): 440–443.

Coupled vibration characteristic analysis of shaft-propeller system of ship based on fluid structure interaction

LI Xiao-jun1, ZHU Han-hua1, XIONG Wei2, WU Ji-dong2
(1. School of Energy and Power Engineering, Wuhan University of Technology, Wuhan 430063, China; 2. Wuchang Shipbuilding Industry Group Company Limited, Wuhan 430060, China)

In order to explore the effect of propeller hydrodynamics and centrifugal force on coupled vibration characteristic of shaft-propeller system, taking Workbench as a tool, modal analysis of shaft-propeller system of ship was completed with the fluid-structure interaction method. The open water performance of propeller was calculated in CFX, the natural frequency and mode shape of shaft-propeller system were analyzed in Ansys taking propeller water pressure and centrifugal force as pre-stress, and the difference of natural frequency of the single model of the shaft, the single model of the propeller and the combination model of shaft-propeller was analyzed. The results show that, the natural frequency of the shaftpropeller system is much lower than the natural frequency of the single model of the shaft and the single model of the propeller; the centrifugal force generated by the rotation of the shaft-propeller system has little effect on the natural frequency; the pre-stress produced by the propeller in the flow field increases the natural frequency of the longitudinal vibration slightly, but the influence is small, which can be ignored in the practical application.

fluid-structure interaction (FSI)；propeller；shaft；Ansys；modal analysis

U664.21

A

1672 – 7649(2017)07 – 0019 – 05

10.3404/j.issn.1672 – 7649.2017.07.004

2016 – 11 – 08；

2016 – 11 – 22

國家自然科學基金重點資助項目(51139005)

李小軍(1994 – )，男，碩士研究生，研究方向為船舶推進技術。