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

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

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.     

Short-term statistics of ice loads on ship bow frames in floe ice fields: Full-scale measurements in the Antarctic ocean

Ship operation and ice loading in floe ice fields have received considerable interest during recent years. There have been several numerical simulators developed by different institutes which can simulate ship navigation through floe ice fields and estimate ship performance and local ice loads. However, public data obtained from full-scale measurement covering comprehensively ship performance and ice loads under various ice thicknesses, concentrations and floe sizes are rare. The 2018/19 Antarctic voyage of the Polar Supply and Research Vessel (PSRV) S.A. Agulhas II gathered considerable data of the ship in floe ice fields under various thicknesses, concentrations, and floe sizes. The aim of this paper is to carry out statistical analysis to seek suitable probability distributions which adequately fit the measured ice load and therefore suitable to be used as parent distributions for long-term estimation. For this aim, three categories of probability distributions, namely standard distributions, truncated distributions and mixture distributions are tested. It is found that truncated distributions can fit the load data better than standard distributions bounded at the threshold. In addition, mixture distributions are shown to have promising features, which fit the data well and are able to separate distribution components. Subsequentially, the well-performed distributions are used as parent distributions to make long-term load estimations. The estimation results demonstrate that long-term estimations are sensitive to the selection of parent distribution, which addresses the importance of finding correct distribution to model short-term ice loads. The data of ten selected cases will be published for the use of other researchers.     

Investigations on the ship-ice impact: Part 1. Experimental methodologies

In the structural design of the Polar Class ships, glancing impact with ice has been considered as the governing load scenario for dimensioning the bow structure. At present, non-linear transient dynamic analysis has been reckoned as a fundamental requirement in the assessment of the structural strength of ships with higher ice classes. Under such requirement, understanding of the dynamic characteristics of the ship-ice impact load appears significant importance. The present study aims to investigate the spatial and temporal variations of the impact loads by simulating the glancing impact events between a polar research vessel and giant ice floes in ice tank. The ship-ice impact loads were recorded by a grid-based tactile sensor attached on the bow area of the model ship. To achieve a reasonable simulation on the design scenario described by the IACS Unified Requirements for Polar Class Ships (UR I), a series of methodological calibration tests were preliminarily carried out to determine the key parameters that should be carefully monitored and controlled, accompanying with thorough discussions on the ship-ice impact process. This paper provides detailed information on the preliminary methodological calibrations and the tactile data processing techniques, including the identification of the ship-ice contact area, the depiction of the ice loading trail and the outline of the spatial distributions of local ice pressures. A companion paper provides detailed analyses and discussions regarding the spatial and temporal variation of the ice impact loads from the formal tests based on the proposed testing procedures obtained by the methodological calibrations.     

船舶与冰排碰撞结构响应研究

冰排的层厚、流速以及特性在很大程度上影响船舶与冰排碰撞的结构响应。本文分别建立船舶和冰排的有限元模型,采用LS-DYNA进行碰撞接触计算,模拟船艏与冰排发生碰撞。主要考虑相同条件下分别改变冰排厚度、运动状态以及物理性质等单个碰撞参数,对比研究不同碰撞工况下船舶的损伤变形、碰撞力等结构响应差异。得出上述因素对船—冰碰撞的影响规律,为提高船舶抗冰载荷设计提供参考。     

破冰船冲撞式破冰模拟及结构损伤分析

提出一种适用于船冰碰撞的海冰材料模型。基于该材料模型,通过非线性有限元模拟冰台撞击刚性墙,验证该模型的正确性,并成功将其运用于某PC3冰级破冰船冲撞式破冰过程数值模拟中。研究破冰船冲撞不同厚度层冰的动响应结果,总结船冰相互作用过程中碰撞力和能量的变化、船体结构损伤特性和破冰能力,研究成果可为破冰船结构设计提供参考。     

Dynamic ice loads for offshore wind support structure design

For offshore wind farms which are planned in sub-arctic regions like the Baltic Sea and Bohai Bay, support structure design has to account for load effects from dynamic ice-structure interaction. There is relatively high uncertainty related to dynamic ice loads as little to no load- and response data of offshore wind turbines exposed to drifting ice exists. In the present study the potential for the development of ice-induced vibrations for an offshore wind turbine on monopile foundation is experimentally investigated. The experiments aimed to reproduce at scale the interaction of an idling and operational 14 MW turbine with ice representative of 50-year return period Southern Baltic Sea conditions. A real-time hybrid test setup was used to allow the incorporation of the specific modal properties of an offshore wind turbine at the ice action point, as well as virtual wind loading. The experiments showed that all known regimes of ice-induced vibrations develop depending on the magnitude of the ice drift speed. At low speed this is intermittent crushing and at intermediate speeds is ‘frequency lock-in’ in the second global bending mode of the turbine. For high ice speeds continuous brittle crushing was found. A new finding is the development of an interaction regime with a strongly amplified non-harmonic first-mode response of the structure, combined with higher modes after moments of global ice failure. The regime develops between speeds where intermittent crushing and frequency lock-in in the second global bending mode develop. The development of this regime can be related to the specific modal properties of the wind turbine, for which the second and third global bending mode can be easily excited at the ice action point. Preliminary numerical simulations with a phenomenological ice model coupled to a full wind turbine model show that intermittent crushing and the new regime result in the largest bending moments for a large part of the support structure. Frequency lock-in and continuous brittle crushing result in significantly smaller bending moments throughout the structure.     

A tri-axial ice model for simulating ice-stiffened panel impact: Experiments and numerical modeling

The interaction of ice with a polar ship is complex, which may involve several ice failure mechanisms such as local crushing, tensile, cracking, bending, shearing, sliding and the ship structures could attain permanent deformations. A reasonable modeling of ice behavior is, therefore, critical to the analysis of ship-ice interactions. This paper reports an experimental and numerical study on the behavior of stiffened panels subject to the impact of a wedge-shaped ice block/indenter. The tested stiffened panel was mounted to a Falling Weight Impact Tester, and measurements were taken to help understand the dynamic responses of the structure and the ice, and also the permanent deformations. Finite element analyses using ABAQUS/EXPLICIT were also performed for the lab tests. The proposed ice model features a multi-surface yield function, empirical failure criteria for a few ice failure modes, and a representation of the remaining load carrying capacity of crushed ice. This ice model is implemented in a user-defined subroutine VUMAT and uses the cohesive element in a numerical solution. The predicted and measured impact force and structural deformation compared very well with the tests, indicating that the proposed ice model coded in VUMAT is reasonable. A series of parametric studies was then carried out to identify the key parameters of this new ice model that would have major effects on the prediction of ice impacts. The paper intends to provide test data that may be useful in understanding the complex ice-ship interaction and to introduce a numerical solution with a new ice model that improves the simulation of high-energy ice-ship impacts.     

Numerical study of ice-induced loads on ship hulls

A numerical model is introduced in this paper to investigate both global and local ice loads on ship hulls. This model is partly based on empirical data, by which the observed phenomena of continuous icebreaking can be well reproduced. In the simulation of a full-scale icebreaking trial, the interdependence between the ice load and the ship’s motion is considered, and the three degree-of-freedom rigid body equations of surge, sway and yaw are solved by numerical integration. The variations in the level ice thickness and in the strength properties of ice can also be taken into account. The simulated ice loads on ship hulls are discussed through two case studies, in which the ship’s performance, the statistics of ice-induced frame loads, and the spatial distribution of ice loads around the hull are analyzed and compared with field measurements. As far as we know the present paper is the first to integrate all the features above. It is hoped that further studies on this numerical model can supplement the field and laboratory measurements in establishing a design basis for the ice-going ships especially for ships navigating in the first-year ice.     

冰区航行船舶冰阻力研究方法综述

船舶在冰区航行时的冰阻力性能一直是国内外关注和研究的重点。冰阻力的研究主要集中在经验方法、数值模拟和实验研究三个方面。同时,由于冰区船舶在航行过程中频繁地与冰层或浮冰产生碰撞,海冰的材料结构和力学性质对冰阻力的研究有重要影响。文章从冰的物理力学特性出发,简要回顾几种重要的冰本构关系模型及其适用性;并从经验方法及试验与数值模拟相结合两个角度,回顾和讨论浮冰区和平整冰区中船舶的冰阻力性能研究进展;最后,基于研究现状提出尚需进一步解决的问题。文章旨在介绍冰区船舶冰阻力性能的研究进展,望能为后续冰阻力研究提供参考。     

水介质对船冰碰撞结构响应的影响

以极地运输船舶艏部作为研究对象,建立基于流固耦合算法的船-水-冰耦合技术对三维船艏与冰体碰撞的结构响应问题进行了研究,结合非线性有限元软件LS-DYNA对比分析了考虑流固耦合(即有水介质)的船-水-冰碰撞模型和不考虑流固耦合(即无水介质)的船-冰碰撞模型与同质量不同速度的冰体发生碰撞下的结构响应问题,揭示了不同碰撞工况下船舶在损伤变形、碰撞力、速度等方面的变化特征及差异,同时阐述了水介质在船冰碰撞中的作用,可以为极地运输船舶的抗冲击结构设计提供参考。     

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.     

冰区海上风电基础的抗冰性能分析

海上风电基础属于典型的柔性结构。由于冰与柔性抗冰结构相互作用的复杂性,长期以来尚未形成基于动冰力响应分析的结构设计。结构抗冰设计中大都是从极端荷载出发,只考虑最大静冰力或最大倾覆力矩。基于对渤海辽东湾柔性抗冰结构的多年监测,发现强烈的冰激振动引起柔性结构的风险性要远大于极端静冰荷载下结构的整体安全问题。为了明确冰区风电基础结构的抗冰性能及抗冰设计的合理性,文章结合基于多年现场冰与结构作用观测及冰荷载的研究成果,明确该类柔性结构与海冰作用形式及其动力特性;提出了柔性抗冰结构设计中应考虑的主要失效模式及评价方法。最后,以渤海某典型风电基础为例,对其抗冰性能进行评价。该文的研究可为寒区风电基础的抗冰设计及安全保障提供合理依据。     

基于SPH法的冰与船舶螺旋桨碰撞数值模拟

在结冰海域航行的船舶螺旋桨会因受到冰块的撞击而出现严重的损坏。该文建立了冰—桨碰撞计算模型,冰模型采用SPH(光滑粒子流体动力学)法,运用非线性有限元软件ANSYS/LS-DYNA数值模拟了螺旋桨与冰在不同速度、位置下碰撞和螺旋桨与不同半径冰碰撞下的动态响应。对比研究了在以上单一碰撞因素下桨的变形、单元最大应力以及碰撞过程中产生的碰撞力等响应差异,定性地得出上述因素对冰—桨碰撞的影响规律,为船舶螺旋桨抗冰性能的结构设计提供参考。     

Fatigue damage induced by vortex-induced vibrations in oscillatory flow

Vortex-induced vibration (VIV) of a flexible cylinder in oscillatory flow was experimentally investigated in an ocean basin. An intermittent VIV was confirmed to have occurred during the tests. The fatigue damage caused by VIV was calculated based on rainflow counting and a standard S–N curve. There are 3 main observations for fatigue damage from VIV in oscillatory flow: 1) the damage varied significantly with the KC number, which is a unique feature for VIV in oscillatory flow. 2) Fatigue damage at small KC number cases was found to be larger compared to those at large KC numbers owing to the fact that number of vortex shedding cycles per half of the motion cycle is low, and damping within half of the motion cycle will hence become low. The fact that vortices from the previous cycle still are active during the next may also contribute to the large response at small KC numbers. 3) ‘Amplitude modulation’ and ‘mode transition’, two specific features for VIV in oscillatory flow, were found to have a strong influence on fatigue. Fatigue damage has also been calculated by an empirical VIV prediction model assuming that all cases have steady flow at an equivalent velocity. Finally, a simplified method for calculating fatigue damage from VIV in oscillatory flow based on steady flow conditions is proposed. A modification factor diagram is presented, but the scope of the present study is too limited to provide a good basis for a general model for this factor. A general model for how to apply results from constant current analysis to predict fatigue in oscillatory flow will therefore need further research.     

Influence of different combinations of impact mass and velocity with identical kinetic energy or momentum on the impact response of RC piles

Marine engineering structures may be subjected to various levels of vessel collisions, causing different degrees of damage to the RC piles. However, the influence of different combinations of impact mass and velocity under the same initial kinetic energy or momentum on the impact behavior of RC piles has not been explored. A numerical investigation of the impact responses of RC pile structures (including plumb piles and inclined pile groups) considering the material strain rate effect is carried out using the software ABAQUS/Explicit in this paper. Parametric studies are conducted based on the validated finite element models of RC piles subjected to horizontal impact. The internal force distributions, as well as the effect of the impact kinetic energy and the momentum (including different combinations of impact mass and initial velocity) on the impact response of RC piles, are discussed. The results indicate that the development path of the bending moment‒overall structural displacement curves in the plastic hinge regions of RC piles are almost identical under both impact and static loading conditions. Compared to momentum, the impact response of RC piles is insensitive to the variation of impact kinetic energy with different combinations of impact mass and velocity. Therefore, the impact kinetic energy is more suitable for evaluating the structural deformation of RC piles after impact. Finally, a performance-based impact-resistance design framework for RC piles is proposed.     

A comparative analysis of the fluid-structure interaction method and the constant added mass method for ice-structure collisions

The constant added mass (CAM) method and the fluid-structure interaction (FSI) method are widely used to simulate ship-ship and ship-ice collisions. In the CAM method, the hydrodynamic effect of the surrounding water is treated as a constant added mass, whereas in the FSI method the surrounding fluid flow is explicitly modelled. The objective of the paper is to compare the two methods and to explain the causes of the differences in the results. We considered collision between a freshwater ice mass and a floating steel structure. For both methods, the numerical simulations were performed with the LS-DYNA software. The behaviour of the ice mass was modelled using an elliptic yield criterion and a strain-based pressure-dependent failure criterion. To ensure realistic ice behaviour, the ice model was calibrated using general trends found in laboratory and in-situ indentation tests with focus on the laboratory-grown ice and the fluid model in the LS-DYNA was verified by comparing the added mass coefficients for a spherical body and a rectangular block with the corresponding WADAM results. To validate and benchmark the numerical simulations, experimental data on ice-structure interactions in water were used, including the acceleration of the floater wall measured with the dynamic motion unit (DMU), the relative velocity between the ice mass and the floater before the impact and some images extracted from video recording of the test. The comparisons indicated that the FSI method yields better results for the motion of the floater, i.e., the acceleration of the floater wall caused by the ice mass’s impact and the relative velocity were in reasonably good agreement with experimental measurements. It was also found that the CAM method was faster but predicted a higher peak contact force and more dissipated energy in the ice mass than in the FSI method.     

海上三体风电运维船阻力和运动性能综合研究

结合海上风电运维船的使用要求和三体船型特点,确定小型高速三体风电运维船船型,并寻求特定海况下阻力性能和运动性能优异的侧体布局方案。基于CFD计算静水阻力和规则波中迎浪增阻与运动响应,并预报不规则波中运动有义值和阻力值。通过比较表明,主、侧体艉部对齐时三体运维船的综合性能最优。从而为海上三体风电运维船的设计和性能预报提供参考。     

Distributed mass/discrete floe model for pack ice rheology computation

A new model, called the distributed mass/discrete floe model, is proposed for performing practical computations of mesoscale pack ice rheology. This model possesses the advantages of both the continuum and the discrete element models: it can express the discrete nature of pack ice, for which it is difficult to use a continuum model, and requires a much shorter computation time than a discrete element model. The pack ice is divided into ice bunches in which the floes, assumed to be distributed uniformly, are modeled as inelastic disks or rectangles floating on the water. The ice interaction forces are formulated from the relationship between the impulse on the bunch and the variation of momentum in the bunch. The ocean flow is calculated simultaneously with the floe movement using a multilayer model. In a circulating water channel, drift tests of physical model floes were performed in order to investigate the characteristics of their motion and interaction with ocean structure models. Near the structure, the floe motion depends on the floe shape. Disk floes show a lateral motion in front of the structure. They flow out around both sides of the structure and the number of floes in front of the structure decreases with the lapse of time. On the other hand, rectangular floes scarcely flow laterally. The number of floes in front of the structure remains constant over time. These experiments indicate that when the motion of pack ice around a structure is simulated, it is important to consider the floe shape. The disk floe motion and the rectangular floe motion can be regarded as extreme cases of pack ice motion. Actual pack ice motion may be between these two extremes. Computations were carried out using the distributed mass/discrete floe (DMDF) model. Simulation results were compared with the circulating water channel experiment results and sea ice motion in the southern part of the Sea of Okhotsk. The DMDF model predicted the circulating water channel drift test results quite closely. The DMDF model results also compared quite well with the sea ice motion.     

海冰对港口作业的影响及应对措施

渤海和北黄海每年冬季都有结冰现象。海冰能够封锁港口、堵塞航道、破坏近海结构物、使航行中的船舶受损,给船舶进出港航行、操纵及靠离码头作业等带来极大困难,对港口运营也有较大的影响,严重情况会导致重大海难事故。从港口建设决策者及运营管理者的角度出发,多方面阐述了海冰对港口作业的影响,并重点论述了港口应对海冰灾害的措施,为冰况区港口能够延长作业天数及安全生产提供参考。