EI、Scopus 收录
中文核心期刊
黄焱, 王建平, 孙剑桥. 基于近场动力学的单晶冰弹性各向异性数值模拟方法. 力学学报, 2022, 54(6): 1641-1650. DOI: 10.6052/0459-1879-22-064
引用本文: 黄焱, 王建平, 孙剑桥. 基于近场动力学的单晶冰弹性各向异性数值模拟方法. 力学学报, 2022, 54(6): 1641-1650. DOI: 10.6052/0459-1879-22-064
Huang Yan, Wang Jianping, Sun Jianqiao. A numerical simulation method for the elastic anisotropy of single crystal ice based on peridynamics. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(6): 1641-1650. DOI: 10.6052/0459-1879-22-064
Citation: Huang Yan, Wang Jianping, Sun Jianqiao. A numerical simulation method for the elastic anisotropy of single crystal ice based on peridynamics. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(6): 1641-1650. DOI: 10.6052/0459-1879-22-064

基于近场动力学的单晶冰弹性各向异性数值模拟方法

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

  • 摘要: 天然冰材料在变形与破坏行为上的各向异性特征是冰与结构相互作用中产生复杂载荷过程的关键诱因, 而天然冰各向异性的根源则在于单晶冰的各向异性. 目前, 学术界针对单晶冰各向异性的数值模拟方法研究仍较为缺乏. 为了准确再现天然冰材料的特殊力学性质, 本文基于近场动力学理论, 提出了一种单晶冰弹性各向异性的数值模拟方法. 该方法的核心思想是将单晶冰杨氏模量沿不同加载方向的变化规律引入到近场动力学力密度向量的影响函数中. 以前人实验测试得到的杨氏模量值为参考, 通过开展与C轴呈0°, 45°和90°三个加载方向的单晶冰单轴压缩数值模拟实验, 提出了针对该影响函数的修正和辅助参数标定方法, 最终在15°, 30°, 60°和75°等其他四个加载方向进行了验证. 结果表明: 本文提出的针对影响函数的修正与参数标定方法, 能够较为便捷地找到数值模型杨氏模量与参考杨氏模量相一致的影响函数最优解, 即本文提出的基于影响函数的近场动力学数值模拟方法, 能够合理、准确地模拟单晶冰的弹性各向异性行为. 本文研究成果可为后续多晶冰各向异性数值模拟方法的建立提供基础性参考.

     

    Abstract: 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.

     

/

返回文章
返回