Pub Date : 2025-02-07DOI: 10.1016/j.compgeo.2025.107128
Maozhu Peng , Zhen-Yu Yin , Zefeng Zhou
This paper presents one of the first systematic finite element investigations into the rate effect during plate anchor uplift in clay, with a focus on the influence of soil drainage. Some numerical difficulties in simulating partially/fully drained anchor uplift are first discussed, and a method to address these issues is proposed. The most striking benefit of the proposed method is that it allows the hydro-mechanical coupled anchor-soil interaction to be efficiently simulated, with improved numerical stability, by simply using a non-detachable anchor-soil interface. After validation, the proposed method is used to simulate anchor uplift at different embedment and uplift rates, spanning from drained to fully undrained conditions. The influence of uplift rate on the force–displacement relations and failure mechanisms are discussed. The threshold rates to differentiate drained, partially drained, and undrained conditions are quantified. Anchor penetrations are also simulated under the same conditions. Comparisons with uplift cases indicate that the force–displacement curves in uplift and penetration are only identical in undrained conditions.
{"title":"Numerical investigations into the influence of drainage on the uplift behaviour of plate anchors in clay","authors":"Maozhu Peng , Zhen-Yu Yin , Zefeng Zhou","doi":"10.1016/j.compgeo.2025.107128","DOIUrl":"10.1016/j.compgeo.2025.107128","url":null,"abstract":"<div><div>This paper presents one of the first systematic finite element investigations into the rate effect during plate anchor uplift in clay, with a focus on the influence of soil drainage. Some numerical difficulties in simulating partially/fully drained anchor uplift are first discussed, and a method to address these issues is proposed. The most striking benefit of the proposed method is that it allows the hydro-mechanical coupled anchor-soil interaction to be efficiently simulated, with improved numerical stability, by simply using a non-detachable anchor-soil interface. After validation, the proposed method is used to simulate anchor uplift at different embedment and uplift rates, spanning from drained to fully undrained conditions. The influence of uplift rate on the force–displacement relations and failure mechanisms are discussed. The threshold rates to differentiate drained, partially drained, and undrained conditions are quantified. Anchor penetrations are also simulated under the same conditions. Comparisons with uplift cases indicate that the force–displacement curves in uplift and penetration are only identical in undrained conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107128"},"PeriodicalIF":5.3,"publicationDate":"2025-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143324761","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}
This paper proposes a novel numerical simulation approach for fluid–solid coupling in fractured reservoirs, incorporating the mechanical properties of proppants into the coupling framework. By integrating the Displacement Discontinuity Method (DDM) with a three-dimensional Embedded Discrete Fracture Model (3D EDFM), the proposed method enables coupled simulations of geomechanics and fluid flow. Unlike conventional approaches, this method considers the time-dependent evolution of fracture aperture and accurately captures the mechanical effects of proppants during fracture closure. A comprehensive analysis is conducted to investigate the impacts on fracture dynamics and fluid flow behavior. The model’s reliability is validated through comparisons with analytical solutions and the Extended Finite Element Method (XFEM). Results demonstrate that the time-dependent evolution of fracture apertures significantly influences fluid flow in fractured reservoirs. Specifically, during the production stage, fracture closure results in a sharp decline in oil production rates. Larger proppant-supported apertures effectively mitigate fracture closure, sustain high fracture conductivity, and prolong the production period. Additionally, low-porosity reservoirs are shown to be more sensitive to rock deformation, leading to substantial changes in flow parameters during production. Thus, incorporating fluid–solid coupling is crucial for accurate modeling in low-porosity reservoirs. This study provides valuable theoretical insights for the development of unconventional resources.
{"title":"A novel fluid-solid coupling method for fractured reservoirs: 3D DDM-EDFM integration with proppant mechanics","authors":"Luoyi Huang , Wentao Zhan , Hui Zhao , Guanglong Sheng","doi":"10.1016/j.compgeo.2025.107127","DOIUrl":"10.1016/j.compgeo.2025.107127","url":null,"abstract":"<div><div>This paper proposes a novel numerical simulation approach for fluid–solid coupling in fractured reservoirs, incorporating the mechanical properties of proppants into the coupling framework. By integrating the Displacement Discontinuity Method (DDM) with a three-dimensional Embedded Discrete Fracture Model (3D EDFM), the proposed method enables coupled simulations of geomechanics and fluid flow. Unlike conventional approaches, this method considers the time-dependent evolution of fracture aperture and accurately captures the mechanical effects of proppants during fracture closure. A comprehensive analysis is conducted to investigate the impacts on fracture dynamics and fluid flow behavior. The model’s reliability is validated through comparisons with analytical solutions and the Extended Finite Element Method (XFEM). Results demonstrate that the time-dependent evolution of fracture apertures significantly influences fluid flow in fractured reservoirs. Specifically, during the production stage, fracture closure results in a sharp decline in oil production rates. Larger proppant-supported apertures effectively mitigate fracture closure, sustain high fracture conductivity, and prolong the production period. Additionally, low-porosity reservoirs are shown to be more sensitive to rock deformation, leading to substantial changes in flow parameters during production. Thus, incorporating fluid–solid coupling is crucial for accurate modeling in low-porosity reservoirs. This study provides valuable theoretical insights for the development of unconventional resources.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107127"},"PeriodicalIF":5.3,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143324444","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}
Pub Date : 2025-02-06DOI: 10.1016/j.compgeo.2025.107129
Weiwei Niu , Yuan-Yuan Zheng , Zhen-Yu Yin , Chi Yao , Pengchang Wei
The mechanical behavior of expansive soil in geotechnical engineering is significantly sensitive to loading rates, hydration, confining pressure, etc., where most engineering problems are attributed to the existence of montmorillonite in expansive soil. Here, the hydration, confining pressure, and loading rate effect on the mechanical behavior of montmorillonite were investigated through the triaxial tests and molecular dynamics (MD) simulation method, revealing their fundamental mechanism between the microscale and macroscale. The average basal spacing of hydrated montmorillonite system, the diffusion coefficient and density distribution of interlayer water molecules were calculated for the verification of MD model. The experimental results indicated that the stress–strain relationship of montmorillonite was the strain-hardening type. The failure stress did not increase monotonously with the increase in loading rate, and there were two obvious critical points. The failure stress of the soil sample increased with the increase of the confining pressure, and the decrease of the water content, where their fundamental mechanism between microscale and macroscale were adequately discussed. Furthermore, the stress–strain response, total energy evolution, deformation evolution of atomistic structure, and broken bonds evolution were analyzed to deeply understand the fundamental deformation mechanism at the microscale. The multi-scale studies could effectively examine the macroscopic mechanical behavior of expansive soil and elucidate its microscopic mechanisms.
{"title":"Multiscale mechanical behavior of hydrated expansive soil: Insights from experimental and MD study","authors":"Weiwei Niu , Yuan-Yuan Zheng , Zhen-Yu Yin , Chi Yao , Pengchang Wei","doi":"10.1016/j.compgeo.2025.107129","DOIUrl":"10.1016/j.compgeo.2025.107129","url":null,"abstract":"<div><div>The mechanical behavior of expansive soil in geotechnical engineering is significantly sensitive to loading rates, hydration, confining pressure, etc., where most engineering problems are attributed to the existence of montmorillonite in expansive soil. Here, the hydration, confining pressure, and loading rate effect on the mechanical behavior of montmorillonite were investigated through the triaxial tests and molecular dynamics (MD) simulation method, revealing their fundamental mechanism between the microscale and macroscale. The average basal spacing of hydrated montmorillonite system, the diffusion coefficient and density distribution of interlayer water molecules were calculated for the verification of MD model. The experimental results indicated that the stress–strain relationship of montmorillonite was the strain-hardening type. The failure stress did not increase monotonously with the increase in loading rate, and there were two obvious critical points. The failure stress of the soil sample increased with the increase of the confining pressure, and the decrease of the water content, where their fundamental mechanism between microscale and macroscale were adequately discussed. Furthermore, the stress–strain response, total energy evolution, deformation evolution of atomistic structure, and broken bonds evolution were analyzed to deeply understand the fundamental deformation mechanism at the microscale. The multi-scale studies could effectively examine the macroscopic mechanical behavior of expansive soil and elucidate its microscopic mechanisms.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107129"},"PeriodicalIF":5.3,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143325092","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}
Pub Date : 2025-02-05DOI: 10.1016/j.compgeo.2025.107110
Xiao-Xuan Chen , Pin Zhang , Zhen-Yu Yin
Physics-informed learning has emerged as a promising approach for solving forward and inverse partial differential equations in engineering practice, but selecting an optimal loss function remains unclear and parameter identification for inverse analysis lacks efficiency. Meanwhile, their values for engineering-scale elastoplastic problems have not been deeply investigated. In this research, a comprehensive comparison between the strong-form collocation point method (CPM) and the deep Ritz method (DRM) based loss functions for both forward and inverse analysis is conducted, and a novel exponential acceleration method is proposed to enlarge the search space of unknown parameters for inverse analysis. By applying these methods to linear elasticity and elastoplasticity footing cases, we found that physics-informed learning equipped with DRM-based loss functions shows more excellent accuracy in forwardly computing displacement but poor accuracy in predicting strain and stress. Physics-informed learning with CPM-based loss functions shows more excellent performance in inverse analysis than their forward-solving ability. The exponential acceleration method largely enhances the efficiency of inverse analysis without sacrificing accuracy. These new findings inspire the future application of physics-informed learning to engineering-scale elastoplastic problems.
{"title":"A Comprehensive Investigation of Physics-Informed Learning in Forward and Inverse Analysis of Elastic and Elastoplastic Footing","authors":"Xiao-Xuan Chen , Pin Zhang , Zhen-Yu Yin","doi":"10.1016/j.compgeo.2025.107110","DOIUrl":"10.1016/j.compgeo.2025.107110","url":null,"abstract":"<div><div>Physics-informed learning has emerged as a promising approach for solving forward and inverse partial differential equations in engineering practice, but selecting an optimal loss function remains unclear and parameter identification for inverse analysis lacks efficiency. Meanwhile, their values for engineering-scale elastoplastic problems have not been deeply investigated. In this research, a comprehensive comparison between the strong-form collocation point method (CPM) and the deep Ritz method (DRM) based loss functions for both forward and inverse analysis is conducted, and a novel exponential acceleration method is proposed to enlarge the search space of unknown parameters for inverse analysis. By applying these methods to linear elasticity and elastoplasticity footing cases, we found that physics-informed learning equipped with DRM-based loss functions shows more excellent accuracy in forwardly computing displacement but poor accuracy in predicting strain and stress. Physics-informed learning with CPM-based loss functions shows more excellent performance in inverse analysis than their forward-solving ability. The exponential acceleration method largely enhances the efficiency of inverse analysis without sacrificing accuracy. These new findings inspire the future application of physics-informed learning to engineering-scale elastoplastic problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107110"},"PeriodicalIF":5.3,"publicationDate":"2025-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143324679","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}
Pub Date : 2025-02-05DOI: 10.1016/j.compgeo.2025.107095
Sinem Bozkurt, Ayman Abed, Minna Karstunen
The stabilisation of deep excavations using columnar inclusions introduces three-dimensional (3D) complexities into design calculations. Often the intricate nature of the wall-soil-column interaction is simulated using simplified averaging techniques based on elasticity theory. This paper introduces a numerical technique to accurately describe the nonlinear elastoplastic response of both stabilised and natural soft clay for the case of excavations stabilised with deep-mixed columns. The technique allows mapping the 3D problem into a plane-strain (2D) counterpart by replacing the composite material made of natural clay and a region of overlapping deep-mixed columns with an equivalent homogenised material. The stress–strain response of the homogenised continuum is computed with a volume averaging technique (VAT) based on the volume fraction of each component (i.e. clay and column). The technique is implemented into a 2D finite element code enabling an effective representation of the behaviour of each constituent represented by advanced elastoplastic material models. After presenting the theoretical background and implementation procedure, the proposed method was verified against the results from the 3D calculations. The technique emerges as an efficient tool for the numerical analysis of stabilised deep excavations since it allows for plane strain analysis to yield results akin to computationally expensive 3D analysis. Thus, the method can significantly reduce the computational costs and can facilitate the easier incorporation of sensitivity studies.
{"title":"Homogenisation method for braced excavations stabilised with deep-mixed columns","authors":"Sinem Bozkurt, Ayman Abed, Minna Karstunen","doi":"10.1016/j.compgeo.2025.107095","DOIUrl":"10.1016/j.compgeo.2025.107095","url":null,"abstract":"<div><div>The stabilisation of deep excavations using columnar inclusions introduces three-dimensional (3D) complexities into design calculations. Often the intricate nature of the wall-soil-column interaction is simulated using simplified averaging techniques based on elasticity theory. This paper introduces a numerical technique to accurately describe the nonlinear elastoplastic response of both stabilised and natural soft clay for the case of excavations stabilised with deep-mixed columns. The technique allows mapping the 3D problem into a plane-strain (2D) counterpart by replacing the composite material made of natural clay and a region of overlapping deep-mixed columns with an equivalent homogenised material. The stress–strain response of the homogenised continuum is computed with a volume averaging technique (VAT) based on the volume fraction of each component (i.e. clay and column). The technique is implemented into a 2D finite element code enabling an effective representation of the behaviour of each constituent represented by advanced elastoplastic material models. After presenting the theoretical background and implementation procedure, the proposed method was verified against the results from the 3D calculations. The technique emerges as an efficient tool for the numerical analysis of stabilised deep excavations since it allows for plane strain analysis to yield results akin to computationally expensive 3D analysis. Thus, the method can significantly reduce the computational costs and can facilitate the easier incorporation of sensitivity studies.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107095"},"PeriodicalIF":5.3,"publicationDate":"2025-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143395383","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}
Pub Date : 2025-02-05DOI: 10.1016/j.compgeo.2025.107124
Xin-Dong Wei, Gao-Feng Zhao
Lattice Spring Methods (LSM) have gained increasing attention in recent years due to their effectiveness in handling fracturing and nonlinear problems. However, their ability to address general contact issues in geomechanics remains limited. Directly handling contact using discrete geometry often results in local locking, a common challenge in particle-based methods like DEM and LSM. This work introduces a level set approach to overcome this limitation. Block and granular geometries are represented by a level set function, which is also employed to manage contact interactions between blocks and granules. A specialized contact algorithm is developed for the lattice model using a signed distance function (SDF) that accounts for lattice discretization. The new LS-4D-LSM demonstrates its capability to solve complex contact problems through various numerical examples, including both deformation and fracturing of blocks. Given that contact plays a key role in many geomechanics tests, incorporating realistic contact behavior into these analyses leads to a more accurate fit with experimentally observed data. This suggests that the digital twin concept could be a powerful tool for analyzing experimental data in geomechanics. The proposed method is also applicable to other numerical approaches, such as DEM, enhancing their ability to handle deformable and fragmentable blocks or granules.
{"title":"Four-dimensional lattice spring method using a level set approach for contact problems in geomechanics","authors":"Xin-Dong Wei, Gao-Feng Zhao","doi":"10.1016/j.compgeo.2025.107124","DOIUrl":"10.1016/j.compgeo.2025.107124","url":null,"abstract":"<div><div>Lattice Spring Methods (LSM) have gained increasing attention in recent years due to their effectiveness in handling fracturing and nonlinear problems. However, their ability to address general contact issues in geomechanics remains limited. Directly handling contact using discrete geometry often results in local locking, a common challenge in particle-based methods like DEM and LSM. This work introduces a level set approach to overcome this limitation. Block and granular geometries are represented by a level set function, which is also employed to manage contact interactions between blocks and granules. A specialized contact algorithm is developed for the lattice model using a signed distance function (SDF) that accounts for lattice discretization. The new LS-4D-LSM demonstrates its capability to solve complex contact problems through various numerical examples, including both deformation and fracturing of blocks. Given that contact plays a key role in many geomechanics tests, incorporating realistic contact behavior into these analyses leads to a more accurate fit with experimentally observed data. This suggests that the digital twin concept could be a powerful tool for analyzing experimental data in geomechanics. The proposed method is also applicable to other numerical approaches, such as DEM, enhancing their ability to handle deformable and fragmentable blocks or granules.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"181 ","pages":"Article 107124"},"PeriodicalIF":5.3,"publicationDate":"2025-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143357316","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}
Pub Date : 2025-02-04DOI: 10.1016/j.compgeo.2025.107099
Linfeng Zhang , Guorong Wang , Zhiyuan Li , Lin Zhong , Qiang Fu , Jiwei Wu , Mengdi Wu , Dongmei Liang , Yuhang Zheng
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.
{"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}
Pub Date : 2025-02-04DOI: 10.1016/j.compgeo.2025.107113
Hao Chen, Shiwei Zhao, Jidong Zhao
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.
{"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}
Pub Date : 2025-02-01DOI: 10.1016/j.compgeo.2025.107114
Weilie Zou , Qiuyang Pei , Xilin Xia , Zhong Han
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.
{"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}
Pub Date : 2025-02-01DOI: 10.1016/j.compgeo.2025.107106
Penghao Zhang, Yunxuan Cui, Kurt Douglas, Chongmin Song, Adrian R. Russell
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.
{"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}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号