Author(s): Jiabao Li; Qing Wang; Yiheng Zhang; Chongyang Jing; Xuesong Yang

Linked Author(s):

Keywords: peridynamics; elastic anisotropy; non-spherical influence function; single crystal; auxiliary coefficients

Abstract:

To accurately describe the elastic anisotropic behavior of ice Ih single crystals, in this paper, based on peridynamics (PD), the non-spherical influence function was adopted to extend the existing method into the three-dimensional (3D) space. Firstly, by rotating the coordinate system, the transformation relationship between the global and local coordinate systems was obtained, so that the mapping relationship between the crystal orientation and Young’s modulus can be established. Subsequently, the auxiliary coefficients were introduced into the mapping relationship as a non-spherical influence function and incorporated into the PD equilibrium equation to establish a numerical model for uniaxial compression of single crystal in PD. Finally, the Young’s modulus of a specific oriented single crystal model was obtained through the stress-strain relationship, and the auxiliary coefficients were iteratively calibrated by com- paring it with the theoretical value. After calibrating the auxiliary coefficients, it was found that the maximum error of the Young’s modulus represented by the 3D and 2D single crystal numerical model at each angle did not exceed 5% and 3%, respectively. All other errors were minimal, confirming that this method can be applied to characterize the elastic anisotropy of single crystal in space. Subsequently, a 2D columnar ice tensile model was established to investigate the effect of fracture criteria on tensile strength. The strength ranges corresponding to different fracture criteria were obtained. By comparing with previous experimental results, the effectiveness of the proposed elastic anisotropy characterization model was verified, which can be used for the construction of polycrystalline models

DOI: https://doi.org/10.5281/zenodo.14529930

Year: 2024