Computers and Geotechnics最新文献_第4页

Coupled mechanisms of load transfer and fracture in bonded granular materials using CT image-based discrete element method 基于CT图像离散元法的粘结颗粒材料载荷传递与断裂耦合机制研究

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-01-23 DOI: 10.1016/j.compgeo.2025.107891

Hu Yang , Yiik Diew Wong , Liyan Shan , Lingwen Li

Bonded granular materials (BGMs) consist of the skeleton of granular aggregate particles and a cementitious agent between the particles. This study presents a damage analysis framework for BGMs using a computed tomography (CT) image-based discrete element method. X-ray CT scanning and image processing technology are employed to construct in situ discrete element models, categorizing particle bonding into interface bonding and mortar paste bonding. The force chains are decomposed into compressive, tensile, and shear chains, and their evolution is analyzed using complex network theory. The study investigates the coupled mechanisms between force chain evolution and bonding failure in BGMs under loading. Key findings include the critical role of fracture resistance differences between interface bonding and mortar paste bonding in determining BGM strength. The bonding-failure rate, particularly at peak load, reliably indicates material strength, with interface bonding failure surpassing mortar paste bonding failure. Force chain evolution shows rapid concentration in the loading zone, forming a stable backbone structure despite post-peak degradation. Shear and tensile force chain networks evolve with crack propagation, and force transmission paths are reorganized in the final phase. The critical points of force chain network parameters align with macroscopic load response and fracture evolution, thereby offering insights into coupled load transfer and fracture mechanisms. The proposed framework not only advances the understanding of force chain dynamics but also supports damage prediction and structural optimization by providing a comprehensive tool to track materials behavior under loading.

粘结颗粒材料（BGMs）由颗粒骨料骨架和颗粒之间的胶凝剂组成。本研究提出了一种基于计算机断层扫描（CT）图像的离散元方法的bgm损伤分析框架。采用x射线CT扫描和图像处理技术构建原位离散单元模型，将颗粒粘结分为界面粘结和砂浆粘结。将力链分解为压缩链、拉伸链和剪切链，并利用复杂网络理论分析了它们的演化过程。研究了应力链演化与粘结破坏的耦合机制。主要发现包括界面粘结和砂浆粘结之间的抗断裂性差异在决定BGM强度中的关键作用。粘结失败率，特别是在峰值荷载下，可靠地反映了材料的强度，界面粘结破坏超过砂浆粘结破坏。力链演化在加载区迅速集中，峰后退化形成稳定的骨干结构。剪切和拉伸力链网络随着裂纹扩展而演化，力传递路径在最后阶段重新组织。力链网络参数的临界点与宏观载荷响应和断裂演化一致，从而深入了解载荷耦合传递和断裂机制。提出的框架不仅促进了对力链动力学的理解，而且通过提供一个全面的工具来跟踪材料在载荷下的行为，支持损伤预测和结构优化。

{"title":"Coupled mechanisms of load transfer and fracture in bonded granular materials using CT image-based discrete element method","authors":"Hu Yang , Yiik Diew Wong , Liyan Shan , Lingwen Li","doi":"10.1016/j.compgeo.2025.107891","DOIUrl":"10.1016/j.compgeo.2025.107891","url":null,"abstract":"<div><div>Bonded granular materials (BGMs) consist of the skeleton of granular aggregate particles and a cementitious agent between the particles. This study presents a damage analysis framework for BGMs using a computed tomography (CT) image-based discrete element method. X-ray CT scanning and image processing technology are employed to construct in situ discrete element models, categorizing particle bonding into interface bonding and mortar paste bonding. The force chains are decomposed into compressive, tensile, and shear chains, and their evolution is analyzed using complex network theory. The study investigates the coupled mechanisms between force chain evolution and bonding failure in BGMs under loading. Key findings include the critical role of fracture resistance differences between interface bonding and mortar paste bonding in determining BGM strength. The bonding-failure rate, particularly at peak load, reliably indicates material strength, with interface bonding failure surpassing mortar paste bonding failure. Force chain evolution shows rapid concentration in the loading zone, forming a stable backbone structure despite post-peak degradation. Shear and tensile force chain networks evolve with crack propagation, and force transmission paths are reorganized in the final phase. The critical points of force chain network parameters align with macroscopic load response and fracture evolution, thereby offering insights into coupled load transfer and fracture mechanisms. The proposed framework not only advances the understanding of force chain dynamics but also supports damage prediction and structural optimization by providing a comprehensive tool to track materials behavior under loading.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107891"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026151","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}

引用次数: 0

Tunnel excavation and swelling analysis of expansive bedrock with multiphysics elasto-plastic model capable of describing different swelling behavior due to exchangeable cation species 基于多物理场弹塑性模型的膨胀基岩隧道开挖及溶胀分析

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-01-30 DOI: 10.1016/j.compgeo.2026.107956

Keitaro Hoshi, Shotaro Yamada, Yuta Abe, Takashi Kyoya

Swelling of smectite-bearing bedrock can cause severe tunnel deformation, depending on the type of exchangeable cation present in the interlayer structure. This study proposes an extended expansive bedrock model capable of capturing distinct swelling behaviors induced by different cation species. The model incorporates a double-layer repulsive force, formulated based on Stern theory, into a previously developed finite elastoplastic framework. Finite element analyses of tunnel excavation and subsequent swelling were performed using the proposed model. The results indicate that yielding of the bedrock skeleton acts as a trigger for accelerated swelling deformation, and that the swelling behavior is strongly influenced by the type of exchangeable cation: in sodium-type smectite, pronounced swelling occurred primarily at the tunnel invert, whereas calcium- and potassium-type smectites exhibited only minor expansion. The analysis also investigated the mechanical interaction between the expansive bedrock and an invert concrete layer. Under the assumed conditions, compressive axial stresses exceeding 20 MPa developed in the invert, suggesting that the swelling pressure can surpass the compressive strength of ordinary unreinforced concrete. These findings elucidate the fundamental mechanism of tunnel invert deformation, highlighting the distinct swelling behaviors associated with various exchangeable cation species, clarifying the multiscale and multiphysics interactions between electrochemical processes in the interlaminar region and the elastoplastic response of the surrounding rock mass, and quantitatively demonstrating the mitigating effect of the invert on swelling-induced tunnel deformation.

含蒙脱石基岩的膨胀会导致严重的隧道变形，这取决于层间结构中存在的交换阳离子的类型。本研究提出了一种扩展的膨胀基岩模型，该模型能够捕捉不同阳离子引起的不同膨胀行为。该模型将基于斯特恩理论的双层斥力纳入到先前开发的有限弹塑性框架中。利用该模型对隧道开挖及后续膨胀进行了有限元分析。结果表明，基岩骨架的屈服是加速膨胀变形的触发因素，并且膨胀行为受交换阳离子类型的强烈影响：在钠型蒙脱石中，明显的膨胀主要发生在隧道仰拱处，而钙和钾型蒙脱石仅表现出轻微的膨胀。分析还探讨了膨胀基岩与倒拱混凝土层之间的力学相互作用。在假设条件下，仰拱出现了超过20 MPa的轴压应力，表明膨胀压力可以超过普通无筋混凝土的抗压强度。这些发现阐明了隧道反转变形的基本机制，突出了与不同交换阳离子相关的独特膨胀行为，阐明了层间区电化学过程与围岩弹塑性响应之间的多尺度和多物理场相互作用，并定量地证明了反转对膨胀引起的隧道变形的缓解作用。

{"title":"Tunnel excavation and swelling analysis of expansive bedrock with multiphysics elasto-plastic model capable of describing different swelling behavior due to exchangeable cation species","authors":"Keitaro Hoshi, Shotaro Yamada, Yuta Abe, Takashi Kyoya","doi":"10.1016/j.compgeo.2026.107956","DOIUrl":"10.1016/j.compgeo.2026.107956","url":null,"abstract":"<div><div>Swelling of smectite-bearing bedrock can cause severe tunnel deformation, depending on the type of exchangeable cation present in the interlayer structure. This study proposes an extended expansive bedrock model capable of capturing distinct swelling behaviors induced by different cation species. The model incorporates a double-layer repulsive force, formulated based on Stern theory, into a previously developed finite elastoplastic framework. Finite element analyses of tunnel excavation and subsequent swelling were performed using the proposed model. The results indicate that yielding of the bedrock skeleton acts as a trigger for accelerated swelling deformation, and that the swelling behavior is strongly influenced by the type of exchangeable cation: in sodium-type smectite, pronounced swelling occurred primarily at the tunnel invert, whereas calcium- and potassium-type smectites exhibited only minor expansion. The analysis also investigated the mechanical interaction between the expansive bedrock and an invert concrete layer. Under the assumed conditions, compressive axial stresses exceeding 20 MPa developed in the invert, suggesting that the swelling pressure can surpass the compressive strength of ordinary unreinforced concrete. These findings elucidate the fundamental mechanism of tunnel invert deformation, highlighting the distinct swelling behaviors associated with various exchangeable cation species, clarifying the multiscale and multiphysics interactions between electrochemical processes in the interlaminar region and the elastoplastic response of the surrounding rock mass, and quantitatively demonstrating the mitigating effect of the invert on swelling-induced tunnel deformation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107956"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080666","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}

引用次数: 0

Time-dependent deformations in deep tunnels: Insights into uncertainty and variability of rheological behavior 深隧道中随时间变化的变形：流变行为的不确定性和可变性

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-01-31 DOI: 10.1016/j.compgeo.2026.107942

Milad Zaheri , Pierpaolo Oreste , Masoud Ranjbarnia , Elham Mahmoudi

The uncertainty in determining rock mass properties significantly impacts tunnel stability. Additionally, squeezing conditions worsen tunnel stability, causing the tunnel to gradually converge over time. This paper addresses this issue and investigates the long-term behavior of deep tunnels using both visco-elastic and visco-elasto-plastic models. This study also includes risk-based analyses to offer a quantitative tool for engineering decision-making. Initially, an analytical method is introduced to calculate tunnel convergence in a visco-elastic rock mass. The uncertainty of key parameters that significantly affect tunnel behavior is also considered. Using MATLAB, the probability distributions of tunnel wall deformations over time are determined. The results indicate that, except in one case, the long-term tunnel convergence follows a right-skewed Gamma distribution, especially with a low GSI in both visco-elastic and visco-elasto-plastic models. This suggests that deterministic methods may not be reliable for ensuring the safety of long-term tunnel designs.

岩体性质的不确定性对隧道的稳定性影响很大。此外，挤压条件使隧道稳定性恶化，导致隧道随时间逐渐收敛。本文解决了这一问题，并采用粘弹性和粘弹塑性模型研究了深埋隧道的长期行为。本研究还包括基于风险的分析，为工程决策提供定量工具。首先介绍了粘弹性岩体中隧道收敛的解析计算方法。还考虑了对隧道性能有重要影响的关键参数的不确定性。利用MATLAB软件，确定了隧道壁变形随时间的概率分布。结果表明，除了一种情况外，长期隧道收敛遵循右偏Gamma分布，特别是粘弹性和粘弹塑性模型的低GSI。这表明，确定性方法可能不可靠，以确保长期隧道设计的安全。

{"title":"Time-dependent deformations in deep tunnels: Insights into uncertainty and variability of rheological behavior","authors":"Milad Zaheri , Pierpaolo Oreste , Masoud Ranjbarnia , Elham Mahmoudi","doi":"10.1016/j.compgeo.2026.107942","DOIUrl":"10.1016/j.compgeo.2026.107942","url":null,"abstract":"<div><div>The uncertainty in determining rock mass properties significantly impacts tunnel stability. Additionally, squeezing conditions worsen tunnel stability, causing the tunnel to gradually converge over time. This paper addresses this issue and investigates the long-term behavior of deep tunnels using both visco-elastic and visco-elasto-plastic models. This study also includes risk-based analyses to offer a quantitative tool for engineering decision-making. Initially, an analytical method is introduced to calculate tunnel convergence in a visco-elastic rock mass. The uncertainty of key parameters that significantly affect tunnel behavior is also considered. Using MATLAB, the probability distributions of tunnel wall deformations over time are determined. The results indicate that, except in one case, the long-term tunnel convergence follows a right-skewed Gamma distribution, especially with a low GSI in both visco-elastic and visco-elasto-plastic models. This suggests that deterministic methods may not be reliable for ensuring the safety of long-term tunnel designs.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107942"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080663","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}

引用次数: 0

The performance of quadratic finite-discrete element method (qFDEM) and its potential advantages 二次有限离散元法（qFDEM）的性能及其潜在优势

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-01-28 DOI: 10.1016/j.compgeo.2026.107925

Zhonghao Li , Xiaofeng Li , Haibo Li , Qi Zhao , Giovanni Grasselli

The combined finite–discrete element method (FDEM) has proven effective for simulating crack initiation, propagation, and coalescence in brittle solids. However, existing FDEM frameworks remain limited to constant-strain elements, leading to restricted capability in representing complex stress fields, pronounced sensitivity to shear and volumetric locking, and a strong tendency toward numerical dispersion in dynamic problems. To overcome these limitations, this study develops a high order element-based framework incorporating a novel quadratic cohesive element to enhance model accuracy and continuity. The proposed quadratic cohesive element ensures uniform traction distribution along edges, avoiding the mid-node stress concentrations that typically lead to mesh incompatibility and artificial strength reduction. Three quasi-static loading tests and one wave propagation test are performed to compare quadratic and linear models. The results show that the quadratic model consistently outperforms the linear counterpart in stress path, crack propagation, and mitigating numerical dispersion. In quasi-static loading, the new quadratic model exhibits a lower error in stress, predicts a more precise crack initiation load, and provides more reliable crack path predictions compared with previous models. In dynamic conditions, it can effectively mitigate the numerical dispersion of high-frequency wave components that low-order elements struggle with and provide more stable wave propagation simulations. Moreover, the quadratic elements FDEM framework offers an economical alternative for enhancing the fidelity of crack simulations: compared to mesh refinement, quadratic elements achieve comparable accuracy in crack initiation load prediction with only 50–60% of the computational cost.

有限-离散元组合方法（FDEM）已被证明是模拟脆性固体裂纹萌生、扩展和合并的有效方法。然而，现有的FDEM框架仍然局限于恒应变单元，导致表征复杂应力场的能力有限，对剪切和体积锁定的敏感性明显，并且在动态问题中有很强的数值离散倾向。为了克服这些限制，本研究开发了一个基于高阶元素的框架，其中包含了一个新的二次内聚元素，以提高模型的准确性和连续性。所提出的二次内聚单元确保沿边缘均匀的牵引力分布，避免了通常导致网格不兼容和人工强度降低的中间节点应力集中。进行了3次准静态加载试验和1次波传播试验，对二次模型和线性模型进行了比较。结果表明，二次型模型在应力路径、裂纹扩展和减轻数值弥散方面均优于线性模型。在准静态加载下，新二次元模型的应力误差更小，预测裂纹起裂载荷更精确，裂纹路径预测更可靠。在动态条件下，它可以有效地缓解低阶元与高频波分量的数值频散，提供更稳定的波传播模拟。此外，二次元FDEM框架为提高裂纹模拟的保真度提供了一种经济的替代方案：与网格细化相比，二次元在裂纹起裂载荷预测方面达到了相当的精度，而计算成本仅为50-60%。

{"title":"The performance of quadratic finite-discrete element method (qFDEM) and its potential advantages","authors":"Zhonghao Li , Xiaofeng Li , Haibo Li , Qi Zhao , Giovanni Grasselli","doi":"10.1016/j.compgeo.2026.107925","DOIUrl":"10.1016/j.compgeo.2026.107925","url":null,"abstract":"<div><div>The combined finite–discrete element method (FDEM) has proven effective for simulating crack initiation, propagation, and coalescence in brittle solids. However, existing FDEM frameworks remain limited to constant-strain elements, leading to restricted capability in representing complex stress fields, pronounced sensitivity to shear and volumetric locking, and a strong tendency toward numerical dispersion in dynamic problems. To overcome these limitations, this study develops a high order element-based framework incorporating a novel quadratic cohesive element to enhance model accuracy and continuity. The proposed quadratic cohesive element ensures uniform traction distribution along edges, avoiding the mid-node stress concentrations that typically lead to mesh incompatibility and artificial strength reduction. Three quasi-static loading tests and one wave propagation test are performed to compare quadratic and linear models. The results show that the quadratic model consistently outperforms the linear counterpart in stress path, crack propagation, and mitigating numerical dispersion. In quasi-static loading, the new quadratic model exhibits a lower error in stress, predicts a more precise crack initiation load, and provides more reliable crack path predictions compared with previous models. In dynamic conditions, it can effectively mitigate the numerical dispersion of high-frequency wave components that low-order elements struggle with and provide more stable wave propagation simulations. Moreover, the quadratic elements FDEM framework offers an economical alternative for enhancing the fidelity of crack simulations: compared to mesh refinement, quadratic elements achieve comparable accuracy in crack initiation load prediction with only 50–60% of the computational cost.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107925"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080670","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}

引用次数: 0

Hydro-mechanical–chemical modelling of solute–consolidation coupling in near-saturated soils with entrapped bubbles 含气泡近饱和土壤溶质-固结耦合的水-力学-化学模拟

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-12 DOI: 10.1016/j.compgeo.2026.107967

Bolin Wang, Dong-Sheng Jeng

Long-term solute transport in deformable, low-permeability porous media plays a crucial role in various geo-environmental and environmental problems. This study presents a coupled three-dimensional (3D) hydro-mechanical–chemical (HMC) model tailored for unsaturated conditions. Two distinct chemical mechanisms are incorporated: osmotic pressure governing chemical–hydraulic feedback, and chemically induced strain contributing to the mechanical response. Numerical simulations show that osmotic effects dominate solute migration by amplifying advective transport and modifying pore-pressure variation, whereas chemically induced deformation plays a secondary role under small strain condition. Furthermore, the pore-radius-dependent interfacial effect primarily modulates the early consolidation response by altering the pore-pressure field, with smaller characteristic pore radii producing steeper hydraulic gradients and stronger advective solute transport. Sensitivity analyses identify threshold values of solute flux, mechanical loading and osmotic efficiency beyond which chemical effects become significant, while under milder conditions conventional hydro–mechanical (HM) models provide reliable approximations of the full HMC responses. The proposed framework provides a robust basis for predicting long-term chemo-hydraulic interactions in unsaturated soils subjected to coupled environmental loading.

溶质在可变形低渗透多孔介质中的长期运移在各种地质环境和环境问题中起着至关重要的作用。本研究提出了一个为非饱和条件量身定制的耦合三维（3D）水-机械-化学（HMC）模型。两种不同的化学机制被纳入其中：控制化学-液压反馈的渗透压，以及促进机械响应的化学诱导应变。数值模拟表明，渗透效应通过放大平流输运和改变孔隙压力变化主导溶质迁移，而在小应变条件下，化学诱导变形起次要作用。此外，孔隙半径相关的界面效应主要通过改变孔隙压力场来调节早期固结响应，较小的特征孔隙半径产生更陡的水力梯度和更强的平流溶质输运。敏感性分析确定溶质通量、机械负荷和渗透效率的阈值，超过这些阈值，化学效应就会变得显著，而在较温和的条件下，传统的水-机械（HM）模型提供了完整的HMC响应的可靠近似值。所提出的框架为预测耦合环境荷载作用下非饱和土的长期化学-水力相互作用提供了坚实的基础。

{"title":"Hydro-mechanical–chemical modelling of solute–consolidation coupling in near-saturated soils with entrapped bubbles","authors":"Bolin Wang, Dong-Sheng Jeng","doi":"10.1016/j.compgeo.2026.107967","DOIUrl":"10.1016/j.compgeo.2026.107967","url":null,"abstract":"<div><div>Long-term solute transport in deformable, low-permeability porous media plays a crucial role in various geo-environmental and environmental problems. This study presents a coupled three-dimensional (3D) hydro-mechanical–chemical (HMC) model tailored for unsaturated conditions. Two distinct chemical mechanisms are incorporated: osmotic pressure governing chemical–hydraulic feedback, and chemically induced strain contributing to the mechanical response. Numerical simulations show that osmotic effects dominate solute migration by amplifying advective transport and modifying pore-pressure variation, whereas chemically induced deformation plays a secondary role under small strain condition. Furthermore, the pore-radius-dependent interfacial effect primarily modulates the early consolidation response by altering the pore-pressure field, with smaller characteristic pore radii producing steeper hydraulic gradients and stronger advective solute transport. Sensitivity analyses identify threshold values of solute flux, mechanical loading and osmotic efficiency beyond which chemical effects become significant, while under milder conditions conventional hydro–mechanical (HM) models provide reliable approximations of the full HMC responses. The proposed framework provides a robust basis for predicting long-term chemo-hydraulic interactions in unsaturated soils subjected to coupled environmental loading.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107967"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146174286","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}

引用次数: 0

Modeling non-Newtonian fluid–solid flows containing non-spherical particles by the SPH-DEM coupling model 用SPH-DEM耦合模型模拟含非球形颗粒的非牛顿流固流动

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-12 DOI: 10.1016/j.compgeo.2026.107990

Jiahang Du , Yihang Yu , Yuqi Guo , Huaqing Ma , Yongzhi Zhao

Given the ubiquity and significance of fluid–solid systems composed of non-Newtonian fluids and non-spherical particles in nature and industry, the deep investigation of their interaction dynamics is crucial for geological studies and engineering applications. In this paper, a resolved coupling model based on Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) is proposed to simulate the interaction between non-Newtonian fluids and non-spherical particles. The non-Newtonian fluid is modelled using SPH with the power-law model, while the dynamics of the non-spherical particles is simulated using DEM incorporated with the super-ellipsoid and polyhedral particle models. A modified boundary repulsive force model, combined with the fixed ghost particle method, is developed to handle boundary conditions. For super-ellipsoid and polyhedron boundaries, an improved approach is proposed to generate uniformly distributed, body-fitted boundary dummy particles on their surfaces, thereby ensuring computational accuracy near boundaries and in fluid–solid interactions. To rigorously validate the framework, four benchmark cases are conducted, and the resulting numerical simulations are subsequently compared with experiments in this study. The SPH and DEM methods are first independently validated using a 3D non-Newtonian dam-break and a dry granular column collapse, respectively. The coupled framework is then assessed through two fluid–solid interaction cases: cube entry into a non-Newtonian fluid and complex dam-break scenarios involving mixtures of non-Newtonian fluids and non-spherical particles. The simulated flow behaviors (e.g., the cube penetration depth and downstream propagation of fluid–solid mixture) agree well with the corresponding experiments, validating the effectiveness of the proposed boundary treatment method and the accuracy of the framework to simulate the fluid–solid interaction systems involving non-Newtonian fluids and non-spherical particles.

鉴于由非牛顿流体和非球形颗粒组成的流固系统在自然界和工业中的普遍存在和重要意义，深入研究它们的相互作用动力学对于地质研究和工程应用至关重要。本文提出了一种基于光滑粒子流体力学（SPH）和离散元法（DEM）的解析耦合模型，用于模拟非牛顿流体与非球形粒子之间的相互作用。非牛顿流体采用幂律模型SPH建模，非球形颗粒采用超椭球和多面体粒子模型相结合的DEM模拟动力学。提出了一种修正的边界斥力模型，结合固定鬼粒子法处理边界条件。对于超椭球体和多面体边界，提出了一种改进的方法，在其表面生成均匀分布的拟体边界虚拟粒子，从而保证了边界附近和流固相互作用的计算精度。为了严格验证该框架，进行了四个基准案例，并将所得数值模拟结果与本研究的实验进行了比较。SPH和DEM方法首先分别使用三维非牛顿溃坝和干颗粒柱崩塌进行独立验证。然后通过两种流固相互作用情况来评估耦合框架：立方体进入非牛顿流体和涉及非牛顿流体和非球形颗粒混合物的复杂溃坝情景。模拟的流动行为（如立方体穿透深度和流固混合物的下游传播）与相应的实验结果吻合良好，验证了所提出的边界处理方法的有效性和框架模拟非牛顿流体和非球形颗粒流固相互作用系统的准确性。

{"title":"Modeling non-Newtonian fluid–solid flows containing non-spherical particles by the SPH-DEM coupling model","authors":"Jiahang Du , Yihang Yu , Yuqi Guo , Huaqing Ma , Yongzhi Zhao","doi":"10.1016/j.compgeo.2026.107990","DOIUrl":"10.1016/j.compgeo.2026.107990","url":null,"abstract":"<div><div>Given the ubiquity and significance of fluid–solid systems composed of non-Newtonian fluids and non-spherical particles in nature and industry, the deep investigation of their interaction dynamics is crucial for geological studies and engineering applications. In this paper, a resolved coupling model based on Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) is proposed to simulate the interaction between non-Newtonian fluids and non-spherical particles. The non-Newtonian fluid is modelled using SPH with the power-law model, while the dynamics of the non-spherical particles is simulated using DEM incorporated with the super-ellipsoid and polyhedral particle models. A modified boundary repulsive force model, combined with the fixed ghost particle method, is developed to handle boundary conditions. For super-ellipsoid and polyhedron boundaries, an improved approach is proposed to generate uniformly distributed, body-fitted boundary dummy particles on their surfaces, thereby ensuring computational accuracy near boundaries and in fluid–solid interactions. To rigorously validate the framework, four benchmark cases are conducted, and the resulting numerical simulations are subsequently compared with experiments in this study. The SPH and DEM methods are first independently validated using a 3D non-Newtonian dam-break and a dry granular column collapse, respectively. The coupled framework is then assessed through two fluid–solid interaction cases: cube entry into a non-Newtonian fluid and complex dam-break scenarios involving mixtures of non-Newtonian fluids and non-spherical particles. The simulated flow behaviors (<em>e</em>.<em>g</em>., the cube penetration depth and downstream propagation of fluid–solid mixture) agree well with the corresponding experiments, validating the effectiveness of the proposed boundary treatment method and the accuracy of the framework to simulate the fluid–solid interaction systems involving non-Newtonian fluids and non-spherical particles.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107990"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146174285","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}

引用次数: 0

Simulation of brittle particle breakage using a spherical harmonic-based discrete element method 基于球面谐波的离散元法模拟脆性颗粒破碎

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-11 DOI: 10.1016/j.compgeo.2026.107972

Jiabao Gao , Daosheng Ling , Fubin Tu

Particle breakage significantly influences the strength and deformation behavior of materials in geotechnical engineering and powder processing. However, conventional discrete element methods remain limited in capturing irregular particle morphology and complex breakage mechanisms. This study develops a 3D breakage model grounded in the spherical harmonic-based discrete element method (SH-DEM). A dual-center definition is introduced to couple geometric representation with dynamic response, forming a breakage identification method applicable to 3D irregular particles. Fracture planes are identified using the Boussinesq–Cerruti analytical solution together with Dijkstra path search, with predictions validated against finite element analyses and experimental data. Numerical single-particle compression simulations reveal an inverse correlation between fragment number and the major-fragment volume fraction. They further demonstrate strong morphology-dependent breakage. Spherical and plate-like particles exhibit multifragment breakage with abundant fines, while blade-like and rod-like particles fail along high-curvature regions, resulting in through-particle separation. These findings underscore the central influence of particle morphology on brittle breakage mechanisms.

在岩土工程和粉体加工中，颗粒破碎对材料的强度和变形行为有重要影响。然而，传统的离散元方法在捕获不规则颗粒形态和复杂的断裂机制方面仍然有限。本文建立了基于球面谐波离散元法（SH-DEM）的三维断裂模型。引入双中心定义，将几何表示与动态响应耦合，形成了一种适用于三维不规则颗粒的破碎识别方法。利用Boussinesq-Cerruti解析解和Dijkstra路径搜索识别裂缝面，并通过有限元分析和实验数据验证预测结果。单颗粒压缩数值模拟表明，碎片数与主碎片体积分数呈负相关关系。它们进一步证明了强烈的形态依赖性断裂。球形和片状颗粒表现为多碎片破碎，细粒丰富，而叶片状和棒状颗粒沿高曲率区域破坏，导致颗粒间分离。这些发现强调了颗粒形态对脆性断裂机制的核心影响。

{"title":"Simulation of brittle particle breakage using a spherical harmonic-based discrete element method","authors":"Jiabao Gao , Daosheng Ling , Fubin Tu","doi":"10.1016/j.compgeo.2026.107972","DOIUrl":"10.1016/j.compgeo.2026.107972","url":null,"abstract":"<div><div>Particle breakage significantly influences the strength and deformation behavior of materials in geotechnical engineering and powder processing. However, conventional discrete element methods remain limited in capturing irregular particle morphology and complex breakage mechanisms. This study develops a 3D breakage model grounded in the spherical harmonic-based discrete element method (SH-DEM). A dual-center definition is introduced to couple geometric representation with dynamic response, forming a breakage identification method applicable to 3D irregular particles. Fracture planes are identified using the Boussinesq–Cerruti analytical solution together with Dijkstra path search, with predictions validated against finite element analyses and experimental data. Numerical single-particle compression simulations reveal an inverse correlation between fragment number and the major-fragment volume fraction. They further demonstrate strong morphology-dependent breakage. Spherical and plate-like particles exhibit multifragment breakage with abundant fines, while blade-like and rod-like particles fail along high-curvature regions, resulting in through-particle separation. These findings underscore the central influence of particle morphology on brittle breakage mechanisms.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107972"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146174284","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}

引用次数: 0

A phase field fracture model for rock creep: Theoretical framework and engineering applications 岩石蠕变相场断裂模型：理论框架与工程应用

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-18 DOI: 10.1016/j.compgeo.2026.108007

Rui Liu , Shuwei Zhou , Shikang Qin , Chengkai Zhang , Meng Xing

This study proposes a double phase field framework for modeling rock creep fracture, which incorporates the tensile and compressive strain induced fracture and considers the coupling between rock creep and fracture. This framework enables the separation of tensile and compressive–shear cracks through a double phase field model, which subsequently influences the rock’s constitutive matrix. Furthermore, it accounts for the mutual coupling between creep and fracture mechanisms. Experimental observations confirm that the proposed model aligns with existing tests in three aspects: creep behavior, failure patterns, and creep-induced failure. Studies on specimens containing single and double prefabricated flaws reveal that: creep induces the development of compressive cracks; under identical loading pressure, creep rates do not alter crack propagation directions; the proposed phase-field fracture model effectively captures crack evolution during rock fracturing. Further validation against engineering practices indicates that the model demonstrates promising performance in characterizing surrounding rock damage in underground engineering.

本文提出了一种模拟岩石蠕变断裂的双相场框架，该框架考虑了岩石蠕变与断裂之间的耦合，同时考虑了拉应变和压应变引起的断裂。该框架通过双相场模型实现了拉伸和压缩-剪切裂缝的分离，从而影响了岩石的本构矩阵。此外，它还解释了蠕变与断裂机制之间的相互耦合。实验观察证实，所提出的模型在三个方面与现有的试验一致：蠕变行为、破坏模式和蠕变诱发的破坏。对含单、双预制裂纹试件的研究表明：蠕变诱发压缩裂纹的发展；在相同加载压力下，蠕变速率不改变裂纹扩展方向；所提出的相场裂缝模型有效地捕捉了岩石破裂过程中的裂缝演化过程。进一步的工程实践验证表明，该模型对地下工程围岩损伤的表征具有良好的效果。

{"title":"A phase field fracture model for rock creep: Theoretical framework and engineering applications","authors":"Rui Liu , Shuwei Zhou , Shikang Qin , Chengkai Zhang , Meng Xing","doi":"10.1016/j.compgeo.2026.108007","DOIUrl":"10.1016/j.compgeo.2026.108007","url":null,"abstract":"<div><div>This study proposes a double phase field framework for modeling rock creep fracture, which incorporates the tensile and compressive strain induced fracture and considers the coupling between rock creep and fracture. This framework enables the separation of tensile and compressive–shear cracks through a double phase field model, which subsequently influences the rock’s constitutive matrix. Furthermore, it accounts for the mutual coupling between creep and fracture mechanisms. Experimental observations confirm that the proposed model aligns with existing tests in three aspects: creep behavior, failure patterns, and creep-induced failure. Studies on specimens containing single and double prefabricated flaws reveal that: creep induces the development of compressive cracks; under identical loading pressure, creep rates do not alter crack propagation directions; the proposed phase-field fracture model effectively captures crack evolution during rock fracturing. Further validation against engineering practices indicates that the model demonstrates promising performance in characterizing surrounding rock damage in underground engineering.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 108007"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147385400","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}

引用次数: 0

A three-dimensional arbitrary grid material point method for large deformation problems with geometrically complex boundaries 具有几何复杂边界的大变形问题的三维任意网格材料点法

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-14 DOI: 10.1016/j.compgeo.2026.107993

Hongyu Ma, Ruopu Zhou, Xiong Zhang

The Material Point Method (MPM) is a hybrid “mesh-particle” approach for large deformation simulations that avoids the high computational costs and boundary challenges of traditional meshfree particle methods. However, bottlenecks arise when modeling geometrically complex boundaries, especially in three dimensions. In this article, a three-dimensional arbitrary grid material point method (3D-AGMPM) is proposed to construct complex boundaries using polyhedral grid cells. This novel method establishes a unified framework for various cell types. In the particle-to-grid mapping procedure, the 3D-Wachspress basis functions are introduced as shape functions for arbitrary grid cells, while an improved hash-cell based particle localization algorithm is proposed to enhance computational efficiency. Nonlinear frictional boundary conditions are proposed in a trial-correction framework to model the frictional interaction between the body and complex geometric surfaces. In addition, a contact algorithm for polyhedral cells is proposed to handle contacts among multiple bodies within the arbitrary grid. Several benchmarks demonstrate the effectiveness of the proposed 3D-AGMPM for handling problems involving complex geometric boundaries. Furthermore, practical case studies, such as the Wangjiayan landslide, highlight its robustness and potential in addressing large-scale, complex engineering problems. Due to the simple implementation, flexibility, and high efficiency, the proposed 3D-AGMPM shows promise as a powerful tool for solving large deformation problems with geometrically complex boundaries.

材料点法（MPM）是一种用于大变形模拟的混合“网格-粒子”方法，它避免了传统的无网格粒子方法的高计算成本和边界挑战。然而，在几何上复杂的边界建模时，特别是在三维空间中，就会出现瓶颈。本文提出了一种利用多面体网格单元构建复杂边界的三维任意网格材料点法（3D-AGMPM）。这种新方法为各种细胞类型建立了统一的框架。在粒子到网格的映射过程中，引入3D-Wachspress基函数作为任意网格单元的形状函数，提出了一种改进的基于哈希单元的粒子定位算法，提高了计算效率。在试验-修正框架中提出了非线性摩擦边界条件来模拟物体与复杂几何表面之间的摩擦相互作用。此外，提出了一种多面体单元接触算法，用于处理任意网格内多个体之间的接触。几个基准测试证明了所提出的3D-AGMPM在处理涉及复杂几何边界的问题方面的有效性。此外，实际案例研究，如王家岩滑坡，突出了其在解决大规模复杂工程问题方面的稳健性和潜力。由于实现简单、灵活和高效，所提出的3D-AGMPM有望成为解决具有几何复杂边界的大变形问题的有力工具。

{"title":"A three-dimensional arbitrary grid material point method for large deformation problems with geometrically complex boundaries","authors":"Hongyu Ma, Ruopu Zhou, Xiong Zhang","doi":"10.1016/j.compgeo.2026.107993","DOIUrl":"10.1016/j.compgeo.2026.107993","url":null,"abstract":"<div><div>The Material Point Method (MPM) is a hybrid “mesh-particle” approach for large deformation simulations that avoids the high computational costs and boundary challenges of traditional meshfree particle methods. However, bottlenecks arise when modeling geometrically complex boundaries, especially in three dimensions. In this article, a three-dimensional arbitrary grid material point method (3D-AGMPM) is proposed to construct complex boundaries using polyhedral grid cells. This novel method establishes a unified framework for various cell types. In the particle-to-grid mapping procedure, the 3D-Wachspress basis functions are introduced as shape functions for arbitrary grid cells, while an improved hash-cell based particle localization algorithm is proposed to enhance computational efficiency. Nonlinear frictional boundary conditions are proposed in a trial-correction framework to model the frictional interaction between the body and complex geometric surfaces. In addition, a contact algorithm for polyhedral cells is proposed to handle contacts among multiple bodies within the arbitrary grid. Several benchmarks demonstrate the effectiveness of the proposed 3D-AGMPM for handling problems involving complex geometric boundaries. Furthermore, practical case studies, such as the Wangjiayan landslide, highlight its robustness and potential in addressing large-scale, complex engineering problems. Due to the simple implementation, flexibility, and high efficiency, the proposed 3D-AGMPM shows promise as a powerful tool for solving large deformation problems with geometrically complex boundaries.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107993"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147385403","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}

引用次数: 0

Interface behavior of concrete cut-off walls: a thermodynamics-based generalized elasto-plastic-damage model 混凝土防渗墙界面行为：基于热力学的广义弹塑性损伤模型

IF 6.2 1区工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS

Computers and Geotechnics

Pub Date : 2026-05-01 Epub Date: 2026-02-20 DOI: 10.1016/j.compgeo.2026.107966

Gang Wang , Mingke Liao , Zhichao Zhang , Wei Jin

A generalized elasto-plastic-damage model unified for stiff and flexible interface behaviors in concrete cut-off walls is developed in this paper based on interface thermodynamics. The state-dependent couplings among elasticity, plasticity, and damage in coarse interface behavior are accounted for within the thermodynamic framework, with considerations of the interface roughness. The concept of interface contact fabric is proposed to take into account the effect of contact normal orientations at meso scale, which evolves gradually under shearing to allow the mobilization of higher shear strengths. It is also coupled with the interface damage that induces degradations of both the contact fabric and the interface roughness. This leads to significant strain softening even in contractive stiff interfaces. The proposed model is validated by predicting a series of interface shear tests implemented for three groups of concrete and concrete-mud interfaces. It is shown that the model well captures the effects of interface roughness on the diverse behaviors of shear hardening/softening and shear dilation/contraction in both stiff and flexible interfaces, with the interface damage and the contact fabric evolution as two basic underlying mechanisms.

基于界面热力学，建立了混凝土防渗墙刚柔界面行为统一的广义弹塑性损伤模型。在热力学框架内，考虑了界面粗糙度，考虑了粗界面行为中弹性、塑性和损伤之间的状态依赖耦合。提出了界面接触织物的概念，考虑了接触法向在细观尺度上的影响，这种影响在剪切作用下逐渐演变，从而可以调动更高的剪切强度。它还与界面损伤相结合，导致接触织物和界面粗糙度的退化。这导致显著的应变软化，即使在收缩刚性界面。通过预测三组混凝土和混凝土-泥浆界面的一系列界面剪切试验，验证了该模型的有效性。结果表明，该模型较好地反映了界面粗糙度对刚性和柔性界面剪切硬化/软化和剪切扩张/收缩等多种行为的影响，界面损伤和接触组织演化是两种基本的潜在机制。

{"title":"Interface behavior of concrete cut-off walls: a thermodynamics-based generalized elasto-plastic-damage model","authors":"Gang Wang , Mingke Liao , Zhichao Zhang , Wei Jin","doi":"10.1016/j.compgeo.2026.107966","DOIUrl":"10.1016/j.compgeo.2026.107966","url":null,"abstract":"<div><div>A generalized elasto-plastic-damage model unified for stiff and flexible interface behaviors in concrete cut-off walls is developed in this paper based on interface thermodynamics. The state-dependent couplings among elasticity, plasticity, and damage in coarse interface behavior are accounted for within the thermodynamic framework, with considerations of the interface roughness. The concept of interface contact fabric is proposed to take into account the effect of contact normal orientations at meso scale, which evolves gradually under shearing to allow the mobilization of higher shear strengths. It is also coupled with the interface damage that induces degradations of both the contact fabric and the interface roughness. This leads to significant strain softening even in contractive stiff interfaces. The proposed model is validated by predicting a series of interface shear tests implemented for three groups of concrete and concrete-mud interfaces. It is shown that the model well captures the effects of interface roughness on the diverse behaviors of shear hardening/softening and shear dilation/contraction in both stiff and flexible interfaces, with the interface damage and the contact fabric evolution as two basic underlying mechanisms.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107966"},"PeriodicalIF":6.2,"publicationDate":"2026-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147385469","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}

引用次数: 0