Linfeng Zhang , Guorong Wang , Zhiyuan Li , Lin Zhong , Qiang Fu , Jiwei Wu , Mengdi Wu , Dongmei Liang , Yuhang Zheng
{"title":"A fully Coupled thermal-hydrodynamic–chemical numerical model for simulating gas hydrate-bearing sediments dissociation based on peridynamic differential operator method","authors":"Linfeng Zhang , Guorong Wang , Zhiyuan Li , Lin Zhong , Qiang Fu , Jiwei Wu , Mengdi Wu , Dongmei Liang , Yuhang Zheng","doi":"10.1016/j.compgeo.2025.107099","DOIUrl":"10.1016/j.compgeo.2025.107099","url":null,"abstract":"<div><div>During methane hydrate-bearing sediments (MHBS) extraction multiphase flow with heat transfer transitioning MHBS phase followed by pore dynamic evolution. Which involves a complex thermo-hydrodynamic (THC) coupled process. Effective theoretical and fully integrated numerical models are essential for evaluating the economic feasibility of gas hydrate production and comprehending the fundamental physical concepts. This study presents a novel approach to constructing a fully integrated THC numerical model by employing a peridynamic differential operator (PDDO). The dimensionless governing equations of Darcy’s law and the heat source approach are transformed into a non-local integral form. The Euler forward difference method is used for time integration. The accuracy and reliability of the newly developed PDDO THC model were tested by comparing it with experimental examples from other popular simulators, which provides an alternative to explicit modeling of the phase change process involving multiphase flow with heat transfer and may find comprehensive and useful applications to a variety of industrial and geophysical processes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107099"},"PeriodicalIF":5.3,"publicationDate":"2025-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143171743","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 sparse-memory-encoding GPU-MPM framework for large-scale simulations of granular flows","authors":"Hao Chen, Shiwei Zhao, Jidong Zhao","doi":"10.1016/j.compgeo.2025.107113","DOIUrl":"10.1016/j.compgeo.2025.107113","url":null,"abstract":"<div><div>The Material Point Method (MPM) is increasingly recognized as an effective tool for simulating complex granular flows. While GPU computing has been widely used in MPM applications for large-scale problems, its heavy reliance on contiguous memory distribution can significantly hinder efficiency and limit simulation capabilities due to memory capacity constraints. This study presents a sparse-memory-encoding framework that incorporates advanced algorithms to address these limitations in large-scale simulations. We introduce a novel algorithm for atomic-free dual mapping between material points and nodes, in conjunction with warp-wise particle-to-grid mappings organized within a block-cell-material point hierarchy. Moreover, the framework features an efficient memory shift algorithm that optimizes memory usage for material properties. This optimization enables the seamless integration of commonly used material constitutive models, including elastic, elastoplastic, and hyper-plastic models, as well as various iteration schemes such as “update stress first”, “update stress last”, and “modified update stress last” within a cohesive framework. Furthermore, the framework accommodates incorporating diverse boundary conditions, such as Dirichlet, Neumann, and arbitrary-shaped rigid body contact, thus broadening its applicability to real-world engineering challenges, including landslides. The framework can effectively and efficiently handle large-scale, high-fidelity simulations of granular flows.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107113"},"PeriodicalIF":5.3,"publicationDate":"2025-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143171740","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":"Experimental and numerical investigations of undrained stress relaxation behavior considering freeze-thaw-drying-wetting effects for saturated clay","authors":"Weilie Zou , Qiuyang Pei , Xilin Xia , Zhong Han","doi":"10.1016/j.compgeo.2025.107114","DOIUrl":"10.1016/j.compgeo.2025.107114","url":null,"abstract":"<div><div>The effects of freeze–thaw-drying-wetting cycles on undrained stress relaxation behavior are studied from the perspectives of physical experiments and numerical simulations using the proposed multi-set material point method (MPM) algorithm. The effectiveness of the established MPM model is validated by comparing the simulation results with experimental data. To gain insights into the microscopic characteristics of undrained stress relaxation, the evaluation focuses on the variations in elastic Green-Lagrangian strain, water pressure, water particle velocity, and solid particle space. The deviator stress relaxation and water pressure evolution are correlated with the decrease in elastic strain and spatial variation in water velocity, respectively. The stress relaxation process breaks the interparticle bonds and results in structural weakening. Freeze-thaw-drying-wetting cycles reduce the relaxed deviator stress within different relaxation stages and promote variation in water pressure, instability in water velocity field, and adjustment of solid particles.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107114"},"PeriodicalIF":5.3,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143171741","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}
Penghao Zhang, Yunxuan Cui, Kurt Douglas, Chongmin Song, Adrian R. Russell
{"title":"Phase field fracture modeling of cohesive-frictional materials like concrete and rock using the scaled boundary finite element method","authors":"Penghao Zhang, Yunxuan Cui, Kurt Douglas, Chongmin Song, Adrian R. Russell","doi":"10.1016/j.compgeo.2025.107106","DOIUrl":"10.1016/j.compgeo.2025.107106","url":null,"abstract":"<div><div>Understanding and predicting the mechanisms and behaviors of crack propagations in engineering structures made of rock and concrete is important when producing reliable designs. This study proposes a phase field crack model suitable for cohesive-frictional materials such as rock and concrete. Novelty lies in introducing and coupling the multiaxial strength criterion for cohesive-frictional materials and the micro-damage evolution law within the fracture process zone within the basic phase field method. Another novel feature is the decomposition of volumetric and deviatoric components of the stiffness matrices in the precomputation phase of the scaled boundary finite element method. The combination of this numerical technique with the proposed constitutive model enables the effective and efficient simulation of compression-shear cracks. Novelty also lies in integrating the proposed model with scaling theory. This keeps the computational cost for large-scale problems within an acceptable range. The technique we adopt to achieve this advancement involves making the simulation results insensitive to the characteristic length. This improvement allows the characteristic length for large-scale problems to be uniquely determined and larger than the value used at the laboratory scale. This enables the use of coarser meshes, reducing the computational resources needed for simulating large-scale problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107106"},"PeriodicalIF":5.3,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104975","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}
Zhi Yong Ai, He Wei Kuang, Zi Kun Ye, Qing Song Lu
{"title":"Interaction between piles and layered fractional viscoelastic soils considering groundwater level","authors":"Zhi Yong Ai, He Wei Kuang, Zi Kun Ye, Qing Song Lu","doi":"10.1016/j.compgeo.2025.107119","DOIUrl":"10.1016/j.compgeo.2025.107119","url":null,"abstract":"<div><div>Based on the time-varying solution of layered fractional viscoelastic soils considering groundwater level, the interaction between soils and piles is studied by a finite element–boundary element coupling method. The single pile is regarded as a one-dimensional compressible bar to establish its stiffness matrix equation by the finite element method. Then, the fundamental solutions of layered fractional viscoelastic soils considering groundwater level are used as kernel functions in the boundary element method. According to the displacement continuity condition on the pile-soil interface, the interaction between single pile and soils is solved. Furthermore, by incorporating the pile-pile interaction as well as the balance and coordination conditions of the rigid cap, the solution of the single pile is extended to the pile group. Finally, the correctness of the presented method and program is validated, and several examples are designed to explore the effects of groundwater level, soil properties, pile parameters, and soil stratification on the pile-soil interaction. The analysis results show that, as the groundwater level declines, the displacement of the pile at the initial moment increases, and the reaction forces of center and side piles increase. Moreover, the increase of fractional order can accelerate the rate of pile displacement, reduce the axial force on the upper pile body, and decrease the reaction force at the top of the center and side piles.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107119"},"PeriodicalIF":5.3,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143171744","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}
Sen Zheng , Weihua Li , Yexin Wan , Zhe Yang , Sainan Zhu
{"title":"Wave propagation in an ocean site considering fractional viscoelastic constitution of porous seabed","authors":"Sen Zheng , Weihua Li , Yexin Wan , Zhe Yang , Sainan Zhu","doi":"10.1016/j.compgeo.2025.107098","DOIUrl":"10.1016/j.compgeo.2025.107098","url":null,"abstract":"<div><div>Wave propagation in an ocean site is an essential research topic in various scientific fields, such as offshore geotechnical engineering, ocean seismology, and underwater acoustics. Previous studies have considered the seabed soil as elastic or poroelastic, ignoring the viscoelastic characteristics of its solid skeleton. Based on the fractional-derivative viscoelastic theory and the modified Biot theory, considering the flow-independent viscosity related to solid skeleton, this paper proposes a generalized viscoelastic wave equation for a fluid-saturated porous medium. The equation has a flexible mathematical form to describe soil rheological properties more accurately through fractional order. On this basis, the total wave field equation of an ocean site, modeled as the fluid–poroviscoelastic–solid media, is established. Then an analytical solution for wave propagation in an ocean site subjected to obliquely incident P and SV waves is obtained, and its degeneration and extension are studied. The proposed method is comprehensively validated through experiment, analytical, and numerical methods. Finally, a parameter analysis is performed to investigate the effects of water depth, seabed properties (including viscoelastic parameters, fractional order and permeability), and incident angle on the seismic response of a poroviscoelastic seabed.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107098"},"PeriodicalIF":5.3,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172643","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}
Shao-Lin Ding , Jia-Jun Pan , Yanli Wang , Han Xu , Dian-Qing Li , Xin Liu
{"title":"Developing a digital twin for dam safety management","authors":"Shao-Lin Ding , Jia-Jun Pan , Yanli Wang , Han Xu , Dian-Qing Li , Xin Liu","doi":"10.1016/j.compgeo.2025.107120","DOIUrl":"10.1016/j.compgeo.2025.107120","url":null,"abstract":"<div><div>A dynamic, reliable, and even automatic evaluation of dam performance (e.g., deformation) is essential for safety management of earth and rockfill dams. The digital twin has emerged as a promising tool for smart dam safety management. Many countries have launched projects to develop digital twins of dams, which represent a living model that can continuously learn from monitoring data from the physical counterpart and produce a complete and timely description of the dam’s performance. A practical framework is proposed in this study to develop digital twin of an earth or rockfill dam for predicting its mechanical responses. The proposed framework utilizes both a physics-based dam model and monitoring data to enhance the model’s performance through Bayesian updating. It is illustrated by an operational digital twin project, i.e., the Danjiangkou Digital Twin Project for developing the digital twin of the core-wall rockfill dam at the right bank. The proposed method enables three novel features for dam safety management, namely, real-time simulation, future forecast, and scenario projection of dam performance. The results showed the updated dam model predicted the dam crest settlements accurately with an RMSE as small as 4.20 mm, verifying the effectiveness of the proposed framework.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107120"},"PeriodicalIF":5.3,"publicationDate":"2025-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104944","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}
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}