冰与斜面结构作用过程的离散元模拟与分析

doi:10.16183/j.cnki.jsjtu.2021.182

冰与斜面结构作用过程的离散元模拟与分析

王延林¹, 郭麒¹, 孙珊珊^,², 魏思浩¹, 许宁³

1.大连理工大学海洋科学与技术学院,辽宁盘锦 124221

2.大连理工大学物理学院,辽宁盘锦 124221

3.国家海洋环境监测中心,辽宁大连 116023

Discrete Element Simulation and Analysis of Ice-Inclined Structure Interaction

WANG Yanlin¹, GUO Qi¹, SUN Shanshan^,², WEI Sihao¹, XU Ning³

1. School of Ocean Science and Technology, Dalian University of Technology, Panjin 124221, Liaoning, China

2. School of Physics, Dalian University of Technology, Panjin 124221, Liaoning, China

3. National Marine Environment Monitoring Center, Dalian 116023, Liaoning, China

通讯作者: 孙珊珊,女,工程师,电话(Tel.):0427-2631888;E-mail:sunss@dlut.edu.cn.

责任编辑: 陈晓燕

收稿日期: 2021-05-26

基金资助:

国家自然科学基金(41806221)
国家自然科学基金(41676087)
中央高校基本科研业务费(DUT20JC27)
中央高校基本科研业务费(DUT21JC34)

Received: 2021-05-26

作者简介 About authors

王延林(1981-),男,辽宁省大连市人,副教授,从事海洋工程问题研究.

摘要

斜面结构是冰区海洋结构的重要形式之一.斜面倾角变化会改变海冰的主要破坏模式,使结构所受冰荷载的峰值有较大变化.为了更好地模拟平整冰破碎的随机性特点,构建了具有黏结-破碎功能的非规则分布的扩展圆盘单元海冰模型,并基于该模型模拟了平整冰与斜面结构作用的动力过程.通过将数值模拟所得冰力峰值与实测冰力峰值进行比较,对该模型进行了验证,分析了斜面倾角对冰荷载及海冰破坏模式的影响规律.结果发现:数值模拟的冰荷载随结构倾角的变化趋势与二维理论模型计算的变化趋势基本一致,斜面倾角的增加会导致弯曲破坏所占比例减小,并造成冰荷载峰值及发生概率都呈增大趋势,斜面倾角是影响海冰破坏模式变化以及冰荷载大小的重要因素.该研究工作可以为海冰的离散元数值模拟及斜面海洋工程结构的抗冰设计提供借鉴和参考.

关键词： 离散单元法; 黏结-破碎模型; 斜面结构; 冰荷载; 破坏模式

Abstract

Inclined structure is an important marine structure in the iced area. The change of the inclined angle will change the main failure mode of sea ice and affect the peak ice force acted on the structure. In order to simulate the random breaking characteristics of level ice, an irregular distributed dilated disk element model with bond-break function is constructed, and the dynamic process of interaction between the level ice and the inclined structure is simulated based on this model, which is verified by comparing the peak ice forces obtained by numerical simulation with the peak ice forces measured in the field. The influence law of the inclined angle on the ice force and the ice failure mode is analyzed. It shows that the variation of ice load with the changes of inclined angles simulated by the numerical method is basically consistent with the variation calculated by the two-dimensional theoretical model. With the increase of the incline angle, the proportion of bending failure decreases and the peak ice force and its occurrence probability increase. The inclined angle is an important factor in the change of sea ice failure modes and the peak ice force. This paper can be used as a reference for discrete element numerical simulation of sea ice and ice-resist design of inclined marine structures.

Keywords： discrete element method (DEM); bond-break model; inclined structure; ice force; failure mode

PDF (4028KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

王延林, 郭麒, 孙珊珊, 魏思浩, 许宁. 冰与斜面结构作用过程的离散元模拟与分析[J]. 上海交通大学学报, 2022, 56(9): 1168-1175 doi:10.16183/j.cnki.jsjtu.2021.182

WANG Yanlin, GUO Qi, SUN Shanshan, WEI Sihao, XU Ning. Discrete Element Simulation and Analysis of Ice-Inclined Structure Interaction[J]. Journal of Shanghai Jiaotong University, 2022, 56(9): 1168-1175 doi:10.16183/j.cnki.jsjtu.2021.182

根据冰块破碎形式的不同,可以将冰区海洋工程结构分为直立结构和斜面结构两种形式.其中,冰与直立结构作用主要发生挤压破坏,而冰与斜面结构作用主要发生弯曲破坏.相同条件下海冰发生弯曲破坏产生的冰荷载峰值要远小于海冰发生挤压破坏时产生的冰荷载峰值.例如,Timco等^[1]对波弗特海域沉箱结构的冰荷载进行了现场实测,通过对实测数据的分析发现,海冰发生挤压破坏时产生的冰荷载约为弯曲破坏时产生冰荷载的6倍.海冰发生弯曲破坏的概率又与斜面结构的倾角有关,如白旭等^[2]基于模型实验和现场实测数据,对船体斜面与海冰作用的过程进行研究,总结出破坏模式与倾角的关系:当斜面倾角小于70° 时,海冰以弯曲破坏为主;当斜面倾角在70°~82° 之间时,挤压破坏和弯曲破坏同时存在;当斜面倾角大于82° 时,海冰主要发生挤压破坏.因此,斜面结构倾角是冰区海洋工程结构设计时需要考虑的重要因素之一.

针对海冰问题,数值模拟方法是重要研究手段之一,其中离散单元法(DEM)广泛用于海冰问题的模拟.如Hopkins等^[3]采用圆盘单元模拟了浮冰在海浪作用下与海洋结构相互作用过程,通过将模拟结果与试验数据比较,验证了模型的合理性.李紫麟等^[4]基于圆盘单元模拟了船舶在碎冰区航行所受冰荷载的问题.上述研究都是基于DEM模拟破碎的浮冰与海洋工程结构的作用.为了同时模拟海冰与结构作用破碎的物理过程,有学者在离散单元法的基础上引入虚拟的“黏结带”概念,构造出具有黏结-破碎功能的离散单元模型,如球体离散单元^[5-6].但是这些工作都是基于均匀的离散单元,通过规则排布离散单元来构建海冰,这种模型冰破碎形式固定,不能很好地模拟海冰由于气泡、杂质及冻结不均匀等因素造成的海冰物理力学的随机性、非线性及各向异性的特性.因此,有必要构建一种更能体现海冰这种物理力学特性的海冰模型.

综上,本文通过非统一尺寸且初始分布随机的圆盘单元进行黏结来构建平整冰模型.这种平整冰模型相较于规则分布的离散单元模型冰,更能反映海冰破碎时的随机性、非线性及各向异性等特性.本文采用这种非规则排布的离散元海冰模型模拟平整冰与斜面结构作用并发生破碎的动力过程,分析斜面倾角对冰荷载的影响以及海冰破坏模式与斜面倾角间的关系.所构建的数值模型及相关研究内容可以为海冰的离散元数值模拟研究以及斜面海洋工程结构的抗冰设计提供参考.

1 离散元模型

1.1 扩展圆盘单元

参考自然界莲叶冰的形态,Hopkins^[7]构造了扩展圆盘单元模型,并用这种圆盘单元模拟了浮冰在海浪作用下与海洋结构相互作用过程,通过将模拟结果与试验数据比较,验证了模型的合理性.这种扩展圆盘单元是通过对一个平面圆盘和球体做闵可夫斯基和运算得到的,可以简单地描述为在二维圆盘表面上的每个点都扩充一个球体,这样就得到了一个表面光滑并具有一定厚度的三维圆盘模型^[7-8].这种扩展圆盘单元的直径等于平面圆盘直径和球体直径之和,单元的厚度是球体的直径.根据圆盘的几何性质可以把圆盘单元间的接触情况分为平面-曲面、曲面-曲面和平面-平面接触这3种类型.根据接触单元间的相对位置判断其接触类型,进而计算得出单元间的变形量,再根据弹簧-黏壶模型计算出单元间的接触力^[9].圆盘单元受到的外力除了单元间的接触力以外,还有浮力、拖曳力及附加质量等^[10].

1.2 非规则分布的圆盘单元及其黏结-破碎模型

海冰是由纯冰、卤水及空气组成的复杂物质^[11],具有非匀质的特点.数值模拟中为了更贴近海冰的非匀质特点,采用初始分布为非规则、非统一尺寸的圆盘单元进行黏结来构造平整冰,如图1所示.采用这种非统一尺寸扩展圆盘单元构建的平整冰可以真实模拟海冰非匀质的力学性质,使海冰破碎更加随机,更符合真实的海冰破碎过程.平整冰与斜面结构作用时,海冰破碎长度与冰厚之比约为7^[12],而本文算例中最薄的平整冰厚度是0.5 m,即海冰的破碎尺寸约为3.5 m.平整冰的破碎是黏结带的破碎,圆盘单元不破坏,因此圆盘单元的尺寸要小于破碎尺寸,这样圆盘单元尺寸对破碎尺寸的影响较小.本文数值模型中选择的圆盘直径范围为1.59~2.65 m.平整冰模型构建时,首先计算每对相邻圆盘单元间的最短距离,当单元间的距离小于冰厚的一半时进行黏结,否则不黏结.依据这样的判断准则,圆盘单元将与其直接相邻并距离较近的单元黏结.如图1中的1号单元,与其相邻的2~8号单元直接相邻且距离均小于0.5h(h为冰厚),因此将它们黏结在一起.基于这样的算法遍历所有圆盘单元,最终将计算域中的所有满足条件的圆盘单元黏结在一起,形成具有一定孔隙的平整冰模型.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 具有黏结-破碎功能的非规则分布圆盘单元冰模型

Fig.1 Irregularly arranged disk discrete element ice model with bond-break function

黏结模型中的黏结带是虚拟的,可以传递力和力矩,不考虑其质量,如图2所示.在黏结带的作用下原本离散的圆盘单元被黏结在一起形成一个整体.黏结带的横截面如图2(b)所示,图中,r为扩展球体的半径;r_d为二维圆盘的半径.参数r和r_d均取黏结带中尺寸较小的圆盘单元.本文假定应力最大值位于图2(b)中圆点标识的位置,即横截面长轴和短轴与单元边缘的交界点共计4个点,即 $A_{1}$ 、 $A'_{1}$ 、 $A_{2}$ 、 $A'_{2}$ . 其中, $A_{1}$ 、 $A'_{1}$ 是横截面边缘与长轴的两个交点,这两个点上的最大正应力和切应力分别为 $σ_{1}$ 和 $τ_{1}$ , $A_{2}$ 、 $A'_{2}$ 上的最大正应力和切应力分别为 $σ_{2}$ 和 $τ_{2}$ .图2(c)为黏结带受力简化图,图中t1、t2分别为黏结带横截面的长轴和短轴; $F_{b, n}$ 为黏结带横截面上的法向力; $F_{b, t 1} F_{b, t 2}$ 分别代表黏结带横截面上沿长轴、短轴方向力的分量; $M_{b, n}$ 代表绕法向轴n扭转的力矩; $M_{b, t 1}$ 、 $M_{b, t 2}$ 分别为绕长轴 $t_{1}$ 、短轴 $t_{2}$ 扭转的力矩.根据材料力学的相关理论,可以得到黏结带横截面的边缘处的4个点 $A_{1}$ 、 $A'_{1}$ 、 $A_{2}$ 、 $A'_{2}$ 上的正应力 $σ_{1}$ 、 $σ_{2}$ 和切应力 $τ_{1}$ 、 $τ_{2}$ 的计算公式如下:

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 圆盘间的黏结带模型

Fig.2 Model of band between two dilated disks

$\begin{matrix} \begin{array}{l} σ_{1} = \frac{|F_{b, n}|}{A} + \frac{|M_{b, n}|}{I_{t 1}} y_{m a x, t 1} \\ σ_{2} = \frac{|F_{b, n}|}{A} + \frac{|M_{b, t 2}|}{I_{t 2}} y_{m a x, t 2} \end{array}\} \end{matrix}$

(1)

$\begin{array}{l} τ_{1} = \frac{|F_{b, t 1}|}{A} + \frac{|M_{b, n}|}{I_{n}} y_{m a x, t 1} \\ τ_{2} = \frac{|F_{b, t 2}|}{A} + \frac{|M_{b, n}|}{I_{n}} y_{m a x, t 2} \end{array}\}$

(2)

式中: A为黏结带横截面面积; $y_{m a x, t 1} = r + r_{d}$ 为横截面中心O到 $A_{1}$ 或 $A'_{1}$ 的距离; $y_{m a x, t 2} = r_{d}$ 为横截面中心与交接点 $A_{2}$ 、 $A'_{2}$ 间的距离;I_t1、I_t2及I_n为黏结带横截面3个方向的惯性矩,其计算公式为

$I_{t 1} = \frac{π}{4} r^{4} + \frac{4}{3} r^{3} r_{d}$

(3)

$I_{t 2} = \frac{4}{3} r^{4} + 4 r^{3} r_{d} + 4 r^{2} r_{d}^{2} + \frac{4}{3} r r_{d}^{3}$

(4)

$I_{n} = (\frac{π}{4} + \frac{4}{3}) r^{4} + \frac{16}{3} r^{3} r_{d} + 4 r^{2} r_{d}^{2} + \frac{4}{3} r r_{d}^{3}$

(5)

真实海冰的破坏形式是复杂多样的,本文将海冰的破坏模式简单地分为3种,即挤压破坏、弯曲破坏及剪切破坏.根据正应力的正负可以判别力是拉应力还是压应力,将应力值分别与拉应力极限或者压应力极限对比判断黏结带是发生弯曲破坏还是挤压破坏,将计算得到的剪应力值与切应力极限对比判断黏结带是否发生了剪切破坏.

1.3 冰与斜面结构作用的离散元模型

计算域的宽度为30 m,宽度方向上采用循环边界来模拟大范围的斜面结构.平整冰的长度为 86 m,计算域长度200 m,水深15 m,斜面结构的倾角为α,如图3所示.圆盘单元的运动采用局部坐标系与整体坐标系,局部坐标系Oxyz固定在每个圆盘单元的质心;整体坐标系OXYZ固定于计算域内.局部坐标系和整体坐标系间的转换通过四元数方法实现^[4].为避免计算结果的偶然性,本文对冰排前沿部分海冰单元赋予垂直方向的初始微扰动.每个算例进行3次数值模拟,对应3种不同的微扰动初始速度 $v_{i} = (i - 1) \times 10^{- 12} \frac{m}{s}$ ,其中i=1,2,3,表示第i次数值模拟时施加在冰排前沿的初始速度微扰动.海冰的断裂强度参考了文献[13-14]的数值模拟工作,其余参数参考文献[1]在原型监测中的工况,相关参数值如表1所示.平整冰在海流的驱动下撞向斜面结构并发生破碎,随着模拟的进行,破碎的冰单元越来越多地在斜面结构前产生堆积.数值模拟中冰排与结构作用过程截图如图4所示.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 计算示意图(m)

Fig.3 Schematic diagram of calculation (m)

表1 数值模拟使用的参数

Tab.1 Parameters used in numerical simulation

参数	取值	参数	取值
水深H/m	15	海冰与结构间的摩擦因数μ_is	0.1
计算域长度L/m	200	冰单元间摩擦因数μ_ii	0.1
计算域宽度B/m	30	圆盘直径d/m	1.59~2.65
冰速U/(m·s^-1)	0.05	单元间挤压强度σ_b/MPa	1
重力加速度g/(m·s^-2)	9.81	单元间弯曲强度σ_f/MPa	0.5
冰弹性模量E/GPa	4	单元间剪切强度σ_s/MPa	0.6
冰泊松比ν	0.3	水密度ρ_w/(kg·m^-3)	1010
冰密度ρ_i/ (kg·m^-3)	900	斜面倾角α/ (°)	30~82
冰厚h/m	1.25

新窗口打开| 下载CSV

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 平整冰与斜面结构相互作用过程截图

Fig.4 Screenshot of interaction process between ice sheet and inclined structure

1.4 模型验证

根据现场实测数据^[1]的工况,本文对0.5、0.875 及1.25 m共3个冰厚分别进行3次模拟试验,并增加0.6、0.7、0.8、0.9、1.0及1.1 m共6种冰厚工况,并将这几个工况的冰荷载最大值与实测数据进行比较.图5所示为当冰厚为1.25 m时3次模拟的冰力时程.图5(a)~5(c)是冰力时程图,图中t为模拟时间,F为冰荷载, F^P为峰值冰荷载.冰力时程曲线有分明的波峰和波谷,可以看出作用过程中冰力的加、卸载过程.对比可以发现峰值冰荷载的大小以及出现的时间阶段在一定范围内具有随机性,但是3次试验结果得到的冰荷载峰值的最大值差异不大.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 冰厚为1.25 m时3次模拟计算结果

Fig.5 Simulation results at h=1.25 m

将数值模拟中得到的冰荷载值除以斜面结构的宽度,得到单位宽度的结构所受冰荷载情况,如图6所示.其中,黑色实线由实测数据中海冰发生挤压破坏的数据进行线性拟合得到;黑色虚线由实测数据中海冰发生弯曲破坏时的数据进行线性拟合得到;红色实线由实测数据中海冰发生混合破坏模式时的数据进行线性拟合得到.比较发现,数值模拟的冰力峰值基本在混合破坏拟合直线偏下的区域中.数值模拟的结果均比挤压破坏的冰力小、比弯曲破坏的冰力大,并在混合破坏冰力范围内.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 数值模拟数据与实测数据比较

Fig.6 Comparison of numerical simulation data and measured data

2 斜面倾角影响分析

2.1 斜面倾角对冰荷载的影响

基于该模型,本文模拟了海冰与不同倾角的斜面结构作用的动力过程.海冰与斜面结构作用的受力分析如图7所示.根据M-C理论有 $f_{s} = μ_{i s} f_{N}$ , $f_{s}$ 为碎冰上爬时与斜面结构产生的摩擦力, $f_{N}$ 为斜面对碎冰的支反力.将力分解到水平方向和竖直方向可知,当竖直方向分量为0时可以认为海冰仅发生挤压破坏^[2],即 $μ_{i s} s i n α - c o s α = 0$ .极区海冰与海洋工程结构间的摩擦因数范围是0.04~0.14^[15],当等式成立时倾角的范围是 82° ~88°.本文数值模拟中采用的海冰与结构间的摩擦因数是0.1,故当等式成立时对应的倾角为84.3°.由于仅发生挤压破坏时冰荷载增长迅速,所以为保守起见,本文取斜面最大倾角为 82°.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 海冰与斜面结构相互作用简图

Fig.7 Diagram of sea ice acting on an inclined structure

综上,设定斜面倾角工况分别为30°、40°、50°、60°、70°、75°、80° 和 82°.设定冰速为0.05 m/s,冰厚为1.25 m,冰排长度为86 m,模拟时长为 1285 s,部分冰荷载时程如图8所示.比较发现,冰荷载大小(整体变化趋势)随着斜面倾角的增加而增加.根据图8中冰荷载整体的变化规律,可将结构所受冰荷载分为两个阶段,即冰力发展阶段和冰力稳定阶段.在冰力发展阶段,冰力均值逐步增加;在冰力稳定阶段,冰力均值变化较小.当斜面倾角较小时,冰荷载有较长的发展阶段,如图8(a)所示;而当斜面倾角较大时,冰荷载很快进入稳定阶段,如图8(c)所示.当斜面倾角较小时,海冰会在斜面上爬升并发生堆积,因此冰力发展阶段主要是由海冰在斜面上爬坡堆积的重力作用造成的,当堆积稳定后,冰力就进入稳定阶段.而当斜面倾角较大时,海冰不容易发生爬坡堆积,故冰力会较快进入稳定阶段.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 不同斜面倾角下斜面结构的冰荷载数值模拟结果

Fig.8 Numerical simulation results of ice force at different inclined angles

平整冰与斜面结构作用后,海冰发生破碎后会在斜面结构前发生堆积,图9所示为数值模拟的堆积情况截图.当斜面倾角较小时,冰块更容易沿斜面上爬,上爬高度与斜面角度呈反比.因此当斜面角度较小时,形成稳定堆积高度的时间长.而当斜面倾角较大时,冰块不易沿斜面上爬,其爬坡高度有限,因此海冰在斜面结构前的堆积很快达到平衡.当斜面倾角较小时,斜面结构前碎冰的堆积爬坡现象明显,因此由爬坡堆积的冰块重力所造成的荷载对结构的影响也较大.海冰对斜面结构的作用力可以分为两部分,一是海冰爬坡堆积所产生的重力,二是海冰与斜面结构作用发生断裂破碎而产生的动冰力^[16].

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 不同斜面倾角工况800 s时碎冰堆积情况

Fig.9 Broken ice accumulation in 800 s at different inclined angles

冰荷载理论研究中,根据二维理论模型,在作用初期,即没有大量碎冰堆积形成时,其冰力的计算公式为^[17]

$F = \frac{F_{h}}{D} = σ_{b} {(\frac{ρ_{w} g h^{5}}{E})}^{0.25} 0.68 \frac{s i n α + μ_{i s} c o s α}{c o s α - μ_{i s} s i n α} + z_{m a x} h ρ_{i} g (s i n α + μ_{i s} c o s α) \times (\frac{s i n α + μ_{i s} c o s α}{c o s α - μ_{i s} s i n α} + \frac{c o s α}{s i n α})$

(6)

$z_{m a x} = 0.68 σ_{f} {(\frac{h}{E})}^{0.25} {(\frac{1}{ρ_{i} g})}^{0.75} \frac{1}{s i n α + μ_{i s} c o s α}$

(7)

式中:F_h为水平冰力;D为平整冰与结构的接触长度;z_max为海冰沿斜面上爬的最大高度^[14].

为了验证数值模拟结果的合理性,将本文数值模拟结果与上述理论值计算值进行比较,如图10所示.考虑到理论公式的实用范围是在海冰作用初期,图中增加了数值模拟初期时的冰力峰值采点.对比发现,离散元模型计算得到的初期冰荷载变化趋势与由理论公式得到的变化趋势基本一致,特别是当斜面倾角为82° 时,数值模拟结果与理论值几乎重合.数值模拟结果中,全程冰力峰值较大的原因是数值模拟计算的冰荷载包含了碎冰堆积与爬坡时的重力分力,而理论计算没有考虑碎冰堆积的重力影响,故数值模拟计算的结果比理论公式的计算值大.如图10所示,当斜面倾角小于70° 时,海冰的破坏模式以弯曲破坏为主;当斜面倾角不小于70° 时,海冰的挤压破坏比率超过弯曲破坏的比率.

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 冰荷载随斜面倾角变化趋势

Fig.10 Ice forces versus inclined angles

2.2 海冰破坏模式分析

在冰与结构作用过程中,海冰的主要破坏模式有挤压破坏、剪切破坏及弯曲破坏^[17-18].为了进一步探究斜面倾角对海冰破坏模式的影响,并验证前文关于冰荷载变化的分析,本文统计出不同斜面倾角工况下不同破坏模式发生的总次数,并计算得到不同破坏模式所占比率(C),如图11所示.可见,当斜面倾角为30° 时,海冰发生弯曲破坏的比率为85%.随着斜面倾角的增加,挤压破坏模式的发生概率逐渐增大,弯曲破坏模式所占比率逐渐下降.当斜面倾角为70° 时,海冰发生挤压破坏的概率超过发生弯曲破坏的概率;而当斜面倾角为80° 时,海冰主要发生挤压破坏.结合图10对比发现,挤压破坏使冰力峰值急剧增加.当斜面倾角小于70° 时,海冰主要发生弯曲破坏且冰力峰值增加相对缓慢,因此冰区斜面结构的倾角设计应尽量小于 70°,同时需综合考虑建造成本以及冰堆积等因素.

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 各破坏模式占比随斜面倾角变化趋势

Fig.11 Change of proportion of each failure mode with dip angle of inclined structure

3 结论

为了更好地表现海冰物理力学性能的随机性、非线性及各向异性的特点,本文构建了非规则分布且具有黏结-破碎功能的扩展圆盘单元海冰模型,并基于该模型模拟了海冰与不同斜面倾角结构作用并发生破碎的动力过程,获得以下结论:

(1) 基于本文提出的离散元模型模拟的冰力时程曲线有分明的波峰和波谷,可以看出作用过程中冰力的加卸载过程,且冰荷载的峰值大小及出现的时间具有一定随机性.表明该模型可以准确模拟出冰力的加卸载过程以及海冰力学性能的随机性.

(2) 斜面倾角对海冰破坏模式及冰荷载峰值影响较大,通过计算粗略认为,当斜面倾角超过70° 后,海冰发生挤压破坏的概率将超过弯曲破坏,结构所受冰荷载大小也将急剧增大.因此斜面结构的倾角设计应尽量小于 70°.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

TIMCO

G W

, JOHNSTON

Ice loads on the caisson structures in the Canadian Beaufort Sea

[J]. Cold Regions Science and Technology, 2004, 38(2/3): 185-209.

DOI:10.1016/j.coldregions.2003.10.007 URL [本文引用: 2]

[2]

白旭, 周利, 陶冶.

破冰结构角度对整冰失效模式的影响分析

[J]. 极地研究, 2018, 30(4): 406-410.

[本文引用: 2]

BAI

, ZHOU

, TAO

Slope angle of icebreaking structure and failure modes in level ice

[J]. Chinese Journal of Polar Research, 2018, 30(4): 406-410.

DOI:10.13679/j.jdyj.20180008 [本文引用: 2]

<p>Understanding of ice failure modes is very important for the estimation of the forces needed in icebreaking. This paper analyzed the mechanics of crushing during ice-structure interactions. The ice breaking force component can affect the ice failure mode, and its variation with the dimension of the icebreaking structure was investigated to further confirm the relationships between friction, slope angle and ice failure mode. On the basis of a range of friction coefficients of sea ice, it is concluded that bending failure occurs when the slope angle of the structure is below 70o. Crushing and bending failure occurs simultaneously when the slope angle lies between 70° and 82°. When the slope angle is above 82°, only crushing failure occurs. A model to derive the ice breaking force under different failure modes was presented.</p>

[3]

HOPKINS

M A

, SHEN

H H

Simulation of pancake-ice dynamics in a wave field

[J]. Annals of Glaciology, 2001, 33: 355-360.

DOI:10.3189/172756401781818527 URL [本文引用: 1]

[4]

李紫麟, 刘煜, 孙珊珊, 等.

船舶在碎冰区航行的离散元模型及冰载荷分析

[J]. 力学学报, 2013, 45(6): 868-877.

DOI:10.6052/0459-1879-13-020 [本文引用: 2]

采用离散元模型对碎冰区浮冰与船舶结构的相互作用进行了数值研究。碎冰由三维圆盘单元构成,并考虑其在海流作用下的浮力、拖曳力和附加质量。船体结构由一系列三角形单元组合构造。通过海冰与船体单元间的接触判断和接触力计算,确定海冰与船体结构之间的相互作用。采用以上离散单元模型对不同冰况（冰速、冰厚、冰块尺寸和密集度）以及航速条件下,海冰对船体的动力作用过程进行了数值分析,对比分析了以上因素对船体冰载荷的影响,可为冰区船舶的安全运行和结构设计提供一定的借鉴作用。

Zilin

, LIU

, SUN

Shanshan

, et al.

Analysis of ship maneuvering performances and ice loads on ship hull with discrete element model in broken-ice fields

[J]. Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(6): 868-877.

DOI:10.6052/0459-1879-13-020 [本文引用: 2]

With the development of navigation channels in Arctic regions and oil/gas exploitations in cold fields, the investigation of ice load on ship hull is required urgently in offshore engineering under various ice and maneuvering conditions. In this paper, the discrete element method (DEM) is adopted to simulate the interactions between drifting ice floes and a moving ship. The pancake ice floes are modeled with three-dimensional discrete elements considering the buoyancy, drag force and added mass of current. The ship hull is constructed with triangle elements. The interaction between ice floe and ship hull is determined through element contact detection and impact force calculation. The influences of ice conditions (current velocities and directions, ice thicknesses, concentrations and ice floe sizes) and navigational speed on the dynamic ice forces of ship hull are examined. This work is of reference value to the ship structure design and the navigation security in ice-covered fields reseaches.

[5]

狄少丞

. 基于GPU并行算法的海洋平台及船舶结构冰荷载的离散元分析[D]. 大连: 大连理工大学, 2015.

[本文引用: 1]

Shaocheng

. Discrete element simulation of ice load on offshore platform and ship hull based on GPU parallel algorithm[D]. Dalian: Dalian University of Technology, 2015.

[本文引用: 1]

[6]

龙雪, 宋础, 季顺迎, 等.

锥角对锥体结构抗冰性能影响的离散元分析

[J]. 海洋工程, 2018, 36(6): 92-100.

[本文引用: 1]

LONG

Xue

, SONG

Chu

, JI

Shunying

, et al.

Influence of cone angle on anti-icing performance of conical structure with numerical simulations of discrete element method

[J]. The Ocean Engineering, 2018, 36(6): 92-100.

DOI:10.1016/j.oceaneng.2008.10.006 URL [本文引用: 1]

[7]

HOPKINS

M A

Discrete element modeling with dilated particles

[J]. Engineering Computations, 2004, 21(2/3/4): 422-430.

DOI:10.1108/02644400410519866 URL [本文引用: 2]

[8]

SUN

S S

, SHEN

H H

Simulation of pancake ice load on a circular cylinder in a wave and current field

[J]. Cold Regions Science and Technology, 2012, 78: 31-39.

DOI:10.1016/j.coldregions.2012.02.003 URL [本文引用: 1]

[9]

王延林, 张亚伟, 许宁, 等.

核电取水口海冰堆积问题的数值模拟

[J]. 哈尔滨工程大学学报, 2019, 40(10): 1716-1722.

[本文引用: 1]

WANG

Yanlin

, ZHANG

Yawei

, XU

Ning

, et al.

Numerical simulation of sea ice pile-up on nuclear water intake structures

[J]. Journal of Harbin Engineering University, 2019, 40(10): 1716-1722.

[本文引用: 1]

[10]

SUN

S S

, SHEN

H H

Is the wave-induced impact load from pancake ice important for offshore structures

[C]// The 21st IAHR International Symposium on Ice. Dalian, China: IAHR, 2012.

[本文引用: 1]

[11]

丁德文

. 工程海冰学概论[M]. 北京: 海洋出版社, 1999.

[本文引用: 1]

DING

Dewen

. The generality of engineering glaciology of seas[M]. Beijing: Ocean Press, 1999.

[本文引用: 1]

[12]

屈衍, 岳前进, BRISTER

冰与锥体作用破碎周期及破碎长度分析

[J]. 冰川冻土, 2003, 25(Sup.2): 326-329.

[本文引用: 1]

Yan

, YUE

Qianjin

, BRISTER

Analysis of ice sheet break period and length on the ice-breaking cone

[J]. Journal of Glaciology and Geocryology, 2003, 25(Sup.2): 326-329.

[本文引用: 1]

[13]

RANTA

, POLOJÄRVI

, TUHKURI

Scatter and error estimates in ice loads-Results from virtual experiments

[J]. Cold Regions Science and Technology, 2018, 148: 1-12.

DOI:10.1016/j.coldregions.2018.01.002 URL

[14]

PAAVILAINEN

, TUHKURI

Parameter effects on simulated ice rubbling forces on a wide sloping structure

[J]. Cold Regions Science and Technology, 2012, 81: 1-10.

DOI:10.1016/j.coldregions.2012.04.005 URL [本文引用: 1]

[15]

TIMCO

G W

, WEEKS

W F

A review of the engineering properties of sea ice

[J]. Cold Regions Science and Technology, 2010, 60(2): 107-129.

DOI:10.1016/j.coldregions.2009.10.003 URL [本文引用: 1]

[16]

Croasdale

K R

, Cammaert

A B

, Metge

A method for the calculation of sheet ice loads on sloping structures

[C]// IAHR 94 Proceedings of the 12th International Symposium on Ice. Trondheim, Norway: IAHR, 1994: 874-885.

[本文引用: 1]

[17]

LUBBAD

, LØSET

A numerical model for real-time simulation of ship-ice interaction

[J]. Cold Regions Science and Technology, 2011, 65(2): 111-127.

DOI:10.1016/j.coldregions.2010.09.004 URL [本文引用: 2]

[18]

ZHOU

, CHUANG

Z J

, BAI

Ice forces acting on towed ship in level ice with straight drift. Part II: Numerical simulation

[J]. International Journal of Naval Architecture and Ocean Engineering, 2018, 10(2): 119-128.

DOI:10.1016/j.ijnaoe.2017.06.001 URL [本文引用: 1]

Ice loads on the caisson structures in the Canadian Beaufort Sea

2004

... 根据冰块破碎形式的不同,可以将冰区海洋工程结构分为直立结构和斜面结构两种形式.其中,冰与直立结构作用主要发生挤压破坏,而冰与斜面结构作用主要发生弯曲破坏.相同条件下海冰发生弯曲破坏产生的冰荷载峰值要远小于海冰发生挤压破坏时产生的冰荷载峰值.例如,Timco等^[1]对波弗特海域沉箱结构的冰荷载进行了现场实测,通过对实测数据的分析发现,海冰发生挤压破坏时产生的冰荷载约为弯曲破坏时产生冰荷载的6倍.海冰发生弯曲破坏的概率又与斜面结构的倾角有关,如白旭等^[2]基于模型实验和现场实测数据,对船体斜面与海冰作用的过程进行研究,总结出破坏模式与倾角的关系:当斜面倾角小于70° 时,海冰以弯曲破坏为主;当斜面倾角在70°~82° 之间时,挤压破坏和弯曲破坏同时存在;当斜面倾角大于82° 时,海冰主要发生挤压破坏.因此,斜面结构倾角是冰区海洋工程结构设计时需要考虑的重要因素之一. ...

... 根据现场实测数据^[1]的工况,本文对0.5、0.875 及1.25 m共3个冰厚分别进行3次模拟试验,并增加0.6、0.7、0.8、0.9、1.0及1.1 m共6种冰厚工况,并将这几个工况的冰荷载最大值与实测数据进行比较.图5所示为当冰厚为1.25 m时3次模拟的冰力时程.图5(a)~5(c)是冰力时程图,图中t为模拟时间,F为冰荷载, F^P为峰值冰荷载.冰力时程曲线有分明的波峰和波谷,可以看出作用过程中冰力的加、卸载过程.对比可以发现峰值冰荷载的大小以及出现的时间阶段在一定范围内具有随机性,但是3次试验结果得到的冰荷载峰值的最大值差异不大. ...

破冰结构角度对整冰失效模式的影响分析

2018

... 基于该模型,本文模拟了海冰与不同倾角的斜面结构作用的动力过程.海冰与斜面结构作用的受力分析如图7所示.根据M-C理论有

$f_{s} = μ_{i s} f_{N}$

$f_{s}$

为碎冰上爬时与斜面结构产生的摩擦力,

$f_{N}$

为斜面对碎冰的支反力.将力分解到水平方向和竖直方向可知,当竖直方向分量为0时可以认为海冰仅发生挤压破坏^[2],即

$μ_{i s} s i n α - c o s α = 0$

.极区海冰与海洋工程结构间的摩擦因数范围是0.04~0.14^[15],当等式成立时倾角的范围是 82° ~88°.本文数值模拟中采用的海冰与结构间的摩擦因数是0.1,故当等式成立时对应的倾角为84.3°.由于仅发生挤压破坏时冰荷载增长迅速,所以为保守起见,本文取斜面最大倾角为 82°. ...

破冰结构角度对整冰失效模式的影响分析

2018

... 基于该模型,本文模拟了海冰与不同倾角的斜面结构作用的动力过程.海冰与斜面结构作用的受力分析如图7所示.根据M-C理论有

$f_{s} = μ_{i s} f_{N}$

$f_{s}$

为碎冰上爬时与斜面结构产生的摩擦力,

$f_{N}$

为斜面对碎冰的支反力.将力分解到水平方向和竖直方向可知,当竖直方向分量为0时可以认为海冰仅发生挤压破坏^[2],即

$μ_{i s} s i n α - c o s α = 0$

Simulation of pancake-ice dynamics in a wave field

2001

... 针对海冰问题,数值模拟方法是重要研究手段之一,其中离散单元法(DEM)广泛用于海冰问题的模拟.如Hopkins等^[3]采用圆盘单元模拟了浮冰在海浪作用下与海洋结构相互作用过程,通过将模拟结果与试验数据比较,验证了模型的合理性.李紫麟等^[4]基于圆盘单元模拟了船舶在碎冰区航行所受冰荷载的问题.上述研究都是基于DEM模拟破碎的浮冰与海洋工程结构的作用.为了同时模拟海冰与结构作用破碎的物理过程,有学者在离散单元法的基础上引入虚拟的“黏结带”概念,构造出具有黏结-破碎功能的离散单元模型,如球体离散单元^[5-6].但是这些工作都是基于均匀的离散单元,通过规则排布离散单元来构建海冰,这种模型冰破碎形式固定,不能很好地模拟海冰由于气泡、杂质及冻结不均匀等因素造成的海冰物理力学的随机性、非线性及各向异性的特性.因此,有必要构建一种更能体现海冰这种物理力学特性的海冰模型. ...

船舶在碎冰区航行的离散元模型及冰载荷分析

2013

... 计算域的宽度为30 m,宽度方向上采用循环边界来模拟大范围的斜面结构.平整冰的长度为 86 m,计算域长度200 m,水深15 m,斜面结构的倾角为α,如图3所示.圆盘单元的运动采用局部坐标系与整体坐标系,局部坐标系Oxyz固定在每个圆盘单元的质心;整体坐标系OXYZ固定于计算域内.局部坐标系和整体坐标系间的转换通过四元数方法实现^[4].为避免计算结果的偶然性,本文对冰排前沿部分海冰单元赋予垂直方向的初始微扰动.每个算例进行3次数值模拟,对应3种不同的微扰动初始速度

$v_{i} = (i - 1) \times 10^{- 12} \frac{m}{s}$

,其中i=1,2,3,表示第i次数值模拟时施加在冰排前沿的初始速度微扰动.海冰的断裂强度参考了文献[13-14]的数值模拟工作,其余参数参考文献[1]在原型监测中的工况,相关参数值如表1所示.平整冰在海流的驱动下撞向斜面结构并发生破碎,随着模拟的进行,破碎的冰单元越来越多地在斜面结构前产生堆积.数值模拟中冰排与结构作用过程截图如图4所示. ...

船舶在碎冰区航行的离散元模型及冰载荷分析

2013

$v_{i} = (i - 1) \times 10^{- 12} \frac{m}{s}$

2015

锥角对锥体结构抗冰性能影响的离散元分析

2018

锥角对锥体结构抗冰性能影响的离散元分析

2018

Discrete element modeling with dilated particles

2004

... 参考自然界莲叶冰的形态,Hopkins^[7]构造了扩展圆盘单元模型,并用这种圆盘单元模拟了浮冰在海浪作用下与海洋结构相互作用过程,通过将模拟结果与试验数据比较,验证了模型的合理性.这种扩展圆盘单元是通过对一个平面圆盘和球体做闵可夫斯基和运算得到的,可以简单地描述为在二维圆盘表面上的每个点都扩充一个球体,这样就得到了一个表面光滑并具有一定厚度的三维圆盘模型^[7-8].这种扩展圆盘单元的直径等于平面圆盘直径和球体直径之和,单元的厚度是球体的直径.根据圆盘的几何性质可以把圆盘单元间的接触情况分为平面-曲面、曲面-曲面和平面-平面接触这3种类型.根据接触单元间的相对位置判断其接触类型,进而计算得出单元间的变形量,再根据弹簧-黏壶模型计算出单元间的接触力^[9].圆盘单元受到的外力除了单元间的接触力以外,还有浮力、拖曳力及附加质量等^[10]. ...

... [7-8].这种扩展圆盘单元的直径等于平面圆盘直径和球体直径之和,单元的厚度是球体的直径.根据圆盘的几何性质可以把圆盘单元间的接触情况分为平面-曲面、曲面-曲面和平面-平面接触这3种类型.根据接触单元间的相对位置判断其接触类型,进而计算得出单元间的变形量,再根据弹簧-黏壶模型计算出单元间的接触力^[9].圆盘单元受到的外力除了单元间的接触力以外,还有浮力、拖曳力及附加质量等^[10]. ...

Simulation of pancake ice load on a circular cylinder in a wave and current field

2012

核电取水口海冰堆积问题的数值模拟

2019

核电取水口海冰堆积问题的数值模拟

2019

Is the wave-induced impact load from pancake ice important for offshore structures

2012

1999

... 海冰是由纯冰、卤水及空气组成的复杂物质^[11],具有非匀质的特点.数值模拟中为了更贴近海冰的非匀质特点,采用初始分布为非规则、非统一尺寸的圆盘单元进行黏结来构造平整冰,如图1所示.采用这种非统一尺寸扩展圆盘单元构建的平整冰可以真实模拟海冰非匀质的力学性质,使海冰破碎更加随机,更符合真实的海冰破碎过程.平整冰与斜面结构作用时,海冰破碎长度与冰厚之比约为7^[12],而本文算例中最薄的平整冰厚度是0.5 m,即海冰的破碎尺寸约为3.5 m.平整冰的破碎是黏结带的破碎,圆盘单元不破坏,因此圆盘单元的尺寸要小于破碎尺寸,这样圆盘单元尺寸对破碎尺寸的影响较小.本文数值模型中选择的圆盘直径范围为1.59~2.65 m.平整冰模型构建时,首先计算每对相邻圆盘单元间的最短距离,当单元间的距离小于冰厚的一半时进行黏结,否则不黏结.依据这样的判断准则,圆盘单元将与其直接相邻并距离较近的单元黏结.如图1中的1号单元,与其相邻的2~8号单元直接相邻且距离均小于0.5h(h为冰厚),因此将它们黏结在一起.基于这样的算法遍历所有圆盘单元,最终将计算域中的所有满足条件的圆盘单元黏结在一起,形成具有一定孔隙的平整冰模型. ...

1999

冰与锥体作用破碎周期及破碎长度分析

2003

冰与锥体作用破碎周期及破碎长度分析

2003

Scatter and error estimates in ice loads-Results from virtual experiments

2018

Parameter effects on simulated ice rubbling forces on a wide sloping structure

2012

... 式中:F_h为水平冰力;D为平整冰与结构的接触长度;z_max为海冰沿斜面上爬的最大高度^[14]. ...

A review of the engineering properties of sea ice

2010

... 基于该模型,本文模拟了海冰与不同倾角的斜面结构作用的动力过程.海冰与斜面结构作用的受力分析如图7所示.根据M-C理论有

$f_{s} = μ_{i s} f_{N}$

$f_{s}$

为碎冰上爬时与斜面结构产生的摩擦力,

$f_{N}$

为斜面对碎冰的支反力.将力分解到水平方向和竖直方向可知,当竖直方向分量为0时可以认为海冰仅发生挤压破坏^[2],即

$μ_{i s} s i n α - c o s α = 0$

A method for the calculation of sheet ice loads on sloping structures

1994

... 平整冰与斜面结构作用后,海冰发生破碎后会在斜面结构前发生堆积,图9所示为数值模拟的堆积情况截图.当斜面倾角较小时,冰块更容易沿斜面上爬,上爬高度与斜面角度呈反比.因此当斜面角度较小时,形成稳定堆积高度的时间长.而当斜面倾角较大时,冰块不易沿斜面上爬,其爬坡高度有限,因此海冰在斜面结构前的堆积很快达到平衡.当斜面倾角较小时,斜面结构前碎冰的堆积爬坡现象明显,因此由爬坡堆积的冰块重力所造成的荷载对结构的影响也较大.海冰对斜面结构的作用力可以分为两部分,一是海冰爬坡堆积所产生的重力,二是海冰与斜面结构作用发生断裂破碎而产生的动冰力^[16]. ...

A numerical model for real-time simulation of ship-ice interaction

2011

... 冰荷载理论研究中,根据二维理论模型,在作用初期,即没有大量碎冰堆积形成时,其冰力的计算公式为^[17] ...

... 在冰与结构作用过程中,海冰的主要破坏模式有挤压破坏、剪切破坏及弯曲破坏^[17-18].为了进一步探究斜面倾角对海冰破坏模式的影响,并验证前文关于冰荷载变化的分析,本文统计出不同斜面倾角工况下不同破坏模式发生的总次数,并计算得到不同破坏模式所占比率(C),如图11所示.可见,当斜面倾角为30° 时,海冰发生弯曲破坏的比率为85%.随着斜面倾角的增加,挤压破坏模式的发生概率逐渐增大,弯曲破坏模式所占比率逐渐下降.当斜面倾角为70° 时,海冰发生挤压破坏的概率超过发生弯曲破坏的概率;而当斜面倾角为80° 时,海冰主要发生挤压破坏.结合图10对比发现,挤压破坏使冰力峰值急剧增加.当斜面倾角小于70° 时,海冰主要发生弯曲破坏且冰力峰值增加相对缓慢,因此冰区斜面结构的倾角设计应尽量小于 70°,同时需综合考虑建造成本以及冰堆积等因素. ...

Ice forces acting on towed ship in level ice with straight drift. Part II: Numerical simulation

2018

〈

〉