Revisiting the thermodynamics of hardening plasticity for unsaturated soils

A thermodynamically consistent extension of the constitutive equations of saturated soils to unsaturated conditions is often worked out through the use of a unique ‘effective’ interstitial pressure, accounting equivalently for the pressures of the saturating fluids acting separately on the internal solid walls of the pore network. The natural candidate for this effective interstitial pressure is the space averaged interstitial pressure. In contrast experimental observations have revealed that, at least, a pair of stress state variables was needed for a suitable framework to describe stress–strain–strength behaviour of unsaturated soils. The thermodynamics analysis presented here shows that the most general approach to the behaviour of unsaturated soils actually requires three stress state variables: the suction, which is required to describe the invasion of the soil by the liquid water phase through the retention curve; two effective stresses, which are required to describe the soil deformation at water saturation held constant. However a simple assumption related to the plastic flow rule leads to the final need of only a Bishop-like effective stress to formulate the stress–strain constitutive equation describing the soil deformation, while the retention properties still involve the suction and possibly the deformation. Commonly accepted models for unsaturated soils, that is the Barcelona Basic Model and any approach based on the use of an effective averaged interstitial pressure, appear as special extreme cases of the thermodynamic formulation proposed here.     

Review of fundamental principles in modelling unsaturated soil behaviour

An unsaturated soil is a state of the soil. All soils can be partially saturated with water. Therefore, constitutive models for soils should ideally represent the soil behaviour over entire ranges of possible pore pressure and stress values and allow arbitrary stress and hydraulic paths within these ranges. The last two decades or so have seen significant advances in modelling unsaturated soil behaviour. This paper presents a review of constitutive models for unsaturated soils. In particular, it focuses on the fundamental principles that govern the volume change, shear strength, yield stress, water retention and hydro-mechanical coupling. Alternative forms of these principles are critically examined in terms of their predictive capacity for experimental data, the consistency between these principles and the continuity between saturated and unsaturated states.     

非饱和土弹塑性模型参数的试验确定及有限元法

简要介绍广泛应用的非饱和土Alonso模型的原理与发展,并通过3组共7个试样的非饱和土三轴试验,得到模型的11个参数。在饱和土Cam本构关系的基础上,推导出一个基于非饱和土Alonso模型的应力-应变增量方程的计算公式,该方程与饱和土的本构方程形式相同。编制了能够考虑净应力和吸力二者对土体硬化规律影响的有限元程序,从而为非饱和土弹塑性计算提供一条途径,便于进一步应用于实际工程。     

Formulation of a three‐dimensional constitutive model for unsaturated soils incorporating mechanical–water retention couplings

Wheeler, Sharma and Buisson proposed an elasto‐plastic constitutive model for unsaturated soils that couples the mechanical and water retention behaviours. The model was formulated for isotropic stress states and adopts the mean Bishop's stress and modified suction as stress state variables. This paper deals with the extension of this constitutive model to general three‐dimensional stress conditions, proposing the generalized stress–strain relationships required for the numerical integration of the constitutive model. A characteristic of the original model is the consideration of a number of elasto‐plastic mechanisms to describe the complex behaviour of unsaturated soils. This work presents the three‐dimensional formulation of these coupled irreversible mechanisms in a generalized way including anisotropic loading. The paper also compares the results from the model with published experiments performed under different loading conditions. The response of the model is very satisfactory in terms of both mechanical and water retention behaviours. Copyright © 2013 John Wiley & Sons, Ltd.     

A density‐dependent elastoplastic hydro‐mechanical model for unsaturated compacted soils

This paper presents a three‐dimensional elastoplastic constitutive model for predicting the hydraulic and mechanical behaviour of unsaturated soils. It is based on experimental results obtained from a series of controlled‐suction triaxial tests on unsaturated compacted clay with different initial densities. Hydraulic hysteresis in the water‐retention behaviour is modelled as an elastoplastic process, with the elastic part modelled by a series of scanning curves and the elastoplastic part modelled by the main drying and wetting curves. The effect of void ratio on the water‐retention behaviour is studied using data obtained from controlled‐suction wetting–drying cyclic tests on unsaturated compacted clay with different initial densities. The effect of the degree of saturation on the stress–strain‐strength behaviour and the effect of void ratio on the water‐retention behaviour are considered in the model, as is the effect of suction on the hydraulic and mechanical behaviour. The initial density dependency of the compacted soil behaviour is modelled by experimental relationships between the initial density and the corresponding yield stress and, thereby, between the initial density and the normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure and yield surfaces in the deviatoric stress plane are given by the Matsuoka–Nakai criterion. Model predictions of the stress–strain and water‐retention behaviour are compared with those obtained from triaxial tests with different initial densities under isotropic compression, triaxial compression and triaxial extension, with or without variation in suction. The comparisons indicate that the model accurately predicts the hydraulic and mechanical behaviour of unsaturated compacted soils with different initial densities using the same material constant. Copyright © 2006 John Wiley & Sons, Ltd.     

On the use of effective stress in three-dimensional hydro-mechanical coupled model

In the last decades, a number of hydro-mechanical elastoplastic constitutive models for unsaturated soils have been proposed. Those models couple the hydraulic and mechanical behaviour of unsaturated soils, and take into account the effects of the degree of saturation on the stress–strain behaviour and the effects of deformation on the soil–water characteristic response with a simple reversible part for the hysteresis. In addition, the influence of the suction on the stress–strain behaviour is considered. However, until now, few models predict the stress–strain and soil–water characteristic responses of unsaturated soils in a fully three-dimensional Finite Element code. This paper presents the predictions of an unsaturated soil model in a Three-dimensional Framework, and develops a study on the effect of partial saturation on the stability of shallow foundation resting on unsaturated silty soil. Qualitative predictions of the constitutive model show that incorporating a special formulation for the effective stress into an elastoplastic coupled hydro-mechanical model opens a full range of possibilities in modelling unsaturated soil behaviour.     

ACMEG‐TS: A constitutive model for unsaturated soils under non‐isothermal conditions

This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo‐mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto‐plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo‐hydro‐mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non‐isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi‐mechanism plasticity, whereas water retention characteristics are described using elasto‐plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 John Wiley & Sons, Ltd.     

Coupled modelling of hydro‐mechanical behaviour of unsaturated compacted expansive soils

The Barcelona basic model cannot predict the mechanical behaviour of unsaturated expansive soils, whereas the Barcelona expansive model (BExM) can only predict the stress–strain behaviour of unsaturated expansive soils without the water‐retention behaviour being incorporated. Moreover, the micro‐parameters and the coupling function between micro‐structural and macro‐structural strains in the BExM are difficult to determine. Experimental data show that the compression curves for non‐expansive soils under constant suctions are shifted towards higher void ratios with increasing suction, whereas the opposite is true for expansive soils. According to the observed water‐retention behaviour of unsaturated expansive soils, the air‐entry value increases with density, and the relationship between the degree of saturation and void ratio is linear at constant suction. According to the above observation, an elastoplastic constitutive model is developed for predicting the hydraulic and mechanical behaviour of unsaturated expansive soils, based on the existing hydro‐mechanical model for non‐expansive unsaturated soil. The model takes into consideration the effect of degree of saturation on the mechanical behaviour and that of void ratio on the water‐retention behaviour. The concept of equivalent void ratio curve is introduced to distinguish the plastic potential curve from the yield curve. The model predictions are compared with the test results of an unsaturated expansive soil, including swelling tests under constant net stress, isotropic compression tests and triaxial shear tests under constant suction. The comparison indicates that the model offers great potential for quantitatively predicting the hydraulic and mechanical behaviour of unsaturated expansive soils. Copyright © 2011 John Wiley & Sons, Ltd.     

吸力历史对非饱和土力学性质的影响

现在被广泛公认的由Fredlund提出的非饱和土力学的双参数理论,即净应力和吸力为非饱和土的应力状态变量,不能直接考虑吸力历史及其饱和度对非饱和土的应力-应变关系和强度的影响。非饱和土三轴试验结果表明,即使净应力和吸力相同的条件下,经过干-湿循环试样与未经过干-湿循环试样的应力比-应变关系和强度是不相同的。在其他条件相同时,经历过干-湿循环的试样比未经过干-湿循环试样的应力比-应变关系要高、强度大和体变小。经过干-湿循环试样的饱和度低而强度高,主要是由于经过先期较高的吸力,相当于受过较大的前期有效压力,使试样成为超固结土。更多不同吸力历史的对比试验有待于进一步研究,以便为非饱和土的水力-力学特性耦合弹塑性本构模型定量地表示上述非饱和土的性质提供基础性试验数据。     

Effect of degree of saturation on mechanical behaviour of unsaturated soils and its elastoplastic simulation

The two stress-state variable approach has been widely used in interpreting unsaturated soil behaviour. However this approach cannot take into account the effect of degree of saturation or water contents on the stress–strain behaviour and strength of unsaturated soils. The triaxial test results presented in this paper show that even if the same path of net stress and suction is followed, the stress–strain relation and strength are different due to different degrees of saturation. When other conditions are the same, the higher the degree of saturation for the soil sample is, the higher the stress ratio corresponding to a given axial strain will be. This effect can be modeled by using an elasto-plastic constitutive model coupling hydraulic and mechanical behaviour of unsaturated soils. Comparisons between the predicted and measured results are presented, which demonstrate that the model can quantitatively simulate the influence of the degree of saturation on stress–strain behaviour and strength of unsaturated soils.     

A bounding surface plasticity model for unsaturated soil at small strains

Most existing hydromechanical models for unsaturated soils are not able to fully capture the nonlinearity of stress–strain curves at small strains (less than 1%). They cannot therefore, for example, accurately predict ground movements and the performance of many earth structures under working conditions. To tackle this problem, a state‐dependent bounding surface plasticity model has been newly developed. Particularly, the degradation of shear modulus with strain at small strains ranging from 0.001% to 1% is focused. The proposed model is formulated in terms of mean average skeleton stress, deviator stress, suction, specific volume and degree of saturation. Void ratio‐dependent hydraulic hysteresis is coupled with the stress–strain behaviour. Different from other elastoplastic models for unsaturated soils, plastic strains are allowed inside bounding surfaces. In this paper, details of model formulations and calibration procedures of model parameters are presented. To evaluate the capability of the new model, it is applied to simulate a series of triaxial compression tests on compacted unsaturated silt at various suctions. Effects of suction, drying and wetting as well as net stress on unsaturated soil behaviour are well captured. The model shows good predictions of the degradation of shear modulus with strain over a wide range of strains from 0.001% to 1%. Copyright © 2015 John Wiley & Sons, Ltd.     

A unified bounding surface plasticity model for unsaturated soils

A unified constitutive model for unsaturated soils is presented in a critical state framework using the concepts of effective stress and bounding surface plasticity theory. Consideration is given to the effects of unsaturation and particle crushing in the definition of the critical state. A simple isotropic elastic rule is adopted. A loading surface and a bounding surface of the same shape are defined using simple and versatile functions. The bounding surface and elastic rules lead to the existence of a limiting isotropic compression line, towards which the stress trajectories of all isotropic compression load paths approach. A non‐associated flow rule of the same general form is assumed for all soil types. Isotropic hardening/softening occurs due to changes in plastic volumetric strains as well as suction for some unsaturated soils, enabling the phenomenon of volumetric collapse upon wetting to be accounted for. The model is used to simulate the stress–strain behaviour observed in unsaturated speswhite kaolin subjected to three triaxial test load paths. The fit between simulation and experiment is improved compared to that of other constitutive models developed using conventional Cam‐Clay‐based plasticity theory and calibrated using the same set of data. Also, the model is used to simulate to a high degree of accuracy the stress–strain behaviour observed in unsaturated Kurnell sand subjected to two triaxial test load paths and the oedometric compression load path. For oedometric compression theoretical simulations indicate that the suction was not sufficiently large to cause samples to separate from the confining ring. Copyright © 2005 John Wiley & Sons, Ltd.     

A kinematic hardening based model for unsaturated soils considering different hydraulic conditions

At present, several of the existing elastoplastic constitutive models are adapted for describing the stress–strain behavior of unsaturated soils. However, most of them present certain limitations in this field. These limitations can be related to the basic model and/or added unsaturated state variables and formulations. In this regard, inability to model the hydro‐mechanical behavior in constant water (CW) conditions is an example of these limitations. In this paper, an advanced version of CJS model is selected for adaptation to the unsaturated states. Adaptation to unsaturated states is achieved in the framework of effective stress approach. Effective stress equation and unsaturated state variables are selected based on the recent research existing in the literature. The developed model is capable of describing the complex behavior of unsaturated soil in the CW condition in addition to predicting the behavior at failure and post–failure, nonlinear elastoplastic behavior at low levels of stress and strain (by selecting a very small elastic domain), as well as wetting and collapse behaviors. In order to validate the model, results of triaxial tests in CD and CW conditions are used. The validation results indicate the good capability of the proposed model. Behavior of the unsaturated soils during wetting is an important issue. For this reason, the model is also evaluated based on the results of wetting and collapse triaxial tests. A comparison between the tests and simulation results shows that the model is able to predict the soil behavior under the wetting path. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling the mechanical behaviour of unsaturated soils using a genetic algorithm-based neural network

Modelling the mechanical behaviour of unsaturated soils has been the subject of many research works in the past few decades. A number of constitutive models have been developed to describe the complex behaviour of unsaturated soils. Despite the significant advances in the constitutive theories for unsaturated soils, none of the existing models can completely describe the various aspects of the real behaviour of unsaturated soils. In this paper, a new unified approach is presented, based on the integration of a neural network and a genetic algorithm, for the modelling of unsaturated soils. In the proposed approach, a genetic algorithm was used to optimise the weights of the neural network. A three-layer sequential architecture was chosen for the neural network. The network had eight input neurons, five neurons in the hidden layer and three neurons in the output layer. The eight input neurons represented the initial gravimetric water content, initial dry density, degree of saturation, net mean stress with respect to pore-air pressure, axial strain, deviatoric stress, soil suction and volumetric strain, and the three neurons in the output layer represented the deviatoric stress, suction and volumetric strain at the end of each increment. The network was trained and tested using a database that included results from a comprehensive set of triaxial tests on unsaturated soils from the literature. The predictions of the proposed model were compared with the experimental results. The comparison of the results indicates that the proposed approach was accurate and robust in representing the mechanical behaviour of unsaturated soils.     

A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation

Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities.     

A unified constitutive model for unsaturated soil under monotonic and cyclic loading

This study presents a sophisticated elastoplastic constitutive model for unsaturated soil using Bishop-type skeleton stress and degree of saturation as state variables in the framework of critical state soil mechanism. The model is proposed in order to describe the coupled hydromechanical behavior of unsaturated soil irrespective of what kind of the loadings or the drainage conditions may be. At the same time, a water retention characteristic curve considering the influence of deformation on degree of saturation is also proposed. In the model, the superloading and subloading concepts are introduced to consider the influences of overconsolidation and structure on deformation and strength of soils. The proposed model only employs nine parameters, among which five parameters are the same as those used in Cam-Clay model. The other four parameters have the clear physical meanings and can be easily determined by conventional soil tests. The capability and accuracy of the proposed model have been validated carefully through a series of laboratory tests such as isotropic loading tests and triaxial monotonic and cyclic compression tests under different mechanical and hydraulic conditions.

     

非饱和土的本构模型研究

引用平均土骨架应力的概念,研究推导出非饱和土的刚度参数随吸力变化而变化的关系式,进而推导得到用平均土骨架应力表述的非饱和土LC屈服面函数以及硬化规律。从土力学原理推导,得到土样由于在净应力和吸力作用下产生体积变形引起土样饱和度变化的关系式。由平均土骨架应力推广,得到三轴应力状态的椭圆屈服函数,这一非饱和土本构模型的优点在于考虑了应力作用后土样饱和度的变化,通过对已有试验数据的初步验证,表明提出的非饱和土本构模型的合理性和适用性。     

A three‐phase model for unsaturated soils

A comprehensive framework to define the constitutive behaviour of unsaturated soils is developed within the theory of mixtures applied to three‐phase porous media. Each of the three phases is endowed with its own strain and stress. Elastic and elastic–plastic constitutive equations are developed. Particular emphasis is laid on the interactions between the phases both in the elastic and plastic regimes. Nevertheless, the clear structure of the constitutive equations requires a minimal number of material parameters. Their identification is provided: in particular, it incorporates directly the soil–water characteristic curve. Crucial to the formulation is an appropriate definition of the effective stress. The coupled influence of this effective stress and of suction makes it possible to describe qualitatively many of the characteristic features observed in experiments, e.g. for normally consolidated soils, a plastic behaviour up to air entry followed by an elastic behaviour at increasing suctions, and, on the way back, an elastic behaviour, unless compression is applied in which case plastic collapse occurs. Copyright © 2000 John Wiley & Sons, Ltd.     

New basis for the constitutive modelling of aggregated soils

Natural and compacted soils are usually characterized by aggregation of particles. The mechanical behaviour of these materials depends on soil structure. The oedometric compression tests performed on aggregated samples presented here showed that these materials exhibit a yield limit depending not only on stress history and stress state but also on soil structure. Evidence is provided using the neutron tomography technique. These results revealed that soil structure modification occurs together with plastic deformations. The experimental results are used to propose a new state parameter to quantify the soil structure. Based on pore-scale experimental observations, an evolution law for this parameter is proposed as a function of associated plastic strains. Considering both soil fabric and inter-particle bonding effects, a new yield limit depending on stress state, stress history and soil structure is introduced for the aggregated soils. Accordingly, a new constitutive framework consistent with strain hardening plasticity is proposed to consider soil structure effects in the modelling of aggregated soils.