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

{"title":"Computational Analysis of the Vibrational Response of Lunar Regolith Made of Nonspherical Grains","authors":"Lewis Traquair,Peter Watson,Ali Anwar,Sebastien Vincent Bonnieu,Marcello Lappa","doi":"10.1002/nag.70419","DOIUrl":"https://doi.org/10.1002/nag.70419","url":null,"abstract":"ABSTRACT Vibration‐based systems have shown that appropriately oriented and tuned mechanical oscillations can disrupt regolith interlocking and promote flow. Building on this evidence, the present study investigates vibration‐controlled granular discharge as a potential ISRU (In Situ Resource Utilization) strategy, with the aim of supporting predictable material processing through physically consistent numerical modeling. The equations governing the dynamics of nonspherical particles are presented accordingly and used to simulate the behavior of a fixed mass of granular material flowing through a hopper device in the framework of the discrete element method (DEM). The results show that even for monodisperse distributions vibration‐induced fluidization does not follow simple monotonic trends. The mass flow rate initially increases with the forcing frequency f , reaching a plateau in the range of 150–200 Hz, beyond which further frequency increases produce no significant enhancement. At the same time, higher values of the dimensionless acceleration Γ generally reduce the mass flow rate and cause an increase in the frequency threshold above which vibration‐driven discharge exceeds gravity‐driven flow. Moreover, for fixed, moderately low values of Γ, increasing f leads to the emergence of higher‐order harmonics in the instantaneous mass flow rate. For fixed and relatively high values of f , a larger value of Γ can cause the appearance of subharmonics in the spectrum. In contrast, for relatively low f and high Γ, incommensurate spectral components appear, which reduce the regularity and predictability of the flow response. These findings are supported by a detailed analysis of particle patterning behavior and frequency spectra, which reveals the fine‐scale dynamics underlying the observed flow regimes.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"82 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894678","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 Approach to Long‐Term Stability Analysis in Rock Engineering","authors":"Muhammad Shoaib,Xin‐Dong Wei,Gao‐Feng Zhao,Xifei Deng","doi":"10.1002/nag.70427","DOIUrl":"https://doi.org/10.1002/nag.70427","url":null,"abstract":"ABSTRACT Long‐term stability of rock engineering structures is controlled by progressive strength degradation under sustained loading, yet direct experimental determination of long‐term strength (LTS) is severely constrained by impractically long testing durations and strong material variability. To address this limitation, this study proposes a feasible numerical framework for long‐term stability analysis by linking time‐dependent microscale damage evolution to engineering‐scale strength parameters. The four‐dimensional lattice spring model (4D‐LSM), incorporating creep, plasticity, and progressive bond fracture, is calibrated using conventional creep test data, in which the maximum bond‐fracture threshold governs long‐term failure time. Time effects are introduced through temporal evolution of model input parameters, enabling numerical creep tests to reproduce stress‐time‐to‐failure behavior and predict LTS. The numerically obtained LTS data are fitted with an empirical stress‐time relationship and further converted into time‐dependent equivalent Mohr‐Coulomb cohesion and friction angle. These time‐varying parameters are then implemented in a classical strength reduction finite element framework, and a machine learning model is adopted to establish the correlation between the factor of safety and the degrading strength parameters. The proposed approach facilitates the practical prediction of long‐term stability of rock engineering structures using temporally evolving strength parameters derived from 4D‐LSM simulations, thereby providing a mechanistically sound and computationally efficient tool for the assessment of delayed failure in rock engineering.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"11 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894680","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":"Faulted Rock Slope Stability Under Blasting Excavation Using Copula‐Based Rotated Anisotropic Random Fields","authors":"Jiajun Wu,Chong Yu,Haibo Li,Yongan Ma,Renjie Wu","doi":"10.1002/nag.70429","DOIUrl":"https://doi.org/10.1002/nag.70429","url":null,"abstract":"ABSTRACT Blasting excavation can induce permanent displacement in faulted rock slopes and adversely affect slope stability. Most existing stability assessments use deterministic approaches and therefore do not adequately account for the spatial variability and cross‐correlation of geomechanical properties. In addition, conventional assumptions of transverse anisotropy and the arbitrary selection of copula models may bias blast‐design assessments. This study develops a stochastic framework that integrates rotated anisotropic random fields with copula models to characterize spatially variable rock properties. The effects of copula selection and rotated anisotropy on the statistical distributions of permanent displacement and a displacement‐based risk index are investigated. The framework is demonstrated through an engineering case study and compared with conventional deterministic analysis. The results indicate that the commonly used Gaussian copula tends to overestimate slope stability, whereas strata orientation, scale of fluctuation, and the correlation between cohesion and friction angle significantly affect the response statistics. The case study illustrates the applicability of the framework to comparative stability assessments of blast‐affected faulted rock slopes.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"5 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894681","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":"Mechanical Stability of Granular Assemblies: A Matter of Scale","authors":"Xiaoxiao Wang,Yang Liu,Antoine Wautier,François Nicot","doi":"10.1002/nag.70431","DOIUrl":"https://doi.org/10.1002/nag.70431","url":null,"abstract":"ABSTRACT Granular assemblies are prone to form self‐organized dissipative structures. The microscopic interparticle sliding, entailing frictional dissipation at contacts drives configurational changes that critically influence the mechanical stability of the system. Using the discrete element method (DEM), this study investigates the mechanical stability of granular assemblies based on a mesoscale second‐order work criterion. In dense specimen, the vanishing of the mesoscopic second‐order work exhibits a spatial localization after the stress peak, consistent with the formation of a shear band. By analyzing the evolution of stability indexes, the sub‐domain inside the shear band undergoes a fluctuating regime of stability, which governs the fluctuating stability characteristic of dense specimen at the critical state. For the loose specimen, this fluctuating stability regime maintains throughout loading up to the critical state. The shear band domain and the whole loose specimen were identified as inherent dissipative structures, governing the fluctuating stability regime, where the different loop categories composing the topology of the granular assembly share similar stability mechanisms. Furthermore, the dissipative structure exhibits a similar mean characteristic size of unstable meso‐loops at the critical state regime. The meso‐loops with sufficiently large instability influence range contribute to a loss of stability of the dissipative structures, ultimately controlling macroscopic instability of granular assembly.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"28 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894682","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 Phase‐Field Model for Quasi‐Dynamic Rupture Nucleation and Propagation of In‐Plane Faults","authors":"Fan Fei,Md Shumon Mia,Ahmed E. Elbanna,Jinhyun Choo","doi":"10.1002/nag.70416","DOIUrl":"https://doi.org/10.1002/nag.70416","url":null,"abstract":"ABSTRACT Computational modeling of faulting processes is an essential tool for understanding earthquake mechanics but remains challenging due to the structural and material complexities of fault zones. The phase‐field method has recently enabled unified modeling of fault propagation and off‐fault damage within a fracture‐mechanics‐based framework. However, its capability has been restricted to simplified anti‐plane settings. In this study, we extend the phase‐field method to in‐plane faulting by introducing two key advancements: (i) the incorporation of enhanced fault kinematics and pressure‐dependent shear strength for a more accurate representation of fault behavior, and (ii) a revised fault propagation criterion that explicitly accounts for the coupling between shear strength and normal stress. The proposed formulation is verified against standard discontinuous approaches to quasi‐dynamic fault rupture under in‐plane conditions and validated using experimental observations and numerical data on nucleation and propagation of fault slip surfaces and associated localized damage zones. Simulations incorporating structural complexities and material heterogeneities demonstrate the robustness and versatility of the phase‐field model, establishing it as a powerful tool for investigating the interactions between fault zone properties and earthquake processes.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"492 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894679","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":"Physics‐Guided Neural Regression Modeling of Energy‐Balance Proxies for Hydraulic Fracture Propagation in 1D–3D","authors":"Dennis Delali Kwesi Wayo,Lei Wang,Sonny Irawan,Leonardo Goliatt,Randy Hazlett","doi":"10.1002/nag.70418","DOIUrl":"https://doi.org/10.1002/nag.70418","url":null,"abstract":"ABSTRACT Energy‐balance surrogates can support early hydraulic‐fracture screening before fully coupled simulations are run. This paper defines a physics‐guided neural regression benchmark for controlled energy‐proxy learning across increasing input dimensionality. We generate samples with fracture length , fracture width , pressure gradient , and fluid viscosity with numerical extrema, compute , add noise, and train on . Inputs are (1D), with fixed (2D), and (3D); features are standardized, with target standardization only for 1D/2D. The 1D/2D networks use 4256 ReLU layers, and the 3D network uses 3 128 Swish layers selected from a limited practical sweep. On noise‐free proxy test targets after noisy‐target training, MSEs are 0.0031 (1D), 1.4580 (2D), and 10917.8261 (3D), with NRMSE 0.0004/0.0069/0.0918 and 1.0000/1.0000/0.9916. KGD‐ and PKN‐inspired analytical trends are used as reference controls to separate the proposed attenuation‐driven energy proxy from classical opening‐dominated fracture scalings. Their poor agreement with the proxy shows that the benchmark is a distinct surrogate task rather than a reduced KGD/PKN solution. The result is a reproducible benchmark for assessing dimensional scaling, training stability, and surrogate accuracy in hydraulic‐fracture energy modeling before extension to coupled elasticity–fluid‐flow validation.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"8 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148894721","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 Total Displacement Catastrophe Criterion for Assessing Open‐Pit Slope Stability Under Rainfall Conditions","authors":"Lu He, Jianchuan Ren, Yuqi Shang, Gaofeng Song, Dezhong Kong, Di Wu","doi":"10.1002/nag.70430","DOIUrl":"https://doi.org/10.1002/nag.70430","url":null,"abstract":"The aim of the present study is to address the challenge of accurately assessing the stability of open‐pit mine slopes under rainfall infiltration. This study innovatively combines total displacement theory with the cusp catastrophe model. A slope instability criterion based on total displacement mutation is proposed. The MATLAB platform was used for model validation and comparative analysis. Firstly, a quartic polynomial model was developed to characterize the relationship between the total displacement and the reduction coefficient based on the strength reduction method and total displacement theory. The equilibrium and discriminant equations for the cusp catastrophe model were then derived. Secondly, a typical mine slope in Guizhou Province, China, was selected as a case study to assess the applicability and accuracy of three instability criteria: the plastic zone penetration criterion, the key point displacement catastrophe criterion, and the total displacement catastrophe criterion (TDCC). The results indicate that the slope safety factor calculated by the TDCC is 1.201. This criterion is objective and independent of the monitoring point location, significantly reducing human errors. The criterion's accuracy was further verified by varying the cohesion and the internal friction angle. The results obtained from the TDCC were stable and central, aligning with the variation trends of the other two criteria, thereby confirming its reliability and accuracy. This work provides a quantitative, operational method for identifying instability in slope rainfall‐infiltration stability analysis. It is of theoretical and engineering guidance significance for the prevention and control of slope disasters in open‐pit mines.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"51 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148860892","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 Non‐Isothermal Viscoplastic Framework for Assessing Shaft Resistance in Energy Geostructures","authors":"Saeed Tourchi, Arash Alimardani Lavasan","doi":"10.1002/nag.70426","DOIUrl":"https://doi.org/10.1002/nag.70426","url":null,"abstract":"This study presents a coupled thermo‐hydro‐mechanical (THM) num erical investigation of shaft resistance degradation in geothermal piles embedded in low‐permeability clays. A thermo‐viscoplastic constitutive model is developed to simulate the pile–soil interface behavior under non‐isothermal conditions. The model captures the evolution of excess porewater pressure, effective stress reduction, and their combined impact on shaft bearing capacity during thermal activation. Parametric simulations are carried out for a range of fluid temperatures (21–50) and soil permeabilities (– ), representative of in‐situ conditions. Results show that undrained heating leads to substantial excess porewater pressure buildup, especially at low permeability and high temperature, resulting in significant reductions in zero‐thickness interface element shear strength. A shaft resistance reduction function is proposed to quantify this phenomenon. The numerical predictions are benchmarked against full‐scale energy‐pile measurements, showing close agreement in temperatures, strains, uplift, and stress. The findings underline the critical role of thermal pressurization in energy pile design and provide an enhanced modeling strategy for performance evaluation in fine‐grained soils.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"43 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148860893","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":"Thermal‐Mechanical Solution of Saturated Soft Soil around a Semi‐Permeable Tunnel Considering Internal Heat Effects","authors":"Senlin Xie, Jiajun Niu, Changqing Xia, Anfeng Hu, Liwu Wang, Zhengwei Jiang, Xiangsheng Chen","doi":"10.1002/nag.70423","DOIUrl":"https://doi.org/10.1002/nag.70423","url":null,"abstract":"Urban tunnels operate for decades under complex service conditions where mechanical loading and thermal disturbance are invariably intertwined. In practice, long‐term train operation, ventilation and air‐conditioning systems, energy‐tunnel heat exchange, and even accidental fire scenarios continuously inject heat into tunnel linings and the surrounding saturated ground. These realities create a coupled thermo‐hydro‐mechanical environment around tunnels that is far more complex than what is captured by classical consolidation theories, and the long‐term implications for ground deformation and tunnel serviceability remain insufficiently understood. To address this gap, this study develops an analytical model for the thermo consolidation of soils surrounding tunnels by extending the Terzaghi–Rendulic consolidation framework to incorporate heat diffusion, temperature‐dependent permeability, and soil viscous effects. On the basis of the derived closed‐form solution, the study systematically explores the roles of tunnel boundary temperature, thermal conductivity, tunnel depth and radius, tunnel permeability, and soil viscosity on temperature fields, excess pore‐water pressure evolution, and settlement response. The results reveal that tunnel‐induced thermal effects predominantly reshape pore‐pressure dissipation pathways and consolidation timescales by altering the spatiotemporal distribution of permeability. Geometric parameters act in a unified manner on both heat transfer and consolidation by modifying the effective overburden thickness. Tunnel permeability emerges as a key control on the consolidation rate, whereas the soil viscosity coefficient dictates whether settlements exhibit persistent, rheology‐driven long‐term development. By providing a tractable yet physically enriched analytical framework, this work offers a theoretical basis for evaluating the long‐term performance of tunnels subjected to combined operational thermal loads and sudden surcharges.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"50 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148860894","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, Shui‐Hua Jiang, Jian‐Hong Wan, Guo‐Tao Ma, Mohammad Rezania, Jingjing Meng","doi":"10.1002/nag.70422","DOIUrl":"https://doi.org/10.1002/nag.70422","url":null,"abstract":"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.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"39 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2026-08-27","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}