{"title":"Modelling of slurry infiltration in saturated porous media integrating particle dispersion and blockage","authors":"Jianxin Ning, Maosong Huang, Jian Yu","doi":"10.1016/j.compgeo.2024.106875","DOIUrl":"10.1016/j.compgeo.2024.106875","url":null,"abstract":"<div><div>Existing analytical solutions available for simulating slurry infiltration do not account for the effect of particle dispersion and blockage at the same time. In view of this, a mathematical model of slurry particle migration in a saturated porous media with the convection–dispersion-deposition effect is established, and a semi-analytical solution of the particle transport problem with time discretization is obtained using an integral transformation. The correctness and rationality of the method are verified by comparing the experimental and theoretical results of the one-dimensional particle transport problem in constant pressure injection and constant velocity injection modes. The results of the method are in good agreement with those obtained from commercial finite element analyses. The spatiotemporal distribution of particle concentration, slurry deposition, soil porosity, pore water pressure, flow velocity in soil column can be easily obtained by the method. The pore water pressure calculated using this method is compared with the measured result during slurry shield drilling stops, demonstrating the potential application of this method in slurry engineering. The parametric analysis indicated that increasing the excavation chamber pressure and slurry concentration can accelerate the pore water pressure dissipation and the mud cake formation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106875"},"PeriodicalIF":5.3,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142653000","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method","authors":"Fan Fei , Jinhyun Choo","doi":"10.1016/j.compgeo.2024.106890","DOIUrl":"10.1016/j.compgeo.2024.106890","url":null,"abstract":"<div><div>The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. We transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106890"},"PeriodicalIF":5.3,"publicationDate":"2024-11-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142652999","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"DEM simulation of cone penetration tests in sand in a virtual calibration chamber","authors":"Yue Song , Xiaoqiang Gu , Jing Hu , Xing Zheng","doi":"10.1016/j.compgeo.2024.106900","DOIUrl":"10.1016/j.compgeo.2024.106900","url":null,"abstract":"<div><div>In this paper, cone penetration tests (CPT) in a virtual calibration chamber (VCC) were simulated using the discrete element method (DEM) to investigate the factors influencing the cone penetration tip resistance (<em>q<sub>c</sub></em>) in sand. The microscopic parameters of the Fontainebleau sand were calibrated by comparing the simulated macroscopic mechanical behavior to that from experimental drained triaxial tests. The influences of VCC dimensions and prescribed boundary conditions, penetration rate, and confining stress on the <em>q<sub>c</sub></em> were investigated. The results show that the <em>q<sub>c</sub></em> under different boundary conditions converges as the diameter of the VCC increases. It is recommended that a ratio of chamber diameter to cone diameter larger than 20 should be used to eliminate the potential boundary effects for CPT in VCC for samples of different relative density (<em>D<sub>r</sub></em>). The comparison between the numerical simulations and cavity expansion theory predictions presents significant relationships among critical state parameters (<em>Ψ</em><sub>cri</sub>), peak internal friction angle (<em>φ</em><sub>peak</sub>) and normalized cone tip resistance (<em>Q</em>).</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106900"},"PeriodicalIF":5.3,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660211","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yuyang You , Junsheng Yang , Xiangcou Zheng , Yipeng Xie , Jingkang Lyu , Ashraf S. Osman
{"title":"A 3D material point discretization approach for complex terrain and geological body: Numerical implementation and application","authors":"Yuyang You , Junsheng Yang , Xiangcou Zheng , Yipeng Xie , Jingkang Lyu , Ashraf S. Osman","doi":"10.1016/j.compgeo.2024.106884","DOIUrl":"10.1016/j.compgeo.2024.106884","url":null,"abstract":"<div><div>The construction of three-dimensional numerical models for complex terrains and geological bodies has always posed a challenge in geoengineering. This study introduces an efficient discretization approach that facilitates the establishment of material point models for three-dimensional complex terrain and geological conditions. The proposed approach incorporates three structural components and employs an efficient algorithm for data processing to characterize terrain and geological bodies within the model. This enables the construction of complex, heterogeneous discrete models using elevation and material information as input data. A detailed procedure for implementing the proposed three-dimensional material point discretization approach is outlined. Additionally, robust metrics for discretization deviation and efficiency evaluation are introduced to assess the accuracy and efficiency of the constructed discretized models. The efficacy of the proposed approach is verified and evaluated through both two-dimensional and three-dimensional examples, demonstrating its accuracy and effectiveness.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106884"},"PeriodicalIF":5.3,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660212","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhi Ye , Yongchi Li , Lihua Li , Wenjin Xu , Hong Wu , Yu Miao , Huabei Liu
{"title":"Seismic response of a long shield tunnel crossing through multiple soil deposits","authors":"Zhi Ye , Yongchi Li , Lihua Li , Wenjin Xu , Hong Wu , Yu Miao , Huabei Liu","doi":"10.1016/j.compgeo.2024.106892","DOIUrl":"10.1016/j.compgeo.2024.106892","url":null,"abstract":"<div><div>This study investigated the seismic response of a long shield tunnel crossing multiple soil deposits. The P2PSand model and practical-oriented hysteretic model were used to simulate the cyclic behaviour of liquefiable sand and soft clay, respectively. A refined tunnel model was adopted to simulate the mechanical behaviour of segmental linings. Dynamic analysis was then performed to explore the seismic response of shield tunnels passing through liquefiable sand and soft clay under transverse or longitudinal asynchronous ground motions. The joint displacement and lining force were found to be prominent near the interface between the liquefiable sand and soft clay. The counterintuitive seismic deformation mechanism of tunnel structure was then clarified. Parametric studies were carried out to elucidate the tunnel deformation mechanism. The seismic responses of the tunnel structure under longitudinal ground motions were more prominent than those under transverse ground motions. The earthquake intensities, relative densities and liquefaction potentials of sand, and clay stiffness exerted major influences on the seismic response of tunnel structure, while the interface inclinations and earthquake wave velocities had a slight impact on the tunnel seismic response.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106892"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660209","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Hyperspheric Integral Reliability Method for efficient reliability analysis of geotechnical ultimate limit states","authors":"Ivan Depina","doi":"10.1016/j.compgeo.2024.106861","DOIUrl":"10.1016/j.compgeo.2024.106861","url":null,"abstract":"<div><div>This paper introduces the Hyperspheric Integral Reliability Method (HINT) for efficient reliability analysis of geotechnical ultimate limit states. The method is motivated by the mechanism of the Shear Strength Reduction Method (SSRM), which is often employed to calculate the factor of safety for geotechnical ultimate limit states. HINT exploits the observation that a factor of safety, computed by the SSRM, is not only a pointwise estimate of safety, but also a measure of distance to the failure limit along a radial direction in the space of the shear strength parameters. This observation is utilized to transform the reliability integral in the hyperspheric coordinate system and develop an efficient estimator of failure probability for geotechnical ultimate limit states. HINT was examined on several benchmarking problems, demonstrating stable and highly efficient performance on low to medium dimensional reliability problems. HINT can be also applied to general reliability problems if a radial search mechanism is implemented as in the SSRM. Given that the SSRM is already available in most commercial geotechnical software packages, HINT is perfectly suited to advance geotechnical reliability analyses and make them accessible to a wider set of use-cases.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106861"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660210","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhengzhen Wang , Zhentao Zhang , Guoliang Dai , Yong Zhou , Chengming Cao , Tiantao Su
{"title":"Reliability analysis of slopes reinforced by frame-anchor structures under earthquake conditions","authors":"Zhengzhen Wang , Zhentao Zhang , Guoliang Dai , Yong Zhou , Chengming Cao , Tiantao Su","doi":"10.1016/j.compgeo.2024.106885","DOIUrl":"10.1016/j.compgeo.2024.106885","url":null,"abstract":"<div><div>Earthquakes often cause slope collapse, which poses a serious threat to people’s lives and property. Therefore, it is important to study slope stability under earthquake conditions. The limit element upper bound method was improved and applied to the stability analysis of the slope reinforced by frame-anchor structures under earthquake action. The influences of anchor prestress and seismic force on the stability of slope were analysed comprehensively, and the soil was regarded as a heterogeneous material in its natural state under the action of weathering, deposition and historical stress. The rationality of the slope’s safety factor and failure mode obtained by this method was verified by an example. Based on the Monte Carlo method, considering the random distribution of soil parameters in space and the autocorrelation between any two points, the midpoint method based on Cholesky decomposition was used to simulate the relevant Gaussian random field. The solution program was compiled in MATLAB. The influence of several commonly used autocorrelation function types on the failure probability of slopes was analysed. The influence of the anchor rod prestress, anchor rod inclination angle and seismic force on the slope’s safety factor was discussed, and the failure mode of the slopes under a random field was studied. The results show that the velocity field boundary obtained by the reliability analysis method proposed in this paper is close to the sliding surface obtained by the traditional methods and the deviation of the slope’s safety factors between the proposed method and the traditional methods is small. The probability density functions of the slope’s safety factors under different autocorrelation functions are approximately normally distributed. Increasing the seismic force will increase the probability of slope failure, while the increase in anchor rod prestress and the decrease in anchor rod inclination angle have positive effects on the slope stability. The failure modes of slopes under different random distribution of soil parameters change significantly. The method proposed in this paper, which does not require to assume the failure model, can realistically and intuitively display the actual failure of slopes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106885"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660207","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Micro-fracture mechanism of microwave induced fracturing of basalt based on a novel Electromagnetic–Thermal–Mechanical coupling model","authors":"Jian Ma, Zheng-Wei Li, Wen-Feng Guo, Liang-Xiao Chen","doi":"10.1016/j.compgeo.2024.106874","DOIUrl":"10.1016/j.compgeo.2024.106874","url":null,"abstract":"<div><div>Microwave-assisted rock fracturing is recognized for its efficiency, energy savings, and environmental benefits. Investigating microscopic mechanisms of microwave-induced rock fracturing is essential for predicting the weakening effect on rock. A coupled Electromagnetic–Thermal–Mechanical model based on FEM–DEM was established to describe the response of rock under microwave irradiation. This model employs interpolation algorithms and mineral lattices randomization algorithms to establish a 2D cross-sectional representation of rock. A discrete element calculation method is proposed to synchronize computational time with the experimental time. The model can simulate the multi-physical field response of different rocks under various conditions, making it an effective tool for studying microwave-induced rock fracturing. The effectiveness of the numerical model was validated through open-end microwave-induced fracturing experiments on basalt. Additionally, the study elucidates the micro-fracture mechanism of basalt under microwave irradiation. The results indicate that the direction of crack propagation is influenced by microwave power and boundary effects. The patterns of fracture development between minerals are summarized as follows: Initial fractures primarily result from the rapid heating of microwave-absorbing minerals like enstatite,creating a significant temperature gradient. With increased heating time, heat transfers to highly expansive minerals such as olivine, causing fractures due to localized thermal expansion.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106874"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660206","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Fan Chen , Xiaohui Liu , Junfeng Sun , Hao Xiong , Zhen-Yu Yin , Xiangsheng Chen
{"title":"Micro mechanism and analytical model of stress distribution in loosened soil zone above active trapdoor","authors":"Fan Chen , Xiaohui Liu , Junfeng Sun , Hao Xiong , Zhen-Yu Yin , Xiangsheng Chen","doi":"10.1016/j.compgeo.2024.106841","DOIUrl":"10.1016/j.compgeo.2024.106841","url":null,"abstract":"<div><div>The formation of discontinuities within soil masses due to arching phenomena is a critical issue in geotechnical engineering, profoundly influencing the stability and integrity of underground structures. This study investigated the microscopic mechanisms and associated stress redistribution in the soil mass above an actively moving trapdoor. Utilizing the Discrete Element Method (DEM) model, the influence of varying initial soil packing densities on the induced arching phenomena is examined concerning the soil displacement mechanism and shearing configuration: transfer from inner triangular to vertical bands is observed in dense soil while only vertical shearing bands are observed in initially loose soil. To account for principal stress trajectories under differing initial density conditions, an analytical model was developed to quantify the variations in soil stress. The findings reveal a noteworthy impact of the initial soil void ratio <span><math><mi>e</mi></math></span>, on the pattern evolution and stress redistribution of soil arching, manifesting a passive stress limit coefficient, <span><math><msub><mrow><mi>K</mi></mrow><mrow><mi>p</mi></mrow></msub></math></span>, in denser soils along the central line above the trapdoor. The analytical model demonstrates robust agreement with both numerical and experimental data in soil loosening stress distribution and the load–displacement response of the trapdoor.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106841"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660213","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Finite-strain gradient-extended damage-plastic modeling of rock: Understanding catastrophe from material failure and structural instability","authors":"Xiaofeng Cheng , Xianhui Feng , Chun’an Tang","doi":"10.1016/j.compgeo.2024.106891","DOIUrl":"10.1016/j.compgeo.2024.106891","url":null,"abstract":"<div><div>Rock, a natural geological material, exhibits diverse collapse modes that make predicting catastrophic behavior challenging. Material failure and structural instability each provide independent mechanical explanations for rock catastrophes. However, isolated perspectives often obscure the distinctions and connections between these two critical mechanisms. Here, we propose a finite-strain gradient-extended damage-plastic scheme within a thermodynamically consistent framework to encapsulate the indispensable dual effects of failure and instability. The crack propagation and frictional dissipation of crack clusters in rock materials inspire coupled damage-plastic theory, whereas the large displacement, large rotation and large strain of rock structures motivate the application of finite strain theory. The proposed scheme incorporates nonlocal gradient-enhanced terms to mitigate mesh dependence and is immune to spurious energy dissipation under cyclic loading. Constitutive treatment at finite strain retains the easily achievable features of the small deformation case. The model is validated through laboratory- and engineering-scale simulations, offering insights into the mechanisms of rock catastrophes. Our findings highlight the dual role of material failure and structural instability as interconnected drivers of rock catastrophes, offering a more holistic understanding for effective prediction and mitigation strategies.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"177 ","pages":"Article 106891"},"PeriodicalIF":5.3,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660208","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}