{"title":"Enhancing accuracy and efficiency in cyclic liquefaction modeling: An automatic calibration framework for advanced constitutive models","authors":"Jan Machaček , Sheng Zeng , Mahdi Taiebat","doi":"10.1016/j.compgeo.2025.107208","DOIUrl":"10.1016/j.compgeo.2025.107208","url":null,"abstract":"<div><div>Precise calibration of constitutive models for cyclic liquefaction is essential but often time-consuming and requires significant expertise, limiting broader application in geotechnical practice. This paper introduces an automatic calibration tool designed to streamline the process for advanced constitutive models under both monotonic and cyclic loading. The tool supports various types of monotonic and cyclic laboratory test data, offers multiple choices of suitable comparison planes for error calculation, with a focus to also suit cyclic liquefaction problems, and employs advanced optimization techniques. The calibration follows a two-stage approach: first, optimizing parameters governing monotonic response using monotonic test data; second, refining these and additional parameters with both monotonic and cyclic data. The critical state parameters are fixed throughout, while the elasticity parameters are fixed in the second stage, all within defined bounds. Using this automatic calibration tool and the adapted calibration strategy, extensive element-level test data was used to determine the parameters of the SANISAND-MSf model for a given sand. These calibrated parameters were then used to simulate boundary value problems, including centrifuge tests of liquefiable sand slopes and sheet-pile-supported liquefiable sand deposits, all subjected to base excitations, demonstrating excellent alignment with experimental results. This validation highlights the robustness, reproducibility, and accuracy of the tool to model cyclic liquefaction while significantly reducing the expertise and time required for calibration. This represents a significant advancement toward the broader adoption of advanced constitutive soil models in geotechnical engineering practice.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107208"},"PeriodicalIF":5.3,"publicationDate":"2025-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143697310","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Submarine slope failure in sensitive clay due to buoyancy force of trapped gas","authors":"Zhipeng Zhu , Dong Wang , Wangcheng Zhang","doi":"10.1016/j.compgeo.2025.107220","DOIUrl":"10.1016/j.compgeo.2025.107220","url":null,"abstract":"<div><div>Submarine landslides pose significant threats to coastal communities and offshore infrastructures. Free gas is regarded as a key triggering or pre-conditioning factor for slope failure. Slope instability caused by reduction of the shear strength due to the presence of gas has been investigated widely. An alternative trigger is the buoyancy force of trapped gas within the gas structure, as observed in both in-situ investigations and model tests. In this study, typical simplified initiation scenarios are investigated under plane strain conditions through finite element analyses, where the trapped gas within a gas structure is simplified as a pressure applied to the adjacent soils and increased to trigger slope failures in sensitive clay. Punching shear failures in the overlying layer and catastrophic landslides are identified under different soil conditions. The latter is attributed to the strength degradation induced by the tensile failure under buoyancy force. Empirical criteria for these two failure patterns are proposed. The proposed criteria are then validated to be applicable to three-dimensional conditions, providing a conservative assessment. A case study of the landslide offshore the northern Svalbard margin is conducted using the criteria proposed, finding that the safety factor of slope may be reduced significantly due to the buoyancy force within the gas structure.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107220"},"PeriodicalIF":5.3,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681402","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}
Xiaohua Bao , Zhizao Bao , Lijuan Li , Peng Peng , Jun Shen , Yingpeng Li , Xiangsheng Chen , Feng Zhang
{"title":"A modified VG model for the SWCC of fiber-reinforced sand","authors":"Xiaohua Bao , Zhizao Bao , Lijuan Li , Peng Peng , Jun Shen , Yingpeng Li , Xiangsheng Chen , Feng Zhang","doi":"10.1016/j.compgeo.2025.107215","DOIUrl":"10.1016/j.compgeo.2025.107215","url":null,"abstract":"<div><div>Fibre reinforcement technology has been widely adopted in soil improvement due to its cost-effectiveness, simplicity, and environmental benefits. In many fibre reinforcement projects, the soil is often in an unsaturated state. However, the numerical simulation mechanisms of fibre-reinforced unsaturated soils remain poorly understood. In this study, a Vangenuchten (VG) model considering fibre incorporating fibres was proposed based on the original VG model. This model considering fibre accurately describes the soil water characteristic curve (SWCC) of fibre-reinforced sand (FRS), as verified by water-holding characteristics tests. Then, unsaturated triaxial tests confirmed the applicability of an unsaturated soil elastoplastic constitutive model and a fully coupled soil–water-air finite element-finite difference (FE-FD) method for simulating the mechanical behaviour of unsaturated FRS. Finally, using the SWCC parameters derived from the VG model considering fibres and mechanical parameters from saturated triaxial tests, slope models were established to analyse the stability of both unreinforced and fibre-reinforced slopes. The results show that the interweaving action of fibres within the soil enhances its strength, reduce permeability, and decreases both saturation and pore water pressure, ultimately increasing slope stability. This study provides valuable insights into the SWCC characteristics and the numerical calculation of FRS under unsaturated conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107215"},"PeriodicalIF":5.3,"publicationDate":"2025-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143696888","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}
Yexin Wan , Weihua Li , Sen Zheng , Zhe Yang , Sainan Zhu
{"title":"Analytical analysis on the influence of groundwater level variation on seismic ground motion under plane wave incidence","authors":"Yexin Wan , Weihua Li , Sen Zheng , Zhe Yang , Sainan Zhu","doi":"10.1016/j.compgeo.2025.107212","DOIUrl":"10.1016/j.compgeo.2025.107212","url":null,"abstract":"<div><div>In this study, a seismic analysis model considering the dynamic interactions among unsaturated porous media, saturated porous media, and single-phase solid media is established to analytically evaluate the seismic response of sites with varying groundwater levels. First, the dynamic characteristics of soils above and below the groundwater table are described based on the mixture theory for unsaturated media and the Biot theory framework, and the governing equations for different media are introduced. Then, the wavefields within the site under obliquely incident plane P/SV waves are analyzed, and the analytical solution for seismic ground motion of the unsaturated–saturated soil–bedrock system is derived based on boundary and continuity conditions. The accuracy and reliability of the proposed analytical solution are validated by comparisons with existing reference solutions. The parameter analysis reveals that groundwater level changes, soil saturation, and porosity significantly influence seismic ground motion. Rising groundwater levels typically reduce displacement amplitudes and increase natural frequencies, while increasing porosity weakens soil stiffness, leading to greater surface displacement. The presence of air phase in unsaturated soil notably affects seismic response, with small saturation changes causing significant variations in displacement. This study provides essential insights into evaluating seismic site safety under varying groundwater conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107212"},"PeriodicalIF":5.3,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681393","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}
{"title":"Seepage features and reliability analysis of deep overburden dam foundations considering the spatial variability in permeability","authors":"Ting Wang, Zilong Wang, Yanlong Li, Haopeng Yuan, Junrui Chai, Yunhe Liu","doi":"10.1016/j.compgeo.2025.107217","DOIUrl":"10.1016/j.compgeo.2025.107217","url":null,"abstract":"<div><div>Seepage stability is an important issue for dams constructed on overburden layers. However, less interest has been focused on the influence of the spatial variability of the overburden layers of dam foundations on seepage features and stability. These layers are too deep such that the cut-off walls cannot reach the rock base because of the limitations of the technology. We propose an efficient way to evaluate the effect of spatial variability on seepage features for different lengths of the cut-off walls, and we perform reliability analysis. By combining the Cholesky decomposition method with a single exponential autocorrelation function (SNX), we generate many random fields for probability analysis. We study the influence of spatial variability and the autocorrelation coefficient of hydraulic conductivity, together with different cut-off wall lengths and drainage prism configurations, on the typical maximum hydraulic gradient (<em>i</em><sub>t</sub> at the base of the cut-off wall and <em>i</em><sub>s</sub> at the downstream overburden) and downstream seepage flow (<em>Q</em><sub>e</sub>). The results show that the deterministic seepage calculation underestimate the critical hydraulic gradient, because soil material heterogeneity often makes the seepage channel more focused. The existence of a drainage prism enlarges the critical maximum hydraulic gradient at the bottom of the cut-off wall. A low risk of seepage failure is shown only at the deep cut-off wall configuration but is also greatly affected by spatial variability. This work improves the understanding of how spatial variability influences the seepage reliability of overburden dams and is critical for controlling the seepage safety of dams with cut-off walls built upon deep overburden foundations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107217"},"PeriodicalIF":5.3,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681403","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}
{"title":"Evaluating clay stiffness effects on offshore pile running with the Coupled Eulerian Lagrangian Method","authors":"A.P. Dyson , A. Tolooiyan , K. Gavin","doi":"10.1016/j.compgeo.2025.107185","DOIUrl":"10.1016/j.compgeo.2025.107185","url":null,"abstract":"<div><div>Driven pipe piles are used extensively in coastal and offshore projects. Traditionally piles with diameters of 2–3 m were common in the offshore wind industry, however the diameter of monopiles to support a 10 MW wind turbine is more commonly 10 m. Offshore wind projects are being developed at sites with very low seabed strengths and pipe pile weights are increasing significantly. Self-weight penetration occurs when the pile is first placed on the seabed. A combination of low strength seabed conditions and increased pile self-weight leads to the risk of pile run (uncontrolled self-weight penetration) during installation at some sites. Predicting pile run risk, run velocities and penetration depths is challenging due to inherent rate effects and the large strains involved. While rapid penetration processes can be considered using both analytic methods and Large Deformation Finite Element simulations, the role of soil rigidity is seldom taken into account, despite known implications from static pile assessments. This study uses large deformation simulation with the Coupled Eulerian Lagrangian method to simulate the pile running process for five well-studied fine-grained soils with varying elastic stiffnesses. Results are compared with analytic methods, highlighting the limitations of current predictive techniques in terms of both the end tip and shaft resistance. As a corollary, a linear trend for the final penetration depth with respect to the logarithm of the soil rigidity index is incorporated in an existing analytic code based on results obtained from large deformation simulations.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107185"},"PeriodicalIF":5.3,"publicationDate":"2025-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681399","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Weitao Yang , Liang Xiao , Zheng Tang , Guoxiong Mei
{"title":"An improved theoretical method for assessing tunnel response to pre-excavation dewatering: Time-dependent deflection, internal force, joint opening, and dislocation","authors":"Weitao Yang , Liang Xiao , Zheng Tang , Guoxiong Mei","doi":"10.1016/j.compgeo.2025.107216","DOIUrl":"10.1016/j.compgeo.2025.107216","url":null,"abstract":"<div><div>Adjacent construction activities can cause significant deformation of the existing shield tunnel, yet the time-dependent development of these deformations is rarely documented, hindering the accurate prediction of potential hazards. This study presents a theoretical model based on the displacement input method to characterize the time-dependent response of tunnel induced by pre-excavation dewatering in an unconfined aquifer. The tunnel and subgrade are modelled as a Timoshenko beam and Pasternak foundation, respectively. The greenfield soil displacement is derived by consolidation theory, incorporating dynamic changes in the phreatic surface and the embedding depth of the waterproof curtain. The proposed solution is evaluated by well documented results of model testing on drawdown and finite element analyses on deformation of both soil and tunnel. Parametric assessments of tunnel deformation are conducted, analyzing the time-dependent response and influences of factors including the tunnel’s relative position to the dewatering zone, soil modulus, specific yield, and the embedding depth of the waterproof curtain. Results indicate that accounting for time-dependent effects significantly reduces the overestimation of tunnel deformation prior to excavation. Additionally, higher soil modulus and greater curtain embedding depth decrease final tunnel deformation, while tunnel position and specific yield primarily influence deformation distribution without altering maximum deformation. The findings of the study provide a more accurate basis for designing dewatering strategies and offer improved prediction for existing tunnel deformation caused by adjacent foundation pit projects.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107216"},"PeriodicalIF":5.3,"publicationDate":"2025-03-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681398","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}
Jiangyong Pu , Qinglei Yu , Hongyuan Liu , Yong Zhao , He Wang , Kai Guan , Yongsheng Cao
{"title":"A coupled humidity diffusion and swelling stress damage model with large deformation analysis and application in studying the swelling behaviors of skarn surrounding a tunnel","authors":"Jiangyong Pu , Qinglei Yu , Hongyuan Liu , Yong Zhao , He Wang , Kai Guan , Yongsheng Cao","doi":"10.1016/j.compgeo.2025.107210","DOIUrl":"10.1016/j.compgeo.2025.107210","url":null,"abstract":"<div><div>This study focuses on numerical modeling of the hygroscopic swelling behavior of skarn rock. A coupled humidity diffusion and swelling stress damage model is proposed by integrating Fick’s second law, the humid-elastic theory and elastic damage mechanics. Geometric nonlinear algorithms are introduced into the coupled model to characterize the swelling large deformation behavior of the skarn rock and varying swelling coefficient is also taken into account. The proposed model is subsequently implemented by the finite element method and validated by comparing the simulated results with well-known analytical solutions, as well as experimental results. Taking a tunnel excavated from skarn rock as an example, the proposed model is applied to investigate the swelling behavior of the skarn rock in humid environments and discuss the effects of the dominant parameters on the swelling behavior. The proposed model can reproduce the swelling deformation and damage of the skarn surrounding the tunnel, which demonstrates good agreement with the in situ monitoring results. Moreover, the damage evolution of the skarn rock induced by humidity diffusion exhibits stepped growth due to the accumulation of swelling stress. The proposed coupled model provides a novel tool for the design and analysis of underground engineering constructions in swelling rocks under humid conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107210"},"PeriodicalIF":5.3,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681400","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}
Jose Salomon , Fernando Patino-Ramirez , Catherine O’Sullivan
{"title":"Stress–dilatancy and micromechanics of sand under partially drained conditions","authors":"Jose Salomon , Fernando Patino-Ramirez , Catherine O’Sullivan","doi":"10.1016/j.compgeo.2025.107200","DOIUrl":"10.1016/j.compgeo.2025.107200","url":null,"abstract":"<div><div>The mechanical behaviour of soil subject to shear loading or deformation is typically considered either completely drained or undrained. Under certain conditions, these drained and undrained scenarios can represent boundaries on the allowed volumetric strain. There is growing interest in exploring the response under intermediate conditions where partial drainage is allowed, particularly in the development of new approaches to mitigate the risk of liquefaction induced failure and the design of off-shore structures. This study uses the discrete element method (DEM) to investigate the effect of partial drainage conditions on the mechanical behaviour of spherical assemblies. Samples with different interparticle friction values are isotropically compressed and then subjected to undrained, drained, and partially drained triaxial shearing. The partially drained conditions are simulated in the DEM samples by applying a controlled volumetric strain that is a fraction of the drained volumetric strain. Results on loose samples indicate that allowing drainage enhances peak shear resistance and can also prevent liquefaction. Moreover, dense samples show a substantial increase in shear resistance when small changes in drainage and volumetric strain take place. The peak stress ratio and the stress ratio at the phase transformation point are insensitive to the drainage level. There is a linear correlation between the state parameter and the drainage level at the peak stress ratio and the phase transformation point. This observation could be used to trace partially drained stress-paths and could also aid the development of uncoupled constitutive models that account for drainage effects.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107200"},"PeriodicalIF":5.3,"publicationDate":"2025-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681401","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A theoretical analysis method for stiffened deep cement mixing (SDCM) pile groups under vertical load in layer soils","authors":"Zhiyu Gong, Haoran Ouyang, Guoliang Dai, Xinsheng Chen","doi":"10.1016/j.compgeo.2025.107211","DOIUrl":"10.1016/j.compgeo.2025.107211","url":null,"abstract":"<div><div>A theoretical analysis method is proposed to forecast the vertical bearing behavior of a long-core SDCM pile group. The load-settlement behavior of a long-core SDCM pile group is different from that of a single long-core SDCM pile due to the existence of the pile group effect. In this study, the nonlinear behaviors of the inner core–cemented soil interface and the inner core–soil interface are expressed via exponential models, whereas the nonlinear relationships of the cemented soil–soil interface and the pile base–soil interface are calculated via an elastic–plastic model. Additionally, the soil between piles is considered a medium that generates additional displacement. Based on the above conditions, the interaction between long-core SDCM pile groups was analyzed. This method was first used on a single long-core SDCM pile, and the results were compared with the analytical solutions and FEM results from previous studies; then, the field test results of long-core SDCM pile groups were compared. A reasonable prediction can be achieved via the method proposed in this article. Finally, the law of additional displacement caused by the pile group effect and the optimal solution for the cemented soil coverage size for the long-core SDCM pile group were obtained by analyzing important parameters, including the pile spacing <em>s</em><sub>ij</sub>, the height ratio of the cemented soil to the PHC pipe pile <em>L<sub>c</sub>/L<sub>p</sub></em>, the radius ratio of the cemented soil to the PHC pipe pile <em>R<sub>c</sub>/R<sub>p</sub></em>, and the slenderness ratio <em>L</em>/<em>R</em><sub>c</sub>.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107211"},"PeriodicalIF":5.3,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143642130","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}