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 cr...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.展开更多
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 ...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.展开更多
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...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.展开更多
For ship structural design and good maneuverability in an ice-covered sea, the local and global load of ice cover on ships should be well understood. This paper reviews the extensive work done on ice loads on ships, i...For ship structural design and good maneuverability in an ice-covered sea, the local and global load of ice cover on ships should be well understood. This paper reviews the extensive work done on ice loads on ships, including： （a） Ice pressure and local load determination based on field and model tests; （b） Global ice loads on ships from full-scale field observations, model tests and numerical models under different ice conditions （level ice and pack ice） and ship operations （maneuvering and mooring）. Spe- cial attention is paid to the discrete element simulation of global ice loads on ships; and （c） Analytical solutions and numerical models of impact loads of icebergs on ships for polar navigation. Finally, research potential in these areas is discussed.展开更多
A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed partic...A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.展开更多
The splitting test is a competitive alternative method to study the tensile strength of sea ice owing to its suitability for sampling.However,the approach was questioned to the neglect of local plastic deformation dur...The splitting test is a competitive alternative method to study the tensile strength of sea ice owing to its suitability for sampling.However,the approach was questioned to the neglect of local plastic deformation during the tests.In this study,splitting tests were performed on sea ice,with 32 samples subjected to the regular procedure and 8 samples subjected to the digital image correlation method.The salinity,density,and temperature were measured to determine the total porosity.With the advantage of the digital image correlation method,the full-field deformation of the ice samples could be determined.In the loading direction,the samples mainly deformed at the ice-platen contact area.In the direction vertical to the loading,deformation appears along the central line where the splitting crack occurs.Based on the distribution of the sample deformation,a modified solution was derived to calculate the tensile strength with the maximum load.Based on the modified solution,the tensile strength was further calculated together with the splitting test results.The results show that the tensile strength has a negative correlation with the total porosity,which agrees with previous studies based on uniaxial tension tests.展开更多
The quasi-solid-fluid phase transition of granular materials is closely related to the shear rate and solid concentration in addition to their intrinsic properties. The contact duration and the coordination number are...The quasi-solid-fluid phase transition of granular materials is closely related to the shear rate and solid concentration in addition to their intrinsic properties. The contact duration and the coordination number are two important temporal-spatial parame-ters to describe the granular interaction in phase transition. In this study, characteristics of the contact duration and the coordination number associated with the transition processes are determined using a 3D discrete element model under different shear rates and concentrations. The resulting macroscopic stress and strain-rate relations are discussed. The temporal and spatial parameters provide a linkage between the macroscopic constitutive law and inter- particle micromechanics.展开更多
The quasi-solid-liquid phase transition exists widely in different fields,and attracts more attention due to its instinctive mechanism.The structure of force chains is an important factor to describe the phase transit...The quasi-solid-liquid phase transition exists widely in different fields,and attracts more attention due to its instinctive mechanism.The structure of force chains is an important factor to describe the phase transition properties.In this study,the discrete element model(DEM) is adopted to simulate a simple granular shear flow with period boundary condition on micro scale.The quasi-solid-liquid phase transition is obtained under various volume fractions and shear rates.Based on the DEM results,the probability distribution functions of the inter-particle contact force are obtained in different shear flow phases.The normal,tangential and total contact forces have the same distributions.The distribution can be fitted as the exponential function for the liquid-like phase,and as the Weibull function for the solid-like phase.To describe the progressive evolution of the force distribution in phase transition,we use the Weibull function and Corwin-Ngan function,respectively.Both of them can determine the probability distributions in different phases and the Weibull function shows more reasonable results.Finally,the force distributions are discussed to explain the characteristics of the force chain in the phase transition of granular shear flow.The distribution of the contact force is an indicator to determine the flow phase of granular materials.With the discussions on the statistical properties of the force chain,the phase transition of granular matter can be well understood.展开更多
基金The authors would like to acknowledge the supports by the National Natural Science Foundation of China under contract No.40206004partly by the East-Asia and Pacific Program of US National Science Foundation under contract No.INT-9912246.
文摘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.
基金The Special Funding for National Marine Commonwealth Industry of China under contract Nos 201105016 and 2012418007the National Natural Science Foundation of China under contract No.41176012
文摘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.
基金Special Fund of Marine Commonweal Industry under contact Nos 201105016 and 201205007supported by National Marine Environment Forecasting Centrethe National Natural Science Foundation of China under contact No.41176012
文摘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.
基金supported by the Ocean Public Welfare Scientific Research Project of China (Grant No. 201105016,2012418007)the National Natural Science Foundation of China (Grant No.41176012)the American Bureau of Shipping (ABS)
文摘For ship structural design and good maneuverability in an ice-covered sea, the local and global load of ice cover on ships should be well understood. This paper reviews the extensive work done on ice loads on ships, including： （a） Ice pressure and local load determination based on field and model tests; （b） Global ice loads on ships from full-scale field observations, model tests and numerical models under different ice conditions （level ice and pack ice） and ship operations （maneuvering and mooring）. Spe- cial attention is paid to the discrete element simulation of global ice loads on ships; and （c） Analytical solutions and numerical models of impact loads of icebergs on ships for polar navigation. Finally, research potential in these areas is discussed.
基金The study was supported by the National Natural Science Foundation of China under contract No.10772041the State 0ceamic Administration Key Laboratory for Ploar Science of China under contract No.KP 2007004.
文摘A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.
基金This study was supported financially by the National Key Research and Development Program of China(Grant no.2018YFA0605902)the National Natural Science Foundation of China(Grant no.52101300)+1 种基金the Fundamental Research Funds for the Central Universities(Grant no.DUT21LK03)Joint Scientific Research Fund Project of DBJI(Grant no.ICR2102).
文摘The splitting test is a competitive alternative method to study the tensile strength of sea ice owing to its suitability for sampling.However,the approach was questioned to the neglect of local plastic deformation during the tests.In this study,splitting tests were performed on sea ice,with 32 samples subjected to the regular procedure and 8 samples subjected to the digital image correlation method.The salinity,density,and temperature were measured to determine the total porosity.With the advantage of the digital image correlation method,the full-field deformation of the ice samples could be determined.In the loading direction,the samples mainly deformed at the ice-platen contact area.In the direction vertical to the loading,deformation appears along the central line where the splitting crack occurs.Based on the distribution of the sample deformation,a modified solution was derived to calculate the tensile strength with the maximum load.Based on the modified solution,the tensile strength was further calculated together with the splitting test results.The results show that the tensile strength has a negative correlation with the total porosity,which agrees with previous studies based on uniaxial tension tests.
文摘The quasi-solid-fluid phase transition of granular materials is closely related to the shear rate and solid concentration in addition to their intrinsic properties. The contact duration and the coordination number are two important temporal-spatial parame-ters to describe the granular interaction in phase transition. In this study, characteristics of the contact duration and the coordination number associated with the transition processes are determined using a 3D discrete element model under different shear rates and concentrations. The resulting macroscopic stress and strain-rate relations are discussed. The temporal and spatial parameters provide a linkage between the macroscopic constitutive law and inter- particle micromechanics.
基金supported by the National Basic Research Program of China (Grant No. 2010CB731502)the Fundamental Research Funds forthe Central Universities (Grant No. DUT12YQ02)
文摘The quasi-solid-liquid phase transition exists widely in different fields,and attracts more attention due to its instinctive mechanism.The structure of force chains is an important factor to describe the phase transition properties.In this study,the discrete element model(DEM) is adopted to simulate a simple granular shear flow with period boundary condition on micro scale.The quasi-solid-liquid phase transition is obtained under various volume fractions and shear rates.Based on the DEM results,the probability distribution functions of the inter-particle contact force are obtained in different shear flow phases.The normal,tangential and total contact forces have the same distributions.The distribution can be fitted as the exponential function for the liquid-like phase,and as the Weibull function for the solid-like phase.To describe the progressive evolution of the force distribution in phase transition,we use the Weibull function and Corwin-Ngan function,respectively.Both of them can determine the probability distributions in different phases and the Weibull function shows more reasonable results.Finally,the force distributions are discussed to explain the characteristics of the force chain in the phase transition of granular shear flow.The distribution of the contact force is an indicator to determine the flow phase of granular materials.With the discussions on the statistical properties of the force chain,the phase transition of granular matter can be well understood.
