留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

冰碰载荷下船体板弹塑性动力响应分析

蔡伟 朱凌

蔡伟, 朱凌. 冰碰载荷下船体板弹塑性动力响应分析[J]. 中国舰船研究, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
引用本文: 蔡伟, 朱凌. 冰碰载荷下船体板弹塑性动力响应分析[J]. 中国舰船研究, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
CAI W, ZHU L. Elastic-plastic responses of ship plates under ice impact loading[J]. Chinese Journal of Ship Research, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
Citation: CAI W, ZHU L. Elastic-plastic responses of ship plates under ice impact loading[J]. Chinese Journal of Ship Research, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804

冰碰载荷下船体板弹塑性动力响应分析

doi: 10.19693/j.issn.1673-3185.01804
基金项目: 中央高校基本科研业务费资助项目(2019YB013)
详细信息
    作者简介:

    蔡伟,男,1996年生,博士生。研究方向:冰区船舶与海洋结构安全性。E-mail:Caiwei199696@163.com

    朱凌,男,1962年生,博士,教授,博士生导师。研究方向:冰区船舶与海洋结构安全性与完整性。E-mail:lingzhu@whut.edu.cn

    通讯作者:

    朱凌

  • 中图分类号: U663.2

Elastic-plastic responses of ship plates under ice impact loading

  • 摘要:   目的  北极航道在远洋运输中发挥着越来越重要的作用。极地船舶在北极航行过程中不可避免地会与浮冰发生碰撞,从而使结构产生塑性变形损伤,严重影响船舶的航行安全。为了研究浮冰碰撞载荷作用下船体板所产生的塑性变形损伤和能量吸收情况,开展楔形冰与船体板碰撞模型实验研究。  方法  首先,分析船体板的塑性变形损伤以及冰体破碎失效情况;然后,总结并对比现有不同冰体材料模型下冰撞刚性墙的计算结果,揭示不同冰材料模型下计算结果的差异;最后,采用数值模拟方法研究冰碰撞载荷作用下固支矩形板的塑性变形和碰撞力时间历程,并与模型实验结果进行对比验证。  结果  结果显示,不同冰体材料模型的数值模拟结果其差异性和离散性较大;采用了各向同性弹脆性冰数值材料的计算模型其数值计算结果与实验结果吻合较好;在碰撞条件相同冰体与刚体碰撞时,冰体碰撞下船体板的最终塑性变形值和能量吸收值比刚体碰撞下的小。  结论  研究成果可为极地船舶结构安全评估方法提供参考建议。
  • 图  1  水平冲击实验机

    Figure  1.  Horizontal impact experimental facility

    图  2  冰体与船体板碰撞模型实验装置

    Figure  2.  Model test device for ice-plate collision

    图  3  冰体尺寸

    Figure  3.  Ice model dimensions

    图  4  冰体模型

    Figure  4.  Ice model

    图  5  冰体破坏过程[12]

    Figure  5.  Ice failure process[12]

    图  6  板的塑性变形损伤[12]

    Figure  6.  Plastic deformation damage of the plate[12]

    图  7  圆台冰体冲击刚性墙示意图

    Figure  7.  Schematic of conical ice indenter collision with a rigid wall

    图  8  有限元网格划分

    Figure  8.  Division of finite element mesh

    图  9  不同冰体材料模型的碰撞力–时间曲线

    Figure  9.  Collision force -time curves of different ice material models

    图  10  不同冰体材料模型和规范推荐的压强–接触面积曲线

    Figure  10.  Pressure-area curves recommended by different ice material models and some standards

    图  11  有限元计算模型

    Figure  11.  Finite element model

    图  12  冰体与船体板模型碰撞过程[12]

    Figure  12.  Collision process of ice and plate model[12]

    图  13  碰撞力时间历程[12]

    Figure  13.  Time history of collision force[12]

    图  14  板的最终中心变形图[12]

    Figure  14.  Final deformations of plate center line[12]

    图  15  事故极限状态设计(韧性设计、折中设计和强度设计)

    Figure  15.  Accident limit state design (ductile design, shared-energy design, strength design)

    图  16  冰体与船体板碰撞下的能量耗散

    Figure  16.  Energy dissipate when the ice collides with the plate

    图  17  冰体/刚体碰撞下船体板的中心点变形时间历程[13]

    Figure  17.  Center point's deformation time history of the hull plate under ice/steel wedge collision[13]

    图  18  冰体/刚体碰撞下船体板能量吸收时间历程[13]

    Figure  18.  Energy absorption time history of hull plate under ice/steel wedge collision[13]

    表  1  冰模型的材料参数(*MAT_013)[12]

    Table  1.   Material parameters of ice model (*MAT_013)[12]

    参数数值
    材料密度ρ/(kg•m-3)893
    剪切模量G/GPa2.20
    屈服应力σy/MPa2.04
    塑性硬化模量Etan/GPa4.26
    体积模量K/GPa5.26
    失效压力Pt/MPa−4.00
    塑性失效应变0.002
    下载: 导出CSV
  • ZHU L, JAMES P, ZHANG S M. Statistics and damage assessment of ship grounding[J]. Marine Structures, 2002, 15(4/5): 515–530.
    FERRARI N, RIZZUTO E, PRESTILEO A. Modelling ice characteristics in iceberg-ship collision analyses[M]//Towards Green Marine Technology and Transport. London: CRC Press, 2015: 305–316.
    SU B, RISKA K, MOAN T. Numerical study of ice-induced loads on ship hulls[J]. Marine Structures, 2011, 24(2): 132–152. doi:  10.1016/j.marstruc.2011.02.008
    ABRAHAM J. Plastic response of ship structure subject to ice loading[D]. Canada: Memorial University of Newfoundland, 2008.
    LIU Z H, AMDAHL J, LØSET S. Plasticity based material modelling of ice and its application to ship-iceberg impacts[J]. Cold Regions Science and Technology, 2011, 65(3): 326–334. doi:  10.1016/j.coldregions.2010.10.005
    LIU Z H, AMDAHL J, LØSET S. Integrated numerical analysis of an iceberg collision with a foreship structure[J]. Marine Structures, 2011, 24(4): 377–395. doi:  10.1016/j.marstruc.2011.05.004
    张建, 万正权, 黄进浩, 等. 舷侧板架与冰体碰撞数值仿真及模型试验研究[J]. 船舶力学, 2014, 18(4): 424–433. doi:  10.3969/j.issn.1007-7294.2014.04.010

    ZHANG J, WAN Z Q, HUNAG J H, et al. Research on numerical simulation and model test of collision between side grillage and icebergs[J]. Journal of Ship Mechanics, 2014, 18(4): 424–433 (in Chinese). doi:  10.3969/j.issn.1007-7294.2014.04.010
    GAO Y, HU Z Q, RINGSBERG J W, et al. An elastic-plastic ice material model for ship-iceberg collision simulations[J]. Ocean Engineering, 2015, 102: 27–39. doi:  10.1016/j.oceaneng.2015.04.047
    王健伟, 邹早建. 基于非线性有限元法的船舶−冰层碰撞结构响应研究[J]. 振动与冲击, 2015, 34(23): 125–130.

    WANG J W, ZOU Z J. Ship's structural response during its collision with level ice based on nonlinear finite element method[J]. Journal of Vibration and Shock, 2015, 34(23): 125–130 (in Chinese).
    ZHU L, YAN M J, YU T X, et al. Dynamic analysis of ship plates under ice floes impact[C]//International Conference on Collision and Grounding of Ships and Offshore Structures. Ulsan:[s.n.], 2016.
    闫孟娇, 朱凌. 浮冰碰撞载荷作用下船体板结构动力响应研究[J]. 武汉理工大学学报(交通科学与工程版), 2017, 41(2): 268–272.

    YAN M J, ZHU L. Dynamic analysis of ship plates under ice floes impact[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2017, 41(2): 268–272 (in Chinese).
    闫孟娇. 冰碰载荷作用下船体板的冲击动力响应研究[D]. 武汉: 武汉理工大学, 2017.

    YAN M J. Dynamic response of ship plates impacted by ice floes[D]. Wuhan: Wuhan University of Technology, 2017 (in Chinese).
    ZHU L, QIU X M, CHEN M S, et al. Simplified ship-ice collision numerical simulations[C]//The 26th International Ocean and Polar Engineering Conference. Rhodes, Greece: International Society of Offshore and Polar Engineers, 2016: 1192-1196.
    ZHU L, SHI S Y, YU T X. A new ice load-response model for structural design of ice classed ships[C]//The 25th International Ocean and Polar Engineering Conference. Hawaii, USA:[s.n.], 2015.
    ZHU L, CAI W, CHEN M S, et al. Dynamic analysis of ship plates under repeated ice floes impacts based on a simplified ship-ice collision model[C]//The 28th International Ocean and Polar Engineering Conference. Sapporo, Japan: International Society of Offshore and Polar Engineers, 2018.
    吕保达, 刘敬喜, 解德, 等. 基于UMAT的冰−结构相互作用数值仿真[J]. 中国舰船研究, 2015, 10(1): 39–45. doi:  10.3969/j.issn.1673-3185.2015.01.006

    LYU B D, LIU J X, XIE D, et al. Numerical simulation of ice-structure interaction based on UMAT[J]. Chinese Journal of Ship Research, 2015, 10(1): 39–45 (in Chinese). doi:  10.3969/j.issn.1673-3185.2015.01.006
    GAGNON R E. A numerical model of ice crushing using a foam analogue[J]. Cold Regions Science and Technology, 2011, 65(3): 335–350. doi:  10.1016/j.coldregions.2010.11.004
    KIM H, DALEY C, COLBOURNE B. A numerical model for ice crushing on concave surfaces[J]. Ocean Engineering, 2015, 106: 289–297. doi:  10.1016/j.oceaneng.2015.07.020
    VON BOCK UND POLACH R, EHLERS S. Model scale ice—Part B: numerical model[J]. Cold Regions Science and Technology, 2013, 94: 53–60. doi:  10.1016/j.coldregions.2013.06.009
    International Standardizations Organization. Petroleum and natural gas industries-arctic offshore structures (draft): ISO/CD 19906[S]. 1st edition. Geneva, Switzerland:[s.n.], 2010: 434.
    API. American petroleum institute recommended practice for planning, designing, and constructing structures and pipelines for arctic conditions[S]. Washington, DC: API, 1995.
  • [1] 王少博, 张英俊, 胡鑫2021: 考虑船位预测不确定性的船舶碰撞危险度计算方法 , 中国舰船研究. doi: 10.19693/j.issn.1673-3185.01914
    [2] 张新宇, 曹俊伟, 骆伟, 张敏2019: 双壳船体结构在楔形物撞击下的损伤特性试验及数值仿真分析 , 中国舰船研究, 14, 89-94. doi: 10.19693/j.issn.1673-3185.01452
    [3] 张斐, 张春辉, 张磊, 王志军, 周春桂2019: 多次水下爆炸作用下钢板动态响应数值模拟 , 中国舰船研究, 14, 122-129. doi: 10.19693/j.issn.1673-3185.01431
    [4] 马坤, 胡高源, 梁舒凡2019: 极地船舶积冰计算及其对参数横摇的影响分析 , 中国舰船研究, 14, 61-65. doi: 10.19693/j.issn.1673-3185.01572
    [5] 李春鹏, 张攀, 刘均, 程远胜2018: 空爆载荷下功能梯度泡沫铝夹层板动响应数值仿真 , 中国舰船研究, 13, 77-84. doi: 10.19693/j.issn.1673-3185.01172
    [6] 董良雄, 杨意, 高军凯, 龚雅萍2017: 基于船舶艉轴-油膜-艉部结构系统的碰撞载荷响应研究 , 中国舰船研究, 12, 122-127. doi: 10.3969/j.issn.1673-3185.2017.01.018
    [7] 刘俊峰, 胡志强2017: 船舶碰撞机理三维解析法实现及恢复系数研究 , 中国舰船研究, 12, 84-91. doi: 10.3969/j.issn.1673-3185.2017.02.011
    [8] 段菲, 张利军, 陈鸽, 姜海宁, 张琪2017: 基于多目标优化算法NSGA II的极地穿梭油轮型线设计 , 中国舰船研究, 12, 66-72. doi: 10.3969/j.issn.1673-3185.2017.06.010
    [9] 谢强, 陈海龙, 章继峰2017: 极地航行船舶及海洋平台防冰和除冰技术研究进展 , 中国舰船研究, 12, 45-53. doi: 10.3969/j.issn.1673-3185.2017.01.008
    [10] 吴林杰, 朱锡, 侯海量, 陈长海2016: 冲击载荷下舰船夹芯舱室结构动态响应及吸能特性 , 中国舰船研究, 11, 70-76,96. doi: 10.3969/j.issn.1673-3185.2016.06.011
    [11] 刘敬喜, 崔濛, 龚榆峰2015: 船舶碰撞仿真失效准则比较 , 中国舰船研究, 10, 79-85. doi: 10.3969/j.issn.1673-3185.2015.04.012
    [12] 葛亮, 田正东, 袁利毫2015: 气垫船碰撞冲击动力学响应三维数值模拟 , 中国舰船研究, 10, 57-62. doi: 10.3969/j.issn.1673-3185.2015.03.010
    [13] 孙登成, 郭君, 吴铭, 姬秀滨, 祝祥刚, 陈舸2014: 潜艇水下碰撞特性数值研究 , 中国舰船研究, 9, 105-110. doi: 10.3969/j.issn.1673-3185.2014.01.015
    [14] 胡志强, 朱旻, 岑松2013: 导管架平台立柱结构的抗撞性能分析 , 中国舰船研究, 8, 54-63. doi: 10.3969/j.issn.1673-3185.2013.01.009
    [15] 宋波2011: 船舶能量管理系统设计研究 , 中国舰船研究, 6, 93-97. doi: 10.3969/j.issn.1673-3185.2011.02.020
    [16] 梅志远, 李卓, 吕岩松2011: 单/双壳体潜艇结构耐撞剩余强度特性研究 , 中国舰船研究, 6, 7-10. doi: 10.3969/j.issn.16733185.2011.05.002
    [17] 张艺凡, 宗智, 张文鹏2011: 水下冲击波作用下结构损伤的数值预报误差分析 , 中国舰船研究, 6, 38-44. doi: 10.3969/j.issn.1673-3185.2011.06.008
    [18] 舒雄, 朱锡, 张振华, 牟金磊2010: 水下爆炸气泡能载荷对单层板架塑性变形的影响研究 , 中国舰船研究, 5, 12-15. doi: 10.3969/j.issn.1673-3185.2010.04.003
    [19] 张维俊2009: 船舶结构碰撞相似律的数值研究 , 中国舰船研究, 4, 38-41,61. doi: 10.3969/j.issn.1673-3185.2009.03.008
    [20] 肖桃云, 王涛2008: 三峡库区船舶碰撞能量损失方法的研究 , 中国舰船研究, 3, 25-31. doi: 10.3969/j.issn.1673-3185.2008.05.006
  • 加载中
图(18) / 表 (1)
计量
  • 文章访问数:  48
  • HTML全文浏览量:  17
  • PDF下载量:  822
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-21
  • 修回日期:  2020-01-10
  • 网络出版日期:  2020-07-18
  • 刊出日期:  2020-12-30

冰碰载荷下船体板弹塑性动力响应分析

doi: 10.19693/j.issn.1673-3185.01804
    基金项目:  中央高校基本科研业务费资助项目(2019YB013)
    作者简介:

    蔡伟,男,1996年生,博士生。研究方向:冰区船舶与海洋结构安全性。E-mail:Caiwei199696@163.com

    朱凌,男,1962年生,博士,教授,博士生导师。研究方向:冰区船舶与海洋结构安全性与完整性。E-mail:lingzhu@whut.edu.cn

    通讯作者: 朱凌
  • 中图分类号: U663.2

摘要:   目的  北极航道在远洋运输中发挥着越来越重要的作用。极地船舶在北极航行过程中不可避免地会与浮冰发生碰撞,从而使结构产生塑性变形损伤,严重影响船舶的航行安全。为了研究浮冰碰撞载荷作用下船体板所产生的塑性变形损伤和能量吸收情况,开展楔形冰与船体板碰撞模型实验研究。  方法  首先,分析船体板的塑性变形损伤以及冰体破碎失效情况;然后,总结并对比现有不同冰体材料模型下冰撞刚性墙的计算结果,揭示不同冰材料模型下计算结果的差异;最后,采用数值模拟方法研究冰碰撞载荷作用下固支矩形板的塑性变形和碰撞力时间历程,并与模型实验结果进行对比验证。  结果  结果显示,不同冰体材料模型的数值模拟结果其差异性和离散性较大;采用了各向同性弹脆性冰数值材料的计算模型其数值计算结果与实验结果吻合较好;在碰撞条件相同冰体与刚体碰撞时,冰体碰撞下船体板的最终塑性变形值和能量吸收值比刚体碰撞下的小。  结论  研究成果可为极地船舶结构安全评估方法提供参考建议。

English Abstract

蔡伟, 朱凌. 冰碰载荷下船体板弹塑性动力响应分析[J]. 中国舰船研究, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
引用本文: 蔡伟, 朱凌. 冰碰载荷下船体板弹塑性动力响应分析[J]. 中国舰船研究, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
CAI W, ZHU L. Elastic-plastic responses of ship plates under ice impact loading[J]. Chinese Journal of Ship Research, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
Citation: CAI W, ZHU L. Elastic-plastic responses of ship plates under ice impact loading[J]. Chinese Journal of Ship Research, 2020, 15(6): 155–161 doi:  10.19693/j.issn.1673-3185.01804
    • 随着全球气候变暖,北极冰川加快融化,导致浮冰,特别是冰山逐渐增多,给航行于冰区的船舶带来了很大的困难以及许多无法预测的威胁。特别是,水下体积约占90%的自由漂浮冰山,大多带有尖锐的边角,船舶若不小心触碰到,就有可能导致船身损伤、货物泄漏及环境污染等问题,甚至还会造成沉船事故。此外,北极航道上还大量散布着的浮冰,极地船舶在航行过程中不可避免地会遭受到浮冰的反复碰撞,致使船体结构产生塑性变形损伤,从而给船舶航行带来极大的困难和潜在威胁。根据劳氏船级社每年对全球船体破坏事故的统计,由船舶碰撞和搁浅引起的船舶结构损伤占船舶总事故的1/3[1],其中,浮冰与船舶碰撞事故屡见不鲜。1912年4月14日,被认为是当时全球最大的“永不沉没”的“泰坦尼克”(Titanic)号邮轮,在北大西洋的纽芬兰海域因船体右舷撞到了冰山,致使艏部水线撞开了一条长约93 m的裂口,水密舱划破全部进水,最终导致船舶沉没,伤亡人数约1 517人。1959年,丹麦的“赫脱夫特”号海轮在航行过程中因与冰山相撞,导致船体破裂,船上近百人死亡。1996年,Reduta Ordona号散货船在与冰山相撞后船舶结构严重受损。2015年,俄罗斯大型拖网渔船“远东”号与浮冰相撞,导致船舶沉没,船上56人死亡,13人失踪。2019年1月,“雪龙”号破冰船在执行中国第35次南极考察任务期间与冰山发生碰撞,导致船艏桅杆被撞倒,部分舷墙受损,给考察任务造成了严重影响。在这些事故案例中,极地船舶与浮冰或冰山相撞时,其结构的塑性变形损伤程度评估显得尤为重要。

      国内外学者提出了许多冰体理想化材料模型,例如各向同性弹性断裂模型、弹塑性材料模型,以及可压碎泡沫型材料模型等数值模型,用来分析浮冰碰撞下船体结构的塑性动力响应。Ferrari等[2]采用NORSOK方法分析了船–冰碰撞过程中的能量耗散情况。Su等[3]应用数值模型对作用于船体总体和局部的冰载荷予以了估算。Abraham[4]应用有限元方法研究了横向冰载荷作用下船舶结构的塑性响应,并采用冰载荷作用下的板架结构实验对数值模拟结果进行了验证。Liu等[5-6]对船舶与冰山碰撞的外部及内部动力学问题进行了较系统的研究,提出了考虑船体几何形状构造和摩擦系数对能量损耗影响的计算公式。张健等[7]进行了船–冰碰撞数值仿真,研究了冰体形状、冰体质量、碰撞速度等因素对船–冰碰撞的影响机理。Gao等[8]开展了船舶与冰山碰撞有限元计算,研究了船–冰碰撞过程中的碰撞力和能量耗散状态。王健伟和邹早建[9]对船–冰碰撞进行数值模拟,研究了船舶在不同速度与不同厚度冰层碰撞下的动态结构响应。朱凌等[10-15]开展了船体板与浮冰碰撞的有限元仿真和实验研究,分析了船体板–浮冰碰撞参数对结构动态响应与塑性变形的影响规律,并基于刚性块、实际冰与船体板碰撞数值仿真计算结果,提出了一种可以保守估计局部冰载荷的计算方法,经对比分析,得到了考虑冰体破碎引起的船体板变形损伤折减比与能量吸收折减比,提出了浮冰碰撞下船体板的塑性设计方法。吕保达等[16]基于ABAQUS用户自定义子程序UMAT,采用已有的冰体弹性失效准则,对冰体与结构的相互作用进行了数值模拟。

      目前,国内外针对冰体与船舶结构碰撞实验的研究工作相对较少,主要是以数值模拟方法为主,相关的数值模拟结果缺乏实验验证。因此,本文将在总结前期楔形冰与船体板碰撞实验工作的基础上,首先研究船体板的塑性变形损伤以及冰体破碎失效情况,并验证数值方法的可靠性;然后利用数值模拟方法分析冰碰载荷作用下船体板的塑性变形损伤和能量吸收情况,并与刚体碰撞下船体板的弹塑性响应进行对比,研究刚体碰撞和冰体碰撞下船体板的塑性动力响应差异。

    • 船体结构–冰体碰撞实验中采用的主要设备为武汉理工大学自行设计的冲击与流固耦合实验室水平冲击实验机,如图1所示。该实验装置包括5个部分,分别为基座、夹具系统、冲击车、弧形滑轨以及牵引系统。在碰撞实验过程中,所用到的实验仪器有激光位移传感器、加速度传感器、高速摄影机、M+P数据采集器及TMR数据采集系统等,相关的实验装置和仪器设备如图2所示。在实际情况中,冰体的形状各种各样,本文中的冰体主要简化为了楔形冰,其尺寸和形状分别如图3图4所示。图3中:a1b1b2d1d2为楔形冰的主要几何尺寸参数;h为楔形冰的初始高度(图中参数单位为mm)。冰体被装载在冲击车内,通过滑轨给予冲击速度,以方便实现撞击过程。制冰前,先做好冰体模具,然后逐步加入淡水,在低温下(−15 ℃)冻结成冰。在每次冰体冲击实验中,冰体尺寸被严格控制成图4所示的尺寸,此外,每次的制冰流程和步骤(例如加水时间、加水量、冷冻时间和温度等)都受到严格控制,以保证每次冷冻成形的冰体没有太大差异。冰体水平冲击矩形板的实验过程如下:首先,在实验前测量冰体质量;接着,将冰体前端与试板中心处对齐;然后,拉高小车,让其自由落下,以给予冰体一个初速度V0;最后,采集相关传感器数据,并测量矩形板变形和剩余冰体质量。实验选取的试板尺寸为1 200 mm×400 mm×1 mm,具体的实验过程和制冰细节见文献[11-12]。

      图  1  水平冲击实验机

      Figure 1.  Horizontal impact experimental facility

      图  2  冰体与船体板碰撞模型实验装置

      Figure 2.  Model test device for ice-plate collision

      图  3  冰体尺寸

      Figure 3.  Ice model dimensions

      图  4  冰体模型

      Figure 4.  Ice model

    • 在冰体碰撞实验过程中,通过高速摄影设备记录碰撞过程中冰体的破碎情况,如图5所示。在碰撞过程中,冰体以颗粒的形式断裂剥落,并会出现较深的裂纹。楔形冰由于其顶端较尖锐,在与试板碰撞接触的过程中会逐渐破碎,使得楔形冰与试板的接触面积增大。当楔形冰的速度减为0时,楔形冰和试板开始一起回弹,直至楔形冰与试板分离。碰撞结束后试板的永久变形图如图6所示。在试板与冰体碰撞接触区域,中心的变形最大,两端逐渐变小。

      图  5  冰体破坏过程[12]

      Figure 5.  Ice failure process[12]

      图  6  板的塑性变形损伤[12]

      Figure 6.  Plastic deformation damage of the plate[12]

    • 为了对比现有的一些不同冰材料数值有限元计算模型,例如弹塑性材料模型[5-6]、各向同性弹性断裂模型[10-15]、可压碎泡沫型材料模型[17-18]以及损伤材料模型[19]等,按照上述冰体材料模型对冰体进行定义,利用ANSYS/LS-DYNA商业计算软件,并采用上述冰体数值材料模型对圆台冰体冲击刚性墙进行数值仿真。图7中:d1d2分别为圆台冰体的上、下底面直径,d1=0.2 m,d2=10 m;l0为圆台冰体的初始长度,l0=5 m;x为冲击过程中冰体的破碎长度;y为冲击过程中圆台冰体与刚性墙接触圆面的直径。冲击速度V=10 m/s。有限元网格划分如图8所示,冰体网格单元长度为100 mm。

      图  7  圆台冰体冲击刚性墙示意图

      Figure 7.  Schematic of conical ice indenter collision with a rigid wall

      图  8  有限元网格划分

      Figure 8.  Division of finite element mesh

      刚性墙圆台撞击数值计算的主要目的是,测得一个圆锥形的冰在撞击刚性墙的过程中,撞击面平均压强P与冲击过程中不断变化冰体的撞击接触面积A的关系,可以与P-A标准理论曲线(ISO/CD 19906)[20]等所提出的P-A曲线进行对比。P-A标准理论曲线为

      $${P_{\rm {ISO}} } = 7.4{A^{ - 0.7}}$$ (1)

      此外,美国石油协会API/CSA推荐的压力–面积曲线[21]

      $${P_{\rm {API}}} = 8.1{A^{ - 0.572}}$$ (2)

      图9对比了运用现有一些冰体材料模型计算得到的碰撞力–时间曲线。由图可看出,不同冰体材料模型计算出来的结果差异较显著:材料模型(2)的碰撞时间最短,材料模型(1)的碰撞时间较长;材料模型(2)和材料模型(4)的碰撞力较大。可见材料的碰撞时间长短和碰撞力大小取决于冰体材料的刚度,冰体材料越强,碰撞力越大,碰撞时间越短;反之,则碰撞力越小,碰撞时间越长。此外,基于文献[5, 8]中冰体撞击刚性墙的计算结果以及图9中4种模型的计算结果,图10给出了不同冰体材料模型下计算得到的压强–接触面积曲线,以及ISO和API推荐的压强–接触面积曲线。不同冰体材料模型计算出来的结果离散性较大,但其结果与ISO和API推荐的压强–接触面积的值总体上比较接近。

      图  9  不同冰体材料模型的碰撞力–时间曲线

      Figure 9.  Collision force -time curves of different ice material models

      图  10  不同冰体材料模型和规范推荐的压强–接触面积曲线

      Figure 10.  Pressure-area curves recommended by different ice material models and some standards

    • 在该冰体与船体板碰撞实验模型中,试板的材料参数通过标准的材料拉伸实验得到,冰体的材料参数参考了淡水圆台冰单轴压缩实验数据结果以及现有的各向同性弹脆性有限元数值材料模型(*MAT_013),相关的冰材料参数和试板参数可参见文献[12],其中冰模型的材料参数如表1所示。有限元计算模型如图11所示,其中试板尺寸为1 200 mm×400 mm×1 mm,采用壳单元Shell 163模拟;冰体总质量为85.8 kg,其尺寸如图4所示,采用实体单元Solid 164模拟。楔形冰与试板碰撞的初速度为2.58 m/s,试板四周进行刚性固定,冰体仅保持碰撞速度方向的自由度。楔形冰与试板的接触方式采用面–面接触中的侵蚀接触,即Eroding surface-to-surface contact。楔形冰与船体板模型碰撞过程的有限元计算结果如图12所示,其碰撞力时间历程及碰撞后试板最终塑性变形的实验结果和数值计算结果如图13图14所示,结果显示两者吻合较好。

      表 1  冰模型的材料参数(*MAT_013)[12]

      Table 1.  Material parameters of ice model (*MAT_013)[12]

      参数数值
      材料密度ρ/(kg•m-3)893
      剪切模量G/GPa2.20
      屈服应力σy/MPa2.04
      塑性硬化模量Etan/GPa4.26
      体积模量K/GPa5.26
      失效压力Pt/MPa−4.00
      塑性失效应变0.002

      图  11  有限元计算模型

      Figure 11.  Finite element model

      图  12  冰体与船体板模型碰撞过程[12]

      Figure 12.  Collision process of ice and plate model[12]

      图  13  碰撞力时间历程[12]

      Figure 13.  Time history of collision force[12]

      图  14  板的最终中心变形图[12]

      Figure 14.  Final deformations of plate center line[12]

    • 对船舶结构在进行事故状态(accidental limit state,ALS)设计时,通常有3种设计方法:韧性设计、折中设计和强度设计,如图15所示。当船体结构强度明显弱于冰体强度时,船体结构会通过变形破坏来吸收大部分的碰撞能量,此时的结构设计即为韧性设计;相反地,当船体结构强度明显强于冰体强度时,此时船体结构几乎不发生变形,吸收的能量较少,此时的结构设计为强度设计;介于两者之间的即为折中设计。

      图  15  事故极限状态设计(韧性设计、折中设计和强度设计)

      Figure 15.  Accident limit state design (ductile design, shared-energy design, strength design)

      图16所示为楔形冰与船体板碰撞过程的能量转换算例,即冰体的初始动能Etotal转换为碰撞过程中板变形所吸收的能量Ep和冰体变形破坏所耗散的能量Ei;当碰撞速度变为0,即冰体动能为0时,此时板不再吸收能量,并释放弹性能转换为冰体的回弹动能Er

      图  16  冰体与船体板碰撞下的能量耗散

      Figure 16.  Energy dissipate when the ice collides with the plate

      图17图18给出了冰体与刚体碰撞下船体板的弹塑性动力响应对比。图中:$\varepsilon = {W_ {\rm {c}}}/H$,为板中心点变形与板厚的比值;$\beta = {E_{ {\rm {i - plate}}}}/\sum E$,为板吸收的能量与碰撞总能量的比值。有限元计算模型如图11所示。试板尺寸为250 mm×200 mm×1.65 mm,楔形冰质量为23.3 kg,冰体厚度为100 mm,具体的模型参数和计算前提条件参见文献[13]。从图中可以看出,在保证冰体和刚性块冲头的前段楔形体几何形状、质量和速度相同的条件下,由于冰体在碰撞过程中出现了变形破坏的特征,冰体与刚体碰撞下的矩形板弹塑性动力响应有着较大差异。在碰撞初始过程中,由于冰体破碎特征不明显,2种碰撞情况下板的变形值几乎相等,但随着冰体失效破坏特征的显著增加,两者的变形出现了较明显的差异,刚性块碰撞下板的最大变形和最终塑性变形值比冰体碰撞下板的最大变形和最终塑性变形值要大,但其碰撞时间较短。在能量吸收方面,冰体碰撞下板变形所吸收的能量明显小于刚体碰撞下板变形所吸收的能量。

      图  17  冰体/刚体碰撞下船体板的中心点变形时间历程[13]

      Figure 17.  Center point's deformation time history of the hull plate under ice/steel wedge collision[13]

      图  18  冰体/刚体碰撞下船体板能量吸收时间历程[13]

      Figure 18.  Energy absorption time history of hull plate under ice/steel wedge collision[13]

    • 本文对楔形冰与船体板模型的碰撞试验工作予以了总结,研究了冰体碰撞下船体板的弹塑性动力响应,分析了冰体失效破坏过程和板的塑性变形损伤情况。此外,还采用各向同性弹脆性冰数值材料模型,对冰体碰撞下船体板的弹塑性动力响应进行了有限元仿真,并将其碰撞力时程曲线和塑性变形计算结果与实验结果进行了对比,结果显示两者吻合较好,验证了该数值方法的可靠性。本文还对现有的一些冰材料数值计算模型予以了总结,并对不同冰体材料模型计算出的冰撞击刚性体的计算结果进行了对比分析,结果表明:不同冰材料模型下的计算结果差异较大,离散性也较大,但其数值模型的计算结果与ISO和API推荐的压强–接触面积值总体上比较接近,可以用于船–冰碰撞动力响应研究。最后,分析了相同碰撞条件及冰体和刚体碰撞下船体结构的弹塑性动力响应,研究了两者间板塑性变形和能量吸收的差异。两者结果的对比显示,由于冰体在碰撞过程中会因变形破坏而耗散能量,故在冰体碰撞下船体板的最终塑性变形值和能量吸收值比刚体碰撞下的小。

      本文的冰体模型采用的是淡水冰模型,在今后的研究中,需要将其拓展到海冰模型中。此外,在船体板和浮冰碰撞实验以及数值仿真中,本文忽略了船–冰碰撞时的真实流体环境影响,而在实际情况中,船舶是在水中与浮冰发生碰撞的,因此,复杂船舶结构与浮冰碰撞中的流固耦合问题还有待深入研究。

参考文献 (21)

目录

    /

    返回文章
    返回