{"title":"A Stress‐Driven Double‐Phase–Field Framework for Tensile Fracturing Processes in Transversely Isotropic Rocks","authors":"Weihong Yuan, Yang Zhao, Bingyin Zhang","doi":"10.1002/nag.3830","DOIUrl":"https://doi.org/10.1002/nag.3830","url":null,"abstract":"We present a double‐phase–field framework for tensile fracturing processes in transversely isotropic rocks. Two distinct phase‐field variables are introduced to represent smeared approximations of tensile fractures along the weak bedding planes and through the anisotropic rock matrix, respectively. Driving forces that control fracture propagation in the phase‐field framework are constructed as a stress‐based formula with a recently developed tensile failure criterion that distinguishes the two failure modes in transversely isotropic rocks. For numerical implementation, we adopt a staggered integration scheme and decouple the governing equations so that the displacement field and phase‐field variables can be updated in sequence for a given loading step. The finite element formulation of the proposed framework is introduced in detail in this paper and is implemented in an in‐house finite element code. The numerical implementation is then validated by reproducing the uniaxial tension test results of Lyons sandstone. After that, we conduct simulations on a pre‐notched square plate loaded in tension to demonstrate the features of the proposed framework. Finally, we conduct simulations of three‐point bending tests of Pengshui shale and show that the proposed model can reproduce the force–displacement curves and failure patterns of specimens with different bedding plane orientations observed in laboratory experiments.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142275653","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}
Mengfan Zong, Jing Zhang, Wenbing Wu, Ziye Yu, Yi Zhang, Guoxiong Mei
{"title":"Semi‐Analytical Solution for One‐Dimensional Nonlinear Consolidation of Multilayered Soil Considering Self‐Weight and Boundary Time Effect","authors":"Mengfan Zong, Jing Zhang, Wenbing Wu, Ziye Yu, Yi Zhang, Guoxiong Mei","doi":"10.1002/nag.3839","DOIUrl":"https://doi.org/10.1002/nag.3839","url":null,"abstract":"The self‐weight stress in multilayered soil varies with depth, and traditional consolidation research seldom takes into account the actual distribution of self‐weight stress, resulting in inaccurate calculations of soil consolidation and settlement. This paper presents a semi‐analytical solution for the one‐dimensional nonlinear consolidation of multilayered soil, considering self‐weight, time‐dependent loading, and boundary time effect. The validity of the proposed solution is confirmed through comparison with existing analytical solutions and finite difference solution. Based on the proposed semi‐analytical solution, this study investigates the influence of self‐weight, interface parameter, soil properties, and nonlinear parameters on the consolidation characteristics of multilayered soil. The results indicate that factoring in the true distribution of self‐weight leads to a faster dissipation rate of excess pore water pressure and larger settlement and settlement rate, compared to not considering self‐weight. Both boundary drainage performance and soil nonlinearity have an impact on consolidation. If the boundary drainage capacity is inadequate, the influence of soil nonlinearity on consolidation diminishes.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142245219","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":"Modeling of Drain Consolidation in the Quick Triaxial Test and Its Analytical Solution","authors":"Zhibo Chen, Jungao Zhu, Xinjiang Zheng, Lei Wang","doi":"10.1002/nag.3842","DOIUrl":"https://doi.org/10.1002/nag.3842","url":null,"abstract":"Sand columns have been widely used to accelerate drainage and then improving the mechanical properties of soft soil foundations. The sand column has also been introduced into the triaxial test by researchers, in the center of the cylindrical specimen, to greatly accelerate drainage and consolidation process. The objective of this paper is to evaluate the consolidation properties of the triaxial cylindrical specimen considering the presence of a sand column, and then to propose a consolidation model that simulates the consolidation process of the triaxial test. The consolidation equations were derived considering the drainage of the specimen with a sand column composed of both vertical and double‐radial flows. Then the analytical solution of the model was obtained based on specific initial and boundary conditions. The comparison between the consolidation model and the laboratory tests yielded highly consistent. The case study demonstrated that the proposed consolidation model accurately simulates the evolution of average pore pressure and degree of consolidation in triaxial specimens containing a sand column. The studies on the consolidation parameters showed that there were different effects on the drainage rate for the diameter of specimen, the permeability coefficients of specimen and sand column, as well as the radius of the sand column.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142236620","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":"Micro‐Mechanical Analysis for Residual Stresses and Shakedown of Cohesionless‐Frictional Particulate Materials Under Moving Surface Loads","authors":"Wei Cai, Ping Xu, Runhua Zhang","doi":"10.1002/nag.3837","DOIUrl":"https://doi.org/10.1002/nag.3837","url":null,"abstract":"Residual stresses and shakedown have been successfully presented by two‐dimensional numerical experiments based on the discrete element method (DEM), wherein a cohesionless‐frictional material under moving surface loads was replicated through irregular‐shaped particles. With surface loads below the shakedown limit, both permanent deformations and residual stresses cease to accumulate and the numerical structure shakes down after a number of load passes. Corresponding micro‐mechanical analyses indicate that strong forces and normal forces make a dominant contribution to residual stresses. Besides, averaged magnitudes of interparticle forces and corresponding total contact numbers initially change with load passes, and their final variation trends will differ as the structure shakes down or not. Furthermore, polar distributions of interparticle forces and contacts have been presented, and variations of their preferential orientations were emphasised. Lastly, the fabric tensor and anisotropy of resultant forces were studied, presenting the anisotropy weakening of macro‐stress fields, induced by developments of residual stresses.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142236220","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 Longitudinal Seismic Responses of Circular Tunnel Crossing Fault Zone","authors":"Jie Tang, Manchao He, Hanbing Bian, Yafei Qiao","doi":"10.1002/nag.3841","DOIUrl":"https://doi.org/10.1002/nag.3841","url":null,"abstract":"This paper proposes a simplified analytical solution for longitudinal seismic responses of a circular tunnel crossing a fault zone under longitudinally propagating shear waves. The transmissions and reflections of shear waves at two geological interfaces between the fault zone and intact rock are considered when calculating the free‐field displacement. An improved elastic foundation beam model considering different tangential contact conditions at the tunnel‒rock interface is also adopted. According to the continuous conditions at the two geological interfaces, explicit expressions for the tunnel displacement, bending moment, and shearing force are given. The effectiveness of the proposed analytical solution is validated via numerical simulations, and the importance of accounting for tangential contact conditions at the tunnel‒rock interface is emphasized. Moreover, parametric studies are performed to investigate the effects of the fault zone width, rock conditions, tunnel lining stiffness, tangential contact conditions, and earthquake frequency on the deformation and internal forces of tunnels subjected to seismic waves. This novel analytical solution can be utilized to quickly estimate the longitudinal seismic responses of circular tunnels crossing fault zones subjected to longitudinally propagating shear waves, particularly in the preliminary engineering design, and can be extended to geological conditions with multiple interfaces.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142236180","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 Solutions for a Fully Coupled Hydraulic‐Mechanical‐Chemical Model With Nonlinear Adsorption","authors":"Lin Han, Zhihong Zhang, Jiashu Zhou","doi":"10.1002/nag.3829","DOIUrl":"https://doi.org/10.1002/nag.3829","url":null,"abstract":"Adsorption characteristics play a crucial role in solute transport processes, serving as a fundamental factor for evaluating the performance of clay liners. Nonlinear adsorption isotherms are commonly found with metal ions and organic compounds, which introduce challenges in obtaining analytical solutions for solute transport models. In this study, analytical solutions are proposed for a fully coupled hydraulic‐mechanical‐chemical (HMC) model that accounts for both the Freundlich and Langmuir isotherms. To mitigate the difficulties arising from the variable coefficients, the system of second‐order partial differential equations involving three variables is linearized. The method of separation of variables, theory of integration, and Fourier series are utilized to derive analytical solutions. The analytical method presented can potentially be extended to a broad spectrum of nonlinear adsorption isotherms. The results reveal a 56.5% reduction in solute breakthrough time under the Freundlich isotherm and a remarkable 2.6‐fold extension under the Langmuir isotherm when compared to the linear isotherm. The adsorption constants of the Freundlich and Langmuir isotherms exhibit a positive correlation with breakthrough time, while the exponent of the Freundlich isotherm and the maximal adsorption capacity in the Langmuir isotherm demonstrate a negative association with breakthrough time. This study enhances the precision of solute transport prediction and provides a more scientific assessment of clay liner performance.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142236179","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 on the Instability Mechanism of Expansive Soil Slope With Weak Interlayer Based on Strain Softening","authors":"Shuai Xu, Hanjing Jiang, Yongfu Xu, Aoxun Wang, Shunchao Qi","doi":"10.1002/nag.3834","DOIUrl":"https://doi.org/10.1002/nag.3834","url":null,"abstract":"Expansive soils are widespread in the world and coincide with areas of high human activity. The main cause of deep instability of expansive soil slopes is due to their softening caused by excavation and seepage. By developing a comprehensive numerical model based on the theory of unsaturated soil, this study examines the characteristics of stress and displacement distribution of expansive soil slopes through hydraulic‐mechanical coupled numerical simulation. This study analyzes the evolution patterns of slopes with excavation unloading and seepage of water storage to reveal the mechanisms of deep‐seated instability of expansive soil slopes. The findings demonstrate that: The instability of expansive soil slopes begins at the foot of the slope and propagates along the interlayer, affecting the entire slope. Excavation leads to the softening of the expansive soil interlayer and the transfer of shear stress. During water storage, the weakening of the soil strength results in slope instability along the weak interlayer slip. Softening of the expansive soil interlayer facilitates the redistribution of shear forces in the slope and alters the distribution law of the plastic zone in the deep layer. Overly slowing down the slope leads to significant excavation unloading, which is detrimental to the slope's stability.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142236219","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":"Cover Image, Volume 48, Issue 14","authors":"Kehao Chen, Rui Pang, Bin Xu, Xingliang Wang","doi":"10.1002/nag.3840","DOIUrl":"https://doi.org/10.1002/nag.3840","url":null,"abstract":"<p>The cover image is based on the article <i>Elastoplastic constitutive model for overconsolidated clays with an advanced dilatancy relation</i> by Kehao Chen et al., https://doi.org/10.1002/nag.3803.\u0000\u0000 <figure>\u0000 <div><picture>\u0000 <source></source></picture><p></p>\u0000 </div>\u0000 </figure>\u0000 </p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3840","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142165751","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":"Extended B-spline-based implicit material point method for saturated porous media","authors":"Yuya Yamaguchi, Shuji Moriguchi, Kenjiro Terada","doi":"10.1002/nag.3827","DOIUrl":"10.1002/nag.3827","url":null,"abstract":"<p>The large deformation and fluidization process of a solid–fluid mixture includes significant changes to the temporal scale of the phenomena and the shape and properties of the mixed material. This paper presents an extended B-spline (EBS)-based implicit material point method (EBS-MPM) for the coupled hydromechanical analysis of saturated porous media to enhance the overall versatility of MPM in addressing such diverse phenomena. The proposed method accurately represents phenomena such as high-speed motion in both the quasi-static and dynamic states by employing a full formulation of coupled hydromechanical modeling. The weak imposition of boundary conditions based on Nitsche's method allows representing the boundary conditions independent of the relative position of the particles and computational grid. In addition, it enables dynamic changes in the boundary domain based on the deformation. The robustness of this boundary representation is reinforced using EBS basis functions, which prevent the degradation of the condition number of the system matrices regardless of the position of the boundary domain with respect to the computational grid. Furthermore, a stabilization method based on a variational multiscale method (VMS) approach is employed to provide the flexibility in choosing arbitrary basis functions for spatial discretization, facilitating the effective construction of EBS. Numerical examples including comparisons between a full formulation and a simplified formulation are presented to demonstrate the performance of the developed method under various boundary conditions and loading states across different time scales.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":3.4,"publicationDate":"2024-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3827","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142131006","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":"Development of improved finite element formulations for pile group behavior analysis under cyclic loading","authors":"Jian‐Hong Wan, Shui‐Hua Jiang, Xue‐You Li, Zhilu Chang","doi":"10.1002/nag.3828","DOIUrl":"https://doi.org/10.1002/nag.3828","url":null,"abstract":"The effect of cyclic loading is an essential factor leading to progressive soil strength degradation. Therefore, a comprehensive analysis of the pile‐soil system behavior under cyclic loading is required to ensure the stability of pile group. There is room for improvement in the inherent constraint of the conventional numerical model in terms of approximating the soil resistance distribution along the pile by point loads at element nodes, necessitating a specific element that integrates considerations of pile group effect and cyclic loading within a unified framework. This study aims to develop a newly specific type of element for efficiently predicting nonlinear behavior within the pile‐soil system, addressing simulations involving nonlinear pile‐soil interaction, pile group effect, and cyclic loading. Modified element formulations based on soil stiffness matrices and soil resistance vectors specifically address pile group effect and consider parameters that influence pile behavior under cyclic lateral loading. The numerical solution procedure with Newton‐Raphson iteration allows the calculation of pile responses in geometric and material nonlinear analyses. The validation of the proposed method includes several examples, comparing it with existing numerical solutions and experimental tests of single piles and pile groups under cyclic loading. These comparisons further support the consistency of the proposed method with measured data and validate its accuracy in considering group effect and cyclic loading. The parametric study illustrates the ability of the proposed method to capture cyclic loading parameters while considering the influence of the number and magnitude of load cycles, the cyclic load direction, and the installation methods.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142131007","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}