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.     

A visco‐plastic constitutive model for granular soils modified according to non‐local and gradient approaches

An already available non‐associated elastic–viscoplastic constitutive model with anisotropic strain hardening is modified in order to describe both the constitutive parameter dependency on relative density and the spatio‐temporal evolution of strain localization. To achieve this latter goal, two distinct but similar approaches are introduced: one inspired by the gradient theory and one by the non‐local theory. A one‐dimensional case concerning a simple shear test for a non‐homogeneous infinitely long dense sand specimen is numerically discussed and a finite difference scheme is employed for this purpose. The results obtained by following the two different approaches are critically analysed and compared. Copyright © 2001 John Wiley & Sons, Ltd.     

Continuum large cracking in a rate‐dependent plastic–damage model for cyclic‐loaded concrete structures

Formulation and algorithmic treatment of a rate‐dependent plastic–damage model modified to capture large tensile cracking in cyclic‐loaded concrete structures are presented in detail for a three‐dimensional implementation. The plastic–damage model proposed by Lee and Fenves in 1998 was founded based on isotropic damaged elasticity in combination with isotropic multi‐hardening plasticity to simulate cracking and crushing of concrete under cyclic or dynamic loadings. In order that the model can capture large crack opening displacements, which are inevitable in plain concrete structures, the excessive increase in plastic strain causing unrealistic results in cyclic behaviors is prevented when the tensile plastic–damage variable controlling the evolution of tensile damage is larger than a critical value. In such a condition, the crack opening/closing mechanism becomes similar to discrete cracking. The consistent tangent operator required to accelerate convergence rate is also formulated for the large cracking state including viscoplasticity. The validation and performance of the modified algorithm implemented in a special finite element program is exemplified through several single‐element tests as well as three structural applications. The last example examines the model in the seismic fracture analysis of Koyna dam as a benchmark problem and the resulting crack profile is compared with the available experiment. The numerical experimentations well demonstrate that the developed model whose modification is necessary to properly simulate the cyclic behavior of plain concrete subjected to large tensile strains is robust and reasonably accurate. Copyright © 2012 John Wiley & Sons, Ltd.     

On failure criteria for anisotropic cohesive‐frictional materials

Anisotropic failure criteria are formulated using two different approaches. The first one employs a spatial distribution of strength parameters and defines the failure condition in terms of traction components acting on the critical plane. The second one incorporates a microstructure tensor and the relevant mixed invariants. Both formulations are illustrated by some numerical examples. In particular, the variation of strength with orientation of the sample is examined for a series of uniaxial compression tests. Copyright © 2001 John Wiley & Sons, Ltd.     

Stress‐driven integration strategies and m‐AGC tangent operator for Perzyna viscoplasticity and viscoplastic relaxation: application to geomechanical interfaces

The paper proposes a stress‐driven integration strategy for Perzyna‐type viscoplastic constitutive models, which leads also to a convenient algorithm for viscoplastic relaxation schemes. A generalized trapezoidal rule for the strain increment, combined with a linearized form of the yield function and flow rules, leads to a plasticity‐like compliance operator that can be explicitly inverted to give an algorithmic tangent stiffness tensor also denoted as the m‐AGC tangent operator. This operator is combined with the stress‐prescribed integration scheme, to obtain a natural error indicator that can be used as a convergence criterion of the intra‐step iterations (in physical viscoplasticity), or to a variable time‐step size in viscoplastic relaxation schemes based on a single linear calculation per time step. The proposed schemes have been implemented for an existing zero‐thickness interface constitutive model. Some numerical application examples are presented to illustrate the advantages of the new schemes proposed. Copyright © 2016 John Wiley & Sons, Ltd.     

A rate‐dependent cohesive crack model based on anisotropic damage coupled to plasticity

In quasi‐brittle material the complex process of decohesion between particles in microcracks and localization of the displacement field into macrocracks is limited to a narrow fracture zone, and it is often modelled with cohesive crack models. Since the anisotropic nature of the decohesion process in separation and sliding is essential, it is particularly focused in this paper. Moreover, for cyclic and dynamic loading the unloading, load reversal (including crack closure) and rate dependency are essential features that are included in a new model. The modelling of degradation is based on a ‘localized’ version of anisotropic continuum damage coupled to inelasticity. The concept of strain energy equivalence between the states in the effective and nominal settings is adopted in order to define the free energy of the interface. The proposed fracture criterion is of the Mohr type, with a smooth transition of the failure and kinematics (slip and dilatation) characteristics between tension and shear. The chosen potential, of the Lemaitre‐type, for evolution of the dissipative processes is additively decomposed into plastic and damage parts, and non‐associative constitutive equations are obtained. The constitutive equations are integrated by applying the backward Euler rule and by using Newton iteration. The proposed model is assessed analytically and numerically and a typical calibration procedure for concrete is proposed. Copyright © 2006 John Wiley & Sons, Ltd.     

Stress–strain model for clays with anisotropic void ratio distribution

The porosity of soils is considered to be a directional measure and its distribution is characterized by a functional form. This form has been used to extend the critical state soil mechanics framework to include the effects of structure in soils. A new internal plastic energy dissipation formulation has been proposed to account for fabric arrangement. New expressions for the yield locus, and the plastic stress–strain response of structural soils have been derived. The applicability of the concepts to model the plastic stress–strain behaviour of a number of soils is illustrated. The advantage of the new model is very well identified in modelling the stress–strain behaviour of K₀ consolidated and natural clays. © 1998 John Wiley & Sons, Ltd.     

Localized failure of fibre‐reinforced elastic–plastic materials subjected to plane strain loading

We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd.     

A unified approach to the implicit integration of standard,non‐standard and viscous plasticity models

It is shown how modern concepts to integrate the elasto‐plastic rate equations of standard plasticity via an implicit algorithm can be generalized to plasticity without an explicitly defined yield surface and to overstress‐type models of viscoplasticity, where the stress point can be located outside the loading surface. For completeness, a tangent operator is derived that is consistent with the update algorithm. Copyright © 2002 John Wiley & Sons, Ltd.     

A two‐step Taylor–Galerkin algorithm applied to shock wave propagation in soils

This paper presents a two‐step Taylor–Galerkin algorithm formulated in terms of velocities and stresses which can be applied to solid dynamics problems requiring good resolution of small wavelengths, such as propagation of shocks. The proposed model is both fast, as it uses simple linear elements (triangles in 2D and tetrahedra in 3D), and accurate. It avoids locking and mesh alignment problems, and therefore can be applied to localized failure computations without the limitations exhibited by the classical displacement formulations. Copyright © 2003 John Wiley & Sons, Ltd.     

An anisotropic elastic‐decohesive constitutive relation for sea ice

Thermodynamic growth or melt and mechanical redistribution due to lead opening or ridge formation shape the thickness distribution of the Arctic ice cover and impact the overall strength of pack ice. Specifically, the deformation and strength of ice are not isotropic but vary with the thickness and lead orientation. To reflect these facts, we develop an anisotropic, elastic‐decohesive constitutive model for sea ice together with a model to describe an oriented, ice thickness distribution. The tight connection between the mechanical response and the thickness distribution is an improvement over a previous model that only depended on the average ice thickness. The model describes mechanical responses anisotropically in both the elastic and failure regimes. In the elastic regime, the constitutive relation implicitly reflects strong and weak directions of the pack ice depending on the distribution of thin ice (including open water) and thicker ice (e.g., multi‐year ice or ridges). In the failure regime, the model predicts both failure initiation and the lead orientation. Evolution from initial failure to complete failure when traction‐free crack surfaces are formed is also modeled. Crack or lead width is determined during the evolution. Various examples of failure surfaces are presented to describe the behavior of modeled ice when the thickness distribution varies. The model predictions are also illustrated and compared with previous modeling efforts by examining regions of ice under idealized loading. Copyright © 2015 John Wiley & Sons, Ltd.     

Consolidation analysis of saturated multi‐layered soils with anisotropic permeability caused by a point sink

This paper presents the analytical layer‐element method to analyze the consolidation of saturated multi‐layered soils caused by a point sink by considering the anisotropy of permeability. Starting from the governing equations of the problem, the solutions of displacements and stresses for a single soil layer are obtained in the Laplace–Hankel transformed domain. Then, the analytical layer‐element method is utilized to further derive the solutions for the saturated multi‐layered soils in the transformed domain by combining with the boundary conditions of the soil system and continuity conditions between adjacent layers. The actual solutions in the physical domain can be acquired by the inversion of Laplace–Hankel transform. Numerical results are carried out to show the accuracy and stability of the proposed method and evaluate the influence of sink depth and anisotropic permeability on excess pore pressure and surface settlement. Copyright © 2011 John Wiley & Sons, Ltd.     

A two‐surface thermomechanical model for saturated clays

This paper presents an advanced thermomechanical model – TEAM in the framework of two‐surface plasticity for saturated clays, with emphasis put on some important thermomechanical features of natural clays evidenced experimentally such as the limited thermomechanical elastic zone, the smooth transition from elastic to plastic behavior. Two plastic mechanisms are introduced in the model: one is to reproduce the thermoplasticity involving thermal expansion and contraction observed at high over‐consolidation ratios and the second one describes the temperature effect on the yield behavior. The model adopts additional yield surfaces, namely inner yield surfaces that are associated with the two proposed plastic mechanisms to account for the plastic behavior inside the existing conventional thermomechanical yield surface namely yield surfaces. The general expressions of the yield surfaces and plastic potentials in p′–q–T space are introduced. A progressive plastic hardening mechanism associated with the inner yield surface is defined, enabling the plastic modulus to vary smoothly during thermomechanical loadings inside the yield surfaces. Several tests on natural Boom clay along different thermomechanical loading paths have been simulated by TEAM, and results show its relevance in describing the thermomechanical behavior of saturated clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Homogenization‐based analysis of anisotropic damage in brittle materials with unilateral effect and interactions between microcracks

This paper is devoted to micromechanical modeling of induced anisotropic damage in brittle geomaterials. The formulation of the model is based on a proper homogenization procedure by taking into account unilateral effects and interactions between microcracks. The homogenization procedure is developed in the framework of Eshelby's inclusion solution and Ponte‐Castaneda and Willis (J. Mech. Phys. Solids 1995; 43 :1919–1951) estimate. The homogenization technique is combined with the thermodynamics framework at microscopic level for the determination of damage evolution law. A rigorous crack opening–closure transition condition is established and an energy‐release‐rate‐based damage criterion is proposed. Computational aspects on the implementation of micromechanical model are also discussed. The proposed model is evaluated by comparing numerical predictions with experimental data for various laboratory tests on concrete. Parametric studies on unilateral effects and influences of microcracks interactions are finally performed and analyzed. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.     

A Runge–Kutta,Taylor–Galerkin scheme for hyperbolic systems with source terms. Application to shock wave propagation in viscoplastic geomaterials

This paper presents an alternative formulation of Solid Dynamics problems based on (i) a mathematical model consisting of a system of hyperbolic PDEs where the source term is originated by the viscoplastic strain rate and (ii) a splitting scheme where the two‐step Taylor–Galerkin is used for the advective part of the PDE operator while the sources are integrated using a fourth‐order Runge–Kutta. Use of the splitting scheme results in a higher accuracy than that of the original two‐step Taylor–Galerkin. The scheme performs well when used with linear triangle or tetrahedra for (i) bending‐dominated situations (ii) localized failure under dynamic conditions and keeps the advantages of the two‐step Taylor–Galerkin concerning numerical dispersion and damping of short wavelengths. Copyright © 2006 John Wiley & Sons, Ltd.     

An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

A constitutive model for the dynamic and high‐pressure behaviour of a propellant‐like material: Part II: Model development and applications

This paper presents the second part of a work that aims at developing a mechanical model for the behaviour of propellant‐like materials under high confining pressure and strain rate. The model is to be employed to determine the temperature rise due to mechanical dissipation during a dynamic penetration event. Using the kinematical and thermodynamic background derived in the first part, a viscoelastic–viscoplastic–compaction model is put forward and identified. Viscoelasticity and compaction refer to the elastomeric nature of the material, while viscoplasticity, including implicit damage through dilatancy, reflects its granular nature. Some numerical exercises are performed, in view of determining the major model forces and weaknesses, and of assessing the numerical algorithm robustness. A penetration event is preliminarily simulated, and a temperature rise field predicted. Necessary model improvements are finally discussed, together with several ways of research for a longer term approach. Copyright © 2001 John Wiley & Sons, Ltd.     

Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading

This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd.     

Three‐dimensional finite elements for the analysis of soil contamination using a multiple‐porosity approach

A three‐dimensional finite‐element model of contaminant migration in fissured clays or contaminated sand which includes multiple sources of non‐equilibrium processes is proposed. The conceptual framework can accommodate a regular network of fissures in 1D, 2D or 3D and immobile solutions in the macro‐pores of aggregated topsoils, as well as non‐equilibrium sorption. A Galerkin weighted‐residual statement for the three‐dimensional form of the equations in the Laplace domain is formulated. Equations are discretized using linear and quadratic prism elements. The system of algebraic equations is solved in the Laplace domain and solution is inverted to the time domain numerically. The model is validated and its scope is illustrated through the analysis of three problems: a waste repository deeply buried in fissured clay, a storage tank leaking into sand and a sanitary landfill leaching into fissured clay over a sand aquifer. Copyright © 2005 John Wiley & Sons, Ltd.