Computers and Geotechnics最新文献_第5页

Relating solution tests to pore-scale CaCO3 crystal growth: Numerical simulation based on the phase field method 孔尺度CaCO3晶体生长的相关溶液试验：基于相场法的数值模拟

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

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107947

Guoliang Ma , Yang Xiao , Jinxuan Zhang , Zhichao Zhang

Microbially induced carbonate precipitation (MICP) has been widely studied for its broad application potential in geotechnical engineering. Exploring the pore-scale CaCO₃ precipitation processes is one of the key aspects of fully understanding the mechanism of MICP for enhancing mechanical properties and controlling hydraulic properties of soils and rocks. Additionally, solution tests are typically the first step to understanding the reaction processes of MICP before conducting pore-scale and macro-scale tests. However, most explanations remain qualitative and do not establish a direct link between solution tests, pore-scale processes, and macro-scale performance. In this study, a saturation-dependent kinetic model is developed to describe the chemical reactions involved in MICP, and a phase-field model is employed to simulate the growth of CaCO₃ crystals at the pore scale. The two models are coupled through the CaCO₃ precipitation rate, enabling direct translation of solution test parameters into pore-scale simulations. Model predictions are validated against experimental data from microdroplet tests, showing good agreement in crystal growth dynamics under varying bacterial densities. The results demonstrate the capability of the proposed approach to quantitatively link solution chemistry with pore-scale mineralization processes. This work provides a robust framework for further investigation of MICP-induced pore sealing and particle bridging, contributing to the optimization and design of bio-mediated ground improvement strategies.

微生物诱导碳酸盐沉淀（MICP）因其在岩土工程中的广泛应用潜力而受到广泛的研究。探索孔隙尺度CaCO3的沉淀过程是充分认识MICP增强岩土力学性能和控制水力性能机理的关键方面之一。此外，在进行孔隙尺度和宏观尺度测试之前，溶液测试通常是了解MICP反应过程的第一步。然而，大多数解释仍然是定性的，并没有在溶液测试、孔隙尺度过程和宏观尺度性能之间建立直接联系。在本研究中，建立了一个饱和度相关的动力学模型来描述MICP中涉及的化学反应，并采用相场模型来模拟CaCO3晶体在孔隙尺度上的生长。这两个模型通过CaCO3沉淀速率耦合，可以将溶液测试参数直接转换为孔隙尺度模拟。模型预测与微滴测试的实验数据相对照，显示出在不同细菌密度下晶体生长动力学的良好一致性。结果表明，该方法能够定量地将溶液化学与孔隙尺度的矿化过程联系起来。这项工作为进一步研究micp诱导的孔隙密封和颗粒桥接提供了一个强大的框架，有助于优化和设计生物介导的地面改善策略。

{"title":"Relating solution tests to pore-scale CaCO3 crystal growth: Numerical simulation based on the phase field method","authors":"Guoliang Ma , Yang Xiao , Jinxuan Zhang , Zhichao Zhang","doi":"10.1016/j.compgeo.2026.107947","DOIUrl":"10.1016/j.compgeo.2026.107947","url":null,"abstract":"<div><div>Microbially induced carbonate precipitation (MICP) has been widely studied for its broad application potential in geotechnical engineering. Exploring the pore-scale CaCO<sub>3</sub> precipitation processes is one of the key aspects of fully understanding the mechanism of MICP for enhancing mechanical properties and controlling hydraulic properties of soils and rocks. Additionally, solution tests are typically the first step to understanding the reaction processes of MICP before conducting pore-scale and macro-scale tests. However, most explanations remain qualitative and do not establish a direct link between solution tests, pore-scale processes, and macro-scale performance. In this study, a saturation-dependent kinetic model is developed to describe the chemical reactions involved in MICP, and a phase-field model is employed to simulate the growth of CaCO<sub>3</sub> crystals at the pore scale. The two models are coupled through the CaCO<sub>3</sub> precipitation rate, enabling direct translation of solution test parameters into pore-scale simulations. Model predictions are validated against experimental data from microdroplet tests, showing good agreement in crystal growth dynamics under varying bacterial densities. The results demonstrate the capability of the proposed approach to quantitatively link solution chemistry with pore-scale mineralization processes. This work provides a robust framework for further investigation of MICP-induced pore sealing and particle bridging, contributing to the optimization and design of bio-mediated ground improvement strategies.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107947"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080673","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 graph-based algorithm for the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U) to simulate two-phase flows in naturally fractured porous media using the MPFA-D method on general unstructured meshes 基于图的连续投影嵌入式离散裂缝模型（CpEDFM-U）在一般非结构化网格上采用MPFA-D方法模拟天然裂缝多孔介质中的两相流动

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

Computers and Geotechnics

Pub Date : 2026-01-27 DOI: 10.1016/j.compgeo.2026.107918

Gabriel Martins Cavalcanti Feitosa , Emanoel Rodrigues dos Santos , Pedro Victor Paixão Albuquerque , Artur Castiel Reis de Souza , Darlan Karlo Elisiário de Carvalho , Paulo Roberto Maciel Lyra

Modeling fluid flow in naturally fractured porous media is crucial for applications such as hydrocarbon production and CO2 sequestration. However, accurately simulating these flows remains challenging due to fractures with complex permeability distributions. The Embedded Discrete Fracture Model (EDFM) has been widely used but has limitations in representing fractures acting as barriers, especially in multiphase flows. To overcome these challenges, the Projection-based Embedded Discrete Fracture Model (pEDFM) was developed, offering better handling of fractures with lower permeability than the matrix. However, it can still exhibit discontinuities in fracture projections. To address these limitations, we propose the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U), a graph-based algorithm that guarantees continuous fracture projections in both structured and unstructured 2D meshes. The CpEDFM-U uses Dijkstra’s algorithm to find the shortest path between fracture tips and applies the MPFA-D method for matrix flow and TPFA for fractures. In numerical simulations analyzed, CpEDFM-U outperforms EDFM and pEDFM, demonstrating lower errors and robust performance across different fracture types and mesh resolutions.

模拟天然裂缝多孔介质中的流体流动对于油气生产和二氧化碳封存等应用至关重要。然而，由于具有复杂渗透率分布的裂缝，准确模拟这些流动仍然具有挑战性。嵌入式离散裂缝模型（EDFM）已被广泛应用，但在描述裂缝作为屏障时存在局限性，特别是在多相流中。为了克服这些挑战，开发了基于投影的嵌入式离散裂缝模型（pEDFM），可以更好地处理渗透率较低的裂缝。然而，它仍然可以在裂缝突出处表现出不连续性。为了解决这些限制，我们提出了连续投影嵌入式离散裂缝模型（CpEDFM-U），这是一种基于图的算法，可以保证在结构化和非结构化二维网格中连续的裂缝投影。CpEDFM-U使用Dijkstra算法寻找裂缝尖端之间的最短路径，并对基质流动采用MPFA-D方法，对裂缝采用TPFA方法。在数值模拟分析中，CpEDFM-U优于EDFM和pEDFM，在不同裂缝类型和网格分辨率下显示出更低的误差和更强的性能。

{"title":"A graph-based algorithm for the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U) to simulate two-phase flows in naturally fractured porous media using the MPFA-D method on general unstructured meshes","authors":"Gabriel Martins Cavalcanti Feitosa , Emanoel Rodrigues dos Santos , Pedro Victor Paixão Albuquerque , Artur Castiel Reis de Souza , Darlan Karlo Elisiário de Carvalho , Paulo Roberto Maciel Lyra","doi":"10.1016/j.compgeo.2026.107918","DOIUrl":"10.1016/j.compgeo.2026.107918","url":null,"abstract":"<div><div>Modeling fluid flow in naturally fractured porous media is crucial for applications such as hydrocarbon production and CO2 sequestration. However, accurately simulating these flows remains challenging due to fractures with complex permeability distributions. The Embedded Discrete Fracture Model (EDFM) has been widely used but has limitations in representing fractures acting as barriers, especially in multiphase flows. To overcome these challenges, the Projection-based Embedded Discrete Fracture Model (pEDFM) was developed, offering better handling of fractures with lower permeability than the matrix. However, it can still exhibit discontinuities in fracture projections. To address these limitations, we propose the Continuous-Projection Embedded Discrete Fracture Model (CpEDFM-U), a graph-based algorithm that guarantees continuous fracture projections in both structured and unstructured 2D meshes. The CpEDFM-U uses Dijkstra’s algorithm to find the shortest path between fracture tips and applies the MPFA-D method for matrix flow and TPFA for fractures. In numerical simulations analyzed, CpEDFM-U outperforms EDFM and pEDFM, demonstrating lower errors and robust performance across different fracture types and mesh resolutions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107918"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080672","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

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

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

Computers and Geotechnics

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

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

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

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

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

引用次数: 0

Utilizing physics-informed neural network and geotechnical distance field for solving three-dimensional nonlinear consolidation 利用物理信息神经网络和岩土距离场求解三维非线性固结

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

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107926

Tran-Gia-Khiem Nguyen, Jongmuk Won

Solving three-dimensional (3D) nonlinear consolidation is complex and computationally expensive. This study proposes a framework for solving 3D nonlinear consolidation by utilizing an improved physics-informed neural network with hard constraints coupling with machine learning based geotechnical distance functions for three-dimensional spatial interpolation. The performance of the developed framework was assessed by comparing pore water pressure data between the developed framework and those obtained from COMSOL Multiphysics. In addition, the impact of vertical hydraulic conductibility heterogeneity, compression index, and void ratio on long-term settlement was also evaluated and discussed. It was found that the proposed framework showed a reliable estimation of the 3D distribution of pore water pressure across the 3D domain, achieving results that are comparable to data obtained from COMSOL. In addition, the heterogeneity of hydraulic conductivity can be successfully considered using the developed framework, which enables assessing the long-term settlement of a clay deposit with high uncertainty of hydraulic conductivity. Overall, the developed framework shown in this study can be applied to complex consolidation problems with low computational costs and high accuracy.

求解三维（3D）非线性固结是复杂和计算昂贵的。本研究提出了一个解决三维非线性固结的框架，该框架利用改进的物理信息神经网络与硬约束耦合，基于机器学习的岩土距离函数用于三维空间插值。通过比较开发框架与COMSOL Multiphysics获得的孔隙水压力数据，对开发框架的性能进行了评估。此外，还评价和讨论了竖向导电性非均质性、压缩指数和孔隙比对长期沉降的影响。研究发现，所提出的框架可以可靠地估计孔隙水压力在三维域中的三维分布，其结果与COMSOL获得的数据相当。此外，利用所开发的框架可以成功地考虑水力导电性的非均质性，从而能够评估具有高水力导电性不确定性的粘土沉积物的长期沉降。总的来说，本研究中显示的开发框架可以应用于复杂的固结问题，计算成本低，精度高。

{"title":"Utilizing physics-informed neural network and geotechnical distance field for solving three-dimensional nonlinear consolidation","authors":"Tran-Gia-Khiem Nguyen, Jongmuk Won","doi":"10.1016/j.compgeo.2026.107926","DOIUrl":"10.1016/j.compgeo.2026.107926","url":null,"abstract":"<div><div>Solving three-dimensional (3D) nonlinear consolidation is complex and computationally expensive. This study proposes a framework for solving 3D nonlinear consolidation by utilizing an improved physics-informed neural network with hard constraints coupling with machine learning based geotechnical distance functions for three-dimensional spatial interpolation. The performance of the developed framework was assessed by comparing pore water pressure data between the developed framework and those obtained from COMSOL Multiphysics. In addition, the impact of vertical hydraulic conductibility heterogeneity, compression index, and void ratio on long-term settlement was also evaluated and discussed. It was found that the proposed framework showed a reliable estimation of the 3D distribution of pore water pressure across the 3D domain, achieving results that are comparable to data obtained from COMSOL. In addition, the heterogeneity of hydraulic conductivity can be successfully considered using the developed framework, which enables assessing the long-term settlement of a clay deposit with high uncertainty of hydraulic conductivity. Overall, the developed framework shown in this study can be applied to complex consolidation problems with low computational costs and high accuracy.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107926"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026134","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

High-fidelity digital modeling and comparative morphological analysis of Chang’E-5 and Chang’E-6 lunar grains “嫦娥五号”和“嫦娥六号”月球颗粒高保真数字建模及形态对比分析

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

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107933

Ziwei Tian , Songzheng Yu , Hanyu Wang , Guang Zhang , Yiwei Liu , Yizhuo He , Quan Zheng , Guodong Liu , Xin Liu , Ronghua Pang , Guanghui Liu , Shijing He , Yang Li , Peng Zhang

Using high-resolution X-ray microscopy (XRM), the 3D structures of lunar regolith from the Chang’E-5 and Chang’E-6 missions were reconstructed, and a four-stage CNN-based volumetric segmentation framework was established to enable accurate extraction of tens of thousands of particles from both samples. Quantitative analysis shows that Chang’E-5 particles exhibit more elongated and irregular shapes (Aspect Ratio = 0.612, Elongation Index = 0.727, Flatness Index = 0.845) and higher energy contained in the high-degree spherical harmonic spectrum (0.0346) than Chang’E-6 (Aspect Ratio = 0.623, Elongation Index = 0.732, Flatness Index = 0.854; high-degree spectral energy = 0.0299), indicating rougher surfaces and lower maturity. XRM-derived porosity measurements further reveal a pronounced structural contrast, with porosities of 38.37 % for Chang’E-5 and 54.44 % for Chang’E-6, which leads to a substantial difference in their estimated bearing capacities. These results establish a direct quantitative link between micro-scale grain morphology, macro-scale regolith structure, and surface mechanical behavior, providing critical constraints for landing-site evaluation and lunar surface infrastructure design.

利用高分辨率x射线显微镜（XRM），重建了嫦娥五号和嫦娥六号任务月球风化层的三维结构，并建立了一个基于cnn的四级体积分割框架，以便从两个样本中精确提取数万个颗粒。定量分析表明，与嫦娥6号（长径比= 0.623，延伸指数= 0.732，平整度指数= 0.854，高度光谱能量= 0.0299）相比，嫦娥5号颗粒的形状更拉长、不规则（宽高比= 0.612，延伸指数= 0.727，平整度指数= 0.845），高度球谐波谱能量（0.0346）更高，表面粗糙，成熟度较低。xrm孔隙度测量进一步揭示了明显的结构差异，嫦娥五号和嫦娥六号的孔隙度分别为38.37%和54.44%，这导致了两者估算的承载能力存在很大差异。这些结果在微观尺度的颗粒形态、宏观尺度的风化层结构和表面力学行为之间建立了直接的定量联系，为着陆点评估和月球表面基础设施设计提供了关键约束。

{"title":"High-fidelity digital modeling and comparative morphological analysis of Chang’E-5 and Chang’E-6 lunar grains","authors":"Ziwei Tian , Songzheng Yu , Hanyu Wang , Guang Zhang , Yiwei Liu , Yizhuo He , Quan Zheng , Guodong Liu , Xin Liu , Ronghua Pang , Guanghui Liu , Shijing He , Yang Li , Peng Zhang","doi":"10.1016/j.compgeo.2026.107933","DOIUrl":"10.1016/j.compgeo.2026.107933","url":null,"abstract":"<div><div>Using high-resolution X-ray microscopy (XRM), the 3D structures of lunar regolith from the Chang’E-5 and Chang’E-6 missions were reconstructed, and a four-stage CNN-based volumetric segmentation framework was established to enable accurate extraction of tens of thousands of particles from both samples. Quantitative analysis shows that Chang’E-5 particles exhibit more elongated and irregular shapes (Aspect Ratio = 0.612, Elongation Index = 0.727, Flatness Index = 0.845) and higher energy contained in the high-degree spherical harmonic spectrum (0.0346) than Chang’E-6 (Aspect Ratio = 0.623, Elongation Index = 0.732, Flatness Index = 0.854; high-degree spectral energy = 0.0299), indicating rougher surfaces and lower maturity. XRM-derived porosity measurements further reveal a pronounced structural contrast, with porosities of 38.37 % for Chang’E-5 and 54.44 % for Chang’E-6, which leads to a substantial difference in their estimated bearing capacities. These results establish a direct quantitative link between micro-scale grain morphology, macro-scale regolith structure, and surface mechanical behavior, providing critical constraints for landing-site evaluation and lunar surface infrastructure design.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107933"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026152","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 random simplified analysis method to evaluate braced excavation-induced wall deflection considering spatial variability of soil 考虑土体空间变异性的支挖墙挠度随机简化分析方法

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

Computers and Geotechnics

Pub Date : 2026-01-23 DOI: 10.1016/j.compgeo.2026.107917

Xiaolu Gan , Nianwu Liu , Jin Guo , Wenbo Li , Jianlin Yu , Xiaonan Gong

The wall deflection in braced excavations is inevitably and significantly influenced by the inherent spatial variability of soil parameters, and the applicability of current beam-spring model-based simplified methods becomes severely constricted when considering soil spatial variability. Therefore, a random simplified analysis method for the evaluation of the braced excavation-induced wall deflection in spatially variable soil is proposed in this paper. The deterministic model is established based on an improved beam-spring model that incorporates spatial variations in soil parameters, and its finite difference solution is derived. The spatial variability of soil parameters is characterized by the lognormal random field, which is further discretized using the Karhunen–Loève expansion method. The sparse polynomial chaos expansion is employed to construct an efficient metamodel for rapid random analysis. The effectiveness of the proposed method is demonstrated through a case study of a well-documented braced excavation project. Random analysis results reveal that the variability in the normalized maximum wall deflection initially decreases and then increases, while that in its position increases abruptly before progressively decreasing, due to the coupled effect of the increased excavation depth and the addition of internal supports. The ignorance of the soil spatial variability may result in overestimation of the variability in the normalized deflection and underestimation of that in its position. The presented methodology provides a useful reference for the uncertainty analysis of excavation-induced wall deflection during preliminary design stages.

土体参数固有的空间变异性对支撑式基坑墙体挠度的影响是不可避免的和显著的，在考虑土体空间变异性时，现有基于梁-弹簧模型的简化方法的适用性受到严重限制。因此，本文提出了一种评价空间变土中支撑开挖引起的墙体挠度的随机简化分析方法。在考虑土壤参数空间变化的改进梁-弹簧模型的基础上，建立了确定性模型，并推导了其有限差分解。土壤参数的空间变异性表现为对数正态随机场，并采用karhunen - lo展开法对其进行离散化。利用稀疏多项式混沌展开构造了一个快速随机分析的有效元模型。所提出的方法的有效性通过一个有充分记录的支撑挖掘项目的案例研究来证明。随机分析结果表明，由于开挖深度增加和内支护的耦合作用，归一化最大墙体挠度的变异性先减小后增大，而其位置的变异性先增大后逐渐减小。对土壤空间变异性的忽视可能导致对归一化挠度变异性的高估和对其位置变异性的低估。该方法为初步设计阶段开挖引起的墙体挠度的不确定性分析提供了有益的参考。

{"title":"A random simplified analysis method to evaluate braced excavation-induced wall deflection considering spatial variability of soil","authors":"Xiaolu Gan , Nianwu Liu , Jin Guo , Wenbo Li , Jianlin Yu , Xiaonan Gong","doi":"10.1016/j.compgeo.2026.107917","DOIUrl":"10.1016/j.compgeo.2026.107917","url":null,"abstract":"<div><div>The wall deflection in braced excavations is inevitably and significantly influenced by the inherent spatial variability of soil parameters, and the applicability of current beam-spring model-based simplified methods becomes severely constricted when considering soil spatial variability. Therefore, a random simplified analysis method for the evaluation of the braced excavation-induced wall deflection in spatially variable soil is proposed in this paper. The deterministic model is established based on an improved beam-spring model that incorporates spatial variations in soil parameters, and its finite difference solution is derived. The spatial variability of soil parameters is characterized by the lognormal random field, which is further discretized using the Karhunen–Loève expansion method. The sparse polynomial chaos expansion is employed to construct an efficient metamodel for rapid random analysis. The effectiveness of the proposed method is demonstrated through a case study of a well-documented braced excavation project. Random analysis results reveal that the variability in the normalized maximum wall deflection initially decreases and then increases, while that in its position increases abruptly before progressively decreasing, due to the coupled effect of the increased excavation depth and the addition of internal supports. The ignorance of the soil spatial variability may result in overestimation of the variability in the normalized deflection and underestimation of that in its position. The presented methodology provides a useful reference for the uncertainty analysis of excavation-induced wall deflection during preliminary design stages.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107917"},"PeriodicalIF":6.2,"publicationDate":"2026-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146026135","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

Stochastic simulation of three-dimensional unsaturated flow in water repellent heterogeneous soil 拒水非均质土中三维非饱和流动的随机模拟

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

Computers and Geotechnics

Pub Date : 2026-01-21 DOI: 10.1016/j.compgeo.2026.107938

Evan John Ricketts , Peter John Cleall , Anthony Jefferson , Pierre Kerfriden , Paul Lyons

It has been established that spatial variability of material properties can lead to distinct unstable hydraulic behaviour, and that this is prominent in soils due to their large spatial heterogeneity. This characteristic can lead to non-uniform wetting behaviour and is particularly prominent when the wettability of the medium is also non-uniform. In water repellent soil, its wettability is often spatially varying, such that a network of flow paths is created where fluid can move preferentially, leading to fingered flow patterns. In this study, the development of a model to represent moisture transport in hydrophobic soil is presented. Local spatial variations in material properties are represented by Gaussian random fields as part of a stochastic finite element based model. Key components of the model include an approach to represent the transition region between wettable and non-wettable layers, and the adoption of a suitable saturation–capillary pressure relationship. For wettable soil, this can be achieved with the standard van Genuchten relation. For hydrophobic soil, this is not applicable; thus, an alternative is employed. The model is then validated against field-scale experimental observations by Lipsius and Mooney (2006), which examined the impact of soil heterogeneity on infiltration profiles. The results demonstrate the model’s ability to capture complex flow dynamics in hydrophobic soils, extending the understanding of moisture transport in heterogeneous soils by explicitly modelling the spatial variability of wettability and its impact on soil hydraulic response.

已经确定，材料特性的空间变异性会导致不同的不稳定水力行为，并且由于其巨大的空间异质性，这在土壤中是突出的。这一特性可能导致不均匀的润湿行为，当介质的润湿性也不均匀时，这一特性尤为突出。在疏水性土壤中，其润湿性通常是空间变化的，因此流体可以优先移动的流动路径网络会形成指状流动模式。在这项研究中，提出了一个模型来表示水分在疏水土壤中的输送。材料性能的局部空间变化用高斯随机场表示，作为随机有限元模型的一部分。该模型的关键组成部分包括一种表示可湿层和不可湿层之间过渡区域的方法，以及采用合适的饱和度-毛细管压力关系。对于可湿性土壤，这可以用标准的van Genuchten关系来实现。对于疏水性土壤，这是不适用的；因此，采用了另一种方法。然后，Lipsius和Mooney（2006）根据现场尺度的实验观察验证了该模型，该实验研究了土壤异质性对入渗剖面的影响。结果表明，该模型能够捕捉疏水性土壤中复杂的流动动力学，通过明确模拟润湿性的空间变异性及其对土壤水力响应的影响，扩展了对非均质土壤中水分输送的理解。

{"title":"Stochastic simulation of three-dimensional unsaturated flow in water repellent heterogeneous soil","authors":"Evan John Ricketts , Peter John Cleall , Anthony Jefferson , Pierre Kerfriden , Paul Lyons","doi":"10.1016/j.compgeo.2026.107938","DOIUrl":"10.1016/j.compgeo.2026.107938","url":null,"abstract":"<div><div>It has been established that spatial variability of material properties can lead to distinct unstable hydraulic behaviour, and that this is prominent in soils due to their large spatial heterogeneity. This characteristic can lead to non-uniform wetting behaviour and is particularly prominent when the wettability of the medium is also non-uniform. In water repellent soil, its wettability is often spatially varying, such that a network of flow paths is created where fluid can move preferentially, leading to fingered flow patterns. In this study, the development of a model to represent moisture transport in hydrophobic soil is presented. Local spatial variations in material properties are represented by Gaussian random fields as part of a stochastic finite element based model. Key components of the model include an approach to represent the transition region between wettable and non-wettable layers, and the adoption of a suitable saturation–capillary pressure relationship. For wettable soil, this can be achieved with the standard van Genuchten relation. For hydrophobic soil, this is not applicable; thus, an alternative is employed. The model is then validated against field-scale experimental observations by <span><span>Lipsius and Mooney (2006)</span></span>, which examined the impact of soil heterogeneity on infiltration profiles. The results demonstrate the model’s ability to capture complex flow dynamics in hydrophobic soils, extending the understanding of moisture transport in heterogeneous soils by explicitly modelling the spatial variability of wettability and its impact on soil hydraulic response.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107938"},"PeriodicalIF":6.2,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022685","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

Two-stage bonded-replacement breakage model and study on mechanical properties of calcareous sand considering fracture breakage and corner breakage 考虑断口破坏和角点破坏的钙质砂两阶段粘结置换破坏模型及力学性能研究

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

Computers and Geotechnics

Pub Date : 2026-01-21 DOI: 10.1016/j.compgeo.2026.107934

Hongxiang Tang , Feng Liu , Xiangji Ye , Honghua Zhao , Zongyuan Pan , Feng Zhu

Fracture breakage and corner breakage are crucial physical characteristics of calcareous sand that significantly influences its macroscopic mechanical response. However, existing breakage models based on the Discrete Element Method (DEM) struggle to simultaneously capture these two types of failure modes, and since they require newly generated particles to strictly fit within the space of the original particle after breakage, the inevitable gaps between these new particles inevitably result in volume loss within the particle assembly. This compromises the physical consistency and computational accuracy of the results. To address those limitations, this study proposes a Two-Stage Bonded-Replacement Breakage Model (TS-BRBM) based on regular triangulation to simulate both fracture breakage and corner breakage and conducts a series of triaxial compression tests of calcareous sand. The results show that TS-BRBM can effectively simulate the fracture breakage and corner breakage while maintain the mass conservation. Both fracture breakage and corner breakage have an impact on the peak strength of the test, and fracture breakage has a greater impact on the peak strength of the test than corner breakage. For same parameters, the peak strength of the traditional breakage replacement model is lower due to loss of particle volume during breakage, whereas TS-BRBM is more reasonable as it experiences no such loss. The particle shape has a certain influence on the peak strength, and tests show the peak strength of elongated angular particles is higher than that of spherical angular particles. The sedimentation angles have varying impacts on mechanical properties, with little effect in cases involving fracture breakage and corner breakage, but a more pronounced effect when breakage does not occur. This study facilitates a comprehensive consideration of both the fracture breakage and corner breakage characteristics of calcareous sand in engineering practice.

断裂断裂和角点断裂是钙质砂的重要物理特性，对钙质砂的宏观力学响应有重要影响。然而，现有的基于离散元法（DEM）的破碎模型很难同时捕捉到这两种类型的破坏模式，并且由于它们要求新生成的颗粒在破碎后严格贴合在原始颗粒的空间内，这些新颗粒之间不可避免的间隙不可避免地导致颗粒组装内的体积损失。这损害了结果的物理一致性和计算精度。针对上述局限性，本文提出了一种基于规则三角剖分的两阶段粘结-置换破坏模型（TS-BRBM），模拟裂隙破坏和角点破坏，并对钙质砂进行了一系列三轴压缩试验。结果表明，TS-BRBM能在保持质量守恒的前提下，有效地模拟断裂破坏和边角破坏。断口断裂和角点断裂对试验峰值强度均有影响，且断口断裂对试验峰值强度的影响大于角点断裂。在相同参数下，传统的破碎置换模型由于破碎过程中颗粒体积的损失，其峰值强度较低，而TS-BRBM模型由于没有颗粒体积的损失，其峰值强度更为合理。颗粒形状对峰值强度有一定的影响，试验表明，细长角状颗粒的峰值强度高于球形角状颗粒。沉降角对力学性能的影响各不相同，在断口破碎和断角破碎的情况下影响较小，而在不发生破碎的情况下影响更为显著。本研究有助于在工程实践中综合考虑钙质砂的断裂破坏和角破坏特性。

{"title":"Two-stage bonded-replacement breakage model and study on mechanical properties of calcareous sand considering fracture breakage and corner breakage","authors":"Hongxiang Tang , Feng Liu , Xiangji Ye , Honghua Zhao , Zongyuan Pan , Feng Zhu","doi":"10.1016/j.compgeo.2026.107934","DOIUrl":"10.1016/j.compgeo.2026.107934","url":null,"abstract":"<div><div>Fracture breakage and corner breakage are crucial physical characteristics of calcareous sand that significantly influences its macroscopic mechanical response. However, existing breakage models based on the Discrete Element Method (DEM) struggle to simultaneously capture these two types of failure modes, and since they require newly generated particles to strictly fit within the space of the original particle after breakage, the inevitable gaps between these new particles inevitably result in volume loss within the particle assembly. This compromises the physical consistency and computational accuracy of the results. To address those limitations, this study proposes a Two-Stage Bonded-Replacement Breakage Model (TS-BRBM) based on regular triangulation to simulate both fracture breakage and corner breakage and conducts a series of triaxial compression tests of calcareous sand. The results show that TS-BRBM can effectively simulate the fracture breakage and corner breakage while maintain the mass conservation. Both fracture breakage and corner breakage have an impact on the peak strength of the test, and fracture breakage has a greater impact on the peak strength of the test than corner breakage. For same parameters, the peak strength of the traditional breakage replacement model is lower due to loss of particle volume during breakage, whereas TS-BRBM is more reasonable as it experiences no such loss. The particle shape has a certain influence on the peak strength, and tests show the peak strength of elongated angular particles is higher than that of spherical angular particles. The sedimentation angles have varying impacts on mechanical properties, with little effect in cases involving fracture breakage and corner breakage, but a more pronounced effect when breakage does not occur. This study facilitates a comprehensive consideration of both the fracture breakage and corner breakage characteristics of calcareous sand in engineering practice.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107934"},"PeriodicalIF":6.2,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022686","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-hydro-mechanical modeling of root–soil interaction in unsaturated slopes 非饱和边坡根-土相互作用的热-水-力学模拟

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

Computers and Geotechnics

Pub Date : 2026-01-20 DOI: 10.1016/j.compgeo.2026.107931

Ehsan Badakhshan , Jean Vaunat

Vegetation effects in unsaturated soils are often simplified by adjusting strength parameters such as cohesion or friction angle, but this approach overlooks the progressive nature of root–soil interaction. To overcome this, an extended Barcelona Basic Model (BBM-VEG) is developed within the thermo-hydro-mechanical (THM) framework. The model introduces a strain-dependent reinforcement parameter (Rp_veg), linked to the root mass fraction and activation strain, which dynamically modifies soil stiffness and strength. By coupling suction-dependent elastoplasticity with Bishop’s effective stress, the yield surface evolves with both suction and root effects, capturing strain hardening and softening in vegetated soils. Unlike earlier approaches, BBM-VEG explicitly represents the mobilization and degradation of root reinforcement with strain, enabling more realistic stress–strain and strength predictions. Validation against triaxial and tensile strength tests on rooted soils confirms its ability to reproduce both initial stiffness and peak strength across varying root densities. Sensitivity analyses showed that while the peak reinforcement controls the maximum strength, the activation strain governs when reinforcement develops, with lower values causing earlier stiffening. Finally, simulations of a representative slope subjected to hydraulic (0.001, 0.002, 0.003, and 0.0005 kg/s infiltration) and thermal (30, 45, 60, and 15 °C) cycles demonstrate that vegetation promotes drier conditions by limiting infiltration and enhancing evapotranspiration. While rainfall triggers displacement, total deformation in vegetated slopes is reduced by 65–70 % compared with bare slopes.

非饱和土中植被效应通常通过调整黏聚力或摩擦角等强度参数来简化，但这种方法忽略了根土相互作用的渐进性。为了克服这个问题，在热-水-机械（THM）框架内开发了扩展的巴塞罗那基本模型（BBM-VEG）。该模型引入了一个与根质量分数和激活应变相关的应变相关加固参数（Rpveg），该参数动态修改了土壤的刚度和强度。通过将吸力依赖的弹塑性与Bishop有效应力相耦合，研究了植被土壤屈服面在吸力和根效应下的演化，捕捉了植被土壤的应变硬化和软化过程。与早期的方法不同，BBM-VEG明确地表示了根筋随应变的动员和降解，从而实现了更现实的应力-应变和强度预测。对有根土壤的三轴和拉伸强度试验的验证证实了其在不同根密度下重现初始刚度和峰值强度的能力。敏感性分析表明，峰值配筋控制最大强度，激活应变控制配筋发展，配筋值越低，越早加劲。最后，对一个具有代表性的斜坡进行了水力（0.001、0.002、0.003和0.0005 kg/s入渗）和热（30、45、60和15°C）循环的模拟，结果表明植被通过限制入渗和增强蒸散发来促进干燥条件。当降雨引起位移时，植被边坡的总变形比光秃秃的边坡减少65 - 70%。

{"title":"Thermo-hydro-mechanical modeling of root–soil interaction in unsaturated slopes","authors":"Ehsan Badakhshan , Jean Vaunat","doi":"10.1016/j.compgeo.2026.107931","DOIUrl":"10.1016/j.compgeo.2026.107931","url":null,"abstract":"<div><div>Vegetation effects in unsaturated soils are often simplified by adjusting strength parameters such as cohesion or friction angle, but this approach overlooks the progressive nature of root–soil interaction. To overcome this, an extended Barcelona Basic Model (BBM-VEG) is developed within the thermo-hydro-mechanical (THM) framework. The model introduces a strain-dependent reinforcement parameter (<em>Rp<sub>veg</sub></em>), linked to the root mass fraction and activation strain, which dynamically modifies soil stiffness and strength. By coupling suction-dependent elastoplasticity with Bishop’s effective stress, the yield surface evolves with both suction and root effects, capturing strain hardening and softening in vegetated soils. Unlike earlier approaches, BBM-VEG explicitly represents the mobilization and degradation of root reinforcement with strain, enabling more realistic stress–strain and strength predictions. Validation against triaxial and tensile strength tests on rooted soils confirms its ability to reproduce both initial stiffness and peak strength across varying root densities. Sensitivity analyses showed that while the peak reinforcement controls the maximum strength, the activation strain governs when reinforcement develops, with lower values causing earlier stiffening. Finally, simulations of a representative slope subjected to hydraulic (0.001, 0.002, 0.003, and 0.0005 kg/s infiltration) and thermal (30, 45, 60, and 15 °C) cycles demonstrate that vegetation promotes drier conditions by limiting infiltration and enhancing evapotranspiration. While rainfall triggers displacement, total deformation in vegetated slopes is reduced by 65–70 % compared with bare slopes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107931"},"PeriodicalIF":6.2,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022684","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

Image-based study of SRM seepage considering the surface roughness 考虑表面粗糙度的SRM渗流图像研究

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

Computers and Geotechnics

Pub Date : 2026-01-20 DOI: 10.1016/j.compgeo.2026.107936

Xinying Shen , Lingxi Chu , Jie Wu , Tingshuo Jia

The permeability of soil rock mixture (SRM), which are ubiquitous geomaterials, is fundamentally controlled by their meso-structure. However, the inherent complexity of the structure and the difficulty in generating conformal meshes have historically hindered the efficient, automated development of high-fidelity geometric models required for comprehensive parameter sensitivity analysis. To address this need, this paper propose a parameter sensitivity analysis approach that integrates image-based models with unfitted meshes. The methodology starts with the extraction of rock blocks from real images, followed by the parametric synthesis of digital images that accurately represent SRM meso-structures. Utilizing these constructed images within an established automated analysis workflow by the authors Shen et al. (2025), we systematically investigate how rock block content and size distribution influence seepage behavior. To the best of our knowledge, this study is the first to introduce surface roughness as a key meso-structural parameter. We employed the Savitzky–Golay filter to modulate and the fractal dimension to quantify block surface roughness, thereby quantitatively demonstrating its significant impact on permeability. This research provides valuable insights and a practical framework for automated seepage analysis that incorporates realistic SRM structures.

土石混合体是一种普遍存在的岩土材料，其渗透性从根本上受其细观结构控制。然而，结构固有的复杂性和生成保形网格的难度一直阻碍着综合参数灵敏度分析所需的高保真几何模型的高效、自动化开发。为了解决这一问题，本文提出了一种将基于图像的模型与未拟合网格相结合的参数灵敏度分析方法。该方法首先从真实图像中提取岩石块，然后对精确表示SRM细观结构的数字图像进行参数化合成。利用Shen等人（2025）在已建立的自动化分析工作流中构建的这些图像，我们系统地研究了岩块含量和尺寸分布如何影响渗流行为。据我们所知，这项研究是第一次将表面粗糙度作为一个关键的细观结构参数。我们采用Savitzky-Golay滤波进行调制，采用分形维数量化块体表面粗糙度，从而定量证明其对渗透率的显著影响。该研究为自动化渗流分析提供了有价值的见解和实用框架，并结合了实际的SRM结构。

{"title":"Image-based study of SRM seepage considering the surface roughness","authors":"Xinying Shen , Lingxi Chu , Jie Wu , Tingshuo Jia","doi":"10.1016/j.compgeo.2026.107936","DOIUrl":"10.1016/j.compgeo.2026.107936","url":null,"abstract":"<div><div>The permeability of soil rock mixture (SRM), which are ubiquitous geomaterials, is fundamentally controlled by their meso-structure. However, the inherent complexity of the structure and the difficulty in generating conformal meshes have historically hindered the efficient, automated development of high-fidelity geometric models required for comprehensive parameter sensitivity analysis. To address this need, this paper propose a parameter sensitivity analysis approach that integrates image-based models with unfitted meshes. The methodology starts with the extraction of rock blocks from real images, followed by the parametric synthesis of digital images that accurately represent SRM meso-structures. Utilizing these constructed images within an established automated analysis workflow by the authors Shen et al. (2025), we systematically investigate how rock block content and size distribution influence seepage behavior. To the best of our knowledge, this study is the first to introduce surface roughness as a key meso-structural parameter. We employed the Savitzky–Golay filter to modulate and the fractal dimension to quantify block surface roughness, thereby quantitatively demonstrating its significant impact on permeability. This research provides valuable insights and a practical framework for automated seepage analysis that incorporates realistic SRM structures.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"192 ","pages":"Article 107936"},"PeriodicalIF":6.2,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146022681","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