冰体质量和撞击角度对船首结构碰撞性能的影响

2019-11-12 07:54王甫超刘海冬杜成忠
造船技术 2019年5期
关键词:船首外板船体

张 健 ,王甫超,2*,刘海冬,杜成忠

(1.江苏科技大学 土木工程与建筑学院,江苏 镇江 212003;2.浙江大华技术股份有限公司,浙江 杭州 310000;3.上海外高桥造船有限公司,上海 200137;4.江苏新扬子造船有限公司,江苏 靖江214500)

0 引 言

我国国民经济正持续快速发展,海洋经济已成为我国国民经济的新增长点。开通北极航道将有助于我国减少对传统航道的依赖,减少航运成本,缩短航行时间,降低航运安全风险,进一步推动我国北极战略部署,但海冰与海上结构物相互作用问题、海冰对结构操作和安全问题也日益突出。我国对于极地原油运输船、极地甲板运输船等方向的研究还处于起步阶段,开展对船-冰碰撞的研究不仅有利于保护人员财产安全,防止海洋污染,还能为开辟极地航道提供理论依据。

针对船-冰碰撞对船首结构的影响,前人已做了大量研究。张充霖[1]应用大型有限元软件MSC.Dytran计算渤海地区3种不同形状的海冰与一艘典型散货船首部碰撞时的动态响应,对3种碰撞工况下船首结构的碰撞力大小、构件吸能情况和损伤变形情况进行分析,得出船-冰碰撞下船首构件的失效规律和能量吸收机制受冰体形状的影响。张健等[2]利用数值仿真的方法,通过变化撞击参数(碰撞速度、冰体形状、碰撞位置等)研究极地船舶在撞击冰山时船舶构件的动态响应差异,讨论船舶与冰体碰撞工况下的损伤变形、碰撞力等变化特性,进而科学地得出冰体形状、冰体质量、船速等撞击参数对船-冰碰撞的影响机理。王健伟等[3]通过改变冰层厚度与船舶速度研究破冰船的动态结构响应,分析船首结构的变形损伤、碰撞冰力的大小和结构吸能等特性。

本文利用ANSYS建立一艘16万吨级的油船首部模型,并使用非线性有限元软件LS-DYNA求解分析,研究在不同冰体质量、碰撞角度等撞击参数下的结构响应问题,得到冰载荷下船首结构的损伤变形、碰撞力等方面的变化特征,对分析船舶与冰体的碰撞性能具有参考价值。

1 船-冰有限元模型

1.1 材料模型

为真实模拟碰撞过程中船体材料的特性,参考文献[4],选用线性强化的Cowper-Symonds弹塑性模型作为材料本构模型,其基本参数如表1所示。国内外学者在对船-冰碰撞进行数值仿真研究时,普遍认为很难建立海冰的本构模型,因为海冰材料受自然因素影响较大,比如气候、温度、盐度等。根据收集整理的资料[5],选择国内外学者普遍使用的冰体本构模型,该材料模型不仅可以考虑应变率对von Mises失效应力、切线模量和弹性模量的影响,还能够考虑屈服应力受应变率的影响。冰体材料参数如表2所示。

表1 船体钢的塑性材料参数

表2 冰体材料参数

1.2 有限元模型

选取一艘极地航行油船首部作为研究对象,因冰载荷下船体损伤具有局部性,故有限元模型只建立油船首部结构。油船首部碰撞区域的结构较复杂,包含甲板平台、船体外板、舷侧纵骨、水密舱壁、横框架、桁材等,在建立模型时尽可能使用板单元进行模拟[6]。船首结构有限元模型如图1所示。虽然冰山形状各异,但是船-冰碰撞的接触区域具有局部性,远离碰撞区域的冰体形状对船-冰碰撞仿真计算结果影响不大。借鉴挪威船级社推荐的冰山数值仿真形状[7],在ANSYS中建立正方体冰体模型如图2所示。

图1 船首网格划分图

图2 冰体模型图

1.3 碰撞方案

在船舶实际碰撞过程中,撞击参数具有多样性[8]。从冰体质量、撞击角度研究船舶碰撞结构响应,使用侵蚀接触模拟舷侧与冰山表面的接触。碰撞方案如表3所示。

表3 碰撞方案

所有工况下的船舶初速度都为8 m/s。不同撞击角度下的船-冰碰撞有限元模型如图3所示。

图3 船-冰碰撞模型图

2 计算结果

2.1 冰体质量对船-冰碰撞的影响

图4为船舶在碰撞不同尺寸冰体过程中的碰撞力随时间变化曲线。由图4可知,碰撞过程持续1 s左右,由于在碰撞过程中船体与冰体发生冲撞挤压作用,碰撞力呈现明显的非线性特征,在碰撞过程的不同阶段船体结构出现了不同程度的卸载现象,伴随有多个峰值出现。值得指出的是,某种构件或者冰体的失效或破坏都以碰撞力的卸载反映出来,因此碰撞力会在一定程度上反映构件失效情况和结构强弱程度。对比分析3种碰撞工况,可以得到船冰首次接触的时间基本在0.01 s左右。边长为6 m的冰体在0.98 s时的碰撞力达到最大峰值110 MN,边长为12 m的冰体在0.85 s时的碰撞力达到最大峰值219 MN,边长为24 m的冰体在0.71 s时的碰撞力达到最大峰值271 MN。船-冰碰撞的碰撞力峰值随冰体质量的增大而增大,碰撞力峰值产生的时刻也随之提前。

图4 不同工况下碰撞力时历曲线

为更好地分析船舶以8 m/s的初速度撞击不同尺寸冰体后船体结构损伤的形态与特征,选取碰撞过程中一系列典型的有效塑性应变与变形图。图5和图6为船首外板及舷侧纵骨在1 s时刻的变形损伤和应力分布情况:球鼻艏的碰撞区域受到挤压产生弯曲变形,舷侧纵骨及船首外板被冰体挤压弯曲变形严重;船舶受边长为6 m的冰体撞击后,船体结构损伤变形区域主要分布在与冰体碰撞区域,局部高应力区域分布与变形区域保持一致,舷侧纵骨和船首外板损伤变形不严重;随着冰体质量的增大,舷侧纵骨和船首外板产生严重的撕裂、弯曲及凹陷,船首结构损伤程度逐渐增大。从图7可知,与船首接触区域的冰体也产生了大量的破碎、失效现象,冰体质量越大,冰体破碎程度越小。

图5 不同冰体尺寸下船首外板的变形损伤和应力分布图

图6 不同冰体尺寸下舷侧纵骨的变形损伤和应力分布图

图7 不同工况下冰体变形损伤和应力分布图

2.2 撞击角度对船-冰碰撞的影响

图8给出了不同碰撞角度下船-冰碰撞的碰撞力曲线。当角度为0°时,碰撞力峰值为219 MN;当角度为30°时,碰撞力峰值为182 MN;当角度为60°时,碰撞力峰值为109 MN。当船冰的撞击角度不断增大时,船首与冰体正面的碰撞接触面积不断减小,当增大到一定角度时,船舶肩部与冰体棱角碰撞,由于肩部结构强度不及艏部,船体会被冰体棱角迅速划破,碰撞力峰值逐渐减小,而且产生碰撞力峰值的时刻也逐渐提前。因此,船-冰碰撞的碰撞力峰值随着撞击角度的增大而减小,碰撞力峰值产生的时刻也随之提前。

从图9可以看出不同撞击角度下冰体与船体结构损伤的形态与特征。当碰撞角度为0°时,冰体被撞击并伴随着大量的破碎,球鼻艏的碰撞区域被冰体挤压产生了较明显的变形。当碰撞角度为60°时,船首与冰体正面的碰撞接触面积减小,此时船首外板被冰体棱角划割,由于肩部结构强度不及艏部,船体会被冰体棱角迅速划破,对船体内结构造成破坏,在这种碰撞情形下,船体的损伤比较严重。因此,在冰区航行的船舶与冰山碰撞不可避免时,应避开冰体棱角。从冰体的破碎情况可知,碰撞角度越大,冰体破碎越小,同时冰体与船首的接触位置随碰撞角度的改变而发生变化。

图8 不同角度下碰撞力曲线对比图

图9 不同撞击角度下船首外板和冰体的变形损伤及应力分布图

3 结 论

(1)在船-冰碰撞过程中碰撞力均呈现明显的非线性,碰撞力的波动较连续且碰撞力增长快速。船-冰碰撞的碰撞力峰值随冰体质量的增大而增大,碰撞力峰值产生的时刻也随之提前;船-冰碰撞的碰撞力峰值随着撞击的增大而减小,碰撞力峰值产生的时刻也随之提前。

(2)在船-冰碰撞过程中,船体结构与冰体损伤具有明显的局部性。随着冰体质量增大,船首结构损伤程度逐渐增大,船首外板和舷侧纵骨产生了巨大的撕拉、弯曲和凹陷。随着撞击角度增大,船首结构凹陷程度减小,但会出现冰体的棱角开始划割船首外板的情况。

(3)冰体尺寸越小,撞击角度越小,冰体损伤越严重。

猜你喜欢
船首外板船体
船体行驶过程中的压力监测方法
不同类型层冰载荷作用下船首结构响应研究
侧围外板转角深拉伸起皱缺陷研究
侧围外板尾灯处排料困难的解决方案
超大型FPSO火炬塔及船体基座设计
汽车侧围外板A柱起皱缺陷分析及处理方法
航行状态下三体船砰击的动力响应
船首形状对船-冰碰撞性能的影响研究
船体剖面剪流计算中闭室搜索算法
水下爆炸气泡作用下船体总纵强度估算方法