Lilin Wang , Shaoyang Wang , Lizhong Wang , Yi Hong
{"title":"A dynamic p-y + M-θ model for monopile in soft clay considering failure mechanism under combined actions of wind and earthquake","authors":"Lilin Wang , Shaoyang Wang , Lizhong Wang , Yi Hong","doi":"10.1016/j.compgeo.2025.107074","DOIUrl":"10.1016/j.compgeo.2025.107074","url":null,"abstract":"<div><div>The increasing turbine sizes have necessitated monopile in soft clay to have larger diameter and rigidity, from early design of flexible piles to recent semi-rigid piles, with a future outlook on rigid piles. Existing failure mechanism-based soil-pile interaction model, i.e. <em>p-y + M-θ</em> model, is specifically developed for monopiles under lateral wind loading in non-seismic areas. To date, there is still a lack of failure mechanism-based <em>p-y + M-θ</em> model considering the combined actions of wind loading and seismic loading that is transmitted upward from the pile toe. This study aims to (a) reveal the failure mechanisms of monopile with varying rigidity under combined wind and seismic loading, and (b) to develop a dynamic <em>p-y + M-θ</em> model in accordance with these mechanisms. The first objective is achieved through a series of 3D finite element analyses well-calibrated by centrifuge model tests, which reveal a new mechanism (i.e., translation-shear failure) introduced by seismic loading, as an addition to the three-zone failure mechanism typically observed for a pile solely under wind loading. A dynamic <em>p-y + M-θ</em> model is then developed in light of these failure mechanisms associated with both wind and seismic loadings, with hysteretic damping and frequency-dependent radiation damping specifically introduced to enable dynamic analyses. The new model is validated against numerical analyses on piles subjected to seismic and wind loadings. Compared to the authors’ original <em>p-y + M-θ</em> model, the newly proposed model can better describe dynamic soil-pile interaction in seismically active areas, as it poses two additional simulation capabilities: (a) amplified lateral pile displacement due to the translation-shear failure caused by the seismic movement of whole pile embedment; (b) suppressed structural response due to the radiation damping aroused from the high-frequency seismic movement.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107074"},"PeriodicalIF":5.3,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143171742","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}
Zhenkun Xie , Shili Qiu , Shaojun Li , Quan Jiang , Dingping Xu , Minzong Zheng
{"title":"Excavation damage mechanism of deep buried layered fractured rock mass based on three-dimensional bonded block damage model","authors":"Zhenkun Xie , Shili Qiu , Shaojun Li , Quan Jiang , Dingping Xu , Minzong Zheng","doi":"10.1016/j.compgeo.2025.107101","DOIUrl":"10.1016/j.compgeo.2025.107101","url":null,"abstract":"<div><div>Under complex mineralization, the geological environment of deep mining projects is often accompanied by fractured rock masses. The existence of structural planes and cracks control the mechanical behavior of fractured rock masses. To describe the mechanical response mechanism of deep buried fractured rock mass, a three-dimensional bonded block damage constitutive model (BBDM) is proposed in this paper. Based on the damage characteristics of rock mass, the model will degrade the tensile strength, cohesion, dilation angle, normal and shear stiffness parameters of the joint based on the fracture energy value when the joint is in tension and shear yield state, and make the model eventually degenerate into a pure friction Mohr-Coulomb model under zero cohesion. Meanwhile, taking a deep buried roadway excavation project as the research background, the 610 m main slope excavation process is simulated by using the BBDM. Combined with the field test results, the stress, displacement and joint damage law of the surrounding rock excavation process are analyzed. The results show that in the closer position to the side wall, the potential interlayer fracture damage is larger, and the damage mechanism is mainly tensile damage. With the increase of the distance from the side wall, the damage degree gradually decreases, and the damage mechanism becomes mainly compressive shear damage, and eventually transitions to the state of no damage to the cracks. The research results reveal the damage process and failure mechanism of interlayer fracture in fractured rock bodies, which deepens the understanding of the mechanical response of deeply buried fractured rock masses and is significant for ensuring the stability of surrounding rocks and the safe and efficient production of the mining area.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107101"},"PeriodicalIF":5.3,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172652","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}
Zhaowei Ding , Wang Wu , Chunyu Song , Lingsen Zhao , Shengli Chen
{"title":"Novel semi-analytical model for the transient response of laterally loaded pile considering geometric nonlinear behaviors","authors":"Zhaowei Ding , Wang Wu , Chunyu Song , Lingsen Zhao , Shengli Chen","doi":"10.1016/j.compgeo.2025.107112","DOIUrl":"10.1016/j.compgeo.2025.107112","url":null,"abstract":"<div><div>This paper presents a semi-analytical model based on the radiant stress theory for analyzing the transient response of pile foundations under impulse loading. Geometric nonlinear behaviors at the pile-soil interface, including the sliding and debonding, are properly considered through the introduction and implementation of the mixed boundary conditions. Laplace transform and Durbin inversion algorithm are employed to calculate the transient response of the laterally loaded pile in the time domain. An iterative strategy is proposed to determine the depth range of geometric nonlinearity. Comparison with the results from finite element method confirms the reliability of the semi-analytical model and demonstrates the significance of incorporating the geometric nonlinearity. Neglecting such nonlinear behaviors can lead to an underestimated pile displacement amplitude and a significantly overrated radiation damping, thereby eliminating the rebound phase and potentially resulting in an overdamped response. Furthermore, extensive parametric analyses are conducted to investigate the influences of modulus ratio, impulse duration, and pile slenderness ratio on the transient response of pile. The numerical results show that neglecting geometric nonlinearity tends to diminish the influences of modulus ratio, while high-frequency impulse loading leads to amplified rebound and re-impact phenomena. The semi-analytical model may serve as an efficient and accurate tool for analyzing and optimizing pile foundation design, offering a practical alternative to computationally more intensive numerical methods.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107112"},"PeriodicalIF":5.3,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172645","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":"Numerical investigation on the landslide dam formation in landslide-river interaction","authors":"H.Y. Luo , P. Shen , L.M. Zhang , J. He","doi":"10.1016/j.compgeo.2025.107118","DOIUrl":"10.1016/j.compgeo.2025.107118","url":null,"abstract":"<div><div>Landslides occur in valleys always interact with the rivers. The sliding materials are likely to deposit in the river and form a landslide dam that impedes the river water flow, threatening human lives and properties both upstream and downstream. It is thus essential to investigate the formation mechanisms of landslide dams during landslide-river interaction. In this study, a multi-phase depth-averaged model is adopted to systematically explore the key factors that influence landslide dam formation, including the landslide volume, river discharge, landslide discharge, initial solid concentration and internal friction angle of landslide material. It is found that the landslide volume is always the key on landslide dam formation. The river water plays a critical role on retarding the inertia dynamic of landslide and accelerating the deposition process. One interesting finding is that the river blockage condition is insensitive to both the river and landslide discharge rates while the river flow depth is a more direct factor that controls river blockage. The two material properties on landslide dam formation (i.e., the initial solid concentration and internal friction angle) are represented by the enhanced flow mobility. After a landslide dam forms, the lasting impact of river water leads to the evolution of landslide deposit in the river. The kinetic energy ratio of landslide deposit to river water is linearly correlated with the solid concentration of deposit. The in-depth study on landslide formation mechanisms provides a solid basis for the evaluation of landslide hazard chain and risk mitigation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107118"},"PeriodicalIF":5.3,"publicationDate":"2025-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104972","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":"Undrained cyclic and post-cyclic shear behaviour of sand with varying liquefaction degrees: insights from DEM","authors":"Zhiyong Liu , Xinran Chen , Bo Liu , Jianfeng Xue","doi":"10.1016/j.compgeo.2025.107116","DOIUrl":"10.1016/j.compgeo.2025.107116","url":null,"abstract":"<div><div>The stability of geotechnical structures after an earthquake is primarily determined by the residual strength of surrounding soils that have not fully liquefied. This research employs the discrete element method (DEM) to study the undrained post-cyclic shear behaviour of sand under triaxial conditions, focusing on the effect of varying degrees of liquefaction (LD) simulated by subjecting the samples to different lengths of cyclic loading. Different types of cyclic loading, i.e. symmetric (fully reversal), partially reversal, and non-reversal ones, as well as the effect of sample density, have been considered. The results indicate that the samples under fully or partially reversal cyclic loading eventually liquefied, displaying a cyclic mobility failure mode. In contrast, samples under non-reversal cyclic loading develop plastic strain accumulation (PSA) failure without liquefaction. The post-cyclic shear stiffness of the samples is affected by both LD and the type of cyclic loading. For samples under reversal cyclic loading, the post-cyclic shear stiffness decreases as LD increases. Notably, the liquefied samples (LD = 1) initially exhibit near-zero stiffness during post-liquefaction shear until highly anisotropic force chains are formed along the loading direction, with their buckling leading to stiffness recovery. The length of the low-stiffness stage is influenced by the static shear stress and the relative density of the sample, which determines the rate of anisotropy accumulation during cyclic loading. The onset and completion of stiffness recovery are marked by a peak in anisotropy and an abrupt increase in effective anisotropy, respectively. For samples under non-reversal cyclic loading, the post-cyclic shear stiffness initially decreases with the increase in LD but increases at higher LDs due to the significant anisotropy developed during the cyclic loading stage.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107116"},"PeriodicalIF":5.3,"publicationDate":"2025-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172644","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":"A novel depth-averaged model of landslide over erodible bed using (b, s) coordinates","authors":"Van Khoi Pham , Changhoon Lee , Van Nghi Vu","doi":"10.1016/j.compgeo.2025.107105","DOIUrl":"10.1016/j.compgeo.2025.107105","url":null,"abstract":"<div><div>A novel model of landslides over erodible beds is developed using the nonlinear shallow water equations (NSWE). In the model, the <span><math><mrow><mo>(</mo><mi>b</mi><mo>,</mo><mspace></mspace><mspace></mspace><mi>s</mi><mo>)</mo></mrow></math></span> coordinates, which follow the global coordinate with the <em>z</em>-axis vertically, are used and the debris surface from the reference line (<em>s</em>) is used as a variable rather than the conventional debris depth in terms of conveniences. The erodible beds are employed in the model by including the erosion rate in the bottom boundary condition. A hybrid finite volume-finite difference scheme is applied to discretize the set of governing equations. The present simulations are compared with experimental data and other numerical simulations to figure out the advantage of this model. Finally, the two case studies of Sindonga and Raemian landslides at Umyeon mountain in the year 2011 are simulated for real applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107105"},"PeriodicalIF":5.3,"publicationDate":"2025-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172651","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}
Lei Xing , Wenping Gong , Jinsong Huang , Hongbo Zhang , Baoyin Xing , Lei Wang
{"title":"An improved CFD-DEM coupling method for simulating the steady seepage-induced behaviors of soil-rock mixture slopes","authors":"Lei Xing , Wenping Gong , Jinsong Huang , Hongbo Zhang , Baoyin Xing , Lei Wang","doi":"10.1016/j.compgeo.2025.107069","DOIUrl":"10.1016/j.compgeo.2025.107069","url":null,"abstract":"<div><div>Modeling the seepage-induced progressive failure of soil-rock mixture (SRM) slopes is challenging because of the large deformation and the complexity of stress-seepage coupling. To address these challenges, this study presents a numerical approach for modeling the seepage-induced progressive failure of SRM slopes, where rock blocks within slopes and the effect of cracks on seepage behavior during slope failure are explicitly considered and modeled. Within the proposed method, the traditional unresolved computational fluid dynamics-discrete element method (CFD-DEM) coupling method is first improved by introducing the unstructured mesh and hydraulic boundary condition to model the seepage-induced progressive failure of SRM slopes. The seepage behavior within SRM slopes is modeled with a CFD solver based on finite volume methods, while the particle motion under these interaction forces is simulated with a DEM solver. Note that fluid-particle interaction forces are calculated with empirical equations in the improved CFD-DEM coupling method. CFD-DEM coupling is achieved by exchanging data between the two solvers at each timestep. The proposed numerical method’s effectiveness is illustrated through two seepage problems (within an SRM sample and an SRM slope) and three model tests with different rock contents (in terms of 10%, 20%, and 30%).</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107069"},"PeriodicalIF":5.3,"publicationDate":"2025-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104973","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":"A novel surrogate model for hydro-mechanical coupling in unsaturated soil with incomplete physical constraints","authors":"Charles W.W. Ng, Qianyu Zhou, Qi Zhang","doi":"10.1016/j.compgeo.2025.107091","DOIUrl":"10.1016/j.compgeo.2025.107091","url":null,"abstract":"<div><div>Physics-informed neural networks (PINNs) are increasingly employed for surrogate modelling of soil behaviour. Existing surrogate models for unsaturated soil only account for seepage in rigid soil, neglecting the complex coupling between deformation and seepage in unsaturated soil. This study develops a new surrogate model for hydro-mechanical coupling in unsaturated soil using the PINN approach. Dimensionless governing equations, including mass balance and force balance equations, are derived and adopted for physical constraints. With absence of explicit constitutive relations, this new surrogate model utilises sparse measured data to identify pore water pressure, effective stress and deformation in unsaturated soil. Separate neural networks are employed to facilitate efficient back-propagation for coupled problem involving multiple outputs. The newly developed model is then applied to simulate two cases with sparse measurements in unsaturated soil. The results illustrate that the newly developed surrogate model successfully learns the elasto-plastic constitutive relation of suction-induced volume change from experimental data. Meanwhile, model predictions regarding both water flow and stress distribution align within the 95 % confidence interval of theoretical values, demonstrating interpretability of PINN model. Furthermore, by adhering to physical constraints, the relative error in predicting soil deformation from neural networks significantly reduces from 49 % to less than 10 %. These findings suggest PINN model with separate networks is capable to simulate unsaturated soil considering both deformation and seepage, even with sparse measured data and incomplete physical constraints.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107091"},"PeriodicalIF":5.3,"publicationDate":"2025-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172647","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}
Jianxiong Li , Zhanyuan Zhu , Wen Hua , Yi Yang , Shiming Dong , Tiankui Guo
{"title":"Numerical investigation of fracture competitive propagation mechanisms for temporary plugging staged fracturing (TPSF) in shale gas reservoirs","authors":"Jianxiong Li , Zhanyuan Zhu , Wen Hua , Yi Yang , Shiming Dong , Tiankui Guo","doi":"10.1016/j.compgeo.2025.107109","DOIUrl":"10.1016/j.compgeo.2025.107109","url":null,"abstract":"<div><div>Temporary plugging staged fracturing (TPSF) is a technique rapidly growing for enhancing the uniform growth of multiple fractures in shale gas reservoirs. However, the mechanisms of fracture competition during TPSF remain incompletely understood, primarily due to pre-existing natural fractures and stress interference. This study presents a fully coupled Finite Element Method (FEM) model incorporating the pore pressure cohesive zone method to investigate how fractures compete for growth in TPSF. Moreover, a user-defined perforation element connected to a discrete fracture network model is employed to analyze how fluid flow is partitioned. The results indicate that denser natural fractures (NFs) create more connections for hydraulic fractures (HFs), thereby increasing the complexity of the fracture network. Reduction in cluster spacing can result in distorted fracture morphologies and altered propagation pathways, potentially halting fracture growth. Higher injection rates elevate injection pressure, which may inhibit fracture growth but could aid in penetrating stress interference zones to generate asymmetrical fractures. Fractures tend to propagate towards larger inclination angles, with an inclination angle close to 45° being optimal for enhancing fracture growth and reservoir stimulation. This study offers valuable insights into the competitive propagation mechanism of fractures and provides guidance on optimal field design for TPSF.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107109"},"PeriodicalIF":5.3,"publicationDate":"2025-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172650","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}
Yaoying Liang , Ming Peng , Liu Liu , Siang Huat Goh , Dengyi Wang , Jian Shen
{"title":"Stability analysis of rock slopes considering strata uncertainty using dual-source surface wave inversion with borehole constraints","authors":"Yaoying Liang , Ming Peng , Liu Liu , Siang Huat Goh , Dengyi Wang , Jian Shen","doi":"10.1016/j.compgeo.2025.107068","DOIUrl":"10.1016/j.compgeo.2025.107068","url":null,"abstract":"<div><div>The stability of rock slopes is mainly controlled by subsurface characteristics, such as soft layers and broken zones. However, determining their locations and associated uncertainties is challenging due to their concealment and complexity. This study proposes a framework for stability analysis of rock slopes, which accounts for strata uncertainty quantified through dual-source surface wave inversion with borehole constraints. The shear wave velocity profiles of the slope were obtained by combining active surface wave and passive microtremor using multichannel analysis of surface waves (MASW) method and microtremor survey method (MSM), respectively. The prior constraints for the inversion were obtained from borehole television and ultrasonic testing. Subsequently, the inverted depths of the interfaces were represented as continuous distributions and fitted to various models. The stability of the slope was then analyzed by incorporating strata uncertainty. The proposed framework was applied to a steep road-cut rock slope located in Shandong Province, China. The results show that the proposed framework effectively detects soft layers and quantifies uncertainties in interface depths for probabilistic stability analysis. The dual-source surface wave inversion provides both sufficient detection depth and high resolution of shallow strata. Incorporating the prior borehole-derived constraints, the accuracy and speed of the dual-source surface wave inversion have been greatly enhanced. The interface depths can be quantitatively described by normal distributions, which have been selected based on the lowest Akaike Information Criterion (AIC) value. Three soft layers were identified in the presented case, and variations in the thickness of the soft layers influenced the safety factor.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107068"},"PeriodicalIF":5.3,"publicationDate":"2025-01-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172649","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}