International Journal for Numerical and Analytical Methods in Geomechanics最新文献

{"title":"Semi‐Analytical Solution for a Lined Non‐Circular Tunnel in Viscoelastic Rock","authors":"Hongliang Liu, Xin Gao, Hui Cai, Yuxue Chen, Hongyun Fan","doi":"10.1002/nag.70323","DOIUrl":"https://doi.org/10.1002/nag.70323","url":null,"abstract":"Based on the complex variable method and the corresponding principle of viscoelasticity, viscoelastic solutions for the stress and displacement of a lined non‐circular tunnel subjected to in‐situ stresses and internal water pressure is derived. The basic equations for solving the analytic functions are established according to the stress boundary condition along the inner boundary of the lining and the stress and displacement continuity conditions along the rock‐lining interface. The analytic functions are expressed as Laurent series and the Laplace transformation is performed on the basic equation. Herein, the power series method is applied to obtain the linear equations which are expressed by the analytic function coefficients in the Laplace domain. The stress and displacement solutions of tunnel in Laplace domain can be addressed by solving the equations, and then the viscoelastic solutions are obtained through Laplace Inverse transformation. Subsequently, an example for the horseshoe‐shaped tunnel is performed. The example used the generalized Kelvin model to simulate the rheological properties of surrounding rock mass. The obtained solution is compared with the numerical solution. The influences of the lateral pressure coefficient and the internal water pressure on the stresses and displacements of lining are analyzed.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"47 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-04-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147681791","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Microscopic Particle Interaction Mechanisms in Fluidized Soil Inrush Revealed by DEM‐CFD","authors":"Xiangdong Meng, Wanghua Sui, Chang Zhou, Binglun Gao","doi":"10.1002/nag.70313","DOIUrl":"https://doi.org/10.1002/nag.70313","url":null,"abstract":"Soil inrush is a typical silo‐type geohazard. To investigate particle behavior during this process, DEM‐CFD numerical simulation was conducted to comprehensively analyze the evolution of velocity fields, stress fields, force chains, and coordination numbers. The results indicate that the arching force chains near the inrush point continuously break and re‐form at higher locations as the inrush progresses. The interior of the arch, characterized by rapid particle transport, a sharp stress decrease, and a low coordination number, indicates that the soil mass is highly loosened and experiences continuous particle loss. This loosened instability zone progressively develops upward. Once the instability zone reaches the ground surface, it rapidly expands across the surface, causing surrounding particles to converge toward the center. Based on the evolution of microscopic contact behavior and macroscopic flow characteristics, this study reveals the dynamic transition mechanism of the soil mass. Specifically, five zones with distinct particle‐behavior characteristics are identified during the early stage of soil inrush, and a sixth zone subsequently emerges near the ground surface as surface collapse develops. Finally, a soil state model based on the ellipsoid of motion was established and used to delineate the influence range of a documented inrush accident. This study systematically investigates particle behavior during soil inrush and provides guidance for inrush prevention and mitigation in concealed engineering projects.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"5 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147655937","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Development of Safe Cut‐Out Distance Framework for Continuous Miner Operations in Indian Underground Coal Mines: A Numerical Simulation and Field Studies","authors":"Tamilprasanth M, Sahendra Ram","doi":"10.1002/nag.70317","DOIUrl":"https://doi.org/10.1002/nag.70317","url":null,"abstract":"The continuous miner (CM)‐based coal extraction was introduced two decades ago in Indian coalfields. Safe cut‐out distance (COD), represents the maximum safe and stable span of a fixed‐width gallery excavated in a single pass using CM. The amount of coal that can be cut at once without support defines the successful implementation of CM. The existing guidelines exhibit considerable variation under similar conditions, making it challenging to determine a safe COD. In this research, a simulation study is carried out for 20 selected panels of different coalfields in India, considering strain‐softening material properties for the immediate roof to accurately capture the geo‐mechanical response. The simulation study found that Rock Mass Rating (RMR), Young's modulus, gallery width, and depth of cover are the most influential parameters in determining safe COD. Safe COD is defined considering a maximum displacement of 5 mm before the installation of roof bolts in the numerical modelling. Multivariate regression analysis is carried out on the basis of the measured safe COD in the models under different conditions, and a formulation is developed with an <jats:italic>R</jats:italic> <jats:sup>2</jats:sup> value of 0.86. The developed formulation is validated in the field through strata monitoring observation and found to predict the safe COD consistently compared to the other available design guidelines. However, this formulation applies only in the absence of weak geological disturbances and is unsuitable for areas affected by such conditions. This newly developed design guideline increase productivity and reduce the cost and time involved in designing the COD in the field.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"28 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147655891","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Investigation of the Influence of Interparticle Cohesion During the Shear Behavior of Granular Soil","authors":"Thao Doan, Buddhima Indraratna, Thanh T. Nguyen, Cholachat Rujikiatkamjorn","doi":"10.1002/nag.70319","DOIUrl":"https://doi.org/10.1002/nag.70319","url":null,"abstract":"Interparticle cohesive forces play a crucial role in governing the shear behavior of many soils; however, quantifying this effect remains challenging due to limited microscopic insights available from laboratory experiments. In this context, the current study aims to investigate the influence of microscale cohesion on the shear response using the discrete element method (DEM), where the relevant stress, strain, and dilation characteristics are compared with experimental data before conducting a detailed micromechanical analysis. The results successfully replicate important features of a moderately compacted soil under shearing such as strain softening and contraction‐dilation responses. Interestingly, cohesion‐induced samples exhibit a sustained increase in shear stress even at high shear strain levels; for instance, the residual shear stress in the highly cohesive sample increased by approximately 40% compared to the non‐cohesive sample. This occurs because particle‐scale cohesion enhances interparticle bonding, leading to higher average contact numbers (coordination number CN) and increased skeletal forces. A key novelty of this study lies in establishing a quantitative correlation between conventional macroscale shear parameters and microscale CN under varying loading and cohesion conditions. This proves that the more cohesive the particles, the stronger the bonding effect with higher CN, consequently resulting in higher shear resistance. Moreover, this trend is exacerbated when the rolling friction of particles is introduced as the coupled effect of microscale cohesion and rolling friction strengthens the soil's resistance to shear deformation.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"20 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147655892","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Soil Nail Design Optimization Through Adaptive Slope Stability Analysis Using Improved Limit Equilibrium Method","authors":"Weihang Ouyang,&nbsp;Kai Liu,&nbsp;Si-Wei Liu","doi":"10.1002/nag.70252","DOIUrl":"10.1002/nag.70252","url":null,"abstract":"<p>Design optimization for automatically generating optimal design results is a promising technique for enhancing the efficiency of design processes and outcomes. However, its development for soil nail reinforced slopes is limited since the traditional slope stability analysis using the limit equilibrium method (LEM) becomes relatively time-consuming when the LEM is repetitively employed during design optimization. In this article, an improved LEM for slip surfaces within reinforced slopes is developed, which eliminates the need for dense slice division by only aligning several sampling points using Gaussian quadrature. Unlike conventional design methods, wherein soil nails are represented by predetermined fixed point loads, the proposed improved LEM considers them through equivalent point loads adaptively updated according to the slip surface and nail positions to avoid overly conservative design results. Moreover, an adaptive slope stability analysis (ASSA) is proposed for reinforced slopes, offering an effective evaluation of stability without the need to step through all potential slip surfaces. The effectiveness of the proposed method is validated through several examples, demonstrating substantial computational cost savings of 98.5% compared to directly implementing traditional LEM-based analysis methods in design optimization, as well as material savings of 30% relative to results from manually designing.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2806-2822"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.70252","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146122037","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Plane Strain Consolidation of Saturated Soils in a Finite Rectangular Domain With Partially Pervious Vertical Boundaries","authors":"Kaijun Zhang,&nbsp;Linzhong Li,&nbsp;Xiao Yang,&nbsp;Manlin Liu,&nbsp;Yi Tian","doi":"10.1002/nag.70253","DOIUrl":"10.1002/nag.70253","url":null,"abstract":"<div>\u0000 \u0000 <p>In practical geotechnical engineering, boundary conditions often exhibit partially pervious characteristics. Previous research on the consolidation behavior of saturated media within a finite domain has generally overlooked the effect of boundary perviousness on the consolidation behavior. This study develops a modified analytical model to examine the plane strain consolidation of saturated porous media with a partially pervious vertical boundary. By incorporating the existing generalized solution with the partially pervious vertical boundary condition, the analytical derivation is first carried out in the transformed domain. The corresponding exact solution in the physical domain is subsequently obtained through a Fourier series expansion combined with Crump's numerical Laplace inversion technique. The reliability of the developed analytical solution is verified by comparing the present results with those in the literature under conditions of full perviousness. Using the proposed solution, numerical calculations are then performed to investigate the effects of boundary perviousness parameters and loading width on the consolidation behavior. Furthermore, the evolution of excess pore water pressure, settlement, and horizontal displacement over space and time is analyzed. The results demonstrate the significant influence of partially pervious boundaries on excess pore water pressure dissipation and soil displacement evolution.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2823-2834"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146122038","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Calibration and Stabilization Performance of Second-Generation Negative Poisson's Ratio Bolts in Deep Steeply Inclined Roadways","authors":"Yifei Gong,&nbsp;Ming Huang,&nbsp;Chun Zhu,&nbsp;Yi Tan,&nbsp;Murat Karakus,&nbsp;Zhigang Tao,&nbsp;Wei Sun","doi":"10.1002/nag.70229","DOIUrl":"10.1002/nag.70229","url":null,"abstract":"<div>\u0000 \u0000 <p>As coal mining in China advances to greater depths, deep and steeply inclined layered roadways are increasingly subjected to high-stress environments, leading to pronounced nonuniform deformation and instability. These failure behaviors pose significant challenges to roadwatability control and engineering safety. To address this problem, this study investigates the −1000 meters (m) steeply inclined roadway of the Qishan Mine in Xuzhou and aims to clarify the deformation evolution mechanism under deep high-stress conditions and to evaluate the applicability of second-generation negative Poisson's ratio (2G-NPR) bolts for large-deformation control. Building on previously conducted physical model tests, a 3DEC numerical model was established to reproduce the full deformation–failure process of the surrounding rock. A key contribution of this work is the calibration and numerical implementation of the mechanical characteristics of 2G-NPR bolts, including their high ductility and constant-resistance plateau behavior, enabling realistic representation of NPR reinforcement in discrete-element simulations. The calibrated model was then applied to assess the reinforcing effects of NPR bolts in deep layered rock masses. The results reveal a distinct asymmetric deformation pattern, with the roadway roof and right side identified as critical instability zones. Deformation and failure are dominated by slip and separation along the layered structural planes under high stress. Owing to their constant resistance and large elongation capacity, 2G-NPR bolts substantially suppress large deformations, reducing roadway sidewall displacement by more than 20% compared with traditional bolt support. This study provides new insights into the deformation mechanisms of deep, steeply inclined layered roadways and demonstrates the engineering advantages of 2G-NPR bolts, offering an effective reinforcement strategy for controlling large deformations in deep mining environments.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2771-2789"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146095767","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Simple Bounding Surface Plasticity Model for Overconsolidated Clays: Theory, Validation, and Numerical Implementation","authors":"Kehao Chen,&nbsp;Rui Pang,&nbsp;Bin Xu,&nbsp;Yang Zhou,&nbsp;Long Yu","doi":"10.1002/nag.70231","DOIUrl":"10.1002/nag.70231","url":null,"abstract":"<div>\u0000 \u0000 <p>Constitutive models with concise theoretical formulations and readily measurable parameters are vital for practical geotechnical engineering applications. This study presents a novel constitutive model for overconsolidated (OC) clays by integrating bounding surface plasticity theory into the modified Cam Clay (MCC) framework without introducing any additional model parameters. The proposed model enhances the MCC model in both strength and deformation predictions by (i) introducing the overconsolidation parameter into the dilatancy relation, enabling a more accurate representation of the shear dilatancy behavior of OC clays, and (ii) incorporating the Hvorslev envelope into the plastic modulus interpolation function to capture the strain-softening behavior and peak strength. The model's performance is validated through element-level simulations of compression and extension tests on clays, encompassing a broad range of overconsolidation ratios (OCRs) and stress paths. Additionally, the model is implemented in the ABAQUS finite element platform using an explicit integration scheme with automatic error control. Its practical applicability is demonstrated through the simulation of a centrifuge plate loading test on an OC clay foundation, with numerical results showing strong agreement with experimental data.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2675-2690"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146056197","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Semi-Analytical Study on Thermal Consolidation of Semi-Infinite Saturated Normally-Consolidated Soils Under Impeded Boundary","authors":"Lun Hua,&nbsp;Yue Gui,&nbsp;Yi Tian,&nbsp;Jinxiao Luo,&nbsp;Wenbing Wu,&nbsp;Mingxi Ou","doi":"10.1002/nag.70262","DOIUrl":"10.1002/nag.70262","url":null,"abstract":"<div>\u0000 \u0000 <p>This paper presents a semi-analytical investigation into the thermal consolidation behavior of saturated normally-consolidated (NC) clay within a two-dimensional semi-infinite domain, incorporating both thermal contraction and a realistic impeded boundary condition. A novel thermo-hydro-mechanical (THM) coupling model is established, building upon a constitutive relationship that decomposes thermal strain into three components to accurately capture the irreversible volumetric contraction of saturated NC clay under drained heating. The corresponding semi-analytical solutions for the temperature, excess pore water pressure (EPWP), and displacement under a strip-type thermomechanical load are derived using Laplace-Fourier transformation and validated against existing analytical results. Parametric studies reveal that heating induces complex EPWP distributions, including transient negative pore pressure zones. Soil deformation exhibits an initial heave followed by consolidation settlement, the magnitude of which is primarily controlled by the ultimate temperature increase. The surface drainage condition significantly influences the response, with poorly permeable boundaries promoting heave-dominant deformation. Furthermore, both the magnitude and spatial extent of the thermal load are shown to be critical factors governing the evolution of EPWP and the rate of heat transfer.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2886-2902"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146122034","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Investigation of Nonuniform Creep Behavior and Crest Cracking of an Asphalt–Concrete Core Rockfill Dam After Reservoir Impoundment","authors":"Wen He,&nbsp;Hangyu Mao,&nbsp;Sihong Liu,&nbsp;Liujiang Wang,&nbsp;Chaomin Shen","doi":"10.1002/nag.70251","DOIUrl":"10.1002/nag.70251","url":null,"abstract":"<div>\u0000 \u0000 <p>Variations in environmental humidity induce changes in the moisture content of rockfill materials, thereby influencing their creep behavior. In core rockfill dams, reservoir impoundment alters the environmental humidity conditions of the upstream rockfill materials, which often induces nonuniform creep deformation between the upstream and downstream zones. Excessive nonuniform creep deformation may lead to structural damage of the dam during its subsequent operation. This study presents an in-depth analysis of the impact of nonuniform creep deformation following reservoir impoundment on the long-term safety of an asphalt–concrete core dam exhibiting pre-existing crest cracks. To this end, a classical empirical creep model was modified to incorporate the humidity-dependent creep behavior of rockfill materials and subsequently integrated into a finite-element program for numerical analysis. Model parameters used in the modified creep model were identified using on-site monitoring data through the backpropagation-particle swarm optimization inversion method. The good agreement between the calculated results and existing monitoring data indicates the validity of the proposed numerical simulation scheme. Based on the simulation results, a thorough discussion is presented to clarify the causes of crest cracking and evaluate the safety of the dam. Stress analysis of the core crest pavement reveals that differential saturation levels on either side of the core wall induce tensile stresses in the pavement, which lead to cracking. Further stress analysis of the core wall suggests a low probability of cracking or hydraulic fracturing. Meanwhile, predictive deformation analyses indicate that creep deformation is expected to stabilize approximately 4 years after impoundment.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 6","pages":"2855-2870"},"PeriodicalIF":3.6,"publicationDate":"2026-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146122036","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}