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":"10.1016/j.compgeo.2025.107354","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.3,"publicationDate":"2025-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144116388","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A first streamline-based simulation method within the projection-based embedded discrete fracture model (pEDFM)","authors":"Xiang Rao , Shuqing Guo , Xupeng He , Hyung Kwak , Hussein Hoteit","doi":"10.1016/j.compgeo.2025.107357","DOIUrl":"10.1016/j.compgeo.2025.107357","url":null,"abstract":"<div><div>This work presents the first streamline (SL) simulation method within the framework of a projection-based embedded discrete fracture model (pEDFM). The implemented pEDFM is formulated by combining two-point flux approximation (TPFA) and mimetic finite difference (MFD) schemes. We refer to the proposed numerical framework as SL-based pEDFM using hybrid TPFA-MFD scheme. It is operated in an IMPES manner. Specifically, the hybrid TPFA-MFD scheme is adopted to implicitly solve the pressure equation. The SL tracking workflow within the pEDFM framework is developed and applied to solve the transport equation in parallel along each streamline. We benchmark the proposed method with other existing streamline- or finite-volume-based approaches on various numerical examples. The results show that the proposed method exhibits broader general applicability compared to the SL-based EDFM, higher accuracy and efficiency than the TPFA-based pEDFM, and greater flexibility in mesh generation than the DFM-based approaches. The advantages underscore the great potential of the proposed method to be implemented in field-scale diagnostics and simulations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107357"},"PeriodicalIF":5.3,"publicationDate":"2025-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144116387","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Seismic performance of buried iron pipeline considering spatial variability of soil properties in 2D model","authors":"Benwei Hou, Chenzhao Xu, Qianyi Xu, Junyan Han, Zilan Zhong","doi":"10.1016/j.compgeo.2025.107347","DOIUrl":"10.1016/j.compgeo.2025.107347","url":null,"abstract":"<div><div>Iron pipes connected by bell-spigot joints are utilized in buried pipeline systems for urban water and gas supply networks. The joints are the weak points of buried iron pipelines, which are particularly vulnerable to damage from excessive axial opening during seismic motion. The axial joint opening, resulting from the relative soil displacement surrounding the pipeline, is an important indicator for the seismic response of buried iron pipelines. The spatial variability of soil properties has a significant influence on the seismic response of the site soil, which subsequently affects the seismic response of the buried iron pipeline. In this study, two-dimensional finite element models of a generic site with explicit consideration of random soil properties and random mechanical properties of pipeline joints were established to investigate the seismic response of the site soil and the buried pipeline, respectively. The numerical results show that with consideration of the spatial variability of soil properties, the maximum axial opening of pipeline joints increases by at least 61.7 %, compared to the deterministic case. Moreover, in the case considering the variability of pipeline-soil interactions and joint resistance, the spatial variability of soil properties remains the dominant factor influencing the seismic response of buried iron pipelines.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107347"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144106829","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xin Wei , Dengxing Qu , Zhengrong Zhou , Xinping Li , Yingwei Zhu , Shaohua Hu , Wenhao Li
{"title":"Three-dimensional refined discrete element numerical modeling method and its application for reef limestone based on digital core technology","authors":"Xin Wei , Dengxing Qu , Zhengrong Zhou , Xinping Li , Yingwei Zhu , Shaohua Hu , Wenhao Li","doi":"10.1016/j.compgeo.2025.107362","DOIUrl":"10.1016/j.compgeo.2025.107362","url":null,"abstract":"<div><div>The increasing scale of reef engineering construction has heightened the importance of understanding reef limestone’s mechanical properties in load-bearing strata. High-precision CT scanning provided digital core data, enabling precise characterization of pore-matrix structures through multi-level filtering and adaptive threshold segmentation. To enhance spatial mapping accuracy and computational efficiency in high-resolution modeling, we developed a three-dimensional refined discrete element method incorporating BallTree algorithm and two-level Euclidean distance logical filtering. The model was validated through pore characteristic comparison and mechanical testing before conducting systematic uniaxial compression simulations. Analysis revealed that shallow weakly-cemented compact reef limestone exhibits vesicular and banded pores, demonstrating microscopic heterogeneity alongside macroscopic homogeneity. Under uniaxial loading, the limestone undergoes five distinct stages: pore compaction, linear elasticity, crack propagation, peak failure, and residual strength. Stress transmission occurs preferentially along the 45° direction, governing crack development. The failure process initiates with tensile cracking and evolves into a combined tensile-shear failure mode. Energy analysis indicates that elastic strain energy dominates storage, while sliding friction heat represents the primary dissipation mechanism during failure. This study integrates digital core with numerical simulation to elucidate the relationship between pore structure and failure evolution in reef limestone, offering new perspectives on deformation and failure mechanisms in porous rocks.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107362"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144106702","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jinghua Zhang , Zhihua Yang , Emilio Bilotta , Lucia Mele , Qing Sun , Haitao Yu , Yong Yuan
{"title":"Soil improvement with laponite: effects on soil-structure interaction in liquefiable sands","authors":"Jinghua Zhang , Zhihua Yang , Emilio Bilotta , Lucia Mele , Qing Sun , Haitao Yu , Yong Yuan","doi":"10.1016/j.compgeo.2025.107360","DOIUrl":"10.1016/j.compgeo.2025.107360","url":null,"abstract":"<div><div>This paper investigates the effects of laponite treatment as a mitigation technique against seismically induced soil liquefaction and its influence on soil-structure interaction. Firstly, undrained cyclic triaxial tests are conducted on clean and laponite-treated Hostun sands. The tests confirm that laponite could enhance the resistance of the treated sand against liquefaction and that the constitutive model PM4Sand is able to capture the main aspects of the dynamic behavior of both clean and laponite-treated Hostun sands, by adjusting the relevant constitutive parameters. Then, two-dimensional plane-strain finite element models are constructed using OpenSees and PM4Sand based on the centrifuge tests involving two schemes of laponite treatment. The numerical results match well with the centrifuge test data, thus verifying the validity of the numerical soil-structure systems. Finally, the numerical models are re-configured to test the effects of laponite treatment on the soil-structure interaction of a shallow buried tunnel, a surface structure, and a combination of the two in liquefiable soils.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107360"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144116385","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A bottom-up hierarchical model for MICP-treated soil based on DEM","authors":"Haiyang Zhao , Annan Zhou , Shui-Long Shen , Arul Arulrajah","doi":"10.1016/j.compgeo.2025.107345","DOIUrl":"10.1016/j.compgeo.2025.107345","url":null,"abstract":"<div><div>Microbially induced calcite precipitation (MICP) is a promising soil improvement technique. This study presents a novel discrete element model (DEM) for the MICP-treated soil, featuring a bottom-up hierarchical approach encompassing three different levels. At level I (the CaCO<sub>3</sub>-sand system), the model incorporates realistic representations of CaCO<sub>3</sub> crystals on a soil particle surface, capturing their size, quantity, and distribution. At level II (representative sand pair with CaCO<sub>3</sub> cementation), the model investigates the influence of interparticle gaps, soil particle sizes, and CaCO<sub>3</sub> bond strength and mass fraction on the CaCO<sub>3</sub> cementation properties between soil particles. This is achieved by upscaling the model from level I to a representative pair. The analysis at level II leads to the introduction of generalised equivalent cementation bonds (ECB) that effectively capture the mechanical behaviour of interparticle CaCO<sub>3</sub> crystals and their dependence on sand properties. At level III (MICP-treated soils), large-scale simulations of soil particles are utilised to examine the effects of CaCO<sub>3</sub> mass fraction and soil void ratio on the mechanical behaviours of soil samples, based on the generalised ECB developed at level II. Additionally, to account for the CaCO<sub>3</sub> heterogeneity of MICP-treated soil, random field theory is utilised to generate spatially varied CaCO<sub>3</sub> distribution at this level. The results of DEM analysis indicate that CaCO<sub>3</sub> mass fraction significantly enhances soil strength, while its influence on elastic modulus is less pronounced. Furthermore, a reduction in void ratio increases interparticle bonds, leading to improved mechanical performance. This effect is further amplified by higher CaCO<sub>3</sub> mass fractions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107345"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144106830","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A discrete element-based study of freeze–thaw damage in water-filled fractured rock","authors":"Fengqi Shen , Wenliang Qiu , Kai Xie , Mingming Xing , Haowei Zhu , Shibo Huang","doi":"10.1016/j.compgeo.2025.107343","DOIUrl":"10.1016/j.compgeo.2025.107343","url":null,"abstract":"<div><div>Freeze–thaw damage to fractured rock severely threatens the safety of geotechnical engineering. The existing studies lack numerical modeling of freeze–thaw cycles in fractured rock. In this study, a discrete element model (DEM)-based method for simulating freeze–thaw damage in water-bearing fractured rock is proposed. Frost expansion of pore water and fracture water is simulated using the particle expansion method and effective volume expansion method, respectively. Changes in the water content during freeze–thaw cycles are considered. The proposed model is validated by comparison with laboratory results. The evolution of freeze–thaw damage and cracking behavior are analyzed for different numbers of freeze–thaw cycles and different fracture inclinations. The numerical results reveal that the effect of freeze–thaw damage caused by fracture water is greater than that of damage caused by pore water after a few freeze–thaw cycles. After 40 freeze–thaw cycles, the peak stress of the fractured rock decreases by 67%. The cracking resistance range of the fracture controls the freeze–thaw damage, and for standard samples, the freeze–thaw damage of the 0° fracture inclination rock is the greatest. The fracture inclination first controls the direction of fracture propagation and then shifts to the localized damage zone. This numerical approach provides a meaningful tool for understanding freeze–thaw damage and frost cracking in engineering application in cold regions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107343"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144106701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tianju Wang , Degao Zou , Jingmao Liu , Kai Chen , Xiuyang Zhang , Guoyang Yi
{"title":"Numerical simulation of environmental and seismic loading of 10 MW monopile offshore wind turbines in liquefiable seabed","authors":"Tianju Wang , Degao Zou , Jingmao Liu , Kai Chen , Xiuyang Zhang , Guoyang Yi","doi":"10.1016/j.compgeo.2025.107358","DOIUrl":"10.1016/j.compgeo.2025.107358","url":null,"abstract":"<div><div>With the gradual application of 10 MW offshore wind turbines (OWTs) in the world, it is very important to study the liquefaction response of high power OWTs under multiple fields. In this paper, based on the state-dependent generalized plasticity model, a cross-scale refined analysis model of structure-monopile-liquefied seabed is established, and a method for the analysis of seismic dynamic response of large-diameter monopile offshore wind turbine (MOWT) in liquefiable seabed is developed. The applicability of the proposed method is validated through simulations of existing centrifuge model tests. Using this approach, the dynamic response characteristics of a 10 MW MOWT under combined wind, wave, and seismic loading are systematically investigated. The study reveals the spatial distribution of seabed liquefaction under varying earthquake intensities, elucidates the flow mechanisms, and explores the effects of different loading conditions and relative depth (<em>Rd</em>) of liquefiable sandy soil on the system’s response. The results show that: 1. Seismic loads significantly influence seabed liquefaction, while environmental loads primarily cause cumulative rotation in MOWT. 2. As <em>Rd</em> increases, the rate of pore pressure accumulation slows, and liquefaction depth, soil strain around the pile, MOWT deflection, and rotation angle all increase. 3. Liquefaction results in overall subsidence of the soil within a range of one pile diameter (1D, D = 10.375 m) around the pile, while soil flow and heave within a range of 2D behind the pile are the main contributors to foundation rotational failure.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107358"},"PeriodicalIF":5.3,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144106703","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Earth pressure on rotating-about-top excavation retaining walls adjacent to existing structures","authors":"Xin-yi Huang, Chang Chen, Fu-quan Chen, Yi-xing Kuang","doi":"10.1016/j.compgeo.2025.107335","DOIUrl":"10.1016/j.compgeo.2025.107335","url":null,"abstract":"<div><div>The load transfer mechanism and failure mode of narrow soils induced by rotating about the top of excavation retaining structures (RT displacement mode) are investigated using the Finite Element Limit Analysis (FELA) method with the Hardening Mohr-Coulomb (HMC) model. A logarithmic spiral curve model is innovatively proposed to accurately characterize the evolution of the failure surface in narrow soils behind retaining structures under RT displacement mode based on numerical simulation results. Through analysis of principal stress vector diagrams obtained from FELA simulations, an asymmetric soil arching effect is identified in the upper zone of the narrow soils behind the retaining structure. An optimized differential element method is developed by constructing asymmetric arched differential elements along the deflection of principal stresses at the soil’s limit state, thereby establishing an analytical framework for calculating earth pressure in narrow soils under RT displacement mode. The proposed analytical method is validated through favorable agreement with finite element analysis results. Additionally, sensitivity analyses are performed to evaluate the effects of soil strength parameters, interface friction angles, and aspect ratios on earth pressure distribution, earth pressure coefficients, and the location of the resultant thrust.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107335"},"PeriodicalIF":5.3,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144089350","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Numerical investigation into the formation of pipe flow and the influencing factors in soil-rock mixtures","authors":"Yang Yu , Zihan Xu , Yongyong Yang , Ni An","doi":"10.1016/j.compgeo.2025.107342","DOIUrl":"10.1016/j.compgeo.2025.107342","url":null,"abstract":"<div><div>Soil-rock mixture (S-RM) slopes are widely distributed in mountainous regions, and the occurrence of landslides poses threats to residents and the economy. Hence, it is of great importance to investigate the stability of the S-RM slope. The existence of a pipe drainage system plays a key role in the assessment of slope stability. This work investigates the formation and development of a pipe drainage system in the S-RM during the seepage process. A new numerical approach based on a characteristic geometric model of the S-RM and a multi-physical field finite element coupling model of seepage-erosion-stress is developed, which is further validated by the experimental results. Moreover, the effects of the initial porosity distribution, volumetric block proportion (<em>VBP</em>), and hydraulic gradient on the development of drainage systems are further studied. The results show that the formation of the seepage channel is influenced by both the initial porosity and rock distribution in S-RM. The drainage channels are generated preferentially in the zone with higher initial porosity and the zone where the seepage direction is not blocked by rocks. With the increase of <em>VBP</em>, both the number and width of the drainage channels decrease. A higher hydraulic gradient can improve the formation of the seepage channels in terms of the number of channels as well as the width of channels. This work enhances the understanding of the evolution of pipe drainage systems in S-RM and provides further insights into the potential control measures in landslide treatments.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"185 ","pages":"Article 107342"},"PeriodicalIF":5.3,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144099881","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}