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.     

Thermodynamic‐based model for coupling temperature‐dependent viscoelastic,viscoplastic, and viscodamage constitutive behavior of asphalt mixtures

Based on the continuum damage mechanics, a general and comprehensive thermodynamic‐based framework for coupling the temperature‐dependent viscoelastic, viscoplastic, and viscodamage behaviors of bituminous materials is presented. This general framework derives systematically Schapery‐type nonlinear viscoelasticity, Perzyna‐type viscoplasticity, and a viscodamage model analogous to the Perzyna‐type viscoplasticity. The resulting constitutive equations are implemented in the well‐known finite element code Abaqus via the user material subroutine UMAT. A systematic procedure for identifying the model parameters is discussed. Finally, the model is validated by comparing the model predictions with a comprehensive set of experimental data on hot mix asphalt that include creep‐recovery, creep, uniaxial constant strain rate, and repeated creep‐recovery tests in both tension and compression over a range of temperatures, stress levels, and strain rates. Comparisons between model predictions and experimental measurements show that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

A three‐phase soil model for simulating stress wave propagation due to blast loading

A three‐phase soil model is proposed to simulate stress wave propagation in soil mass to blast loading. The soil is modelled as a three‐phase mass that includes the solid particles, water and air. It is considered as a structure that the solid particles form a skeleton and their voids are filled with water and air. The equation of state (EOS) of the soil is derived. The elastic–plastic theory is adopted to model the constitutive relation of the soil skeleton. The damage of the soil skeleton is also modelled. The Drucker–Prager strength model including the strain rate effect is used to describe the strength of the soil skeleton. The model is implemented into a hydrocode Autodyn. The recorded results obtained by explosion tests in soil are used to validate the proposed model. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of earthquake‐induced slope deformation using SPH method

This paper presents the development, calibration, and validation of a smoothed particle hydrodynamics (SPH) model for the simulation of seismically induced slope deformation under undrained condition. A constitutive model that combines the isotropic strain softening viscoplasticity and the modified Kondner and Zelasko rule is developed and implemented into SPH formulations. The developed SPH model accounts for the effects of wave propagation in the sliding mass, cyclic nonlinear behavior of soil, and progressive reduction in shear strength during sliding, which are not explicitly considered in various Newmark‐type analyses widely used in the current research and practice in geotechnical earthquake engineering. Soil parameters needed for the developed model can be calibrated using typical laboratory shear strength tests, and experimental or empirical shear modulus reduction curve and damping curve. The strain‐rate effects on soil strength are considered. The developed SPH model is validated against a readily available and well‐documented model slope test on a shaking table. The model simulated slope failure mode, acceleration response spectra, and slope deformations are in excellent agreement with the experimental data. It is thus suggested that the developed SPH model may be utilized to reliably simulate earthquake‐induced slope deformations. This paper also indicates that if implemented with appropriate constitutive models, SPH method can be used to model large‐deformation problems with high fidelity. Copyright © 2013 John Wiley & Sons, Ltd.     

Non‐coaxial elasto‐plasticity model and bifurcation prediction of shear banding in sands

Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 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.     

Formulation of a three‐dimensional rate‐dependent constitutive model for chalk and porous rocks

Rate‐dependent behaviour of chalk and other porous rocks has undergone widespread study in geomechanics due to its implications on the performance of engineering structures. We present a rate‐dependent constitutive model for chalk and other porous rocks with several new features. The model formulation is based on a viscoplastic rate‐lines approach in which the axial strain rate depends on the proximity of the stress point to an elliptical reference surface. A non‐associated viscoplastic potential surface and an axial scaling algorithm are used to determine the viscoplastic strain components. The model predicts that axial yields stress varies as a power function of applied axial strain rate, as shown by published laboratory data. Comparisons with published experimental data indicate that the model is capable of reproducing observed rate‐dependent behaviour of chalk under a variety of loading conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

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

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.     

Effects of transport of pore water and material heterogeneity on strain localization of fluid‐saturated gradient‐dependent viscoplastic geomaterial

The purpose of the present paper is to clarify the effects of permeability and initial heterogeneity on the strain localization of fluid‐saturated cohesive soil modelled by a strain gradient‐dependent poro‐viscoplastic constitutive model. The effects of permeability and gradient parameters on the growth rate of the fluctuation were obtained by a linear instability analysis. Deformation behaviour of clay specimens modelled as a viscoplastic model with a second order strain gradient during shear was numerically analysed by a soil–water coupled FEM under both globally undrained and partially drained conditions. It was found that the deformation pattern and the stress–strain curve greatly depend on the permeability, the drainage conditions and the initial non‐homogeneous properties. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

A non‐coaxial constitutive model for sand deformation under rotation of principal stress axes

A constitutive model for the simulation of non‐coaxiality, an aspect of anisotropic behavior of sand subjected to the rotation of the principal stress axes, is presented in this paper. Experimental studies have shown that non‐coaxiality or non‐coincidence of principal plastic strain increments with principal stress axes under loadings involving the rotation of principal stress axes may be considerable. Besides, the rotation of the principal stress axes results in dramatic effects on stiffness and dilatant behavior of sand. Therefore, the consequences of principal stress axes rotation on deformational behavior, dilatancy and soil stiffness must be taken into account in theoretical and practical problems. To this aim, the following steps are taken: (1) A general relationship for flow direction with respect to possibility of non‐coaxial flow is developed. Moreover, special circumstances linking non‐coaxiality to instantaneous interaction between loading and soil fabric are proposed. (2) Proposing novel expressions for plastic modulus and dilatancy function, the model is enforced to provide realistic simulations when sand is subjected to the rotation of the principal stress axes. Finally, with numerous examples and comparisons, the model capabilities are shown under various stress paths and drainage conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

On finite volume method implementation of poro‐elasto‐plasticity soil model

Accurate prediction of the interactions between the nonlinear soil skeleton and the pore fluid under loading plays a vital role in many geotechnical applications. It is therefore important to develop a numerical method that can effectively capture this nonlinear soil‐pore fluid coupling effect. This paper presents the implementation of a new finite volume method code of poro‐elasto‐plasticity soil model. The model is formulated on the basis of Biot's consolidation theory and combined with a perfect plasticity Mohr‐Coulomb constitutive relation. The governing equation system is discretized in a segregated manner, namely, those conventional linear and uncoupled terms are treated implicitly, while those nonlinear and coupled terms are treated explicitly by using any available values from previous time or iteration step. The implicit–explicit discretization leads to a linearized and decoupled algebraic system, which is solved using the fixed‐point iteration method. Upon the convergence of the iterative method, fully nonlinear coupled solutions are obtained. Also explored in this paper is the special way of treating traction boundary in finite volume method compared with FEM. Finally, three numerical test cases are simulated to verify the implementation procedure. It is shown in the simulation results that the implemented solver is capable of and efficient at predicting reasonable soil responses with pore pressure coupling under different loading situations. Copyright © 2015 John Wiley & Sons, Ltd.     

Limitations of suction‐controlled triaxial tests in the characterization of unsaturated soils

Behavior of unsaturated soils is influenced by many factors, and the influences of these factors are usually coupled together. Suction‐controlled triaxial (SCTX) tests are considered to allow researchers to investigate influences of individual variables on unsaturated soils under specified stress path with controls of stresses, pore water, and air pressures. In the past 50 years, SCTX testing method has been established as a standard approach to characterize constitutive behavior of unsaturated soils. Most important concepts for modern unsaturated soil mechanics were developed upon results from the SCTX tests. Among these, one of the most important contributions in the constitutive modeling of elasto‐plastic behavior for unsaturated soils is the Barcelona basic model (BBM) proposed by Alonso et al. in 1990. The BBM successfully explained many features of unsaturated soils and received extensive acceptance. However, the SCTX tests are designed based upon the divide‐and‐conquer approach in which an implicit assumption is used: soil behavior is stress‐path independent. However, it is well‐established that unsaturated soil behavior is elasto‐plastic and stress‐path dependent. It is found that the SCTX tests in fact cannot control the stress path of an unsaturated soil during loading. This incapability, in combination with complicated loading/collapse behavior of unsaturated soils, makes the SCTX tests for characterizing unsaturated soil questionable. This paper discusses the limitations of the SCTX tests in the characterization of unsaturated soils. A possible solution to the problem was proposed based on a newly developed modified state surface approach. The discussions are limited for isotropic conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

A hypoplastic macroelement model for shallow foundations under monotonic and cyclic loads

In this work, a macroelement for shallow foundations on sands has been developed based on the theory of hypoplasticity. The incrementally nonlinear constitutive equations of the macromodel are defined in terms of generalized forces and displacements and are constructed based on the general approach proposed by Niemunis (Extended Hypoplastic Models for Soils. Habilitation Thesis, Bochum University, 2002). A suitable vectorial internal variable—mimicking the concept of intergranular strain introduced by Niemunis and Herle (Mech Cohes Frict Mater 2:279–299, 1997) for continuum hypoplasticity—is employed to provide the model sufficient memory of past displacement history to be able to reproduce the observed behavior under cyclic loading paths. As compared to similar macroelements formulated within the framework of the theory of elastoplasticity, the proposed approach has the advantage of a much simpler mathematical structure, which allows a straightforward implementation in existing structural analysis FE codes. The model performance has then been evaluated by comparing the model predictions with available experimental results from a series of small-scale model tests reported by Nova and Montrasio (Géotechnique 41:243–256, 1991). Overall, the model captures reasonably well the observed response under nonproportional, complex loading paths. A series of simulations reproducing the tests performed by di Prisco et al. (Shallow footing under cyclic loading: experimental behaviour and constitutive modelling. Patron, Bologna, 2003) has also demonstrated the potential of the proposed model to simulate the observed behavior of footings under cyclic loading paths, at least from a qualitative point of view.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.