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

{"title":"Study on Key Parameters of Automatic Roadway Formation Without Pillar by Roof Cutting Based on Floor Failure","authors":"Xingxing Zhang,&nbsp;Yubing Gao,&nbsp;Huacheng Sun,&nbsp;Qingli Liu,&nbsp;Chuanjiu Zhang,&nbsp;Qiukai Gai","doi":"10.1002/nag.70414","DOIUrl":"10.1002/nag.70414","url":null,"abstract":"<div>\u0000 \u0000 <p>Floor water inrush threatens the safe mining of the mine. The automatic roadway formation without pillar by roof cutting (ARFPRC) method can reduce the floor water inrush risk. However, the key parameters for mining above confined aquifers are difficult to determine. In response, this study proposes a design method for determining the key parameters of ARFPRC based on floor failure. First, the unique advantages of ARFPRC in reducing floor water inrush risk is elucidated. Then, numerical simulation studies are conducted to analyze the characteristics of floor damage under different roof cutting parameters, leading to the determination of the optimal parameters. It is also clarified that grouting reinforcement to improve rock properties is beneficial for preventing water inrush. Subsequently, a field test of floor grouting reinforcement was conducted on site. After its effectiveness was confirmed through geophysical detection, the ARFPRC industrial test was carried out. Finally, based on the measured depth of floor failure, it was demonstrated that the designed parameters effectively controlled floor damage. The research findings provide a theoretical foundation and engineering guidance for the application of ARFPRC technology in mining areas with the risk of floor water inrush.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5804-5822"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148768447","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":"Influence of Relative Permeability Formulations on Hysteretic Seepage in Infinite Slopes","authors":"Diana Bianchi,&nbsp;Rossella Bovolenta,&nbsp;Martino Leoni,&nbsp;Domenico Gallipoli","doi":"10.1002/nag.70406","DOIUrl":"10.1002/nag.70406","url":null,"abstract":"<p>The paper investigates the impact of two classical relative permeability formulations, one expressed in terms of pore water pressure and the other in terms of degree of saturation, on the predicted seepage in an infinite slope subjected to surface pore pressure variations. A hysteretic retention law is employed, inducing irreversible changes in hydraulic conductivity during wetting–drying cycles. Two soils with different gradings are analysed and, for each of them, the two relative permeability formulations are matched at a reference point so as to represent the same material. Despite this calibration, the corresponding permeability fields differ significantly over the hysteretic domain, by as much as five orders of magnitude. These discrepancies are mostly the consequence of retention hysteresis and are substantially reduced if a unique, non-hysteretic retention curve is selected. For reference, a slope model with constant permeability, equal to the maximum saturated value, is also evaluated. The magnitude and rate of the slope response to surface perturbations, the number of cycles to restore equilibrium and the characteristics of the steady state regime all depend on the chosen permeability formulation. Slopes with a constant permeability react rapidly and uniformly, whereas slopes with permeability depending on pore water pressure or degree of saturation display slower, more complex and depth-attenuated responses, particularly for the finer soil. Overall, the results highlight the need for consistent calibration of both relative permeability and retention hysteresis, while underscoring the limitations of constant permeability models in capturing unsaturated seepage, even when appropriate hysteretic retention laws are adopted.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5663-5679"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.70406","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148674997","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":"Analytical Method for Longitudinal Response of Shield Tunnels Under Arbitrary Loading Based on the Kerr Foundation Model","authors":"Yukai Zhou,&nbsp;Weiming Huang,&nbsp;Jinchang Wang,&nbsp;Xiaojia Ji,&nbsp;Guobo Wang,&nbsp;Jinpin Chen","doi":"10.1002/nag.70413","DOIUrl":"10.1002/nag.70413","url":null,"abstract":"<div>\u0000 \u0000 <p>The operational safety of existing shield tunnels is increasingly threatened by construction activities in the proximity domain. Accurate prediction of their longitudinal response is thus of great importance. Previous research usually adopted the Winkler or Vlasov foundation models coupled with Euler-Bernoulli beams to establish analytical solutions. However, the solutions were tailored solely to and available for a specific disturbance scenario, which required an extension when encountering diverse ones. In addition, they are insufficient to capture the behavior of the soil continuum. To address these issues, this paper presents a comprehensive model by integrating the more advanced three-parameter Kerr foundation model with Timoshenko beam theory. By introducing the state-space method, a unified solution procedure is developed for predicting the tunnel response. The method's validity and versatility are demonstrated through comparisons with field measurements and existing solutions under three distinct scenarios: new tunnel undercrossing, adjacent pit excavation, and surface surcharge. Finally, a sensitivity analysis for the excavation case examines the influences of soil elastic modulus, pit-tunnel angle, and burial depth. The proposed framework offers a unified and efficient approach applicable to various disturbances, providing an efficient tool for the prediction of shield tunnel deformation.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5789-5803"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148754634","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":"DOE‐Based Sensitivity Analysis for Rational Constitutive Parameter Selection in Coupled Hydro‐Mechanical Simulations of Test in Unsaturated Soils","authors":"Isabela Augusto Silveira,Higor Biondo de Assis,Roger Augusto Rodrigues,Caio Gorla Nogueira","doi":"10.1002/nag.70455","DOIUrl":"https://doi.org/10.1002/nag.70455","url":null,"abstract":"ABSTRACT The Barcelona Basic Model (BBM) is the most well‐known elasto‐plastic constitutive model for describing the behavior of unsaturated soils. However, the selection of BBM parameters from laboratory data remains challenging due to the strong coupling and interdependence among the constitutive parameters governing hydro‐mechanical behavior. In this context, the Design of Experiments (DOE) methodologies provide a systematic approach for quantifying parameter influence and reducing the dimensionality of calibration problems. This paper presents a sensitivity analysis of the parameters used in the coupled hydro‐mechanical numerical simulation of a suction‐controlled oedometer test for three distinct soils, enabled by a computational framework designed to automatically execute and manage many simulations. The analyses were performed using both exhaustive Full Factorial (FF) designs and reduced Plackett–Burman (PB) matrices. The investigated variables included BBM mechanical parameters and hydraulic parameters, while the response variable was the simulated vertical strain during oedometric loading. The FF analyses were used as reference and PB matrices with different numbers of trials were subsequently assessed regarding their ability to reproduce the patterns obtained from the FF analyses with fewer simulations. The results showed that the soil stiffness parameters λ(0) and r , together with the historical variable p 0 * , had the greatest influence on vertical strain across all three soils. In contrast, hydraulic parameters had little to no effect on the response variable in this case. The results demonstrate that DOE methodologies provide an effective and computationally efficient framework for sensitivity assessment and parameter prioritization in coupled hydro‐mechanical simulations involving unsaturated soils.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"32 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"149090378","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 Multiscale Numerical Method for Studying Resaturation Phenomena in the COx Claystone Damaged Zone","authors":"Sarah Abou Chakra,&nbsp;Benoît Bary,&nbsp;Eric Lemarchand,&nbsp;Sylvie Granet,&nbsp;Jean Talandier","doi":"10.1002/nag.70412","DOIUrl":"10.1002/nag.70412","url":null,"abstract":"<div>\u0000 \u0000 <p>In the context of the planned deep geological disposal of nuclear waste in the Callovo-Oxfordian claystone (COx) in France, this research focuses on the long-term behavior of the host rock damaged during tunnel excavation. A general multiscale numerical methodology is developed and applied as a first step toward analyzing the hydro-mechanical response of a representative damaged zone surrounding the excavated drifts during resaturation, while accounting for swelling phenomena. A nonlinear viscoelastic-viscoplastic numerical model incorporating irreversible creep mechanisms is employed to describe the clayey rock at the material scale. The evaluation of the damaged matrix behavior is conducted using 3D representative elementary volumes (REVs) of the fractured host rock, consisting of the clay matrix containing dispersed penny-shaped cracks. A multiscale numerical technique based on a modified Finite Element squared (FE²) method is then used to simulate the response of a macroscopic gallery region, which is discretized into zones sharing similar crack density parameters estimated from experimental data. This approach is intended to provide a first-order assessment of the damaged zone response, and constitutes an initial contribution toward a more comprehensive multiscale modeling framework.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5738-5758"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148687830","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":"Research on the Macro-Micro Energy Evolution Law and Segmented Damage Constitutive Model of Sandstone Under Dry-Wet Cycles","authors":"Senlin Gao,&nbsp;Qingyang Ren,&nbsp;Bin Chen,&nbsp;Songqiang Xiao,&nbsp;Hang Song,&nbsp;Yanping Jia","doi":"10.1002/nag.70399","DOIUrl":"10.1002/nag.70399","url":null,"abstract":"<div>\u0000 \u0000 <p>The periodic fluctuation of reservoir water levels induces dry-wet cycles, deteriorating sandy rock slope stability and potentially triggering collapses. This study investigates the degradation mechanism of Three Gorges Reservoir fine sandstone through uniaxial compression tests and PFC2D simulations under varying dry-wet cycles. Key findings include: (1) Increasing cycles reduce uniaxial compressive strength and elastic modulus (showing “V-shaped” and “N-shaped” degradation trends), while permeability and porosity rise, and P-wave velocity declines. (2) Crack growth is nonlinear, dominated by shear cracks and high-angle microcracks, with force chains aligning with the loading direction. (3) Dissipation energy rate follows a “W-shaped” trend, while elastic energy rate exhibits an “M-shaped” pattern; both energies at crack initiation, damage, and peak stress correlate exponentially with cycle count (<i>K</i>_sd being more sensitive). (4) Dry-wet cycles weaken intergranular bonds, reduce elastic energy storage (<i>U</i>^e), increase dissipated energy (<i>U</i>^d), lower <i>M</i>-value stability, and shift failure from brittle to ductile. (5) A Weibull-based segmented damage model effectively simulates sandstone behavior under cyclic conditions. These insights enhance understanding of reservoir slope stability under hydrological fluctuations.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5823-5846"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148768446","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":"An Interfacial Flow Contact Model for Axisymmetric Consolidation Analysis of Unsaturated-Saturated Soils Under Continuous Drainage Boundaries","authors":"Dansheng Yu,&nbsp;Minjie Wen,&nbsp;Pan Ding,&nbsp;Yiming Zhang,&nbsp;Yuhui Xu,&nbsp;Manyu Yang","doi":"10.1002/nag.70408","DOIUrl":"10.1002/nag.70408","url":null,"abstract":"<div>\u0000 \u0000 <p>The tortuous morphology of fluid channels and inherent fluid viscosity cause significant friction and air resistance at the unsaturated-saturated soil interface. This resistance significantly reduces fluid seepage velocity, leading to an interfacial flow contact resistance effect. To address this, a flow contact transfer coefficient (<span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>R</mi>\u0000 <mi>w</mi>\u0000 </msub>\u0000 <annotation>${{R}_w}$</annotation>\u0000 </semantics></math>) and a flow partition coefficient (<span></span><math>\u0000 <semantics>\u0000 <msub>\u0000 <mi>η</mi>\u0000 <mi>w</mi>\u0000 </msub>\u0000 <annotation>${{eta }_w}$</annotation>\u0000 </semantics></math>) are introduced to establish a general imperfect flow contact model. This model quantitatively describes the resistance encountered by pore water as it seeps across interface channels. Based on this model, semi-analytical solutions for pore water pressure, pore air pressure, and foundation settlement of unsaturated-saturated soil foundations with vertical drains under continuous drainage conditions are derived by employing the Laplace transform and the finite Hankel transform. The consolidation behavior of foundations improved by vertical drains is then systematically investigated. Results indicate that the interfacial flow contact resistance effect significantly impedes pore water discharge, inducing a maximum pore water pressure gradient of 32.2 kPa at the interface. Compared with the complete flow contact model, the proposed model exhibits significant consolidation retardation. Specifically, at a depth of 8m, the initial dissipation of pore water pressure is delayed by approximately 11.1 h, and the complete dissipation period is extended by 46.29 days. Furthermore, under weak radial drainage conditions, the time required to reach a settlement of 183 mm is delayed by approximately 69 days.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5680-5700"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148674999","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":"Probabilistic Evaluation of Heterogeneous Landslide Influence Zones Accounting for Stratigraphic Dips With Borehole Data","authors":"Jian-Ping Li,&nbsp;Shui-Hua Jiang,&nbsp;Jian-Hong Wan,&nbsp;Guo-Tao Ma,&nbsp;Mohammad Rezania,&nbsp;Jingjing Meng","doi":"10.1002/nag.70422","DOIUrl":"10.1002/nag.70422","url":null,"abstract":"<div>\u0000 \u0000 <p>Accurate prediction of landslide runout and hazard zones is crucial for effective disaster risk management. Current studies often overlook complex soil structures by using stationary isotropic or transversely anisotropic unconditional random fields (RFs) in landslide post-failure modeling. Limited attention has been given to using borehole data for enhancing the accuracy of post-failure behavior predictions in slopes. To address these issues, this study presents a novel framework to evaluate landslide hazard zones using conditional RFs. It integrates enhanced Bayesian Updating with Structural Reliability Methods (BUS) to infer soil parameter distributions, identify dips and capture soil nonstationarity. The process of landslide is simulated using the generalized interpolation material point (GIMP) method. Additionally, an automated strategy is proposed for the borehole location selection to ensure that predictions align with actual conditions. Results indicate that this method improves prediction accuracy for runout and influence distances by using borehole data and reduces uncertainty compared to unconditional RFs, thereby simplifying decision-making and lowering control costs. Moreover, fewer boreholes are required for predicting runout distances compared to influence distances. This study highlights the necessity of considering complex heterogeneity and integration of borehole data in landslide risk management.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5947-5964"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148836258","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":"Effects of Ambient Moisture and Large-Strain Dynamic Preloading on the Dynamic Properties of the EAC-1A Regolith Simulant","authors":"Christos Vrettos","doi":"10.1002/nag.70409","DOIUrl":"10.1002/nag.70409","url":null,"abstract":"<p>The EAC-1A lunar mare simulant, a silty sand, was selected by the European Astronaut Center (EAC) for use in the LUNA analogue facility's lunar simulation testbed. Its stress- and strain-dependent shear modulus and damping ratio were determined via resonant column testing on specimens prepared either from ambient moisture material or oven-dried material. Significant differences were observed across the strain range. The tests further examined specimen behaviour during unloading from large to small strains and subsequent reloading. The reloading curves were observed to merge with the virgin response curves at the strain corresponding to load reversal, while a notable degradation of the small-strain modulus due to prior loading was evident. The overall behaviour in terms of shear modulus and damping was captured using simple equations.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5701-5711"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.70409","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148675131","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":"Analytical Methodology to Determine Thermally Induced Soil Deformation From Laboratory Data","authors":"Jibril B. Coulibaly,&nbsp;Alessandro F. Rotta Loria","doi":"10.1002/nag.70407","DOIUrl":"10.1002/nag.70407","url":null,"abstract":"<p>Experimental laboratory investigations underlie most of the knowledge about the mechanics of soils under non-isothermal conditions and their thermally induced deformation. Despite advances, the understanding of thermally induced deformation of soils remains inconsistent, hindering the development of a coherent framework that can advance science and engineering. Considering this challenge, this paper presents an analysis of the established methodology to determine thermally induced soil deformations from laboratory experiments. Established equations are re-derived, verified, and compared, and the resulting predictions are analyzed. Data support that the quantification of thermally induced soil deformations can be affected by volume corrections and calibrations, the integration of conservation equations under non-isothermal conditions, the misuse of thermal expansion coefficients, and the lack of repeatability tests and uncertainty quantification. Quantitative analyses support the conclusion that these aspects contribute to the contradictions in the state-of-the-art. Therefore, this work proposes an updated methodology to experimentally quantify the thermally induced deformation of soils, whose widespread adoption is supported by the provision of an open-source toolbox that implements the proposed changes to the state-of-the-art: The Thermally Induced Deformation Analysis Library, TIDAL.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 16","pages":"5868-5890"},"PeriodicalIF":3.6,"publicationDate":"2026-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.70407","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148768444","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}