Computers and Geotechnics最新文献

{"title":"Multiscale investigation of shear mechanical behaviour in cemented soil–rock mixtures (CSRMs) regulated by matrix–block welding states","authors":"Minghui Ren ,&nbsp;Hai Pu ,&nbsp;Guangsi Zhao ,&nbsp;Ruilin Li ,&nbsp;Lulu Liu ,&nbsp;Runhua Zhang","doi":"10.1016/j.compgeo.2025.107254","DOIUrl":"10.1016/j.compgeo.2025.107254","url":null,"abstract":"<div><div>Soil–rock mixtures (SRMs) are characterized by heterogeneous structural features that lead to multiscale mechanical evolution under varying cementation conditions. However, the shear failure mechanisms of cemented SRMs (CSRMs) remain insufficiently explored in existing studies. In this work, a heterogeneous three-dimensional (3D) discrete element model (DEM) was developed for CSRMs, with parameters meticulously calibrated to examine the role of matrix–block interfaces under different volumetric block proportions (VBPs). At the macroscopic scale, significant influences of the interface state on the peak strength of CSRMs were observed, whereas the residual strength was found to be largely insensitive to the interface cementation properties. Pronounced dilatancy behaviour was identified in the postpeak and residual phases, with a positive correlation with both interface cementation and VBP. Quantitative particle-scale analyses revealed substantial heterogeneity and anisotropy in the contact force network of CSRMs across different components. A highly welded interface was shown to reduce the number of interface cracks at the peak strength state while increasing the proportion of tensile cracks within the interface zone. Furthermore, the welding degree of the interface was found to govern the formation and morphology of shear cracking surfaces at the peak strength state. Nevertheless, a reconstruction method for the shear slip surface was proposed to demonstrate that, at the same VBP, the primary roughness of the slip surfaces remained consistent and was independent of the interface properties. Based on the extended simulations, the peak strength of the weakly welded CSRMs progressively decreased with increasing VBP, whereas further exploration of the enhanced residual strength is needed.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107254"},"PeriodicalIF":5.3,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143817713","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":"Experimental and MPM modelling of widened levee failure under the combined effect of heavy rainfall and high riverine water levels","authors":"Yanhao Zheng,&nbsp;Jinhui Li,&nbsp;Tanfang Zhu,&nbsp;Junru Li","doi":"10.1016/j.compgeo.2025.107259","DOIUrl":"10.1016/j.compgeo.2025.107259","url":null,"abstract":"<div><div>The practice of widening levees to mitigate frequent river flooding is globally prevalent. This paper addresses the pressing issue of sand-filled widened levee failures under the combined effect of heavy rainfall and high riverine water levels, as commonly observed in practice. The primary objective is to illuminate the triggering mechanism and characteristics of such levee failures using the well-designed physical model experiment and Material Point Method (MPM), thus guiding practical implementations. Experimentally, the macro-instability of the levee, manifested as slope failure within the sand-filled widened section, is primarily triggered by changes in the stress regime near the levee toe and continuous creep deformation. Upon failure initiation, the levee slope experiences a progressive failure mode, starting with local sliding, followed by global sliding, and ultimately transitioning into a flow-like behaviour, which characterises the slide-to-flow failure pattern. The slope failure along the interface between the original and new levees is the result of shear deformation rather than the cause. Parametric studies conducted using the calibrated MPM model reveal a critical threshold for the widening width, beyond which the volume of sliding mass and travel angle exhibit no further variation. Increasing the cohesion of the river sand used for levee widening demonstrates the most pronounced improvement in levee stability in the face of the combined effect of intense rainfall and elevated river levels. The MPM-based evaluation of common slope protection measures demonstrates the superior protective benefits of grouting reinforcement and impervious armour layer protection, providing valuable insights for reinforcement strategies in levee engineering applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107259"},"PeriodicalIF":5.3,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143817631","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":"2D orthogonal simulation method for multi-physics field evolution in material-structure heterogeneous rocks under microwave irradiation","authors":"Yueyang Li ,&nbsp;Lei He ,&nbsp;Huaiguang Xiao ,&nbsp;Tienan Wang ,&nbsp;Yanqi Su ,&nbsp;Weiqiang Zhu ,&nbsp;Chee Kiong Soh","doi":"10.1016/j.compgeo.2025.107237","DOIUrl":"10.1016/j.compgeo.2025.107237","url":null,"abstract":"<div><div>Due to the significant physical property differences among various diagenetic minerals, it is crucial to characterize the multi-physics field evolution process of microwaves in heterogeneous rocks. To overcome the issues of insufficient refinement in existing numerical studies, this paper proposes a two-dimensional (2D) orthogonal model considering structural-material heterogeneous rocks for the detailed simulation of the strong coupling process of microwave radiation in rocks. Based on a high-fidelity granite model with texture structure, a strong coupled numerical simulation method is developed from both section and facade perspectives. Furthermore, the electrical-thermal-stress-damage physical responses of rocks to microwave radiation are meticulously calculated. The results indicate that both the section model and the facade model exhibit a stress distribution with central compression and surrounding tension, while the damage evolves from the inside out and from the surface inward, respectively. The change in the proportion of the tensile zone and the transformation of the stress properties of the mineral are the response characteristics of damage. The order of mineral damage is biotite, feldspar and quartz, with damage more likely to occur in the feldspar phase. Additionally, the texture structure and physical property differences in heterogeneous rocks are the main driving forces to rock damage. This study provides a novel simulation method for the microwave interaction with rocks and offers new insights into the evolution of multi-physics fields.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107237"},"PeriodicalIF":5.3,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143817712","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":"On the role of aspect ratio on the geocell-reinforced soils: A coupled continuum-discontinuum study","authors":"Haibo Wang ,&nbsp;Ge Gao ,&nbsp;Mohamed A. Meguid ,&nbsp;Lulu Zhang","doi":"10.1016/j.compgeo.2025.107256","DOIUrl":"10.1016/j.compgeo.2025.107256","url":null,"abstract":"<div><div>Geocells are three-dimensional, interconnected cellular geosynthetics widely used to enhance the overall strength of soils. Their foldable structure can cause variations in pocket shape during installation, depending on the extent of extension. Understanding the impact of these shape variations is essential for optimizing reinforcement efficiency and reducing the associated geocell application costs. The aspect ratio, defined as the ratio of the cell’s transverse (welded) axis to the longitudinal (wall summit) axis, is proposed to evaluate the degree of extension of the most commonly utilized honeycomb-shaped geocell. A coupled continuum-discontinuum numerical method was employed to investigate the behavior of honeycomb-shaped geocell reinforced soils across various aspect ratios under confined compressive loading. The simulation results indicate that a geocell with an aspect ratio of 1.0 exhibits optimal reinforcement efficiency, and whereas reinforcement efficiency decreases as the aspect ratio deviates from 1.0 causing pocket geometries to flatten. The superior performance of rounded geocells is attributed to their enhanced ability to promote load-bearing in strong contact subnetworks. This results in denser packing structures, higher contact force anisotropy from a microscopic perspective, and greater confinement capacity against deformation from a macroscopic perspective.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107256"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808577","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 coupled mathematical model of microbial grouting reinforced seabed considering the response of wave-induced pore pressure and its application","authors":"Yujie Li ,&nbsp;Lizhong Wang ,&nbsp;Changrui Dong ,&nbsp;Guohui Feng ,&nbsp;Xingye Sun ,&nbsp;Zhen Guo","doi":"10.1016/j.compgeo.2025.107235","DOIUrl":"10.1016/j.compgeo.2025.107235","url":null,"abstract":"<div><div>Wave-induced instability of the sandy seabed has become an issue that cannot be ignored in the design of offshore foundations. Microbially induced calcium carbonate precipitation (MICP) is believed as a potential method to reinforce seabed against liquefaction. A bio-chemo-hydro-wave (B-C–H-W) coupled theoretical model of microbial grouting reinforced seabed considering the response of wave-induced pore pressure was developed in this study. The stability of MICP-reinforced seabed under wave actions with flower grouting and spot grouting conditions was also analyzed in detail. Results show that MICP reaction generates local calcium carbonate, which reduces the porosity and permeability coefficient, and enhances the stiffness of seabed, thereby reducing the amplitude of oscillating pore pressure. The mechanism of MICP reaction to resist seabed liquefaction include: increases the weight of seabed, provides a degree of cohesion, and reduces the amplitude of oscillating pore pressure. Different grouting methods determine the distribution pattern of calcium carbonate. Flower grouting has a larger effective protection length than spot grouting at the same conditions. There is a critical distance during double-tube grouting, it is recommended to arrange grouting tubes according to the critical distance and consider a certain safety factor in actual operation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107235"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808573","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 simulation on enhanced coalbed methane extraction by hot flue gas based on the ECM-DFM model","authors":"Ting Liu ,&nbsp;Jianhao Wang ,&nbsp;Baiquan Lin ,&nbsp;Xinhao Wang ,&nbsp;Qiming Huang ,&nbsp;Jianfeng Li","doi":"10.1016/j.compgeo.2025.107253","DOIUrl":"10.1016/j.compgeo.2025.107253","url":null,"abstract":"<div><div>The multi-scale characteristics of the storage and transportation structure of coal seams will lead to abnormal heat and mass transport phenomena. It is necessary to explore the thermo-hydro-mechanical (THM) coupling mechanism of coalbed methane displacement by hot flue gas (CBMDHFG) and the mechanical response behaviour of macro fractures (MFs). Therefore, this paper develops a THM-coupled dual-medium model, which includes both the equivalent continuous medium (ECM) representing the coal matrix and micro fractures and the discrete fracture medium (DFM) representing the MFs. The fracture network randomly generated by the MATLAB code is imported into COMSOL to realize the solution of the model. Based on this, the evolution of multiple fields and the response law of gas flow rate during the CBMDHFG process are studied. The results show that the hot flue gas first migrates in the high-permeability MFs and then evolves towards the surrounding low-permeability ECM system. The high temperature and high pressure of the hot flue gas effectively weaken the compressive stress applied to the MFs in the near field and significantly promote shear dilation and expansion. A high gas injection temperature may induce the compression and closure of MFs outside the high-temperature excitation zone. CO₂ and N₂ mainly contribute to the shear dilation of MFs in the near field and far field respectively. Appropriately increasing the gas injection pressure, temperature, and the proportion of CO₂ in the early stage and gradually increasing the proportion of N₂ in the later stage can achieve the long-term development of gas. Loading high anisotropic stress induces complex THM coupling effects, leading to the inhomogeneity of the aperture of MFs. The side closer to the maximum principal stress is more likely to undergo shear dilation, and the peak gas production is higher. The research results are expected to provide new insights for exploring the enhanced recovery mechanism of CBMDHFG and optimizing the production scheme.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107253"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808574","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 DEM sensitivity study on the effects of contact parameters on triaxial response for the development of a calibration method","authors":"Damon Nguyen ,&nbsp;Yuyan Chen ,&nbsp;Alejandro Martinez","doi":"10.1016/j.compgeo.2025.107241","DOIUrl":"10.1016/j.compgeo.2025.107241","url":null,"abstract":"<div><div>Discrete element modeling (DEM) is a useful tool for linking global responses of granular materials to underlying particle-level interactions. A DEM model capable of capturing realistic soil behavior must be calibrated to a reference dataset, typically consisting of laboratory experiments. Calibration of a DEM model often requires numerous simulations as contact parameters need to be iterated upon until the simulation results satisfactorily replicate the experimentally observed behaviors. This paper presents a sensitivity investigation that examines the effects of the contact parameters on the drained triaxial compression response of a poorly-graded sand. It then introduces a calibration procedure capable of providing contact parameters that satisfactorily reproduce the results of laboratory triaxial results in a few simulations. Results show that friction and rolling resistance coefficients jointly influence the mobilized peak and critical state friction angles, secant shear modulus, maximum dilation rate, total volumetric strain, and strain softening magnitude. These parameters also influence the mode of failure at contacts and the evolution of fabric anisotropy. The influence of μ_r or μ on the triaxial response and particle-level interactions is coupled, becoming more profound as the other parameter is increased. Contact stiffness is shown to influence the shear modulus and volumetric change behavior independently of μ and μ_r. An algorithm that estimates values for μ and μ_r needed to reproduce experimental results is developed using triaxial response parameters from experimental datasets. The performance of the proposed calibration method is demonstrated for three natural sands showing that it provides appropriate calibrated parameters for poorly graded sands with different relative densities and confined with varying effective stress magnitudes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107241"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808570","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":"Study on the torsion and sway responses of a massless rigid foundation excited by spatial obliquely incident P-waves considering the nonlinearity of the surrounding soil","authors":"Gang Xiong,&nbsp;Fangyuan Zhou,&nbsp;Bo Zhou","doi":"10.1016/j.compgeo.2025.107252","DOIUrl":"10.1016/j.compgeo.2025.107252","url":null,"abstract":"<div><div>In this study, the influences of the spatial incidence angles of P-waves and soil nonlinearity on the soil–structure interaction system are investigated. A finite element model of the SSI system is developed using the equivalent nodal force method, viscoelastic artificial boundary, and the revised Davidenkov model to analyze the dynamic responses of a rectangular foundation and its associated structure. The results indicate that the rectangular foundation predominantly exhibits sway responses to spatially oblique incidence P-waves, with a negligible torsion response. Compared with vertically incident P-waves, spatially oblique incident P-waves significantly increase the maximum sway and torsional angles of structural floors. Furthermore, soil nonlinearity increases the sway angle of the foundation and intensifies the maximum torsional and sway responses of the floors. Therefore, soil nonlinearity and the spatial incidence angles of P-waves need to be considered in structural seismic hazard assessments.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107252"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808571","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 analysis of soil arching and critical height in column-supported embankments using the Material Point Method","authors":"Wei He,&nbsp;Kaibo Yang,&nbsp;Pingbao Yin","doi":"10.1016/j.compgeo.2025.107229","DOIUrl":"10.1016/j.compgeo.2025.107229","url":null,"abstract":"<div><div>The arching effect and critical height in column-supported embankments were studied using the Material Point Method (MPM). First, the method was validated using model and field test results. Based on these validations, additional scenarios covering various embankment heights and column spacings were analyzed using MPM. The results showed that the shear pattern evolves from a punching shear pattern for embankments lower than 0.7(s-a) to an inverted general bearing pattern for embankments higher than 1.3(s-a). These observations were further validated through numerous 3D analyses based on the soil profile of a real project. Lastly, additional case studies from China, the 3D MPM analysis results, and previously available data, were compared with existing criteria, leading to the proposal of a modified formula for the critical height of embankment to consistently reduce differential settlement to under 0.1%.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107229"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808572","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 compressive rate-dependent hypoplastic constitutive model and numerical application in SPH for granular soil","authors":"Yaru Lv,&nbsp;Sichao Ding,&nbsp;Sitong Liu,&nbsp;Ding Chen,&nbsp;Yuchen Su,&nbsp;Wenxiong Huang","doi":"10.1016/j.compgeo.2025.107223","DOIUrl":"10.1016/j.compgeo.2025.107223","url":null,"abstract":"<div><div>The mechanical behavior of granular soil is significantly influenced by the strain rate to which it is subjected. However, there is currently no widely applicable constitutive model that effectively captures the rate-dependent behavior of granular soil across the spectrum from quasi-static to extreme dynamic loading conditions. This paper presents the findings from both physical and numerical experiments, demonstrating that particle strength displays three distinct regimes over a strain rate range of 10⁻⁶ to 10⁵ s⁻¹ resulting from the competition between coexisting thermally activated and macro-viscous mechanisms. Given that particle strength is associated with the granular hardness (<em>h</em>_s), this study proposes a compressive rate-dependent hypoplastic constitutive (CRDHC) model by establishing a rate-dependent <em>h</em>_s. Furthermore, the significantly elevated initial stiffness resulting from high strain rates and the concave limit stress surface induced by denser states are addressed through small strain stiffness and stress transformation, respectively. The proposed CRDHC model is validated through a one-dimensional quasi-static test and a Split Hopkinson Pressure Bar (SHPB) test, and it is implemented numerically using the Smoothed Particle Hydrodynamics (SPH) meshfree method. The method incorporates the substepping scheme with error control and stability correction at minimal stress states to establish a comprehensive framework for granular soil under low to high strain rates. Both discharge of silo and high-velocity penetration are simulated to validate the efficiency and accuracy of the framework in capturing the rate-dependent response in large deformation problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107223"},"PeriodicalIF":5.3,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143808576","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}