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.     

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.     

A thermo‐hydro‐mechanical model for unsaturated soils based on the effective stress concept

A simple thermo‐hydro‐mechanical (THM) constitutive model for unsaturated soils is described. The effective stress concept is extended to unsaturated soils with the introduction of a capillary stress. This capillary stress is based on a microstructural model and calculated from attraction forces due to water menisci. The effect of desaturation and the thermal softening phenomenon are modelled with a minimal number of material parameters and based on existing models. THM process is qualitatively and quantitatively modelled by using experimental data and previous work to show the application of the model, including a drying path under mechanical stress with transition between saturated and unsaturated states, a heating path under constant suction and a deviatoric path with imposed suction and temperature. The results show that the present model can simulate the THM behaviour in unsaturated soils in a satisfactory way. Copyright © 2010 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.     

Stress history effects on 1‐D consolidation of soft soils: a rheological model

Stress history plays an important role in controlling the consolidation behavior of soft clays, but few models exist that can provide quantitative estimate of its influence. In this paper, the Gibson–Lo rheological model is used to simulate the coupled processes of drainage and creep of soft soils that takes stress history into account. A hybrid combination of analytical and numerical methods is adopted to solve the governing equations of consolidation with the nonlinear rheological model. The methodology is applied to a saturated soft soil subjected to surface loading. The soil profile is separated into normally consolidated and overconsolidated layers by a boundary that is allowed to move. Comparisons of the model predictions and its simulations are used to evaluate the effects of stress history, model parameters, and loading pattern on consolidation behavior. It is shown that stress history influences the location of the moving boundary, variations of the profiles of excess pore water pressure dissipation, stress and deformation‐based average degrees of consolidation. Parametric studies conducted show that when soil is stiffer, the excess pore water pressure dissipates much more quickly, and thus the soil consolidates much faster especially at the early stages. The results also show that soil viscosity influences the deformation‐based average degree of consolidation at the latter stages. The consolidation process of soil layer under linear loading is shown to lag behind those under instantaneous loading: the longer the loading period is, the smaller the average degrees of consolidation are no matter how they are defined. Copyright © 2012 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.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Active Soils of the Niger Delta in Road Pavement Design and Construction

Active soils existing to irreplaceable extents (both laterally and vertically) within the freshwater zone of the Niger delta are known and have been shown to be highly undesirable in their natural forms for road construction. This is because they swell extensively in the rainy season and shrink proportionately in the dry months together with any structure including road pavement constructed with or placed on them. In this paper, four samples of these soils collected from across the Niger delta were subjected to chemical stabilization with slaked-lime and calcium chloride (CaCl₂) to reduce swelling potential while resulting mixtures (corresponding to tolerable swelling limits) were stabilized with Portland cement to improve strength. First, swell indicator (i.e. 24-h free swell index) was found to diminish with increasing stabilizer content while for combined plot (independent of location/percentage fines), about 10.7% and 6% slaked-lime and CaCl₂ were found to be required to achieve the minimum (2.5% free swell) tolerable, normal or non-swell soil condition. Under this condition and for individual soil materials, addition of cement was found to enhance continual increase in MDD, soaked-CBR and durability index such that the soil becomes better compactible, stronger and more durable with increasing cement content. About 2.6% and 14.6% cement contents were found needed to respectively produce sub-base and base-course materials in roadworks with the latter much more durable than the former. The graphs so obtain could also constitute a complete graphical model for preparing deltaic active soils for use as road pavement materials.     

Forecasting vertical ground surface movement from shrinking/swelling soils with artificial neural networks

Artificial neural networks (ANNs) are used to estimate vertical ground surface movement when soils expand and contract due to changes in soil moisture content caused by changing climate conditions. Several counterpropagation ANN test cases were investigated to map climate data (i.e. temperature and rainfall) to vertical ground surface movement at field sites in Texas and Australia. Three of the four ANN test cases use a historical time series of climate data to forecast ground surface elevation relative to a specified datum. The fourth ANN test case predicts the rate of ground surface movement, and requires post‐processing of the predicted rates to calculate ground surface elevation relative to a specified datum. The counterpropagation network has demonstrated a successful mapping of temperature and rainfall data to vertical ground surface movement at a field site when it is trained with a subset of data from the same field site (test cases 1 and 2). The results of training an ANN on one field site and testing it on another field site (test cases 3 and 4) demonstrate the ability of the ANN to capture trends in vertical ground surface movement. When compared with the predictions from a physics‐based method (shrink test‐water content method) that requires measurements/estimates of changes in soil water content, the ANN‐based predictions (based on climatic changes) captured the trends in the field measurements of shrinking–swelling soil surface movements equally well. These findings are promising and merit further investigation with data from additional field sites. Copyright © 2007 John Wiley & Sons, Ltd.     

Discussing the definition of the second‐order work for unsaturated soils

The volumetric compaction due to wetting processes is a phenomenon observed quite often in unsaturated soils. Under certain circumstances, saturation events can result into a sudden and unexpected collapse of the system. These phenomena are usually referred to as wetting‐induced collapses, without providing any detailed theoretical justification for this terminology. In order to predict in a general fashion the occurrence of coupled instabilities induced by saturation processes, a generalization of the theoretical approaches usually employed for saturated geomaterials is here provided. More specifically, this paper addresses the problem of hydro‐mechanical instability in unsaturated soils from an energy standpoint. For this purpose, an extension of the definition of the second‐order work is here suggested for the case of unsaturated porous media. On the basis of some examples of numerical simulations of laboratory tests, coupled hydro‐mechanical instabilities are then interpreted in the light of this second‐order energy measure. Finally, the implications of the theoretical results here presented are commented from a constitutive modelling perspective. Two possible alternative approaches to formulate incremental coupled constitutive relations are indeed discussed, showing how the onset of hydro‐mechanical instabilities can be predicted using an extended form of Hill's stability criterion. Copyright © 2010 John Wiley & Sons, Ltd.     

A double structure generalized plasticity model for expansive materials

The constitutive model presented in this work is built on a conceptual approach for unsaturated expansive soils in which the fundamental characteristic is the explicit consideration of two pore levels. The distinction between the macro‐ and microstructure provides the opportunity to take into account the dominant phenomena that affect the behaviour of each structural level and the main interactions between them. The microstructure is associated with the active clay minerals, while the macrostructure accounts for the larger‐scale structure of the material. The model has been formulated considering concepts of classical and generalized plasticity theories. The generalized stress–strain rate equations are derived within a framework of multidissipative materials, which provides a consistent and formal approach when there are several sources of energy dissipation. The model is formulated in the space of stresses, suction and temperature; and has been implemented in a finite element code. The approach has been applied to explaining and reproducing the behaviour of expansive soils in a variety of problems for which experimental data are available. Three application cases are presented in this paper. Of particular interest is the modelling of an accidental overheating, that took place in a large‐scale heating test. This test allows the capabilities of the model to be checked when a complex thermo‐hydro‐mechanical (THM) path is followed. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of lime and fly ash on swell,consolidation and shear strength characteristics of expansive clays: a comparative study

Expansive soils swell on absorbing water and shrink on evaporation thereof. Because of this alternate swelling and shrinkage, civil engineering structures founded in them are severely damaged. For counteracting the problems of expansive soils, different innovative techniques were suggested. Stabilization of expansive clays with various additives has also met with considerable success. This paper presents, by comparison, the effect of lime and fly ash on free swell index (FSI), swell potential, swelling pressure, coefficient of consolidation, compression index, secondary consolidation characteristics and shear strength. Lime content (weight of lime/weight of dry soil) was varied as 0%, 2%, 4%?and 6%?and fly ash content (weight of fly ash/weight of dry soil) as 0%, 10%?and 20%. A fly ash content of 20%?showed significant reduction in swell potential, swelling pressure, compression index and secondary consolidation characteristics and resulted in increase in maximum dry density and shear strength. Swell potential and swelling pressure decreased with increase in lime content also. Further, consolidation characteristics improved. Compaction characteristics and unconfined compression strength improved at 4%?lime and reduced at 6%?lime.     

A thermo‐mechanical model for variably saturated soils based on hypoplasticity

The paper presents a mechanical model for non‐isothermal behaviour of unsaturated soils. The model is based on an incrementally non‐linear hypoplastic model for saturated clays and can therefore tackle the non‐linear behaviour of overconsolidated soils. A hypoplastic model for non‐isothermal behaviour of saturated soils was developed and combined with the existing hypoplastic model for unsaturated soils based on the effective stress principle. Features of the soil behaviour that are included into the model, and those that are not, are clearly distinguished. The number of model parameters is kept to a minimum, and they all have a clear physical interpretation, to facilitate the model usefulness for practical applications. The step‐by‐step procedure used for the parameter calibration is described. The model is finally evaluated using a comprehensive set of experimental data for the thermo‐mechanical behaviour of an unsaturated compacted silt. Copyright © 2011 John Wiley & Sons, Ltd.     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

Analysis of the FEBEX multi‐barrier system including thermoplasticity of unsaturated bentonite

Deep geological repository involving a multibarrier system constitutes one of the most promising options for isolating high‐level radioactive waste from the human environment. To certify the efficiency of waste isolation, it is essential to understand the behaviour of confining geomaterial under a variety of environmental conditions. To this end, results from a near‐to‐real experiment, the full‐scale engineered barriers in situ experiment, are studied by means of a thermo–hydro–mechanical finite element approach, including a consistent thermoplastic constitutive model for unsaturated soils. Laboratory tests are simulated to calibrate model parameters. The results of the numerical simulations are compared with sensor measurements and show the ability of the model to reproduce the main behavioural features of the system. The influence of the hysteretic and temperature‐dependent retention of water on the mechanical response is exhibited. Finally, those results are interpreted in the light of thermoplasticity of unsaturated soils, which reveals the highly coupled and non‐linear characters of the processes encountered. Copyright © 2011 John Wiley & Sons, Ltd.     

膨胀土胀缩性评价中有关问题的研究

膨胀土吸水膨胀、失水收缩,产生地表变形,进而导致建筑物破坏。本文就现行膨胀土胀缩性指标试验及地基变形量评价计算中有关问题进行了研究,认为以天然含水量为起始含水量测定土样胀缩率是不合适的,建议以缩限为起点进行膨胀试验,以胀限为起点进行收缩试验,采用加荷膨胀量（而非卸荷膨胀量）参加有荷载下的膨胀率计算。膨胀土地基勘察必须结合场地条件分析可能影响土体含水量变化的因素,对同一工程采用几种方法评价,互相对比、参照。     

A random solid‐porous model to simulate the retention curves of soils

This paper describes a random solid‐porous model capable of simulating the structure of porous materials. To this purpose, the grain and pore size distributions as well as the void ratio of the material are required. Solids and pores are distributed at random in the model's space according to a size strategy. Herein, the model is used to simulate the retention curves of soils. The Laplace equation is used to determine the size of the pores able to saturate or dry during a wetting or drying process, respectively. The continuous path principle is used to define those elements that effectively saturate or dry during these processes. With this procedure, it is possible to simulate the main retention curves as well as the scanning curves during wetting–drying cycles. Some experimental results reported in the international literature have been used to test the model. This model can be enhanced to study the mechanical behavior of unsaturated soils. Copyright © 2012 John Wiley & Sons, Ltd.