Numerical challenges, incorporating non-uniqueness, non-convexity, undefined gradients, and high curvature, of the positive level sets of yield function F > 0 are encountered in stress integration when utilizing the return-mapping algorithm family. These phenomena are illustrated by an assessment of four typical yield functions: modified spatially mobilized plane criterion, Lade criterion, Bigoni-Piccolroaz criterion, and micromechanics-based upscaled Drucker-Prager criterion. One remedy to these issues, named the Hop-to-Hug (H2H) algorithm, is proposed via a convexification enhancement upon the classical cutting-plane algorithm (CPA). The improved robustness of the H2H algorithm is demonstrated through a series of integration tests in one single material point. Furthermore, a constitutive model is implemented with the H2H algorithm into the Abaqus/Standard finite-element platform. Element-level and structure-level analyses are carried out to validate the effectiveness of the H2H algorithm in convergence. All validation analyses manifest that the proposed H2H algorithm can offer enhanced stability over the classical CPA method while maintaining the ease of implementation, in which evaluations of the second-order derivatives of yield function and plastic potential function are circumvented. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

This article presents a micro-structure tensor enhanced elasto-plastic finite element (FE) method to address strength anisotropy in three-dimensional (3D) soil slope stability analysis. The gravity increase method (GIM) is employed to analyze the stability of 3D anisotropic soil slopes. The accuracy of the proposed method is first verified against the data in the literature. We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90 degrees turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions. Based on our numerical results, the end effect significantly impacts the failure mechanism and safety factor. Anisotropy degree notably affects the safety factor, with higher degrees leading to deeper landslides. For concave slopes, they can be approximated by straight slopes with suitable boundary conditions to assess their stability. Furthermore, a case study of the Saint-Alban test embankment A in Quebec, Canada, is provided to demonstrate the applicability of the proposed FE model. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

This study investigates the impact of nearby structures on the cyclic settlement mechanisms of shallow foundations in liquefiable soils using a numerical model based on Biot's porous media theory. The model predicts excess pore water pressure and settlement by coupling equilibrium and continuity equations, solved using an implicit time integration scheme. Soil nonlinearity under cyclic loading is represented using generalized plasticity, boundary surfaces, and non-associated models. Three scenarios are simulated to study the effect of spacing between light and heavy foundations and variation in acceleration intensity. Results show that as spacing between foundations increases, lateral displacement and settlement decrease. Excess pore water pressure generation also decreases with increased foundation spacing. Soil just below the foundation exhibits maximum settlement, decreasing with depth. When input acceleration increases from 0.1 g to 0.15 g and 0.2 g, settlement increases by 40%-55% and 90%-110% respectively for both light and heavy foundations, regardless of spacing. Excess pore water pressure also increases sharply with higher acceleration intensity. The findings highlight the importance of considering foundation-soil-foundation interaction effects in liquefaction-prone urban settings and provide insights for designing resilient shallow foundations. The advanced numerical modeling approach offers engineers a more informed way to mitigate liquefaction risk and build safer, more durable structures in earthquake-prone areas.

The multi-dimensional estimation of frost heave deformation is crucial for predicting soil deformation and the pressure generated during the freezing process. This study offers a comprehensive review and analysis of frost heave characteristics in the heat flow direction and the transverse direction to it. Based on the frost heave test results, an innovative method for calculating the anisotropic parameter has been introduced. This method includes only two parameters: one reflects the more pronounced characteristic of frost heave in the direction of heat flow, and the other represents the sensitivity of anisotropic parameters to constraint stress. Through comparative analysis with experimental results, this method can effectively express the evolution of the anisotropic frost heave with changes of the confining stress in different directions. Then, a combined thermalhydraulic-mechanical coupling model is established, offering a way of applying the improved model. The coupling model can predict significant frost heave under conditions of sufficient water supply and effectively captures the frost heave characteristics under various temperature and stress boundary conditions. This research contributes significantly to predicting frost heave deformation in low-temperature natural gas pipelines and calculating the frozen soil pressure exerted on the pipelines.