液貨晃蕩對雙殼油船碰撞性能的影響研究

吳文鋒，楊雨濱，盧金樹，鄧佳佳，朱發(fā)新

(浙江海洋大學(xué) 海運與港航建筑工程學(xué)院，浙江 舟山 316022)

液貨晃蕩對雙殼油船碰撞性能的影響研究

吳文鋒，楊雨濱，盧金樹，鄧佳佳，朱發(fā)新

(浙江海洋大學(xué) 海運與港航建筑工程學(xué)院，浙江 舟山 316022)

利用ANSYS/LS-DYNA分析載貨情形下的雙殼油船發(fā)生碰撞導(dǎo)致艙內(nèi)液貨發(fā)生晃蕩時，晃蕩載荷對舷側(cè)結(jié)構(gòu)碰撞性能的影響，并將所得結(jié)果與空載狀態(tài)下的碰撞事故進(jìn)行比對，發(fā)現(xiàn)液體晃蕩對碰撞性能的影響主要體現(xiàn)在當(dāng)撞擊船接觸到被撞擊船內(nèi)殼時，艙內(nèi)液體動能較大，對碰撞性能產(chǎn)生顯著影響，艙內(nèi)液貨能有效吸收撞擊船的撞擊能量，使得撞擊速度快速下降，同時由于艙內(nèi)液貨與艙壁之間的耦合作用，使得碰撞力增大，導(dǎo)致內(nèi)殼受損更為嚴(yán)重，造成內(nèi)殼提前破裂。

船舶碰撞；流固耦合；液貨晃蕩；數(shù)值模擬

油船碰撞持續(xù)時間較短，且涵蓋多種非線性問題，無疑增大了探索相關(guān)科學(xué)問題的難度。油船碰撞其本質(zhì)上是一個流固耦合的問題[1]，對于載有液體的雙殼油船而言，這種流固耦合同時存在于船體外部及船體內(nèi)部區(qū)域。

對于載貨油船而言，其具備碰撞持續(xù)時間短、產(chǎn)生的撞擊能量大、艙內(nèi)液貨會在短時間內(nèi)發(fā)生液貨晃蕩等特點。據(jù)Kristjan T的船模試驗表明[2]，液艙載貨發(fā)生船舶碰撞時，液體晃蕩產(chǎn)生的晃蕩沖擊能量會對船舶碰撞性能造成顯著影響，在撞擊能量作用下，被撞擊船船身的晃動會導(dǎo)致艙內(nèi)液體的晃蕩，甚至?xí)a(chǎn)生巨大的瞬間晃蕩沖擊能量，該能量對被撞擊船造成的損傷主要在于對液貨艙艙壁的沖擊作用及晃蕩壓力對艙壁造成的局部損傷。因此，研究液貨晃蕩對雙殼油船碰撞性能的影響對于增強船舶安全和保護(hù)環(huán)境具有十分重要的意義。

以往對油船碰撞的仿真研究中，幾乎未考慮到液艙中液體與艙壁結(jié)構(gòu)的相互耦合作用，即忽視了液貨的影響[3-8]。隨著對船舶碰撞問題的深入研究，有學(xué)者已開始研究裝載液貨的雙殼油船碰撞損傷機理。娜日薩等[9-10]運用數(shù)值分析方法分析了VLCC液艙滿載情形下液貨發(fā)生晃蕩的全過程，并對晃蕩作用下的船體結(jié)構(gòu)強度進(jìn)行評估。張愛年等[11]運用ANSYS/LS-DYNA數(shù)值軟件，采用3種不同的數(shù)值仿真方法對碰撞過程中滿載狀態(tài)下油艙內(nèi)流固耦合作用進(jìn)行研究，根據(jù)研究結(jié)果可知，在實際工程運用中，可采用拉格朗日—歐拉算法，不僅能夠滿足實際工程的精度需求還能大大縮減數(shù)值計算所需時間。張延昌等運用MSC.Dytran有限元軟件分析了VLCC在滿載、壓載工況下的碰撞損傷特性，但缺乏艙內(nèi)液體晃蕩的直觀展示及艙內(nèi)液貨未滿載時的晃蕩特性。Kristjan Tabri等開展一系列船舶碰撞的模型試驗，在試驗過程中未能考慮到碰撞船舶裝載液貨的情形。張日紅運用ANSYS/LS-DYNA針對滿載狀態(tài)下的油船發(fā)生碰撞事故時，分析液體與舷側(cè)結(jié)構(gòu)之間的耦合作用，對碰撞性能開展相關(guān)研究。崔杰等分析80%裝載率狀態(tài)下艙內(nèi)液體對雙層舷側(cè)碰撞性能的影響。

本文結(jié)合以往的工作經(jīng)驗，以雙殼油船為研究對象，應(yīng)用ANSYS/LS-DYNA非線性有限元軟件，利用仿真技術(shù)建立船舶碰撞數(shù)值模型，為凸顯液貨的晃蕩效果，分別對空載和載貨率為30%的2種情景開展數(shù)值計算，通過對碰撞力、結(jié)構(gòu)損傷變形等比較分析，得出液貨晃蕩對雙殼油船碰撞性能的一般影響規(guī)律。

1 流固耦合算法

本文運用任意拉格朗日—歐拉算法(ALE, Arbitrary Lagrange Euler)處理雙殼油船艙內(nèi)液貨與艙壁之間的流固耦合問題。該算法匯集了仿真軟件中拉格朗日算法和歐拉算法的優(yōu)點，能有效處理有限空間內(nèi)的大位移和自身大變形問題。

ALE算法在復(fù)雜的載荷條件下，能夠很好的處理流體與結(jié)構(gòu)之間的相互作用。ALE網(wǎng)格在計算時，首先得到固體網(wǎng)格的變形特征，并將相關(guān)計算結(jié)果導(dǎo)入空間節(jié)點中，隨后根據(jù)固體網(wǎng)格的變形特征重新劃分固體網(wǎng)格，并通過空間網(wǎng)格節(jié)點將變形后的邊界條件分配到新劃分的固體網(wǎng)格中。結(jié)構(gòu)參數(shù)和流體單元之間的傳遞可用關(guān)鍵字*CONSTRAIN-ED_LAGRANGE_IN_SOLID處理。

2 碰撞方案

根據(jù)對歷年油船碰撞事故的統(tǒng)計，相撞船舶原型尺寸選取5萬t雙殼油船和5萬t散貨船。相撞船舶的主要尺寸如表1所示。

本文設(shè)計的碰撞方案為，散貨船以5 m/s的速度，垂直對中撞擊靜止的油船，撞擊位置選取在舷側(cè)最為薄弱的縱骨之間，碰撞方案如圖1所示。

表1 相撞船舶的主要尺寸 m

圖1 碰撞方案

3 模型建立

為縮減計算時間，在建立數(shù)值模型時，對相撞船舶做了相應(yīng)簡化處理。撞擊船船首材料選取彈塑性材料，鋼板厚度與實際船首厚度一致，內(nèi)部構(gòu)件根據(jù)實際船首結(jié)構(gòu)進(jìn)行簡化。將撞擊船除船首外的所有結(jié)構(gòu)簡化成為一個貨艙，可通過改變該貨艙的尺寸及質(zhì)量來確保撞擊船的質(zhì)量和重心位置與實際一致?？紤]到碰撞的局部性，同時結(jié)合船舶重心及質(zhì)量等因素，將雙殼油船簡化為一個液貨艙，該液貨艙的舷側(cè)鋼板厚度、內(nèi)外殼間距等均與實際船舶一致，艙段結(jié)構(gòu)采用塑性動態(tài)模型模擬，如圖2所示。圖3為艙內(nèi)液貨模型。

計算過程中，設(shè)定失效應(yīng)變?yōu)?.1，材料模型采用軟件提供的Cowper-Symonds本構(gòu)方程，并考慮材料應(yīng)變硬化影響。接觸算法采用主從面接觸和自動單面接觸。對可能發(fā)生碰撞的區(qū)域進(jìn)行網(wǎng)格細(xì)化，其余部分采用較粗網(wǎng)格。

在LS-DYNA中，油船艙內(nèi)液貨壓力與體積的關(guān)系可通過關(guān)鍵字*EOS-GURNEISEN來描述，該狀態(tài)方程定義壓縮材料的壓力見式(1)。

+(γ0+aμ)E,

(1)

式中：p為壓力；ρ0為流體初始密度；C為撞擊速度—質(zhì)點速度曲線截距；μ為比體積，其值為ρ/(ρ0-1)；ρ為流體過程中的密度，ρ的值是隨計算過程變化而變化的；γ0為格林愛森常數(shù)；a為γ0的一階體積修正系數(shù)；S1、S2、S3為撞擊速度—質(zhì)點速度曲線的斜率系數(shù)；E為單位體積內(nèi)能。各參數(shù)值見表2。

表2 液貨材料參數(shù)

圖2 相撞船舶數(shù)值模型

圖3 艙內(nèi)液貨模型

4 結(jié)果分析

4.1 艙內(nèi)液貨響應(yīng)

圖4為液貨平均動能隨撞擊船位移變化關(guān)系曲線。圖5為撞擊過程中液貨狀態(tài)的變化。

圖4 液貨平均動能-位移曲線

圖5 液貨狀態(tài)變化圖

觀察圖4曲線，可將液貨動能變化大致分為3個階段。第一階段，被撞船位移在約0.8 m之前，液貨動能幾乎為零。在該階段撞擊船船首已撞擊被撞擊船的外殼，但由于液貨存在黏性，因此存在一個滯后效應(yīng)。此外，被撞擊船的漂移速度較慢，因此液貨相對被撞擊船而言處于靜止?fàn)顟B(tài)，如圖5(a)所示。第二階段，撞擊船位移約在0.8～ 3.1 m之間，在該階段液貨明顯有動能，但該階段液貨流動較為緩慢，動能較小，如圖5(b)所示。在該階段被撞擊船外殼已發(fā)生破裂，船首開始通過雙層殼之間的構(gòu)件向內(nèi)殼進(jìn)發(fā)，此時在力的傳遞作用下，被撞擊船漂移速度開始逐漸增大，艙內(nèi)液貨在黏性力、慣性力以及碰撞力作用下開始緩慢的流動。第三階段在3.1～6.4 m之間，曲線開始急劇上升，液貨動能迅速增大，如圖5(c)所示。在該階段撞擊船船首已觸及內(nèi)殼，內(nèi)殼發(fā)生變形。這表明當(dāng)撞擊船觸及到內(nèi)殼時，艙內(nèi)液貨會發(fā)生劇烈運動。此時導(dǎo)致液貨發(fā)生劇烈運動的原因是流固耦合力。

根據(jù)圖5可以看出，在撞擊船位移為1 m時，液貨狀態(tài)基本沒有變化，說明撞擊船沒有接觸到被撞船的內(nèi)殼板。當(dāng)撞擊船位移為2 m時，液貨狀態(tài)發(fā)生微小變化，說明撞擊船剛觸及到內(nèi)殼板，內(nèi)殼板還沒有發(fā)生變形，液貨運動不劇烈。而撞擊船位移為5 m時，液貨狀態(tài)發(fā)生顯著變化，這時內(nèi)殼板出現(xiàn)明顯變形，液貨劇烈運動。

4.2 碰撞力分析

圖6為船舶空載和載貨情況下碰撞力隨位移變化關(guān)系曲線。該過程可分為3個階段：第一階段在2條曲線第一個交點之前，顯然載貨狀態(tài)的碰撞力要高于空載狀態(tài)的碰撞力，這是由于考慮載貨時，被撞擊船的總質(zhì)量比空載要大。第二階段為2條曲線2個交點之間，此時載貨狀態(tài)的碰撞力要低于空載狀態(tài)的碰撞力。由于經(jīng)歷完第一階段后，載貨狀態(tài)的被撞擊船外殼已然發(fā)生破裂，而空載狀態(tài)的則要至位移約為3 m處才開始破裂，這是由于2類情形下被撞擊船質(zhì)量不一致，導(dǎo)致被撞擊船漂移速度不一致。第三階段，載貨狀態(tài)的碰撞力逐漸高于空載狀態(tài)的碰撞力，此時的碰撞力除固體之間接觸產(chǎn)生之外，還有很大部分是由于流固耦合力的作用。此時撞擊船船首已接觸到被撞擊船內(nèi)殼，艙內(nèi)液貨隨著內(nèi)殼的變形開始劇烈運動，液貨與內(nèi)殼結(jié)構(gòu)之間呈現(xiàn)出非常強烈的的耦合作用。艙內(nèi)液貨受到內(nèi)殼的擠壓變形作用，使得撞擊船運動明顯減緩，這種減緩作用表現(xiàn)為碰撞力的提高。

圖6 碰撞力-位移曲線

4.3 撞擊船速度變化分析

圖7為撞擊船速度隨位移變化關(guān)系曲線。由圖7可知，在船舶碰撞初期，二者曲線較為吻合，在撞擊船前進(jìn)位移到達(dá)3.1 m時，2條曲線開始出現(xiàn)明顯差異，具體體現(xiàn)撞擊船在撞擊載貨狀態(tài)的油船時速度下降更快。這說明，撞擊船在未觸及內(nèi)殼時，艙內(nèi)液貨幾乎不會對撞擊船的運動造成影響；當(dāng)撞擊船接觸到內(nèi)殼后，艙內(nèi)液貨明顯對撞擊船的運動產(chǎn)生阻礙作用，使得相同位移下撞擊船速度減小的更加迅速。

4.4 吸能分析

圖8為2種碰撞情形下被撞擊船吸能隨位移變化曲線。由圖8知，被撞擊船內(nèi)艙與撞擊船位移成正比。該曲線劃分為3個階段。第一階段在位移為0.8 m之前，此時，被撞擊船能量變化主要是依靠外殼及其舷側(cè)肋板結(jié)構(gòu)變形，由于該階段外殼尚未破裂，且損傷模式基本一致，因此2條吸能曲線基本吻合，艙內(nèi)液貨對舷側(cè)外板吸能幾乎不產(chǎn)生影響。第二階段在位移為0.8～3.1 m之間，該階段吸能增加主要是依靠外殼持續(xù)變形而發(fā)生破裂，以及雙層殼之間結(jié)構(gòu)損傷變形。通過圖5(b)可知，此時液貨已開始緩慢流動，為此載貨狀態(tài)船舶吸能稍高于空載狀態(tài)船舶的吸能。第三階段在位移為3.1 m之后。此時撞擊船已觸及內(nèi)殼，致使內(nèi)殼開始發(fā)生變形。該階段，空載狀態(tài)船舶能量變化主要依靠內(nèi)殼變形；而載貨狀態(tài)船舶能量變化除內(nèi)殼變形外，還有內(nèi)殼與液貨的耦合作用，使得液貨動能大幅度增加，從而導(dǎo)致載貨狀態(tài)船舶能量吸收迅速上升且大幅度超過空載狀態(tài)船舶能量吸收。

圖8 被撞擊船內(nèi)能-位移曲線

4.5 損傷變形分析

在碰撞過程中，被撞船舷側(cè)結(jié)構(gòu)可通過自身變形來吸收撞擊能量，能量吸收性能取決于舷側(cè)結(jié)構(gòu)的損傷變形程度，而舷側(cè)結(jié)構(gòu)吸能性能的好壞直接關(guān)乎其抗撞擊能力。

圖9反映了2種碰撞情形下被撞擊船外殼損傷變形，由圖可見，外殼損傷變形幾乎相同。外殼損傷均以膜拉伸為主，且變形區(qū)域及程度幾乎一致。

圖10反映了2種碰撞情形下被撞擊船內(nèi)殼損傷變形情況，由圖10知，碰撞結(jié)束后，內(nèi)殼發(fā)生了以膜拉伸為主的變形，且尚未破裂，但顯然載貨狀態(tài)下的內(nèi)殼變形程度和范圍均比空載狀態(tài)下的內(nèi)殼變形程度和范圍要大，且載貨狀態(tài)下內(nèi)殼破裂時間將提前。造成該現(xiàn)象的原因在于，內(nèi)殼除受到撞擊船撞擊的作用外，還受到液貨晃蕩產(chǎn)生的動壓力，內(nèi)殼變形是液貨壓力與接觸力共同作用的結(jié)果。

(a) 空載 (b) 載貨圖9 外殼損傷變形圖

(a) 空載 (b) 載貨圖10 內(nèi)殼損傷變形圖

5 結(jié)束語

本文以雙殼油船為研究對象，利用ANSYS/LS-DYNA對比分析載貨狀態(tài)和空載狀態(tài)下雙殼油船舷側(cè)結(jié)構(gòu)碰撞性能，得出如下結(jié)論。

1)艙內(nèi)液貨在碰撞前期對碰撞性能影響不大，對舷側(cè)結(jié)構(gòu)損傷變形幾乎不會產(chǎn)生影響。

2)艙內(nèi)液貨晃蕩主要在碰撞后期對雙殼油船碰撞性能產(chǎn)生較為明顯的影響。在內(nèi)殼發(fā)生變形時，艙內(nèi)液貨會發(fā)生劇烈運動，將會產(chǎn)生較大的動壓力，對結(jié)構(gòu)碰撞性能產(chǎn)生顯著影響，艙內(nèi)液貨的存在極大降低了撞擊船運動速度，同時由于流固耦合力，使得碰撞力大幅增加，使得內(nèi)殼變形嚴(yán)重，最終導(dǎo)致內(nèi)殼破裂時間提前。故研究雙殼油船碰撞時，不能忽略艙內(nèi)液貨晃蕩的影響。

[1] Z liu，J Amdahl. A new formulation of the impact mechanics of ship collisions and its application to a ship-iceberg collision [J]. Marine Structures, 2010, 23(3):360-384.

[2]K Jabri,J matu siak.P Varsta. Sloshing interaction in ship collisions-An experimental and Numerical study [J].Ocean Engineering,2009(36)：1366-1376.

[3] Saydam D, Frangopol D M. Performance assessment of damaged ship hulls [J]. Ocean Engineering, 2009,36(17-18):1366-1376.

[4] PT Pedersen. Review and application of ship collision and grounding analysis procedures [J]. Marine Structures, 2010, 23(3):241-262.

[5] 肖曙明,劉敬喜,王娜,等.雙殼船內(nèi)外殼結(jié)構(gòu)耐撞性試驗與仿真研究[J].應(yīng)用力學(xué)學(xué)報,2013,30(2)：257-261.

[6] S Haris, J Amdahl. Analysis of ship-ship collision damage accounting for bow and side deformation interaction [J]. Marine Structures, 2013, 32(7): 18-48.

[7] P Hogstr?m,J W Ringsberg. An extensive study of a ship's survivability after collision-A parameter study of material characteristics, non-linear FEA and damage stability analyses [J]. Marine Structures, 2012, 27(1):1-28.

[8] P T Pedersen ,Y Li. On the global ship hull bending energy in ship collisions [J]. Marine Structures, 2009, 22(1):2-11.

[9] Na Risa, HL Ren，ZL Gao, et al. Simulation of sloshing and structure response in VLCC tanks [J]. Journal of Ship Mechanics, 2007, 11(6):879-887.

[10] 娜日薩.VLCC液艙晃蕩仿真及結(jié)構(gòu)強度評估方法研究[D].哈爾濱,哈爾濱工程大學(xué),2006.

[11] A zhang,K Suzuki. A comparative study of numerical simulations for fluid-structure interaction of liquid-filled tank during ship collision [J].Ocean Engineering，2007,34(5-6)：645-652.

The ANSYS/LS-DYNA was used to study the effect of sloshing load on ship-ship collision which resulted in liquid cargo sloshing.Compared with the light condition,the influence of liquid sloshing on the impact performance is found mainly reflected in the impact when the sticking ship had been attached to the struck ship's inner shell with large kinetic energy of liquid.A significant impact on the performance of the collision shows that the liquid cargo could effectively absorb the impact energy of the striking ship,which made the collide velocity decreased rapidly.At the same time,due to the coupling of liquid cargo compartment and the bulkhead between the collision force,the increased collision damaged the inner shell seriously,resulting the shell ruptured in advance.

ship-ship collision;fluid-structure coupling;liquid sloshig;numerical simulation

浙江省自然科學(xué)基金青年基金(LQ16E090003，LQ14E090001)；舟山市科技局項目(2015C41009)

吳文鋒(1984-)，男，江西上饒人，講師，博士，主要從事船舶安全與防污染研究。

U661.42

10.13352/j.issn.1001-8328.2016.06.003

2016-07-25