Computers and Geotechnics最新文献_第4页

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

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

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

Computers and Geotechnics

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

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

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

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

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

引用次数: 0

On the integration of an ACST-based bounding surface model 基于acst的边界面模型的集成

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

Computers and Geotechnics

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

Srinivas Vivek Bokkisa , Jorge Macedo , Pedro Arduino

Anisotropic critical state theory (ACST) provides a framework for incorporating fabric effects in constitutive models. However, most previous efforts have focused on constitutive aspects with comparatively limited attention to numerical implementations. This study presents a comprehensive assessment of explicit and implicit implementations of the ACST-based bounding surface model, SANISAND-F. Assessments are conducted in terms of stability, accuracy, computational efficiency, and both local and global performance.

In the implicit implementation, the critical importance of accurate gradient calculations is highlighted, introducing a verification procedure that enables quadratic convergence. The explicit and implicit implementations exhibit stability, producing smooth and bounded responses across a wide range of strain increments and numerical tolerances. However, their accuracy differs significantly. The implicit implementation is sensitive to the initial loading state, strain increment, and loading direction, showing minor dependence on the solver tolerance. In contrast, the explicit implementation is influenced by both strain increment and substepping tolerance, and at practical tolerance and strain increment levels, it often outperforms the implicit scheme in accuracy. Regarding efficiency, the explicit implementation proves more efficient at the local integration level. However, at the global level, the implicit implementation with the consistent tangent exhibits a faster rate of convergence in global equilibrium iterations. Nonetheless, the overall computational cost at the global level is not definitive when comparing explicit and implicit schemes; it varies with simulations and loading-specific parameters, as demonstrated through the included boundary-value problems.

各向异性临界状态理论（ACST）为将织物效应纳入本构模型提供了一个框架。然而，大多数以前的努力都集中在本构方面，对数值实现的关注相对有限。本研究对基于acst的边界面模型SANISAND-F的显式和隐式实现进行了全面评估。评估是根据稳定性、准确性、计算效率以及本地和全局性能进行的。在隐式实现中，强调了精确梯度计算的关键重要性，引入了一个验证程序，使二次收敛成为可能。显式和隐式实现表现出稳定性，在广泛的应变增量和数值公差范围内产生平滑和有界的响应。然而，它们的准确性差别很大。隐式实现对初始加载状态、应变增量和加载方向敏感，对求解器公差的依赖较小。相比之下，显式实现受应变增量和步进公差的影响，在实际公差和应变增量水平上，显式实现的精度往往优于隐式方案。关于效率，显式实现证明在本地集成级别更有效。然而，在全局水平上，具有一致切线的隐式实现在全局平衡迭代中表现出更快的收敛速度。尽管如此，当比较显式和隐式方案时，全球一级的总体计算成本并不确定；它随模拟和加载特定参数的变化而变化，正如通过所包含的边值问题所证明的那样。

{"title":"On the integration of an ACST-based bounding surface model","authors":"Srinivas Vivek Bokkisa , Jorge Macedo , Pedro Arduino","doi":"10.1016/j.compgeo.2026.107911","DOIUrl":"10.1016/j.compgeo.2026.107911","url":null,"abstract":"<div><div>Anisotropic critical state theory (ACST) provides a framework for incorporating fabric effects in constitutive models. However, most previous efforts have focused on constitutive aspects with comparatively limited attention to numerical implementations. This study presents a comprehensive assessment of explicit and implicit implementations of the ACST-based bounding surface model, SANISAND-F. Assessments are conducted in terms of stability, accuracy, computational efficiency, and both local and global performance.</div><div>In the implicit implementation, the critical importance of accurate gradient calculations is highlighted, introducing a verification procedure that enables quadratic convergence. The explicit and implicit implementations exhibit stability, producing smooth and bounded responses across a wide range of strain increments and numerical tolerances. However, their accuracy differs significantly. The implicit implementation is sensitive to the initial loading state, strain increment, and loading direction, showing minor dependence on the solver tolerance. In contrast, the explicit implementation is influenced by both strain increment and substepping tolerance, and at practical tolerance and strain increment levels, it often outperforms the implicit scheme in accuracy. Regarding efficiency, the explicit implementation proves more efficient at the local integration level. However, at the global level, the implicit implementation with the consistent tangent exhibits a faster rate of convergence in global equilibrium iterations. Nonetheless, the overall computational cost at the global level is not definitive when comparing explicit and implicit schemes; it varies with simulations and loading-specific parameters, as demonstrated through the included boundary-value problems.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107911"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080667","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

Extension of explicit Runge-Kutta substepping stress integration for viscoplastic model of saturated soils 饱和土粘塑性模型显式Runge-Kutta分段应力积分的推广

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

Computers and Geotechnics

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

Wei CHENG , Zhen-Yu YIN

Stable integration schemes are critically important for rate-dependent constitutive models, serving as a cornerstone for ensuring accuracy, efficiency, and robustness in finite element implementations. This paper investigates the numerical performance of explicit stress integration schemes with adaptive substepping for integrating a newly proposed fractional consistency two-surface viscoplastic model for saturated clays. The incremental stress–strain-strain rate relation of the model can be linearized following the consistency condition of the rate-dependent loading surface and subsequently integrated using four distinct explicit Runge-Kutta substepping integration algorithms (i.e., RK12, RK23, RK34, RK45) with automatic error control and stress drift correction techniques. The overall numerical performance of the algorithms in terms of accuracy and efficiency is evaluated at both the material point level (i.e., isotropic, oedometric, and triaxial compression tests) and the boundary-value problem level (i.e., piezocone penetration and underground gallery excavation), which demonstrates that the RK23 and RK34 algorithms perform excellently in balancing accuracy and computational cost. The proposed algorithms provide a versatile and adaptive framework for integrating time-dependent constitutive equations, particularly those based on the consistency viscoplastic approaches commonly used in advanced rate-dependent modeling, allowing for a wide range of geotechnical engineering applications.

稳定的积分方案对于速率相关的本构模型至关重要，是确保有限元实现中准确性、效率和鲁棒性的基石。本文研究了带自适应分步法的显式应力积分方案对新提出的饱和黏土分数浓度双表面粘塑性模型的数值性能。模型的增量应力-应变-应变率关系可根据速率相关加载面一致性条件进行线性化处理，并采用四种不同的显式龙格-库塔分步积分算法（RK12、RK23、RK34、RK45），结合自动误差控制和应力漂移校正技术进行积分。在材料点水平（即各向同性、测径和三轴压缩试验）和边值问题水平（即压电锥穿透和地下巷道开挖）上对算法在精度和效率方面的整体数值性能进行了评估，表明RK23和RK34算法在平衡精度和计算成本方面表现出色。所提出的算法为积分时变本构方程提供了一个通用的自适应框架，特别是那些基于先进速率相关建模中常用的一致性粘塑性方法的框架，允许广泛的岩土工程应用。

{"title":"Extension of explicit Runge-Kutta substepping stress integration for viscoplastic model of saturated soils","authors":"Wei CHENG , Zhen-Yu YIN","doi":"10.1016/j.compgeo.2026.107937","DOIUrl":"10.1016/j.compgeo.2026.107937","url":null,"abstract":"<div><div>Stable integration schemes are critically important for rate-dependent constitutive models, serving as a cornerstone for ensuring accuracy, efficiency, and robustness in finite element implementations. This paper investigates the numerical performance of explicit stress integration schemes with adaptive substepping for integrating a newly proposed fractional consistency two-surface viscoplastic model for saturated clays. The incremental stress–strain-strain rate relation of the model can be linearized following the consistency condition of the rate-dependent loading surface and subsequently integrated using four distinct explicit Runge-Kutta substepping integration algorithms (i.e., RK12, RK23, RK34, RK45) with automatic error control and stress drift correction techniques. The overall numerical performance of the algorithms in terms of accuracy and efficiency is evaluated at both the material point level (i.e., isotropic, oedometric, and triaxial compression tests) and the boundary-value problem level (i.e., piezocone penetration and underground gallery excavation), which demonstrates that the RK23 and RK34 algorithms perform excellently in balancing accuracy and computational cost. The proposed algorithms provide a versatile and adaptive framework for integrating time-dependent constitutive equations, particularly those based on the consistency viscoplastic approaches commonly used in advanced rate-dependent modeling, allowing for a wide range of geotechnical engineering applications.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"193 ","pages":"Article 107937"},"PeriodicalIF":6.2,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146080669","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