{"title":"Cyclic p-y curve modeling for laterally loaded two-pile group in sand based on strain wedge model","authors":"Fanhuan Zeng, Chong Jiang","doi":"10.1016/j.compgeo.2025.107499","DOIUrl":"10.1016/j.compgeo.2025.107499","url":null,"abstract":"<div><div>The lateral response of pile groups used to support offshore structures differs significantly from that of single piles. The cyclic lateral loads induced by waves, tides, and winds complicate the study of the calculation method for cyclic pile response. A cyclic <em>p-y</em> curve model is developed to address the cyclic lateral response of a two-pile group in this study. The strain wedge model is introduced and modified to account for the shadow effect of a laterally loaded pile group and to derive the backbone <em>p-y</em> curves. The cyclic <em>p-y</em> curves framework for a two-pile group is constructed using the improved collapse and recompression model to represent cyclic pile-soil interaction in sand. A program is developed for the model based on the finite difference method to solve the incremental pile governing equation. The model is evaluated by comparison with the results of existing centrifuge model tests and the discrete element method-finite difference method (DEM-FDM) coupled simulation of this work. The proposed model is effective in predicting the cyclic lateral response of two-pile group. Based on the results, the characteristics of the shadow effect, cyclic <em>p-y</em> curves, and sand movement are discussed.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107499"},"PeriodicalIF":5.3,"publicationDate":"2025-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144657241","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}
Wenjun Ji , Qingshuo Hao , Chongmin Song , Qianhui Liu , Yuzhen Yu
{"title":"A stabilized equal-order linear polygonal and polyhedral scaled boundary finite element formulation for Biot’s consolidation problems in saturated soils","authors":"Wenjun Ji , Qingshuo Hao , Chongmin Song , Qianhui Liu , Yuzhen Yu","doi":"10.1016/j.compgeo.2025.107487","DOIUrl":"10.1016/j.compgeo.2025.107487","url":null,"abstract":"<div><div>This paper addresses the numerical instability of pore pressure in solving Biot’s consolidation problems by developing a stabilized equal-order solution framework based on the scaled boundary finite element method (SBFEM). First, a concise formulation of Laplace generalized shape function is constructed for arbitrary polygonal and polyhedral elements in both two and three dimensions. A novel stabilization strategy is then proposed by introducing a pressure projection scheme into SBFE, resulting in a stabilized low-order formulation that ensures numerical stability while retaining high computational efficiency. Numerical examples demonstrate that the proposed method significantly improves pressure field stability and exhibits strong compatibility with the unified hardening (UH) elastoplastic constitutive model in saturated soils. The approach broadens the applicability of SBFEM to more complex coupled problems, enhances accuracy and robustness, and provides a new theoretical tool for analyzing challenging geotechnical engineering scenarios.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107487"},"PeriodicalIF":5.3,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144657237","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}
Qingzhe Cui , Fei Wu , Liting Wang , Jianfeng Liu , Jie Chen , Cunbao Li , Huiqing Liu , Renbo Gao
{"title":"A rock elasto-visco-plastic creep model fully considering the effect of time-dependent damage","authors":"Qingzhe Cui , Fei Wu , Liting Wang , Jianfeng Liu , Jie Chen , Cunbao Li , Huiqing Liu , Renbo Gao","doi":"10.1016/j.compgeo.2025.107490","DOIUrl":"10.1016/j.compgeo.2025.107490","url":null,"abstract":"<div><div>The establishment of a constitutive model that accurately describes the creep-damage coupling behavior of rock is a theoretical cornerstone for predicting the delayed deformation of surrounding rock. Accurately assessing the impact of time-dependent damage on rock creep characteristics is therefore of paramount importance. In this study, following the Lemaitre strain equivalence principle, time-dependent damage is simultaneously introduced into the Kelvin body, representing the viscoelastic behavior of rock, and the Bingham body, representing the viscoplastic behavior, to establish a rock elasto-visco-plastic creep model (TDEVP) that fully considers the effect of time-dependent damage. The model’s validity is confirmed through comparisons with existing laboratory creep data. The results show that, during the accelerated creep stage, the time-dependent deformation of the viscoelastic component exceeds that of the viscoplastic component, with the viscoelastic strain being up to 2.8 times greater. Furthermore, the TDEVP model was successfully implemented in FLAC3D. Compared to the theoretical model, the secondary development program exhibits stronger feasibility and applicability. Finally, the model was applied to a deeply buried tunnel project. The TDEVP model achieved an <em>L</em>2<sub>P</sub> error of only 0.58 mm and outperformed several classical models across five evaluation metrics, demonstrating its robustness and accuracy. This study addresses key limitations in existing creep-damage models and provides a solid theoretical basis for improving the long-term stability prediction of deep underground engineering.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107490"},"PeriodicalIF":5.3,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144657239","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}
Pengcheng Yu , Jie Li , Yukun Long , Xinyan Peng , Bing Hu , Yingbin Zhang , Dejian Li , Jinmei Wang , Shihong Bai
{"title":"Exploring fracturing and post-failure process of rock mass using modified discontinuous deformation method considering non-linear behavior of joints","authors":"Pengcheng Yu , Jie Li , Yukun Long , Xinyan Peng , Bing Hu , Yingbin Zhang , Dejian Li , Jinmei Wang , Shihong Bai","doi":"10.1016/j.compgeo.2025.107495","DOIUrl":"10.1016/j.compgeo.2025.107495","url":null,"abstract":"<div><div>Discontinuous deformation analysis (DDA) has gained wide acceptance in the field of rock mechanics. Nevertheless, the non-linear behavior of virtual or real joints in DDA can be further considered since the joints in nature commonly presents non-linear features. Therefore, a modified DDA considering non-linear behavior of joints is presented. A cohesive contact model to bond adjacent blocks through virtual joints both in the normal and shear directions is implemented into DDA to capture the fracturing behavior. This cohesive contact model can better describe the softening behavior of rock material after reaching the peak cohesive strength. For real joints due to the existence of real discontinuities such as the existing fractures and the newly propagated cracks, the nonlinear Barton-Bandis model considering an exponentially displacement-dependent joint roughness degradation is adopted. The non-linearity inherent in the Barton-Bandis model can better describe the joint shear behavior. Several validation examples are studied to show the success of implementing the two non-linear models. A series of mixed-mode Brazilian tests with different internal diameters are studied to further show the capability of the modified DDA to predict the rock fracturing behavior. Finally, the modified DDA is applied to simulate the earthquake-induced Donghekou landslide. The fracturing and post-failure process of rock mass can be well predicted.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107495"},"PeriodicalIF":5.3,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144657240","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":"Application of geophysical methods for 2D and 3D geomechanical modeling in geotechnical engineering","authors":"Muhammad Hasan , Lijun Su","doi":"10.1016/j.compgeo.2025.107496","DOIUrl":"10.1016/j.compgeo.2025.107496","url":null,"abstract":"<div><div>The successful execution of deep underground engineering projects depends heavily on comprehensive and accurate geomechanical modeling. Among various geomechanical parameters, the rock mass integrity coefficient (Kv) plays a critical role in assessing the stability and associated risks of subsurface structures. However, accurately determining Kv is often constrained by the high costs, time demands, and limited spatial coverage of conventional borehole-based methods. To address these limitations, this study presents a novel application of the Controlled Source Audio-frequency Magnetotellurics (CSAMT) method to generate 2D and 3D Kv models across a geologically heterogeneous site at depths reaching 1 km. This geophysical approach offers a non-invasive, cost-effective, and efficient alternative for rock mass characterization. The results contribute to a better understanding of complex subsurface conditions, enhance early-collapse risk assessment, and support safer, more economical design of deep underground infrastructures. In regions with limited geomechanical data, this method helps bridge the gap between sparse field measurements and reliable geomechanical modeling, ultimately enabling more accurate and objective geotechnical evaluations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107496"},"PeriodicalIF":5.3,"publicationDate":"2025-07-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144657238","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}
Dongxu Chen , Laigui Wang , Shanyong Wang , Han Gao , Chuang Sun
{"title":"A creep damage model and microscopic damage characteristics of fractured sandstone repaired by microbial induction","authors":"Dongxu Chen , Laigui Wang , Shanyong Wang , Han Gao , Chuang Sun","doi":"10.1016/j.compgeo.2025.107494","DOIUrl":"10.1016/j.compgeo.2025.107494","url":null,"abstract":"<div><div>The Microbially Induced Carbonate Precipitation (MICP) technique provides a novel solution to the long-term stability challenges of fractured sandstone. To evaluate the effectiveness of MICP in controlling the creep behavior of fractured rocks, a series of experiments were conducted to investigate the microscopic mineralization characteristics of MICP and the creep mechanical properties of repaired sandstone. A CRG-MICP numerical model was developed using the Discrete Element Method (DEM) under the assumption of uniformly distributed mineral particles during MICP. The model captures morphological characters of microbial-induced minerals at the microscale while simulating macroscale rock deformation and energy-based damage evolution, establishing a cross-scale framework for rock mechanical properties. This model was further utilized to determine parameters of the creep damage constitutive model and the evolution of microscopic structural damage in the rock. The results indicate that the circularity and diameter of the mineralized particles vary significantly with different MICP solution concentrations. As CaCl<sub>2</sub> concentration increases, the elastic, viscoelastic, and viscoplastic parameters initially rise and subsequently decline, whereas the damage parameter demonstrates an inverse trend (decreasing first and then increasing). The extent of microscopic structural damage in the repaired sandstone similarly follows a pattern of initial reduction followed by augmentation with increasing CaCl<sub>2</sub> concentration. These findings provide critical scientific and engineering insights for controlling rock creep damage.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107494"},"PeriodicalIF":5.3,"publicationDate":"2025-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633092","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":"Heat transfer pathways in unsaturated granular materials: A thermal conductance network model perspective","authors":"Tairu Chen , Wenbin Fei , Guillermo A. Narsilio","doi":"10.1016/j.compgeo.2025.107493","DOIUrl":"10.1016/j.compgeo.2025.107493","url":null,"abstract":"<div><div>Understanding heat transfer in unsaturated granular materials is important to geo-engineering applications, yet the pathways sustaining heat transfer between particles have not been quantitively identified. This work recognised six types of pathways for heat transfer in a granular assembly varying in degree of saturation: (a) particles in physical contact with air bridges, (b) particles in physical contact with both air and water bridges, (c) particles in physical contact with only water bridges, (d) separate particles “connected” by air bridges, (e) separate particles “connected” by both air and water bridges and (f) separate particles “connected” by only water bridges. To quantify each pathway’s contribution to the effective thermal conductivities, thermal conductance network models were developed. In these models, a granular assembly is depicted as a web of nodes (particles) connected by edges (physical and “near contact” via air/water bridges). The thermal conductance of each edge was calculated and used as a weight for that edge to determine the effective thermal conductivity. Results showed that effective thermal conductivities from thermal conductance network models matched the experimental data well. Furthermore, the influence of edge types (b), (e), and (f) on heat transfer is dominant in unsaturated conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107493"},"PeriodicalIF":5.3,"publicationDate":"2025-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633690","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}
Yujie Zhu , Gang Zheng , Haizuo Zhou , Boyang Xia , Xiaoxuan Yu , Yi He , Denis Kalumba
{"title":"Evaluation of settlement profiles of soft clay under embankments","authors":"Yujie Zhu , Gang Zheng , Haizuo Zhou , Boyang Xia , Xiaoxuan Yu , Yi He , Denis Kalumba","doi":"10.1016/j.compgeo.2025.107492","DOIUrl":"10.1016/j.compgeo.2025.107492","url":null,"abstract":"<div><div>Current design procedures for embankment-induced settlement usually assume that the maximum ground settlement occurs beneath the embankment center; however, actual deformation profiles exhibit various characteristics that result in significant discrepancies. In this study, a finite element analysis is implemented to quantify the effects of geometric parameters (the embankment width <em>B</em>, height <em>H</em>, and soft clay thickness <em>H</em><sub>s</sub>) and soil properties (the undrained shear strength ratio <em>S</em><sub>u</sub><em>/σ’</em><sub>v</sub>) on settlement patterns. The results show that the <em>H</em><sub>s</sub> and <em>S</em><sub>u</sub><em>/σ’</em><sub>v</sub> affect the settlement magnitude, whereas the <em>B</em> and <em>H</em> are coupled to the settlement patterns. Three characteristic types (typical, sagged, and transitional patterns) are recognized, with each type corresponding to a unique maximum settlement location. The relationship between the maximum settlement and lateral displacement under the different patterns is established by simplifying the deformation as a triangular distribution; this enables a derivation of the shear stress–induced settlement. A piecewise function describing the profiles of the different patterns is proposed. The mechanism of settlement pattern variation is discussed on the basis of the distribution of incremental settlement induced by shear stress. Classification criteria distinguishing the settlement patterns are then proposed. Finally, validations of the developed method in comparison with field measurements are made.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107492"},"PeriodicalIF":5.3,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144614914","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":"Reliability-based multi-objective optimisation for the design of vegetation planting strategy on landfill cover systems","authors":"Chuanxiang Qu , Charles W.W. Ng , Haowen Guo","doi":"10.1016/j.compgeo.2025.107483","DOIUrl":"10.1016/j.compgeo.2025.107483","url":null,"abstract":"<div><div>Existing design criteria lack specific recommendations for vegetation planting strategies of final landfill cover systems. Factors such as inherent uncertainties in cover material and plant properties are generally ignored in design. Moreover, the current design mainly focuses on ensuring effectiveness while ignoring the total cost and carbon impact. In this study, a reliability-based multi-objective optimisation framework is developed for designing vegetated final cover systems. This framework balances the probability of failure (<em>p</em><sub>f</sub>, the probability that percolation exceeds a recommended limit), total cost, carbon sequestration and emissions of final cover systems. A numerical model integrated with a stress-dependent soil water retention curve and the hydraulic influence of vegetation is validated through soil column tests. The optimal planting strategy, including the selection of root architecture, depth, and plant spacing, is then determined. It is found that a triangular root is more effective in reducing <em>p</em><sub>f</sub> (by up to 3 orders of magnitude) than a uniform root, regardless of plant spacing and root depth. Increasing plant spacing from 0.5 m to 1.5 m in increments of 0.25 m raises <em>p</em><sub>f</sub> by at least 2 orders of magnitude for both root architectures. Conversely, increasing root depth from 0.3 m to 1.0 m in increments of 0.2 m (associated with greater biomass) reduces <em>p</em><sub>f</sub> by at least 77 % with the two root architectures considered. The optimal design, among 32 planting scheme combinations investigated, is a triangular root architecture with a 0.5 m planting spacing and a 0.3 m root depth. This combination can achieve high performance while minimising total cost and net carbon emission (carbon emission minus carbon sequestration).</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107483"},"PeriodicalIF":5.3,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633051","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":"Failure mechanism of rock tunnels subjected to combined reversed fault dislocation and subsequent seismic excitation","authors":"Dong-Mei Zhang , Li Guo , Zhong-Kai Huang , Yi-Yao Shen , Kyriazis Pitilakis","doi":"10.1016/j.compgeo.2025.107485","DOIUrl":"10.1016/j.compgeo.2025.107485","url":null,"abstract":"<div><div>Rock tunnels are particularly vulnerable when crossing seismically active fault zones. While their response to fault dislocation or seismic excitation has been quite well studied, their response under the combined effects of fault dislocation and ground motion remains poorly examined, primarily due to its complexity. To study this challenging problem, the paper proposes a well-constrained numerical approach where a typical case of fault dislocation and subsequent seismic excitation has been modelled for typical tunnel lining and rock conditions. The model incorporates cohesive interface elements and a concrete damage plasticity model to represent fault discontinuities and the nonlinear weakening of tunnel linings. Validation was conducted using centrifuge tests under free-field fault dislocation. Parametric analyses were performed varying the dislocation characteristics of the selected reverse-fault and ground motions. The principle aim is threefold: (a) to propose a methodology to analyse the complex problem of the combined creeping fault dislocations and subsequent seismic excitation, illustrated through a well selected case study; (b) to investigate the complex mechanism and the associated deformations of the tunnel when subjected to the combined loading; (c) to describe an efficient approach for the studied case, with two damage assessment indicators, to assess the vulnerability and the damages of the tunnel when crossing the seismic fault. The study provides a valuable set of insightful data for the assessment of the seismic performance of mountain tunnels in rock conditions crossing seismically active fault zones, contributing to improving the seismic protection of tunnels.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"187 ","pages":"Article 107485"},"PeriodicalIF":5.3,"publicationDate":"2025-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144633050","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}