Processing math: 0%

基于离散元法的极地船舶斜航冰载荷特性分析

杨霄楠, 施海天, 柴威, 吕成, 封培元, 吴达

杨霄楠, 施海天, 柴威, 等. 基于离散元法的极地船舶斜航冰载荷特性分析[J]. 中国舰船研究, 2025, 20(X): 1–9. DOI: 10.19693/j.issn.1673-3185.04269
引用本文: 杨霄楠, 施海天, 柴威, 等. 基于离散元法的极地船舶斜航冰载荷特性分析[J]. 中国舰船研究, 2025, 20(X): 1–9. DOI: 10.19693/j.issn.1673-3185.04269
YANG X N, SHI H T, CHAI W, et al. Study on ice load characteristics of polar ships under oblique sailing based on discrete element method[J]. Chinese Journal of Ship Research, 2025, 20(X): 1–9 (in Chinese). DOI: 10.19693/j.issn.1673-3185.04269
Citation: YANG X N, SHI H T, CHAI W, et al. Study on ice load characteristics of polar ships under oblique sailing based on discrete element method[J]. Chinese Journal of Ship Research, 2025, 20(X): 1–9 (in Chinese). DOI: 10.19693/j.issn.1673-3185.04269
杨霄楠, 施海天, 柴威, 等. 基于离散元法的极地船舶斜航冰载荷特性分析[J]. 中国舰船研究, 2025, 20(X): 1–9. CSTR: 32390.14.j.issn.1673-3185.04269
引用本文: 杨霄楠, 施海天, 柴威, 等. 基于离散元法的极地船舶斜航冰载荷特性分析[J]. 中国舰船研究, 2025, 20(X): 1–9. CSTR: 32390.14.j.issn.1673-3185.04269
YANG X N, SHI H T, CHAI W, et al. Study on ice load characteristics of polar ships under oblique sailing based on discrete element method[J]. Chinese Journal of Ship Research, 2025, 20(X): 1–9 (in Chinese). CSTR: 32390.14.j.issn.1673-3185.04269
Citation: YANG X N, SHI H T, CHAI W, et al. Study on ice load characteristics of polar ships under oblique sailing based on discrete element method[J]. Chinese Journal of Ship Research, 2025, 20(X): 1–9 (in Chinese). CSTR: 32390.14.j.issn.1673-3185.04269

基于离散元法的极地船舶斜航冰载荷特性分析

基金项目: 国家自然科学基金资助项目(52201379);中央高校基本科研业务费专项资金资助项目(3120624109);上海市船舶工程重点实验室开放基金资助项目(SE202305)
详细信息
    作者简介:

    柴威,男,1988年生,博士,教授,博士生导师。研究方向:极地船舶。Email:chaiwei@whut.edu.cn

    通讯作者:

    柴威

  • 中图分类号: U674.81

Study on ice load characteristics of polar ships under oblique sailing based on discrete element method

知识共享许可协议
基于离散元法的极地船舶斜航冰载荷特性分析杨霄楠,采用知识共享署名4.0国际许可协议进行许可。
  • 摘要:
    目的 

    旨在开展斜航作业下船舶冰载荷的准确预报和冰载荷特性研究。

    方法 

    以某极地船型为研究对象,应用离散元法对碎冰区船舶斜航状态下的冰载荷进行预报。首先,针对目标船型建立数值模型,通过对比试验数据验证其准确性;其次,设定不同的斜航角度、航速及冰厚,计算斜航过程中船舶受到的冰载荷;最后,分析不同漂角和冰厚条件下船舶冰载荷的分布特性及变化趋势。

    结果 

    研究显示,随着漂角的增大,船舶所受碎冰阻力及碎冰侧向力呈非线性增长;船模在航速0.6 m/s、海冰密集度70%,冰厚0.014 97 m的条件下,在斜航工况15°漂角下,相比直航工况,碎冰阻力与碎冰侧向力增幅分别达4.25倍和6.04倍;船舶在大于10°漂角工况斜航作业下,提高航速或进入冰厚较大区域需注意碎冰阻力及迎冰侧碎冰侧向力突变可能造成的不利影响。

    结论 

    所做研究可为船舶在极地碎冰区斜航工况下的安全性评估提供数据支持,为斜航工况下极地船舶冰载荷预报与研究提供参考。

    Abstract:
    Objectives 

    As global warming accelerates sea - ice melting, the Arctic region witnesses an increase in ship navigation. The brash ice area, composed of brash ice of various sizes and shapes, is a common operational scenario for polar ships. Understanding the ice - load characteristics of polar ships during oblique navigation in brash ice regions is crucial. It can enhance ship navigation safety in the complex polar marine environment, provide a reference for polar navigation route planning, and fill the gap in the current research that mainly focuses on straight - sailing conditions.

    Methods 

    This study selects a specific polar ship type as the research object and utilizes the Discrete Element Method (DEM) to predict the ice loads on the ship during oblique navigation in the brash ice area. First, a numerical model of the target ship is established. The model parameters include a ship model with a scale ratio of 60, a total length of 2.04 m, a beam of 0.37 m, and a design draft of 0.13 m. The ice particles have a density of 917.0 kg/m³, a Poisson's ratio of 0.3, and other specific properties. The accuracy of the model is verified by comparing with the experimental results of the straight - sailing condition in the literature. Then, different oblique - sailing angles (0° - 15°), speeds (0.6 m/s, 0.7 m/s), and ice thicknesses (0.01167 m, 0.01497 m) are set. The ice - load calculation is carried out based on the momentum conservation equation, angular - momentum conservation equation, and the linear spring contact - force model in the DEM.

    Results 

    The results show that as the drift angle increases, the ice - breaking resistance and lateral force on the ship grow non - linearly. For example, under the conditions of a speed of 0.6 m/s, an ice concentration of 70%, and an ice thickness of 0.01497 m, when the drift angle is 15°, the ice - breaking resistance and lateral force increase by 4.25 times and 6.04 times respectively compared to the straight - sailing condition. In terms of the influence of speed, when the drift angle is in the range of 0° - 10°, the ice - breaking resistance increases slowly, and when it exceeds 10°, it increases significantly. The lateral force also shows a non - linear growth trend, and the influence of speed on the lateral force is more significant than whether the ship is on the ice - facing side. Regarding the influence of ice thickness, when the drift angle is greater than 10°, the ice - breaking resistance and lateral force increase significantly with the increase of ice thickness.

    Conclusions 

    In conclusion, this research provides reliable data support for the safety assessment of ships during oblique navigation in polar brash ice regions. It offers a valuable reference for predicting and studying ice loads on polar ships under such conditions. Ship operators should be cautious when increasing speed or entering thicker ice areas, especially when the drift angle is greater than 10°. This is to avoid potential risks caused by sudden changes in ice - breaking resistance and lateral force, ensuring the safe and stable navigation of polar ships in complex ice - covered waters.

  • 随着全球气候的变暖,海冰加速融化,北极地区船舶航行数量与航行窗口期呈现增长趋势[1]。对于极地航行船舶,由不同大小和形状的破碎冰块组成的浮碎冰区域是其遭遇的典型作业场景之一[2-5]。对浮碎冰区域的船舶冰载荷特性进行研究,有助于提升船舶在极地复杂海洋环境下的航行安全,可为极地航行路线规划提供参考。

    目前,关于碎冰区极地船舶冰载荷预报的研究主要集中在船舶直航工况下,其研究方法主要有经验公式、冰水池试验和数值模拟方法3种。在经验公式方面,Sazonov等[6]提出了一种估算碎冰阻力的公式,该公式是基于早期碎冰航道中预报冰阻力的半经验方法,参考了5种不同的船舶模型试验数据,并对船体参数影响进行了考虑。但该公式更多地关注狭窄航道中碎冰对船舶阻力的影响,其在更广泛的碎冰区域中的应用还有待考证。

    冰水池试验是目前最被认可的一种方法。Jeong等[7]在冰池中开展了IA冰级和IB冰级船舶在碎冰区的阻力试验,在与芬兰−瑞典冰级规则(Finnish-Swedish ice class rules,FSICR)公式进行对比验证其准确性后,开展了不同航道宽度、冰密集度和冰厚对碎冰区阻力影响的试验研究[8]。对于没有冰池试验设施的科研院所,在常规水池中加入非冻结模型冰试验是一种替代选择。Xue等[9]使用非制冷聚丙烯模型冰,开展了船体与碎冰之间的局部冰载荷以及船舶总阻力预报研究。但水池试验存在成本较高以及在冰面上试验次数受限的缺点。

    数值模拟方法是目前研究领域中最常用的方法之一,主要包括有限元法和离散元法等。Yang等[10]基于有限元软件LS–DYNA开展了碎冰阻力预报,并与DuBrovin[11]的经验公式进行了对比。然而,有限元法无法模拟船舶的水阻力,采用该方法开展碎冰区域的船舶冰载荷研究计算精度较低。采用离散元法可以很好地解决碎冰非连续性问题,在碎冰模拟方面具有明显的优势。徐雪松等[12]研究了平板结构在碎冰与水耦合作用下的运动机理,提出了碎冰的4种运动模式:堆积、滑移、翻转和跟随,并探讨了碎冰密集度和平板倾角对平板受碎冰阻力及其运动模式的影响。季顺迎等[13]采用离散元法构建了三维圆盘状碎冰模型,开发了用于模拟海冰与船体结构相互作用的粘弹性接触模型,并对船舶与碎冰之间的相互作用机制进行了系统研究。黄家煜等[14]采用离散元法对碎冰进行建模,分析了不同航速和碎冰环境对船体区域冰压力的影响。Zhong等[15]通过离散元法建立船舶碎冰航行数值模型,研究了不同航速和冰密集度下船体冰压力的分布及其概率特性,并与试验进行了对比验证。王世超等[16]基于离散元法,提出采用HMIGA算法构建碎冰场,并对KCS标准模型进行船型优化,有效降低了碎冰阻力。

    综上所述,国内外学者针对极地船舶碎冰区直航工况冰载荷特性做了大量研究,但斜航作为典型的工况,目前针对该工况下冰载荷特性的研究较少。结合以上现状,本文拟采用离散元法对碎冰区中某极地船在不同漂角下的斜航运动进行数值模拟,系统地开展冰载荷特性预报研究,统计船舶在斜航过程中所受碎冰阻力与侧向力的平均值和标准差,量化分析船舶碎冰区斜航过程中冰载荷的特性规律,为斜航工况下极地船舶冰载荷预报与研究提供参考。

    本文中的数值模拟研究遵循连续性方程和动量守恒方程的要求,使用标准k-ε模型求解湍流问题。利用STAR-CCM+软件离散元法模块的单向耦合和双向耦合模式,可以实现冰与水之间的动量和能量交换。以0.7 m/s航速、0.014 97 m冰厚(实尺度下为0.9 m)、15°漂角工况为例,经验证,双向耦合碎冰阻力平均值为6.80 N,单向耦合碎冰阻力平均值为7.45 N,误差为9.55%,该工况下碎冰的不确定度为9.16%。在相同计算条件下,双向耦合下的计算时间为单向耦合下的2~3倍,因此,为了节省计算效率,采用单向耦合进行计算。

    在船−冰−水耦合作用过程中,模拟碎冰采用基于拉格朗日框架的离散元法,碎冰的运动通常包含平动和旋转2种形式,其运动控制方程基于牛顿第二定律,如下所示[17]

    动量守恒方程:

    m\frac{{{\mathrm{d}}{{\boldsymbol{v}}_{\mathrm{p}}}}}{{{\mathrm{d}}t}} = {{\boldsymbol{F}}_{\mathrm{g}}} + {{\boldsymbol{F}}_{\mathrm{f}}} + {{\boldsymbol{F}}_{\mathrm{c}}} (1)

    角动量守恒方程:

    I\frac{\mathrm{d}{\boldsymbol{\omega }}_{\mathrm{p}}}{\mathrm{d}t}={\boldsymbol{M}}_{\mathrm{f}}+{\boldsymbol{M}}_{\mathrm{c}} (2)

    式中:m为碎冰颗粒的质量;vp为碎冰颗粒的速度;Fg为碎冰颗粒的重力;Ff为碎冰颗粒所受的流体力;Fc为来自其他颗粒或船体的接触力;I为碎冰颗粒的惯性矩;ωp为碎冰颗粒的角速度;Mf为碎冰颗粒所受的扭矩;Mc为来自其他颗粒或船体的扭矩。

    使用线性弹簧模型计算冰−冰或船−冰相互作用的接触力Fc,线性弹簧模型广泛用于碎冰颗粒模型中,以求解接触力[17-19]。在线性弹簧接触模型中,弹簧提供弹性力,阻尼器提供黏性阻尼。结果,接触力Fc可以被分成沿着2个球体接触表面的法向和切向的2个分量,其中法向为沿着2个颗粒中心的连线,切向与之垂直,如图1所示。

    图  1  弹簧−阻尼器接触模型
    Figure  1.  Spring-damper contact model

    法向分量和切向分量可以表示入如下:

    {F_{\mathrm{c}}} = {F_{\mathrm{n}}}{\boldsymbol{n}} + {F_{\mathrm{t}}}{\boldsymbol{\tau}} (3)

    式中:Fn为法向作用力;Ft为切向作用力;nτ分别为法向方向和切向方向的单位向量。

    当计算船−冰接触力时,可以将船舶颗粒视为具有无限直径的颗粒。此时,法向和切向分别垂直于和切向于船体表面。

    本文计算模型选取为某极地船舶,其主要参数如表1所示。围绕船体建立数值水池,在流体域设置中,流体域设置为长7.84L、宽2.40L、高1.00L。如图2所示,在计算区域的高度方向,水区域和空气区域的高度分别为0.543L和0.457L

    表  1  船舶模型主要参数
    Table  1.  Main parameters of ship model
    参数 数值
    缩尺比 60
    总长L/m 2.04
    型宽B/m 0.37
    设计吃水D/m 0.13
    下载: 导出CSV 
    | 显示表格
    图  2  计算域划分图
    Figure  2.  Division of the computational domain

    划分网格时,各漂角网格的尺寸保持不变。通过布尔减运算更新计算域网格。网格划分采用切割体网格。划分完成后,总体网格单元数209万。为确保船体周围流场平滑过渡,设置了6层棱柱层,棱柱层的延伸率为1.2。计算域及船体网格的划分如图3所示。

    图  3  计算域及船体网格划分
    Figure  3.  Division of computational domain and hull mesh

    基于冰块不会发生破碎的假设不涉及海冰模型的弯曲失效。为了节省计算资源,每个冰单元的大小和形状一致,每个碎冰只包含1个正方体单元,碎冰主要参数如表2所示。模拟采用部件喷射器来产生碎冰,喷射方向角度随机的碎冰颗粒。部件喷射器代表从各部件几何形状中提取的一系列喷射点,这些点可以是在部件上定义的点的一个子集,碎冰颗粒从这些点中进行喷射,从而实现造冰过程。

    表  2  碎冰颗粒主要参数
    Table  2.  Main Parameters of the icebroken particle
    参数 数值
    缩尺比 60
    密度/(kg∙m−3) 917.0
    泊松比 0.3
    杨氏模量/Pa 1.67×107
    船−冰摩擦系数 0.138
    冰−冰摩擦系数 0.29
    下载: 导出CSV 
    | 显示表格

    选取文献[10]中的直航工况模型试验结果对本文数值计算模型进行验证。本次验证的 碎冰尺寸与试验中的非冻结冰尺寸一致,碎冰边长为0.067 m(实尺度下为4.0 m),高度为0.014 97 m(实尺度下为0.9 m)。验证选取70%的碎冰密集度,在模型尺度下,分别选取0.5,0.6,0.7和0.8 m/s这4个航速进行敞水阻力与碎冰阻力验证。图4所示为静水及碎冰阻力数值模拟结果与试验数据的对比。由图可知,静水阻力误差小于5%,碎冰阻力误差不超过碎冰的不确定度,二者均具有较高的精度,可以较好地模拟船舶在碎冰区航行时相关的冰载荷特性。

    图  4  静水及碎冰阻力的试验验证
    Figure  4.  Experimental validation of hydrostatic and icebreaking resistance

    本文针对斜航状态下的冰载荷特性进行仿真,为了确保计算域和网格划分的准确性,在保证计算精度和计算效率的情况下,后续的计算均在模型尺度下进行。同时,采用傅汝德及柯西相似定律保证碎冰模型与实际尺度在力学上的相似性[3]

    图5所示,漂角α为船首与x轴的夹角。以航速0.6 m/s、冰厚0.014 97 m工况为例,对船舶在斜航过程中的碎冰整体及局部分布情况进行说明,如图6所示。

    图  5  漂角示意图
    Figure  5.  Schematic diagram of the drift angle
    图  6  斜航状态下的碎冰分布图
    Figure  6.  Distribution of brash ice in oblique sailing

    图6所示,随着漂角的增大,船舶迎冰侧碎冰堆积增加,船−冰接触更加频繁,部分碎冰滑移至船底与船体持续发生相互作用。对比图6(a)和图6(d)可知,随着漂角的增大,船舶尾迹无冰航道的宽度明显增加。

    船舶在0.6 ,0.7 m/s速度下以不同漂角在碎冰区航行时所受的碎冰阻力时历曲线如图7所示。

    图  7  不同速度和漂角下船舶所受碎冰阻力时历曲线
    Figure  7.  Time history curves of brash ice resistance for ships at various drift angles and speeds

    图7可知,碎冰开始与船体发生接触后,由于船−冰碰撞和摩擦的随机性,船舶所受碎冰阻力也呈现出波动性。随着漂角的增大,冰块除了与船体发生摩擦作用外,冰块的翻转、堆积等运动对船体造成的挤压载荷也会持续增大,导致碎冰阻力峰值增加。不同速度下船舶的碎冰阻力平均值及标准差随漂角变化的曲线如图8所示,图中误差棒表示各工况下的阻力均值加上或是减去阻力的标准差。

    图  8  不同漂角下碎冰阻力平均值及标准差的变化曲线
    Figure  8.  Variation curves of the average values and standard deviations of brash ice resistance at different drift angles

    图8所示,船舶在0.6,0.7 m/s速度下的碎冰阻力平均值随着漂角的增大,整体变化趋势一致,均呈非线性增长。当漂角α在0°~10°范围内时,碎冰阻力增长缓慢,接近于线性;当漂角α增至10°~12.5°时,碎冰阻力显著跃升。在小漂角下,碎冰阻力增幅较大,在α = 5°时增幅达28.5%,而在大漂角下其增幅相对较小。在2种航速下,碎冰阻力标准差随着漂角的增大,整体呈增大趋势。在0.6 m/s航速下,当α由7.5°增大至10°时,碎冰不确定度增加,导致标准差迅速增大;在0.7 m/s航速下,当α由10°增加至12.5°时,船−冰接触更加频繁、剧烈,表现为标准差显著增大,船舶碎冰阻力在此突增。当α = 15°时,航速为0.6,0.7 m/s时的碎冰阻力分别为直航工况的4.25倍和3.86倍。因此,斜航状态下船舶所受的碎冰阻力问题值得关注。

    船舶在0.6,0.7 m/s速度下以不同漂角在碎冰区航行时所受的碎冰侧向力时历曲线如图9所示。由图可见,在各漂角下,船舶所受碎冰侧向力均呈现上、下波动的特性,这是因为冰块在与船发生碰撞作用时冰块的运动姿态会发生变化,部分冰块会产生翻转并沿船体表面滑行,通过,还会有部分冰块被波浪排开,导致冰块与船体不断发生加载−卸载,使得碎冰侧向力呈现随机波动的特点。不同漂角下船舶碎冰侧向力平均值的变化曲线如图10所示。

    图  9  不同速度和漂角下船舶所受碎冰侧向力时历曲线
    Figure  9.  Time history curves of brash ice lateral force on ships at various drift angles and speeds
    图  10  船舶碎冰侧向力平均值变化曲线
    Figure  10.  Variation curves of the average value of lateral force of the ship against brash ice

    图10可知,随着漂角的增大,船舶所受碎冰侧向力呈非线性增长趋势。结合图6可知,由于船舶迎冰侧碎冰堆积较多,碎冰−船体接触力更强,因此在相同航速下,迎冰侧的侧向力整体大于非迎冰侧。相较于船舶是否迎冰,航速对船舶所受碎冰侧向力的影响更显著,在0.7 m/s航速下非迎冰侧的侧向力大于0.6 m/s航速下迎冰侧的侧向力。

    在0.6 m/s速度下,α>10°时侧向力显著增大;而在0.7 m/s速度下,α从5°增至7.5°时侧向力略微下降,但在α = 7.5°后也快速上升。不同速度下漂角对侧向力的影响差异明显:当α = 0°时,速度对侧向力的影响较小;当α = 10°时增幅显著,达81.3%。因此,当船舶在α = 10°及以上工况下增加速度时,需关注侧向力突增可能带来的航行风险。

    船舶在冰厚为0.011 67 和0.014 97 m时以不同的漂角在碎冰区航行时所受碎冰阻力时历曲线如图11所示。在不同冰厚下船舶碎冰阻力平均值随漂角的变化如图12所示。

    图  11  不同冰厚和漂角下船舶所受碎冰阻力时历曲线
    Figure  11.  Time history curves of brash ice resistance on ships at various drift angles and ice thicknesses
    图  12  不同漂角下碎冰阻力平均值及标准差变化曲线
    Figure  12.  Variation curves of the average values and standard deviations of brash ice resistance at different drift angles

    图11图12可知,随着漂角的增大,船舶碎冰阻力平均值呈非线性增长;碎冰阻力标准差随着漂角的增大,总体呈增大趋势。当冰厚为0.014 97 m时,α>10°时碎冰阻力显著增加;在较小的冰厚下,小块的碎冰更加不稳定,碎冰阻力标准差更大。当冰厚为0.011 67 m时,α<7.5°时阻力增长较缓,而在α>7.5°时近似呈线性增长。在相同漂角下,不同冰厚对碎冰阻力的影响差异明显:漂角α超过10°后,伴随着冰厚的增加,碎冰阻力显著跃升,最高增幅可达42.9%。因此,需关注较大漂角以及更厚碎冰区对船舶航行的影响。在不同冰厚条件下,当α = 7.5°时2种冰厚标准差均出现了下降情况。

    船舶在冰厚为0.011 67,0.014 97 m时以不同的漂角在碎冰区航行时所受碎冰侧向力时历曲线如图13所示。在不同冰厚下,船舶碎冰侧向力平均值随漂角的变化如图14所示。

    图  13  不同冰厚和漂角下船舶所受碎冰侧向力时历曲线
    Figure  13.  Time history curves of brash ice lateral force on ships at various drift angles and ice thicknesses
    图  14  船舶碎冰侧向力平均值变化曲线
    Figure  14.  Variation curves of the average lateral force of the ship against brash ice

    图13图14可知,随着漂角的增大,船舶两侧的碎冰侧向力峰值随机波动并增加;船舶迎冰侧的侧向力整体大于非迎冰侧。当α<10°时,在2种冰厚条件下碎冰迎冰侧的侧向力随漂角变化的趋势相似,均在α = 7.5°时出现下降;当α>10°时,0.014 97 m冰厚下船舶所受碎冰侧向力的增幅更大。随着冰厚的增加,α = 15°时碎冰侧向力的增幅最大,达24.4%。因此,在大漂角工况且冰厚增大时,需关注迎冰侧碎冰侧向力突然增加时对航行的影响。

    本文采用离散元法建立了极地船舶碎冰区斜航数值模型,研究了斜航工况下船舶在碎冰区的冰载荷特性,主要你得到如下结论:

    1)对比直航工况,当漂角α增大时,船舶尾迹无冰航道宽度增加,迎冰侧碎冰堆积增多,与船体的接触作用增强,船舶所受碎冰阻力及碎冰侧向力均显著增加。在航速0.6 m/s、海冰密集度70%、冰厚0.014 97 m条件下,当α = 15°时船舶碎冰阻力及碎冰侧向力分别为直航工况电费4.25倍和6.04倍。

    2)在大漂角工况下,航速的增加使船舶受到的碎冰阻力和碎冰侧向力显著提高,当α = 5°时,碎冰阻力的增幅为28.5%,当α = 10°时,碎冰侧向力的增幅达81.3%。因此,当船舶在碎冰区以较高航速航行,在部分漂角条件下需谨慎操作,以防止突发碎冰阻力或碎冰侧向力增加带来的潜在风险。

    3)在大漂角工况下,冰厚的增加显著影响船舶所受碎冰阻力及碎冰侧向力,当α = 12.5°时,碎冰阻力的增幅为42.9%;当α = 15°时,碎冰侧向力增幅为24.4%。因此,在大漂角和厚冰环境下,应关注碎冰阻力和碎冰侧向力的变化,以确保航行的平稳与安全。

  • 图  1   弹簧−阻尼器接触模型

    Figure  1.   Spring-damper contact model

    图  2   计算域划分图

    Figure  2.   Division of the computational domain

    图  3   计算域及船体网格划分

    Figure  3.   Division of computational domain and hull mesh

    图  4   静水及碎冰阻力的试验验证

    Figure  4.   Experimental validation of hydrostatic and icebreaking resistance

    图  5   漂角示意图

    Figure  5.   Schematic diagram of the drift angle

    图  6   斜航状态下的碎冰分布图

    Figure  6.   Distribution of brash ice in oblique sailing

    图  7   不同速度和漂角下船舶所受碎冰阻力时历曲线

    Figure  7.   Time history curves of brash ice resistance for ships at various drift angles and speeds

    图  8   不同漂角下碎冰阻力平均值及标准差的变化曲线

    Figure  8.   Variation curves of the average values and standard deviations of brash ice resistance at different drift angles

    图  9   不同速度和漂角下船舶所受碎冰侧向力时历曲线

    Figure  9.   Time history curves of brash ice lateral force on ships at various drift angles and speeds

    图  10   船舶碎冰侧向力平均值变化曲线

    Figure  10.   Variation curves of the average value of lateral force of the ship against brash ice

    图  11   不同冰厚和漂角下船舶所受碎冰阻力时历曲线

    Figure  11.   Time history curves of brash ice resistance on ships at various drift angles and ice thicknesses

    图  12   不同漂角下碎冰阻力平均值及标准差变化曲线

    Figure  12.   Variation curves of the average values and standard deviations of brash ice resistance at different drift angles

    图  13   不同冰厚和漂角下船舶所受碎冰侧向力时历曲线

    Figure  13.   Time history curves of brash ice lateral force on ships at various drift angles and ice thicknesses

    图  14   船舶碎冰侧向力平均值变化曲线

    Figure  14.   Variation curves of the average lateral force of the ship against brash ice

    表  1   船舶模型主要参数

    Table  1   Main parameters of ship model

    参数 数值
    缩尺比 60
    总长L/m 2.04
    型宽B/m 0.37
    设计吃水D/m 0.13
    下载: 导出CSV

    表  2   碎冰颗粒主要参数

    Table  2   Main Parameters of the icebroken particle

    参数 数值
    缩尺比 60
    密度/(kg∙m−3) 917.0
    泊松比 0.3
    杨氏模量/Pa 1.67×107
    船−冰摩擦系数 0.138
    冰−冰摩擦系数 0.29
    下载: 导出CSV
  • [1]

    NG A K Y, ANDREWS J, BABB D, et al. Implications of climate change for shipping: opening the Arctic seas[J]. WIREs Climate Change, 2018, 9(2): e507. doi: 10.1002/wcc.507

    [2]

    CHAI W, LEIRA B J, HØYLAND K V, et al. Statistics of thickness and strength of first-year ice along the Northern Sea Route[J]. Journal of Marine Science and Technology, 2021, 26(2): 331–343. doi: 10.1007/s00773-020-00742-5

    [3]

    LUO W Z, JIANG D P, WU T C, et al. Numerical simulation of an ice-strengthened bulk carrier in brash ice channel[J]. Ocean Engineering, 2020, 196: 106830. doi: 10.1016/j.oceaneng.2019.106830

    [4] 何林, 柴威, 余象鹏, 等. 极地船舶冰载荷随机特性研究综述[J]. 船舶力学, 2023, 27(7): 1109–1117. doi: 10.3969/j.issn.1007-7294.2023.07.014

    HE L, CHAI W, YU X P, et al. Review on random nature of ice loads for Arctic ships[J]. Journal of Ship Mechanics, 2023, 27(7): 1109–1117 (in Chinese). doi: 10.3969/j.issn.1007-7294.2023.07.014

    [5]

    SAZONOV K, DOBRODEEV A. Ice resistance assessment for a large size vessel running in a narrow ice channel behind an icebreaker[J]. Journal of Marine Science and Application, 2021, 20(3): 446–455. doi: 10.1007/s11804-021-00226-x

    [6]

    SAZONOV K, DOBRODEEV A. Ice resistance assessment for a large size vessel running in a narrow ice channel behind an icebreaker[J]. Journal of Marine Science and Application, 2021, 20(3): 446-455. doi: 10.1007/s11804-021-00226-x.

    [7]

    JEONG S Y, JANG J, KANG K J, et al. Implementation of ship performance test in brash ice channel[J]. Ocean Engineering, 2017, 140: 57–65. doi: 10.1016/j.oceaneng.2017.05.008

    [8]

    JEONG S Y, KIM H S. Study of ship resistance characteristics in pack ice fields[C]//25th International Conference on Port and Ocean Engineering under Arctic Conditions. Delft, Netherlands: POAC, 2019: 9-13.

    [9]

    XUE Y Z, ZHONG K, NI B Y, et al. A combined experimental and numerical approach to predict ship resistance and power demand in broken ice[J]. Ocean Engineering, 2024, 292: 116476. doi: 10.1016/j.oceaneng.2023.116476

    [10]

    YANG B Y, ZHANG G Y, HUANG Z G, et al. Numerical simulation of the ice resistance in pack ice conditions[J]. International Journal of Computational Methods, 2020, 17(1): 1844005. doi: 10.1142/S021987621844005X

    [11]

    DUBROVIN O V, ALEKSANDROV M, MOOR R. Calculation of brash ice resistance based on model testing[M]. Translation Department of Naval Architecture and Marine Engineering, University of Michigan, 1970. (查阅网上资料, 未找到本条文献的英文信息, 请确认)

    [12] 徐雪松, 李鸿远, 吴嘉蒙, 等. 平板−碎冰−水耦合作用下的碎冰运动机理研究[J]. 中国造船, 2023, 64(2): 95–103. doi: 10.3969/j.issn.1000-4882.2023.02.009

    XU X S, LI H Y, WU J M, et al. Research on movement mechanism of crushed ice coupled with plate and water[J]. Shipbuilding of China, 2023, 64(2): 95–103 (in Chinese). doi: 10.3969/j.issn.1000-4882.2023.02.009

    [13] 季顺迎, 李紫麟, 李春花, 等. 碎冰区海冰与船舶结构相互作用的离散元分析[J]. 应用力学学报, 2013, 30(4): 520–526. doi: 10.11776/cjam.30.04.D032

    JI S Y, LI Z L, LI C H, et al. Analysis of interaction between ice floe and ship hull with discrete element method in broken-ice field[J]. Chinese Journal of Applied Mechanics, 2013, 30(4): 520–526 (in Chinese). doi: 10.11776/cjam.30.04.D032

    [14] 黄家煜, 周利, 刘仁伟, 等. 航道碎冰条件下极地船舶冰压力分布特性[J]. 中国舰船研究, 2024, 19(2): 53–61. doi: 10.19693/j.issn.1673-3185.03276

    HUANG J Y, ZHOU L, LIU R W, et al. Distribution characteristics of polar ship ice pressure in pack ice channel[J]. Chinese Journal of Ship Research, 2024, 19(2): 53–61 (in Chinese). doi: 10.19693/j.issn.1673-3185.03276

    [15]

    ZHONG K, NI B Y, LI Z Y, et al. Direct measurements and CFD simulations on ice-induced hull pressure of a ship in floe ice fields[J]. Ocean Engineering, 2023, 272: 113523. doi: 10.1016/j.oceaneng.2022.113523

    [16] 王世超, 刘刚. 考虑碎冰阻力和静水阻力的高效船型多目标优化方法[J]. 中国舰船研究, 2024, 19(6): 97–107. doi: 10.19693/j.issn.1673-3185.03998

    WANG S C, LIU G. Efficient ship hull multi-objective optimization method considering ice resistance and calm water resistance[J]. Chinese Journal of Ship Research, 2024, 19(6): 97–107 (in Chinese). doi: 10.19693/j.issn.1673-3185.03998

    [17]

    Siemens PLM Software. Simcenter STAR-CCM+® Documentation[M]. Plano: Siemens Digital Industries Software, 2020. (查阅网上资料, 未找到本条文献, 请确认)

    [18]

    JOHNSON K L. Contact mechanics[M]. Cambridge: Cambridge University Press, 1987.

    [19]

    DI RENZO A, DI MAIO F P. Comparison of contact-force models for the simulation of collisions in DEM-based granular flow codes[J]. Chemical Engineering Science, 2004, 59(3): 525–541. doi: 10.1016/j.ces.2003.09.037

图(14)  /  表(2)
计量
  • 文章访问数:  84
  • HTML全文浏览量:  31
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-11-12
  • 修回日期:  2024-12-06
  • 网络出版日期:  2024-12-18

目录

/

返回文章
返回