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.     

Non‐coaxiality and stress–dilatancy rule in granular materials: FE investigation within micro‐polar hypoplasticity

This paper deals with FE investigations of shear localization in dilatant granular bodies. The calculations were carried out with a hypoplastic constitutive law enhanced by micro‐polar terms to properly model the shear zone evolution. The behaviour of an initially medium dense sand specimen with very smooth and very rough horizontal boundaries was analyzed during a plane strain compression test. A stochastic distribution of the initial void ratio was assumed to be spatially correlated. Attention was focused on the non‐coaxiality of the directions of the principal strain increments and principal stresses in the shear zone and on the stress–dilatancy rule. Copyright © 2008 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.     

A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation

Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities.     

A single numerically efficient equation for approximating the Mohr–Coulomb and the Matsuoka–Nakai failure criteria with rounded edges and apex

This paper presents a novel formulation for defining soil failure. It plots in the principal stress space as a surface with the shape ranging between an approximation of the Matsuoka–Nakai and of the Mohr–Coulomb criteria depending on the value of a single parameter. The new function can be used as a replacement of the original equations of these well‐established criteria for implementing in a program for numerical analyses, and it is particularly effective for approximating the Matsuoka–Nakai criterion. Both the Mohr–Coulomb and the Matsuoka–Nakai failure criteria present numerical difficulties during implementation and also at run‐time. In the case of the Matsuoka–Nakai, the new formulation plots in the first octant only, whereas the original criterion plots in all octants, which causes severe convergence problems particularly for those Gauss points with low stress state, such as those on the side of a shallow footing. When the shape parameter is set to reproduce the Mohr–Coulomb failure criterion, on the other hand, the new formulation plots as a pyramid with rounded edges. Moreover, as the new function is at least of class C², the second derivatives are continuous, thus ensuring quadratic convergence of the Newton's method used within the integration scheme of the constitutive law. The proposed formulation can also provide both sharp and rounded apex of the surface at the origin of the stress space by setting accordingly one additional parameter. Copyright © 2013 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.