Computers and Geotechnics最新文献

{"title":"Modelling the lateral behavior of trunk-root-soil systems with representative branching patterns using a combined p-y framework","authors":"Fan He ,&nbsp;Teng Liang ,&nbsp;Rui Zhao ,&nbsp;Tianfeng Zhu ,&nbsp;Yu Zhao ,&nbsp;Anthony Kwan Leung ,&nbsp;Liangtong Zhan","doi":"10.1016/j.compgeo.2025.107654","DOIUrl":"10.1016/j.compgeo.2025.107654","url":null,"abstract":"<div><div>The Beam-on-Non-linear-Winkler-foundation (BNWF) model using <em>p</em>-<em>y</em> formulation is widely used on laterally loaded structures (e.g., piles) due to its computational efficiency. However, its applicability to branched structural systems (e.g., trunk-root-soil systems) remains highly uncertain given their intricate branching geometries and proclivity for substantial deflections, which contrast with conventional pile structures. This work tries to address this issue through developing a combined <em>p</em>-<em>y</em> framework. Initially, physical model tests were conducted using four 3D-printed representative branched root architectures—tap, herringbone, dichotomous, and T-shape patterns—under dry and saturated sand conditions, thereby providing a comprehensive dataset for numerical model development. Subsequently, the <em>p</em>-<em>y</em> framework was established by considering the large deflection response of parent roots and the push-in/uplift behavior of branch roots during tree overturning. For modeling large deflection response, two types of <em>p</em>-<em>y</em> formulations were critically evaluated: the commonly used Reese <em>p</em>-<em>y</em> curve within previous studies and the hyperbolic <em>p</em>-<em>y</em> curve. For push-in/uplift behavior, the deep flow mechanism and a bilinear softening trend were considered, respectively. The proposed <em>p-y</em> framework was integrated into an efficient hybrid beam-FE model to predict the lateral response of the trunk-root-soil systems. Model validation demonstrated that the Reese <em>p</em>-<em>y</em> curve inadequately captures large deflection responses due to its inherent assumptions. Conversely, the hyperbolic <em>p</em>-<em>y</em> curve effectively predicts lateral resistance and root behavior in branched root systems. Compared to other analytical models, the proposed model enhances the capability to predict the overturning behavior of tree structures and enables the modeling of complex 3D root system geometries and root material variability, offering a potentially effective method for evaluating tree stability in practical applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107654"},"PeriodicalIF":6.2,"publicationDate":"2025-09-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145109585","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":"High-Fidelity Morphological and Mineralogical Reconstruction of Chang’E-5 Lunar Regolith Grains","authors":"Quan Zheng ,&nbsp;Ziwei Tian ,&nbsp;Songzheng Yu ,&nbsp;Ronghua Pang ,&nbsp;Guang Zhang ,&nbsp;Guanghui Liu ,&nbsp;Yiwei Liu ,&nbsp;Yang Li ,&nbsp;Xin Liu ,&nbsp;Shijing He ,&nbsp;Ran Niu ,&nbsp;Peng Zhang","doi":"10.1016/j.compgeo.2025.107661","DOIUrl":"10.1016/j.compgeo.2025.107661","url":null,"abstract":"<div><div>The complex structure and physicochemical properties of lunar soil particles form a critical foundation for future lunar exploration and in-situ resource utilization. To address the distortion introduced by oversimplified particle morphology and composition assumptions in conventional modeling approaches, we propose a high-fidelity morphological and mineralogical reconstruction method of Chang’E-5 lunar regolith grains using micro-CT characterization and advanced image processing. Twenty-one representative lunar particles were selected and subjected to high-resolution X-ray computed tomography (X-CT) scanning for comprehensive characterization of their 3D internal structures. A segmentation workflow was developed to achieve precise particle identification and mesh model generation, integrating unsupervised clustering, morphological optimization, and voxel-to-mesh conversion. Scanning electron microscopy and energy-dispersive spectroscopy were employed to identify mineral phases and reconstruct their spatial distribution within the particles. Based on these characterizations, models applicable to the discrete element method and finite element method were constructed. This work forms a comprehensive digital lunar particle dataset encompassing optical images, X-CT 3D images, mesh models, discrete element and finite element models. This integrated approach enables high-fidelity representation of morphology, internal structure, and mineralogical composition of particles, significantly enhancing modeling accuracy and simulation scalability compared to previous methods. The digital model repository established in this study provides standardized and extensible data resources to support mechanical simulations, in-situ resource evaluation, and the design of lunar detection devices.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107661"},"PeriodicalIF":6.2,"publicationDate":"2025-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107014","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 modeling of columnar jointed basalt with a dual irregular joint network and its unloading response under true triaxial conditions","authors":"Xiangcheng Que ,&nbsp;Shu Zhu ,&nbsp;Zhende Zhu ,&nbsp;Jian Zhang ,&nbsp;Zihao Niu","doi":"10.1016/j.compgeo.2025.107660","DOIUrl":"10.1016/j.compgeo.2025.107660","url":null,"abstract":"<div><div>The irregular columnar structure and internal hidden joints of Baihetan columnar jointed basalt (CJB) lead to its prominent anisotropy, and developing an appropriate model and accurately understanding its unloading responses are crucial for engineering safety. Based on the natural structure of Baihetan CJB, a dual irregular joint network numerical model, including both columnar joints and internal hidden joints, was proposed to investigate its anisotropic mechanical characteristics and damage failure mechanisms under true triaxial loading and unloading conditions. The results demonstrate that the dual irregular joint model with an irregularity coefficient of 0.47 effectively captures the cross-sectional geometry and internal structure features of the actual CJB. After a two-stage parameter calibration, the numerical and laboratory uniaxial test results are relatively consistent, and cracks exhibit four typical propagation modes. Under true triaxial loading and unloading conditions, both the deformation and strength behaviors of CJB show significant anisotropy. The unloading of minimum principal stress weakens the mechanical properties of CJB and increases its anisotropic level. The failure of the model with 15° inclination angle under unloading condition helps explain the column fragmentation observed at the side walls after excavation, which aligns with stress-structure controlled failure characteristics. The progressive damage of CJB involves three stages, and unloading the minimum principal stress or reducing the restriction of intermediate principal stress will accelerate the internal damage rate. The research results provide valuable guidance for the rational design and safe construction of related projects.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107660"},"PeriodicalIF":6.2,"publicationDate":"2025-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107011","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":"Comparative analysis of helical pile uplift behaviour in sand using three different constitutive models","authors":"Mohsen Misaghian,&nbsp;Pan Hu,&nbsp;Chin Leo,&nbsp;Samanthika Liyanapathirana","doi":"10.1016/j.compgeo.2025.107651","DOIUrl":"10.1016/j.compgeo.2025.107651","url":null,"abstract":"<div><div>This research conducted a comparative study to compare the performance of three constitutive models, i.e., the Mohr-Coulomb model, the Hypoplasticity with intergranular strain model, and the Sanisand model, in assessing the behaviour of helical piles subjected to uplift loading in dry dense sand. Triaxial tests were initially simulated using three models to calibrate their parameters and examine their predictive capabilities. Subsequently, finite element simulations incorporating the calibrated models were developed to replicate a comprehensive series of centrifuge uplift tests. These tests include single-helix helical piles under both shallow and deep embedment conditions, as well as multi-helix helical piles with varying helix spacing ratios. The performance of each model was evaluated based on element-level response, load–displacement behaviour, failure mechanisms, void ratio evolution, and uplift capacity predictions. Comparative analysis of the numerical results against centrifuge test data revealed that the advanced constitutive models, i.e., Hypoplasticity and Sanisand models, provided a more realistic representation of helical pile behaviour under uplift loading conditions. This research concludes that these advanced models offer higher reliability for designing and analysing helical piles, which are increasingly recognised as an efficient foundation solution for offshore structures.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107651"},"PeriodicalIF":6.2,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107012","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 three-dimensional grouting model of slurry diffusion-particle deposition-permeability evolution in discrete fracture network","authors":"Jie Hu,&nbsp;Gao Liang Liu,&nbsp;Jia Qing Chen,&nbsp;Yun Min Chen","doi":"10.1016/j.compgeo.2025.107658","DOIUrl":"10.1016/j.compgeo.2025.107658","url":null,"abstract":"<div><div>Accurately characterizing slurry behavior in discrete fracture network remains challenging due to the complexity of fracture geometries, the deformability of apertures, and the nonlinear rheology of slurry. This study develops a novel three-dimensional grouting model that integrates slurry diffusion, particle deposition, and permeability evolution in discrete fracture network. The evolving composition of the slurry is described by the volume fractions of deposited (<em>S</em>) and suspended (<em>C</em>) particles. A spatial–temporal viscosity field is incorporated to capture the time-dependent changes in slurry properties as it diffuses along the flow paths. Moreover, fracture apertures deform dynamically under the influence of slurry pressure and particle deposition, establishing a bidirectional feedback mechanism between the flow field and fracture geometry. The accuracy of the numerical model is validated by comparing the predicted and experimental results of permeability evolution in both single fracture and fracture network. The results show that permeability decreases hysteresis arises from delayed particle deposition, while increased fracture intersections accelerate the transformation from suspended to deposited particles. Sensitivity analyses reveal that slurry diffusion is primarily governed by slurry properties such as grouting pressure and water–cement ratio, whereas deposition and sealing behavior are strongly influenced by fracture aperture and network density. Finally, the model facilitates the determination of optimal grouting duration under various conditions, offering a reliable reference for designing effective grouting strategies in complex fracture networks.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107658"},"PeriodicalIF":6.2,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107013","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":"Dimensionless physics-informed neural networks for remote-excitation SH-wave propagation in truncated domains","authors":"Hamid Taghavi Ganji,&nbsp;Elnaz Seylabi","doi":"10.1016/j.compgeo.2025.107612","DOIUrl":"10.1016/j.compgeo.2025.107612","url":null,"abstract":"<div><div>Physics-informed neural networks (PINNs) are explored for the first time as a mesh-free solver for two-dimensional shear–horizontal (SH) wave propagation in truncated domains excited by remote seismic sources. Building on the Domain Reduction Method (DRM), we derive a jump (substructuring) condition that transfers free-field displacements and tractions onto an interior interface and combine it with Lysmer–Kuhlemeyer (LK) absorbing and transmitting boundary conditions. A dimensionless scaling, extended here from one to two dimensions, removes heuristic loss-weighting and markedly improves training robustness. Three truncation configurations of increasing sophistication are benchmarked against finite-element (FEM) solutions for (i) a synthetic toy model, (ii) a field-scale homogeneous basin, and (iii) a heterogeneous basin with depth-varying shear-wave velocity. Across all cases, the dimensionless PINN reconstructs displacement fields with normalized errors below 9%, while the most effective configuration — separating the jump interface from the absorbing boundary by half a characteristic wavelength — reduces FEM truncation error to <span><math><mo>&lt;</mo></math></span> 4% and total PINN error to <span><math><mo>&lt;</mo></math></span> 8%. The proposed framework, therefore, offers an accurate and flexible alternative for large-scale site response and soil–structure interaction analyses subjected to remote excitation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107612"},"PeriodicalIF":6.2,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107009","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":"Dynamic stability of piles with fractional damping foundation","authors":"Mohammadmehdi Shahroudi,&nbsp;Yanglin Gong,&nbsp;Jian Deng","doi":"10.1016/j.compgeo.2025.107653","DOIUrl":"10.1016/j.compgeo.2025.107653","url":null,"abstract":"<div><div>This paper investigates the dynamic stability and vibration response of a slender pile embedded in soil, both theoretically and numerically. The pile is modeled as a column resting on a Winkler-type elastic foundation with fractional damping characteristics. Investigation of the pile loaded axially leads to a fractional Mathieu differential equation of motion. To analyze stability, the Bolotin method is employed to derive approximate instability boundaries. Additionally, a numerical approach based on block-pulse functions is developed to construct detailed instability diagrams, which also serve to validate the theoretical results obtained using the Bolotin method. It is found that higher order instability regions generally need higher order infinite determinants. Finally, a numerical case study is presented to explore the effects of key parameters—such as the fractional order, foundation stiffness, damping coefficient, and static and dynamic load components—on the pile’s dynamic stability. The dynamic stability of a pile subjected to a real earthquake is determined.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107653"},"PeriodicalIF":6.2,"publicationDate":"2025-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107008","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":"An extended quasi-bond method with isoparametric implementation for fracture problems in rock-like materials","authors":"Ren-Jie Zhang ,&nbsp;Qi-Zhi Zhu ,&nbsp;Wei-Jian Li ,&nbsp;Xing-Guang Zhao","doi":"10.1016/j.compgeo.2025.107611","DOIUrl":"10.1016/j.compgeo.2025.107611","url":null,"abstract":"<div><div>Efficient and accurate prediction of damage and fracture remains one of main challenges in solid mechanics and computational mechanics. While the original quasi-bond method (QBM) has demonstrated its applicability in simulating solid deformation and characterizing crack nucleation and propagation, it remains the issue of fixed Poisson’s ratio and its numerical implementation still relies on computationally intensive calculation of geometric intersection points between quasi-bonds and element edges (sides) and lacks a standardized method for assembling the system’s stiffness matrix. In this paper, the effects of shear deformation are first incorporated into the QBM theory in order to remove the limitation of Poisson’s ratio. An isoparametric quasi-bond method is developed by introducing an integral transformation technique, which establishes a unified approach of constructing a system of quasi-bonds for both 2D and 3D problems. When comparing with the conventional operation, the present technique allows reducing computational cost of the initialization process while achieving higher numerical accuracy with a lower bond density. The proposed QBM exhibits lower matrix bandwidth compared to peridynamics approach, and eliminates the need for calculating continuous damage field as required in phase-field method, thereby combining computational efficiency with simplicity. Furthermore, the isoparametric formulation leads to an efficient and natural QB–FE coupling, which provides a seamless interface for integration into commercial finite element programs. Various validations through typical benchmark tests confirm computational efficiency of the proposed approach in solving both elasticity and elastic fracture problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107611"},"PeriodicalIF":6.2,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145106988","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":"Modelling cliff collapse and run-out with the material point method","authors":"Sam J.V. Sutcliffe ,&nbsp;William M. Coombs ,&nbsp;Wangcheng Zhang ,&nbsp;Ravindra Duddu ,&nbsp;Charles E. Augarde","doi":"10.1016/j.compgeo.2025.107608","DOIUrl":"10.1016/j.compgeo.2025.107608","url":null,"abstract":"<div><div>This paper proposes a method for modelling the combined collapse and large deformation run-out of cliffs using the Material Point Method (MPM). An explicit dynamic MPM formulation is derived with features key for simulating ductile–brittle fracture dynamics over varying timescales. An isotropic non-local plastic-damage scheme is presented that may be used to simulate crack initiation and growth leading to brittle collapse of a cliff. The scheme is applicable to both quasi-static and dynamic cases. The behaviour of the proposed plastic-damage model is compared to experimental data utilising numerical shear-box tests. A numerical case study is carried out on a real-world chalk cliff to show the model’s applicability to rock-slope mechanics in both failure initiation and post-failure run-out.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107608"},"PeriodicalIF":6.2,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145106989","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":"Fractional derivative viscoelastic models based on L2-1σ formula: modelling and numerical application","authors":"Pan Ding ,&nbsp;Riqing Xu ,&nbsp;Lixing Wang ,&nbsp;Zhiran Gao ,&nbsp;Xiaonan Ge ,&nbsp;Minjie Wen","doi":"10.1016/j.compgeo.2025.107633","DOIUrl":"10.1016/j.compgeo.2025.107633","url":null,"abstract":"<div><div>This study develops two three-dimensional fractional derivative viscoelastic models and introduces a novel <em>L</em>_2-1σ difference formula for numerical implementation. In deriving the fractional derivative models, a spring-pot element, defined by Caputo fractional calculus, replaces the dashpot or spring in the classical Kelvin-Voigt model. Using the Laplace transform approach, the creep compliance and stress relaxation modulus of the fractional derivative model are derived. Additionally, the models are extended into a generalized three-dimensional form, and a novel <em>L</em>_2-1σ difference algorithm is employed to discretize the constitutive equations. This method has been shown to provide superior computational efficiency and accuracy compared to the GL algorithm, with a convergence order of 3-<em>α</em>. Subsequently, one of the models was compiled and implemented through the UMAT interface of ABAQUS, facilitating both element-level and large-scale project simulations. The numerical predictions from the model were compared with the results of triaxial creep tests conducted on Hangzhou clay and Hong Kong marine clay. The results confirm a reasonable degree of agreement between the predicted and experimental results, validating the accuracy and feasibility of both the model and the <em>L</em>_2-1σ difference algorithm. Furthermore, when applied to a practical case involving the long-term deformation of an artificial island construction project, the predictions from this model exhibit strong consistency with the measured data. This work provides valuable insights into understanding and predicting the time-dependent behavior of viscoelastic materials.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"189 ","pages":"Article 107633"},"PeriodicalIF":6.2,"publicationDate":"2025-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145107010","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}