Kinematic hardening plasticity formulation of small strain behaviour of soils

An objective of this paper is to demonstrate that the small strain model developed by the authors can be incorporated into the conventional kinematic hardening plasticity framework to predict pre‐failure defor mations. The constitutive model described in this paper is constituted by three elliptical yield surfaces in triaxial stress space. Two inner surfaces are rotated ellipses of the same shape, representing the boundaries of the linear elastic and small strain regions, while the third surface is the modified Cam clay large‐scale yield surface. Within the linear elastic region, the soil behaviour is elastic with cross‐coupling between the shear and volumetric stress–strain components. Within the small strain region, the soil behaviour is elasto‐plastic, described by the kinematic hardening rule with an infinite number of loading surfaces defined by the incremental energy criterion. Within the large‐scale yield surface, the soil behaviour is elasto‐plastic, described by kinematic and isotropic hardening of the small strain region boundary. Since the yield surfaces have different shapes, the uniqueness of the plastic loading condition imposes a restriction on the ratio between their semi‐diameters. The model requires 12 parameters, which can be determined from a single consolidated undrained triaxial compression test. Copyright © 2000 John Wiley & Sons, Ltd.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

A multi‐mechanism constitutive model for plastic adaption under cyclic loading

As is well known, granular soils under cyclic loading dissipate a large amount of energy and accumulate large irreversible strains. Usually, with time, this second effect reduces and the accumulation rate decreases with the number of cycles until obtaining a sort of ideal stationary cyclic state at which ratcheting disappears. In this paper, only this ideal state is taken into consideration and simulated by means of a multi‐mechanism constitutive model for plastic adaptation. For this purpose, the concept of cycle is discussed, many different categories of cyclic stress/strain paths are considered and some theoretical issues concerning both the flow and the strain‐hardening rules are tackled. Even though the paper focuses on soil behaviour, the conclusions can be extended to all materials exhibiting ratcheting due to volumetric behaviour.Copyright © 2013 John Wiley & Sons, Ltd.     

Directional response of a reconstituted fine‐grained soil—Part II: performance of different constitutive models

In this paper, the performance of different advanced constitutive models for soils is evaluated with respect to the experimentally observed behaviour of a soft reconstituted clay subject to a wide range of loading directions, see (presented in the companion paper). The models considered include a three‐surface kinematic hardening elastoplastic model; the CLoE hypoplastic model; a recently proposed K‐hypoplastic model for clays, and an enhanced version of the same model incorporating the concept of intergranular strain. A clear qualitative picture of the relative performance of the different models as a function of the loading direction is obtained by means of the incremental strain response envelopes. The definition of suitable error measures allows to obtain further quantitative information in this respect. For the particular initial conditions and loading programme considered in this study, the kinematic hardening and the enhanced K‐hypoplastic models appear to provide the best performance overall. Copyright © 2006 John Wiley & Sons, Ltd.     

Fibre reinforced sands: from experiments to modelling and beyond

Based on hypotheses derived directly from experimental observations of the triaxial behaviour, a constitutive model for fibre reinforced sands is built in this paper. Both the sand matrix and the fibres obey their own constitutive law, whereas their contributions are superimposed using a volumetric homogenization procedure. The Severn‐Trent sand model, which combines well‐known concepts such as critical state theory, Mohr‐Coulomb like strength criterion, bounding surface plasticity and kinematic hardening, is adopted for the sand matrix. Although the fibres are treated as discrete forces with defined orientation, an equivalent continuum stress for the fibre phase is derived to allow the superposition of effects of sand and fibres. The fibres are considered as purely tensile elements following a linear elastic constitutive rule. The strain in the fibres is expressed as a fraction of the strain in the reinforced sample so that imperfect bonding is assumed at the sand‐fibre interface. Only those fibres oriented within the tensile strain domain of the sample can mobilize tensile stress—the orientation of fibres is one of the key ingredients to capture the anisotropic behaviour of fibre reinforced soil that is observed for triaxial compression and extension loading. A further mechanism of partition of the volume of voids between the fibres and the sand matrix is introduced and shown to be fundamental for the simulation of the volumetric behaviour of fibre‐reinforced soils. Copyright © 2012 John Wiley & Sons, Ltd.     

A theoretical framework for constructing elastic/plastic constitutive models of triaxial tests

Modern ideas of thermomechanics are used to develop families of models describing the elastic/plastic behaviour of cohesionless soils deforming under triaxial conditions. Once the form of the free energy and dissipation potential functions have been specified, the corresponding yield loci, flow rules, isotropic and kinematic hardening rules as well as the elasticity law are deduced in a systematic manner. The families contain the classical linear frictional (Coulomb type) models and the classical critical state models as special cases. The generalized models discussed here include non‐associated flow rules, shear as well as volumetric hardening, anisotropic responses and rotational yield loci. The various parameters needed to describe the models can be interpreted in terms of ratio of the plastic work, which is dissipated, to that which is stored. Non‐associated behaviour is found to occur whenever this division between dissipated and stored work is not equal. Micro‐level interpretations of stored plastic work are discussed. The models automatically satisfy the laws of thermodynamics, and there is no need to invoke any stability postulates. Some classical forms of the peak‐strength/dilatancy relationship are established theoretically. Some representative drained and undrained paths are computed. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonlinear SEM numerical analyses of dry dense sand specimens under rapid and dynamic loading

The paper mainly concerns the mechanical response of 2D dry dense sand specimens under shock loading. The problem is numerically analysed by means of a SEM dynamic code, within which an already conceived non‐local viscoplastic constitutive model characterized by a non‐associated flow rule and by an anisotropic strain hardening has been implemented. In particular the strain localization and time dependency of the material mechanical response are taken into consideration. Both rapid/static loading and dynamic histories are numerically simulated. In the first case, the time dependency of the material mechanical response can be captured by neglecting inertial effects, while in the second one the two factors are superimposed and the propagation of the stress waves within the specimen is discussed. The interest of these analyses derives from the fact that the diffusion phenomenon takes place within a specimen already localized. Copyright © 2006 John Wiley & Sons, Ltd.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

基于Drucker-Prager准则的岩石弹塑性损伤本构模型研究

大多数岩石材料软化本构模型在硬化函数中引入塑性内变量来表示材料的硬化/软化性质,但并不能反映岩石微裂隙损伤对材料力学性能的影响及单轴拉伸和压缩所表现的初始屈服强度f0与屈服极限fu的差异。基于D-P准则同时考虑塑性软化及损伤软化,建立岩石类材料的弹塑性本构关系及其数值算法。塑性屈服函数采用Borja等的应力张量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;由于岩石损伤软化是微裂隙扩展所导致的体积膨胀引起的,因此,提出用体积应变表征岩石损伤变量的演化,并用回映隐式积分算法编制了岩石的弹塑性损伤本构程序。对单轴压缩及拉伸荷载作用下的岩石材料试验进行数值模拟,结果表明,所提出的岩石弹塑性损伤本构模型可以较好地符合岩石材料的力学特性。     

A general approach to model interfaces using existing soil constitutive models application to hypoplasticity

The modelling of interfaces is important for the holistic simulation of geotechnical structures (e.g. piles, tunnels and geogrids). For this reason, advanced constitutive interface models and numerical techniques are needed. There are few user-friendly models, and these are rarely implemented. In this paper, a new approach for advanced interface models is proposed. This is based on the assumption that the fully rough interface can be modelled considering simple-shear behaviour at the interface. A 3D soil model is used as a constitutive driver for a frictional subroutine. This minimises the effort required, and advanced interface models are available with less effort. Two different hypoplastic models are used with the new approach. The approach was verified for several aspects (e.g. mesh size dependence), and the volumetric behaviour was studied. The user-friendliness and absence of additional parameters led to more realistic simulation results. The proposed method can be extended to other modelling techniques and will improve the modelling of contacts in soil-structure interaction analysis.     

Constitutive model for cyclic behaviour of cohesionless sands

This paper presents a constitutive model for describing the stress-strain response of sands under cyclic loading. The model, formulated using the critical state theory within the bounding surface plasticity framework, is an upgraded version of an existing model developed for monotonic behaviour of cohesionless sands. With modification of the hardening law, plastic volumetric strain increment and unloading plastic modulus, the original model was modified to simulate cyclic loading. The proposed model was validated against triaxial cyclic loading tests for Fuji River sand, Toyoura sand and Nigata sand. Comparison between the measured and predicted results suggests that the proposed modified model can capture the main features of cohesionless sands under drained and undrained cyclic loading.     

Probabilistic yielding and cyclic behavior of geomaterials

In this paper, the novel concept of probabilistic yielding is used for 1‐D cyclic simulation of the constitutive behavior of geomaterials. Fokker–Planck–Kolmogorov equation‐based probabilistic elastic–plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening‐type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic‐perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response. Copyright © 2010 John Wiley & Sons, Ltd.     

基于旋转运动硬化理论的修正剑桥模型的改进

修正剑桥模型是最早建立和得到广泛承认的经典土体弹塑性模型之一,但不能模拟应力路径转折时土体的应力-应变特性以及应力引起的各向异性。将旋转运动硬化理论引入到剑桥模型中,给出了椭圆屈服面的旋转运动硬化机制,在不增加任何模型参数的情况下,把等向硬化的修正剑桥模型扩展为旋转运动硬化模型。扩展的新模型既保留了单调加载时的等向硬化的特性,又能反映应力路径转折时土体的本构特性与应力诱发的各向异性,初步验证了模型的有效性。     

Visualization of material stiffness in geomechanics analysis

This paper presents novel visualization techniques to simplify representation of the fourth‐order material stiffness tensor as a set of three‐dimensional geometric objects. Stiffness visualization aids in understanding the complex stiffness characteristics of highly non‐linear constitutive models including modelled material anisotropy and loading path dependent stiffness variation. Stiffness visualization is relevant for understanding the relationship of material stiffness to global behaviour in the analysis of a boundary value problem. The spherical pulse stiffness visualization method, developed in the acoustics field, is extended to visualize stiffness of geomaterials using three three‐dimensional objects. This method is limited to relatively simple constitutive models with symmetric stiffness matrices insensitive to loading magnitude and direction. A strain dependent stiffness visualization method is developed that allows the examination of material stiffness for a range of loading directions and is suitable for highly non‐linear and path dependent material models. The proposed stiffness visualization can be represented as 3‐D, 2‐D and 1‐D objects. The visualization technique is used to represent material stiffness and its evolution during simulated soil laboratory tests and deep excavation construction. Copyright © 2005 John Wiley & Sons, Ltd.     

On generalization of constitutive models from two dimensions to three dimensions

In this paper, a study is made of the generalization of constitutive models for geomaterials from two‐dimensional stress and strain states to three‐dimensional stress and strain states. Existing methods of model generalization are reviewed and their deficiencies are highlighted. A new method is proposed based on geometries of the model imprints on two normal planes. Using the proposed method, various three‐dimensional failure criterions suitable for geomaterials are implemented directly into a two‐dimensional model and the generalized model is identical to its original form for the axially symmetric condition. To demonstrate the application of the proposed method, the Modified Cam Clay model is extended using the Matsuoka–Nakai failure criterion. Simulations of soil behaviour for loading in the principal stress space are presented and analysed. Copyright © 2008 John Wiley & Sons, Ltd.     

Liquefaction and cyclic mobility model for saturated granular media

A new constitutive law for the behaviour of undrained sand subjected to dynamic loading is presented. The proposed model works for small and large strain ranges and incorporates contractive and dilative properties of the sand into the unified numerical scheme. These features allow to correctly predict liquefaction and cyclic mobility phenomena for different initial relative densities of the soil. The model has been calibrated as an element test, by using cyclic simple shear data reported in the literature. For the contractive sand behaviour a well‐known endochronic densification model has been used, whereas a plastic model with a new non‐associative flow rule is applied when the sand tends to dilate. Both dilatancy and flow rule are based on a new state parameter, associated to the stiffness degradation of the material as the shaking goes on. Also, the function that represents the rearrangement memory of the soil takes a zero value when the material dilates, in order to easily model the change in the internal structure. Proceeding along this kind of approach, liquefaction and cyclic mobility are modelled with the same constitutive law, within the framework of a bi‐dimensional FEM coupled algorithm developed in the paper. For calibration purposes, the behaviour of the soil in a cyclic simple shear test has been simulated, in order to estimate the influence of permeability, frequency of loading, and homogeneity of the shear stress field on the laboratory data. Copyright © 2006 John Wiley & Sons, Ltd.     

NGI‐ADP: Anisotropic shear strength model for clay

Many geotechnical problems involve undrained behavior of clay and the capacity in undrained loading. Most constitutive models used today are effective stress based and only indirectly obtain values for the undrained shear strength. To match the design profiles of undrained shear strengths, in active (A), direct simple shear (D) and passive (P) modes of loading are complicated. This paper presents the elastoplastic constitutive model NGI‐ADP which is based on the undrained shear strength approach with direct input of shear strengths. Consequently, exact match with design undrained shear strengths profiles is obtained and the well‐known anisotropy of undrained shear strength and stiffness is accounted for in the constitutive model. A non‐linear stress path‐dependent hardening relationship is used, defined from direct input of failure strains in the three directions of shearing represented by triaxial compression, direct simple shear and triaxial extension. With its clear input parameters the model has significant advantages for design analysis of undrained problems. The constitutive model is implemented, into finite element codes, with an implicit integration scheme. Its performance is demonstrated by a finite element analysis of a bearing capacity problem. Copyright © 2011 John Wiley & Sons, Ltd.     

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.