A viscoplastic subloading soil model for rate‐dependent cyclic anisotropic structured behaviour

This paper presents a new purely viscoplastic soil model based on the subloading surface concept with a mobile centre of homothety, enabling the occurrence of viscoplastic strains inside the yield surface and avoiding the abrupt change in stiffness of the traditional overstress viscoplastic models. This is required for overconsolidated soils. The model is formulated to reproduce the soil rate‐dependent behaviour under cyclic loading (changes in loading direction) and incorporates both initial and induced anisotropy, as well as destructuring. The model shows good qualitative response to some imposed three‐dimensional stress paths under quasi‐inviscid (elastoplastic) behaviour. Some of the main time‐dependent aspects of soil behaviour that the model is capable of reproducing were also illustrated. The capability of the model to adequately reproduce the results from an undrained triaxial test performed on stiff overconsolidated clays from the Lisbon region (Formação de Benfica), with an unloading–reloading deviatoric stress cycle at constant mean stress, that incorporates a series of staggered fast loading and creep stages, was evaluated. The model was able to reproduce well the main observed aspects of the time‐dependent stress–strain response and pore pressure evolution of a stiff overconsolidated clay under complex loading. The revised and generalised viscoplastic subloading surface concept is viable and can be applied to a consistent extension to viscoplasticity, including in the interior of the yield surface, of existing elastoplastic models formulated for soils and other materials. Copyright © 2016 John Wiley & Sons, Ltd.     

Anisotropic viscoplastic modeling of rate‐dependent behavior of clay

The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K₀‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluation of an anisotropic elastoplastic–viscoplastic bounding surface model for clays

An anisotropic time-dependent bounding surface model for clays is developed by generalizing a previous time-independent model that adopts a flexible bounding surface. It is based on the framework for coupled elastoplasticity–viscoplasticity for clays and Perzyna’s overstress theory. Three viscoplastic parameters were introduced and explained in detail. The model was validated against undrained creep tests for both isotropically and anisotropically consolidated clays, undrained and drained stress relaxation tests on some undisturbed clays, and undrained triaxial tests with varying strain rates on natural Hong Kong marine deposit clay. The general agreement between the model simulations and test results was satisfactory. The varying effects of lower-level parameters were discussed on the undrained multistage stress relaxation response for normally consolidated soils which had been ignored in literature. The flexibility of the model in capturing the shear strengths, which is the unique feature of the current model, was shown in the simulations of time-dependent triaxial tests on Taipei silty clay. All the simulations show that the proposed model is a relatively practical model considering both anisotropy and time dependency of clays.     

Nonstationary flow surface theory for modeling the viscoplastic behaviors of soils

This paper presents a three-dimensional elastic viscoplastic model that can describe the time-dependent behaviors of soft clays. The constitutive model is formulated based on the nonstationary flow surface theory and incorporates new developments, including (i) an improved definition of the nonstationary flow surface that is capable of capturing the stress–strain behaviors under different loading paths, (ii) a unique stress–strain—viscoplastic-strain-rate equation that is able to explicitly describe the nonstationary flow surface, and (iii) a final stable state concept that identifies the final equilibrium state at the end of creep and stress relaxation, which is also used to simplify the loading criteria. The consistency condition is validated for the proposed model, and the viscoplastic multipliers are calculated by solving the consistency equations. The model performance is investigated and validated via simulation of both oedometer and triaxial tests. The numerical results demonstrate that the proposed model is able to reproduce the main viscoplastic behaviors of soils, including creep, undrained creep rupture, stress relaxation, rate effect and accumulated effect.     

各向异性弹黏塑性模型在ABAQUS中的研发

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明：(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。     

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 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.     

A simple hypoplastic constitutive model for sand

The paper presents a hypoplastic constitutive model for the three-dimensional non-linear stress–strain and dilatant volume change behaviour of sand. The model is developed without recourse to the concept in elastoplasticity theory such as yield surface, plastic potential and decomposition into elastic and plastic parts. Benefited from the non-linear tensorial functions available from the representation theorem the model possesses simple mathematical formulation and contains only four material parameters, which can be easily identified with triaxial compression tests. Comparison of the predictions with the experimental results shows that the model is capable of capturing the salient behaviour of sand under monotonic loading and is applicable to both drained and undrained conditions.     

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.     

循环荷载下非饱和结构性黏土的弹塑性双面模型

非饱和黏土的结构性能够显著影响其力学特性。基于非饱和土经典模型BBM（Barcelona basic model）和一种可描述循环塑性的硬化法则,引入体积破损率的作为标准土体结构破损的参数,建立了一个描述常吸力下非饱和结构性黏土静态及动态力学特性的弹塑性双面模型。模型在应力空间中包含与重塑非饱和土屈服面几何相似的结构性边界面和加载面,采用径向映射法则和可移动的记忆中心原理,通过结构性边界面和加载面在应力空间中的演化来反映循环加载过程中非饱和结构性黏土的循环塑性特征和结构损伤过程。通过与相关非饱和黏土控制吸力试验数据的比较,表明该模型能够较好地反映静态加载下非饱和结构性黏土的力学特性,而模型预测的循环荷载下的应力?应变特征也具有一定的合理性。     

A critical state model for overconsolidated structured clays

This paper presents a generalised critical state model with the bounding surface theory for simulating the stress–strain behaviour of overconsolidated structured clays. The model is formulated based on the framework of the Structured Cam Clay (SCC) model and is designated as the Modified Structured Cam Clay with Bounding Surface Theory (MSCC-B) model. The hardening and destructuring processes for structured clays in the overconsolidated state can be described by the proposed model. The image stress point defined by the radial mapping rule is used to determine the plastic hardening modulus, which varies along loading paths. A new proposed parameter h, which depends on the material characteristics, is introduced into the plastic hardening modulus equation to take the soil behaviour into account in the overconsolidated state. The MSCC-B model is finally evaluated in light of the model performance by comparisons with the measured data of both naturally and artificially structured clays under compression and shearing tests. From the comparisons, it is found that the MSCC-B model gives reasonable good simulations of mechanical response of structured clays in both drained and undrained conditions. With its simplicity and performance, the MSCC-B model is regarded as a practical geotechnical model for implementation in numerical analysis.     

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.     

Undrained Triaxial Behavior of Cement Treated Marine Clay

Clays treated with lower cement contents often exhibit behaviour similar to stiff clays with planar failure surface under triaxial compression. In the present work the behaviour of a marine clay treated with 5 % cement, subjected to undrained triaxial compression tests is studied. The pre-consolidation pressure of the cemented clay due to the cementation bonding is observed to be very high. It is attempted to model the behaviour of cement treated clay using a bounding surface plasticity formulation as the plastic behaviour of the cemented clays within the yield surface has to be considered. The effect of cementation is included in the model as the pre-consolidation pressure obtained from consolidation tests. The tensile strength due to cementation bonds is included in the equation of the bounding surface. Simulations of the undrained triaxial compression tests on cemented clays are carried out and the results are validated with the experimental results.     

On double shearing in frictional materials

This paper evaluates the mechanical behaviour of yielding frictional geomaterials. The general Double Shearing model describes this behaviour. Non‐coaxiality of stress and plastic strain increments for plane strain conditions forms an important part of this model. The model is based on a micro‐mechanical and macro‐mechanical formulation. The stress–dilatancy theory in the model combines the mechanical behaviour on both scales. It is shown that the general Double Shearing formulation comprises other Double Shearing models. These models differ in the relation between the mobilized friction and dilatancy and in non‐coaxiality. In order to describe reversible and irreversible deformations the general Double Shearing model is extended with elasticity. The failure of soil masses is controlled by shear mechanisms. These shear mechanisms are determined by the conditions along the shear band. The shear stress ratio of a shear band depends on the orientation of the stress in the shear band. There is a difference between the peak strength and the residual strength in the shear band. While peak stress depends on strength properties only, the residual strength depends upon the yield conditions and the plastic deformation mechanisms and is generally considerably lower than the maximum strength. It is shown that non‐coaxial models give non‐unique solutions for the shear stress ratio on the shear band. The Double Shearing model is applied to various failure problems of soils such as the direct simple shear test, the biaxial test, infinite slopes, interfaces and for the calculation of the undrained shear strength. Copyright © 2006 John Wiley & Sons, Ltd.     

A viscoplastic model including anisotropic damage for the time dependent behaviour of rock

This paper presents a new constitutive model for the time dependent mechanical behaviour of rock which takes into account both viscoplastic behaviour and evolution of damage with respect to time. This model is built by associating a viscoplastic constitutive law to the damage theory. The main characteristics of this model are the account of a viscoplastic volumetric strain (i.e. contractancy and dilatancy) as well as the anisotropy of damage. The latter is described by a second rank tensor. Using this model, it is possible to predict delayed rupture by determining time to failure, in creep tests for example. The identification of the model parameters is based on experiments such as creep tests, relaxation tests and quasi‐static tests. The physical meaning of these parameters is discussed and comparisons with lab tests are presented. The ability of the model to reproduce the delayed failure observed in tertiary creep is demonstrated as well as the sensitivity of the mechanical response to the rate of loading. The model could be used to simulate the evolution of the excavated damage zone around underground openings. Copyright © 2005 John Wiley & Sons, Ltd.     

Creep behaviour of intact and remoulded fibrous peat

This paper presents the creep behaviour of intact and remoulded specimens of fibrous peat obtained from a field site near Anzac, Alberta, Canada. The creep behaviour was investigated by means of long-term drained and undrained triaxial tests. The development of volumetric, axial, and undrained axial strain and strain rate during drained and undrained creep tests under variable stress conditions is presented. The stress _– strain _– strain rate (p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\)) relationship is found to be unique for different stress and loading durations. The p′–ε _v–\(\dot{\varepsilon }_{\text{v}}\) relationship is analysed and represented by creep isotaches. The applicability of different creep models developed for normally consolidated clay is discussed and applied to define the development of creep strain in fibrous peat under varying isotropic and deviator stresses. The secondary consolidation coefficient for evaluating the volumetric strain rate of peat is found to be applicable with some limits. The drained creep behaviour of remoulded peat specimens differs from the behaviour shown by Shelby tube specimens, whereas the undrained creep behaviour in remoulded and Shelby tube specimens is similar.     

A Strain-Rate Dependent Clay Constitutive Model with Parametric Sensitivity and Uncertainty Quantification

This paper presents a strain-rate dependent plastic constitutive model for clays. Based on the concepts of critical-state soil mechanics and bounding surface plasticity theory, the model reproduces the mechanical response of clays under triaxial and simple shear loading conditions. The model parameters are determined for Boston Blue Clay, London Clay and Kaolin Clay, and the performance of the model in simulating the mechanical response of these clays is demonstrated for low to medium strain rates. The sensitivity of each model parameter is checked by perturbing the calibrated values by ±20 %. Subsequently, a probabilistic analysis using Monte Carlo simulations is performed by treating the model parameters as random variables and the impact of the statistics of the parameters on the undrained shear strength is investigated.     

A plasticity model for the mechanical behaviour of anisotropically consolidated clay

An incremental plasticity theory has been developed to describe the mechanical behaviour of anisotropically consolidated clays. The theory removes some of the shortcomings of the existing critical state models by incorporating the effects of the initial anisotropy due to a known depositional stress history and its subsequent alteration during further plastic deformation under a general stress system. From the extensive comparisons with the reported results in the literature, it is established that the model is satisfactory in predicting the various aspects of drained as well as undrained behaviour of K₀-consolidated clays.