The behavior of soft soils distributed in coastal areas usually exhibits obvious time-dependent behavior after loading. To reasonably describe the stress-strain relationship of soft soils, this paper establishes a viscoelastic-viscoplastic small-strain constitutive model based on the component model and the hardening soil model with small-strain stiffness (HSS model). First, the Perzyna's viscoplastic flow rule and the modified Hardin-Drnevich model are introduced to derive a one-dimensional incremental Nishihara constitutive equation. Next, the flexibility coefficient matrix is utilized to extend the one-dimensional model to three-dimensional conditions. Then, by combining the HSS elastoplastic theory with the component model, the viscoelastic-viscoplastic small-strain constitutive model is subsequently established. To implement the proposed model for numerical analysis, the corresponding UMAT subroutine is developed using Fortran. After comparing the results of numerical simulations with those of existing literature, the reliability of the constitutive model and the program written in this paper is verified. Finally, numerical examples are designed to further analyze the effects of small-strain parameters and viscoelastic-viscoplastic parameters on the time-dependent behavior of soft soils.

Investigation of thermal effects on the strain rate-dependent properties of compacted bentonite is crucial for the long-term safety assessment of deep geological repository for disposal of high-level radioactive waste. In the present work, cylindrical GMZ01 bentonite specimens were compacted with suctioncontrolled by the vapor equilibrium technique. Then, a series of temperature- and suction-controlled stepwise constant rate of strain (CRS) tests was performed and the rate-dependent compressibility behavior of the highly compacted GMZ01 bentonite was investigated. The plastic compressibility parameter l, the elastic compressibility parameter k, the yield stress p0, as well as the viscous parameter a were determined. Results indicate that l, k and a decrease and p0 increases as suction increases. Upon heating, parameters l, a and p0 decrease. It is also found that p0 increases linearly with increasing CRS in a double-logarithm coordinate. Based on the experimental results, a viscosity parameter a(s, T) was fitted to capture the effects of suction s and temperature Ton the relationship between yield stress and strain rate. Then, an elastic-thermo-viscoplastic model for unsaturated soils was developed to describe the thermal effects on the rate-dependent behavior of highly compacted GMZ01 bentonite. Validation showed that the calculated results agreed well to the measured ones. (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/).

A viscoplastic constitutive model is developed to describe the viscoplastic behavior of unsaturated soils. The proposed model accounts for the stain-rate and suction effects on yield by adopting an unsaturated isotach concept. The nonstationary flow surface theory (NSFS) is applied for modeling the viscoplastic behavior, with the yield surface which can evolve with the viscoplastic strain, viscoplastic strain rate and suction. Meanwhile, the progressively hardening concept is adopted for reproducing the viscoplastic behavior of soil at an over-consolidated state. For the validation, a series of loading conditions are considered based on the data from the literature. Results show that the proposed model is able to reproduce the main viscoplastic behaviors of unsaturated highly compacted Gaomiaozi (GMZ) bentonite, Glenroy silt and Qianjiangping landslide (QL) soil, including CRS compression tests, rate-dependent triaxial shear tests, and triaxial creep tests.

The paper is dedicated to developing a comprehensive analysis method of the criteria for defining the compressible thickness critical for estimating long-term settlements in buildings and structures situated on soft soils, focusing on their creep behavior. This study introduces an engineering method grounded on the criterion of soil's undrained condition within the mass, considering both elastic and residual deformations through equivalent creep deformations. Unlike previous methodologies, the proposed method facilitates the assessment of long-term settlements by incorporating creep effects over time, employing undrained shear strength for both normally consolidated and overconsolidated soils. The method enables settlement calculations based on static-sounding data, enhancing predictions' accuracy and reliability. This research endeavors to broaden the application of numerical and analytical calculations in real-world practices, employing elastoviscoplastic soil models to design structures on weak foundations.

In order to evaluate the performance of deep geological disposal of radioactive waste, an underground research laboratory (URL) was constructed by Andra in the Callovo-Oxfordian (COx) claystone formation at the Meuse/Haute-Marne (MHM). The construction of URL induced the excavation damage of host formations, and the ventilation in the galleries desaturated the host formation close to the gallery wall. Moreover, it is expected that the mechanical behaviour of COx claystone is time -dependent. This study presents a constitutive model developed to describe the viscoplastic behaviour of unsaturated and damaged COx claystone. In this model, the unsaturation effect is considered by adopting the Bishop effective stress and the van Genuchten (VG) water retention model. In terms of the viscoplastic behaviour, the nonstationary flow surface (NSFS) theory for unsaturated soils is used with consideration of the coupled effects of strain rate and suction on the yield stress. A progressive hardening law is adopted. Meanwhile, a non -associated flow rule is used, which is similar to that in Barcelona basic model (BBM). In addition, to describe the damage effect induced by suction change and viscoplastic loading, a damage function is defined based on the crack volume proportion. This damage function contains two variables: unsaturated effective stress and viscoplastic volumetric strain, with the related parameters determined based on the mercury intrusion porosimetry (MIP) tests. For the model validation, different tests on COx claystone under different loading paths are simulated. Comparisons between experimental and simulated results indicated that the present model is able to well describe the viscoplastic behaviour of damaged COx claystone, including swelling/shrinkage, triaxial extension and compression, and triaxial creep. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting 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/).