Establishmen and vertification of an improved similitude law for dynamic ice force model test

To study the similitude laws for dynamic ice force model test which is a hot problem in ice mechanics.In this paper,a new modified similitude law for dynamic ice force model test is proposed based on ice-induced virbration theory and two reasonable assumptions.And the new law is vertified with an ice-induced pier virbration example.The numerical calculation results show that,firstly,the error of frequency factor between test value and theory value is small as well as that of displacement and velocity;Secondly,the time scale completely meet the similitude relation,the waveforms anastomotic well,and the phenomenon of offset and omit do not occur;Thirdly,as the damping ratio between prototype and model is nearly equal,the improved similitude law will still be applied.Therefore,the new modified similitude law is feasible to study dynamic ice force model test and can provide reference for the ice mechanics.

Dynamics of Sea Ice and Spectrum of Ice Force

Spectrum and self-excite characters are the two significant characters of the dynamics of sea ice. The spectrum character of sea ice is mainly shown by the spectrum of ice force. The spectrum character of the sea ice is its intrinsic attributes. When the spectrum of ice force from the dynamic response of ice and structure interaction are evaluated, the effect of dynamic character of the structure must be eliminated. In this paper, the ice force spectrum at Bohai Bay and Liaodong Bay is evaluated from the displacement and strain responses of a single degree and a multi-degree freedom structure. The evaluated ice force spectrum can be used to define the spectrum character of ice in the analysis of ice induced vibration.

A viscoelastic-plastic constitutive model with Mohr-Coulomb yielding criterion for sea ice dynamics

A new viscoelastic-plastic （VEP） constitutive model for sea ice dynamics was developed based on continuum mechanics. This model consists of four components： Kelvin-Vogit viscoelastic model, Mohr-Coulomb yielding criterion, associated normality flow rule for plastic rehololgy, and hydrostatic pressure. The numerical simulations for ice motion in an idealized rectangular basin were made using smoothed particle hydrodynamics （SPH） method, and compared with the analytical solution as well as those based on the modified viscous plastic（VP） model and static ice jam theory. These simulations show that the new VEP model can simulate ice dynamics accurately. The new constitutive model was further applied to simulate ice dynamics of the Bohai Sea and compared with the traditional VP, and modified VP models. The results of the VEP model are compared better with the satellite remote images, and the simulated ice conditions in the JZ20-2 oil platform area were more reasonable.

A modified 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.

Mechanical and numerical models for sea ice dynamics on small-meso scale

On small-meso scale, the sea ice dynamic characteristics are quite different from that on large scale. To model the sea ice dynamics on small-meso scale, a new elastic-viscous-plastic （EVP） constitutive model and a hybrid Lagrangian- Eulerian （HLE） numerical method are developed based on continuum theory. While a modified discrete element model （DEM） is introduced to model the ice cover at discrete state. With the EVP constitutive model, the numerical simulation for ice ridging in an idealized rectangular basin is carried out and the results are comparable with the analytical solution of jam theory. Adopting the HLE numerical model, the sea ice dynamic process is simulated in a vortex wind field. The furthering application of DEM is discussed in details for modeling the discrete distribution of sea ice. With this study, the mechanical and numerical models for sea ice dynamics can be improved with high precision and computational efficiency.

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.

Numerical simulation of the dynamic characteristics of flow fields under ice

Atmospheric and oceanic drag are the main environmental forces controlling sea ice drift. Oceanic drag includes the form drag generated by water pressure gradients on the side of ice floes or on ice ridges, and the skin friction generated by viscous flow on the bottom of ice floes. In this study, we carried out a two-dimensional numerical simulation using FLUENT software to investigate the characteristics of dynamic flow under ice with a smooth undersurface. We studied water drag and flow field distribution below the ice under different conditions of ice draft and flow velocity, and the results agreed well with data from laboratory-based physical modeling tests, demonstrating the ability of the numerical model to reproduce the dynamic interactions between sea ice and the flow field. The degree of distortion in the flow field caused by ice increased as the ice draft increased. Vortexes occurred in the wake field of the floe, and the centers of the vortexes moved away from the ice with increasing ice draft. The simulated drag of water on ice showed a clear linear relationship with the square of the flow velocity.

Model Test Study on Ice-Induced Vibrations of Compliant Multi-Cone Structures

For the study on the ice-induced vibration of a compliant mono-cone structure,a series of model tests were performed from 2004 to 2006.In these tests,the ice sheet before the compliant conical structure was found to be failed in two-time breaking.Based on this important finding,model tests study of the ice force on a compliant multi-cone structure were performed from 2006 to 2007.In these tests,the ice sheet broke before each single cone non-simultaneously.The exciting energy of the total ice force was found to be in a wide range of frequencies,and the structure can be easily excited with nonlinear resonance.

Classification Criteria of Narrow/Wide Ice-ResistantConical Structures Based on Direct Measurements

Ice-induced structural vibration generally decreases with an increase in structural width at the waterline. Definitions of wide/narrow ice-resistant conical structures, according to ice-induced vibration, are directly related to structure width, sea ice parameters, and clearing modes of broken ice. This paper proposes three clearing modes for broken ice acting on conical structures: complete clearing, temporary ice pile up, and ice pile up. In this paper, sea ice clearing modes and the formation requirements of dynamic ice force are analyzed to explore criteria determining wide/narrow ice-resistant conical structures. According to the direct measurement data of typical prototype structures, quantitative criteria of the ratio of a cone width at waterline(D) to sea ice thickness(h) is proposed. If the ratio is less than 30(narrow conical structure), broken ice is completely cleared and a dynamic ice force is produced; however, if the ratio is larger than 50(wide conical structure), the front stacking of broken ice or dynamic ice force will not occur.

DRUCKER-PRAGER YIELD CRITERIA IN VISCOELASTIC-PLASTIC CONSTITUTIVE MODEL FOR THE STUDY OF SEA ICE DYNAMICS

Based on the characteristics of sea ice drifting and ridging at meso-small scale, the Drucker-Prager （D-P） yield criteria was introduced into the Viscoelastic-Plastic （VEP） constitutive model for the study of sea ice dynamics. In this model, the Kelvin-Vogit viscoelastic model was adopted in the elastic stage, and the associated normal flow rule was used in the plastic stage. Using the VEP model, the sea ice ridging process was simulated in an idealized rectangular basin, and the simulation results show that the simulated ice ridge thickness is consistent with the analytical solution. Moreover, the VEP model with the D-P yield criteria was also applied for the sea ice simulation of Bohai Sea, and the ice thickness, concentration, velocity, and ice stress were obtained in 48 h. The simulated thickness distributions agree well with the satellite images. The singular problem in the Mohr^7oulomb （M-C） yield criteria was overcome by the D-P yield criteria, and the computational efficiency was also improved. In the numerical simulations described above, the smoothed particle hydrodynamics was applied.