冰激振动及冰荷载动力特性的试验研究

回顾了解决海冰引起结构振动分析方法的要求，通过室内冰力模型试验得出了冰激振动的重要参数－冰力频率，幅值的特性。研究表明，在低冰速时，冰力频率与冰速对冰厚之比呈线性关系；高冰速时，冰力步与冰速对冰厚之比呈现恒定性，文章还根据试验数据以及统计分析，推荐了冰力幅值计算公式和取值范围。     

正倒锥体结构上冰力谱分析实验研究

利用研制的新型非冻结可破碎合成模型冰材料，进行了冰排与正倒锥体结构物作用的物理模拟试验研究，试验中观察记录了冰排的破碎过程并测定了作用于锥体结构物上的冰力，通过快速傅立叶变换方法由所测得的冰力时程曲线得到冰力频谱，分析并讨论了冰排速度，本结构物的水线直径，锥角以及冰排弯曲破坏方向等对冰力频谱对冰力频谱因子的影响。     

基于黏聚单元法的抗冰海洋平台与层冰相互作用数值模拟

通过商业有限元软件LS-DYNA采用黏聚单元法模拟层冰的断裂问题,在IIHR冰池模型试验相同缩尺比参数下对KULLUK抗冰平台受到的冰力进行数值计算,探究冰厚与冰速等冰况参数对平台冰力的影响。计算结果表明,在极地作业平台主要遭遇北极海域冰况范围下,相同冰厚情况时冰速变化对倒锥形结构平台的横向冰力影响不大,而同冰速时冰力随着冰厚的增加接近线性增加。同时将计算所得结果与同条件下IIHR冰池模型试验数据、ARCTEC、ACL非冻结模型冰试验数据、实船数据和通过ISO19906规范中提供的针对大锥面结构受到层冰冰力的经验公式计算所得数据等进行了对比,验证了选用冰参数的可靠性以及该数值模型在模拟大尺度宽锥形平台遭受层冰作用问题中的适用性。     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

The Shielding Effect of Multi-Pile Structures on Ice Force

The shielding effect of the front pile-row on the ice force acting on the hack pile-row is studied by ice force model tests. In the tests, the front pile-row is designed to model jacket legs and the hack pile-row to model the water resisting pipe-phalanx within the jacket. The shielding factor for ice force corresponding to different conditions are given in this paper. The research indicates that there are many factors, including the longitudinal and lateral spacing between the front and hack pile-row, ice attacking angle and the ratio of pile diameter to ice thickness, that influence the shielding effect on ice force.     

极地吊舱推进船舶回转性能的离散元分析

吊舱推进装置在极地船舶中被广泛应用,其转舵模块可以带动螺旋桨摆动从而产生任意方向的推进力,使船舶操纵更为灵活。提出了一种吊舱推进船舶冰区操纵的离散元方法,对具有吊舱推进装置的冰区船舶破冰过程进行了数值模拟。以“雪龙2”号破冰船为研究对象,计算分析了船舶定速直航时的冰阻力,并通过与Lindqvist经验公式的对比验证了冰阻力计算的合理性。在此基础上进一步对船舶在不同冰厚、不同吊舱转向角下的回转破冰运动进行了离散元模拟,分析了回转半径与船长比值随冰厚的变化规律。计算结果表明：船舶的回转性能随冰厚的增大而降低,并随吊舱转向角的增大而显著提高。     

锥体导管架海洋平台冰激振动的DEM-FEM耦合分析及高性能算法

在寒区海域,冰荷载是影响海洋平台安全运行的主要环境荷载之一,由其引起的冰激振动给平台结构及其上部设备带来了严重危害。为分析不同冰况下平台的振动响应,本文建立了导管架海洋平台冰激振动的离散元（DEM）-有限元（FEM）耦合模型。采用具有黏接-破碎性能的等粒径球体离散单元对海冰的破碎特性进行模拟,通过由梁单元构建的海洋平台有限元模型获得结构的振动响应。在离散元与有限元的接触区域中实现了两个模型间计算参数的传递。为提高该耦合模型的计算效率和规模,发展了基于动力子结构方法的DEM-FEM耦合模型。为验证该耦合模型的有效性和可行性,将不同冰况下得到的冰荷载与ISO19906和JTS 144-1-2010标准进行了对比。结果表明,计算得到的冰荷载与标准相近,且冰厚与冰荷载呈二次非线性关系。同时,从冰速和冰厚两方面对比了渤海四桩腿JZ20-2 MUQ锥体导管架平台冰激振动加速度的数值结果和现场实测数据,发现冰速与振动加速度呈线性关系,冰厚与振动加速度呈二次非线性关系,并且振动加速度与冰速和冰厚平方的乘积呈线性关系。     

冰和直立结构相互作用的静冰力可靠性分析

我国渤海每年约有3个月的结冰期,冰与海洋结构的相互作用远较波浪为甚,因此有必要进行冰与海洋结构相互作用的可靠性分析.首先,采用热力学的方法根据每年实测的环境参数对海冰平整冰厚进行数值模拟计算,然后根据计算结果对不同重现期的极值冰厚采用极值统计法给出了极值分布,并与其他文献作比较.其次,以应变速率为海冰强度的主要影响参数,给出了海冰强度的概率分布.最后,以Schwarz公式为静冰力计算公式,结合得到的冰厚和冰强度的概率分布,采用映射变换法进行不同重现期静冰力的可靠度分析计算,并用蒙特卡罗重点抽样法进行模拟.     

海冰动力学的混合拉格朗日-欧拉数值方法

综合考虑欧拉坐标下有限差分法(FDM)在海冰动力学计算中的效率,以及拉格朗日坐标下光滑质点流体动力学方法(SPH)对海冰流变行为的精确模拟,本文发展了一种海冰动力学的混合拉格朗日-欧拉(HLE)数值方法。该方法首先在拉格朗日坐标下将海冰离散为若干个具有厚度、密集度的海冰质点,并由这些海冰质点通过Gauss函数对欧拉网格上的海冰参量进行积分插值;然后,在欧拉坐标下对海冰动量方程进行差分计算以确定各网格节点的海冰速度,并由此采用Gauss函数积分插值出拉格朗日坐标下各海冰质点的速度分布;最后,通过对海冰质点运动和分布的计算,确定出各海冰质点的位置、厚度和密集度等参量。采用该HLE方法对规则区域内的海冰堆积过程和涡动风场作用下的海冰动力演化趋势进行了数值试验;最后,采用该HLE方法对渤海海冰的动力过程进行了72h数值模拟,其计算结果与卫星遥感图像和现场观测资料吻合较好。以上计算结果均表明该HLE方法在海冰动力学数值模拟中具有较高的计算效率和模拟精度,可用于海冰动力过程的数值模拟。     

渤海模型冰运动特征的现场观测研究

Sea ice drift is mainly controlled by ocean currents, local wind, and internal ice stress. Information on sea ice motion, especially in situ synchronous observation of an ice velocity, a current velocity, and a wind speed, is of great significance to identify ice drift characteristics. A sea ice substitute, the so-called "modelled ice", which is made by polypropylene material with a density similar to Bohai Sea ice, is used to complete a free drift experiment in the open sea. The trajectories of isolated modelled ice, currents and wind in the Bohai Sea during non-frozen and frozen periods are obtained. The results show that the currents play a major role while the wind plays a minor role in the free drift of isolated modelled ice when the wind is mild in the Bohai Sea. The modelled ice drift is significantly affected by the ocean current and wind based on the ice–current–wind relationship established by a multiple linear regression. The modelled ice velocity calculated by the multiple linear regression is close to that of the in situ observation, the magnitude of the error between the calculated and observed ice velocities is less than12.05%, and the velocity direction error is less than 6.21°. Thus, the ice velocity can be estimated based on the observed current velocity and wind speed when the in situ observed ice velocity is missing. And the modelled ice of same thickness with a smaller density is more sensitive to the current velocity and the wind speed changes. In addition, the modelled ice drift characteristics are shown to be close to those of the real sea ice, which indicates that the modelled ice can be used as a good substitute of real ice for in situ observation of the free ice drift in the open sea, which helps solve time availability, safety and logistics problems related to in situ observation on real ice.     

海冰与直立结构相互作用的破坏模式研究

海冰在结构前的破坏模式以及产生的冰载荷与结构尺度和海冰厚度密切相关。采用海冰离散元方法（DEM）模拟平整冰与直立结构相互作用过程中的海冰破坏模式及整体冰载荷。该离散元方法的计算参数通过与Norstr?msgrund灯塔的现场实测数据对比进行了可靠性验证。在此基础上,对不同宽厚比（结构宽度与海冰厚度的比值）工况下平整冰与直立结构作用过程中的海冰破坏模式和整体冰压力进行了离散元分析。模拟结果表明：当宽厚比＜10时,海冰破坏模式主要为挤压破坏;当10≤宽厚比＜30时,海冰混合破坏发生;当宽厚比≥30时,海冰屈曲破坏发生。从海冰断裂长度与平均冰压力两个方面进一步说明了海冰破坏模式的转变过程。最后,构建了直立结构极值冰压力计算公式。研究成果可为寒区直立结构的整体设计提供参考。     

海上风电单桩式支撑结构冰激振动及参数敏感性分析

针对渤海湾某风电场的海上固定式风机支撑结构,采用适用于大直径单桩结构的PSI曲线模拟桩土相互作用,并采用SACS软件建立支撑结构的动力分析模型。首先对支撑结构进行模态分析;其次考虑海冰结构的随机振动作用模式,根据适用于渤海湾的随机冰力谱构造随机冰载荷时程曲线,基于半耦合的时域方法,采用SACS软件对支撑结构进行冰激振动分析,输出塔筒顶部加速度、单桩基底剪力及倾覆力矩等响应参数的时程曲线和响应功率谱;最后针对冰厚、冰速和海冰强度等海冰参数对支撑结构的冰激振动进行敏感性分析。研究结果表明,在随机振动模式下,冰载荷及结构动力响应对冰厚和海冰强度较为敏感,在进行冰激振动分析时应合理确定冰厚和海冰强度等参数。     

Amodified discrete element model for sea ice dynamics

Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.     

Reliability Analysis of Ice-Induced Fatigue and Damage in Offshore Engineering Structures

—In Bohai Gulf,offshore and other installations have collapsed by sea ice due to the fatigueand fracture of the main supporting components in the ice environments.In this paper presented are someresults on fatigue reliability of these structures in the Gulf by investigating the distributions of iceparameters such as its floating direction and speed,sheet thickness,compressive strength,ice forces on thestructures,and hot spot stress in the structure.The low temperature,ice breaking modes and componentfatigue failure modes are also taken into account in the analysis of the fatigue reliability of the offshorestructures experiencing both random ice loading and low temperatures.The results could be applied to thedesign and operation of offshore platforms in the Bohai Gulf.     

Discrete elementmodeling of ice loads on ship hulls in broken ice fields

Ice loads on a ship hull affect the safety of the hull structure and the ship maneuvering performance in ice-covered regions. A discrete element method （DEM） is used to simulate the interaction between drifting ice floes and a moving ship. The pancake ice floes are modelled with three-dimensional （3-D） dilated disk elements considering the buoyancy, drag force and additional mass induced by the current. The ship hull is modelled with 3D disks with overlaps. Ice loads on the ship hull are determined through the contact detection between ice floe element and ship hull element and the contact force calculation. The influences of different ice conditions （current velocities and directions, ice thicknesses, concentrations and ice floe sizes） and ship speeds are also examined on the dynamic ice force. The simulated results are compared qualitatively well with the existing field data and other numerical results. This work can be helpful in the shil3 structure design and the navigation securitv in ice-covered fields.     

辽东湾JZ20—2海域海冰参数的概率分布

基于辽东湾JZ2 0 2海域 1996 2 0 0 0年 4个冬季的海冰定点观测资料和海冰数值模拟结果 ,对该海域的平整冰厚、冰速、冰向和压缩强度等海冰参数进行了概率分析 ,确定了各自的分布参数 ,并对冰速和冰向进行了联合概率分析。结果表明 :冰厚和冰速分别服从对数正态分布和瑞利分布 ,海冰压缩强度服从正态分布。计算结果可用于JZ2 0 2海域海洋结构可靠性设计和疲劳累积损伤分析的海冰参数 ,也可作为其邻近海域的参考资料     

Model tests on performance of offshore wind turbine with suction caisson foundation in sand

Abstract

The performance of steel caisson during and after installation with different penetration velocities in medium dense sand is presented. The applied jacking forces, the amount of formed soil heave and bearing capacity were measured in the model tests. The influence of penetration velocities on jacking forces, soil heave and bearing capacity were also discussed in detail. The results indicated that the jacking forces for caisson in medium dense sands were significantly affected by the penetration velocity. The larger the penetration velocity, the more soil flowed into the caisson cavity during installation. This will lead to larger inner shaft resistance and in turn more jacking forces required for the same penetration depth. The height of soil heave during installation increases with penetration velocity. The m value calculated by the ratio of the volumes of the soil heave to that of the penetrated caisson wall can be used to evaluate the soil heave. The larger the applied velocity, the larger the m value and larger bearing capacity of caisson after installation. The relationship between the m value and penetration velocity can be used to control the soil heave for a steel caisson with a wall thickness to external diameter ratio of 4.2% in medium dense sand by jacking method.     

渤海导管架平台桩柱冰压力随机过程模型及其参数确定

首先分析了冰压力随机过程的特征及其产生机理,从而按冰的破坏形式将冰压力随机过程分为挤压破坏过程、屈曲破坏过程和弯曲破坏过程.从渤海海冰对平合桩柱作用大量实测冰压力曲线中选取37条典型的冰挤压破坏时程曲线,经谱分析确定了相应冰压力随机过程的谱密度及其谱参数;进一步经统计分析,建立了谱参数与冰厚的关系;在此基础上,建立了绕圆桩柱的冰压力随机场模型,并利用大量实测数据,经统计回归分析,确定了场函数的具体表达形式,由此用冰压力随机场可以容易地确定圆桩柱的总冰压力随机过程及其谱密度.     

Wave–ice interactions in the marginal ice zone. Part 2: Numerical implementation and sensitivity studies along 1D transects of the ocean surface

The theoretical foundation of a wave–ice interaction model is reported in Part 1 of this study. The model incorporates attenuation of ocean surface waves by sea ice floes and the concomitant breaking of the floes by waves that determines the structure of the marginal ice zone (MIZ). A numerical implementation of the method is presented here. Convergence of the numerical method is demonstrated, as temporal and spatial grids are refined. A semi-analytical method, which does not require time-stepping, is also developed to validate the numerical results, when dispersion is neglected. The wave energy lost during ice breakage is parameterized, as part of the numerical method. Sensitivity studies are conducted in relation to the energy loss and also dispersive effects, the choice of the attenuation model, the properties of the wave field, and sea ice properties such as concentration, thickness and breaking strain. Example simulations intended to represent conditions in the Fram Strait in 2007, which exploit reanalyzed wave and ice model data, are shown to conclude the results section. These are compared to estimates of MIZ widths based on a concentration criteria, and obtained from remotely-sensed passive microwave images.     

渤海海冰现场监测的数字图像技术及其应用

在渤海冰区油气开发中,海冰给平台结构、油气运输和施工作业带来很大影响.油气作业区海冰参数精确、连续、实时的现场监测对分析油气开发的可靠性、检验海冰数值模式、校正海冰卫星遥感数据具有重要意义.针对渤海油气作业区的海冰运动和分布特性,通过数字图像技术对海冰的厚度、运动速度和密集度三个参数的提取进行了算法开发和软件研制.在2...