Computers and Geotechnics最新文献_第7页

A grain-based phase-field method for simulating compressive-fatigue fracture in heterogeneous brittle rocks 模拟非均质脆性岩石压缩疲劳断裂的相场方法

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

Computers and Geotechnics

Pub Date : 2025-12-19 DOI: 10.1016/j.compgeo.2025.107848

Ming Wu , Zaobao Liu , Zhan Yu

The numerical prediction of compressive-fatigue fracture in heterogeneous brittle rocks under cyclic loading presents a formidable challenge in computational geomechanics. To address this, we introduce a novel mixed-phase field model, grounded in a grain-based methodology, for simulating the fatigue damage and fracture processes in such materials. The proposed framework incorporates two distinct fatigue degradation functions, specifically designed to capture the evolution of both Mode I (tensile) and Mode II (shear) cracks under cyclic conditions. The model is validated through its successful application to standard benchmark problems, including cyclic tests on a single-edge notched specimen and a homogeneous rock specimen embedding a pre-defined weak inclusion. Moreover, we develop a comprehensive numerical model to replicate the complete failure process of a heterogeneous brittle sandstone subjected to cyclic compressive loading. The simulations capture the characteristic stress–strain hysteresis and the progressive initiation and propagation of complex crack networks. Our results indicate that macroscopic damage accumulation commences only after the applied stress exceeds a critical fatigue threshold, beyond which progressive degradation ensues. A key insight from parametric studies reveals that a reduction in the Mode II critical energy release rate (G_c_II) significantly accelerates the initial decay of the tensile fatigue degradation function under a constant stress amplitude. The findings of this study substantiate the extension of the phase-field method to the critical problem of fatigue life prediction in heterogeneous brittle rocks, offering a robust tool for advancing the understanding of their long-term mechanical behavior.

循环荷载作用下非均质脆性岩石压缩疲劳断裂的数值预测是计算地质力学领域的一个巨大挑战。为了解决这个问题，我们引入了一种新的混合相场模型，该模型以基于晶粒的方法为基础，用于模拟此类材料的疲劳损伤和断裂过程。提出的框架包含两种不同的疲劳退化函数，专门用于捕获循环条件下I型（拉伸）和II型（剪切）裂纹的演变。该模型成功应用于标准基准问题，包括单边缺口试件和嵌入预定义弱包裹体的均匀岩石试件的循环试验。此外，我们建立了一个全面的数值模型来复制非均质脆性砂岩在循环压缩载荷下的完整破坏过程。模拟捕获了典型的应力-应变滞后和复杂裂纹网络的渐进起裂和扩展。我们的研究结果表明，宏观损伤积累仅在外加应力超过临界疲劳阈值后才开始，超过该阈值后会发生渐进式退化。参数化研究的一个关键发现表明，在恒定应力幅值下，II型临界能量释放率（GcII）的降低显著加速了拉伸疲劳退化函数的初始衰减。本研究的发现证实了相场方法在非均质脆性岩石疲劳寿命预测这一关键问题上的推广，为促进对其长期力学行为的理解提供了强有力的工具。

{"title":"A grain-based phase-field method for simulating compressive-fatigue fracture in heterogeneous brittle rocks","authors":"Ming Wu , Zaobao Liu , Zhan Yu","doi":"10.1016/j.compgeo.2025.107848","DOIUrl":"10.1016/j.compgeo.2025.107848","url":null,"abstract":"<div><div>The numerical prediction of compressive-fatigue fracture in heterogeneous brittle rocks under cyclic loading presents a formidable challenge in computational geomechanics. To address this, we introduce a novel mixed-phase field model, grounded in a grain-based methodology, for simulating the fatigue damage and fracture processes in such materials. The proposed framework incorporates two distinct fatigue degradation functions, specifically designed to capture the evolution of both Mode I (tensile) and Mode II (shear) cracks under cyclic conditions. The model is validated through its successful application to standard benchmark problems, including cyclic tests on a single-edge notched specimen and a homogeneous rock specimen embedding a pre-defined weak inclusion. Moreover, we develop a comprehensive numerical model to replicate the complete failure process of a heterogeneous brittle sandstone subjected to cyclic compressive loading. The simulations capture the characteristic stress–strain hysteresis and the progressive initiation and propagation of complex crack networks. Our results indicate that macroscopic damage accumulation commences only after the applied stress exceeds a critical fatigue threshold, beyond which progressive degradation ensues. A key insight from parametric studies reveals that a reduction in the Mode II critical energy release rate (<em>G</em><sub>c</sub><em><sub>II</sub></em>) significantly accelerates the initial decay of the tensile fatigue degradation function under a constant stress amplitude. The findings of this study substantiate the extension of the phase-field method to the critical problem of fatigue life prediction in heterogeneous brittle rocks, offering a robust tool for advancing the understanding of their long-term mechanical behavior.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107848"},"PeriodicalIF":6.2,"publicationDate":"2025-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796888","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 state-dependent generalized plasticity constitutive model for coarse-grained soil-structure interfaces 基于状态的粗粒土-结构界面广义塑性本构模型

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

Computers and Geotechnics

Pub Date : 2025-12-17 DOI: 10.1016/j.compgeo.2025.107844

Rui Pang , Yisen Wang , Kehao Chen , Yang Zhou , Bin Xu

Understanding and modeling the mechanical response of coarse-grained soil-structure interfaces is fundamental to geotechnical engineering, particularly for maintaining the safety and functionality of critical structures, including dams and retaining systems. To characterize the nonlinear mechanical response of these interfaces under monotonic and cyclic loading, this study proposes a state-dependent elastoplastic constitutive model formulated within a generalized plasticity framework incorporating critical state soil mechanics. The shear dilatancy formulation was refined using directional parameters to distinguish between loading and unloading paths, thereby capturing cyclic volumetric contraction and dilation. Furthermore, a new plastic modulus formulation was proposed to describe the stiffness degradation and recovery behavior under cyclic loading. These enhancements enable unified modeling of monotonic and cyclic interface behaviors, including strain hardening/softening, cumulative contraction, and dilation-contraction transitions. To validate the model’s robustness and applicability, three typical interface types—rockfill-concrete, gravel-steel, and cushion layer-extrusion wall—are simulated and compared with experimental data. Results show that the proposed model effectively captures the interface response under various initial densities, normal stresses, and loading histories, providing a reliable tool for constitutive modeling and numerical simulation of complex soil-structure interactions.

理解和模拟粗粒度土壤-结构界面的力学响应是岩土工程的基础，特别是对于维护关键结构的安全和功能，包括水坝和挡土墙系统。为了描述这些界面在单调和循环荷载下的非线性力学响应，本研究提出了一个包含临界状态土力学的广义塑性框架内的状态相关弹塑性本构模型。剪切剪胀公式使用方向参数进行细化，以区分加载和卸载路径，从而捕获循环体积收缩和膨胀。此外，提出了一种新的塑性模量公式来描述循环加载下的刚度退化和恢复行为。这些增强功能使单调和循环界面行为的统一建模成为可能，包括应变硬化/软化、累积收缩和膨胀-收缩转变。为验证模型的鲁棒性和适用性，对堆石料-混凝土、碎石-钢和垫层-挤压墙三种典型界面类型进行了模拟，并与实验数据进行了对比。结果表明，该模型能有效地捕捉不同初始密度、正应力和加载历史下的界面响应，为复杂土-结构相互作用的本构建模和数值模拟提供了可靠的工具。

{"title":"A state-dependent generalized plasticity constitutive model for coarse-grained soil-structure interfaces","authors":"Rui Pang , Yisen Wang , Kehao Chen , Yang Zhou , Bin Xu","doi":"10.1016/j.compgeo.2025.107844","DOIUrl":"10.1016/j.compgeo.2025.107844","url":null,"abstract":"<div><div>Understanding and modeling the mechanical response of coarse-grained soil-structure interfaces is fundamental to geotechnical engineering, particularly for maintaining the safety and functionality of critical structures, including dams and retaining systems. To characterize the nonlinear mechanical response of these interfaces under monotonic and cyclic loading, this study proposes a state-dependent elastoplastic constitutive model formulated within a generalized plasticity framework incorporating critical state soil mechanics. The shear dilatancy formulation was refined using directional parameters to distinguish between loading and unloading paths, thereby capturing cyclic volumetric contraction and dilation. Furthermore, a new plastic modulus formulation was proposed to describe the stiffness degradation and recovery behavior under cyclic loading. These enhancements enable unified modeling of monotonic and cyclic interface behaviors, including strain hardening/softening, cumulative contraction, and dilation-contraction transitions. To validate the model’s robustness and applicability, three typical interface types—rockfill-concrete, gravel-steel, and cushion layer-extrusion wall—are simulated and compared with experimental data. Results show that the proposed model effectively captures the interface response under various initial densities, normal stresses, and loading histories, providing a reliable tool for constitutive modeling and numerical simulation of complex soil-structure interactions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107844"},"PeriodicalIF":6.2,"publicationDate":"2025-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797497","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

Advancing internal erosion analysis through three-dimensional FEM simulation: insights from the Agly river dike 通过三维有限元模拟推进内部侵蚀分析：来自Agly河堤防的见解

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

Computers and Geotechnics

Pub Date : 2025-12-17 DOI: 10.1016/j.compgeo.2025.107850

Zezhi Deng , Laurence Girolami , Nadia Benahmed , Pierre Philippe , Stéphane Bonelli , Gang Wang

Internal erosion is a primary cause of degradation and failure in hydraulic structures and natural deposits. While the finite element method (FEM) is widely adopted in field-scale investigations, most existing studies are predominantly limited to two-dimensional (2D) cross-sectional analyses. This simplification proves particularly inadequate when examining complex stratigraphic configurations or specific erosion phenomena including concentrated leaks and sand-boils. This study comprehensively investigates a river dike in southern France. By integrating geophysical imaging data, we developed both a large-scale three-dimensional (3D) geometry model and a series of representative 2D cross-sectional models. A transient FEM that couples multi-field and multi-phase suffusion processes was employed to simulate seepage and fines migration within the dike foundation under realistic flood conditions. Comparative results demonstrate non-negligible discrepancies between 2D and 3D simulations. The 2D analysis demonstrates significant underestimation of suffusion effects by neglecting transverse seepage, while also failing to capture complex flow dynamics inherent in heterogeneous stratigraphic conditions. In contrast, the 3D simulation results aligned closely with field observations, intuitively elucidating the formation mechanisms of sinkholes, leaks, and sand boils while providing clear insights into their complex spatial distribution patterns. The 3D analysis further reveals the complex seepage and suffusion evolution governed by stratigraphic morphology and highlights the interplay between subsurface suffusion and surface signatures of internal erosion. These findings underscore the superiority and necessity of 3D modeling in internal erosion studies for sites featuring complex geological stratification.

内部侵蚀是水工建筑物和自然沉积物退化和破坏的主要原因。虽然有限单元法（FEM）在现场调查中被广泛采用，但大多数现有研究主要局限于二维（2D）截面分析。在研究复杂的地层配置或特定的侵蚀现象（包括集中泄漏和砂沸腾）时，这种简化被证明是特别不充分的。这项研究全面调查了法国南部的一条河堤。通过整合地球物理成像数据，我们建立了大尺度三维（3D）几何模型和一系列具有代表性的二维截面模型。采用耦合多场多相渗流过程的瞬态有限元法，模拟了现实洪水条件下堤基内的渗流和细粒运移。对比结果表明，2D和3D模拟之间存在不可忽略的差异。二维分析表明，由于忽略横向渗流，严重低估了渗透效应，同时也未能捕捉到非均质地层条件下固有的复杂流动动力学。相比之下，3D模拟结果与现场观测结果密切相关，直观地阐明了天坑、泄漏和砂沸腾的形成机制，同时清晰地了解了它们复杂的空间分布模式。三维分析进一步揭示了受地层形态控制的复杂渗流和渗透演化，并强调了地下渗透与地表内部侵蚀特征之间的相互作用。这些发现强调了三维建模在复杂地质分层遗址内部侵蚀研究中的优越性和必要性。

{"title":"Advancing internal erosion analysis through three-dimensional FEM simulation: insights from the Agly river dike","authors":"Zezhi Deng , Laurence Girolami , Nadia Benahmed , Pierre Philippe , Stéphane Bonelli , Gang Wang","doi":"10.1016/j.compgeo.2025.107850","DOIUrl":"10.1016/j.compgeo.2025.107850","url":null,"abstract":"<div><div>Internal erosion is a primary cause of degradation and failure in hydraulic structures and natural deposits. While the finite element method (FEM) is widely adopted in field-scale investigations, most existing studies are predominantly limited to two-dimensional (2D) cross-sectional analyses. This simplification proves particularly inadequate when examining complex stratigraphic configurations or specific erosion phenomena including concentrated leaks and sand-boils. This study comprehensively investigates a river dike in southern France. By integrating geophysical imaging data, we developed both a large-scale three-dimensional (3D) geometry model and a series of representative 2D cross-sectional models. A transient FEM that couples multi-field and multi-phase suffusion processes was employed to simulate seepage and fines migration within the dike foundation under realistic flood conditions. Comparative results demonstrate non-negligible discrepancies between 2D and 3D simulations. The 2D analysis demonstrates significant underestimation of suffusion effects by neglecting transverse seepage, while also failing to capture complex flow dynamics inherent in heterogeneous stratigraphic conditions. In contrast, the 3D simulation results aligned closely with field observations, intuitively elucidating the formation mechanisms of sinkholes, leaks, and sand boils while providing clear insights into their complex spatial distribution patterns. The 3D analysis further reveals the complex seepage and suffusion evolution governed by stratigraphic morphology and highlights the interplay between subsurface suffusion and surface signatures of internal erosion. These findings underscore the superiority and necessity of 3D modeling in internal erosion studies for sites featuring complex geological stratification.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107850"},"PeriodicalIF":6.2,"publicationDate":"2025-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796886","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 multiscale anisotropic damage model for discontinuous rock materials under triaxial compression 三轴压缩下非连续岩体多尺度各向异性损伤模型

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

Computers and Geotechnics

Pub Date : 2025-12-16 DOI: 10.1016/j.compgeo.2025.107845

Jin Zhang , Ying Xia , Ke Ren , Shu Zhu , Qi-Zhi Zhu , Jianfu Shao

This study proposes a thermodynamically consistent constitutive modeling framework for the unified characterization of anisotropic damage in quasi-brittle rock materials containing structural discontinuities, such as cracks and bedding planes. The model integrates an energy release rate-driven damage evolution law with a confinement-sensitive Mohr–Coulomb-type friction criterion, establishing a mechanistic link between microscale inelastic processes and the emergent macroscopic degradation behavior. Anisotropy arising from the orientation and distribution of discontinuities is rigorously captured through directionally dependent second-order damage tensors. Furthermore, a dual-phase decomposition of the Helmholtz free energy function into elastic and dissipative components ensures strict thermodynamic admissibility and internal consistency. The predictive capability of the model is evaluated through three representative applications: bedded yellow sandstone, interlayered sandstone, and fractured rock samples with pre-existing cracks. The results demonstrate the model’s robustness in capturing nonlinear stress–strain responses, strength anisotropy, confinement-induced dilatancy transitions, and failure mode evolution across varying structural configurations. The framework also elucidates the causal relationships between microstructural evolution and macroscopic mechanical behavior, thereby providing a unified approach for modeling damage-induced anisotropy in fractured and layered rock systems under complex loading conditions.

本研究提出了一个热力学一致的本构模型框架，用于统一表征含有裂缝和层理面等结构不连续的准脆性岩石材料的各向异性损伤。该模型将能量释放率驱动的损伤演化规律与约束敏感的mohr - coulomb型摩擦准则相结合，建立了微观尺度非弹性过程与紧急宏观退化行为之间的机制联系。通过方向相关的二阶损伤张量严格捕获了由不连续点的方向和分布引起的各向异性。此外，亥姆霍兹自由能函数的双相分解为弹性和耗散分量，确保了严格的热力学容许性和内部一致性。该模型的预测能力通过三种典型应用进行了评估：层状黄砂岩、层间砂岩和带有预先存在裂缝的破碎岩石样品。结果表明，该模型在捕获非线性应力-应变响应、强度各向异性、约束诱导的剪胀转变和不同结构构型的破坏模式演化方面具有鲁棒性。该框架还阐明了微观结构演化与宏观力学行为之间的因果关系，从而为复杂加载条件下裂隙和层状岩石系统损伤诱导的各向异性建模提供了统一的方法。

{"title":"A multiscale anisotropic damage model for discontinuous rock materials under triaxial compression","authors":"Jin Zhang , Ying Xia , Ke Ren , Shu Zhu , Qi-Zhi Zhu , Jianfu Shao","doi":"10.1016/j.compgeo.2025.107845","DOIUrl":"10.1016/j.compgeo.2025.107845","url":null,"abstract":"<div><div>This study proposes a thermodynamically consistent constitutive modeling framework for the unified characterization of anisotropic damage in quasi-brittle rock materials containing structural discontinuities, such as cracks and bedding planes. The model integrates an energy release rate-driven damage evolution law with a confinement-sensitive Mohr–Coulomb-type friction criterion, establishing a mechanistic link between microscale inelastic processes and the emergent macroscopic degradation behavior. Anisotropy arising from the orientation and distribution of discontinuities is rigorously captured through directionally dependent second-order damage tensors. Furthermore, a dual-phase decomposition of the Helmholtz free energy function into elastic and dissipative components ensures strict thermodynamic admissibility and internal consistency. The predictive capability of the model is evaluated through three representative applications: bedded yellow sandstone, interlayered sandstone, and fractured rock samples with pre-existing cracks. The results demonstrate the model’s robustness in capturing nonlinear stress–strain responses, strength anisotropy, confinement-induced dilatancy transitions, and failure mode evolution across varying structural configurations. The framework also elucidates the causal relationships between microstructural evolution and macroscopic mechanical behavior, thereby providing a unified approach for modeling damage-induced anisotropy in fractured and layered rock systems under complex loading conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107845"},"PeriodicalIF":6.2,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797494","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

Investigation of morphology-dependent breakage behavior in graded aggregates under traffic loading via DEM 基于DEM的交通荷载作用下级配骨料形态依赖破碎行为研究

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

Computers and Geotechnics

Pub Date : 2025-12-16 DOI: 10.1016/j.compgeo.2025.107843

Xianpu Xiao , Hongchuang Zhang , Kang Xie , Zhibo Cheng , Qian Zhang , Taifeng Li , Jianying Zheng

The angular breakage (angular chipping and milling) within graded aggregates significantly weakens the mechanical performance of high-speed railway subgrades and exacerbates mud-pumping. Nevertheless, conventional breakage simulation mainly focuses on bulk breakage, making accurate simulation of angular breakage a significant challenge. To address this issue, an efficient morphology-dependent angular breakage algorithm is proposed within the discrete element method (DEM) framework. Specifically, complex particle morphology of convex particles is precisely quantified using equivalent ellipsoids derived from the inertia tensor. The resulting shape parameter is incorporated into a contact strain criterion and stiffness degradation model to simulate angular breakage. Furthermore, the accuracy of the algorithm is validated by single-particle breakage tests and X-CT assisted cycle loading experiments. Finally, the morphology-dependent breakage mechanism of graded aggregates is revealed through the angular breakage algorithm. The evolution of macro- and meso-indicators reveals a morphology-dependent breakage path comprising sliding limitation, rotation dominance, contact force concentration, and angular breakage. This leads to a critical convex particle morphology (i.e., aspect ratio A_r = 0.8 and angularity index A_i = 0.7) in graded aggregates at which angular breakage is most sensitive. This study provides new insights into the morphology-dependent breakage of graded aggregates in high-speed railway subgrades.

级配骨料的角破碎（角碎裂和角磨）显著削弱了高速铁路路基的力学性能，加剧了抽泥。然而，传统的破碎模拟主要集中在体破碎上，这使得角破碎的精确模拟成为一个重大挑战。为了解决这一问题，在离散元法（DEM）框架下提出了一种高效的基于形态的角断裂算法。具体地说，利用由惯性张量导出的等效椭球精确量化了凸粒子的复杂粒子形态。将得到的形状参数纳入接触应变准则和刚度退化模型来模拟角断裂。通过单颗粒破碎试验和X-CT辅助循环加载试验验证了算法的准确性。最后，通过角破碎算法揭示了级配骨料的形态依赖破碎机制。宏观和中观指标的演化揭示了一个由滑动限制、旋转优势、接触力集中和角度破坏组成的形态依赖的破坏路径。这导致了一个关键的凸颗粒形态（即长径比Ar = 0.8和角度指数Ai = 0.7）在分级聚集体中，角断裂是最敏感的。该研究为高速铁路路基级配骨料的形态依赖性破碎提供了新的认识。

{"title":"Investigation of morphology-dependent breakage behavior in graded aggregates under traffic loading via DEM","authors":"Xianpu Xiao , Hongchuang Zhang , Kang Xie , Zhibo Cheng , Qian Zhang , Taifeng Li , Jianying Zheng","doi":"10.1016/j.compgeo.2025.107843","DOIUrl":"10.1016/j.compgeo.2025.107843","url":null,"abstract":"<div><div>The angular breakage (angular chipping and milling) within graded aggregates significantly weakens the mechanical performance of high-speed railway subgrades and exacerbates mud-pumping. Nevertheless, conventional breakage simulation mainly focuses on bulk breakage, making accurate simulation of angular breakage a significant challenge. To address this issue, an efficient morphology-dependent angular breakage algorithm is proposed within the discrete element method (DEM) framework. Specifically, complex particle morphology of convex particles is precisely quantified using equivalent ellipsoids derived from the inertia tensor. The resulting shape parameter is incorporated into a contact strain criterion and stiffness degradation model to simulate angular breakage. Furthermore, the accuracy of the algorithm is validated by single-particle breakage tests and X-CT assisted cycle loading experiments. Finally, the morphology-dependent breakage mechanism of graded aggregates is revealed through the angular breakage algorithm. The evolution of macro- and <em>meso</em>-indicators reveals a morphology-dependent breakage path comprising sliding limitation, rotation dominance, contact force concentration, and angular breakage. This leads to a critical convex particle morphology (i.e., aspect ratio <em>A<sub>r</sub></em> = 0.8 and angularity index <em>A<sub>i</sub></em> = 0.7) in graded aggregates at which angular breakage is most sensitive. This study provides new insights into the morphology-dependent breakage of graded aggregates in high-speed railway subgrades.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107843"},"PeriodicalIF":6.2,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797496","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

Probabilistic modeling of microwave thermal damage at rough rock-cement interfaces using a Bayesian-ANN framework 基于贝叶斯-神经网络框架的粗糙岩石-水泥界面微波热损伤概率建模

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

Computers and Geotechnics

Pub Date : 2025-12-16 DOI: 10.1016/j.compgeo.2025.107809

Jia-Le Li , Fuxin Rui , Yanlong Zheng , Gao-Feng Zhao , Kostas Senetakis

We examine the thermal damage behavior at the interface of two-phase materials (rock-cement) under microwave exposure, conducting experiments and developing a probabilistic-based methodology. The evolution of the deformation field at the rock-cement interface under varying microwave exposure durations are analyzed using digital image correlation, with a particular focus on the influence of surface morphology on interfacial strength. Three roughness characterization schemes are evaluated, and a 3D roughness metric considering embedding depth and direction is selected to correlate surface morphology with interfacial strength. Owing to inherent uncertainties, a probabilistic framework is subsequently proposed based on the Bayesian approach to identify the optimal thermal damage model for rough rock-cement interfaces. The model incorporates temperature, surface roughness, and interfacial strength. The results demonstrate that a quadratic polynomial of the temperature term, combined with the roughness index, yields the highest evidence value and is therefore selected as the most suitable model to describe the damage behavior. Subsequently, the corresponding optimal parameter distributions and the confidence interval are inferred using the Markov Chain Monte Carlo (MCMC) algorithm within the selected model. By integrating the Bayesian approach with uncertainty derived from synthetic data generated by an artificial neural network (ANN), the robustness and reliability of the selected damage model are further validated through the probability value. Overall, the proposed framework leverages the strengths of the Bayesian approach in handling uncertainty, while integrating data generated by ANN into the inference process, providing a more comprehensive and robust approach to determining the thermal damage model in terms of limited data.

我们研究了微波照射下两相材料（岩石-水泥）界面的热损伤行为，进行了实验并开发了一种基于概率的方法。利用数字图像相关分析了不同微波辐照时间下岩石-水泥界面变形场的演化，重点研究了表面形貌对界面强度的影响。评估了三种粗糙度表征方案，并选择了考虑嵌入深度和方向的三维粗糙度度量来将表面形貌与界面强度相关联。由于岩石-水泥粗糙界面存在固有的不确定性，提出了基于贝叶斯方法的概率框架来识别岩石-水泥粗糙界面的最优热损伤模型。该模型包含温度、表面粗糙度和界面强度。结果表明，温度项的二次多项式与粗糙度指数相结合产生的证据值最高，因此被选为最适合描述损伤行为的模型。随后，在选择的模型中，使用马尔可夫链蒙特卡罗（MCMC）算法推断出相应的最优参数分布和置信区间。将贝叶斯方法与人工神经网络生成的合成数据的不确定性相结合，通过概率值进一步验证所选损伤模型的鲁棒性和可靠性。总体而言，所提出的框架利用了贝叶斯方法在处理不确定性方面的优势，同时将人工神经网络生成的数据集成到推理过程中，提供了一种更全面、更稳健的方法来根据有限的数据确定热损伤模型。

{"title":"Probabilistic modeling of microwave thermal damage at rough rock-cement interfaces using a Bayesian-ANN framework","authors":"Jia-Le Li , Fuxin Rui , Yanlong Zheng , Gao-Feng Zhao , Kostas Senetakis","doi":"10.1016/j.compgeo.2025.107809","DOIUrl":"10.1016/j.compgeo.2025.107809","url":null,"abstract":"<div><div>We examine the thermal damage behavior at the interface of two-phase materials (rock-cement) under microwave exposure, conducting experiments and developing a probabilistic-based methodology. The evolution of the deformation field at the rock-cement interface under varying microwave exposure durations are analyzed using digital image correlation, with a particular focus on the influence of surface morphology on interfacial strength. Three roughness characterization schemes are evaluated, and a 3D roughness metric considering embedding depth and direction is selected to correlate surface morphology with interfacial strength. Owing to inherent uncertainties, a probabilistic framework is subsequently proposed based on the Bayesian approach to identify the optimal thermal damage model for rough rock-cement interfaces. The model incorporates temperature, surface roughness, and interfacial strength. The results demonstrate that a quadratic polynomial of the temperature term, combined with the roughness index, yields the highest evidence value and is therefore selected as the most suitable model to describe the damage behavior. Subsequently, the corresponding optimal parameter distributions and the confidence interval are inferred using the Markov Chain Monte Carlo (MCMC) algorithm within the selected model. By integrating the Bayesian approach with uncertainty derived from synthetic data generated by an artificial neural network (ANN), the robustness and reliability of the selected damage model are further validated through the probability value. Overall, the proposed framework leverages the strengths of the Bayesian approach in handling uncertainty, while integrating data generated by ANN into the inference process, providing a more comprehensive and robust approach to determining the thermal damage model in terms of limited data.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107809"},"PeriodicalIF":6.2,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797500","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

Thermo-plastic response of energy piles during long-term monotonic cooling 长期单调冷却过程中能量桩的热塑性响应

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

Computers and Geotechnics

Pub Date : 2025-12-15 DOI: 10.1016/j.compgeo.2025.107839

Mouadh Rafai , Merita Tafili , Yuepeng Dong , Philip J. Vardon

Realistic numerical modeling of energy piles in soft soil requires advanced constitutive relationships capable of capturing the inherent thermo-plastic behavior of the surrounding ground. In this study, a newly developed rate-dependent, thermo-plastic constitutive model, called AVISA-T, is employed within the Plaxis Finite Element (FE) code to simulate the response of a well-instrumented energy pile embedded in multilayered soft soils subjected to thermo-mechanical loading. Following material parameter calibration and model prediction validation using non-isothermal laboratory tests on the soils surrounding the pile, the model was employed to simulate full-scale in-situ tests. In these simulations, the pile was initially subjected to either 0 % or 60 % of its bearing capacity and then exposed to continuous cooling over a period of up to three months. The AVISA-T model effectively reproduces the development of contractive and expansive strains, as well as compressive and tensile stresses that coexist along the pile shaft, including the accumulation of residual strains and stresses. In the absence of axial mechanical load, both residual contractive and expansive strains were observed, accompanied by irreversible uplift of the pile head, primarily attributed to non-uniform, unrecovered temperature changes. Moreover, under higher mechanical loading, the model captures dragdown effects resulting from thermal shrinkage of the surrounding soil, which contributes to the accumulation of permanent strains, stresses, and settlements. A comparison between simulations using the common Modified Cam Clay (MCC) model, the AVISA model without thermal effects and the AVISA-T model highlights the importance of using models including thermal plasticity for engineering practice.

软土中能源桩的真实数值模拟需要先进的本构关系，能够捕捉周围地基固有的热塑性行为。在这项研究中，一种新开发的速率相关的热塑性本构模型，称为AVISA-T，在Plaxis有限元（FE）代码中使用，以模拟嵌入多层软土中仪器完备的能源桩在热机械载荷下的响应。通过对桩周土体进行非等温室内试验，对材料参数进行校正并对模型进行预测验证后，采用该模型进行了全尺寸现场模拟试验。在这些模拟中，桩最初受到0%或60%的承载力，然后暴露在长达三个月的持续冷却中。AVISA-T模型有效地再现了沿桩身的收缩应变和膨胀应变的发展，以及同时存在的压应力和拉应力，包括残余应变和应力的积累。在不存在轴向机械载荷的情况下，观察到残余收缩应变和膨胀应变，并伴有不可逆的桩顶隆起，主要原因是温度变化不均匀，无法恢复。此外，在较高的机械载荷下，该模型捕获了由周围土壤的热收缩引起的拖拽效应，这有助于永久应变、应力和沉降的积累。通过对常用的修正Cam Clay （MCC）模型、不含热效应的AVISA模型和AVISA- t模型的仿真比较，强调了在工程实践中使用包括热塑性在内的模型的重要性。

{"title":"Thermo-plastic response of energy piles during long-term monotonic cooling","authors":"Mouadh Rafai , Merita Tafili , Yuepeng Dong , Philip J. Vardon","doi":"10.1016/j.compgeo.2025.107839","DOIUrl":"10.1016/j.compgeo.2025.107839","url":null,"abstract":"<div><div>Realistic numerical modeling of energy piles in soft soil requires advanced constitutive relationships capable of capturing the inherent thermo-plastic behavior of the surrounding ground. In this study, a newly developed rate-dependent, thermo-plastic constitutive model, called AVISA-T, is employed within the Plaxis Finite Element (FE) code to simulate the response of a well-instrumented energy pile embedded in multilayered soft soils subjected to thermo-mechanical loading. Following material parameter calibration and model prediction validation using non-isothermal laboratory tests on the soils surrounding the pile, the model was employed to simulate full-scale in-situ tests. In these simulations, the pile was initially subjected to either 0 % or 60 % of its bearing capacity and then exposed to continuous cooling over a period of up to three months. The AVISA-T model effectively reproduces the development of contractive and expansive strains, as well as compressive and tensile stresses that coexist along the pile shaft, including the accumulation of residual strains and stresses. In the absence of axial mechanical load, both residual contractive and expansive strains were observed, accompanied by irreversible uplift of the pile head, primarily attributed to non-uniform, unrecovered temperature changes. Moreover, under higher mechanical loading, the model captures dragdown effects resulting from thermal shrinkage of the surrounding soil, which contributes to the accumulation of permanent strains, stresses, and settlements. A comparison between simulations using the common Modified Cam Clay (MCC) model, the AVISA model without thermal effects and the AVISA-T model highlights the importance of using models including thermal plasticity for engineering practice.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107839"},"PeriodicalIF":6.2,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796890","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 novel fluid-solid coupling method for simulating fracture propagation induced by detonation gas 一种新的流固耦合方法模拟爆轰气体诱导的断裂扩展

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

Computers and Geotechnics

Pub Date : 2025-12-15 DOI: 10.1016/j.compgeo.2025.107840

Huhu Wan, Qinglei Yu, Yuheng Wang, Yufei Jia, Xuerui Yang, Jiangyong Pu

Detonation gas plays a vital role in driving fracture propagation during the quasi-static stage of blasting. However, its transient nature and strong coupling with fractures make it difficult to observe directly. To capture this process, a numerical method that couples FiPy and PFC3D is developed to simulate gas flow and dynamic loading along evolving cracks. The model employs Python API to achieve bidirectional fluid–solid coupling: gas pressure fields are solved using FiPy, while crack propagation and stress responses are tracked by PFC3D. A fracture-feedback mechanism is incorporated to dynamically update permeability and loading within evolving fracture networks. The proposed method was validated against laboratory-scale blast tests on a polycarbonate specimen, demonstrating its accuracy in reproducing both stress wave propagation and gas-induced fracturing. On this basis, a parametric study was conducted to investigate rock mass responses under varying initial gas pressures. Results show that gas pressure has a decisive influence on fracture extent and failure mode. Tensile cracking remains dominant across all cases, while increasing pressure leads to more extensive damage, with fractures exhibiting multidirectional growth and branching. Compared to conventional approaches, the developed model provides a more realistic representation of gas-driven crack propagation, offering new insights for understanding rock fracturing mechanisms and optimizing blast design.

在准静态爆破阶段，爆轰气体对裂缝扩展起着至关重要的推动作用。然而，由于其瞬态性和与裂缝的强耦合性，使其难以直接观测。为了捕捉这一过程，开发了一种耦合FiPy和PFC3D的数值方法来模拟沿裂纹演化的气体流动和动载荷。模型采用Python API实现双向流固耦合，气体压力场采用FiPy求解，裂纹扩展和应力响应采用PFC3D跟踪。采用裂缝反馈机制，在不断变化的裂缝网络中动态更新渗透率和载荷。该方法在实验室规模的聚碳酸酯试样上进行了爆炸试验，验证了其在再现应力波传播和气致压裂方面的准确性。在此基础上，对不同初始气体压力下的岩体响应进行了参数化研究。结果表明，气体压力对断裂程度和破坏模式有决定性影响。在所有情况下，拉伸裂缝仍占主导地位，而压力增加会导致更广泛的破坏，裂缝呈现多向生长和分支。与传统方法相比，所开发的模型能够更真实地描述气致裂缝扩展，为理解岩石破裂机制和优化爆破设计提供了新的见解。

{"title":"A novel fluid-solid coupling method for simulating fracture propagation induced by detonation gas","authors":"Huhu Wan, Qinglei Yu, Yuheng Wang, Yufei Jia, Xuerui Yang, Jiangyong Pu","doi":"10.1016/j.compgeo.2025.107840","DOIUrl":"10.1016/j.compgeo.2025.107840","url":null,"abstract":"<div><div>Detonation gas plays a vital role in driving fracture propagation during the quasi-static stage of blasting. However, its transient nature and strong coupling with fractures make it difficult to observe directly. To capture this process, a numerical method that couples FiPy and PFC3D is developed to simulate gas flow and dynamic loading along evolving cracks. The model employs Python API to achieve bidirectional fluid–solid coupling: gas pressure fields are solved using FiPy, while crack propagation and stress responses are tracked by PFC3D. A fracture-feedback mechanism is incorporated to dynamically update permeability and loading within evolving fracture networks. The proposed method was validated against laboratory-scale blast tests on a polycarbonate specimen, demonstrating its accuracy in reproducing both stress wave propagation and gas-induced fracturing. On this basis, a parametric study was conducted to investigate rock mass responses under varying initial gas pressures. Results show that gas pressure has a decisive influence on fracture extent and failure mode. Tensile cracking remains dominant across all cases, while increasing pressure leads to more extensive damage, with fractures exhibiting multidirectional growth and branching. Compared to conventional approaches, the developed model provides a more realistic representation of gas-driven crack propagation, offering new insights for understanding rock fracturing mechanisms and optimizing blast design.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107840"},"PeriodicalIF":6.2,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796884","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 responses under reverse faulting: A semi-theoretical approach accounting for fault plane effects in soil-tunnel Interactions 反向断层作用下的隧道响应：考虑土-隧道相互作用中断层面效应的半理论方法

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

Computers and Geotechnics

Pub Date : 2025-12-15 DOI: 10.1016/j.compgeo.2025.107817

Guanxiong ZENG , Yasuo SAWAMURA , Ping GENG , Kiyoshi KISHIDA

Increasing numbers of underground infrastructures have been constructed or planned in seismically active and fault-prone regions, resulting in many inevitably intersecting active fault zones. Conventional theoretical approaches have typically simplified the active fault-tunnel system using a beam-foundation model, in which the surrounding soil is represented by equivalent soil springs with uniform stiffness. However, the present study demonstrated that, governed by the tunnel-fault intersection angle and interactions between the hanging wall and footwall, the stiffness of the equivalent soil springs varies significantly near the tunnel-fault intersection, showing a “hook effect”. To capture this effect, spatially varying stiffness distributions were derived through numerical methods and incorporated into a semi-theoretical framework. The active fault-tunnel system was divided into three equivalent sub-systems, and governing equations were formulated and coupled based on the Pasternak elastic foundation beam theory and classical elasticity solutions. Validation of the FE results confirmed the model’s accuracy across diverse geological conditions, with a greatly improved performance over the conventional approach. Furthermore, the proposed semi-theoretical framework was extended to yield two-dimensional results, providing not only consistency with the results of FE simulations, but also intuitive visualizations of tunnel responses under reverse faulting.

越来越多的地下基础设施在地震活跃区和断层易发区建设或规划，导致许多活动断裂带不可避免地相交。传统的理论方法通常使用梁基础模型来简化活动断层-隧道系统，在梁基础模型中，周围的土壤用均匀刚度的等效土弹簧来表示。然而，本研究表明，在隧道-断层相交角和上下盘相互作用的控制下，隧道-断层相交附近等效土弹簧刚度变化显著，表现为“钩效应”。为了捕捉这种效应，通过数值方法推导出空间变化的刚度分布，并将其纳入半理论框架。基于帕斯捷尔纳克弹性地基梁理论和经典弹性解，将活动断洞系统划分为3个等效子系统，建立了控制方程并进行了耦合。有限元结果的验证证实了该模型在不同地质条件下的准确性，比传统方法的性能有了很大的提高。此外，将所提出的半理论框架扩展到二维结果，不仅与有限元模拟结果一致，而且可以直观地显示逆断层作用下隧道的响应。

{"title":"Tunnel responses under reverse faulting: A semi-theoretical approach accounting for fault plane effects in soil-tunnel Interactions","authors":"Guanxiong ZENG , Yasuo SAWAMURA , Ping GENG , Kiyoshi KISHIDA","doi":"10.1016/j.compgeo.2025.107817","DOIUrl":"10.1016/j.compgeo.2025.107817","url":null,"abstract":"<div><div>Increasing numbers of underground infrastructures have been constructed or planned in seismically active and fault-prone regions, resulting in many inevitably intersecting active fault zones. Conventional theoretical approaches have typically simplified the active fault-tunnel system using a beam-foundation model, in which the surrounding soil is represented by equivalent soil springs with uniform stiffness. However, the present study demonstrated that, governed by the tunnel-fault intersection angle and interactions between the hanging wall and footwall, the stiffness of the equivalent soil springs varies significantly near the tunnel-fault intersection, showing a “hook effect”. To capture this effect, spatially varying stiffness distributions were derived through numerical methods and incorporated into a semi-theoretical framework. The active fault-tunnel system was divided into three equivalent sub-systems, and governing equations were formulated and coupled based on the Pasternak elastic foundation beam theory and classical elasticity solutions. Validation of the FE results confirmed the model’s accuracy across diverse geological conditions, with a greatly improved performance over the conventional approach. Furthermore, the proposed semi-theoretical framework was extended to yield two-dimensional results, providing not only consistency with the results of FE simulations, but also intuitive visualizations of tunnel responses under reverse faulting.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107817"},"PeriodicalIF":6.2,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145797498","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

Integrating machine learning and pore-scale physical properties for enhanced permeability prediction of porous rocks 结合机器学习和孔隙尺度物理性质增强多孔岩石渗透率预测

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

Computers and Geotechnics

Pub Date : 2025-12-15 DOI: 10.1016/j.compgeo.2025.107835

Xinyue Zhu , Hongzhuo Fan , Yun Jia , Jianfu Shao

Accurate permeability prediction is critical in subsurface applications. This study developed a hybrid machine-learning framework that integrates image-derived and physical features to improve permeability prediction. Three-dimensional random-field porous-media models were generated using the Stochastic Pore Network Analysis and Modeling (SPAM) method, and two-dimensional images were extracted along the principal axes and processed with a convolutional neural networks (CNNs) to obtain features used as inputs for regression. The corresponding permeabilities were computed via the lattice Boltzmann method (LBM) and used as regression targets. We evaluated support vector regression (SVR) and random forests regression(RFR), SVR performed best and was adopted as the primary regressor for subsequent task. To enhance physical interpretability and generalization, porosity

(ϕ)

and specific surface area

(S_{v})

were incorporated into the image-feature set. Results showed that adding these physical descriptors significantly improved performance on both the training and the test sets, yielding higher explained variance and lower prediction error. Overall, combining image-derived features with physically interpretable descriptors provides an efficient, physically grounded approach to permeability prediction in complex porous materials.

准确的渗透率预测在地下应用中至关重要。该研究开发了一种混合机器学习框架，该框架集成了图像衍生和物理特征，以提高渗透率预测。采用随机孔隙网络分析与建模（SPAM）方法生成三维随机场多孔介质模型，并沿主轴提取二维图像，利用卷积神经网络（cnn）进行处理，获得特征作为回归输入。通过晶格玻尔兹曼方法（LBM）计算相应的渗透率，并将其作为回归目标。我们评估了支持向量回归（SVR）和随机森林回归（RFR），其中SVR表现最好，并被用作后续任务的主要回归因子。为了增强物理可解释性和泛化性，孔隙度（ϕ）和比表面积（Sv）被纳入图像特征集。结果表明，添加这些物理描述符显著提高了训练集和测试集的性能，产生了更高的解释方差和更低的预测误差。总的来说，将图像衍生特征与物理可解释描述符相结合，为复杂多孔材料的渗透率预测提供了一种有效的物理基础方法。

{"title":"Integrating machine learning and pore-scale physical properties for enhanced permeability prediction of porous rocks","authors":"Xinyue Zhu , Hongzhuo Fan , Yun Jia , Jianfu Shao","doi":"10.1016/j.compgeo.2025.107835","DOIUrl":"10.1016/j.compgeo.2025.107835","url":null,"abstract":"<div><div>Accurate permeability prediction is critical in subsurface applications. This study developed a hybrid machine-learning framework that integrates image-derived and physical features to improve permeability prediction. Three-dimensional random-field porous-media models were generated using the Stochastic Pore Network Analysis and Modeling (SPAM) method, and two-dimensional images were extracted along the principal axes and processed with a convolutional neural networks (CNNs) to obtain features used as inputs for regression. The corresponding permeabilities were computed via the lattice Boltzmann method (LBM) and used as regression targets. We evaluated support vector regression (SVR) and random forests regression(RFR), SVR performed best and was adopted as the primary regressor for subsequent task. To enhance physical interpretability and generalization, porosity <span><math><mrow><mo>(</mo><mi>ϕ</mi><mo>)</mo></mrow></math></span> and specific surface area <span><math><mrow><mo>(</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>)</mo></mrow></math></span> were incorporated into the image-feature set. Results showed that adding these physical descriptors significantly improved performance on both the training and the test sets, yielding higher explained variance and lower prediction error. Overall, combining image-derived features with physically interpretable descriptors provides an efficient, physically grounded approach to permeability prediction in complex porous materials.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"191 ","pages":"Article 107835"},"PeriodicalIF":6.2,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145796885","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