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

{"title":"Poromechanical Solution for One‐Dimensional Large Strain Consolidation of Modified Cam‐Clay Soil","authors":"Sheng‐Li Chen, Hai‐Sui Yu, Younane N. Abousleiman, Christopher E. Kees","doi":"10.1002/nag.3924","DOIUrl":"https://doi.org/10.1002/nag.3924","url":null,"abstract":"A theoretical model describing the one‐dimensional large strain consolidation of the modified Cam‐Clay soil is presented in this paper. The model is based on the Lagrangian formulation and is capable of featuring the variability of soil compressibility (inherently due to the direct incorporation of the specific Cam‐Clay plasticity model) and permeability, as well as the impact of the overconsolidation ratio (OCR). The derivation starts from the establishment of the incremental stress–strain relations for both purely elastic and elastoplastic deformations under one‐dimensional compression conditions, and thereafter the coefficients of compressibility/volume change that are essential to the consolidation analysis. The governing partial differential equation is then neatly deduced in conjunction with the continuity and equilibrium conditions for the soil, with the vertical effective stress being the privileged unknown to be solved for. Subsequently, a semi‐analytical solution to the developed rigorous poroelastoplastic large strain consolidation model is obtained and verified with the ABAQUS finite element numerical results. Parametric analyses are finally provided to investigate in detail the influences of the soil overconsolidation ratio, large strain configuration, and the variability of the soil permeability on the calculated one‐dimensional consolidation response.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"56 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142934885","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":"Stress Distribution Around Arbitrarily Shaped Shallow Buried Tunnels in Transversely Isotropic Rock Mass","authors":"Zhi Yong Ai, Qi Liang, Zi Kun Ye, Ke Xin Hu","doi":"10.1002/nag.3936","DOIUrl":"https://doi.org/10.1002/nag.3936","url":null,"abstract":"In this paper, the stress solutions around shallow buried tunnels with arbitrary shapes in transversely isotropic rock mass are derived. First, the hybrid penalty function method is employed to derive the mapping function from an arbitrarily shaped tunnel to a unit circle in the complex plane. Then the complex function method and the Schwartz alternating method are used to derive the solution of the studied problem. A MATLAB program is developed on the basis of the presented theory. Thereafter, the finite element results of ABAQUS are compared with those obtained by the presented solutions, which confirms the theoretical and computational accuracy of the proposed approach, and the impacts of transverse isotropy, tunnel depth, and tunnel shape on the surrounding rock stresses are discussed.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"20 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142924996","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 Porothermoelastodynamic Modeling of Stress Wave Through a Fluid‐Saturated Porous Cylinder","authors":"Chao Liu","doi":"10.1002/nag.3934","DOIUrl":"https://doi.org/10.1002/nag.3934","url":null,"abstract":"Analytical porothermoelastodynamic (PTED) solutions are rare in the literature. The responses of fluid‐saturated porous materials subject to coupled mechanisms of loading frequency, fluid flow, stress, and temperature are unclear. In this paper, we use the PTED theory and derive the analytical solutions of pore pressure, temperature, stress, force, and displacement for an isotropic fluid‐saturated porous cylinder subject to a harmonic vibration. The coupled partial differential equations among pore pressure, displacement, and temperature are decoupled by matrix diagonalization and solved by further introducing a potential function and separation of variables. The PTED solution reproduces the poroelastodynamic (PED) one by easing the thermal effect. A demonstration example shows that the coupled mechanisms among pore pressure, stress, and displacement are highly frequency dependent. The thermal effect is more pronounced at low frequencies than at high frequencies when the inertial impact is more significant. Pore pressure is almost uniform in the ‐direction at low frequencies and becomes nonuniform at high frequencies for both PTED and PED cases. Displacements exhibit linear behavior at low frequencies and become nonlinear at high frequencies. Thermal stress and expansion significantly impact the pore pressure and displacement. A brief sensitivity analysis shows that pore pressure responds linearly and monotonically with the increase of the volumetric thermal expansion coefficient of the solid matrix at low frequencies and becomes nonlinear and nonmonotonic at high frequencies. The volumetric thermal expansion coefficient of the solid matrix has a minor effect on the vertical displacement and significantly influences the radial displacement.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"46 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142908282","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 Fatigue Cohesive Law‐Embedded Finite‐Discrete Element Method for Pulsed Hydraulic Fracture Simulation","authors":"Xuanchun Wei, Lei Wang, Yancao Li, Jingtao Ding, Zhennan Zhang","doi":"10.1002/nag.3935","DOIUrl":"https://doi.org/10.1002/nag.3935","url":null,"abstract":"The pulsed hydraulic fracture (PHF) is a stimulation technique of reservoir, which can lower breakdown pressure by generating fatigue fracture. In this study, a fatigue cohesive law is proposed and embedded into the finite‐discrete element method (FDEM) to describe the fatigue failure of the interface under cyclic loading. The fatigue is assumed to result from the accumulation of plastic deformation, whose increment is related to the traction variation range of each cycle and the total plastic deformation. The proposed fatigue cohesive law is validated by the uniaxial compression and mixed‐mode three‐point bending test simulation. Then this fatigue cohesive law is embedded into a fully hydraulic–mechanical coupled FDEM to simulate the PHF. The influence of loading scheme, flow rate, frequency, viscosity and natural fracture density on PHF behaviors is discussed. The results suggest that both the pressure‐ and injection rate–controlled PHF can reduce the breakdown pressure. The flow rate, frequency, and viscosity have a great impact on the performance of PHF. The natural fractures surrounding a hydraulic fracture (HF) can be gradually activated under cyclic injection. The activated natural fractures contribute to the complicated HF network. These results indicate that the proposed fatigue cohesive model can effectively simulate the fatigue fracture of rock and HF propagation under cyclic injection.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"6 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142908263","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":"Phase-Field Modeling of Fracture Under Compression and Confinement in Anisotropic Geomaterials","authors":"Maryam Hakimzadeh,&nbsp;Carlos Mora-Corral,&nbsp;Noel Walkington,&nbsp;Giuseppe Buscarnera,&nbsp;Kaushik Dayal","doi":"10.1002/nag.3933","DOIUrl":"10.1002/nag.3933","url":null,"abstract":"<p>Strongly anisotropic geomaterials, such as layered shales, have been observed to undergo fracture under compressive loading. This paper applies a phase-field fracture model to study this fracture process. While phase-field fracture models have several advantages—primarily that the fracture path is not predetermined but arises naturally from the evolution of a smooth non-singular damage field—they provide unphysical predictions when the stress state is complex and includes compression that can cause crack faces to contact.</p><p>Building on a recently developed phase-field model that accounts for compressive traction across the crack face, this paper extends the model to the setting of anisotropic fracture. The key features of the model include the following: (1) a homogenized anisotropic elastic response and strongly anisotropic model for the work to fracture; (2) an effective damage response that accounts consistently for compressive traction across the crack face, that is derived from the anisotropic elastic response; (3) a regularized crack normal field that overcomes the shortcomings of the isotropic setting, and enables the correct crack response, both across and transverse to the crack face.</p><p>To test the model, we first compare the predictions to phase-field fracture evolution calculations in a fully resolved layered specimen with spatial inhomogeneity, and show that it captures the overall patterns of crack growth. We then apply the model to previously reported experimental observations of fracture evolution in laboratory specimens of shales under compression with confinement, and find that it predicts well the observed crack patterns in a broad range of loading conditions. We further apply the model to predict the growth of wing cracks under compression and confinement. Prior approaches to simulate wing cracks have treated the initial cracks as an external boundary, which makes them difficult to apply to general settings. Here, the effective crack response model enables us to treat the initial crack simply as a nonsingular damaged zone within the computational domain, thereby allowing for easy and general computations.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1319-1335"},"PeriodicalIF":3.4,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3933","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142908283","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":"Laboratory-Scale Free Fall Cone Penetrometer Test on Marine Clay: A Numerical Investigation Using the Generalized Interpolation Material Point Method","authors":"Debasis Mohapatra,&nbsp;Saeideh Mohammadi,&nbsp;Maarit Saresma,&nbsp;Joonas J. Virtasalo,&nbsp;Wojciech T. Sołowski","doi":"10.1002/nag.3929","DOIUrl":"10.1002/nag.3929","url":null,"abstract":"<p>This paper presents a series of laboratory free-fall cone penetrometer (FFCP) tests conducted on marine clay samples collected from the Gulf of Finland in the Baltic Sea. Subsequently, these tests are replicated numerically with the generalized interpolation material point method (GIMP) simulations. First, the paper gives laboratory-scale FFCP experiment results used for the validation of the numerical framework. In these experiments, a small-scale model of a FFCP was dropped from various heights into a natural marine clay soil sample and recorded using a high-speed camera. The tests were supplemented with a laboratory test program to determine the geotechnical properties of the clay used in the experiments. Following image processing, the tests provided data for numerical simulations: displacement, velocity, acceleration, and reaction force curves associated with the FFCP during the penetration process. The GIMP simulations shown in the paper replicate the process of penetration of the FFCP into the marine clay. The simulations used a strain-rate dependent Tresca constitutive model, extended with strain softening that replicates the reduction of the undrained shear strength due to destructuration, an important feature of the material. The numerical simulations replicate the experiments well. The study examines the effect of cone penetrometer roughness, impact velocity, mesh density, strain rate, and strain softening on the cone penetrometer penetration process. The simulation results indicate that the presented framework can replicate the dynamic penetration process on soft and sensitive clay very well.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1299-1318"},"PeriodicalIF":3.4,"publicationDate":"2024-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3929","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142908281","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":"Computational Large-Deformation-Plasticity Periporomechanics for Localization and Instability in Deformable Porous Media","authors":"Xiaoyu Song,&nbsp;Hossein Pashazad,&nbsp;Andrew J. Whittle","doi":"10.1002/nag.3920","DOIUrl":"10.1002/nag.3920","url":null,"abstract":"<div>\u0000 \u0000 <p>In this article, we formulate a computational large-deformation-plasticity (LDP) periporomechanics (PPM) paradigm through a multiplicative decomposition of the deformation gradient following the notion of an intermediate stress-free configuration. PPM is a nonlocal meshless formulation of poromechanics for deformable porous media through integral equations in which a porous material is represented by mixed material points with nonlocal poromechanical interactions. Advanced constitutive models can be readily integrated within the PPM framework. In this paper, we implement a linearly elastoplastic model with Drucker–Prager yield and post-peak strain softening (loss of cohesion). This is accomplished using the multiplicative decomposition of the nonlocal deformation gradient and the return mapping algorithm for LDP. The paper presents a series of numerical examples that illustrate the capabilities of PPM to simulate the development of shear bands, large plastic deformations, and progressive slope failure mechanisms. We also demonstrate that the PPM results are robust and stable to the material point density (grid spacing). We illustrate the complex retrogressive failure observed in sensitive St. Monique clay that was triggered by toe erosion. The PPM analysis captures the distribution of horst and graben structures that were observed in the failed clay mass.</p></div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1278-1298"},"PeriodicalIF":3.4,"publicationDate":"2024-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142904990","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":"Multiscale Stochastic Modeling of Backward Erosion Piping Initiation, From Grain Kinetics to Weibull Statistics. Part I: Analytical Derivations","authors":"Zhijie Wang,&nbsp;Caglar Oskay,&nbsp;Alessandro Fascetti","doi":"10.1002/nag.3931","DOIUrl":"10.1002/nag.3931","url":null,"abstract":"<p>Backward erosion piping (BEP) is a significant contributor to failures in global flood protection infrastructure, yet it remains among the least understood geotechnical phenomena, particularly concerning the fundamental mechanisms driving its initiation. This study focuses on the development of a novel stochastic framework for the prediction of critical hydraulic gradients causing BEP initiation. The novelty of the study lies in the following: (1) the development of a grain-scale probabilistic model based on fundamental mechanisms by means of the theory of rate processes, (2) quantification of the influence of soil variability on BEP initiation probability by introducing an initiation probability function, and (3) an analytical framework reconciling grain kinetics of BEP initiation with the Weibull distribution. A particle-scale BEP initiation probabilistic model is first established based on fundamental grain kinetics under seepage flow by using the theory of rate processes. To investigate how soil variability influences initiation, a stochastic dual random lattice modeling framework is exercised, complemented by direct x-ray computed tomography measurements of soil variability conducted on sand samples. The analytical probabilistic model for BEP initiation closely aligns with the Weibull distribution, also demonstrating that soil variability influences both the scale and shape parameters of the distribution. This work establishes the linkage between probability of BEP initiation as described by the theory of rate processes and phenomenological Weibull statistics. Findings presented herein bring the potential to develop a multiscale probabilistic framework by means of Weibull statistics for evaluating the probability of BEP initiation at multiple scales.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1262-1277"},"PeriodicalIF":3.4,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3931","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142888310","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":"Multiscale Stochastic Modeling of Backward Erosion Piping Initiation, From Grain Kinetics to Weibull Statistics. Part II: Model Validation and Applications","authors":"Zhijie Wang,&nbsp;Caglar Oskay,&nbsp;Alessandro Fascetti","doi":"10.1002/nag.3930","DOIUrl":"10.1002/nag.3930","url":null,"abstract":"<p>Backward erosion piping (BEP) is a leading internal erosion mechanism for flood protection system failures. A model capable of predicting critical hydraulic conditions for BEP initiation at multiple scales while also incorporating soil variability is a pressing need. This study formulates and validates a novel multiscale probabilistic BEP initiation framework with incorporation of soil variability. The framework is based on a grain-scale probabilistic model and the weakest link theory, and the theory of rate processes. The multiscale framework proposed herein is validated through a wide range of available experimental data from independent sources, encompassing tests performed at multiple scales. Following calibration with small-scale experimental data, the model demonstrates accurate prediction of critical hydraulic gradients at larger scales (3–6 orders of magnitude difference), including the ability to capture the grain size dependence of BEP initiation and providing uncertainty estimates. A systematic analysis is performed to uncover the effects of different soil properties on multiscale critical hydraulic conditions.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1247-1261"},"PeriodicalIF":3.4,"publicationDate":"2024-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3930","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142888188","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":"Analysis and Optimisation of Obstacle-Crossing Performance of Electric Shovel Based on DEM-MBD Coupling Method","authors":"Zeren Chen,&nbsp;Wei Guan,&nbsp;Ruibin Li,&nbsp;Guang Li,&nbsp;Duomei Xue,&nbsp;Zhengbin Liu,&nbsp;Guoqiang Wang","doi":"10.1002/nag.3927","DOIUrl":"10.1002/nag.3927","url":null,"abstract":"<div>\u0000 \u0000 <p>To study and enhance the obstacle-crossing performance of the electric shovel, an obstacle-crossing model that employs a coupling methodology integrating the discrete element method (DEM) and multi-body dynamics (MBD) is constructed. Secondly, the influence of grouser height (GH), track velocity (TV), slope inclination (SI) and slope height (SH) on obstacle-crossing performance is investigated through DEM-MBD simulation, with the objective of obtaining an obstacle-crossing surrogate model through the Kriging method and Box-Behnken experimental design. On this basis, two optimisation solutions for the obstacle-crossing performance of the electric shovel are proposed based on a genetic algorithm (GA), and the corresponding obstacle-crossing performances are analysed. The results demonstrate that the coupling effect between SI and SH exerts a considerable influence on the ground pressure coefficient (GPC), power and disturbance potential energy (DPE). When the optimal TV and GH are set at 0.1 m/s and 9.38 mm, the GPC, power and disturbance kinetic energy (DKE) are observed to diminish to varying degrees, thereby indicating that the obstacle-crossing performance of the electric shovel has been enhanced.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"49 4","pages":"1232-1246"},"PeriodicalIF":3.4,"publicationDate":"2024-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142886728","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}