A NUMERICAL SIMULATION METHOD FOR THE ELASTIC ANISOTROPY OF SINGLE CRYSTAL ICE BASED ON PERIDYNAMICS

The anisotropy in the deformation and failure of natural ice, originating from the anisotropy of single crystal ice, is the main reason of the complex loading process during the ice-structure interaction. However, studies on the numerical simulation method of the anisotropy of single crystal ice are still rare in the academic community. To simulate such mechanical property of ice, a numerical simulation method for the elastic anisotropy of single crystal ice is proposed in this paper based on the theory of peridynamics. In the present method, the variation of Young's modulus of single crystal ice along different loading directions with respect to the c-axis obtained from published experimental results is adopted in the influence function of the force density vector in the state-based peridynamic model. Based on the numerical simulations of the uniaxial compression of single crystal ice along the loading directions of 0°, 45° and 90° with respect to the c-axis, correction method for the influence function of the peridynamic model, as well as the calibration procedure for the related auxiliary parameters, are proposed in this paper. Furthermore, validations of the Young's modulus for other loading directions including 15°, 30°, 60° and 75° are made, and good agreement has been achieved according to the comparison between the numerical and the reference experimental values. The results show that the correction method and calibration procedure presented in this paper can efficiently find the optimal solution to the influence function for the consistency between the Young's modulus of the numerical model and the reference Young's modulus from the published experiments, which indicates that the proposed numerical method based on peridynamic theory can sensibly simulate the elastic anisotropy behavior of single crystal ice. The main findings in this paper can provide basic reference for the future development of the numerical simulation method for the anisotropy of polycrystalline ice.