自由拋落式救生艇落水沖擊響應及可靠性分析

蘇石川，戴成龍，史 俊，張 文，張 力，高 強，張順東

(1.江蘇科技大學能源與動力工程學院，江蘇鎮江212003)

(2.鎮江四洋柴油機制造有限公司，江蘇鎮江212003)

自由拋落式救生艇落水沖擊響應及可靠性分析

蘇石川1，戴成龍1，史 俊2，張 文1，張 力1，高 強1，張順東1

(1.江蘇科技大學能源與動力工程學院，江蘇鎮江212003)

(2.鎮江四洋柴油機制造有限公司，江蘇鎮江212003)

自由拋落式救生艇入水過程開始時，艇體與水面接觸瞬間將產生巨大的沖擊力脈沖，這一脈沖力有可能造成救生艇內部結構失效.在救生艇入水過程中，存在固、液、氣三相的耦合，其過程較為復雜.文中以某救生艇為研究對象，主機配置為380 J，應用有限元分析方法對其入水沖擊過程進行數值仿真，分析自由拋落后沖擊載荷對艇體的沖擊響應.最終得出:最大變形發生在入水時間t=0.35 s，即救生艇開始向水面浮升的時，最大變形量為34.3 mm;最大應力發生在t=0.035 s，為95 MPa，應力變化在材料承受范圍內，艇身不會被破壞.

艇體變形;沖擊載荷;自由拋落式救生艇

自由拋落式救生艇以入水速度快、快速動作、避免近體危險等特點，至今已廣泛應用在船舶、海洋平臺等海上設備.針對救生艇結構安全性和可靠性方面的研究，目前國內外已取得了可觀的進步和收獲，主要是通過實驗測量或有限元分析軟件對艇身結構和強度進行分析［1-2］.

在救生艇入水過程開始時，艇體與水面接觸瞬間將產生巨大的沖擊力脈沖，很有可能造成救生艇內部結構失效.因此，文中通過對救生艇入水沖擊過程進行數值仿真模擬，分析其沖擊過程中艇體變形及應力變化，為救生艇的可靠性設計提供參考.

1 滿載救生艇基本參數

文中以某380J船用自由拋落式救生艇為研究對象，該自拋艇的主要技術參數如表1.

表1 救生艇基本參數Table 1 Basic parameters of the lifeboat

取救生艇拋落角度為30°，入水角度為53°，下落高度為20 m，為便于計算，假定自拋艇是規則物體，底面積近似為S=6.7 m×2.55 m=17.1 m2，滿載總質量m=6 475 kg.

根據文獻［3］中的計算方法，經計算得救生艇入水時豎直方向分速度V2H=18.4 m/s.

當救生艇完全浸入水面，即入水深度為3.2 m時，垂向速度VH=11.6 m/s.

當速度減為零時，所用入水時間t=0.15 s，此時間段也是作用在救生艇上的最大受力時間.該過程入水深度h=0.91 m，則救生艇總入水深度H= 4.11 m［3-4］，該過程入水位移為0.91 m，則救生艇總入水深度H=4.11 m，ρ為水密度，取ρ=103kg/m3.

2 模型建立與處理

2.1 模型建立與網格劃分

根據已有相關圖紙和數據，進行艇體三維建模.再用HyperMesh有限元劃分軟件對模型進行網格劃分.艇體網格劃分模型如圖1.隨后將畫好網格的救生艇模型導入有限元計算分析軟件Patran中.

在救生艇入水沖擊仿真過程中需對流體(空氣、水)部分建立歐拉區域，采用8節點六面體網格單元進行劃分.由于艇體十分復雜，在不影響反映艇體變形整體情況的條件下，為節省計算量，采用節點數量較少，救生艇及流體的節點數和單元數見表2.

補水泵站主廠房抗浮穩定及基底應力均滿足規范要求，井筒及其內部板梁的配筋滿足強度要求。連通洞結構設計符合規范和已有工程經驗，結構的穩定滿足規范要求，配筋滿足強度要求。

圖1 網格劃分Fig.1 Meshing map

表2 節點及單元數Table 2 Number of nodes and elements

最后調整救生艇的入水角度，完成最終的有限元網格分析模型，將其導入有限元分析前處理軟件中.

文中所研究的救生艇用聚酯纖維玻璃鋼的一般性能常數如表3.

表3 聚酯纖維玻璃鋼一般性能Table 3 General performance of polyester fiber glass steel

救生艇玻璃鋼材料采用Von Mises屈服模型，屈服應力取290 MPa［5］.

流體域主要包括空氣域和水域，流體域歐拉單元初始條件的定義內容主要為比內能和密度.流體域中，空氣和水的狀態方程參數值如表4.表中:ρ為密度;γ為空氣比熱比;E為單位質量的內能;a1為水的體積彈性模量.

表4 流體域狀態方程參數Table 4 State parameters of fluid field equation

救生艇入水沖擊仿真計算的初始條件主要包括救生艇結構與流體間的流-固耦合邊界、救生艇初速度、歐拉單元初始狀態等條件.在仿真過程主要模擬救生艇拋落過程的第4階段，即救生艇在接觸水面到法向加速度達到最大時刻的沖擊過程.由上文計算可知，救生艇拋落的水平傾斜角為35°，接觸水面時刻的傾斜角度為53°，拋落高度為20 m.

隨后對艇體模型進行速度加載，并定義流-固耦合邊界.

為提高計算效率，文中只分析救生艇開始接觸水面到加速度達到最大時的沖擊過程.初始，救生艇最低端距水面距離為0.05 m，分析計算時間定義為0.36 s，初始時間步長為1×e-7s.定義初始條件的完整有限元計算分析模型如圖2所示［6-8］.艇身并且應力分布逐漸趨于均勻;由圖d)可以看出:沖擊過程結束后，由于能量不能及時傳遞到艇體外界，應力變化存在一定的殘余響應.當艇體將要浮出水面時，即當t=0.35 s時，沖擊殘余應力已均勻分布于整個艇體，只在救生艇的拐角結構等處出現局部應力集中的情況.

圖2 救生艇入水初始邊界條件Fig.2 Initial diving boundary conditions of lifeboat

3 艇體形變及應力分析

3.1 等效應力分析

救生艇以53°入水角度跌入水中，救生艇與水面接觸的瞬間，艇首的等效應力迅速增大，并將沖擊能量以應變、應力形式傳遞到艇體及動力裝置.沖擊過程不同時刻t的艇體應力變化如圖3.

由圖a)可以看出:當沖擊開始時，艇首與水面正面沖擊，沖擊部位的應力值迅速增大.隨著救生艇著水區域的增大，艇首的最大應力值也隨之迅速增大，并且沖擊能量以應力、應變的形式，傳遞到整個艇身;由圖b)可以看出:當沖擊時間t=0.058 s時，艇首出現比較明顯的應力集中情況并達到最大值;由圖c)可以看出:當t=0.1 s時，救生艇大部分艇身已浸入水中，沖擊載荷已由艇首傳遞到整個

圖3 艇體不同時刻應力分布云圖Fig.3 Hull stress distribution nephogram of different time

圖4 救生艇最大應力變化曲線Fig.4 Maximum stress change curve of lifeboat

沖擊過程中，最大應力的大小及位置隨著沖擊時刻的變化而不斷變化.其沖擊過程中，救生艇最大應力隨時間的變化情況如圖4.從圖3，4可以看到，最大應力σ分布范圍主要集中在艇首與水面接觸的初始區域內，在t=0.035 s時應力分布出現的最大值為95 MPa，發生在救生艇首部與水面初始接觸位置.由救生艇玻璃鋼材料特性可知，此時的應力小于玻璃鋼的最大安全應力120 MPa，所以認為救生艇受沖擊后不會發生破壞狀況.

3.2 艇體變形分析

圖5為救生艇在不同沖擊時刻艇體的變形云圖，從救生艇的變形圖中分析可知，最大變形發生在t=0.35 s時，即救生艇開始向水面浮升的時刻，最大位移為33.3 mm，發生在救生艇尾部.不同沖擊時刻救生艇出現不同的變形響應.由圖a)可以看出:沖擊過程開始瞬間，艇首最先與水面接觸，沖擊能量瞬間集中于救生艇首部，造成了如圖a)所示的變形響應.隨后沖擊響應沿艇首迅速向艇體傳遞，所以艇首的沖擊變形首先迅速增大，然后稍有減小并趨于穩定.救生艇入水初始過程中，首部相當于受到水的約束作用，致使救生艇尾部受到巨大慣性力及扭轉力矩作用，扭轉力矩使救生艇尾部產生彎曲變形，救生艇的彎曲變形要遠大于沖擊載荷所造成的沖擊變形.當救生艇接近入水深度最低點時，救生艇已完全浸入水中，此時扭轉力矩作用消失，尾部變形達到最大值，如圖b).

圖5 艇體應變響應云圖Fig.5 Nephogram of hull strain response

分別取救生艇艇首、艇側、艇尾及頂部觀的測點，對救生艇入水沖擊過程中的應力變化情況進行分析，其觀測點節點編號分別為node 81、node 5 310、node 3 876和node 6 181.圖6為各個觀測點沖擊過程中的位移時程曲線.從圖中可以看出:艇首觀測點(node 81)最先與水面接觸，所以位移D迅速增大并首先達到最大值，之后隨著時間t的變化先減小而后略有增加并趨于穩定值.這種現象發生的主要原因是因為玻璃鋼具有較高的韌性.距離沖擊位置較遠的艇側觀測點(node 5 310)與頂部觀測點(node 6 181)由于沒有與水面直接撞擊，所以沖擊變形響應相對較為緩慢.救生艇艇尾觀測點(node 3 876)的沖擊位移響應隨時間變化較為平緩.不同觀測點的位移響應各不相同，這主要是因為救生艇在入水過程中，與救生艇結構、入水姿態以及沖擊能量在救生艇內部的傳遞和轉化過程有關.

圖6 不同觀測點位移變化曲線Fig.6 Observation points maximum displacement change curve

4 結論

1)艇體在沖擊過程開始時，最大應力迅速增大，并在t=0.035 s時應力分布出現最大值為95 MPa，艇首觀測點變形量于t=0.35 s時達到最大值為33.3 mm，之后隨著時間的變化先減小而后略有增加并趨于穩定值.

2)艇側觀測點與頂部觀測點沖擊變形響應相對較為緩慢.救生艇艇尾觀測點的沖擊位置響應隨時間變化較為平緩.

References)

［1］Re A S，Veitch B.A comparison of three types of evacuation system［J］.Transactions-society of Naval Architects and Marine Engineers，2007(115):119-139.

［2］Bradley Elliott.Computer simulation of a freefall evacuation system in a range of initial conditions［R］.Canada:NRC Research Press，2003.

［3］ 黃春平，陸軍，秦嘉岷.自由降落救生艇試驗過程探討［J］.造船技術，2012(3):37-40.Huang Chunping，Lu Jun，Qin Jiamin.Research on free-fall lifeboats test process［J］.Marine Technology，2012(3):37-40.(in Chinese)

［4］ 賀玉英，楊金芳.船舶自拋救生艇［J］.船舶，1996 (5):30-33.

［5］黃燕玲，郭強波，廖文，等.小型玻璃鋼開敞艇結構強度分析［J］.船海工程，2013，42(1):21-22.Huang Yanling，Guo Qiangbo，Liao Wen，et al.Direct calculation of structural strength of a small FRP yacht［J］.Ship＆Ocean Engineering，2013，42(1): 21-22.(in Chinese)

［6］Von K T.The impact of seaplane floats during landing［R］.Washington USA:National Advisory Committee for Aeronautics，NACA TN 321，1929:1-8.

［7］Wagner V H.Phenomena associated with impacts and sliding on liquid surfaces［J］.Z Angew Math Mech，1932(12):193-215.

［8］Shiffman M，Spencer D C.The force of impact on a cone striking a water surface［J］.Commun Pure Appl Mech，1951(4):379-417.

［9］蘇石川，王憲淼，孫江，等.自由拋落式救生艇推進系統落水沖擊響應計算分析［J］.中國造船，2013，54(3):114-122. Su Shichuan，Wang Xianmiao，Sun Jiang，et al.Simulation and analysis of response of propulsion system of free-fall lifeboat to water entry impact［J］.Shipbuilding of China，2013，54(3):114-122.(in Chinese)

［10］ 蘇石川，張未軍，王永，等.外載荷作用下船舶推進系統非線性動力學分析［J］.中國造船，2012，53 (1):79-87.Su Shichuan，Zhang Weijun，Wang Yong，et al.Calculation and analysis of non-linear dynamics of propulsion system with external load［J］.Shipbuilding of China，2012，53(1):79-87.(in Chinese)

［11］Det Norske Veritas.Design of free fall lifeboats［S］.Oslo Norway:Det Norske Veritas，2010.

［12］Karim M M，Iqbal K S，Khondoker M R H，et al.Influence of falling height on the behavior of skid launching free-fall lifeboat in regular waves［J］.Journal of Applied Fluid Mechanics，2011，4(1):77-88.

［13］ 宋長福.民機機身機構入水沖擊問題數值仿真研究［D］.上海:上海交通大學，2011.

［14］ 王軍.魚雷入水沖擊動力學仿真研究［D］.云南昆明:昆明理工大學，2010.

(責任編輯:貢洪殿)

Analysis of water entry impact response of free fall lifeboats and its reliability

Su Shichuan1，Dai Chenglong1，Shi Jun2，Zhang Wen1，Zhang Li1，Gao Qiang1，Zhang Shundong1

(1.School of Energy and Power Engineering，Jiangsu University of Science and Technology，Zhenjiang Jiansu 212003，China) (2.Zhenjiang Siyang Diesel Engine Manufacturing Co.Ltd.，Zhenjiang Jiangsu 212003，China)

When the lifeboat falls into water，the impact pulse created is so great that the structure of the boat may be destroyed.In the process，there exists the coupling of gas，liquid and solid phases，so the process may be very complex to be analyzed.A lifeboat whose main engine is 380 J is used in this study as the research object.The finite element analysis method is used to simulate the lifeboat water impact process.Finally，the conclusion is made that when t=0.35 s，the biggest displacement is 34.3 mm.When t=0.035 s，the biggest stress is 95 MPa.The biggest stress is within the material stress range and the body of the boat will not be destroyed.

ship deformation;impact load;free fall lifeboat

U662.3

A

1673-4807(2015)06-0511-05

10.3969/j.issn.1673-4807.2015.06.001

2015-06-08

資金項目:鎮江市科技支撐基金資助項目(GY2014047)

蘇石川(1963—)，男，教授，研究方向為動力裝置性能優化及可靠性.E-mail:15606103818@163.com

蘇石川，戴成龍，史俊，等.自由拋落式救生艇落水沖擊響應及可靠性分析［J］.江蘇科技大學學報(自然科學版)，2015，29(6):511-515.