岩石三维裂纹扩展问题的近场动力学数值模拟  

Numerical Simulation of 3D Crack Propagation in Rock by Peridynamics Approach

作　　者：崔昊 郑宏 李春光[3] 韩月 CUI Hao;ZHENG Hong;LI Chunguang;HAN Yue(Power China Huadong Engineering Corporation Limited,Hangzhou 311122,China;Faculty of Architecture,Civil and Transportation Engineering,Beijing University of Technology,Beijing 100124,China;State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Wuhan 430071,China)

机构地区：[1]中国电建集团华东勘测设计研究院有限公司,浙江杭州311122 [2]北京工业大学城市建设学部,北京100124 [3]中国科学院武汉岩土力学研究所,湖北武汉430071

出　　处：《同济大学学报（自然科学版）》2022年第9期1223-1231,共9页Journal of Tongji University:Natural Science

摘　　要：近场动力学(PD)方法基于键的形式处理断裂问题,对裂纹尖端可自动追踪,因而在求解三维裂纹扩展问题时具备极大的优势,但该方法中也存在数值震荡与边界误差等问题。为解决上述问题,首先应用无网格伽辽金弱形式框架对非常规态基近场动力学方法进行了探讨;随后引入近场动力学微分算子(PDDO)近似,并对比分析了该近似与重构核粒子(RKPM)近似间的差异,提出了具备更高数值精度的RKPM-PD耦合算法,并给出了该算法的隐式迭代流程。最后,通过若干数值算例验证了该算法在求解三维裂纹扩展问题上的有效性。The peridynamics method has great advantages in solving crack propagation problems due to its automatic tracking of crack tips by breaking bonds,but it also faces the problems such as numerical oscillations and boundary errors.In order to solve the above defects,this paper first discussed the NOSB-PD method with the meshless Galerkin weak-form framework. Next,it introduced the peridynamics differential operator approximation and fully compared and analyzed the difference between the PDDO approximation and the reproducing kernel particle method(RKPM)approximation.Then,it proposed the RKPM-PD coupling algorithm with a higher numerical accuracy and gave the implicit iterative process of the RKPM-PD method. Finally,it verified the effectiveness of the algorithm in solving three-dimensional crack propagation problems by using several numerical examples.

关 键 词：近场动力学 无网格方法 伽辽金法 三维裂纹扩展 

分 类 号：O343[理学—固体力学]