Suifeng Wang, Hideaki Yasuhara, Li Zhuang, Xianyu Zhao, Liping Zhang, Tao Wang
{"title":"A Novel Hydro‐Grain‐Texture Model to Unveil the Impact of Mineral Grain Anisotropy on Fluid‐Driven Cracking Processes in Crystalline Rock","authors":"Suifeng Wang, Hideaki Yasuhara, Li Zhuang, Xianyu Zhao, Liping Zhang, Tao Wang","doi":"10.1002/nag.3888","DOIUrl":"https://doi.org/10.1002/nag.3888","url":null,"abstract":"The anisotropy at the grain scale significantly impacts cracking behavior of crystalline rocks. However, the anisotropy of mineral structure, especially the grain shape and orientation has been inadequately addressed in studies on hydraulic fracturing. To bridge this gap, this paper introduces a coupled hydro‐grain‐texture model (HGTM) based on discrete element model (DEM) that investigates the influence of grain shape and orientation on fluid‐driven cracking processes in crystalline rock. The HGTM can consider the different mineral grain shapes and orientations by changing the aspect ratio and rotating coordinate axes. Our studies covered six distinct in‐situ stresses, three grain shapes, and five grain orientations. Initially, we present a comprehensive examination of the microcracking processes of hydraulic fracturing. Then the influences of in‐situ stress, grain shape, and grain orientation on cracking processes were studied. The results underscore that both mineral grain and in‐situ stress interplay to influence the hydraulic fracturing of the crystalline rocks. The proposed HGTM can well mimic the propagation process of hydraulic fracturing by comparing with the experimental results and the results reveal that hydraulic fracturing in crystalline rocks is a highly complex process. This research clarifies the complex interplay between grain texture and hydraulic fracturing, offering invaluable insights for optimizing stimulation practices.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142597495","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":"Formulation and Implicit Numerical Integration of a Kinematic Hardening Model for Unsaturated Soils","authors":"Lluís Monforte, Mohamed Rouainia","doi":"10.1002/nag.3878","DOIUrl":"https://doi.org/10.1002/nag.3878","url":null,"abstract":"Currently, our understanding of material‐scale deterioration resulting from meteorologically induced variations in pore water pressure and its significant impact on infrastructure slopes is limited. To bridge this knowledge gap, we have developed an extended kinematic hardening constitutive model for unsaturated soils that refines our understanding of weather‐driven deterioration mechanisms in heterogeneous clay soils. This model has the capability of predicting the irrecoverable degradation of strength and stiffness that has been shown to occur when soils undergo wetting and drying cycles. The model is equipped with a fully coupled and hysteretic water retention curve and a hysteretic loading–collapse curve and has the capability to predict the irrecoverable degradation of strength and stiffness that occurs during cyclic loading of soils. Here, we employ a fully implicit stress integration technique and give particular emphasis to deriving a consistent tangent operator, which includes the linearisation of the retention curve. The proposed algorithm is evaluated for efficiency and performance by simulating various stress and strain‐driven triaxial paths, and the accuracy of the integration technique is evaluated through the use of convergence curves.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142597485","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}
Xiao‐Hong Wang, Meng‐Chen Yue, Zhi‐Feng Liu, Wei‐Dong Cao, Yong Wang, Jun Hu, Chang‐Hao Xiao, Yao‐Yong Li
{"title":"Numerical Strategy on the Grid Orientation Effect in the Simulation for Two‐Phase Flow in Porous Media by Using the Adaptive Artificial Viscosity Method","authors":"Xiao‐Hong Wang, Meng‐Chen Yue, Zhi‐Feng Liu, Wei‐Dong Cao, Yong Wang, Jun Hu, Chang‐Hao Xiao, Yao‐Yong Li","doi":"10.1002/nag.3886","DOIUrl":"https://doi.org/10.1002/nag.3886","url":null,"abstract":"In the context of numerical simulators for multiphase flow in porous media, there exists a long‐standing issue known as the grid orientation effect (GOE), wherein different numerical solutions can be obtained when considering grids with different orientations under certain unfavorable conditions. The GOE is relevant to the instability near displacement fronts. If numerical oscillations accompanied by sharp fronts are not adequately suppressed, the GOE occurs. To reduce or even eliminate the GOE, we propose augmenting adaptive artificial viscosity in the process of solving the saturation equation. It has been demonstrated that appropriate artificial viscosity can effectively reduce or even eliminate the GOE. The proposed numerical method can be easily applied in practical engineering problems.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142598204","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}
Vahab Sarfarazi, Hadi Haeri, Jinwei Fu, Amir Namdarmanesh, Gholamreza Saeedi, Naser Golsanami
{"title":"Experimental and Numerical Indentation Tests to Study the Crack Growth Properties Caused by Various Indenters Under High Temperature","authors":"Vahab Sarfarazi, Hadi Haeri, Jinwei Fu, Amir Namdarmanesh, Gholamreza Saeedi, Naser Golsanami","doi":"10.1002/nag.3885","DOIUrl":"https://doi.org/10.1002/nag.3885","url":null,"abstract":"This paper investigates the influence of indenter shape and rock texture on the crack growth properties and rock hardness. For this purpose, two different rock specimens such as basalt and marble with various textures were prepared and tested by Vickers indentation hardness device with rhombus indenter shape under two different temperatures of 35°C and 100°C. Concurrent with the experimental test, Vickers indentation simulations have been done on three calibrated rock models with nine different indenter shapes. The tensile strength of marble was 8 MPa, while basalt had a tensile strength of 10.8 MPa. Regarding compressive strength, marble exhibited 71 MPa, whereas basalt had a compressive strength of 153 MPa. Marble and basalt had elastic moduli of 45 and 95 GPa, respectively. The physical loading was applied vertically at a rate of 0.004 mm/min. The rock's texture and temperature significantly impact Vickers indentation hardness. Penetration by the Vickers indenter creates an elastic‐plastic stress field, leading to radial cracks due to exceeding the critical stress level. Ring cracks are caused by bulging of the test material and nested cracks directly under the indentation due to high shear and bending stresses in the region. The indentation shape affected the extent of the damage zone below it, leading to increased crack growth with smaller indentation diameters. The radial fracture number increased when the cross‐section changed from circular to quadrilateral shape. The crack growth and damage zone area increased with higher rock temperature due to increased rock brittleness.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142589024","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}
Xiangyun Zhou, Qingkai Zhang, De'an Sun, You Gao, Minjie Wen, Yunzhi Tan
{"title":"Analytical Model for Heat Transfer Around Energy Piles in Layered Soil With Interfacial Thermal Resistance by Integral Transform Method","authors":"Xiangyun Zhou, Qingkai Zhang, De'an Sun, You Gao, Minjie Wen, Yunzhi Tan","doi":"10.1002/nag.3883","DOIUrl":"https://doi.org/10.1002/nag.3883","url":null,"abstract":"Energy piles are commonly deployed in vertically layered geological conditions due to the geological structure and pile foundation backfill. The imperfect contact between adjacent soil layers results in resistance to heat transfer at the interface, known as the interfacial thermal resistance effect. In this paper, the energy pile was simplified as a finite‐length solid cylindrical heat source, and an analytical model was established for layered heat transfer of energy piles considering the interfacial thermal resistance effect. The Laplace‐domain solutions to the temperatures in the layered ground were derived by using the finite Hankel and Laplace transforms. The Crump method was subsequently employed to numerically invert Laplace‐domain solutions to the time‐domain solutions. The proposed model was validated by comparing with an analytical solution of a homogeneous model and COMSOL numerical solution. These solutions were used to analyze the temperature response around energy piles considering interfacial thermal resistance. Finally, a parametric study was performed to explore the effects of interfacial thermal resistance and other thermal properties of the soil layer on the layered heat transfer of energy piles.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142589025","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}
Jie Jiang, Yonglin Ai, Lijun Chen, Wencheng Chai, Mingxi Chen, Xiaoduo Ou
{"title":"Lateral Response Analysis of a Large‐Diameter Pile Under Combined Horizontal Dynamic and Axial Static Loads in Nonhomogeneous Soil","authors":"Jie Jiang, Yonglin Ai, Lijun Chen, Wencheng Chai, Mingxi Chen, Xiaoduo Ou","doi":"10.1002/nag.3884","DOIUrl":"https://doi.org/10.1002/nag.3884","url":null,"abstract":"The large diameter piles are widely used in structures such as offshore wind turbines due to their superior lateral load‐bearing capacity. To explore the lateral response of a large‐diameter pile under combined horizontal dynamic and axial static loads in nonhomogeneous soil, a simplified analytical model of the pile–soil interaction is developed. This model represents the pile as a Timoshenko beam resting on the Pasternak foundation, incorporating the double‐shear effect by considering both pile and soil shear. The governing matrix equations for the pile elements are derived from the principle of virtual work. Further, the pile's lateral deformations and internal forces are obtained using the modified finite beam element method (FBEM) and then validated through existing analytical solutions. Finally, the contribution of various properties of pile, soil, and applied load to the pile's lateral vibration response are performed. It is found that both pile and soil shear effects significantly impact the lateral dynamic response of a large‐diameter pile. Additionally, in nonhomogeneous soil, decreasing surface soil strength and dimensionless frequency lead to increased lateral displacements and bending moments of the pile, which are significantly affected by the <jats:italic>P</jats:italic>‐Δ effect under increasing axial load.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142580360","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}
Zhongxian Liu, Shuo Zhu, Alexander H. D. Cheng, Zhenen Huang
{"title":"Nonsingular Indirect Boundary Element Method and Multistrategy Particle Swarm Optimization Algorithm Applied to 3D Subsurface Cavity Inversion","authors":"Zhongxian Liu, Shuo Zhu, Alexander H. D. Cheng, Zhenen Huang","doi":"10.1002/nag.3879","DOIUrl":"https://doi.org/10.1002/nag.3879","url":null,"abstract":"The inversion of subsurface geological structures is a crucial approach for gaining insights into the internal composition of the earth. In this paper, we propose a novel inversion method combining the nonsingular indirect boundary element method (IBEM) with the multistrategy particle swarm optimization (MSPSO) algorithm, tailored for accurately inverting 3D subsurface cavities. Leveraging the semi‐analytical nature of IBEM offers advantages such as dimensionality reduction, automatic fulfillment of radiation conditions at infinity, and high computational accuracy. Furthermore, to augment global optimization and local search capabilities, an MSPSO algorithm is introduced. Employing multiple optimization strategies enhances particle diversity, accelerates algorithm convergence, and mitigates the risk of local optima. Through the consideration of subsurface cavities with varying parameters, this method quickly identifies the approximate location of the cavity within a wide search range. The final results demonstrate that the proposed method can simultaneously and accurately invert the 3D spatial position, size, and orientation of the cavity.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-11-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142589030","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}
Hamza Mhamdi Alaoui, Chaofa Zhao, Wenbo Niu, Pierre‐Yves Hicher
{"title":"Modeling the Swelling Behavior of Clayey Geomaterials Across Scales: Advances and Challenges","authors":"Hamza Mhamdi Alaoui, Chaofa Zhao, Wenbo Niu, Pierre‐Yves Hicher","doi":"10.1002/nag.3881","DOIUrl":"https://doi.org/10.1002/nag.3881","url":null,"abstract":"Most soils and rocks contain varying fractions of clay minerals within their solid matrix. These geomaterials can exhibit a significant swelling potential toward chemo‐thermo‐hydromechanical loadings. Several multiscale modeling techniques have been developed to ascertain their swelling behavior across various scales, with molecular dynamics (MD), micromechanics‐based approaches, and double‐porosity models being the most common. MD simulation is a computational technique that applies Newton's second law of motion to depict the movement of particles within a granular system. Micromechanics‐based approaches upscale the poro‐elasticity law from the clay layer level to the sample scale through homogenization. Dual‐porosity models are generally based on elasto‐plasticity, incorporating different hydro‐mechanical laws at two distinct scales. These models have been extensively used, particularly for clayey soils and bentonites, though their application to clayey rocks has not been reported in the literature. Although their significant contribution to the understanding of clay swelling behavior, these techniques have been insufficiently reviewed, compared, and discussed mutually in the literature. This paper aims to provide a cross‐look on these multiscale approaches by presenting the theoretical background of existing formulations, highlighting breakthrough results, discussing major differences and current challenges, and proposing future perspectives.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142555783","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 Brittle Failure in Rock Slopes Using Semi‐Lagrangian Nonlocal General Particle Dynamics","authors":"Peng Yin, Xiao‐Ping Zhou, Jinhu Pan","doi":"10.1002/nag.3882","DOIUrl":"https://doi.org/10.1002/nag.3882","url":null,"abstract":"The nonlocal general particle dynamics (NGPD) has been successfully developed to model crack propagation and large deformation problems. In this paper, the semi‐Lagrangian nonlocal general particle dynamics (SL‐NGPD) is proposed to solve brittle failure in rock slopes. In SL‐NGPD, the interaction between particles due to deformation is calculated in the initial configuration, while the friction contact interaction from discontinuities is calculated in the current configuration. The Van der Waals force model is utilized for friction contact. The bond‐level energy‐based failure criterion is developed to predict tensile/compressive‐shear mix‐mode cracks. The artificial viscosity and damage correction are used to enhance the numerical stability and accuracy when modeling brittle failure. The SL‐NGPD paradigm is numerically implemented through adaptive dynamic relaxation and predictor–corrector schemes for stable numerical solutions. The stability and accuracy of SL‐NGPD are verified by simulating compression tests. Thereafter, the crack coalescence patterns of double‐flaw specimens are investigated to understand the triggering failure mechanism of jointed rock slopes. Finally, the progressive failure process of the rock slope with step‐path joints is simulated to demonstrate its validity and robustness in modeling brittle failure in rockslides. The numerical results illustrate that the proposed SL‐NGPD is promising and performant for analyzing brittle failure problems in geotechnical engineering.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-10-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142555782","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}
Raúl Castro, Diego Oyarzo, René Gómez, Kimie Suzuki, Miguel Cifuentes
{"title":"Coupling Geomechanical and Gravity Flow Models to Obtain More Representative Flow Simulations and Air‐Gap Risk Identification in Caving Mining","authors":"Raúl Castro, Diego Oyarzo, René Gómez, Kimie Suzuki, Miguel Cifuentes","doi":"10.1002/nag.3880","DOIUrl":"https://doi.org/10.1002/nag.3880","url":null,"abstract":"The extraction and propagation of caving are complex phenomena involving the breaking of the rock mass, the formation of a column of broken material, and the extraction from the column base. Geomechanical modeling in cave mining commonly uses approaches to model the rock mass as a continuous material, while discontinuous modeling is frequently used for the column of broken material. However, it remains complex to include all mechanisms in a single model. Therefore, to achieve a better representation of ore breakage and extraction in caving mining, this work couples FLAC<jats:sup>3D</jats:sup>, a continuous finite volume tool, with FlowSim, a discrete tool based on cellular automata, to determine the air gap volume. The methodology first defines the height of caving propagation and the cave back with a tool that models solid rock mass in a continuous manner, which are used to constrain the cellular automata tool that simulates the flow of broken material. The results show that unidirectional FLAC<jats:sup>3D</jats:sup>‐FlowSim coupling reproduces the generation of cave backs and air gaps in the propagation of caving, rendering the methodology valuable for preliminary estimation of air volumes over fragmented material and the generation of supportive data to control the caving process.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":null,"pages":null},"PeriodicalIF":4.0,"publicationDate":"2024-10-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142490788","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}