Dilated-polyhedron-based DEM analysis of the ice resistance on ship hulls in escort operations in level ice

The ice resistance on ships in escort operations in level ice are investigated using the discrete element method (DEM). A dilated polyhedron—generated by the Minkowski sum of a sphere and a polyhedron—is employed in the DEM; this dilated polyhedron-based DEM (DPDEM) is adopted to simulate the ship–ice interaction, wherein the contact force and bond-failure criterion are considered for the collision and fracture of sea ice, respectively. A three-point bending test was simulated with DPDEM, and a field test was conducted in the Bohai Sea to validate the DEM results. Further, a parametric analysis of flexural strength was conducted to identify the parameters involved in the bond-failure criterion. The ice resistance on icebreakers and cargo ships in level ice are simulated using DPDEM. The simulated ice resistances are compared with the Lindqvist and Riska formulas and the model test, which proves the validity of the DEM simulation. The interaction between ships and level ice is simulated parametrically to investigate the ice resistance on cargo ships with and without the icebreaker escort. Influencing factors such as ship speed, ice thickness, and ship breadth were examined to investigate the ice resistance on the escorted cargo ship. Analysis and change rules of the ice resistance on cargos affected by those factors were given.     

Research on short-term dynamic ice cases for dynamic analysis of ice-resistant jacket platform in the Bohai Gulf

Jacket platforms in the Bohai Gulf are often victims of random ice-induced vibrations. Their design originally did not account for dynamic ice loading on the jacket structures and significant ice-induced vibrations have been observed during winters. So it is urgent to find out as to how to consider the dynamic ice loading in the design stage of an ice-resistant jacket platform in the Bohai Gulf. In this paper, based on the state-of-art developments of dynamic ice load and ice–structure interaction, the authors proposed two different approaches to determine the short-term dynamic ice cases for dynamic analysis of ice-resistant jacket platform in the Bohai Gulf, named as the failure probability-based approach and the expected loss-based approach. Considering the variation of the ice environments and the variability of ice-resistant structure's properties, the random ice spectrum and the pseudo-excitation method (PEM) are employed to improve the efficiency of the procedure. Through finite element method (FEM) modeling of the real ice-resistant jacket platforms in the Bohai Gulf, such as MSW, MUQ and NW, the characteristics of the two proposed approaches have been investigated and compared. This paper provides an open and universal framework to assist designers and owners for determining the short-term dynamic ice cases in the Bohai Gulf, in which more information can be considered in the procedure to update and modify the results if available.     

On maximum forces exerted by floating ice on a structure due to constrained thermal expansion of ice

The problem of interaction between a floating ice cover and an engineering structure is considered, in which the ice–structure contact forces are caused by an increase in ice temperature due to solar radiation in situations, when the lateral thermal expansion of ice is constrained. The focus is on the determination of the maximum thermally-induced horizontal force exerted on a structure wall, assuming that the magnitude of this force is bound by the smallest force capable of fracturing the ice cover due to its buckling. The ice cover is modelled as a rectangular plate of uniform thickness, with its four edges being constrained by vertical rigid walls, and it is assumed that ice deforms, and eventually fails, by the mechanism of viscous creep buckling. The plate is subjected to in-plane axial compressive stresses developing in ice to prevent its thermal expansion due to solar heating, and is transversely (vertically) bent by the forces caused by the reaction of underlying water. The floating ice is treated as a material whose elastic and viscous properties depend on temperature and the ice porosity, and therefore they vary with time and the depth of ice. The results of numerical simulations, conducted for a variety of the ice plate horizontal dimensions, thicknesses and daytime temperature-change scenarios, illustrate the evolution of the plate deflection surface prior to its failure, and show the time variation of the maximum forces exerted by ice on a structure wall as functions of the ice thickness and maximum daytime temperature rise at the top surface of ice.     

Primary productivity and export fluxes on the Canadian shelf of the Beaufort Sea: A modelling study

We present a coupled sea ice–ocean-biological (including ice algae) model in the Arctic Ocean. The 1D model was developed and implemented on the Canadian Beaufort Sea shelf to examine the importance of different physical processes in controlling the timing and magnitude of primary production and biogenic particle export over an annual cycle (1987). Our results show that the snow and sea ice cover melt and/or break-up controls the timing of the phytoplankton bloom but primary producers (ice algae and phytoplankton) on the outer shelf are essentially nutrient limited. The total annual primary production (22.7 to 27.7 g-C m^? 2) is thus controlled by nutrient “pre-conditioning” in the previous fall and winter and by the depth of wind mixing that is controlled in part by the supply of fresh water at the end of spring (ice melt or runoff). The spring bloom represents about 40% of the total annual primary production and occurs in a period of the year when sampling is often lacking. Time interpolation of observed values to obtain total annual primary production, as done in many studies, was shown to lead to an underestimation of the actual production. Our simulated ratios of export to primary production vary between 0.42 and 0.44.     

Experiments on the resistance of a large transport vessel navigating in the Arctic region in pack ice conditions

In this study, we carried out model tests to investigate the ice failure process and the resistance experienced by a transport vessel navigating in the Arctic region in pack ice conditions. We tested different navigation velocities, ice plate sizes, and ice concentrations. During the tests, we closely observed several phenomena, including the modes of interaction of the ice ship and the moving and failure modes of ice. We also measured the vessel resistances under different conditions. The test results indicate that the navigation velocity is a significant determinant of the moving and failure modes of ice. Moreover, vessel resistance is remarkably dependent on the ice concentration and navigation velocity. The variances of the mean and maximum resistance are also compared and discussed in detail.     

基于SPH方法的连续式破冰数值模拟

根据海冰的力学性质,采用光滑粒子流体动力学(SPH)方法建立海冰本构模型,模拟冰锥实验并与实验结果对比,对冰材料模型准确性进行验证.利用显式动力分析软件LS-DYNA模拟破冰船连续式破冰,得到了海冰的损伤变形和破冰阻力,同时研究了船舶航速和海冰厚度等参数对冰载荷的影响.分析研究表明:应用SPH方法建立的海冰模型能够真实的反映海冰的力学性能,对连续式破冰过程实现较为准确的模拟,为船冰碰撞的数值模拟提供了新思路.     

冰区运输船极地海域航行阻力的模型试验研究（英文）

为了考查冰区运输船在北极结冰海域航行时所受到的航行阻力,进行了一系列室内物理模型试验。试验中船体模型通过一个单向测力传感器与主拖车上的刚性拖曳臂相连。试验包含三种主要的冰条件,即船体独立航行的平整冰条件、在破冰船引导下航行的航道碎冰条件和独立航行的浮冰条件。平整冰及航道碎冰条件下对多种船体航行速度进行考察,而浮冰条件下则针对不同的冰覆盖率展开研究。试验中对冰-船相互作用模式、冰排在船体不同位置处的破坏进程以及碎冰沿船体的运动状态进行了详尽的观测,并对船体的航行阻力进行了有效的测量。此外,文中还对航行阻力均值与极值之间的差异进行了详尽的讨论与分析。结果表明,冰条件与船体航行速度是影响船体航行阻力的重要因素。     

Measuring the thicknesses of the freshwater-layer plume and sea ice in the land-fast ice region of the Mackenzie Delta using helicopter-borne sensors

Helicopter-borne sensors have been used since the early 1990s to monitor ice properties in support of winter marine transportation along the east coast of Canada. The observations are used in ice chart production and to validate ice hazard identification algorithms using satellite advanced synthetic aperture radar (ASAR) imagery. In this study we evaluated the sensors' additional capability to monitor the freshwater plume characteristic beneath land-fast ice. During the Canadian Arctic Shelf Exchange Study (CASES) data were collected over the Mackenzie Delta in the southern Beaufort Sea where a buoyant river plume exists. Results showed that the electromagnetic–laser system could describe not only the ice properties but also the horizontal distribution of the freshwater plume depths that decreased in depth stepwise offshore as the flow of the buoyant plume was restricted by a series of ridge-rubble fields running parallel to the coast. Relative to the 2 m mean ice thickness, the plume layer depth varied from zero under mobile offshore pack ice to 3 m inshore of the third set of ridge-rubble fields.     

冰区航行船层冰作用下的结构响应

船冰碰撞是一个复杂的动力学过程,如何得到碰撞中的冰载荷一直是船舶碰撞研究领域的热点之一.本文分别建立300000 t冰区航行船和层冰的有限元模型,基于弹塑性理论及非线性有限元理论,利用MSC.Dytran对其碰撞进行数值仿真,模拟了船首及层冰的接触碰撞过程,最终得到船冰碰撞过程中的碰撞力、船首结构响应、船航速变化及能量耗散等参数.分析船冰碰撞过程中的碰撞机理及特性,并对冰区航行船特别是其船壳提出结构加强的建议,为设计冰区航行船提供一定的参考.     

中国新型极地考察船冰载荷预估研究

本文针对中国新型极地考察船所受冰载荷进行了研究,其中包含了平整冰阻力和浮冰碰撞力.对于船舶设计者和建造者来说,船舶冰载荷的预估非常重要.本文采用经验公式方法计算了该船在平整冰中航行的冰阻力,并且与模型试验结果进行了对比.结果显示,计算和试验结果中的冰阻力都随船舶航行速度的增大而增大,经验公式方法可以预测出合理的平整冰阻力.通过计算得到了该船的性能曲线,即该船在不同厚度平整冰中航行所能达到的速度.此外,本文还考虑了该船与圆形浮冰之间的三维斜向碰撞问题,采用解析方法评估了浮冰对船舶的撞击力,研究了撞击位置、法向框架角度以及浮冰尺寸对碰撞力的影响.基于计算结果,本文就冰载荷的预测进行了讨论并提出了一些建议.     

考虑水介质作用的船冰碰撞解耦方法及载荷预报

考虑水介质在船冰碰撞过程中水压变化的影响,是模拟水介质中船-冰碰撞的关键.本文开展数值模拟的动网格技术研究,基于RANS方程和改进的双相流流体体积函数VOF,实时预报船冰靠拢全过程中船体表面、冰体表面的压力变化情况.得到船体表面、冰体表面的P(V)函数.进而将水介质对船冰碰撞载荷的影响简化为对船体、冰体碰撞面预加载荷来实现,达到解耦目的.     

新型极地科考艇上冰过程的强度分析及可行性研究

新型极地科考船舶、装备的研发,对于极地的探索、开发起着至关重要的效果。主要查询水冰两栖艇的国内外发展现状,根据冰刀原理的启发,设计了这种新型水冰两栖艇及其爬冰方式。通过探讨新型水冰两栖艇的爬冰系统及其设计特点,运用材料失效、瞬态动力学等理论,分析验证水冰两栖艇爬冰的最佳形式。并使用ANSYS Workbench软件建模,验证其能满足爬冰的强度要求。     

Analysis of the contribution of ice algae to the ice-covered ecosystem in Lake Saroma by means of a coupled ice–ocean ecosystem model

A new coupled ice–ocean ecosystem model that links the pelagic and ice ecosystems was used to clarify the role of ice algae in ice-covered ocean ecosystems. The model was applied to Lake Saroma (Hokkaido, Japan) in 1992. Comparison of the model's results with observational data confirmed that the model reproduced the behavior of the ecosystem with acceptable accuracy during the period from winter to spring. The primary production of the ice algae is effectively transported into the pelagic system by means of physical releasing effects: brine convection, ice melting and freezing, and diffusion generated at the bottom of the ice. Ice algae released from the ice are rapidly exported because of their high sinking speed and the shallow depth of Lake Saroma. For this reason, the zooplankton in Lake Saroma cannot graze these released algae. However, zooplankton actively graze the ice algae living along the bottom of the ice. These results show that, before their release, ice algae play an important role as a food source for overwintering zooplankton. A sensitivity analysis revealed a positive correlation between the sinking speed of the released ice algae and the magnitude of the spring bloom by pelagic phytoplankton, and that the time when secondary production becomes active is an important factor in the linkage between these two algal populations.     

冰载荷下船舶螺旋桨强度的有限元分析

论文研究了冰载荷下船舶螺旋桨的强度问题。基于冰载荷经验公式和五种典型计算工况,利用有限元法进行了冰载荷下螺旋桨的强度分析。建立了一种冰/桨叶碰撞计算模型,对碰撞过程中接触力的变化及螺旋桨的动态响应进行了研究。计算结果可供螺旋桨结构设计时参考。     

Ports in Proximity – Competition and Coordination among Adjacent Seaports,edited by Theo Notteboom,Cesar Ducruet,and Peter de Langen

The summer minimum extent of Arctic sea ice shrank drastically in these years, and the opportunity on Arctic international shipping emerged. The Northern Sea Route (NSR), formerly blocked by permanent ice, was completely ice-free in September in the past 3 years. Because this route is much shorter than conventional Asia–Europe shipping lane, many maritime countries have paid attention to exploit the enormous potential of the Arctic Ocean from economical consideration. This study measured the economical advantage of the seasonal NSR by calculating the shipping cost saved.     

冰载荷下的船舶运动建模

随着极地航道的开辟以及极地海洋能源开发的需要,海冰与船体的相互作用以及极地船舶航行与作业安全越来越受到关注。本文主要针对破冰过程的几个典型阶段,对冰载荷对船舶的作用力进行了分析与计算,并考察了破冰船在不同的海冰密集度、冰厚以及船舶不同航速情况下的不同影响。     

Numerical estimation of ship resistance in broken ice and investigation on the effect of floe geometry

In this paper, a numerical model based on the non-smooth discrete element method is presented to investigate the effect of ice floe shape on ship resistance under low-concentration broken ice condition where ice fields consist of relatively small and unbreakable floes. The accuracy of the numerical model is validated by the comparison with a series of experiments conducted by the authors with artificial ice floes made from polypropylene (PP) material. The mean ice resistance estimated by the numerical model is in fairly good agreement with both the experimental results and those from the existing semi-empirical formulas. Additionally, the floe shape effect captured in the experiment, i.e. “rectangle-like” floes results in higher resistance than “elliptic-like” floes, is also predicted by the numerical model. Then, the validated numerical model is used to further investigate the effect of floe shape on ship resistance. It has been found that the rectangle floe results in higher ice resistance than other studied floe shapes. The ratio between the maximum and minimum caliper diameters of floe shows little influence on the ice resistance of ship. Finally, the effect of ice thickness on the ice resistance is also discussed.     

Distribution of diatoms in the Northeast Water Polynya, Greenland

This study formed part of the Northeast Water project (NEW project) which dealt with physical, geophysical and biological processes in the Northeast Water Polynya off Northeast Greenland. This was part of the International Arctic Polynya Programme (IAPP). The diatom composition of the water masses, sea ice and melt ponds was analysed to show the relationship between ice and the water column near the ice with regard to the origin and fate of the cells in the ice and melt ponds. Fragilariopsis oceanica, Fragiliria sp. I and Chaetoceros socialis usually dominated the phytoplankton, while the ice and melt pond samples showed a wide range of assemblages, with different single-celled pennates and two undescribed species, Navicula sp. 1 and Nitzschia sp. 1 often dominant. Planktonic algae in sea ice can be released into the water column during ice break-up and melt, thus contributing to the spring bloom in the water column, if the timing of the release and the species composition are correct. The number of different ice algal assemblages supports the theory that cells originated from the water column, the benthos and freshwater. In addition, differential growth in the sea ice or melt ponds often altered the relative abundance of species in comparison with what is usually found in their original habitat. However, many of the cells in the ice and melt ponds were dead (empty frustules), making it difficult to determine whether the cells had actually lived in these habitats.     

Numerical simulation of level ice impact on landing craft bow considering the transverse isotropy of Baltic Sea ice based on XFEM

Ice bending is a major failure mechanism of level ice when ships and marine structures interact with level ice. This paper aims to investigate the ice bending and ice load when level ice collides on ships and marine structures using numerical simulation method, and compare the numerical results with field test. The fracture of ice is simulated with extended finite element method (XFEM), and cohesive zone concept is used to describe the crack propagation. In order to consider the characteristics of S2 columnar ice, a transversely isotropic elastic material model is used for the ice bulk elements, and a transversely isotropic Tsai-Wu failure criterion is adopted to predict the initiation of cracks. A well-controlled field test of a landing craft bow colliding with level ice in Baltic Sea is simulated to verify the numerical scheme. The ice plate's continuous deformation, crack initiation and crack propagation at different impact velocities and angles are simulated and the results are discussed. In the simulation, the bending crack emerges at the midline of the top surface of ice plate, then propagates towards free boundary, and finally a circumferential crack forms. It is found that with the impact velocity increases, the bending load increases and the fracture size (perpendicular distance from the crack to the contact edge) decreases. And as the angle between the landing craft bow and vertical direction increases, the bending load and the fracture size decrease. The simulated results corresponds well with the field test. The competition between the circumferential crack and radial crack is also found in the simulation and will be discussed in this paper. The results show that this method well simulates the bending of level ice and predict the ice load, and provides a good approach for investigating the mechanism of different forms of level ice fracture.     

渤黄海冰情分析与对策思考

2009年冬季,我国渤海和黄海北部出现近30年同期最为严重的海冰灾害。海冰袭来,沿海大小港口均遭海冰袭击,部分港口被迫封港,海冰对港口设施、海上交通运输和生产作业造成了严重影响。文中分析了渤黄海海冰的特点和影响,并结合我国港口和海上交通运输特点,提出了应对措施和建议。