首页 > 最新文献

Computers and Geotechnics最新文献

英文 中文
Micromechanical mechanisms of damage evolution and failure transition in sandstone: insights from physically informed digital rock models 砂岩损伤演化和破坏转变的微观力学机制:来自物理信息的数字岩石模型的见解
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-01-02 DOI: 10.1016/j.compgeo.2025.107892
Jiaxin Feng , Xu Yang , Gao Li , Yi Zhang , Hongtao Li , Mubai Duan
Sandstone is widespread in the upper crust and underpins many subsurface engineering and energy applications. Understanding its deformation and failure under stress conditions therefore remains essential. However, the micromechanical mechanisms governing damage evolution and failure transitions under confining pressure remain poorly understood, largely owing to limited experimental reproducibility and intrinsic sample heterogeneity. Discrete element method (DEM)-based digital rock simulation provides a promising way to interrogate micromechanical damage processes, yet contact-parameter calibration often relies on trial-and-error, yielding non-unique, scenario-specific parameter sets and limiting reproducibility and transferability. This study reconstructs a DEM-based digital rock using the mineralogical and microstructural features of sandstone from the Yanchang Formation in the Ordos Basin. Mineral-scale elastic modulus and fracture toughness were determined through nanoindentation and combined with semi-circular bend simulations of individual mineral phases to calibrate the micromechanical parameters of the digital rock. The calibrated parameters were validated against the Hoek-Brown strength criterion and mechanical responses under varying loading conditions, confirming their robustness and applicability. Using the physically informed calibrated parameter set, we performed pseudo-triaxial compression simulations of sandstone under different confining pressures. Spatiotemporal tracking of microcracks shows that higher confining pressure postpones intergranular tensile cracking and amplifies strain incompatibility at interfaces between soft and stiff minerals, which triggers intragranular tensile cracking inside stiff grains. Force chain analysis indicates a confinement-driven alignment and densification of load paths that homogenizes stress transfer and delays instability. Joint microcrack and acoustic emission statistics document a transition from tensile-dominated failure at low confinement to shear-dominated failure at higher confinement. The physically informed calibrated parameters provide a unified framework for digital rock simulations under similar conditions, demonstrating their potential for capturing the mechanisms of damage evolution and failure in rock under confining pressure constraints.
砂岩广泛存在于地壳上部,是许多地下工程和能源应用的基础。因此,了解其在应力条件下的变形和破坏仍然是必要的。然而,控制围压下损伤演化和破坏转变的微观力学机制仍然知之甚少,这主要是由于有限的实验可重复性和固有的样品异质性。基于离散元方法(DEM)的数字岩石模拟为研究微力学损伤过程提供了一种很有前途的方法,但接触参数校准通常依赖于试错,产生的参数集不独特,且场景特定,限制了可重复性和可转移性。利用鄂尔多斯盆地延长组砂岩的矿物学和微观结构特征,重建了基于dem的数字岩石。通过纳米压痕测定矿物尺度的弹性模量和断裂韧性,并结合单个矿物相的半圆弯曲模拟,标定数字岩石的微观力学参数。根据Hoek-Brown强度准则和不同载荷条件下的力学响应对校准参数进行了验证,验证了其鲁棒性和适用性。利用物理信息校准参数集,我们对不同围压下的砂岩进行了伪三轴压缩模拟。对微裂纹的时空追踪表明,较高的围压延缓了晶间拉伸开裂,放大了软、硬矿物界面处的应变不相容,从而引发了硬颗粒内部的晶内拉伸开裂。力链分析表明,约束驱动的载荷路径对齐和致密化,均匀化应力传递和延迟不稳定性。节理微裂纹和声发射统计数据表明,节理从低约束条件下以拉伸为主的破坏转变为高约束条件下以剪切为主的破坏。物理信息校准参数为类似条件下的数字岩石模拟提供了统一的框架,展示了它们在捕获围压约束下岩石损伤演化和破坏机制方面的潜力。
{"title":"Micromechanical mechanisms of damage evolution and failure transition in sandstone: insights from physically informed digital rock models","authors":"Jiaxin Feng ,&nbsp;Xu Yang ,&nbsp;Gao Li ,&nbsp;Yi Zhang ,&nbsp;Hongtao Li ,&nbsp;Mubai Duan","doi":"10.1016/j.compgeo.2025.107892","DOIUrl":"10.1016/j.compgeo.2025.107892","url":null,"abstract":"<div><div>Sandstone is widespread in the upper crust and underpins many subsurface engineering and energy applications. Understanding its deformation and failure under stress conditions therefore remains essential. However, the micromechanical mechanisms governing damage evolution and failure transitions under confining pressure remain poorly understood, largely owing to limited experimental reproducibility and intrinsic sample heterogeneity. Discrete element method (DEM)-based digital rock simulation provides a promising way to interrogate micromechanical damage processes, yet contact-parameter calibration often relies on trial-and-error, yielding non-unique, scenario-specific parameter sets and limiting reproducibility and transferability. This study reconstructs a DEM-based digital rock using the mineralogical and microstructural features of sandstone from the Yanchang Formation in the Ordos Basin. Mineral-scale elastic modulus and fracture toughness were determined through nanoindentation and combined with semi-circular bend simulations of individual mineral phases to calibrate the micromechanical parameters of the digital rock. The calibrated parameters were validated against the Hoek-Brown strength criterion and mechanical responses under varying loading conditions, confirming their robustness and applicability. Using the physically informed calibrated parameter set, we performed pseudo-triaxial compression simulations of sandstone under different confining pressures. Spatiotemporal tracking of microcracks shows that higher confining pressure postpones intergranular tensile cracking and amplifies strain incompatibility at interfaces between soft and stiff minerals, which triggers intragranular tensile cracking inside stiff grains. Force chain analysis indicates a confinement-driven alignment and densification of load paths that homogenizes stress transfer and delays instability. Joint microcrack and acoustic emission statistics document a transition from tensile-dominated failure at low confinement to shear-dominated failure at higher confinement. The physically informed calibrated parameters provide a unified framework for digital rock simulations under similar conditions, demonstrating their potential for capturing the mechanisms of damage evolution and failure in rock under confining pressure constraints.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107892"},"PeriodicalIF":6.2,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883614","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
An improved bond-associated peridynamics for modelling elastoplastic deformation and progressive fracture in geomaterials 用于模拟岩土材料弹塑性变形和渐进断裂的改进粘结相关周动力学
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-01-02 DOI: 10.1016/j.compgeo.2025.107868
Yixin Li , Zhen-Yu Yin , Luyu Wang , Xueyu Geng
Geomaterials are pressure-dependent and exhibit non-linear elastoplastic deformation and failure behaviours, including strain localisation and complex cracking. To study these phenomena, this paper develops an improved non-ordinary state-based peridynamic (NOSB-PD) framework for simulating the elastoplastic response and progressive failure of geomaterials. A bond-associated (BA) stabilisation approach is incorporated to eliminate zero-energy deformation modes, while a divergence-form peridynamic differential operator (PDDO) is employed to remove residual stresses near boundaries. Using the principle of virtual work, a new equation of motion is derived that naturally incorporates both zero and non-zero traction boundary conditions. Furthermore, a modified staggered algorithm is introduced to mitigate spurious plastic yielding arising from surface effects and inaccuracies in stress history. The improved BA-NOSB-PD framework enables direct implementation of elastoplastic constitutive laws and mixed-mode failure criteria. Model accuracy and convergence performance are verified against analytical solutions. The numerical simulations confirm the model’s capability to capture elastoplastic behaviour, including stress redistribution, shear band development, and plasticity-induced crack propagation. Moreover, the mixed-mode fracture mechanisms arising from different pre-existing fissure configurations are thoroughly analysed. The results demonstrate that the proposed model accurately captures shear band formation and fracture coalescence, exhibiting good numerical stability and convergence. The improved PD model provides a robust and feasible framework for investigating elastoplastic behaviour and the underlying failure mechanisms in geomaterials.
岩土材料具有压力依赖性,并表现出非线性弹塑性变形和破坏行为,包括应变局部化和复杂开裂。为了研究这些现象,本文开发了一种改进的非普通基于状态的周动力学(NOSB-PD)框架,用于模拟岩土材料的弹塑性响应和渐进破坏。采用键相关(BA)稳定方法消除零能量变形模式,采用发散形式的周动力学微分算子(PDDO)消除边界附近的残余应力。利用虚功原理,推导了一个包含零和非零牵引力边界条件的新运动方程。此外,还引入了一种改进的交错算法,以减轻由于表面效应和应力历史不准确引起的伪塑性屈服。改进的BA-NOSB-PD框架可以直接实现弹塑性本构律和混合模式破坏准则。通过解析解验证了模型的精度和收敛性。数值模拟证实了该模型能够捕捉弹塑性行为,包括应力重分布、剪切带发展和塑性诱导裂纹扩展。此外,还深入分析了不同既有裂缝形态导致的混合模式断裂机制。结果表明,该模型能较准确地反映剪切带的形成和裂缝的聚并,具有较好的数值稳定性和收敛性。改进的PD模型为研究岩土材料的弹塑性行为和潜在的破坏机制提供了一个鲁棒和可行的框架。
{"title":"An improved bond-associated peridynamics for modelling elastoplastic deformation and progressive fracture in geomaterials","authors":"Yixin Li ,&nbsp;Zhen-Yu Yin ,&nbsp;Luyu Wang ,&nbsp;Xueyu Geng","doi":"10.1016/j.compgeo.2025.107868","DOIUrl":"10.1016/j.compgeo.2025.107868","url":null,"abstract":"<div><div>Geomaterials are pressure-dependent and exhibit non-linear elastoplastic deformation and failure behaviours, including strain localisation and complex cracking. To study these phenomena, this paper develops an improved non-ordinary state-based peridynamic (NOSB-PD) framework for simulating the elastoplastic response and progressive failure of geomaterials. A bond-associated (BA) stabilisation approach is incorporated to eliminate zero-energy deformation modes, while a divergence-form peridynamic differential operator (PDDO) is employed to remove residual stresses near boundaries. Using the principle of virtual work, a new equation of motion is derived that naturally incorporates both zero and non-zero traction boundary conditions. Furthermore, a modified staggered algorithm is introduced to mitigate spurious plastic yielding arising from surface effects and inaccuracies in stress history. The improved BA-NOSB-PD framework enables direct implementation of elastoplastic constitutive laws and mixed-mode failure criteria. Model accuracy and convergence performance are verified against analytical solutions. The numerical simulations confirm the model’s capability to capture elastoplastic behaviour, including stress redistribution, shear band development, and plasticity-induced crack propagation. Moreover, the mixed-mode fracture mechanisms arising from different pre-existing fissure configurations are thoroughly analysed. The results demonstrate that the proposed model accurately captures shear band formation and fracture coalescence, exhibiting good numerical stability and convergence. The improved PD model provides a robust and feasible framework for investigating elastoplastic behaviour and the underlying failure mechanisms in geomaterials.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107868"},"PeriodicalIF":6.2,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883638","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 coupled strain-damage framework for rock deformation and dynamic damage under cyclic loading based on mesoscale elastoplastic elements 基于细观弹塑性单元的循环加载下岩石变形与动力损伤耦合应变-损伤框架
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2026-01-02 DOI: 10.1016/j.compgeo.2025.107879
Haohao Zhang , Rui Chen , Yazhou Wang , Huijie Bi , Chaojun Fan , Hui Peng
In this study, a coupled strain-damage framework based on mesoscale elastoplastic elements is proposed to numerically investigate rock deformation and dynamic damage (from the initiation of mesoscale damage to macroscopic failure) under cyclic loading. Building on mesoscale elastic-damage mechanics, elements are upgraded to elastoplastic counterparts to capture plastic strain accumulation and strength degradation during deformation. Damage is formulated as a strain- and time-dependent field variable and solved together with the solid-mechanics field through a staggered coupling scheme. The framework is evaluated by cyclic uniaxial compression and cyclic Brazilian splitting tests. The simulations reproduce the complete deformation–failure process and key laboratory-observed features, including hysteresis, irreversible strain growth, and fatigue failure. Quantitative validation shows that key mechanical parameters are reproduced with relative errors below 5% (elastic modulus: 4.6%, UCS: 0.9%, UTS: 0.4%, and failure angle: 3.2%). The results indicate that the ultimate failure strain under cyclic loading is governed by its monotonic counterpart, and they quantify the influences of specimen heterogeneity (homogeneity index m), maximum stress ratio, and loading frequency on fatigue life. The proposed approach provides a physics-grounded and computationally tractable basis for analysing rock fatigue and failure in geoengineering settings subjected to dynamic loading.
本文提出了一种基于细观弹塑性单元的应变-损伤耦合框架,对循环加载下岩石的变形和动态损伤(从细观损伤开始到宏观破坏)进行了数值研究。在中尺度弹性损伤力学的基础上,单元升级为弹塑性对应体,以捕获变形过程中的塑性应变积累和强度退化。将损伤表述为与应变和时间相关的场变量,并通过交错耦合格式与固体力学场一起求解。通过循环单轴压缩和循环巴西劈裂试验对框架进行了评价。模拟再现了完整的变形破坏过程和实验室观察到的关键特征,包括迟滞、不可逆应变增长和疲劳破坏。定量验证表明,关键力学参数的再现相对误差小于5%(弹性模量4.6%,UCS 0.9%, UTS 0.4%,失效角3.2%)。结果表明,循环加载下的极限破坏应变受单调应变控制,并量化了试样非均质性(均匀性指数m)、最大应力比和加载频率对疲劳寿命的影响。所提出的方法为分析地球工程环境下受动载荷影响的岩石疲劳和破坏提供了物理基础和易于计算的基础。
{"title":"A coupled strain-damage framework for rock deformation and dynamic damage under cyclic loading based on mesoscale elastoplastic elements","authors":"Haohao Zhang ,&nbsp;Rui Chen ,&nbsp;Yazhou Wang ,&nbsp;Huijie Bi ,&nbsp;Chaojun Fan ,&nbsp;Hui Peng","doi":"10.1016/j.compgeo.2025.107879","DOIUrl":"10.1016/j.compgeo.2025.107879","url":null,"abstract":"<div><div>In this study, a coupled strain-damage framework based on mesoscale elastoplastic elements is proposed to numerically investigate rock deformation and dynamic damage (from the initiation of mesoscale damage to macroscopic failure) under cyclic loading. Building on mesoscale elastic-damage mechanics, elements are upgraded to elastoplastic counterparts to capture plastic strain accumulation and strength degradation during deformation. Damage is formulated as a strain- and time-dependent field variable and solved together with the solid-mechanics field through a staggered coupling scheme. The framework is evaluated by cyclic uniaxial compression and cyclic Brazilian splitting tests. The simulations reproduce the complete deformation–failure process and key laboratory-observed features, including hysteresis, irreversible strain growth, and fatigue failure. Quantitative validation shows that key mechanical parameters are reproduced with relative errors below 5% (elastic modulus: 4.6%, UCS: 0.9%, UTS: 0.4%, and failure angle: 3.2%). The results indicate that the ultimate failure strain under cyclic loading is governed by its monotonic counterpart, and they quantify the influences of specimen heterogeneity (homogeneity index <em>m</em>), maximum stress ratio, and loading frequency on fatigue life. The proposed approach provides a physics-grounded and computationally tractable basis for analysing rock fatigue and failure in geoengineering settings subjected to dynamic loading.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107879"},"PeriodicalIF":6.2,"publicationDate":"2026-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883616","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
Adaptive FEM-SPH for anchored slope failure with fractal dimension crack analysis 基于分维裂纹分析的自适应有限元- sph法锚固边坡破坏
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107872
Zhengnan Tu , Qiang Xie , Jun Feng , Hexing Zhang , Yucheng Chen
Anchored slopes are widely used in geotechnical engineering, yet their failure mechanisms under ultimate loading-particularly crack evolution and slip-surface formation-remain insufficiently understood. This study develops an adaptive FEM-SPH coupling framework governed by a maximum-principal-strain criterion and validates it against a laboratory pull-out test on a single BFRP anchor. The method automatically converts excessively distorted finite elements into SPH particles while preserving stress and damage continuity, thus combining FEM accuracy in the small-strain regime with SPH robustness at large deformations. Validation confirms its capability to reproduce the measured load–displacement response and shear-damage evolution with high fidelity. A prototype slope from the Yunnan high-speed railway project was simulated under a unified crest-displacement loading of 500 mm for four anchoring layouts. The cross-arranged scheme achieves the highest peak contact stress and residual capacity due to multidirectional constraint, while gently inclined anchors mobilize higher axial forces and steeply inclined anchors exhibit faster crack-complexity growth. The box-counting fractal dimension D effectively quantifies the spatial evolution of crack networks, providing a geometric measure of damage complexity and its progression under large deformation, to be interpreted together with conventional mechanical responses. The D-S evolution exhibits a three-stage, nonlinear development of crack-network complexity and reveals layout-dependent differences in the growth rate and saturation of cracking, with the cross-arranged layout, within the present parameter space, combining the slowest complexity growth and a smooth approach to saturation with the highest global bearing capacity. From an engineering standpoint, the findings provide comparative references for representative anchorage layouts in large-deformation marly slopes.
锚固边坡在岩土工程中有着广泛的应用,但其在极限荷载作用下的破坏机制,特别是裂缝演化和滑面形成机制,目前还没有得到充分的了解。本研究开发了一个由最大主应变准则控制的自适应FEM-SPH耦合框架,并对单个BFRP锚进行了实验室拉拔试验。该方法在保持应力和损伤连续性的同时,将过度变形的有限元自动转换为SPH粒子,从而将小应变状态下的FEM精度与大变形状态下的SPH鲁棒性相结合。验证证实其能够高保真地再现所测载荷-位移响应和剪切-损伤演变。以云南高铁工程原型边坡为例,对4种锚固布置方式在500 mm统一顶移荷载下进行了模拟。由于多向约束,交叉布置方案的峰值接触应力和残余承载力最高,缓倾斜锚杆动员的轴向力更高,急倾斜锚杆的裂纹复杂性增长更快。盒计数分形维数D有效地量化了裂纹网络的空间演化,提供了大变形下损伤复杂性及其进展的几何度量,可以与常规力学响应一起解释。D-S演化表现为裂缝网络复杂性的三阶段非线性发展,并揭示了裂缝扩展速率和饱和度的分布差异,在当前参数空间内,以交叉布置的布局,结合了最慢的复杂性增长和最高的全局承载能力的平滑饱和。从工程角度看,研究结果可为具有代表性的大变形泥质边坡锚固布置提供比较参考。
{"title":"Adaptive FEM-SPH for anchored slope failure with fractal dimension crack analysis","authors":"Zhengnan Tu ,&nbsp;Qiang Xie ,&nbsp;Jun Feng ,&nbsp;Hexing Zhang ,&nbsp;Yucheng Chen","doi":"10.1016/j.compgeo.2025.107872","DOIUrl":"10.1016/j.compgeo.2025.107872","url":null,"abstract":"<div><div>Anchored slopes are widely used in geotechnical engineering, yet their failure mechanisms under ultimate loading-particularly crack evolution and slip-surface formation-remain insufficiently understood. This study develops an adaptive FEM-SPH coupling framework governed by a maximum-principal-strain criterion and validates it against a laboratory pull-out test on a single BFRP anchor. The method automatically converts excessively distorted finite elements into SPH particles while preserving stress and damage continuity, thus combining FEM accuracy in the small-strain regime with SPH robustness at large deformations. Validation confirms its capability to reproduce the measured load–displacement response and shear-damage evolution with high fidelity. A prototype slope from the Yunnan high-speed railway project was simulated under a unified crest-displacement loading of 500 mm for four anchoring layouts. The cross-arranged scheme achieves the highest peak contact stress and residual capacity due to multidirectional constraint, while gently inclined anchors mobilize higher axial forces and steeply inclined anchors exhibit faster crack-complexity growth. The box-counting fractal dimension <em>D</em> effectively quantifies the spatial evolution of crack networks, providing a geometric measure of damage complexity and its progression under large deformation, to be interpreted together with conventional mechanical responses. The <em>D-S</em> evolution exhibits a three-stage, nonlinear development of crack-network complexity and reveals layout-dependent differences in the growth rate and saturation of cracking, with the cross-arranged layout, within the present parameter space, combining the slowest complexity growth and a smooth approach to saturation with the highest global bearing capacity. From an engineering standpoint, the findings provide comparative references for representative anchorage layouts in large-deformation marly slopes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107872"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883708","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
Three-dimensional fracture modelling of heterogeneous and anisotropic rock material under true triaxial compression 真三轴压缩条件下非均质和各向异性岩石材料三维裂缝建模
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107862
Hui Li, Shanyong Wang
Underground excavation inevitably disturbs natural in-situ stress equilibrium of rocks, leading to localized deformation, fracture initiation, and potential large-scale failure. Understanding the failure mechanisms of rock under complex stress states is of great importance for ensuring the safety and stability of underground structures. To this end, a mesoscale numerical approach combining Weibull random fields with a hybrid phase-field model is adopted to investigate the failure mechanisms of heterogeneous rock materials under true triaxial compression. The isotropic model is further extended to anisotropic formulation by incorporating an anisotropic crack surface density function to account for directional fracture behaviour. First, the isotropic model is calibrated against experimental data by simulating triaxial compression of intact cuboid sandstone specimens, with the effects of confining stress and mesoscale heterogeneity systematically analysed. Next, the influence of pre-existing notches on 3D crack propagation is studied using the calibrated model. Finally, the anisotropic model is used to simulate the failure of intact and notched laminated rocks. It is found that rock fracture behaviour is governed by the combined effects of confining stress, mesoscale heterogeneity, and notch geometry. The incorporation of anisotropic fracture toughness enables directional control of crack propagation along bedding planes, successfully capturing joint-induced compression–shear failure. The developed approach provides a robust tool for simulating 3D failure in heterogeneous and anisotropic rocks under complex loading conditions.
地下开挖不可避免地破坏岩石的天然地应力平衡,导致局部变形、起裂,并可能发生大规模破坏。了解岩石在复杂应力状态下的破坏机制,对保证地下结构的安全稳定具有重要意义。为此,采用威布尔随机场与混合相场模型相结合的中尺度数值方法研究非均质岩石材料在真三轴压缩下的破坏机制。通过加入各向异性裂缝表面密度函数来解释定向断裂行为,各向同性模型进一步扩展为各向异性公式。首先,通过模拟完整长方体砂岩试件的三轴压缩,根据实验数据对各向同性模型进行校正,系统分析围应力和中尺度非均质性对模型的影响;其次,利用校正后的模型研究了预先存在的缺口对三维裂纹扩展的影响。最后,采用各向异性模型对完整和缺口层状岩石的破坏进行了模拟。研究发现,岩石断裂行为受围应力、中尺度非均质性和缺口几何形状的综合影响。各向异性断裂韧性的结合使裂缝沿顺层面扩展的定向控制成为可能,成功捕获节理引起的压剪破坏。所开发的方法为模拟复杂加载条件下非均质和各向异性岩石的三维破坏提供了一个强大的工具。
{"title":"Three-dimensional fracture modelling of heterogeneous and anisotropic rock material under true triaxial compression","authors":"Hui Li,&nbsp;Shanyong Wang","doi":"10.1016/j.compgeo.2025.107862","DOIUrl":"10.1016/j.compgeo.2025.107862","url":null,"abstract":"<div><div>Underground excavation inevitably disturbs natural in-situ stress equilibrium of rocks, leading to localized deformation, fracture initiation, and potential large-scale failure. Understanding the failure mechanisms of rock under complex stress states is of great importance for ensuring the safety and stability of underground structures. To this end, a mesoscale numerical approach combining Weibull random fields with a hybrid phase-field model is adopted to investigate the failure mechanisms of heterogeneous rock materials under true triaxial compression. The isotropic model is further extended to anisotropic formulation by incorporating an anisotropic crack surface density function to account for directional fracture behaviour. First, the isotropic model is calibrated against experimental data by simulating triaxial compression of intact cuboid sandstone specimens, with the effects of confining stress and mesoscale heterogeneity systematically analysed. Next, the influence of pre-existing notches on 3D crack propagation is studied using the calibrated model. Finally, the anisotropic model is used to simulate the failure of intact and notched laminated rocks. It is found that rock fracture behaviour is governed by the combined effects of confining stress, mesoscale heterogeneity, and notch geometry. The incorporation of anisotropic fracture toughness enables directional control of crack propagation along bedding planes, successfully capturing joint-induced compression–shear failure. The developed approach provides a robust tool for simulating 3D failure in heterogeneous and anisotropic rocks under complex loading conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107862"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883612","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
Transient responses of multilayered transversely isotropic half-space subjected to a rectangular oblique impact loading 多层横向各向同性半空间在矩形斜冲击载荷作用下的瞬态响应
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107875
Ye Yang , Mincai Jia , Qun Yan
To address the oblique impact effect in engineering, this paper develops a theoretical solution for non-axisymmetric transient problems based on a three-dimensional Cartesian coordinate system. By combining the dynamic differential equations with displacement potential functions, the expressions for the displacement responses of single-layered soil are derived via the double Fourier-Laplace transformation. The global stiffness matrix for a multilayered half-space is assembled from the analytically derived layer-element stiffness matrices. Transient displacement and stress responses in the time domain are then obtained through inverse transform algorithms. The accuracy of the proposed method is validated by comparing the solutions with the calculation results from finite difference method. Furthermore, the effects of transversely isotropic parameters, loading embedment depth, and impact inclination angle on the dynamic response of the soil are studied using the presented approach.
针对工程中存在的斜冲击效应,提出了一种基于三维直角坐标系的非轴对称瞬态问题的理论解。将动力微分方程与位移势函数相结合,通过双傅里叶-拉普拉斯变换导出了单层土的位移响应表达式。将解析导出的层单元刚度矩阵组合成多层半空间的整体刚度矩阵。然后通过逆变换算法得到时域内的瞬态位移和应力响应。通过与有限差分法计算结果的比较,验证了所提方法的准确性。在此基础上,研究了横向各向同性参数、荷载埋置深度和冲击倾角对土体动力响应的影响。
{"title":"Transient responses of multilayered transversely isotropic half-space subjected to a rectangular oblique impact loading","authors":"Ye Yang ,&nbsp;Mincai Jia ,&nbsp;Qun Yan","doi":"10.1016/j.compgeo.2025.107875","DOIUrl":"10.1016/j.compgeo.2025.107875","url":null,"abstract":"<div><div>To address the oblique impact effect in engineering, this paper develops a theoretical solution for non-axisymmetric transient problems based on a three-dimensional Cartesian coordinate system. By combining the dynamic differential equations with displacement potential functions, the expressions for the displacement responses of single-layered soil are derived via the double Fourier-Laplace transformation. The global stiffness matrix for a multilayered half-space is assembled from the analytically derived layer-element stiffness matrices. Transient displacement and stress responses in the time domain are then obtained through inverse transform algorithms. The accuracy of the proposed method is validated by comparing the solutions with the calculation results from finite difference method. Furthermore, the effects of transversely isotropic parameters, loading embedment depth, and impact inclination angle on the dynamic response of the soil are studied using the presented approach.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107875"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883639","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 Eulerian multiphase model for collapse and segregation of bidisperse granular columns 双分散颗粒柱崩塌和偏析的欧拉多相模型
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107877
Lei Niu , Liyan Zhang , Shuchen Cai , Kang He , Xin Chen
Granular collapse (landslide) is generally destructive to marine infrastructure and ecosystems. Although particle size distribution exerts a strong influence on collapse behavior, the specific mechanisms through which polydispersity affects granular collapse are not yet fully understood. In this study, a Eulerian multiphase model is developed to simulate the collapse of the bidisperse granular column. A bidisperse constitutive law that combines kinetic theory of granular flow (KTGF) and frictional dilation/contraction rheology is used to close the granular stresses. The KTGF incorporates anisotropic particle collisions and interstitial fluid drag, and the adopted drag model remains effective across dilute to dense particle volume fractions. The fluid turbulence model further accounts for density stratification arising from non-uniformity in the particle volume fraction. Both monodisperse and bidisperse granular collapse experiments were performed, and the predicted flow fields, runout distances, and segregation patterns show good agreement with the experimental observations. The results indicate that a higher proportion of coarse particles leads to a greater total runout distance and more pronounced segregation during bidisperse column collapse. Moreover, the choice of frictional viscosity model significantly influences the morphological evolution, underscoring the importance of accurately representing dilation and contraction effects. Overall, this study develops a multiphase framework for analyzing bidisperse granular collapse, providing valuable insights for hazard assessment and mitigation.
颗粒崩塌(滑坡)对海洋基础设施和生态系统具有普遍的破坏性。尽管粒径分布对颗粒坍塌行为有很大影响,但多分散性影响颗粒坍塌的具体机制尚未完全了解。本文建立了欧拉多相模型来模拟双分散颗粒柱的崩塌。结合颗粒流动动力学理论(KTGF)和摩擦扩张/收缩流变学的双分散本构律来关闭颗粒应力。KTGF考虑了各向异性颗粒碰撞和间隙流体阻力,所采用的阻力模型在稀颗粒体积分数和致密颗粒体积分数之间仍然有效。流体湍流模型进一步解释了由颗粒体积分数的不均匀性引起的密度分层。进行了单分散和双分散颗粒坍塌实验,预测的流场、跳动距离和偏析模式与实验结果吻合较好。结果表明:粗颗粒比例越高,双分散柱坍塌过程中总跳动距离越大,偏析越明显;此外,摩擦粘度模型的选择显著影响形态演变,强调了准确表征扩张和收缩效应的重要性。总的来说,本研究开发了一个多阶段框架来分析双分散颗粒坍塌,为危害评估和缓解提供了有价值的见解。
{"title":"A Eulerian multiphase model for collapse and segregation of bidisperse granular columns","authors":"Lei Niu ,&nbsp;Liyan Zhang ,&nbsp;Shuchen Cai ,&nbsp;Kang He ,&nbsp;Xin Chen","doi":"10.1016/j.compgeo.2025.107877","DOIUrl":"10.1016/j.compgeo.2025.107877","url":null,"abstract":"<div><div>Granular collapse (landslide) is generally destructive to marine infrastructure and ecosystems. Although particle size distribution exerts a strong influence on collapse behavior, the specific mechanisms through which polydispersity affects granular collapse are not yet fully understood. In this study, a Eulerian multiphase model is developed to simulate the collapse of the bidisperse granular column. A bidisperse constitutive law that combines kinetic theory of granular flow (KTGF) and frictional dilation/contraction rheology is used to close the granular stresses. The KTGF incorporates anisotropic particle collisions and interstitial fluid drag, and the adopted drag model remains effective across dilute to dense particle volume fractions. The fluid turbulence model further accounts for density stratification arising from non-uniformity in the particle volume fraction. Both monodisperse and bidisperse granular collapse experiments were performed, and the predicted flow fields, runout distances, and segregation patterns show good agreement with the experimental observations. The results indicate that a higher proportion of coarse particles leads to a greater total runout distance and more pronounced segregation during bidisperse column collapse. Moreover, the choice of frictional viscosity model significantly influences the morphological evolution, underscoring the importance of accurately representing dilation and contraction effects. Overall, this study develops a multiphase framework for analyzing bidisperse granular collapse, providing valuable insights for hazard assessment and mitigation.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107877"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883611","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 GSA-ML hybrid framework combined key geological parameters selection for deformation prediction of fractured rock slopes 结合关键地质参数选择的GSA-ML混合框架用于裂隙岩质边坡变形预测
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107876
Wenxiu Wang , Huanling Wang , Yizhe Wu , Chenglong Yang , Wei Bao , Yuxuan Liu
The deformation prediction of fractured rock slopes is a key challenge in hydropower engineering due to the heterogeneity and stochastic distribution of fracture systems. Traditional approaches often struggle to identify dominant fracture parameters and maintain predictive accuracy while preserving geological interpretability. This study proposes a hybrid feature selection and prediction framework that integrates Discrete Fracture Network (DFN) modeling, Finite Difference Method (FDM), Global Sensitivity Analysis (GSA), and Extreme Gradient Boosting (XGBoost). Eighteen fracture-related geometric and mechanical parameters were initially considered, with maximum slope displacement from DFN-FDM simulations as the output. Sensitivity ranking was performed using Pearson correlation, Sobol variance-based index, and PAWN density-based method (Probabilistic Analysis of Whisker Numbers). The results indicate that 11 parameters dominate slope deformation, reflecting both geometric configuration and mechanical heterogeneity. Based on ranked features, XGBoost surrogate models were established and evaluated through 10-fold cross-validation. The GSA-XGBoost achieves the lowest error, outperforming both the original model and PCA-based (Principal Component Analysis) surrogate model, while reducing dimensionality and retaining geological interpretability. Application to a fractured slope at BDa Hydropower Station demonstrates the framework’s capacity for parameter prioritization and reliable deformation prediction. This approach provides practical guidance for slope stability evaluation and support design in fractured rock masses.
由于裂隙系统的非均质性和随机性,裂隙边坡的变形预测是水电工程中的一个关键问题。传统方法往往难以确定主要裂缝参数,并在保持地质可解释性的同时保持预测精度。本研究提出了一个混合特征选择和预测框架,该框架集成了离散裂缝网络(DFN)建模、有限差分法(FDM)、全局灵敏度分析(GSA)和极端梯度增强(XGBoost)。最初考虑了18个与裂缝相关的几何和力学参数,并将DFN-FDM模拟的最大边坡位移作为输出。采用Pearson相关性、基于Sobol方差的指数和基于PAWN密度的方法(须数概率分析)进行敏感性排序。结果表明,11个参数主导了边坡变形,反映了边坡的几何形态和力学非均质性。基于排序特征,建立XGBoost代理模型,并通过10次交叉验证对模型进行评价。GSA-XGBoost实现了最低的误差,优于原始模型和基于pca(主成分分析)的代理模型,同时降低了维数并保留了地质可解释性。通过对BDa水电站断裂边坡的应用,验证了该框架的参数优选能力和可靠的变形预测能力。该方法对裂隙岩体边坡稳定性评价和支护设计具有实际指导意义。
{"title":"A GSA-ML hybrid framework combined key geological parameters selection for deformation prediction of fractured rock slopes","authors":"Wenxiu Wang ,&nbsp;Huanling Wang ,&nbsp;Yizhe Wu ,&nbsp;Chenglong Yang ,&nbsp;Wei Bao ,&nbsp;Yuxuan Liu","doi":"10.1016/j.compgeo.2025.107876","DOIUrl":"10.1016/j.compgeo.2025.107876","url":null,"abstract":"<div><div>The deformation prediction of fractured rock slopes is a key challenge in hydropower engineering due to the heterogeneity and stochastic distribution of fracture systems. Traditional approaches often struggle to identify dominant fracture parameters and maintain predictive accuracy while preserving geological interpretability. This study proposes a hybrid feature selection and prediction framework that integrates Discrete Fracture Network (DFN) modeling, Finite Difference Method (FDM), Global Sensitivity Analysis (GSA), and Extreme Gradient Boosting (XGBoost). Eighteen fracture-related geometric and mechanical parameters were initially considered, with maximum slope displacement from DFN-FDM simulations as the output. Sensitivity ranking was performed using Pearson correlation, Sobol variance-based index, and PAWN density-based method (Probabilistic Analysis of Whisker Numbers). The results indicate that 11 parameters dominate slope deformation, reflecting both geometric configuration and mechanical heterogeneity. Based on ranked features, XGBoost surrogate models were established and evaluated through 10-fold cross-validation. The GSA-XGBoost achieves the lowest error, outperforming both the original model and PCA-based (Principal Component Analysis) surrogate model, while reducing dimensionality and retaining geological interpretability. Application to a fractured slope at BDa Hydropower Station demonstrates the framework’s capacity for parameter prioritization and reliable deformation prediction. This approach provides practical guidance for slope stability evaluation and support design in fractured rock masses.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107876"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883641","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
Improved p-y curve modeling for large-diameter monopiles in asymmetric scour conditions in sandy soils 沙质土非对称冲刷条件下大直径单桩的改进p-y曲线建模
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107880
Ben Wu, Siau Chen Chian
Scour is a critical factor affecting the bearing capacity of monopile foundations for offshore wind turbines (OWTs). Influenced by factors such as flow direction, turbulence structures and sediment heterogeneity, the geometry of scour holes often exhibit asymmetry. This study investigates the effects of geometric asymmetry of scour holes on the bearing capacity of large-diameter monopiles through the finite element method (FEM) in sandy soil. A self-developed Python Application Programming Interface (Python-API) program is implemented to automatically extract p-y curves from FEM results. The proposed FEM approach is validated through comparison with existing centrifuge test data. A series of parametric analyses is subsequently conducted, revealing that scour hole asymmetry has a minor influence on monopile deformation and the location of rotation center, but has a pronounced impact on the ultimate soil resistance. Among all geometric parameters, the ultimate soil resistance is most sensitive to variations in depth asymmetry, whereas the effects of width and slope asymmetry are negligible. Under the combined influence of pile diameter, the depth asymmetry continues to exhibit a strong correlation with ultimate soil resistance. Findings from the parametric analyses are adopted to improve the existing p-y curve model, which is subsequently validated against existing centrifuge tests and analytical solutions.
冲刷是影响海上风力发电机组单桩基础承载力的重要因素。受水流方向、湍流结构和泥沙非均质性等因素的影响,冲刷孔的几何形状往往呈现不对称性。采用有限元方法研究了砂质土中冲刷孔几何不对称对大直径单桩承载力的影响。实现了自行开发的Python应用编程接口(Python- api)程序,可自动从有限元结果中提取p-y曲线。通过与现有离心机试验数据的对比,验证了所提出的有限元方法的有效性。随后进行的一系列参数分析表明,冲刷孔不对称对单桩变形和旋转中心位置的影响较小,但对极限土抗力的影响较大。在所有几何参数中,土壤极限阻力对深度不对称最为敏感,而宽度和坡度不对称的影响可以忽略不计。在桩径的综合影响下,桩深不对称继续与极限土阻力表现出较强的相关性。采用参数分析的结果来改进现有的p-y曲线模型,随后根据现有的离心机试验和分析解进行验证。
{"title":"Improved p-y curve modeling for large-diameter monopiles in asymmetric scour conditions in sandy soils","authors":"Ben Wu,&nbsp;Siau Chen Chian","doi":"10.1016/j.compgeo.2025.107880","DOIUrl":"10.1016/j.compgeo.2025.107880","url":null,"abstract":"<div><div>Scour is a critical factor affecting the bearing capacity of monopile foundations for offshore wind turbines (OWTs). Influenced by factors such as flow direction, turbulence structures and sediment heterogeneity, the geometry of scour holes often exhibit asymmetry. This study investigates the effects of geometric asymmetry of scour holes on the bearing capacity of large-diameter monopiles through the finite element method (FEM) in sandy soil. A self-developed Python Application Programming Interface (Python-API) program is implemented to automatically extract <em>p</em>-<em>y</em> curves from FEM results. The proposed FEM approach is validated through comparison with existing centrifuge test data. A series of parametric analyses is subsequently conducted, revealing that scour hole asymmetry has a minor influence on monopile deformation and the location of rotation center, but has a pronounced impact on the ultimate soil resistance. Among all geometric parameters, the ultimate soil resistance is most sensitive to variations in depth asymmetry, whereas the effects of width and slope asymmetry are negligible. Under the combined influence of pile diameter, the depth asymmetry continues to exhibit a strong correlation with ultimate soil resistance. Findings from the parametric analyses are adopted to improve the existing <em>p</em>-<em>y</em> curve model, which is subsequently validated against existing centrifuge tests and analytical solutions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107880"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883640","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
Nonlocal binary-medium constitutive model for quasibrittle materials: A framework with implicit gradient enhancement 准脆性材料的非局部二元介质本构模型:隐式梯度增强框架
IF 6.2 1区 工程技术 Q1 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-31 DOI: 10.1016/j.compgeo.2025.107885
Yanbin Chen , Enlong Liu , Ziyin Cai
To address the pathological mesh dependency inherent in strain-softening simulations of quasibrittle materials, this paper presents a novel implicit gradient-enhanced computational framework for nonlocal binary-medium models. The core innovation lies in the synergistic integration of gradient-enhanced regularization with the multiscale homogenization theory underpinning binary-medium models. This approach preserves the sound physical foundation of the constitutive theory while effectively remedying non-objectivity in localized failure analysis through a Helmholtz-type partial differential equation. Within this regularized framework, a new nonlocal binary-medium constitutive model is developed. It comprises a porous elastic medium governed by Hashin–Shtrikman upper bounds and an elastic–perfectly plastic frictional medium with an innovative yield criterion, homogenized via an extended Mori–Tanaka scheme. Critically, to capture complex breakage modes, the model introduces an equivalent strain measure driven by the Unified Strength Theory. This novel feature enables the model to account for the influences of hydrostatic pressure and the stress Lode angle, which are crucial aspects often overlooked in previous models. Numerical investigations, implemented via a Bubnov–Galerkin scheme, confirm that the gradient-enhanced framework ensures mesh-objective solutions during strain-softening regimes. The model successfully replicates a comprehensive suite of quasibrittle behaviors, including multi-stage stress–strain responses, stiffness degradation, tension-compression asymmetry, and sensitivities to both confining pressure and the Lode angle. This work establishes a theoretically rigorous and computationally robust paradigm for modeling breakage progression, significantly advancing the predictive accuracy for engineering applications involving localized failure and material instability.
为了解决准脆性材料应变软化模拟中固有的病态网格依赖,本文提出了一种新的隐式梯度增强计算框架,用于非局部二元介质模型。核心创新在于梯度增强正则化与支撑二元介质模型的多尺度均匀化理论的协同集成。该方法既保留了本构理论的良好物理基础,又有效地弥补了局部破坏分析中通过亥姆霍兹型偏微分方程进行的非客观性问题。在此正则化框架下,提出了一种新的非局部二元介质本构模型。它包括一个由Hashin-Shtrikman上界控制的多孔弹性介质和一个具有创新屈服准则的弹塑性摩擦介质,通过扩展的Mori-Tanaka方案进行均匀化。关键是,为了捕捉复杂的破坏模式,该模型引入了由统一强度理论驱动的等效应变测量。这一新颖的特征使模型能够考虑静水压力和应力Lode角的影响,这是以前模型中经常忽略的关键方面。通过Bubnov-Galerkin格式实施的数值研究证实,梯度增强框架确保了应变软化过程中的网格目标解。该模型成功地复制了一套全面的准脆性行为,包括多阶段应力-应变响应、刚度退化、拉压不对称以及对围压和Lode角的敏感性。这项工作建立了一个理论上严谨、计算上稳健的断裂过程模型,显著提高了涉及局部破坏和材料不稳定的工程应用的预测精度。
{"title":"Nonlocal binary-medium constitutive model for quasibrittle materials: A framework with implicit gradient enhancement","authors":"Yanbin Chen ,&nbsp;Enlong Liu ,&nbsp;Ziyin Cai","doi":"10.1016/j.compgeo.2025.107885","DOIUrl":"10.1016/j.compgeo.2025.107885","url":null,"abstract":"<div><div>To address the pathological mesh dependency inherent in strain-softening simulations of quasibrittle materials, this paper presents a novel implicit gradient-enhanced computational framework for nonlocal binary-medium models. The core innovation lies in the synergistic integration of gradient-enhanced regularization with the multiscale homogenization theory underpinning binary-medium models. This approach preserves the sound physical foundation of the constitutive theory while effectively remedying non-objectivity in localized failure analysis through a Helmholtz-type partial differential equation. Within this regularized framework, a new nonlocal binary-medium constitutive model is developed. It comprises a porous elastic medium governed by Hashin–Shtrikman upper bounds and an elastic–perfectly plastic frictional medium with an innovative yield criterion, homogenized via an extended Mori–Tanaka scheme. Critically, to capture complex breakage modes, the model introduces an equivalent strain measure driven by the Unified Strength Theory. This novel feature enables the model to account for the influences of hydrostatic pressure and the stress Lode angle, which are crucial aspects often overlooked in previous models. Numerical investigations, implemented via a Bubnov–Galerkin scheme, confirm that the gradient-enhanced framework ensures mesh-objective solutions during strain-softening regimes. The model successfully replicates a comprehensive suite of quasibrittle behaviors, including multi-stage stress–strain responses, stiffness degradation, tension-compression asymmetry, and sensitivities to both confining pressure and the Lode angle. This work establishes a theoretically rigorous and computationally robust paradigm for modeling breakage progression, significantly advancing the predictive accuracy for engineering applications involving localized failure and material instability.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107885"},"PeriodicalIF":6.2,"publicationDate":"2025-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145883613","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
期刊
Computers and Geotechnics
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1