Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Limitations of suction‐controlled triaxial tests in the characterization of unsaturated soils

Behavior of unsaturated soils is influenced by many factors, and the influences of these factors are usually coupled together. Suction‐controlled triaxial (SCTX) tests are considered to allow researchers to investigate influences of individual variables on unsaturated soils under specified stress path with controls of stresses, pore water, and air pressures. In the past 50 years, SCTX testing method has been established as a standard approach to characterize constitutive behavior of unsaturated soils. Most important concepts for modern unsaturated soil mechanics were developed upon results from the SCTX tests. Among these, one of the most important contributions in the constitutive modeling of elasto‐plastic behavior for unsaturated soils is the Barcelona basic model (BBM) proposed by Alonso et al. in 1990. The BBM successfully explained many features of unsaturated soils and received extensive acceptance. However, the SCTX tests are designed based upon the divide‐and‐conquer approach in which an implicit assumption is used: soil behavior is stress‐path independent. However, it is well‐established that unsaturated soil behavior is elasto‐plastic and stress‐path dependent. It is found that the SCTX tests in fact cannot control the stress path of an unsaturated soil during loading. This incapability, in combination with complicated loading/collapse behavior of unsaturated soils, makes the SCTX tests for characterizing unsaturated soil questionable. This paper discusses the limitations of the SCTX tests in the characterization of unsaturated soils. A possible solution to the problem was proposed based on a newly developed modified state surface approach. The discussions are limited for isotropic conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastoplastic model with loading memory surfaces (LMS) for monotonic and cyclic behaviour of geomaterials

A new elastoplastic model called loading memory surface based on the critical state concept and the multi‐surface framework is proposed for geomaterials. The model uses a hypoelastic formulation and two plastic mechanisms. The formulations of the model are made in three‐dimensional stress–strain space and work under both monotonic and cyclic loadings. A newly introduced formalism makes it possible to obtain the cyclic response directly from the monotonic loading one. This formalism gives a three‐dimensional generalization of the well‐known Masing rule. The model has been validated against test results of Hostun sand under several conditions: monotonic and cyclic, drained and undrained, tests in compression and in extension, and at different confining pressures and different densities. Copyright © 2014 John Wiley & Sons, Ltd.     

A numerical study of mechanical behavior of a cement paste under mechanical loading and chemical leaching

Based on relevant experimental data of a petroleum cement paste under mechanical loading and chemical leaching, an elastic‐plastic model is first proposed by taking into account plastic shearing and pore collapse. The degradation of mechanical properties induced by the chemical leaching is characterized by a chemical damage variable which is defined as the increase of porosity. Both elastic and plastic properties of the cement paste are affected by the chemical damage. The proposed model is calibrated from and applied to describe mechanical responses in triaxial compression tests respectively on sound and fully leached samples. In the second part, a phenomenological chemical model is defined to establish the relationship between porosity change and calcium dissolution process. The dissolution kinetics is governed by a diffusion law taking into account the variation of diffusion coefficient with calcium concentration. The chemical model is coupled with the mechanical model, and both are applied to describe mechanical response of cement paste samples subjected to progressive chemical leaching and compressive stresses. Comparisons between experimental data and numerical results are presented.     

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.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

A bounding surface model for anisotropically overconsolidated clay incorporating thermodynamics admissible rotational hardening rule

A bounding surface model is formulated to simulate the behavior of clays that are subject to an anisotropic consolidation stress history. Conventional rotational hardening is revisited from the perspective of thermodynamics. As the free energy cannot be accumulated infinitely upon critical state failure, the deviatoric back stress must vanish. This requires the rotated yield surface to be turned back to eventually align on the hydrostatic axis in the stress plane. Noting that most of the previous propositions violate this restriction, an innovative rotational hardening rule is formulated that is thermodynamically admissible. The bounding surface framework that employs the modified yield surface is applied to simulate elastoplastic deformations for overconsolidated clays, with which the overprediction of strength on the “dry” side can be greatly improved with reasonable results. Other important features, including contractive or dilative response and hardening or softening behavior, can also be well-captured. It has been shown that the model can simulate three types of reconstituted clays that are sheared with initial conditions over a wide range of anisotropic consolidation stress ratios and overconsolidation ratios under both triaxial undrained and drained conditions. Limitations and potential improvement of the model regarding the fabric anisotropy at critical state have been discussed.     

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.     

Enhancement of backtracking search algorithm for identifying soil parameters

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.