Diffuse failure in geomaterials: Experiments,theory and modelling

This paper presents a synthesis of the works performed by various teams from France, Italy and Canada around the question of second‐order work criterion. Because of the non‐associative character of geomaterials plastic strains, it is now recognized that a whole bifurcation domain exists in the stress space with various possible modes of failure. In a first part these failure modes are observed in lab experimental tests and in discrete element modelling. Then a theoretical study of second‐order work allows to establish a link with the kinetic energy, giving a basis to explain the transition from a prefailure (quasi)static regime to a postfailure dynamic regime. Eventually the main features of geomaterials failure are obtained by applying second‐order work criterion to five different constitutive rate‐independent models—three being phenomenological and two micromechanical. As a whole this paper tries to gather together all the elements for a proper understanding and use of second‐order work criterion in geomechanics. Copyright © 2010 John Wiley & Sons, Ltd.     

Material stability analysis of rock joints

For prediction of rockfalls, the failure of rock joints is studied. Considering these failures as constitutive instabilities, a second‐order work criterion is used because it explains all divergence instabilities (flutter instabilities are excluded). The bifurcation domain and the loading directions of instabilities, which fulfill the criterion, are determined for any piecewise linear constitutive relation. The instability of rock joints appears to be ruled by coupling features of the behavior (e.g., dilatancy). Depending on the loading parameters, instabilities can lead to failure, even before the plastic limit criterion. Results for two given constitutive relations illustrate the approach. Some given loading paths are especially considered. Constant volume (undrained) shear and τ‐constant paths are stable or not depending on the link between the deviatoric stress and strain along undrained paths, as found for soils. Some unstable loading paths are illustrated. Along these paths, failure before the plastic limit criterion is possible. The corresponding failure rules are determined. Copyright © 2012 John Wiley & Sons, Ltd.     

Second‐order work analysis for granular materials using a multiscale approach

It has been established that the second‐order work criterion is a general necessary condition for all instabilities by divergence in rate‐independent granular materials. The relation between the values of discrete second‐order work at the intergranular contact point level and its global macroscopic value is recalled at the beginning of this paper. Then, the basic purpose of the paper is tackled by an analysis of the main features of second‐order work criterion in relation with the granular microstructure. For that, it is considered a novel micromechanical model (the so‐called ‘H‐microdirectional model’), which has the property to involve three scales: grain scale, mesoscale with a specific granular configuration and continuum mechanics macroscale. Eventually, these exhibited features (a bifurcation stress domain and some instability cones) are qualitatively compared with the ones provided by direct numerical simulations issued from a discrete element model. The ultimate goal is to analyse what happens at the granular scale, when the macrosecond‐order work is vanishing at the macrolevel. Copyright © 2013 John Wiley & Sons, Ltd.     

Bifurcation and generalized mixed loading conditions in geomaterials

The concept of bifurcation has been widely discussed, based, for example, on homogeneous tests in soil mechanics. Essentially two‐dimensional (axisymmetric) conditions have been considered. This paper aims at deriving a bifurcation criterion in geomechanics, independently of the problem's dimension. This criterion is related to the vanishing of the determinant of the symmetric part of the material's (or the discrete system's) constitutive (or stiffness) matrix. The derivation is essentially based on the notion of loading parameters (controlling the loading). Basically, a bifurcation occurs when the existence of a unique solution to the quasistatic problem, involving a given set of loading parameters, is lost. Interestingly, the criterion is shown to be independent of the choice of loading parameters. As a possible extension, the context of structure mechanics is considered, and the close connection between both soil and structure analysis is discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Diffuse instabilities with transition to localization in loose granular materials

This paper is concerned with diffuse and other ensuing failure modes in geomaterials when tested under homogeneous states of shearing in various loading programs and drainage conditions. Material instability is indeed the basic property that accounts for the instability of an initially homogeneous deformation field leading to diffuse failure and strain localization in geomaterials. The former is normally characterized by a runaway type of failure accompanied with a sudden and violent collapse of the material in the absence of any localization phenomena. Against this backdrop, we present a brief overview of material instability in elastoplastic solids where one finds a rich source of theoretical concepts including bifurcation, strain localization, diffuse failure and second‐order work, as well as a considerable body of experiments. Some compelling laboratory experimental studies of material instability with focus to diffuse failure are then presented and interpreted based on the second‐order work. Finally, various material instability analyses using an elastoplastic constitutive and a general finite element analysis of the above‐mentioned laboratory experimental tests are presented as a boundary value problem. It is shown that instability can be captured from otherwise uniform stress, density and hydraulic states, whereas uniform deviatoric loads are being applied on the external boundaries of a specimen. Although the numerical simulations reproduce well the laboratory experimental results, they also highlight the hierarchy of failure modes where localization phenomena emerge in the post‐bifurcation regime as a result of a degradation of homogeneity starting from a diffuse mode signalled by a zero second‐order work. Copyright © 2012 John Wiley & Sons, Ltd.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

Predicting the initiation of static liquefaction of cross‐anisotropic sands under multiaxial stress conditions

Experimental evidence has shown that the liquefaction instability of sands can be affected by its material density, stress state, and inherent anisotropy. In order to predict the initiation of the static liquefaction of inherent cross‐anisotropic sands under multidimensional stress conditions, a rational constitutive model is needed. An elastoplasticity model able to capture the influences of intermediate principal stress ratio (b  = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) and loading direction on stress–strain relationships and volumetric properties was proposed. The yield function was formulated to be controlled by Lode angle, loading direction, and material state; the stress–dilatancy was a material state‐dependent function. After using the existing drained hollow cylinder tests to validate the proposed model, this model was used to simulate the existing undrained hollow cylinder tests. During this simulation, the second‐order work criterion was used to determine the initiation of static liquefaction. The results showed that an increase in both the intermediate principal stress ratio and the loading angle induces a decrease in the second‐order work. Static liquefaction is initiated more easily at a stress state with a large intermediate principal stress ratio and a large loading angle, and the mobilized friction angle at the instability points decreases with the intermediate principal stress ratio and the loading angle. Copyright © 2017 John Wiley & Sons, Ltd.     

Semi‐analytical far field model for three‐dimensional finite‐element analysis

A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd.     

Spectral analysis of nonlocal regularization in two‐dimensional finite element models

Strain‐softening in geomaterials often leads to ill‐posed boundary‐valued problems (BVP), which cannot be solved with finite element methods without introducing some kind of regularization such as nonlocal plasticity. Hereafter we propose to apply spectral analysis for testing the performance of nonlocal plasticity in regularizing ill‐posed BVP and producing mesh‐independent solutions when local plasticity usually fails. The spectral analysis consists of examining the eigenvalues and eigenvectors of the global tangential stiffness matrix of the incremental equilibrium equations. Based on spectral analysis, we propose a criterion for passing or failing the test of constitutive regularization in the context of BVP. If the eigenvalues of the tangential operator are all positive then the regularization succeeds, otherwise it fails and may not prevent artificial mesh‐dependent solutions from appearing. The approach is illustrated in the particular case of a biaxial compression with strain‐softening plasticity. In this particular case, local softening plasticity is found to produce negative eigenvalues in the tangential stiffness matrix, which indicates ill‐posed BVP. In contrast, nonlocal softening plasticity always produces positive eigenvalues, which regularizes ill‐posed BVP. The dominant eigenvectors, which generate localized deformation patterns, have a bandwidth independent of mesh size, provided that the mesh is fine enough to capture localization. These mesh‐independent eigenmodes explain why nonlocal plasticity produces numerical solutions that are mesh‐independent. Copyright © 2011 John Wiley & Sons, Ltd.     

Solid–fluid transition modelling in geomaterials and application to a mudflow interacting with an obstacle

Given the contrasting behaviour observed for geomaterials, for example, during landslides of the flow type, this contribution proposes an original constitutive model, which associates both an elasto‐plastic relation and a Bingham viscous law linked by a mechanical transition criterion. This last is defined as the second‐order work sign for each material point, which is a general criterion for divergence instabilities. Finite element method with Lagrangian integration points is chosen as a framework for implementing the new model because of its well‐known ability to deal with both solid and fluid behaviours in large deformation processes. A first boundary model considering a sample of initially stable soil, a slope and an obstacle is performed. The results show the power of the constitutive model because the consistent evolution of initiation, propagation and arrest of the mudflow is described. A parametric study is led on various plastic and viscous parameters to determine their influence on the flow development and arrest. Finally, forces against the obstacle are compared with good agreement with those of other authors for the same geometry and a pure viscous behaviour. Copyright © 2014 John Wiley & Sons, Ltd.     

Diffuse and localized failure modes: Two competing mechanisms

The concept of failure is one of the most debated in soil mechanics, for two reasons essentially. First, this is a crucial issue in the engineering of structures and geotechnical project design. Second, this is still a challenging academic issue mobilizing significant scientific interest in the development of a unique framework to describe the different failure modes. In this respect, this paper revisits the localized failure mode, replacing the well‐known Rice criterion within the wider context of bifurcation. Considering a micro‐mechanical model, the main theoretical results are covered. In particular, it is established that localized failure is a particular case of failures observed within the so‐called bifurcation domain: the incremental strain within the localization band is associated with a vanishing value of the second‐order work. Copyright © 2010 John Wiley & Sons, Ltd.     

A stability criterion for elasto‐viscoplastic constitutive relationships

In this paper, the onset of mechanical instability in time‐sensitive elasto‐viscoplastic solids is theoretically analyzed at the constitutive level and associated with the occurrence of ‘spontaneous accelerations’ under stationary external perturbations. For this purpose, a second‐order form of Perzyna's constitutive equations is first derived by time differentiation, and a sufficient stability condition is identified for general mixed loading programs. These loading conditions are in fact the most general in both laboratory tests and real boundary value problems, where a combination of certain stress and strain components is known/prescribed. The theoretical analysis leads to find precise stability limits in terms of material hardening modulus. In the case of constitutive relationships with isotropic strain‐hardening, no instabilities are possible while the hardening modulus is larger than the so‐called ‘controllability modulus’ defined for (inviscid) elasto‐plastic materials. It is also shown that the current stress/strain rate may also directly influence the occurrence of elasto‐viscoplastic instability, which is at variance with elasto‐plastic inviscid media. Copyright © 2015 John Wiley & Sons, Ltd.     

Analytical study of mine closure behaviour in a poro‐visco‐elastic medium

This paper is interested in the hydro‐mechanical behaviour of an underground cavity abandoned at the end of its service life. It is an extension of a previous study that accounted for a poro‐elastic behaviour of the rock mass (Int. J. Comput. Geomech. 2007; DOI: 10.1016/j.compgeo.2007.11.003 ). Deterioration of the lining support with time leads to the transfer of the loading from the exterior massif to the interior backfill. The in situ material has a poro‐visco‐elastic constitutive behaviour while the backfill is poro‐elastic, both saturated with water. This loading transfer is accompanied by an inward cavity convergence, thereby compressing the backfill, and induces an outward water flow. This leads to a complex space–time evolution of pore pressures, displacements and stresses, which is not always intuitive. In its general setting, a semi‐explicit solution to this problem is developed, using Laplace transform, the inversion being performed numerically. Analytical inversion leading to a quasi‐explicit solution in the time domain is possible by identifying the characteristic creep and relaxation times of volumetric strains with those of the deviatoric strains, on the basis of a parametric study. A few numerical examples are given to illustrate the hydro‐mechanical behaviour of the cavity and highlight the influence of key parameters (e.g. stiffness of backfill, lining deterioration rate, etc.). Further studies accounting for more general material behaviours for the backfill and external ground are ongoing. Copyright © 2008 John Wiley & Sons, Ltd.     

Finite element modelling of thick plates on two‐parameter elastic foundation

This paper is intended to give some information about how to build a model necessary for bending analysis of rectangular and circular plates resting on a two‐parameter elastic foundation, subjected to combined loading and permitting various types of boundary conditions. The formulation of the problem takes into account the shear deformation of the plate and the surrounding interaction effect outside the plate. The numerical model based on an 18‐node zero‐thickness isoparametric interface element interacting with a thick Reissner–Mindlin plate element with three degrees of freedom at each of the nine nodes, which enforce C⁰ continuity requirements for the displacements and rotations of the midsurface, is proposed. Stiffness matrices of a special interface element are superimposed on the global stiffness matrix to represent the stiffening elastic foundation under and beyond the plate. Some numerical examples are given to illustrate the advantages of the method presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Bifurcation and second‐order work in geomaterials

In this paper, the ability of a material rate‐independent system to evolve toward another mechanical state from an equilibrium configuration, with no change in the control parameters, is investigated. From a mechanical point of view, this means that the system can spontaneously develop kinetic energy with no external disturbance from an equilibrium state, which corresponds to a particular case of bifurcation. The existence of both conjugate incremental strain and stress such that the second‐order work vanishes is established as a necessary and sufficient condition for the appearance of this bifurcation phenomenon. It is proved that this fundamental result is independent of the constitutive relation of the rate‐independent material considered. Then the case of homogeneous loading paths is investigated, and, as an illustration, the subsequent results are applied to interpret the well‐known liquefaction observed under isochoric triaxial loading conditions with loose granular materials. Copyright © 2006 John Wiley & Sons, Ltd.     

Visualization of material stiffness in geomechanics analysis

This paper presents novel visualization techniques to simplify representation of the fourth‐order material stiffness tensor as a set of three‐dimensional geometric objects. Stiffness visualization aids in understanding the complex stiffness characteristics of highly non‐linear constitutive models including modelled material anisotropy and loading path dependent stiffness variation. Stiffness visualization is relevant for understanding the relationship of material stiffness to global behaviour in the analysis of a boundary value problem. The spherical pulse stiffness visualization method, developed in the acoustics field, is extended to visualize stiffness of geomaterials using three three‐dimensional objects. This method is limited to relatively simple constitutive models with symmetric stiffness matrices insensitive to loading magnitude and direction. A strain dependent stiffness visualization method is developed that allows the examination of material stiffness for a range of loading directions and is suitable for highly non‐linear and path dependent material models. The proposed stiffness visualization can be represented as 3‐D, 2‐D and 1‐D objects. The visualization technique is used to represent material stiffness and its evolution during simulated soil laboratory tests and deep excavation construction. Copyright © 2005 John Wiley & Sons, Ltd.     

Diffuse bifurcations engraving diverse shear bands in granular materials

Shear bands with characteristic spatial patterns observed in an experiment for a cubic or parallelepiped specimen of dry dense sand were simulated by numerical bifurcation analysis using the Cam‐clay plasticity model. By incorporating the subloading surface concept into the plasticity model, the model became capable of reproducing hardening/softening and contractive/dilative behavior observed in the experiment. The model was reformulated to be compatible with the multiplicative hyperelasto‐plasticity for finite strains. This enhanced constitutive model was implemented into a finite‐element code reinforced by a stress updating algorithm based on the return‐mapping scheme, and by an efficient numerical procedure to compute critical eigenvectors of elastoplastic tangent stiffness matrix at bifurcation points. The emergence of diamond‐ and column‐like diffuse bifurcation modes breaking uniformity of the materials, followed by the evolution of shear bands through strain localization, was observed in the analysis. In the bifurcation analysis of plane strain compression test, unexpected bifurcation modes, which broke out‐of‐plane uniformity and led to 3‐dimensional diamond‐like patterns, were detected. Diffuse bifurcations, which were difficult to observe by experiments, have thus been found as a catalyst creating diverse shear band patterns.     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.     

Discussing the definition of the second‐order work for unsaturated soils

The volumetric compaction due to wetting processes is a phenomenon observed quite often in unsaturated soils. Under certain circumstances, saturation events can result into a sudden and unexpected collapse of the system. These phenomena are usually referred to as wetting‐induced collapses, without providing any detailed theoretical justification for this terminology. In order to predict in a general fashion the occurrence of coupled instabilities induced by saturation processes, a generalization of the theoretical approaches usually employed for saturated geomaterials is here provided. More specifically, this paper addresses the problem of hydro‐mechanical instability in unsaturated soils from an energy standpoint. For this purpose, an extension of the definition of the second‐order work is here suggested for the case of unsaturated porous media. On the basis of some examples of numerical simulations of laboratory tests, coupled hydro‐mechanical instabilities are then interpreted in the light of this second‐order energy measure. Finally, the implications of the theoretical results here presented are commented from a constitutive modelling perspective. Two possible alternative approaches to formulate incremental coupled constitutive relations are indeed discussed, showing how the onset of hydro‐mechanical instabilities can be predicted using an extended form of Hill's stability criterion. Copyright © 2010 John Wiley & Sons, Ltd.