{"title":"Numerical Modelling Ice Lens Formation and Frost Heave in Unsaturated Soils","authors":"Zili Wang, Jidong Teng, Satoshi Nishimura, Sheng Zhang, Daichao Sheng","doi":"10.1002/nag.70091","DOIUrl":"https://doi.org/10.1002/nag.70091","url":null,"abstract":"Frost heave is a typical thermo‐hydro‐mechanical coupling process, which can lead to significant threats to the infrastructures in cold regions. The initiation and growth of ice lenses is the core issue to understand the process of frost heave. But this process has not been well modelled in previous studies. It is hard to accurately model the formation of ice lenses and the contribution of the vapour phase during the freezing process. This study presents a novel frost heave model that accounts for the effects of vapour migration in unsaturated soils, which is implemented in a C++ simulation programme with an interactive user interface. The model reveals periodic frost heave variations due to cooling gradients and highlights the impact of vapour on ice lens formation under low moisture conditions. It also demonstrates the relative importance of liquid and vapour fluxes with moisture content, with vapour migration playing a key role in frost heave in coarse‐grained soils. These new findings provide new insights into frost heave mechanisms and challenge traditional ice lens formation criteria. In addition, the model effectively simplifies boundary conditions and improves computational stability and efficiency. These advancements can improve frost heave prediction and deepen the understanding of soil freezing mechanisms, offering valuable insights for infrastructure applications.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"18 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145203160","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":"Intelligent Inversion of Multi‐Stratum Parameters in Shield Tunnels and Reliability Analysis of Tunnel Deformation Under Spatial Variability Conditions","authors":"Liu Yang, Jie Yang, Lin Cheng, Chunhui Ma","doi":"10.1002/nag.70094","DOIUrl":"https://doi.org/10.1002/nag.70094","url":null,"abstract":"Owing to deposition, weathering, and historical loading variations, the mechanical properties of underground rock and soil masses demonstrate significant spatial variability and stratified distribution. This study investigates the influence of multi‐layer soil spatial variability on ground settlement and tunnel reliability during shield tunnel construction by developing a refined stochastic finite element model. The CPSO‐TLOOA‐Stacking hybrid intelligent algorithm optimizes the inversion of multi‐stratum mechanical parameters based on the measured surface settlement data from tunnel engineering and the Conditional Tabular GAN (CTGAN) data extension framework. Utilizing the Karhunen–Loève (K‐L) series expansion method and random field theory, a joint analysis framework of stochastic finite element and probability statistics is constructed to evaluate the impact of spatial random field parameters of different soil layers on formation deformation and failure probability. Coupled with the Hamiltonian Monte Carlo‐Subset Simulation algorithm, the reliability of tunnel deformation under conditions of cross‐correlated random fields with multiple surrounding rock parameters is effectively assessed. The results indicate that the <jats:italic>R<jats:sup>2</jats:sup></jats:italic> value of the expanded dataset fitted by the CPSO‐TLOOA‐Stacking hybrid intelligent algorithm is 99.46%, and the relative error between the dataset and the measured value is 0.7%. The Hamiltonian Monte Carlo‐Subset Simulation algorithm significantly enhances the calculation efficiency of tunnel deformation reliability and provides valuable guidance for shield tunnel construction and design.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"6 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145203161","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":"Damping Mode Adaptability and Benchmarks in Explicit Formulation of Discontinuous Deformation Analysis for Typical Geotechnical Scenarios","authors":"Qinglong Deng, Fei Zheng, Yu‐Yong Jiao, Zhi Li","doi":"10.1002/nag.70096","DOIUrl":"https://doi.org/10.1002/nag.70096","url":null,"abstract":"Appropriate choice of both damping mode and the associated damping coefficient significantly influence the stability, accuracy, and efficiency of the simulation results for geotechnical problems. However, in numerical analysis with explicit discontinuous deformation analysis (EDDA), appropriate damping settings for different geotechnical scenarios are rarely discussed. To seek optimized damping in EDDA for various scenarios, four types of damping, including viscous damping, contact damping, unbalanced force damping, and Rayleigh damping, are incorporated in the EDDA framework. The performance of these four types of damping for basic and multi‐block systems is compared through three representative scenarios, including quasi‐static contact, collision, and wave propagation. This study verified the adaptability of the four damping models in certain geotechnical scenarios and provided recommendations for mode selection and parameters setting based on findings from the benchmarking tests. This study on damping adaptability facilitates the application of EDDA in geotechnical problems towards a more accurate and more efficient simulation.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"39 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145195025","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":"Benchmarking Contact Detection Algorithms Used in Polyhedral Particle System","authors":"Yuval Keissar, Michael Gardner, Nicholas Sitar","doi":"10.1002/nag.70088","DOIUrl":"https://doi.org/10.1002/nag.70088","url":null,"abstract":"A critical assessment of contact detection algorithms routinely used for simulating convex polyhedra in the Discrete Element Method is presented herein. Specifically, we focus on accuracy and computational efficiency and discuss the advantages and limitations of four different algorithms: the coupled Gilbert–Johnson–Keerthi – Expanding Polytope Algorithm (GJK‐EPA), Fast Common Plane coupled with the Intersection Polygon Method (FCP‐IPM), GJK‐EPA coupled with IPM (GJK‐EPA‐IPM), and the Linear Programming Algorithm (LP). These algorithms were implemented within the same open source software framework to allow an objective assessment of their performance. Although in this case the Discrete Element Method (DEM) is used as the particulate solver, the characteristics of the selected contact detection algorithms are independent of this choice. The results show that the GJK‐EPA algorithm, although often preferred for its speed, lacks the necessary accuracy to reproduce the physics of even a very simple dynamic problem. The same issue arises with all the other algorithms except the modified FCP‐IP algorithm. These results raise serious concerns about the veracity of results of dynamic analyses involving simulations involving large assemblies of polyhedral particles in which any of these algorithms have been or are being used. Unfortunately, in many situations, an empirical assessment of the behavior of the assembly will have masked the inherent inaccuracy of the results.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"97 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145188377","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}
Chuanhai Zhan, Ming Sun, Zhiqing Zhang, Hao Liu, Yunpeng Zhang, M. Hesham El Naggar, Meijuan Xu, Wenbing Wu
{"title":"Torsional Vibration of a Pipe Pile in Arbitrary Layered Saturated Soil Based on the Additional Mass Model","authors":"Chuanhai Zhan, Ming Sun, Zhiqing Zhang, Hao Liu, Yunpeng Zhang, M. Hesham El Naggar, Meijuan Xu, Wenbing Wu","doi":"10.1002/nag.70093","DOIUrl":"https://doi.org/10.1002/nag.70093","url":null,"abstract":"This study systematically investigates the torsional vibration characteristics of pipe piles embedded in arbitrary layered saturated soils using an additional mass model to account for soil plug effects. A coupled pile‐soil torsional vibration model is developed by integrating Biot's poroelastic theory and the additional mass model, considering viscous damping in both the pile shaft and surrounding soil. Analytical solutions for the torsional dynamic complex impedance at the pile head in the frequency domain are derived using Laplace transforms, separation of variables, and impedance function recursion. The degenerate form of the solution is validated against existing results for solid piles in saturated soils, confirming its accuracy. Parametric analyses are conducted to evaluate the effects of pile geometry, soil properties, and soil plug characteristics on the torsional dynamic stiffness and damping of the pile head, providing insights into the influence mechanisms of key factors in layered saturated soil systems.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"57 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145188376","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}
Piao Miao, Fusheng Zha, Qiao Wang, Shijin Feng, Hongxin Chen, Lin Mu
{"title":"Innovative Sealing Approaches for Karst Flowing Groundwater: Numerical Analysis of Secondary Grouting Mechanism After Grout‐Filled Bag Sealed","authors":"Piao Miao, Fusheng Zha, Qiao Wang, Shijin Feng, Hongxin Chen, Lin Mu","doi":"10.1002/nag.70092","DOIUrl":"https://doi.org/10.1002/nag.70092","url":null,"abstract":"Groundwater pollution in karst regions poses a serious threat to the ecological environment and public health. Effective sealing of karst conduits with flowing water still remains a significant engineering challenge at present. This study introduces an innovative sealing technique combining grout‐filled bag plugging with secondary grouting. Mechanism of grout filling and diffusion during secondary grouting in flowing water at the gap between grout‐filled bag and bedrock boundary of the karst conduit was analyzed through computational fluid dynamics simulations. Effect of water flow velocities, grouting rates, and grout‐filled bag dimensions was investigated. Our results indicate that water flow velocity is a pivotal factor in secondary grout sealing efficacy. Grouting effectiveness is significantly impacted when water flow velocity exceeds 0.1 m/s. However, increasing grouting speed can counteract this effect. When grouting speed exceeds water flow velocity, the grout filling ratio increases. Grout filling ratio reaches its maximum when gap sizes range from 0.7 to 0.8 m. Under identical grout flow velocities, the differences in viscosity and density between water and grout create uneven pressure distribution. This causes water to compress and displace the injected grout, reducing the grout filling ratio within the gap. Notably, when grouting speed surpasses or equals the flowing water rate, a dynamic sealing mass forms upstream of the bag, providing additional pathways for grout migration and significantly bolstering the sealing effect. These insights are instrumental for the construction design of large leakage conduit sealing, especially within karst regions.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"155 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145188378","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 Efficient Hybrid Simulation Method to Quantitatively Analyze the Flexible Wire‐Ring Net Rockfall Protection Systems","authors":"Yuntao Jin, Zhixiang Yu, Liping Guo, Linxu Liao, Lijun Zhang, Liru Luo","doi":"10.1002/nag.70090","DOIUrl":"https://doi.org/10.1002/nag.70090","url":null,"abstract":"Flexible wire‐ring net systems are commonly utilized in rockfall protection engineering on geologically unstable slopes. Due to the complex nonlinear behaviors between discrete wire rings, significant computational resources are required for contact searching between nodes in numerical models. Coupling computations of large‐scale systems consisting of flexible structures, rockfalls, and actual terrain still present challenges. In this study, a hybrid simulation method is proposed, combining the high‐precision coordination model with the efficient membrane equivalent method for wire‐ring nets based on an established partitioning principle. Hybrid transition boundaries were built to accurately reproduce impact responses while significantly improving computational efficiency. The hybrid method was validated by establishing a flexible barrier and an attenuator system compared with full‐scale test results, with 11.9% impact force error for the flexible barrier and 8.4% residual kinetic energy error for the attenuator system. This method enables the prediction of rockfall trajectories, energy evolutions, and dynamic responses of structural components. Compared to pure circular beam models, the overall computational time for the attenuator system exceeding 2300 m<jats:sup>2</jats:sup> is reduced by 97.0%, facilitating quantitative design of large‐scale flexible protection systems in situ terrain.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"2 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145153779","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":"Analytical Solution for Stress and Displacement of the Arbitrary‐Shaped Shallow Tunnel Excavated in Orthotropic Rock","authors":"Yulin Zhou, Ning Zhang, Aizhong Lu","doi":"10.1002/nag.70086","DOIUrl":"https://doi.org/10.1002/nag.70086","url":null,"abstract":"Excavating tunnels in anisotropic rock induces greater stress concentrations at the excavation boundary compared to isotropic rock. However, existing analytical solutions for the shallow tunnel are founded on the simplifying assumption of isotropic rock masses. To clarify the deformation mechanism and mechanical behavior of the shallow tunnel, we proposed an analytical method for solving stress and displacement of the arbitrary‐shaped shallow tunnel excavated in orthotropic rock mass, incorporating the effects of body forces and anisotropy of rock mass. Proposing the specific forms of the analytical functions for the shallow tunnel in anisotropic rock, which can reflect body forces and the anisotropy of the rock mass. A linear equation system, derived from stress boundary conditions at the surface and tunnel excavation, is solved by the boundary collocation method. In the solution process, conformal transformation techniques are employed to construct three polar coordinate systems, which aid in resolving boundary conditions. Subsequently, we analyzed the effects of tunnel depth and anisotropic parameters on the stress and displacement of a horseshoe‐shaped tunnel and verified the correctness of the results through ANSYS software.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"96 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145181114","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":"Cross‐Scale Simulation Study on the Wetting Mechanical Properties and Deformation Characteristics of Rockfill Materials","authors":"Hui Chen, Donghai Liu","doi":"10.1002/nag.70073","DOIUrl":"https://doi.org/10.1002/nag.70073","url":null,"abstract":"The cross‐scale approach based on the finite element method (FEM) and discrete element method (DEM) was employed to investigate the effects of different particle shapes and stress conditions on axial deviatoric stress during wetting. A DEM modeling framework was proposed to simulate the wetting tests considering different particle shapes. The Effect of varying particle shapes on the stress‐strain and volumetric strain behavior of rockfill materials during the wetting process was investigated. Moreover, the study analyzed macroscopic characteristics such as axial and volumetric wetting deformation under triaxial test conditions. Furthermore, the mesoscopic mechanisms influencing the wetting behavior of rockfill materials were investigated using internal sample features such as Euler angles, coordination number, and contact fabric distribution. Finally, a wetting constitutive model incorporating particle shape characterization indices was proposed and validated across multiple scales. This study provides a novel approach for the accurate evaluation and safety control of the structural behavior of high rockfill dams.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"3 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145140657","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 Novel Compressive Sensing Approach for Optimizing Automated Monitoring of Excavation‐Induced Horizontal Displacements","authors":"Cheng Chen, Yang Lyv, Liang‐Tong Zhan, Xin‐Jiang Wei, Xing‐Wang Liu, Guan‐Nian Chen","doi":"10.1002/nag.70089","DOIUrl":"https://doi.org/10.1002/nag.70089","url":null,"abstract":"As excavation projects become increasingly complex, the demand for automated monitoring systems has risen due to their ability to provide continuous, real‐time data collection. However, traditional methods using automated inclinometers are cost‐prohibitive because they require a high number of sensors for accurate data acquisition. This study proposes a novel approach based on compressive sensing (CS) theory to interpret excavation‐induced horizontal displacement profiles using data from a reduced number of sensors. Validation with 19,311 displacement profiles from a 30.2‐m deep excavation project in Hangzhou, China, demonstrated the robustness of the method, achieving a maximum root mean square error (RMSE) of 6.42 mm (a 7.9% relative error for a maximum displacement of 80.9 mm), while reducing sensor deployment costs by a factor of 22 compared to traditional inclinometer techniques. The CS‐based approach consistently outperformed traditional regression models and proved effective across various sensor spacing scenarios. An analysis of 405,531 simulated cases provided an RMSE envelope, allowing engineers to balance accuracy and budget constraints when selecting sensor spacing. Additionally, comparative studies of sensor placement schemes revealed that while uniform spacing resulted in lower RMSE values and superior overall reconstruction, non‐uniform spacing more effectively captured maximum horizontal displacements, offering a cost‐efficient solution for applications that prioritize critical displacement monitoring.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"16 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145133437","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}