Computers and Geotechnics最新文献

CFD-DEM simulation of sand liquefaction with non-spherical particles and inherent anisotropic effects 含非球形颗粒及固有各向异性的砂土液化CFD-DEM模拟

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

Computers and Geotechnics

Pub Date : 2026-01-31 DOI: 10.1016/j.compgeo.2026.107929

Hongmei Gao , Wenhao Xu , Yinqiang Liu , Zhifu Shen , Xinlei Zhang , Zhihua Wang

The evaluation of sand liquefaction has long faced two major technical bottlenecks. Firstly, conventional centrifuge tests and finite element numerical simulations struggle to precisely control granular deposition anisotropy (e.g., deposition angle) and accurately characterize the interactions between fluid and non-spherical particles. Secondly, due to insufficient control of dynamic similarity, the Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) coupled methods encounter significant computational efficiency challenges in the large-scale site simulations. To address these issues, this study innovatively proposes an improved CFD-DEM coupling framework, achieving methodological integration and parameter optimization in two key aspects: (1) incorporation of a non-spherical particle model to accurately characterize the directional effects of particle shape on fluid resistance; and (2) through refined adjustment of key parameter matching relationships including fluid viscosity, coupling forces, and particle Reynolds number, enabling equivalent simulation of high-gravity models while strictly maintaining physical consistency, thereby significantly improving computational efficiency. Within this framework, periodic boundary conditions were effectively employed to eliminate rigid boundary interference and achieve high-precision control of initial fabric anisotropy. Using this methodological system, the study successfully reproduced the liquefaction response differences in the sand layers with three deposition angles (0°, 45°, and 90°). It reveals that deposition angle exerts significant control on the soil liquefaction resistance: horizontally deposited (0°) sand layers demonstrate the optimal anti-liquefaction capacity due to their stable force chain network structure, while vertically deposited (90°) sand layers exhibit the highest liquefaction susceptibility owing to rapid particle suspension (suspension coefficient β_t→1.0) and pronounced pore compression effects. The findings offer some micro-mechanistic insights for seismic liquefaction risk assessment in the sites with natural deposition anisotropy.

砂土液化评价长期以来面临两大技术瓶颈。首先，传统的离心机试验和有限元数值模拟难以精确控制颗粒沉积的各向异性（如沉积角度）和准确表征流体与非球形颗粒之间的相互作用。其次，计算流体动力学-离散元法（CFD-DEM）耦合方法由于对动力相似度控制不足，在大尺度现场模拟中计算效率面临较大挑战。为了解决这些问题，本研究创新性地提出了一种改进的CFD-DEM耦合框架，在两个关键方面实现了方法集成和参数优化：(1)结合非球形颗粒模型，准确表征颗粒形状对流体阻力的定向影响；(2)通过对流体粘度、耦合力、粒子雷诺数等关键参数匹配关系的精细调整，在严格保持物理一致性的前提下，实现了高重力模型的等效模拟，显著提高了计算效率。在此框架下，利用周期边界条件有效地消除了刚性边界干扰，实现了织物初始各向异性的高精度控制。利用该方法系统，该研究成功地再现了三种沉积角度（0°、45°和90°）下砂层液化响应的差异。结果表明，沉积角度对土壤液化阻力具有显著的控制作用，水平沉积（0°）砂层具有稳定的力链网络结构，具有最佳的抗液化能力，而垂直沉积（90°）砂层由于颗粒快速悬浮（悬浮系数βt→1.0）和明显的孔隙压缩效应，具有最高的液化敏感性。研究结果为自然沉积各向异性场地的地震液化风险评估提供了一些微观机制见解。

{"title":"CFD-DEM simulation of sand liquefaction with non-spherical particles and inherent anisotropic effects","authors":"Hongmei Gao , Wenhao Xu , Yinqiang Liu , Zhifu Shen , Xinlei Zhang , Zhihua Wang","doi":"10.1016/j.compgeo.2026.107929","DOIUrl":"10.1016/j.compgeo.2026.107929","url":null,"abstract":"<div><div>The evaluation of sand liquefaction has long faced two major technical bottlenecks. Firstly, conventional centrifuge tests and finite element numerical simulations struggle to precisely control granular deposition anisotropy (e.g., deposition angle) and accurately characterize the interactions between fluid and non-spherical particles. Secondly, due to insufficient control of dynamic similarity, the Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) coupled methods encounter significant computational efficiency challenges in the large-scale site simulations. To address these issues, this study innovatively proposes an improved CFD-DEM coupling framework, achieving methodological integration and parameter optimization in two key aspects: (1) incorporation of a non-spherical particle model to accurately characterize the directional effects of particle shape on fluid resistance; and (2) through refined adjustment of key parameter matching relationships including fluid viscosity, coupling forces, and particle Reynolds number, enabling equivalent simulation of high-gravity models while strictly maintaining physical consistency, thereby significantly improving computational efficiency. Within this framework, periodic boundary conditions were effectively employed to eliminate rigid boundary interference and achieve high-precision control of initial fabric anisotropy. Using this methodological system, the study successfully reproduced the liquefaction response differences in the sand layers with three deposition angles (0°, 45°, and 90°). It reveals that deposition angle exerts significant control on the soil liquefaction resistance: horizontally deposited (0°) sand layers demonstrate the optimal anti-liquefaction capacity due to their stable force chain network structure, while vertically deposited (90°) sand layers exhibit the highest liquefaction susceptibility owing to rapid particle suspension (suspension coefficient <em>β<sub>t</sub></em>→1.0) and pronounced pore compression effects. The findings offer some micro-mechanistic insights for seismic liquefaction risk assessment in the sites with natural deposition anisotropy.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107929"},"PeriodicalIF":6.2,"publicationDate":"2026-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080664","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

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

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

Computers and Geotechnics

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

Milad Zaheri , Pierpaolo Oreste , Masoud Ranjbarnia , Elham Mahmoudi

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

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

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

引用次数: 0

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

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

Computers and Geotechnics

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

Keitaro Hoshi, Shotaro Yamada, Yuta Abe, Takashi Kyoya

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

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

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

引用次数: 0

TKLE-BPINN: A Bayesian physics-informed inversion framework for high-dimensional parameter identification in geotechnical subsurface systems TKLE-BPINN：用于岩土地下系统高维参数识别的贝叶斯物理信息反演框架

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

Computers and Geotechnics

Pub Date : 2026-01-30 DOI: 10.1016/j.compgeo.2026.107957

Zhenjie Tang, Li He

Characterizing spatially variable hydraulic properties in geotechnical subsurface systems is a high-dimensional and ill-posed inverse problem, particularly in unsaturated flow conditions. We introduce TKLE-BPINN, a novel Bayesian framework that integrates Bayesian physics-informed neural networks (B-PINNs) with truncated Karhunen–Loève expansion (KLE) to address these challenges efficiently. By representing unknown fields with latent KLE coefficients and employing Ensemble Kalman Inversion (EKI) for posterior inference, TKLE-BPINN achieves accurate parameter estimation and robust uncertainty quantification. We validate the framework on three subsurface flow and transport test cases: (i) a linear diffusion–reaction problem with unknown diffusivity, (ii) a nonlinear diffusion–reaction problem with unknown reaction coefficients, and (iii) a multi-parameter inversion of diffusivity and contaminant source locations, where TKLE-BPINN consistently outperforms standard B-PINNs in accuracy and stability. We further apply TKLE-BPINN to a realistic hydrogeotechnical scenario, inverting saturated hydraulic conductivity in nonlinear unsaturated flow for heterogeneous soils, yielding reliable parameter estimates and uncertainty quantification essential for subsurface flow modeling. These results indicate that TKLE-BPINN provides a flexible and robust alternative for geotechnical subsurface characterization.

岩土地下系统中空间可变水力特性的表征是一个高维的不适定反问题，特别是在非饱和流动条件下。我们介绍了TKLE-BPINN，这是一种新的贝叶斯框架，它将贝叶斯物理信息神经网络（b - pinn）与截断karhunen - lo扩展（KLE）相结合，以有效地解决这些挑战。TKLE-BPINN通过用潜在KLE系数表示未知域，并采用集合卡尔曼反演（EKI）进行后验推理，实现了精确的参数估计和鲁棒的不确定性量化。我们在三个地下流动和运输测试案例中验证了该框架：(i)具有未知扩散系数的线性扩散反应问题，（ii）具有未知反应系数的非线性扩散反应问题，以及（iii）具有扩散系数和污染源位置的多参数反演，其中TKLE-BPINN在准确性和稳定性方面始终优于标准b - pinn。我们进一步将TKLE-BPINN应用于现实的水文岩土工程场景，反演非均质土壤非线性非饱和流中的饱和水力传导性，为地下流建模提供可靠的参数估计和不确定性量化。这些结果表明，TKLE-BPINN为岩土工程地下表征提供了一种灵活而可靠的替代方法。

{"title":"TKLE-BPINN: A Bayesian physics-informed inversion framework for high-dimensional parameter identification in geotechnical subsurface systems","authors":"Zhenjie Tang, Li He","doi":"10.1016/j.compgeo.2026.107957","DOIUrl":"10.1016/j.compgeo.2026.107957","url":null,"abstract":"<div><div>Characterizing spatially variable hydraulic properties in geotechnical subsurface systems is a high-dimensional and ill-posed inverse problem, particularly in unsaturated flow conditions. We introduce TKLE-BPINN, a novel Bayesian framework that integrates Bayesian physics-informed neural networks (B-PINNs) with truncated Karhunen–Loève expansion (KLE) to address these challenges efficiently. By representing unknown fields with latent KLE coefficients and employing Ensemble Kalman Inversion (EKI) for posterior inference, TKLE-BPINN achieves accurate parameter estimation and robust uncertainty quantification. We validate the framework on three subsurface flow and transport test cases: (i) a linear diffusion–reaction problem with unknown diffusivity, (ii) a nonlinear diffusion–reaction problem with unknown reaction coefficients, and (iii) a multi-parameter inversion of diffusivity and contaminant source locations, where TKLE-BPINN consistently outperforms standard B-PINNs in accuracy and stability. We further apply TKLE-BPINN to a realistic hydrogeotechnical scenario, inverting saturated hydraulic conductivity in nonlinear unsaturated flow for heterogeneous soils, yielding reliable parameter estimates and uncertainty quantification essential for subsurface flow modeling. These results indicate that TKLE-BPINN provides a flexible and robust alternative for geotechnical subsurface characterization.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107957"},"PeriodicalIF":6.2,"publicationDate":"2026-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080662","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

CFD-DEM investigation on particle clogging in porous media at different scales: From single pore to digital rock 不同尺度多孔介质中颗粒堵塞的CFD-DEM研究：从单孔到数字岩石

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

Computers and Geotechnics

Pub Date : 2026-01-30 DOI: 10.1016/j.compgeo.2026.107951

Yuanping Li , Ruyang Yu , Xiaolong Yin , Huanquan Pan , Bin Gong , Lifeng Chen , Jingwei Huang

Particle clogging in porous media is a critical phenomenon with significant implications for geotechnical engineering, subsurface flow, and underground carbon storage. However, existing studies are limited to single-scale analysis, lacking multi-scale insights into the connections between simplified models and real porous media. To address these limitations, this study conducts a multi-scale investigation spanning single pores, homogeneous porous media, and digital rocks to elucidate clogging mechanisms. In the single-pore, three clogging regimes were identified: non-clogging, unstable clogging, and stable clogging. The critical thresholds of the three regimes were determined. In homogeneous porous medium, the clogging process could undergo three stages. Firstly, particles only blocked the dominant flow paths. Secondly, the new preferential pathways experience an increase in velocity and particle flux, leading to secondary clogging. Thirdly, the primary and secondary paths were mostly blocked, and a network of blockages was formed. The three clogging modes were summarized: selective channel clogging, localized bridging, and network-scale blockage, which are controlled by particle-to-throat size ratio, flow velocity and particle concentration. Network-scale blockages result in the most significant decline in permeability, while selective channel clogging leads to the least. In digital rocks, clogging exhibited the similar clogging behaviors observed in homogeneous porous media, but it exhibited distinct permeability loss due to the heterogeneity of realistic pore structures. This multi-scale study quantifies the regulatory effects of key parameters on clogging across scales, and identifies scale-specific patterns and mechanisms.

多孔介质中的颗粒堵塞是岩土工程、地下流动和地下碳储存的重要现象。然而，现有的研究仅限于单尺度分析，缺乏对简化模型与真实多孔介质之间联系的多尺度洞察。为了解决这些局限性，本研究开展了跨单孔、均质多孔介质和数字岩石的多尺度研究，以阐明堵塞机制。在单孔中，确定了三种堵塞状态：非堵塞、不稳定堵塞和稳定堵塞。确定了这三种制度的临界阈值。在均匀多孔介质中，堵塞过程可分为三个阶段。首先，颗粒只阻塞了主要的流动路径。其次，新的优先通道经历速度和粒子通量的增加，导致二次堵塞。三是主次路径多被阻断，形成阻塞网络。总结了三种堵塞模式：选择性通道堵塞、局部桥接和网络级堵塞，这些堵塞模式受颗粒与喉道尺寸比、流速和颗粒浓度的控制。网络规模堵塞导致渗透率下降最显著，而选择性通道堵塞导致的渗透率下降最小。在数字岩石中，堵塞表现出与均匀多孔介质相似的堵塞行为，但由于真实孔隙结构的非均质性，其渗透率损失明显。这项多尺度研究量化了关键参数对跨尺度堵塞的调节作用，并确定了尺度特定的模式和机制。

{"title":"CFD-DEM investigation on particle clogging in porous media at different scales: From single pore to digital rock","authors":"Yuanping Li , Ruyang Yu , Xiaolong Yin , Huanquan Pan , Bin Gong , Lifeng Chen , Jingwei Huang","doi":"10.1016/j.compgeo.2026.107951","DOIUrl":"10.1016/j.compgeo.2026.107951","url":null,"abstract":"<div><div>Particle clogging in porous media is a critical phenomenon with significant implications for geotechnical engineering, subsurface flow, and underground carbon storage. However, existing studies are limited to single-scale analysis, lacking multi-scale insights into the connections between simplified models and real porous media. To address these limitations, this study conducts a multi-scale investigation spanning single pores, homogeneous porous media, and digital rocks to elucidate clogging mechanisms. In the single-pore, three clogging regimes were identified: non-clogging, unstable clogging, and stable clogging. The critical thresholds of the three regimes were determined. In homogeneous porous medium, the clogging process could undergo three stages. Firstly, particles only blocked the dominant flow paths. Secondly, the new preferential pathways experience an increase in velocity and particle flux, leading to secondary clogging. Thirdly, the primary and secondary paths were mostly blocked, and a network of blockages was formed. The three clogging modes were summarized: selective channel clogging, localized bridging, and network-scale blockage, which are controlled by particle-to-throat size ratio, flow velocity and particle concentration. Network-scale blockages result in the most significant decline in permeability, while selective channel clogging leads to the least. In digital rocks, clogging exhibited the similar clogging behaviors observed in homogeneous porous media, but it exhibited distinct permeability loss due to the heterogeneity of realistic pore structures. This multi-scale study quantifies the regulatory effects of key parameters on clogging across scales, and identifies scale-specific patterns and mechanisms.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107951"},"PeriodicalIF":6.2,"publicationDate":"2026-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080722","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

Physics-informed neural networks (PINNs) for dynamic pile-soil interaction problems 桩-土动力相互作用问题的物理信息神经网络

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

Computers and Geotechnics

Pub Date : 2026-01-29 DOI: 10.1016/j.compgeo.2026.107941

Juntao Wu , Weikai Zhao , Kuihua Wang , M. Hesham El Naggar

Traditional analytical and numerical methods for dynamic pile-soil interaction (PSI) have limitations in solving the problem’s partial differential equations (PDEs). To address some of these limitations, a transfer learning enhanced progressive multi-physics-informed neural networks (TLP-mPINNs) framework is developed in this paper. The improved PINNs model utilizes customized independent deep neural networks (DNNs) for pile and soil to adapt to different material domains. First, a pre-trained pile vibration network is transferred to a more sophisticated pile-soil coupled vibration network to accelerate the training process. Second, the parameters of pile and soil DNNs are progressively trained through alternating learning to improve the model stability. The predictions of the TLP-mPINNs are verified using existing solutions of the pile-half space soil coupled vibration system. The validated model is then employed to investigate the dynamic PSI problem under different cases. The results demonstrate that the improved PINN framework can provide a new solution for the dynamic PSI problems, and highlight its potential as an optimization methodology for applying PINN to complex multi-media coupled vibration problems.

传统的桩土动力相互作用的解析和数值方法在求解桩土动力相互作用的偏微分方程时存在局限性。为了解决这些限制，本文开发了一种迁移学习增强的渐进式多物理信息神经网络（TLP-mPINNs）框架。改进的PINNs模型利用自定义的桩土独立深度神经网络（dnn）来适应不同的材料域。首先，将预训练好的桩振动网络转化为更复杂的桩-土耦合振动网络，加快训练过程；其次，通过交替学习逐步训练桩土深度神经网络参数，提高模型稳定性；利用已有的桩-半空间土耦合振动系统解，验证了TLP-mPINNs的预测结果。然后利用验证的模型对不同情况下的动态PSI问题进行了研究。结果表明，改进的PINN框架可以为动态PSI问题提供新的解决方案，并突出了其作为将PINN应用于复杂多媒体耦合振动问题的优化方法的潜力。

{"title":"Physics-informed neural networks (PINNs) for dynamic pile-soil interaction problems","authors":"Juntao Wu , Weikai Zhao , Kuihua Wang , M. Hesham El Naggar","doi":"10.1016/j.compgeo.2026.107941","DOIUrl":"10.1016/j.compgeo.2026.107941","url":null,"abstract":"<div><div>Traditional analytical and numerical methods for dynamic pile-soil interaction (PSI) have limitations in solving the problem’s partial differential equations (PDEs). To address some of these limitations, a transfer learning enhanced progressive multi-physics-informed neural networks (TLP-mPINNs) framework is developed in this paper. The improved PINNs model utilizes customized independent deep neural networks (DNNs) for pile and soil to adapt to different material domains. First, a pre-trained pile vibration network is transferred to a more sophisticated pile-soil coupled vibration network to accelerate the training process. Second, the parameters of pile and soil DNNs are progressively trained through alternating learning to improve the model stability. The predictions of the TLP-mPINNs are verified using existing solutions of the pile-half space soil coupled vibration system. The validated model is then employed to investigate the dynamic PSI problem under different cases. The results demonstrate that the improved PINN framework can provide a new solution for the dynamic PSI problems, and highlight its potential as an optimization methodology for applying PINN to complex multi-media coupled vibration problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107941"},"PeriodicalIF":6.2,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080665","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

Stress integration algorithm of the small-strain Unified Hardening model for soils based on multistage Homotopy continuation method 基于多阶段同伦延拓法的土体小应变统一硬化模型应力积分算法

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

Computers and Geotechnics

Pub Date : 2026-01-29 DOI: 10.1016/j.compgeo.2025.107860

Yangping Yao , Yu Tian , Zijun Yao , Dechun Lu , Xiuli Du

The small-strain Unified Hardening (SSUH) model considers the high initial stiffness and rapid stiffness degradation of soils during the loading process. When its stress integration is implemented by an implicit algorithm in the finite element analysis, a key issue lies in how to solve the highly nonlinear constitutive equations. This paper proposes multistage Homotopy continuation method (MHCM), which progressively optimizes the initial guess by several stages of Homotopic deformation, to ensure the successful solution of the constitutive equations using Newton-Raphson iteration. An adaptive Homotopic deformation rate is introduced to prevent the subsequent stage of Homotopic deformation from repeating the previous failure. Compared with the original single-stage Homotopy continuation method, MHCM improves the convergence and efficiency without compromising the accuracy. Based on the SSUH model and the proposed stress integration algorithm, the ground displacement induced by the excavation of Crossrail tunnels can be reasonably predicted.

小应变统一硬化（SSUH）模型考虑了土体在加载过程中的高初始刚度和快速刚度退化。在有限元分析中采用隐式算法进行应力积分时，关键问题在于如何求解高度非线性的本构方程。本文提出了多阶段同伦延拓法（MHCM），通过多阶段同伦变形逐步优化初始猜测，以保证牛顿-拉夫森迭代求解本构方程的成功。引入自适应同伦变形速率，防止同伦变形后续阶段重复前一阶段的破坏。与原来的单阶段同伦延拓方法相比，MHCM在不影响精度的前提下提高了收敛性和效率。基于SSUH模型和提出的应力积分算法，可以合理预测横贯铁路隧道开挖引起的地面位移。

{"title":"Stress integration algorithm of the small-strain Unified Hardening model for soils based on multistage Homotopy continuation method","authors":"Yangping Yao , Yu Tian , Zijun Yao , Dechun Lu , Xiuli Du","doi":"10.1016/j.compgeo.2025.107860","DOIUrl":"10.1016/j.compgeo.2025.107860","url":null,"abstract":"<div><div>The small-strain Unified Hardening (SSUH) model considers the high initial stiffness and rapid stiffness degradation of soils during the loading process. When its stress integration is implemented by an implicit algorithm in the finite element analysis, a key issue lies in how to solve the highly nonlinear constitutive equations. This paper proposes multistage Homotopy continuation method (MHCM), which progressively optimizes the initial guess by several stages of Homotopic deformation, to ensure the successful solution of the constitutive equations using Newton-Raphson iteration. An adaptive Homotopic deformation rate is introduced to prevent the subsequent stage of Homotopic deformation from repeating the previous failure. Compared with the original single-stage Homotopy continuation method, MHCM improves the convergence and efficiency without compromising the accuracy. Based on the SSUH model and the proposed stress integration algorithm, the ground displacement induced by the excavation of Crossrail tunnels can be reasonably predicted.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107860"},"PeriodicalIF":6.2,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080671","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 modified nonlinear seepage-stress coupling model of heterogeneous pore-fracture dual medium for engineering geomaterials 工程岩土材料非均质孔裂双介质非线性渗流-应力耦合修正模型

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

Computers and Geotechnics

Pub Date : 2026-01-29 DOI: 10.1016/j.compgeo.2026.107939

Wentao Yang, Shaobo Jin, Qiuhua Rao, Wei Yi, Zelin Liu

The seepage-stress coupling effect is critical to the stability of engineering structures, but existing models have certain limitations: the equivalent continuum model neglects seepage differences between pores and fractures; the discrete fracture network model overlooks pore permeability and involves complex calculations; and the dual medium model fails to account for non-Darcy seepage effects and the spatiotemporal evolution of physical parameters. Consequently, solving the seepage-stress coupling problem of spatially heterogeneous pore-fracture dual media remains challenging. In this study, we incorporate two key effects into the modeling framework: the non-Darcy seepage and the nonlinear seepage-stress coupling (including the spatiotemporal evolution of seepage parameters, solid deformation parameters, and effective stress coefficients). The model is developed via self-programming in COMSOL Multiphysics, and its reliability is validated through lab experiments. Results show that both effects suppress the dissipation of fluid pressure, the consolidation rate, and the final settlement of the soil. The degree of such suppression is governed by fracture permeability, pore permeability, and their ratio. Additionally, the practicality of this model is verified through a 3D engineering case. The model can be further extended to a multi-scale fracture network model (explicitly characterizing main fractures and implicitly describing randomly distributed secondary fractures and pores) with consideration of permeability anisotropy.

渗流-应力耦合效应对工程结构的稳定性至关重要，但现有模型存在一定的局限性：等效连续介质模型忽略了孔隙与裂隙之间的渗流差异；离散裂缝网络模型忽略了孔隙渗透率，计算复杂；双介质模型不能考虑非达西渗流效应和物理参数的时空演化。因此，解决空间非均质孔隙-破裂双重介质的渗流-应力耦合问题仍然是一个挑战。在本研究中，我们将非达西渗流和非线性渗流-应力耦合（包括渗流参数、固体变形参数和有效应力系数的时空演变）两个关键效应纳入建模框架。在COMSOL Multiphysics中通过自编程开发了该模型，并通过实验验证了模型的可靠性。结果表明，这两种作用都抑制了流体压力的消散、固结速率和土体的最终沉降。这种抑制程度受裂缝渗透率、孔隙渗透率及其比值的影响。并通过一个三维工程实例验证了该模型的实用性。该模型可进一步扩展为考虑渗透率各向异性的多尺度裂缝网络模型（显式表征主裂缝，隐式描述随机分布的次裂缝和孔隙）。

{"title":"A modified nonlinear seepage-stress coupling model of heterogeneous pore-fracture dual medium for engineering geomaterials","authors":"Wentao Yang, Shaobo Jin, Qiuhua Rao, Wei Yi, Zelin Liu","doi":"10.1016/j.compgeo.2026.107939","DOIUrl":"10.1016/j.compgeo.2026.107939","url":null,"abstract":"<div><div>The seepage-stress coupling effect is critical to the stability of engineering structures, but existing models have certain limitations: the equivalent continuum model neglects seepage differences between pores and fractures; the discrete fracture network model overlooks pore permeability and involves complex calculations; and the dual medium model fails to account for non-Darcy seepage effects and the spatiotemporal evolution of physical parameters. Consequently, solving the seepage-stress coupling problem of spatially heterogeneous pore-fracture dual media remains challenging. In this study, we incorporate two key effects into the modeling framework: the non-Darcy seepage and the nonlinear seepage-stress coupling (including the spatiotemporal evolution of seepage parameters, solid deformation parameters, and effective stress coefficients). The model is developed via self-programming in COMSOL Multiphysics, and its reliability is validated through lab experiments. Results show that both effects suppress the dissipation of fluid pressure, the consolidation rate, and the final settlement of the soil. The degree of such suppression is governed by fracture permeability, pore permeability, and their ratio. Additionally, the practicality of this model is verified through a 3D engineering case. The model can be further extended to a multi-scale fracture network model (explicitly characterizing main fractures and implicitly describing randomly distributed secondary fractures and pores) with consideration of permeability anisotropy.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107939"},"PeriodicalIF":6.2,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080661","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

Multiphysics modelling of millimetre-wave ablation of geological materials 地质物质毫米波消融的多物理场模拟

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

Computers and Geotechnics

Pub Date : 2026-01-29 DOI: 10.1016/j.compgeo.2026.107946

Nandan Gokhale , Candace Gilet , Franck Monmont , Nikos Nikiforakis

This study presents a multiphysics model and the corresponding algorithm for the numerical simulation of millimetre-wave ablation of geological materials. Simulations of this process are particularly challenging due to the low thermal conductivity and limited melt mobility of rock, as well as the high latent heat of evaporation and melt-to-vapour density ratio (exceeding 2500:1), which induce strong velocity divergence effects at the melt–vapour interface. Specifically, the model introduced in this work addresses the low-intensity, volumetric heating of rocks, taking place over significantly larger spatial and temporal scales (approximately 300 and 100 times greater, respectively) compared to similar processes encountered in laser drilling of metals. A comprehensive description is provided of the mathematical formulation, beam model, and numerical algorithm employed to perform robust, fully three-dimensional simulations lasting up to 25 s. These demanding computations are made tractable through MPI parallelisation and hierarchical Adaptive Mesh Refinement. Following the presentation of the model and solution methodology, an experiment is introduced which was devised specifically to produce data for model calibration and validation purposes. Comparison of the numerical results against the experimental reference test case and a parametric study of varying beam intensity demonstrate that the model can make good predictions of the penetration depth and the volume of material removed. The bounds and limitations of the model are also explored through the parametric study, with desirable future extensions identified to improve the accuracy of the model predictions.

提出了地质物质毫米波烧蚀数值模拟的多物理场模型和相应的算法。由于岩石的低导热系数和有限的熔体迁移率，以及蒸发潜热和熔体-蒸汽密度比（超过2500:1）的高潜热，在熔体-蒸汽界面上引起强烈的速度发散效应，因此模拟这一过程特别具有挑战性。具体来说，本工作中引入的模型解决了岩石的低强度体积加热问题，与金属激光钻孔中遇到的类似过程相比，这种加热发生在更大的空间和时间尺度上（分别约为300倍和100倍）。一个全面的描述提供了数学公式，梁模型，和数值算法用于执行鲁棒，全三维模拟持续长达25秒。通过MPI并行化和分层自适应网格细化，这些苛刻的计算变得易于处理。在模型和解决方法的介绍之后，介绍了一个实验，该实验是专门为模型校准和验证目的而设计的。将数值计算结果与实验参考测试案例进行比较，并对不同光束强度的参数化研究表明，该模型能较好地预测穿透深度和去除材料的体积。通过参数化研究还探讨了模型的界限和局限性，并确定了理想的未来扩展，以提高模型预测的准确性。

{"title":"Multiphysics modelling of millimetre-wave ablation of geological materials","authors":"Nandan Gokhale , Candace Gilet , Franck Monmont , Nikos Nikiforakis","doi":"10.1016/j.compgeo.2026.107946","DOIUrl":"10.1016/j.compgeo.2026.107946","url":null,"abstract":"<div><div>This study presents a multiphysics model and the corresponding algorithm for the numerical simulation of millimetre-wave ablation of geological materials. Simulations of this process are particularly challenging due to the low thermal conductivity and limited melt mobility of rock, as well as the high latent heat of evaporation and melt-to-vapour density ratio (exceeding 2500:1), which induce strong velocity divergence effects at the melt–vapour interface. Specifically, the model introduced in this work addresses the low-intensity, volumetric heating of rocks, taking place over significantly larger spatial and temporal scales (approximately 300 and 100 times greater, respectively) compared to similar processes encountered in laser drilling of metals. A comprehensive description is provided of the mathematical formulation, beam model, and numerical algorithm employed to perform robust, fully three-dimensional simulations lasting up to 25 s. These demanding computations are made tractable through MPI parallelisation and hierarchical Adaptive Mesh Refinement. Following the presentation of the model and solution methodology, an experiment is introduced which was devised specifically to produce data for model calibration and validation purposes. Comparison of the numerical results against the experimental reference test case and a parametric study of varying beam intensity demonstrate that the model can make good predictions of the penetration depth and the volume of material removed. The bounds and limitations of the model are also explored through the parametric study, with desirable future extensions identified to improve the accuracy of the model predictions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107946"},"PeriodicalIF":6.2,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080668","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 seepage–deformation analysis of the dynamics of embankment with elastoplasticity based on the full formulation 基于全公式的弹塑性路堤渗流-变形耦合动力学分析

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

Computers and Geotechnics

Pub Date : 2026-01-28 DOI: 10.1016/j.compgeo.2026.107950

Jiawei Xu , Ryosuke Uzuoka , Kyohei Ueda , Yoshikazu Tanaka

The dynamics of embankment considering elastoplasticity is investigated using the coupled seepage–deformation finite element analysis with the full Biot formation in the u–v–p format, where the solid displacement, relative fluid velocity respect to solid, and pore fluid pressure are taken as the primary variables. The seismic response of embankment is first evaluated using the centrifuge experiment, based on which the typical seepage and deformation characteristics of embankment are investigated. The validation against the centrifuge experimental result demonstrates the capability of the coupled finite element analysis using the full formulation and elastoplasticity to predict embankment responses during seismic loading such as the horizontal acceleration, surface settlement, and pore water pressure, thus proving a robust tool to investigate the porous media dynamics. Regarding the dynamics of embankment during seismic loading with various combinations of soil permeability and loading frequency, the predicated embankment responses such as solid acceleration, pore fluid pressure, and soil deformation based on the full analysis tend to show more significant difference in comparison with those predicted by the simplified analysis that neglects the relative fluid acceleration respect to solid acceleration when the soil permeability or loading frequency increases to a high level. Based on the numerical simulation with elastoplasticity, the difference in various dynamic responses of the embankment especially the soil deformation using the full and simplified analysis approaches can mainly be divided into two distinct zones in the

k / f

–

f^{2}

space (

k

and

f

are the permeability and frequency ratios), where the significant difference is found in scenarios with

k / f

or

f

generally larger than 100 m or 100 Hz.

采用u-v-p格式的渗流-变形耦合有限元分析方法，以固体位移、相对于固体的流体相对速度和孔隙流体压力为主要变量，对考虑弹塑性的路堤动力进行了研究。首先利用离心试验对路堤的地震反应进行了评价，在此基础上研究了路堤的典型渗流和变形特征。通过对离心机实验结果的验证，证明了使用完整公式和弹塑性的耦合有限元分析能够预测地震荷载下路堤的反应，如水平加速度、地表沉降和孔隙水压力，从而证明了一个研究多孔介质动力学的强大工具。针对不同土渗透性和加载频率组合下的地震荷载作用下路堤的动力特性，预测路堤的固体加速度、孔隙流体压力、当土体渗透性或加载频率增加到较高水平时，基于完整分析的土体变形与忽略流体相对加速度与固体相对加速度的简化分析预测的土体变形差异更显著。基于弹塑性数值模拟，在k/f - f2空间（k和f分别为渗透比和频率比）上，路堤各种动力响应特别是土体变形的差异主要可以划分为两个不同的区域，其中k/f - f一般大于100 m或100 Hz的情况下差异显著。

{"title":"Coupled seepage–deformation analysis of the dynamics of embankment with elastoplasticity based on the full formulation","authors":"Jiawei Xu , Ryosuke Uzuoka , Kyohei Ueda , Yoshikazu Tanaka","doi":"10.1016/j.compgeo.2026.107950","DOIUrl":"10.1016/j.compgeo.2026.107950","url":null,"abstract":"<div><div>The dynamics of embankment considering elastoplasticity is investigated using the coupled seepage–deformation finite element analysis with the full Biot formation in the <strong><em>u</em></strong>–<strong><em>v</em></strong>–<em>p</em> format, where the solid displacement, relative fluid velocity respect to solid, and pore fluid pressure are taken as the primary variables. The seismic response of embankment is first evaluated using the centrifuge experiment, based on which the typical seepage and deformation characteristics of embankment are investigated. The validation against the centrifuge experimental result demonstrates the capability of the coupled finite element analysis using the full formulation and elastoplasticity to predict embankment responses during seismic loading such as the horizontal acceleration, surface settlement, and pore water pressure, thus proving a robust tool to investigate the porous media dynamics. Regarding the dynamics of embankment during seismic loading with various combinations of soil permeability and loading frequency, the predicated embankment responses such as solid acceleration, pore fluid pressure, and soil deformation based on the full analysis tend to show more significant difference in comparison with those predicted by the simplified analysis that neglects the relative fluid acceleration respect to solid acceleration when the soil permeability or loading frequency increases to a high level. Based on the numerical simulation with elastoplasticity, the difference in various dynamic responses of the embankment especially the soil deformation using the full and simplified analysis approaches can mainly be divided into two distinct zones in the <span><math><mrow><mi>k</mi><mo>/</mo><mi>f</mi></mrow></math></span>–<span><math><msup><mrow><mi>f</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span> space (<span><math><mi>k</mi></math></span> and <span><math><mi>f</mi></math></span> are the permeability and frequency ratios), where the significant difference is found in scenarios with <span><math><mrow><mi>k</mi><mo>/</mo><mi>f</mi></mrow></math></span> or <span><math><mi>f</mi></math></span> generally larger than 100 m or 100 Hz.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107950"},"PeriodicalIF":6.2,"publicationDate":"2026-01-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146081268","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