Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-16DOI: 10.1016/j.compgeo.2026.107914
Gerardo Grelle
{"title":"A Dimensionless Geomechanical Core to Inform ML-AI Prediction of Seismic-Induced Landslides","authors":"Gerardo Grelle","doi":"10.1016/j.compgeo.2026.107914","DOIUrl":"10.1016/j.compgeo.2026.107914","url":null,"abstract":"<div><div>The Geomechanical Core (GMC) is a novel, scalable, and modular computational framework designed as the central engine component of hybrid AI systems for predicting and mapping seismic-induced landslides. GMC is founded on three pillars: a generalized failure model, a dimensionless constitutive model, and a serial modular architecture. The generalized failure model is a physically consistent adaptive function that distributes the initial shear strength over a normalized logarithmic shear-surface template, accounting for potential strength reductions caused by pre-seismic straining. At the heart of the engine, the constitutive model reproduces nonlinear elastoplastic stress–strain behavior through the relationship between normalized seismic load and straining (ductility ratio). This dimensionless formulation integrates the Hardening Soil (HS) approach for slope deformation and the Limit Equilibrium (LE) method for triggering and sliding, enabling the computation of both reversible and irreversible excess pore-water pressures under seismic excitation. The GMC chain of modules can be configured to represent different sectors of a landslide mass, where interacting boxes exchange static stress for initialization and dynamically update deformation states during shaking. GMC demonstrates high computational efficiency and aims to bridge the gap between large-scale predictive modeling and advanced constitutive-based numerical simulations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107914"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145976133","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-07DOI: 10.1016/j.compgeo.2025.107887
Chen-Xi Tong , Xiao-Yi Qi , Hong-Wei Liu , Gang-Hai Huang , Hai-Chao Li , Sheng Zhang
{"title":"A rapid packing method for Sphere Discontinuous Deformation Analysis (SDDA) simulation of confined compression of crushable ballast","authors":"Chen-Xi Tong , Xiao-Yi Qi , Hong-Wei Liu , Gang-Hai Huang , Hai-Chao Li , Sheng Zhang","doi":"10.1016/j.compgeo.2025.107887","DOIUrl":"10.1016/j.compgeo.2025.107887","url":null,"abstract":"<div><div>The computational efficiency of modelling irregular particles has always been a challenge, particularly when particle breakage needs to be considered. This paper first proposes a novel particle packing method for ballast samples based on the Monte Carlo, dichotomy, and random iteration methods for reducing particle contact, minimizing the gravitational potential energy, and optimizing the spatial position of ballast particles. The packing method is then applied to one-dimensional compression tests of ballast samples within the Sphere Discontinuous Deformation Analysis (SDDA) framework. The differences in the crushing behaviour of ballast particles with varying shapes under various contact conditions are investigated. The results indicate a distinct three-stage crushing process for the ballast samples. In the initial compaction stage, particle rearrangement and corner damage without penetrating fracture are observed. A sharp force increases with significant fluctuations, and penetrating fractures, including primary, secondary, and tertiary types, are observed during the rapid development stage. The sample reaches a structurally stable state due to the cessation of particle fracture. Furthermore, the crushing behaviour of ballast particles is strongly influenced by inter-particle interactions, leading to characteristics that differ from single-particle simulations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107887"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924207","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-06DOI: 10.1016/j.compgeo.2025.107881
Run Shi, Jian Zhao
{"title":"Unveiling rock heterogeneity-driven multiphysics field coupling mechanisms under microwave radiation by 3D numerical study","authors":"Run Shi, Jian Zhao","doi":"10.1016/j.compgeo.2025.107881","DOIUrl":"10.1016/j.compgeo.2025.107881","url":null,"abstract":"<div><div>Microwave-assisted hard rock breakage technique can promote the excavation efficiency in a sustainable way. While beneficial, its complex and unclear multiphysics coupling mechanism limits its application and precise control. We qualitatively and quantitively analyse the involving Electromagnetic(E)-Thermal(T)-Mechanical(M) coupling mechanism from macro to <em>meso</em>-perspectives via a novel 3D high-fidelity heterogeneous rock numerical study. The results show that rock heterogeneity has a significant disturbing effect on the multiphysics field, with the maximum disturbance degree reaching 345.92 % in electromagnetic loss density, 12.48 % in temperature, and 104.96 % in first principal strain. Furthermore, microwave power intensity can amplify the disturbance effect of rock heterogeneity, as quantified by the Disturbance Amplification Factor (DAF), which under high intensity 6 kW microwave radiation for 1 min amplified the effect 6 times of electromagnetic field, 3.22 ∼ 4.44 times of thermal field, 1.98 ∼ 2.31 times of mechanical field, compared with the 1 kW for 6 min condition. This amplification effect decays along the ETM pathway, where microwave power nonlinearly intensifies multifield disturbances initially driven by mineral dielectric disparity. The multifield presents a 3D high core area coupling relation from the macroscopic level, and the spatio-temporal anchoring effect of mineral components under microwave radiation was dissected through time evolution analysis from the mesoscopic level. This approach reveals the multifield coupling laws and highlights the rock heterogeneity effects, which provide a foundation for target control strategies for future microwave rock breakage.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107881"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924555","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-06DOI: 10.1016/j.compgeo.2025.107884
Shuqing Wang, Hong Zheng
{"title":"Numerical manifold method for the transient HMC fully coupled model in triple-layer composite liners","authors":"Shuqing Wang, Hong Zheng","doi":"10.1016/j.compgeo.2025.107884","DOIUrl":"10.1016/j.compgeo.2025.107884","url":null,"abstract":"<div><div>A triple-layer liner composed of a geomembrane (GMB), a geosynthetic clay liner (GCL), and a compacted clay liner (CCL) is widely adopted in landfill liner systems to impede leachate migration, resulting in a typically inhomogeneous problem. In this study, a fully coupled three-field (<strong><em>u</em></strong><em>-p-c</em>) formulation is established for GMB/GCL/CCL composite liners. Then, a numerical manifold method (NMM) is developed to address the drawbacks of the finite element method in weakly discontinuous porous media, including the need for interface-fitted meshes, the difficulty in capturing derivative discontinuities, and the element distortion caused by cover cutting. Nonetheless, this also raises the issue of how to enforce the interface continuity of approximations to displacement (<strong><em>u</em></strong>), pore pressure (<em>p</em>), and pollutant concentration (<em>c</em>). Unlike the situation where only a physical field is involved, in which interface continuity can be approximated by either the penalty method or Lagrange multiplier methods, it is difficult to find a proper set of penalties for the coupled variables <strong><em>u</em></strong>, <em>p</em> and <em>c</em>. In this study, this difficulty is overcome by constructing the approximations of <strong><em>u</em></strong>, <em>p</em> and <em>c</em> such that they exactly satisfy the interface continuity conditions. As a result, it avoids the cover-cutting induced by interfaces as well as the need for penalty parameters or Lagrange multipliers. Numerical examples demonstrate that the proposed method can accurately simulate solute migration in composite liners and reveal solute migration characteristics, providing valuable guidance for the improvement of liner materials.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107884"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924597","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2025-12-31DOI: 10.1016/j.compgeo.2025.107875
Ye Yang , Mincai Jia , Qun Yan
{"title":"Transient responses of multilayered transversely isotropic half-space subjected to a rectangular oblique impact loading","authors":"Ye Yang , Mincai Jia , Qun Yan","doi":"10.1016/j.compgeo.2025.107875","DOIUrl":"10.1016/j.compgeo.2025.107875","url":null,"abstract":"<div><div>To address the oblique impact effect in engineering, this paper develops a theoretical solution for non-axisymmetric transient problems based on a three-dimensional Cartesian coordinate system. By combining the dynamic differential equations with displacement potential functions, the expressions for the displacement responses of single-layered soil are derived via the double Fourier-Laplace transformation. The global stiffness matrix for a multilayered half-space is assembled from the analytically derived layer-element stiffness matrices. Transient displacement and stress responses in the time domain are then obtained through inverse transform algorithms. The accuracy of the proposed method is validated by comparing the solutions with the calculation results from finite difference method. Furthermore, the effects of transversely isotropic parameters, loading embedment depth, and impact inclination angle on the dynamic response of the soil are studied using the presented approach.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107875"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883639","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2025-12-31DOI: 10.1016/j.compgeo.2025.107877
Lei Niu , Liyan Zhang , Shuchen Cai , Kang He , Xin Chen
{"title":"A Eulerian multiphase model for collapse and segregation of bidisperse granular columns","authors":"Lei Niu , Liyan Zhang , Shuchen Cai , Kang He , Xin Chen","doi":"10.1016/j.compgeo.2025.107877","DOIUrl":"10.1016/j.compgeo.2025.107877","url":null,"abstract":"<div><div>Granular collapse (landslide) is generally destructive to marine infrastructure and ecosystems. Although particle size distribution exerts a strong influence on collapse behavior, the specific mechanisms through which polydispersity affects granular collapse are not yet fully understood. In this study, a Eulerian multiphase model is developed to simulate the collapse of the bidisperse granular column. A bidisperse constitutive law that combines kinetic theory of granular flow (KTGF) and frictional dilation/contraction rheology is used to close the granular stresses. The KTGF incorporates anisotropic particle collisions and interstitial fluid drag, and the adopted drag model remains effective across dilute to dense particle volume fractions. The fluid turbulence model further accounts for density stratification arising from non-uniformity in the particle volume fraction. Both monodisperse and bidisperse granular collapse experiments were performed, and the predicted flow fields, runout distances, and segregation patterns show good agreement with the experimental observations. The results indicate that a higher proportion of coarse particles leads to a greater total runout distance and more pronounced segregation during bidisperse column collapse. Moreover, the choice of frictional viscosity model significantly influences the morphological evolution, underscoring the importance of accurately representing dilation and contraction effects. Overall, this study develops a multiphase framework for analyzing bidisperse granular collapse, providing valuable insights for hazard assessment and mitigation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107877"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883611","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-20DOI: 10.1016/j.compgeo.2026.107916
Yongzheng Zhang , Jidong Zhao , Gang Wang , Xuxu Yang , Huilong Ren , Timon Rabczuk , Yunpeng Zhao
{"title":"Dual-horizon peridynamic modeling of thermally induced fracture in anisotropic materials: A variational energy-based approach","authors":"Yongzheng Zhang , Jidong Zhao , Gang Wang , Xuxu Yang , Huilong Ren , Timon Rabczuk , Yunpeng Zhao","doi":"10.1016/j.compgeo.2026.107916","DOIUrl":"10.1016/j.compgeo.2026.107916","url":null,"abstract":"<div><div>Fracturing phenomena driven by temperature variations present substantial challenges across geotechnical engineering applications. Accurate computational representation of such behavior demands robust numerical architectures that can reliably capture discontinuous crack evolution within materials characterized by pronounced directional dependencies. This paper introduces a novel approach that integrates dual-horizon non-ordinary state-based peridynamics (DH-NOSBPD) with a variational damage model for simulating thermo-mechanical fracture in anisotropic media. The proposed framework overcomes the limitations of conventional peridynamic methods in representing anisotropic thermal and mechanical coupling while eliminating numerical instabilities inherent in bond-breaking criteria. A staggered coupling strategy is employed to synchronize thermal and mechanical field updates, incorporating anisotropic constitutive relationships for both heat conduction and stress–strain behavior. The variational damage model introduces a history-dependent scalar damage field derived from strain energy density, thereby circumventing the spurious energy release and mesh dependence associated with abrupt bond deletion. This approach yields physically consistent crack evolution. Numerical examples validate the framework’s accuracy in anisotropic heat transfer, mechanical deformation, and complex fracture patterns under combined thermo-mechanical loading. The model demonstrates superior stability and predictive capability compared to conventional bond-breaking approaches.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107916"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022678","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-20DOI: 10.1016/j.compgeo.2026.107935
Yecheng Qian, Lu Jing
{"title":"A GPU accelerated LBM-DEM framework for non-spherical convex particles and its application in permeability tests","authors":"Yecheng Qian, Lu Jing","doi":"10.1016/j.compgeo.2026.107935","DOIUrl":"10.1016/j.compgeo.2026.107935","url":null,"abstract":"<div><div>Particle shape plays an important role in particle–fluid systems, but accurate numerical tools for simulating the interactions between non-spherical particles and the fluid are still limited. One of the challenges is the efficient and accurate evaluation of the solid volume fraction in fluid nodes for non-spherical particles. Here, we develop a GPU accelerated computational framework that couples the lattice Boltzmann method (LBM) for the fluid phase with the spheropolyhedra-based discrete element method (DEM) for particles. Using the immersed moving boundary method to resolve fluid particle interactions, we introduce a novel sub-grid and surface-normal approach to precisely calculate the solid volume fraction for each fluid node. Compared to the traditional approximate polyhedron method in LBM-DEM, this approach achieves higher accuracy while maintaining computational efficiency. The improved method is validated by volume calculation of a given geometry and drag coefficient of individual Platonic polyhedral particles under varying Reynolds numbers. Furthermore, the fully coupled LBM–IMB–spheropolyhedra DEM framework is validated through a single free settling non-spherical particle case in a viscous fluid. This framework is then applied to predict the permeability of random packings of various non-spherical particles, and the predicted values are found to be consistent with the classic Kozeny–Carman and Ergun equations. Finally, to provide a more general prediction model accounting for the particle shape, the Ergun equation is modified with the sphericity of non-spherical particles, which shows good agreement with our results and the data from the literature.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107935"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022679","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-20DOI: 10.1016/j.compgeo.2026.107932
Chunlei Li , Long Yu , Yu Xin , Yunrui Han
{"title":"Quantifying the scale effects on shallow foundation bearing capacity induced by shear band on sensitive marine clays","authors":"Chunlei Li , Long Yu , Yu Xin , Yunrui Han","doi":"10.1016/j.compgeo.2026.107932","DOIUrl":"10.1016/j.compgeo.2026.107932","url":null,"abstract":"<div><div>Mapping from small-scale model tests to prototype design remains challenging in geotechnics, primarily due to scale effects induced by progressive failure in sensitive marine clays. While progressive failure is often negligible in small-scale models due to uniform mobilized strain within the failure soil, it becomes critical at the prototype scale. This study investigates these scale effects using a nonlocal Remeshing and Interpolation Technique with Small Strain (RITSS) method, incorporating a strain-softening constitutive model based on the simple modified elastic-perfectly plastic Tresca criterion. Simulations of plan-strain biaxial tests and shallow foundation penetration successfully capture shear band development. The results reveal that the extent of the strain-softening zone depends on the ratio of the model size to the intrinsic internal length of the clay, approximately the shear band thickness. A clear transition from localized to diffuse softening is observed as the model size decreases. Based on this finding, novel rescaling coefficients are proposed, incorporating the ratio of internal length to footing size and a clay ductility parameter. These coefficients enable more reliable prediction of the ultimate bearing capacity when mapping from reduced-scale models to prototypes in strain-softening clays.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107932"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022680","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}
Computers and GeotechnicsPub Date : 2026-04-01Epub Date: 2026-01-05DOI: 10.1016/j.compgeo.2025.107893
Zhijian Qiu , Qiwei Jin , Qingwei Wu , Muhammad Zayed , Ahmed Ebeido , Yewei Zheng
{"title":"An explainable resilience-informed framework for surrogate modeling and multi-objective optimization of embankments under seismic loading","authors":"Zhijian Qiu , Qiwei Jin , Qingwei Wu , Muhammad Zayed , Ahmed Ebeido , Yewei Zheng","doi":"10.1016/j.compgeo.2025.107893","DOIUrl":"10.1016/j.compgeo.2025.107893","url":null,"abstract":"<div><div>Ensuring seismic resilience of transportation earth structures is critical for maintaining lifeline functionality following major earthquakes, with interpretability and data-driven modeling being key to achieving maximum resilience and rapid demand prediction. This study presents an explainable resilience-informed framework that integrates finite element (FE) simulations, probabilistic demand modeling, and machine learning-based surrogate modeling to rapidly assess and optimize the seismic resilience of earth embankments. Within the framework, FE results from 1,000 embankment configurations subjected to 100 recorded ground motions are employed to train an explainable XGBoost model that accurately captures the nonlinear effects of key parameters on seismic response. In this regard, the trained surrogate model further facilitates efficient derivation of seismic fragility and resilience curves, quantifying both performance degradation and post-earthquake recovery. Consequently, a resilience-informed multi-objective optimization is performed to identify optimal geometric configurations of earth embankments that maximize seismic resilience while minimizing both the lateral and vertical deformations. Representative scenarios show that reducing embankment height and flattening the embankment slope significantly enhance seismic resilience for low-strength materials, minimizing the need for costly reinforcement or material enhancement. Overall, the developed explainable framework provides a transparent, data-driven, and physics-consistent approach for rapid prediction and optimization of equivalent resilient transportation earth structures.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107893"},"PeriodicalIF":6.2,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145924606","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}