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‐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.     

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.     

Time-dependent behaviour of deep clays

The time-dependent behaviour of deep saturated clays is related both to the effects of hydraulic diffusion and of viscosity. In this paper, we present first, by means of an experimental approach, the main features of the mechanical properties of deep clays. Then, we study the effects of these specific behaviours on the time-dependent behaviour of underground structures. We present theoretical and numerical investigations of the effects of pore pressure diffusion resulting from a tunnel excavation in a poroplastic or poroviscoplastic medium. Finite-element calculations show that the time-dependent convergence of the tunnel wall is non negligible. Then, we focus on the behaviour of a cylindrical thermal source buried in a deep clay.

Coupling effects between thermal, hydraulic and mechanical behaviours are very important in soft and low permeable deep clays (saturated compressible clays with high porosity). We show in particular that the displacements and the stresses are very sensitive to the constitutive model. Irreversible behaviour may be traduced by appearance of residual stresses within the rockmass.     

A new bounding surface model for thermal cyclic behaviour

To accurately predict soil volume changes under thermal cycles is of great importance for analysing the performance of many earth structures such as the energy pile and energy storage system. Most of the existing thermo‐mechanical models focus on soil behaviour under monotonic thermal loading only, and they are not able to capture soil volume changes under thermal cycles. In this study, a constitutive model is proposed to simulate volume changes of saturated soil subjected to cyclic heating and cooling. Two surfaces are defined and used: a bounding surface and a memory surface. The bounding surface and memory surface are mainly controlled by the preconsolidation pressure (a function of plastic volumetric strain) and the maximum stress experienced by the soil, respectively. Under thermal cycles, the distance of the two surfaces and plastic modulus increase with an accumulation of plastic strain. By adopting the double surface concept, a new elastoplastic model is derived from an existing single bounding surface thermo‐mechanical model. Comparisons between model predictions and experimental results reveal that the proposed model is able to capture soil volume changes under thermal cycles well. The plastic strain accumulates under thermal cycles, but at a decreasing rate, until stabilization. Copyright © 2017 John Wiley & Sons, Ltd.     

A coupled chemo‐thermo‐hygro‐mechanical model of concrete at high temperature and failure analysis

A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Experimental study and constitutive modelling of elasto‐plastic damage in heat‐treated mortar

This study investigates the effect of a heat‐treatment upon the thermo‐mechanical behaviour of a model cement‐based material, i.e. a normalized mortar, with a (w/c) ratio of 0.5. First, a whole set of varied experimental results is provided, in order to either identify or validate a thermo‐mechanical constitutive model, presented in the second paper part. Experimental responses of both hydraulic and mechanical behaviour are given after different heating/cooling cycling levels (105, 200, 300, 400^°C). The reference state, used for comparison purposes, is taken after mass stabilization at 60^°C. Typical uniaxial compression tests are provided, and original triaxial deviatoric compressive test responses are also given. Hydraulic behaviour is identified simultaneously to triaxial deviatoric compressive loading through gas permeability K_gas assessment. K_gas is well correlated with volumetric strain evolution: gas permeability increases hugely when ε_v testifies of a dilatant material behaviour, instead of contractile from the test start. Finally, the thermo‐mechanical model, based on a thermodynamics approach, is identified using the experimental results on uniaxial and triaxial deviatoric compression. It is also positively validated at residual state for triaxial deviatoric compression, but also by using a different stress path in lateral extension, which is at the origin of noticeable plasticity. Copyright © 2009 John Wiley & Sons, Ltd.     

Gas flow in anisotropic claystone: modelling triaxial experiments

Selected gas pulse tests on initially saturated claystone samples under isotropic confinement pressure are simulated using a 3D thermo‐hydro‐mechanical code. The constitutive model considers the hydro‐mechanical anisotropy of argillaceous rocks. A cross‐anisotropic linear elastic law is adopted for the mechanical behaviour. Elements for a proper modelling of gas flow along preferential paths include an embedded fracture permeability model. Rock permeability and its retention curve depend on strains through a fracture aperture. The hydraulic and mechanical behaviours have a common anisotropic structure. Small‐scale heterogeneity is considered to enhance the initiation of flow through preferential paths, following the direction of the bedding planes. The numerical simulations were performed considering two different bedding orientations, parallel and normal to the imposed flow in the test. Simulations are in agreement with recorded upstream and downstream pressures in the tests. The evolution of fluid pressures, degree of saturation, element permeability and stress paths are presented for each case analysed. This information provides a good insight into the mechanisms of gas transport. Different flow patterns are obtained depending on bedding orientation, and the results provide an explanation for the results obtained in the tests. Copyright © 2012 John Wiley & Sons, Ltd.     

Implicit integration of a mixed isotropic–kinematic hardening plasticity model for structured clays

In recent years, a number of constitutive models have been proposed to describe mathematically the mechanical response of natural clays. Some of these models are characterized by complex formulations, often leading to non‐trivial problems in their numerical integration in finite elements codes. The paper describes a fully implicit stress‐point algorithm for the numerical integration of a single‐surface mixed isotropic–kinematic hardening plasticity model for structured clays. The formulation of the model stems from a compromise between its capability of reproducing the larger number of features characterizing the behaviour of structured clays and the possibility of developing a robust integration algorithm for its implementation in a finite elements code. The model is characterized by an ellipsoid‐shaped yield function, inside which a stress‐dependent reversible stiffness is accounted for by a non‐linear hyperelastic formulation. The isotropic part of the hardening law extends the standard Cam‐Clay one to include plastic strain‐driven softening due to bond degradation, while the kinematic hardening part controls the evolution of the position of the yield surface in the stress space. The proposed algorithm allows the consistent linearization of the constitutive equations guaranteeing the quadratic rate of asymptotic convergence in the global‐level Newton–Raphson iterative procedure. The accuracy and the convergence properties of the proposed algorithm are evaluated with reference to the numerical simulations of single element tests and the analysis of a typical geotechnical boundary value problem. Copyright © 2007 John Wiley & Sons, Ltd.     

Global and local sensitivity analyses for coupled thermo–hydro–mechanical problems

This paper focuses on the sensitivity analysis for coupled thermo–hydro–mechanical problems employing both local and global sensitivity methods. A derivative‐based method is used in the local sensitivity approach, whereas the random balance designs method is used for the global sensitivity analysis. The main goal is to investigate the effect of uncertainties in the constitutive parameters on the results from nonlinear coupled thermo–hydro–mechanical analyses of unsaturated soil behavior whose modeling generally involves large sets of constitutive relations. Knowing the parameter sensitivity allows to qualitatively assess the validity of the results obtained by computational simulations of high‐risk situations, for example, emerging nuclear waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.     

A double hardening thermo-mechanical constitutive model for overconsolidated clays

On the basis of a double hardening model for clays and available experimental results, a new thermo-elasto-plastic constitutive model for saturated clays is proposed to describe the effects of temperature and overconsolidation ratio on the mechanical properties of saturated clays. Two hardening parameters are introduced: s_c^￠ {\sigma}_{{\rm c}}^{\prime} and α. The proposed model is then applied to simulate the relevant important features of saturated clays with different overconsolidation ratios under different temperature and loading conditions. The model predictions are compared with available experimental results to demonstrate its accuracy and usefulness.     

A water retention model accounting for void ratio changes in double porosity clays

This paper presents a constitutive model that predicts the water retention behaviour of compacted clays with evolving bimodal pore size distributions. In line with previous research, the model differentiates between the water present inside the saturated pores of the clay aggregates (the microstructure) and the water present inside the pores between clay aggregates (the macrostructure). A new formulation is then introduced to account for the effect of the macrostructural porosity changes on the retention behaviour of the soil, which results in a consistent evolution of the air-entry value of suction with volumetric deformations. Data from wetting tests on three different active clays (i.e. MX-80 bentonite, FEBEX bentonite, and Boom clay), subjected to distinct mechanical restraints, were used to formulate, calibrate, and validate the proposed model. Results from free swelling tests were also modelled by using both the proposed double porosity model and a published single porosity model, which confirmed the improvement in the predictions of degree of saturation by the present approach. The proposed retention model might be applied, for example, to the simulation of the hydromechanical behaviour of engineered bentonite barriers in underground nuclear waste repositories, where compacted active clays are subjected to changes of both suction and porosity structure under restrained volume conditions.

     

考虑渐近状态特性的超固结土本构模型

基于超固结土的UH模型,结合渐近状态概念,将适用于饱和砂土的渐近状态本构模型推广为考虑渐近状态特性的超固结土本构模型,用于描述土的超固结性和渐近状态特性。通过采用变换应力方法实现了模型的三维化。新提出的本构模型简单、参数较少,与剑桥模型相比,仅增加了一个材料参数,即伏斯列夫面斜率。该模型能够合理反映超固结土的硬化、软化、临界状态、剪缩、剪胀、应力路径依赖性、渐近状态等特性。     

Non‐isothermal plasticity model for cyclic behaviour of soils

On the one hand, it has been observed that liquefaction‐induced shear deformation of soils accumulates in a cycle‐by‐cycle pattern. On the other hand, it is known that heating could induce plastic hardening. This study deals with the constitutive modelling of the effect that heat may have on the cyclic mechanical properties of cohesive soils, a relatively new area of interest in soil mechanics. In this paper, after a presentation of the thermo‐mechanical framework, a non‐isothermal plasticity cyclic model formulation is presented and discussed. The model calibration is described based on data from laboratory sample tests. It includes numerical simulations of triaxial shear tests at various constant temperatures. Then, the model predictions are compared with experimental results and discussed in the final section. Both drained and undrained loading conditions are considered. The proposed constitutive model shows good ability to capture the characteristic features of behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling chemo‐hydro‐mechanical behaviour of unsaturated clays: a feasibility study

Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.     

On incremental non‐linearity in granular media: phenomenological and multi‐scale views (Part I)

On the basis of fundamental constitutive laws such as elasticity, perfect plasticity, and pure viscosity, many elasto‐viscoplastic constitutive relations have been developed since the 1970s through phenomenological approaches. In addition, a few more recent micro‐mechanical models based on multi‐scale approaches are now able to describe the main macroscopic features of the mechanical behaviour of granular media. The purpose of this paper is to compare a phenomenological constitutive relation and a micro‐mechanical model with respect to a basic issue regularly raised about granular assemblies: the incrementally non‐linear character of their behaviour. It is shown that both phenomenological and micro‐mechanical models exhibit an incremental non‐linearity. In addition, the multi‐scale approach reveals that the macroscopic incremental non‐linearity could stem from the change in the regime of local contacts between particles (from plastic regime to elastic regime) in terms of the incremental macroscopic loading direction. Copyright © 2005 John Wiley & Sons, Ltd.     

非饱和土化学-塑性耦合本构行为的数值模拟

基于Hueckel提出的饱和黏土化学-塑性本构模型和Gallipoli提出的非饱和土弹塑性本构模型,提出了一个新的非饱和多孔介质的化学-塑性本构模型,并建立了该模型的隐式积分算法,算法中考虑了化学软化和非饱和吸力的影响。在已有的非饱和多孔介质有限元分析程序平台上进行了程序研发,对孔隙水中化学污染物浓度变化对非饱和土力学行为的影响进行数值模拟,使所研制的程序能够进行岩土工程问题的化学-力学耦合非线性分析。