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.     

Explicit solutions for the instantaneous undrained contraction of hollow cylinders and spheres in porous elastoplastic medium

In this article we present closed‐form solutions for the undrained variations in stress, pore pressure, deformation and displacement inside hollow cylinders and hollow spheres subjected to uniform mechanical pressure instantaneously applied to their external and internal boundary surfaces. The material is assumed to be a saturated porous medium obeying a Mohr–Coulomb model failure criterion, exhibiting dilatant plastic deformation according to a non‐associated flow rule which accounts for isotropically strain hardening or softening. The instantaneous response of a porous medium submitted to an instantaneous loading is undrained, i.e. without any fluid mass exchange. The short‐term equilibrium problem to be solved is now formally identical to a problem of elastoplasticity where the constitutive equations involve the undrained elastic moduli and particular equivalent plastic parameters. The response of the model is presented (i) for extension and compression undrained triaxial tests, and (ii) for unloading problems of hollow cylinders and spheres through the use of appropriately developed closed‐form solutions. Numerical results are presented for a plastic clay stone with strain hardening and an argilite with strain softening. The effects of plastic dilation, of the strain softening law and also of geometry of the cavity on the behaviour of the porous medium have been underlined. Analytical solutions provide valuable benchmarks enabling various numerical methods in undrained conditions with a finite boundary to be verified. Copyright © 2002 John Wiley & Sons, Ltd.     

Relationship between small and large strain solutions for general cavity expansion problems in elasto-plastic soils

This paper presents a closed-form relationship between small and finite strain cavity expansion solutions. Its derivation is based on the non-linearly elastic–perfectly plastic cylindrical (or spherical) problem considering a general Mohr’s criterion and constant plastic dilatancy. It is shown, however, that it is sufficiently accurate for general expansion problems not obeying plane-strain rotationally (or spherically) symmetric conditions and involving strain-hardening/softening constitutive behaviour. Therefore, this relationship quantifies the error stemming from the computational assumption of small deformations and provides a simple and efficient way of accounting for geometric non-linearity based entirely on conventional computational methods: ‘self-correction’ of small strain analyses results.     

UNDRAINED CAVITY EXPANSIONS IN CRITICAL STATE SOILS

Boundary value problems for hardening/softening soils, such as Cam-Clay, usually require the extensive use of finite element methods. Here analytical and semi-analytical solutions for the undrained expansion of cylindrical and spherical cavities in critical state soils are presented. The strain is finite, the initial cavity radius is arbitrary and the procedure applicable to any isotropically hardening materials. In all cases only simple quadratures are involved, and in the case of the original Cam-Clay a complete analytical solution can be found. In addition to providing models of the behaviour of displacement piles and pressuremeters these results also provide valuable benchmark solutions for verifying various numerical methods.     

Strain‐softening analysis of a spherical cavity

This paper studies the excavation of a spherical cavity subjected to hydrostatic initial stresses in the infinite homogeneous and isotropic rock mass with strain‐softening Mohr–Coulomb (M‐C) and Hoek–Brown (H‐B) behaviors. Numerical solutions of the spherical cavity are obtained and the application to determining stress–strain curve of strain‐softening M‐C and H‐B rock mass is studied. A closed‐form solution for the elastic–brittle–plastic medium is introduced first, and then a numerical procedure that simplifies the strain‐softening process into a series of brittle–plastic ones is presented. The approach is validated against the facts that the strain‐softening process evolves into a brittle–plastic one when the softening slope is very steep, whereas it evolves into an elasto‐plastic one when the softening slope approaches zero. Numerical solutions for the prediction of displacements and stresses around the spherical cavity in the strain‐softening M‐C and H‐B rock mass are presented. On the basis of the analysis of the spherical cavity in strain‐softening rock mass, the stress–strain relationship at an infinitesimal cube around the cavity is obtained and discussed with different evolution laws for the strength parameters considered. Copyright © 2011 John Wiley & Sons, Ltd.     

Stress and pore pressure changes in clay during and after the expansion of a cylindrical cavity

The disturbance of a clay mass, due to either the installation of a driven pile or the expansion of a pressuremeter membrane, is often modelled as a cylindrical cavity expansion. In addition, it is usual (and convenient) to assume that the expansion occurs under conditions of plane strain. For this problem a method of analysis is presented which considers the soil to be a saturated two-phase material with a pore fluid which flows according to Darcy's Law. Non-linearity in material behaviour is permitted as long as the effective stress–strain law can be written in an incremental or rate form. The use of a consolidation analysis allows the changes in effective stress and pore pressure to be determined at any stage during both the cavity expansion and the subsequent period of reconsolidation. Expansions may occur at any prescribed rate, including the very fast (undrained) and the very slow (fully drained) case. The technique is illustrated by considering the expansion of a cavity in two different types of elastoplastic soil. It is shown how these solutions may be used to model the disturbance of the soil due to pile driving.     

Finite strain elastoplastic solutions for the undrained ground response curve in tunnelling

The instantaneous response of saturated low permeability grounds to tunnel excavation is important for deformations and stability close to the tunnel face. It is characterized by zero volume change in combination with the development of excess pore pressures. In tunnelling through poor quality ground under great depth of cover and high in situ pore pressure, heavily squeezing conditions (characterized by very large convergences) may occur soon after excavation. This paper presents exact finite strain analytical solutions for the undrained ground response around cylindrical and spherical openings that are unloaded from uniform and isotropic initial stress states, on the basis of the Modified Cam Clay (MCC) model and the Mohr–Coulomb (MC) model. The solution for a Drucker–Prager material is also given as it requires only a very small modification to the MC solution. The so‐called ground response curve, that is, the relationship between the support pressure and the cavity wall displacement, is derived in closed form for the MC model. The solution for the MCC problem is semi‐analytical in that it uses the trapezium rule for the computation of a definite integral. The influence of the significant parameters of the problem on the predicted deformation behaviour is shown by means of dimensionless charts. Finally, the practical usefulness of the solutions presented is illustrated by applying them to the breccia zones of the planned Gibraltar Strait tunnel – an extreme case of weak, low permeability ground under high pore pressure. The solutions can serve as a trustworthy benchmark for numerical procedures that incorporate material and geometric nonlinearities. Copyright © 2014 John Wiley & Sons, Ltd.     

Advanced constitutive modelling of shotcrete: Model formulation and calibration

Shotcrete is one of the main support elements for tunnels constructed according to the principles of the New Austrian Tunnelling Method (NATM). In this paper, a sophisticated constitutive model for the complex mechanical and time-dependent behaviour of shotcrete is presented, which is formulated within the framework of elasto-plasticity. Two independent yield surfaces govern the material behaviour in compression and tension under multi-axial loading conditions, which is further controlled by non-linear plastic strain hardening and softening rules following uniaxial stress–strain curves. The main novelty of the model is the introduction of normalised hardening and softening parameters which enable a more realistic simulation of the shotcrete loading history. In addition, cracking of the shotcrete is treated within the smeared crack concept by applying a fracture energy approach as a regularisation technique, which eliminates the dependency of the results on the size of the elements in a finite element mesh. Furthermore, the model takes into account the gradual change of the main material properties of the shotcrete due to cement hydration during the first 28 days after spraying. Other important time effects such as creep, shrinkage and hydration temperature induced deformations at early ages are considered in the current model formulation. This paper aims to describe in detail the developed structure of the constitutive model and model calibration.     

Arbitrary Lagrangian Eulerian based finite element analysis of cone penetration in soft clay

This paper presents a numerical model for the analysis of cone penetration in soft clay based on the finite element method. The constitutive behaviour of the soil is modelled by modifying an elastic, perfectly-plastic soil model obeying Von-Mises yield criterion to take into account the strain-softening, rate dependent behaviour of soft clay. Since this is a problem involving large soil deformations, the analysis is carried out using an Arbitrary Lagrangian Eulerian method where the quality of the mesh is preserved during penetration. The variation of cone resistance is examined with various parameters such as rigidity index of the soil, in situ stress anisotropy and roughness at the cone–soil interface, which influence the penetration resistance of the cone. A theoretical correlation has been developed incorporating these parameters and the results have been compared with previous correlations based on the cavity expansion theory, finite element method and strain path method. With the increase in strain-softening, relative brittleness of the soil increases and the penetration resistance is significantly reduced. With the rising strain-rate dependency, penetration resistance increases but this increase is independent of the degree of brittleness of the soil.     

Elastic‐plastic solutions for expanding cavities embedded in two different cohesive‐frictional materials

An analytical solution of cavity expansion in two different concentric regions of soil is developed and investigated in this paper. The cavity is embedded within a soil with finite radial dimension and surrounded by a second soil, which extends to infinity. Large‐strain quasi‐static expansion of both spherical and cylindrical cavities in elastic‐plastic soils is considered. A non‐associated Mohr–Coulomb yield criterion is used for both soils. Closed‐form solutions are derived, which provide the stress and strain fields during the expansion of the cavity from an initial to a final radius. The analytical solution is validated against finite element simulations, and the effect of varying geometric and material parameters is studied. The influence of the two different soils during cavity expansion is discussed by using pressure–expansion curves and by studying the development of plastic regions within the soils. The analytical method may be applied to various geotechnical problems, which involve aspects of soil layering, such as cone penetration test interpretation, ground‐freezing around shafts, tunnelling, and mining. © 2014 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons Ltd.     

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.     

Drained cavity expansion in sands exhibiting particle crushing

The expansion of cylindrical and spherical cavities in sands is modelled using similarity solutions. The conventional Mohr–Coulomb failure criterion and the state parameter sand behaviour model, which enables hardening–softening, are used in the analysis. The sand state is defined in terms of a new critical state line, designed to account for the three different modes of compressive deformation observed in sands across a wide range of stresses including particle rearrangement, particle crushing and pseudoelastic deformation. Solutions are generated for cavities expanded from zero and finite radii and are compared to those solutions where a conventional critical state line has been used. It is shown that for initial states typical of real quartz sand deposits, pseudoelastic deformation does not occur around an expanding cavity. Particle crushing does occur at these states and causes a reduction in the stress surrounding the cavity. This has major implications when using cavity expansion theory to interpret the cone penetration test and pressuremeter test. Copyright © 2002 John Wiley & Sons, Ltd.     

Undrained expansion of a cylindrical cavity in clays with fabric anisotropy: theoretical solution

This paper presents a novel, exact, semi-analytical solution for the quasi-static undrained expansion of a cylindrical cavity in soft soils with fabric anisotropy. This is the first theoretical solution of the undrained expansion of a cylindrical cavity under plane strain conditions for soft soils with anisotropic behaviour of plastic nature. The solution is rigorously developed in detail, introducing a new stress invariant to deal with the soil fabric. The semi-analytical solution requires numerical evaluation of a system of six first-order ordinary differential equations. The results agree with finite element analyses and show the influence of anisotropic plastic behaviour. The effective stresses at critical state are constant, and they may be analytically related to the undrained shear strength. The initial vertical cross-anisotropy caused by soil deposition changes towards a radial cross-anisotropy after cavity expansion. The analysis of the stress paths shows that proper modelling of anisotropic plastic behaviour involves modelling not only the initial fabric anisotropy but also its evolution with plastic straining.     

A large deformation analysis for the assessment of failure induced deformations of slopes in strain softening materials

This paper presents a numerical method that can be used to evaluate the post-failure deformation of slopes and embankments. The method is based on a large deformation finite element analysis employing the updated Lagrangian formulation. It can simulate the changing geometry of slopes during failure. An extended Mohr–Coulomb constitutive model represents the strain softening behaviour of slope material. The results of a series of analyses show that the post-failure deformation of slopes is a function of the rate of strength reduction and also the stiffness of the slope material. These two factors influence the initiation of progressive failure of slopes. The validity of the method is evaluated by simulation of a test embankment failure and comparison of the predicted results with the observed record of the failure.     

Liquefaction Analysis Using Viscoplastic Kinematic Hardening Constitutive Model

The dynamic response analysis combined with the generalized return-mapping algorithm is applied to the integration algorithms of viscoplastic constitutive relations including the effect of the shear band. The kinematic hardening model based on modified and extended soil model with isotropic strain-hardening–softening is employed. In this paper, the TESRA (temporary effect of strain rate and acceleration) model is employed for the nonlinear viscosity of sand. The constitutive equations of rate-dependent plasticity originally proposed by Duvaut–Lions are employed as the base of the solutions. Liquefaction of a buried pipe is analyzed by finite element method by employing the above mentioned constitutive relations and the calculated results are compared with experimental results. The dynamic response analysis is applied to the solutions of the problems. The kinematic hardening–softening viscoplastic constitutive relations for geomaterials are promising for the predictions of cumulative deformations and liquefaction of the buried pipe. A great deal of experimental results indicate that the stress is a unique function of irreversible strain and its rate.     

Modelling large plastic deformations of cohesive soils using sequential limit analysis

This paper introduces sequential limit analysis (SLA) as a method for modelling large plastic deformations of purely cohesive materials such as undrained clay. The method involves solving a series of consecutive small‐deformation plastic collapse problems using finite element limit analysis, thus ensuring high levels of accuracy, efficiency, and robustness. The techniques needed to develop an SLA implementation for two‐dimensional (plane strain) problems are described in detail, including model geometry updating routines, treatment of rigid bodies, interfaces and boundaries, and periodic remeshing and interpolation of field variables. A simple total stress‐based constitutive model is used to account for strain softening and strain rate effects. Extensive verifications and validations are performed using analytical solutions and physical model test results, comparing both collapse loads and failure mechanisms, to demonstrate the effectiveness of the SLA approach. Additional solution quality checks on the bracketing discrepancy between lower‐bound and upper‐bound limit analysis solutions, and on the incompressibility of the rigid‐plastic material, are also presented.     

A Physical and Numerical Investigation on the Stability of Shallow Tunnels in Strain Softening Media

Summary The behaviour up to failure of shallow underground openings is discussed on the basis of some laboratory, small-scale model tests and of finite element simulation. The experimental results are first illustrated. They were obtained from two-dimensional (plane strain) and three-dimensional tunnel models tested under standard gravity conditions. Then, the phenomenon of strain localisation that characterizes the medium surrounding the model tunnels is discussed, recalling two alternative approaches for its numerical interpretation. On this basis, a finite element procedure for strain softening analyses is outlined and applied to the simulation of the tests in both two- and three-dimensional conditions. The comparison between experimental and numerical results leads to some conclusions on the influence of strain localisation on the overall behaviour of shallow tunnels and on the stability of their headings.     

Probabilistic elasto-plasticity: solution and verification in 1D

In this paper, a solution is presented for evolution of probability density function (PDF) of elastic–plastic stress–strain relationship for material models with uncertain parameters. Developments in this paper are based on already derived general formulation presented in the companion paper. The solution presented is then specialized to a specific Drucker–Prager elastic–plastic material model. Three numerical problems are used to illustrate the developed solution. The stress–strain response (1D) is given as a PDF of stress as a function of strain. The presentation of the stress–strain response through the PDF differs significantly from the traditional presentation of such results, which are represented by a single, unique curve in stress–strain space. In addition to that the numerical solutions are verified against closed form solutions where available (elastic). In cases where the closed form solution does not exist (elastic–plastic), Monte Carlo simulations are used for verification.     

Finite deformation analysis of geomaterials

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.