Three-dimensional peridynamics based on matrix operation and its application in rock mass compression failure simulation

IF 5.3 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics Pub Date : 2025-05-23 DOI:10.1016/j.compgeo.2025.107354

Junyan Yang , Chenglu Gao , Zongqing Zhou , Daosheng Zhang , Fanlin Bu , Xiaochu Chen , Jinming Song

{"title":"Three-dimensional peridynamics based on matrix operation and its application in rock mass compression failure simulation","authors":"Junyan Yang , Chenglu Gao , Zongqing Zhou , Daosheng Zhang , Fanlin Bu , Xiaochu Chen , Jinming Song","doi":"10.1016/j.compgeo.2025.107354","DOIUrl":null,"url":null,"abstract":"<div><div>Discontinuous structural planes such as joints and cracks strongly influence the stability of surrounding rock in underground engineering. To simulate the failure behaviour of jointed and fractured rock masses under tension and compression complex loads, a peridynamics strength reduction constitutive model of jointed and fractured rock mass is proposed. Furthermore, the Weibull distribution function is introduced into the basic motion equations of peridynamics. This approach solves the problems of material point overlap and penetration in the simulations of the nonuniform damage and compression process of rock mass materials in peridynamics. By transforming the traditional peridynamics solution method into a matrix operation format, the efficiency of the peridynamics solution is improved significantly. Compared with the conventional successive retrieval iterative solution method, the proposed method can improve the calculation efficiency by approximately 30 %. Finally, numerical simulations of the failure process of intact rock masses, jointed rock masses and fractured rock masses are carried out. The results of the numerical simulations are compared with the results of previous laboratory tests and numerical simulations to validate the peridynamics strength-reduced constitutive model and the matrix operation method. The research in this paper is highly important for surrounding rock stability analysis.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107354"},"PeriodicalIF":5.3000,"publicationDate":"2025-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computers and Geotechnics","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0266352X25003039","RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}

引用次数: 0

Abstract

Discontinuous structural planes such as joints and cracks strongly influence the stability of surrounding rock in underground engineering. To simulate the failure behaviour of jointed and fractured rock masses under tension and compression complex loads, a peridynamics strength reduction constitutive model of jointed and fractured rock mass is proposed. Furthermore, the Weibull distribution function is introduced into the basic motion equations of peridynamics. This approach solves the problems of material point overlap and penetration in the simulations of the nonuniform damage and compression process of rock mass materials in peridynamics. By transforming the traditional peridynamics solution method into a matrix operation format, the efficiency of the peridynamics solution is improved significantly. Compared with the conventional successive retrieval iterative solution method, the proposed method can improve the calculation efficiency by approximately 30 %. Finally, numerical simulations of the failure process of intact rock masses, jointed rock masses and fractured rock masses are carried out. The results of the numerical simulations are compared with the results of previous laboratory tests and numerical simulations to validate the peridynamics strength-reduced constitutive model and the matrix operation method. The research in this paper is highly important for surrounding rock stability analysis.

查看原文本刊更多论文

基于矩阵运算的三维周动力学及其在岩体压缩破坏模拟中的应用

在地下工程中，节理、裂缝等不连续结构面对围岩的稳定性影响很大。为了模拟节理裂隙岩体在拉压复合载荷作用下的破坏行为，提出了节理裂隙岩体的周动力强度折减本构模型。在此基础上，将威布尔分布函数引入周期动力学的基本运动方程。该方法解决了围动力模拟岩体材料非均匀损伤与压缩过程中存在的质点重叠和侵彻问题。将传统的周期动力学求解方法转化为矩阵运算形式，显著提高了周期动力学求解效率。与传统的逐次检索迭代求解方法相比，该方法可将计算效率提高约30%。最后，对完整岩体、节理岩体和裂隙岩体的破坏过程进行了数值模拟。将数值模拟结果与室内试验结果和数值模拟结果进行对比，验证了围动力强度降维本构模型和矩阵运算方法的有效性。本文的研究对围岩稳定性分析具有重要意义。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

求助全文

约1分钟内获得全文求助全文

来源期刊

Computers and Geotechnics 地学-地球科学综合

CiteScore

9.10

自引率

15.10%

发文量

438

审稿时长

45 days

期刊介绍： The use of computers is firmly established in geotechnical engineering and continues to grow rapidly in both engineering practice and academe. The development of advanced numerical techniques and constitutive modeling, in conjunction with rapid developments in computer hardware, enables problems to be tackled that were unthinkable even a few years ago. Computers and Geotechnics provides an up-to-date reference for engineers and researchers engaged in computer aided analysis and research in geotechnical engineering. The journal is intended for an expeditious dissemination of advanced computer applications across a broad range of geotechnical topics. Contributions on advances in numerical algorithms, computer implementation of new constitutive models and probabilistic methods are especially encouraged.