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

{"title":"Use of Advanced Constitutive Models for the Mechanical Behavior of Soft Soils With Diatoms From Bogotá (Colombia)","authors":"Cristhian Mendoza, Márcio Muniz de Farias","doi":"10.1002/nag.3863","DOIUrl":"https://doi.org/10.1002/nag.3863","url":null,"abstract":"Most constitutive models did not initially consider special behaviors in some soils with singular characteristics (e.g., soft soils with diatom content). For example, at first, these models did not consider the effect of soil structure and viscosity. However, in the last decades, these variables have been incorporated into several constitutive models to describe the mechanical behavior of the soil in its natural state. Structure and viscosity laws that adequately reproduce the soil behavior had to be developed to include these variables. This paper compares the mechanical behavior of soft soils in Bogotá with different constitutive models. Bogotá’s soft soils are lacustrine deposits with a high content of diatoms in their structure. Natural soil samples with intact structures show a high‐water content, which can be higher than 300%, liquid limits of up to 400%, void ratios higher than five, and friction angles of almost 40°. In addition, the model validations were made through the simulations of triaxial tests in compression and shear paths. Modified Cam Clay (MCC), hypoplastic (HP), and subloading Cam Clay (SCC) were the constitutive models used. Two models are based on an elastoplastic framework, and the third uses a HP framework. Several lessons were learned from the simulations regarding the strengths and weaknesses of the models compared to the tests carried out. Finally, the extensive discussion revolves around determining the most suitable model for simulating the mechanical behavior of soft soils containing diatoms in Bogotá.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"24 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142397744","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 Semi‐Analytical Method for Simulating Near‐Field Antiplane Wave Propagation in Layered Fluid‐Saturated Porous Media","authors":"Liang Li, Man Wang, Hongyun Jiao, Xiuli Du, Peixin Shi","doi":"10.1002/nag.3859","DOIUrl":"https://doi.org/10.1002/nag.3859","url":null,"abstract":"A semi‐analytical method for the near‐field antiplane wave propagation analysis in the layered fluid‐saturated porous media (FSPM) is proposed based on the Biot <jats:italic>u</jats:italic>–<jats:italic>U</jats:italic> dynamic formulation. The wave propagation equations of the FSPM are decoupled by the variable‐separating method. The thin‐layer element method (TLEM) is applied to discretize the infinite domain and construct the consistent artificial boundary condition. The finite element method (FEM) is adopted for the space discretization of the finite domain and the numerical solution of the dynamic response. The proposed method is validated by the comparison of the numerical results of this method with those in the published references and acquired from the remote artificial boundary. Subsequently, this method is applied to investigate typical near‐field antiplane wave propagation problems in the FSPM. Parametric sensitivity investigations are also executed to explore the impact of mechanical parameters, including permeability coefficients, porosity, and shear modulus of the solid phase, on the dynamic response of the FSPM. The study results confirm the efficacy and efficiency of the proposed method in the near‐field antiplane wave propagation analysis in the FSPM.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"13 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142385626","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 Consolidation of Soft Ground With Impervious Columns Considering Non‐Darcian Flow","authors":"Kuo Li, Mengmeng Lu, Jinxin Sun","doi":"10.1002/nag.3857","DOIUrl":"https://doi.org/10.1002/nag.3857","url":null,"abstract":"Cohesive material columns have been extensively used in foundation improvement projects to enhance the bearing capacity of composite foundations and mitigate post‐construction settlement. However, the low permeability of cohesive material columns restricts the dissipation of pore water primarily through the top surface of the foundation, potentially resulting in longer drainage paths compared to foundations treated with granular material columns or vertical drains. Moreover, the impact of non‐Darcian flow within soils on consolidation behavior becomes increasingly pronounced as the drainage path increases. Consequently, a novel analytical model for the consolidation of impervious column‐assisted foundations is established, which can incorporate the seepage model accounting for the initial hydraulic gradient. The accuracy and reasonableness of the obtained solution are then validated by conducting a comparative analysis with existing models and through a detailed case study. Furthermore, a parametric analysis is conducted to delve into the influence of several crucial factors on the consolidation performance. The findings demonstrate that non‐Darcian flow has a greater influence on composite foundations compared to natural foundations. Additionally, the threshold value of the well‐diameter ratio decreases with the increase in the initial hydraulic gradient. Finally, the final seepage front remains at a shallower position when the column–soil modulus ratio becomes larger, and the influence of non‐Darcian flow on the consolidation rate becomes more pronounced.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"27 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142385428","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":"Evaluation of the Spatial Variability of the Mechanical Properties of Rocks Using Non‐Iterative Green's Function Approach and the FOSM Method","authors":"Leonardo C. Mesquita, Elisa D. Sotelino, Matheus L. Peres","doi":"10.1002/nag.3861","DOIUrl":"https://doi.org/10.1002/nag.3861","url":null,"abstract":"The present work proposes a new version of the Green‐FOSM (first‐order second moment) method, which eliminates the iterative calculation process of the original version and, simultaneously, solves the convergence problems related to the mechanical properties of rocks that form the geological formation. In this calculation scheme, the iterative process is eliminated by using a matrix that correlates the nodal displacement vector with the strain vector. Considering the same computational resources, this non‐iterative version of the Green‐FOSM method is up to 200 times faster than the original iterative process. In addition, it allows analyzing problems with more than 10,000 random variables, value that in the original method is less than 3000. To demonstrate its validity, the proposed method is applied to two hypothetical models subjected to different fluid extraction processes. For all the different levels of correlation and spatial variability, the statistical results obtained by the proposed method agree well with the results obtained via Monte Carlo Simulation (MCS). The relationship between CPU times demonstrates that the proposed method is at least 50 times faster than MCS. In the end, the non‐iterative Green‐FOSM method is used to obtain the displacement, strain, and stress fields of a geological section constructed from a seismic image of Brazilian pre‐salt oil region. The results found show that, depending on the levels of spatial variability, the analyzed fields can assume values up to 30.6% higher or lower than the values obtained deterministically.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"29 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2024-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142384212","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":"Efficient Random Field Generation With Rotational Anisotropy for Probabilistic SPH Analysis of Slope Failure","authors":"Zhonghui Bi,&nbsp;Wei Wu,&nbsp;Liaojun Zhang,&nbsp;Chong Peng","doi":"10.1002/nag.3858","DOIUrl":"10.1002/nag.3858","url":null,"abstract":"<p>Due to geological processes such as sedimentation, tectonic movement, and backfilling, natural soil often exhibits characteristics of rotated anisotropy. Recent studies have shown the significant impact of rotated anisotropy on slope stability. However, little research has explored how this rotated anisotropy affects the large deformations occurring after slope failure. Therefore, this study integrates rotated random field theory with smoothed particle hydrodynamics (SPH) to investigate its influence on post-failure slope behavior. Focusing on a typical slope scenario, this research utilizes graphics processing unit (GPU)–accelerated covariance matrix decomposition (CMD) method to create rotated anisotropy random fields and applies the SPH framework for analysis. It examines the influence of rotated anisotropy angles and the cross-correlation between cohesion and internal friction angle on landslides. The results indicate that the rotational anisotropy of the slope significantly influences post-failure behavior. When the rotation angle is close to the slope surface, it tends to amplify both the magnitude and variability of slope failure. Furthermore, the study evaluates the efficiency of generating these random fields and emphasizes the substantial computational speed improvements achieved with GPU acceleration. These findings offer a robust approach for probabilistic analysis of slope large deformations considering rotated anisotropy. They provide a theoretical foundation for accurately assessing the risk of slope collapse, holding significant practical implications for geotechnical engineering.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4520-4538"},"PeriodicalIF":3.4,"publicationDate":"2024-10-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3858","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142384210","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":"Implementation of the Zienkiewicz–Pande Model into a Four-Dimensional Lattice Spring Model for Plasticity and Fracture","authors":"Xin-Dong Wei,&nbsp;Zhe Li,&nbsp;Gao-Feng Zhao","doi":"10.1002/nag.3860","DOIUrl":"10.1002/nag.3860","url":null,"abstract":"<div>\u0000 \u0000 <p>Plasticity and fracture problems have always been hot topics in numerical methods. In this work, a universal implementation procedure for the elasto-plastic constitutive model is developed in the four-dimensional lattice spring model (4D-LSM), in which the Jaumann stress rate is incorporated to exclude the influence of the rigid rotation in the particle stress, expanding the ability of 4D-LSM to deal with large elastic deformation problems by its own to large plastic deformation problems. As an example, the Zienkiewicz–Pande (ZP) constitutive model is implemented. Several numerical examples are carried out to check the performance of the implemented model. Through a comparison with analytical solutions, available experimental data, and other numerical results, the stability of the developed plastic framework and the correctness of the stress calculation scheme are verified. Meanwhile, numerical results show that the developed code is capable of solving elasto-plastic large deformation problems. With the advantage of 4D-LSM in handling fracture problems, the ability of the embedded model to solve plastic fracture problems is verified with a simple maximum deformation failure criterion.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4506-4519"},"PeriodicalIF":3.4,"publicationDate":"2024-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142369109","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":"3D Semi-Analytical and Numerical Upper-Bound Homogenisation Approaches to the Out-of-Plane Strength Domain of a Running-Bond Masonry Wall","authors":"Elodie Donval,&nbsp;Ghazi Hassen,&nbsp;Duc Toan Pham,&nbsp;Patrick de Buhan,&nbsp;Martin Vigroux","doi":"10.1002/nag.3851","DOIUrl":"10.1002/nag.3851","url":null,"abstract":"<p>The present contribution proposes a new semi-analytical homogenisation approach to determine a running-bond masonry wall's in- and out-of-plane strength domain based on the yield design framework. The main novelty of such an approach is that it does not rely on simplifying assumptions such as infinitely thin joints or plane stress state within the blocks, by making use of 3D virtual failure mechanisms in the kinematic approach. The new semi-analytical approach is then compared to a state-of-the-art numerical implementation of the kinematic approach of yield design, relying on semi-definite programming. Several comparisons show a good agreement between the semi-analytical and the numerical approaches and outline the computational efficiency of the semi-analytical approach as well as the fact that it is very well suited for engineering design purposes. Both proposed approaches are then compared to existing approaches based on the limit analysis or yield design framework.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4469-4490"},"PeriodicalIF":3.4,"publicationDate":"2024-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3851","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142369084","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":"Microscopic Thermo-Mechanical Properties and Phase Transition of Bulk Ice-Ih","authors":"Pengchang Wei,&nbsp;Weiwei Niu,&nbsp;Chi Yao,&nbsp;Zhenyu He,&nbsp;Yuan-Yuan Zheng,&nbsp;Wei Ma","doi":"10.1002/nag.3856","DOIUrl":"10.1002/nag.3856","url":null,"abstract":"<div>\u0000 \u0000 <p>The ice–water phase transition of bulk ice could develop with varying temperatures and external loads, significantly affecting its mechanical properties. The coupling effect of temperature and shear loads on the thermo-mechanical properties of bulk ice and its phase transition evolution is poorly understood, especially at the nanoscale. In this study, molecular dynamics (MD) simulation method was employed to investigate the thermo-mechanical behaviours of bulk ice-Ih system at the microscale under various temperatures (73–270 K) and shear paths, where its phase transition, elastic properties, structure deformation mechanism and structural anisotropy were discussed. The simulation results show that (1) the shear modulus, shear strength and ultimate shear strain of bulk ice-Ih system could linearly decrease with rising temperature, aligning with previous studies. (2) Two types of failure modes from bulk ice-Ih system were founded, such as solid–liquid phase co-existence at 73–225 K and liquid phase at 250–270 K. (3) Ice melting into water was attributed to the fracture of hydrogen bond during shear process. (4) Compared to vertical shearing (<i>XZ</i> (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>11</mn>\u0000 <mover>\u0000 <mn>2</mn>\u0000 <mo>¯</mo>\u0000 </mover>\u0000 <mn>0</mn>\u0000 </mrow>\u0000 <annotation>$11bar{2}0$</annotation>\u0000 </semantics></math>) and <i>YZ</i> (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>01</mn>\u0000 <mover>\u0000 <mn>1</mn>\u0000 <mo>¯</mo>\u0000 </mover>\u0000 <mn>0</mn>\u0000 </mrow>\u0000 <annotation>$01bar{1}0$</annotation>\u0000 </semantics></math>)) directions, the mechanical response along the horizontal shearing (<i>XY</i> (0001)) direction was most sensitive to temperature effect.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4491-4505"},"PeriodicalIF":3.4,"publicationDate":"2024-10-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142374122","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 2D Electro-Osmotic Consolidation of Unsaturated Soil With Non-linear Voltage Distribution","authors":"Xudong Zhao,&nbsp;Jie Min,&nbsp;Shaolin Ding,&nbsp;Yang Liu,&nbsp;Jiaxin Liao,&nbsp;Shuai Zhang","doi":"10.1002/nag.3854","DOIUrl":"10.1002/nag.3854","url":null,"abstract":"<div>\u0000 \u0000 <p>Existing solutions for electro-osmotic consolidation assume a linear voltage distribution, which is inconsistent with the experimental findings. The present study introduces a novel two-dimensional electro-osmotic consolidation model for unsaturated soils, which considers the influence of non-linear voltage distribution. The closed-form solution is derived by employing the eigenfunction expansion method and the Laplace transform technique. The accuracy of the analytical solutions is validated through the implementation of finite element simulations. The findings from the parametric studies indicate that the excess pore water pressure (EPWP) observed in electro-osmotic consolidation is influenced by the distribution of voltage. The dissipation rate of EPWP is observed to be higher when subjected to non-linear voltage conditions compared to linear voltage conditions. Moreover, the impact of non-linear voltage distribution becomes more pronounced in unsaturated soil characterised by higher electro-osmosis conductivity and a lower ratio of <i>k_x/k_y</i>. In contrast, the excess pore air pressure (EPAP) remains unaffected by the voltage distribution.</p>\u0000 </div>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4457-4468"},"PeriodicalIF":3.4,"publicationDate":"2024-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142329002","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 Parameter Identification for a High-Cycle Accumulation Model of Sand With Enhancement of Cuckoo Search Algorithm","authors":"Shao-Heng He,&nbsp;Zhen-Yu Yin,&nbsp;Yifei Sun,&nbsp;Zhi Ding","doi":"10.1002/nag.3838","DOIUrl":"10.1002/nag.3838","url":null,"abstract":"<p>This study presents a novel approach of intelligent parameter identification (IPI) for a high-cycle accumulation (HCA) model of sand, which reduces the subjective errors on manual parameter calibration and makes the use of the HCA model more accessible. The technique is based on optimization theory and adopts the cuckoo search algorithm (CSA). To improve search ability and convergence speed of CSA, several enhancements are implemented. First, the improved CSA (ICSA) incorporates quasi-opposition learning to expand the search space and replaces the original search strategy with a Cauchy random walk to enhance global search ability. Second, an adaptive scaling factor is introduced in the algorithm's control parameters to achieve a better balance between exploration speed and accuracy. Third, a dynamic inertia weight is used to balance the search between global and local spaces when generating new nest positions after abandoning old ones. The performance of the ICSA-based IPI approach is evaluated by comparing it with the original CSA-based IPI and manual calibration in determining the HCA model parameters. A comprehensive analysis is also conducted to assess the effectiveness and superiority of each improvement strategy introduced in the ICSA over the original CSA. All comparisons demonstrate that the proposed ICSA-based IPI method is more powerful and efficient in finding optimal parameters.</p>","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"48 18","pages":"4410-4427"},"PeriodicalIF":3.4,"publicationDate":"2024-09-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/nag.3838","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142328665","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}