Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. 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.     

Evaluation of a constitutive model for clays and sands: Part I – sand behaviour

This paper evaluates the performance of a generalized effective stress soil model for predicting the rate independent behaviour of freshly deposited sands, while a companion paper describes model capabilities for clays and silts. Most material parameters can be obtained from standard laboratory data, including hydrostatic or one‐dimensional compression, drained and undrained triaxial shear testing. A compilation of data on compression behaviour allows for estimation of compression parameters when this type of data is not available. Extensive comparisons of model predictions with measured data from undrained triaxial shear tests shows that the model gives excellent predictions of the transition from dilative to contractive shear response as the confining pressure and/or the initial formation void ratio increases. A parametric study of drained response shows that the model describes realistically the variation of peak friction angle and dilation rate as a function of confining pressure and density when compared with an empirical correlation valid for many sands. The proposed formulation predicts a unique critical state locus for both drained and undrained triaxial testing which is non‐linear over the entire range of stresses and is in excellent agreement with recent experimental data. Overall, the model provides excellent predictions of the stress–strain–strength relationships over a wide range of confining pressures and formation densities. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical studies of shear banding in interface shear tests using a new strain calculation method

Strain localization is closely associated with the stress–strain behaviour of an interphase system subject to quasi‐static direct interface shear, especially after peak stress state is reached. This behaviour is important because it is closely related to deformations experienced by geotechnical composite structures. This paper presents a study using two‐dimensional discrete element method (DEM) simulations on the strain localization of an idealized interphase system composed of densely packed spherical particles in contact with rough manufactured surfaces. The manufactured surface is made up of regular or irregular triangular asperities with varying slopes. A new simple method of strain calculation is used in this study to generate strain field inside a simulated direct interface shear box. This method accounts for particle rotation and captures strain localization features at high resolution. Results show that strain localization begins with the onset of non‐linear stress–strain behaviour. A distinct but discontinuous shear band emerges above the rough surface just before the peak stress state, which becomes more expansive and coherent with post‐peak strain softening. It is found that the shear bands developed by surfaces with smaller roughness are much thinner than those developed by surfaces with greater roughness. The maximum thickness of the intense shear zone is observed to be about 8–10 median particle diameters. The shear band orientations, which are mainly dominated by the rough boundary surface, are parallel with the zero extension direction, which are horizontally oriented. Published in 2007 by John Wiley & Sons, Ltd.     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Elasto‐plastic model for cement‐treated sand

This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd.     

Influence of Lode angle‐dependent failure criteria on shear localization analysis in sand

Geotechnical experiments show that Lode angle‐dependent constitutive formulations are appropriate to describe the failure of geomaterials. In the present study, we have adopted one such class of failure criteria along with a versatile constitutive relationship to theoretically analyze the effects of Lode angle on localized shear deformation or shear band formation in loose sand for both drained and undrained conditions. We determine the variation in the possible stress states for shear localization due to the introduction of Lode angle by considering the localized deformation as a bifurcation problem. Further, similar bifurcation analysis is performed for the stress states along a specific loading path, namely, plane strain compression at the constitutive level. In addition, the plane strain compression tests have been simulated as a boundary value finite element problem to see how Lode angle affects the post‐localization response. Results show that the inclusion of a Lode angle parameter within the failure criterion has considerable effects on the onset, plastic strain, and propagation of shear localization in loose sand specimens. For drained condition, we notice early inception of shear localization and multiple band formation when the Lode angle‐dependent failure criterion is used. Undrained localization characteristics, however, found to be independent of Lode angle consideration.     

A laboratory study of anisotropic geomaterials incorporating recent micromechanical understanding

This paper presents an experimental investigation revisiting the anisotropic stress–strain–strength behaviour of geomaterials in drained monotonic shear using hollow cylinder apparatus. The test programme has been designed to cover the effect of material anisotropy, preshearing, material density and intermediate principal stress on the behaviour of Leighton Buzzard sand. Experiments have also been performed on glass beads to understand the effect of particle shape. This paper explains phenomenological observations based on recently acquired understanding in micromechanics, with attention focused on strength anisotropy and deformation non-coaxiality, i.e. non-coincidence between the principal stress direction and the principal strain rate direction. The test results demonstrate that the effects of initial anisotropy produced during sample preparation are significant. The stress–strain–strength behaviour of the specimen shows strong dependence on the principal stress direction. Preloading history, material density and particle shape are also found to be influential. In particular, it was found that non-coaxiality is more significant in presheared specimens. The observations on the strength anisotropy and deformation non-coaxiality were explained based on the stress–force–fabric relationship. It was observed that intermediate principal stress parameter b(b = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) has a significant effect on the non-coaxiality of sand. The lower the b-value, the higher the degree of non-coaxiality is induced. Visual inspection of shear band formed at the end of HCA testing has also been presented. The inclinations of the shear bands at different loading directions can be predicted well by taking account of the relative direction of the mobilized planes to the bedding plane.     

Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.     

Strain softening of K0-consolidated Changi sand under plane-strain conditions

Experimental results are presented in this paper to study the strain softening behaviour of a marine dredged sand under plane-strain conditions. K₀ consolidated drained and undrained tests were conducted using a new plane-strain apparatus to characterize the strain softening behaviour of the sand under plane-strain conditions. For medium dense specimens, strain softening and shear bands were observed to occur under both drained and undrained conditions. For very loose specimens, no shear bands were observed and critical states were reached within the homogeneous deformation region in both drained and undrained tests. Strain softening was observed to occur at small strain for very loose specimens under undrained conditions. Two types of strain softening, the homogenous softening and banding softening, were identified and the conditions for strain softening were established. The results obtained from this study were compared with the studies by Han and Vardoulakis (Géotechnique 41(1):49–78, 1991), Finno et al. (J Geotech Eng ASCE 122(6):462–473, 1996, Géotechnique 47(1):149–165, 1997) and Mokni and Desrues (Mech Cohes-Frict Mat 4:419–441, 1998).     

Ground response curves for rock masses exhibiting strain‐softening behaviour

A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Visco‐elastic load transfer models for axially loaded piles

Viscoelastic or creep behaviour can have a significant influence on the load transfer (t–z) response at the pile–soil interface, and thus on the pile load settlement relationship. Many experimental and theoretical models for pile load transfer behaviour have been presented. However, none of these has led to a closed‐form expression which captures both non‐linearity and viscoelastic behaviour of the soil. In this paper, non‐linear viscoelastic shaft and base load transfer (t–z) models are presented, based on integration of a generalized viscoelastic stress–strain model for the soil. The resulting shaft model is verified through published field and laboratory test data. With these models, the previous closed‐form solutions evolved for a pile in a non‐homogeneous media have been readily extended to account for visco‐elastic response. For 1‐step loading case, the closed‐form predictions have been verified extensively with previous more rigorous numerical analysis, and with the new GASPILE program analysis. Parametric studies on two kinds of commonly encountered loading: step loading, ramp (linear increase followed by sustained) loading have been performed. Two examples of the prediction of the effects of creep on the load settlement relationship by the solutions and the program GASPILE, have been presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of boundary conditions and partial drainage on cyclic simple shear test results—a numerical study

Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density D_r=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd.     

Shear Strength Response of a Geotextile-Reinforced Chlef Sand: A Laboratory Study

During the last earthquake that occurred in Chlef (El Asnam 1980, Algeria), a significant decrease in the shear strength has caused major damages to several civil and hydraulic structures (earth dams, embankments, bridges, slopes and buildings), especially for the saturated sandy soil of the areas near Chlef valley. This paper presents a laboratory study of drained compression triaxial tests conducted on sandy soil reinforced with horizontal layers of geotextile, in order to study the influence of geotextile layer characteristics both on shear stress–strain and on volumetric change–strain. Tests were carried out on medium and dense sand. The experimental programme includes some drained compression tests performed on reinforced sand samples, for different values of the geotextile layers number (N _g), of confining pressure (\( \sigma_{\text{c}}^{\prime } \)) and relative density (D _r). The test results have shown that the contribution of the geotextile at low values of the axial strain (ε ₁) is negligible, for higher values of (ε ₁); geotextile induces a quasi-linear increase in the deviator stress (q) and leads to an increase in the volume contractiveness within the reinforced samples. A negligible influence of geotextile layers number (N _g) on the stress–strain behaviour and the volumetric change has been shown, when normalized with N _g. The results indicate that the contribution of geotextile to the stress–strain mobilization increases with increasing confining pressure, while its contribution to the volume contraction decreases with the increase in the confining pressure.     

A Bounding Surface Plasticity Model for Intact Rock Exhibiting Size-Dependent Behaviour

A new constitutive model for intact rock is presented recognising that rock strength, stiffness and stress–strain behaviour are affected by the size of the rock being subjected to loading. The model is formulated using bounding surface plasticity theory. It is validated against a new and extensive set of unconfined compression and triaxial compression test results for Gosford sandstone. The samples tested had diameters ranging from 19 to 145 mm and length-to-diameter ratios of 2. The model captures the continuous nonlinear stress–strain behaviour from initial loading, through peak strength to large shear strains, including transition from brittle to ductile behaviour. The size dependency was accounted for through a unified size effect law applied to the unconfined compressive strength—a key model input parameter. The unconfined compressive strength increases with sample size before peaking and then decreasing with further increasing sample size. Inside the constitutive model two hardening laws act simultaneously, each driven by plastic shear strains. The elasticity is stress level dependent. Simple linear loading and bounding surfaces are adopted, defined using the Mohr–Coulomb criterion, along with a non-associated flow rule. The model simulates well the stress–strain behaviour of Gosford sandstone at confining pressures ranging from 0 to 30 MPa for the variety of sample sizes considered.     

Kinematic modelling of shear band localization using discrete finite elements

Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 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.     

On the volumetric deformation of reconstituted soils

This paper reviews the phenomenon of volumetric hardening, which is a common feature of the mechanical behaviour of many geo‐materials. Three different material idealizations have been proposed to describe this hardening, and the paper contains the corresponding mathematical formulation. These idealizations vary in their complexity and hence their ability to capture different aspects of real material behaviour. Any of the three postulates can be implemented into most constitutive models. As a demonstration of their capabilities, the postulates have been implemented into the well‐known modified Cam Clay model, and computations are made with the resulting new constitutive models. It is seen that the new models can successfully represent important features of soil behaviour such as plastic yielding associated with loading inside the current virgin yield surface, the loosening or densifying of granular soils caused by shearing, and the accumulation of both volumetric and distortional deformation caused by repeated drained loading over a large number of cycles. Copyright © 2000 John Wiley & Sons, Ltd.