Elasto‐plastic model for cement‐treated sand

This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd.     

Elasto‐plastic modelling of stress‐path‐dependent behaviour of interfaces

A single‐surface elasto‐plastic model developed by Desai and his coworkers is used to predict the behaviour of an interface between sand and a steel plate. The loading in the experiments and in their predictions followed various stress and displacement paths. The results of predictions of the two‐ and three‐dimensional behaviour of the interface under both constant normal stress and constant normal stiffness conditions are presented. The predictions are compared with their corresponding experimental results. The model parameters were determined on the basis of 2‐D conventional experiments under the condition of constant normal stress and they were used in the prediction of the interface behaviour in various stress paths. There is, in general, a good agreement between the predicted and experimental results. Copyright © 2000 John Wiley & Sons, Ltd.     

Template elastic‐plastic computations in geomechanics

In this paper we present a new approach to computations in elasto‐plastic geomechanics. The approach is based on the object oriented design philosophy and observations on similarity of most incremental elastic–plastic material models. This new approach to elastic–plastic computations in geomechanics allows for creation of template material models. The analysis of template material models will in turn allow for an easy implementation of other elastic–plastic material models based on the object oriented design principles. Furthermore we present some illustrative implementation details. Finally we present analysis results that emphasize features of template elastic–plastic computations in geomechanics. Copyright © 2002 John Wiley & Sons, Ltd.     

Non‐normality and induced plastic anisotropy under fractional plastic flow rule: a numerical study

In this paper, an implementation of fractional plastic flow rule in the framework of implicit and explicit procedures is under consideration. The fractional plastic flow rule is obtained from a generalisation of the classical plastic flow rule utilising fractional calculus. The key feature of this new concept is that in general, the non‐associative flow is obtained without necessity of additional potential assumption. If needed, the model can cover the anisotropy induced by plastic deformation. Illustrative examples showing the unusual flexibility of this model are also presented. Copyright © 2015 John Wiley & Sons, Ltd.     

Bifurcation and instability modelling by a multimechanism elasto‐plastic model

The bifurcation and instability conditions in geomechanics are closely related to the elasto‐plastic behaviour. In this paper the potential of a multimechanism elasto‐plastic model to predict various modes of failure is examined. First, a brief overview for the essential aspects of the constitutive model and the development of the elasto‐plastic constitutive matrix for this model are presented. Then, numerical simulations of different drained and undrained paths in the axisymmetric and plane‐strain conditions for the Hostun sand are illustrated. These examples confirm the capacity of the model to reproduce instability and strain localization phenomena. The obtained response is in agreement with experimental observations, theoretical developments and numerical analyses existing in the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Poro‐elasto‐plastic response of an unconsolidated formation confined with stiff seal rocks under radial injection

In this article, we evaluate geomechanics of fluid injection from a fully penetrating vertical well into an unconsolidated formation confined with stiff seal rocks. The coupled behavior of an isotropic, homogeneous sand layer is studied under injection pressures that are high enough to induce plasticity yet not fracturing. Propagation of the significant influence zone surrounding the injection borehole, quantified by the extent of the plastic domain in the elasto‐plastic model, is examined for the first time. First, a new fully coupled axisymmetric numerical model is developed. A comprehensive assessment is performed on pore pressures, stresses/strains, and failure planes during the entire transient period of an injection cycle. Results anticipate existence of five distinctive zones in terms of plasticity state: liquefaction at wellbore; two inner plastic domains surrounding the wellbore, where failure occurs along two planes and major principal stress is in vertical direction; remaining of the plastic domain, where formation fails along one plane and major principal stress is in radial direction; and a non‐plastic region. Failure mechanism at the wellbore is found to be shear followed by liquefaction. Next, a novel methodology is proposed based on which new weakly coupled poro‐elasto‐plastic analytical solutions are derived for all three stress/strain components. Unlike previous studies, extension of the plastic zone is obtained as a function of injection pressure, incorporating plasticity effects on the subsequent elastic domain. Solutions, proven to be a good approximation of numerical simulations, offer a huge advantage as the run time of coupled numerical simulations is considerably long. Copyright © 2016 John Wiley & Sons, Ltd.     

Large stress reversals in true triaxial tests on cross‐anisotropic sand

The results of a series of true triaxial tests with stress paths involving large reversals under 3D conditions are presented. These tests were performed on medium dense Santa Monica Beach sand to provide experimental evidence for the rotational kinematic hardening hypothesis presented in a companion paper and to provide stress–strain and volume change relations for experiments with 3D stress paths and large stress reversals to be predicted by the rotational kinematic hardening model. The experimental equipment and the testing procedures are briefly explained followed by a presentation of the experimental results and their sensitivity to unknown causes as well as effects of cross‐anisotropy on the sand behavior. The stress paths are presented in a σ₃′‐plane and in an octahedral plane and the directions of experimental strain increment vectors are compared with those obtained from the rotational kinematic hardening model. Copyright © 2008 John Wiley & Sons, Ltd.     

A thermo‐hydro‐mechanical model for unsaturated soils based on the effective stress concept

A simple thermo‐hydro‐mechanical (THM) constitutive model for unsaturated soils is described. The effective stress concept is extended to unsaturated soils with the introduction of a capillary stress. This capillary stress is based on a microstructural model and calculated from attraction forces due to water menisci. The effect of desaturation and the thermal softening phenomenon are modelled with a minimal number of material parameters and based on existing models. THM process is qualitatively and quantitatively modelled by using experimental data and previous work to show the application of the model, including a drying path under mechanical stress with transition between saturated and unsaturated states, a heating path under constant suction and a deviatoric path with imposed suction and temperature. The results show that the present model can simulate the THM behaviour in unsaturated soils in a satisfactory way. Copyright © 2010 John Wiley & Sons, Ltd.     

Three dimensional elasto‐plastic interface for embedded beam elements with interaction surface for the analysis of lateral loading of piles

This paper presents a numerical formulation for a three dimensional elasto‐plastic interface, which can be coupled with an embedded beam element in order to model its non‐linear interaction with the surrounding solid medium. The formulation is herein implemented for lateral loading of piles but is able to represent soil‐pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The interface is formulated in order to capture localized material plasticity in the soil surrounding the pile within the range of small to moderate lateral displacements. The interface is formulated following two different approaches: (i) in terms of beam degrees of freedoms; and (ii) considering the displacement field of the solid domain. Each of these alternatives has its own advantages and shortcomings, which are discussed in this paper. The paper presents a comparison of the results obtained by means of the present formulation and by other well‐established analysis methods and test results published in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media

This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd.     

盾构偏航引起的地表位移预测

对于盾构隧道的设计和施工，解析方法预测地表沉降，具有使用简便、物理意义明确的特点。在Sagaseta提出的基于由地层损失引起的地表沉降理论基础上，结合工程实际情况，提出了模拟盾构推进过程的三维地层损失模式，并推导了相应的地表位移计算解析公式，对实际工程具有参考价值。     

Unsaturated slope stability analysis with steady infiltration or evaporation using elasto‐plastic finite elements

The paper presents results of unsaturated slope stability analyses using elasto‐plastic finite elements in conjunction with a novel analytical formulation for the suction stress above the water table. The suction stress formula requires four parameters, three for the soil type and one for the steady infiltration (or evaporation) due to environmental effects. The suction stress approach enables the analysis to proceed in the context of classical effective stress, while maintaining the advantages of a general non‐linear finite element approach in which no advance assumptions need to be made about the shape or location of the critical failure surface. The results show the extent to which suctions above the water table can increase the factor of safety of a slope for a variety of different soil types and infiltration rates. All stability analyses that include the effects of suction stresses are contrasted with more traditional approaches in which water pressures above the water table are ignored. Copyright © 2005 John Wiley & Sons, Ltd.     

A non‐coaxial and unbalanced plastic flow rule for geomaterials and its application to the shear band orientation prediction

In this paper, the non‐coaxial relation between the principal plastic strain increments and the principal stresses, which results from the internal friction in geomaterials, is analyzed, and the phenomenon of the unbalanced development of plastic flow in two conjugate directions is discussed. A non‐coaxial, unbalanced plastic flow model for Coulomb frictional materials is developed and used to determine the orientation of shear band in geomaterials. It is shown that the unbalanced index r of plastic flow has important effect on the orientation of the shear band, and the orientation determined by the conventional plastic flow theory is only a special case of the proposed model when r=0. This result soundly explains the reason that the geomaterials with the same internal friction angle and dilatancy angle can have very different shear band orientations. In addition, the difference between the intrinsic and apparent dilatancy angles is analyzed, and it is emphasized that the dilatancy angle commonly used in practice is indeed the apparent dilatancy angle. Copyright © 2010 John Wiley & Sons, Ltd.     

Elasto‐plastic analysis of block structures through a homogenization method

The paper describes the development and numerical implementation of a constitutive relationship for modeling the elasto‐plastic behavior of block structures with periodic texture, regarded at a macroscopic scale as homogenized anisotropic media. The macroscopic model is shown to retain memory of the mechanical characteristics of the joints and of the shape of the blocks. The overall mechanical properties display anisotropy and singularities in the yield surface, arising from the discrete nature of the block structure and the geometrical arrangement of the units. The model is formulated in the framework of multi‐surface plasticity. It is implemented in an finite element (FE) code by means of two different algorithms: an implicit return mapping scheme and a minimization algorithm directly derived from the Haar–Karman principle. The model is validated against analytical and experimental results: the comparison between the homogenized continuum and the original block assembly shows a good agreement in terms of ultimate inelastic behavior, when the size of the block is small as compared with that of the whole assembly. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of earthquake‐induced slope deformation using SPH method

This paper presents the development, calibration, and validation of a smoothed particle hydrodynamics (SPH) model for the simulation of seismically induced slope deformation under undrained condition. A constitutive model that combines the isotropic strain softening viscoplasticity and the modified Kondner and Zelasko rule is developed and implemented into SPH formulations. The developed SPH model accounts for the effects of wave propagation in the sliding mass, cyclic nonlinear behavior of soil, and progressive reduction in shear strength during sliding, which are not explicitly considered in various Newmark‐type analyses widely used in the current research and practice in geotechnical earthquake engineering. Soil parameters needed for the developed model can be calibrated using typical laboratory shear strength tests, and experimental or empirical shear modulus reduction curve and damping curve. The strain‐rate effects on soil strength are considered. The developed SPH model is validated against a readily available and well‐documented model slope test on a shaking table. The model simulated slope failure mode, acceleration response spectra, and slope deformations are in excellent agreement with the experimental data. It is thus suggested that the developed SPH model may be utilized to reliably simulate earthquake‐induced slope deformations. This paper also indicates that if implemented with appropriate constitutive models, SPH method can be used to model large‐deformation problems with high fidelity. Copyright © 2013 John Wiley & Sons, Ltd.     

A multi‐mechanism constitutive model for plastic adaption under cyclic loading

As is well known, granular soils under cyclic loading dissipate a large amount of energy and accumulate large irreversible strains. Usually, with time, this second effect reduces and the accumulation rate decreases with the number of cycles until obtaining a sort of ideal stationary cyclic state at which ratcheting disappears. In this paper, only this ideal state is taken into consideration and simulated by means of a multi‐mechanism constitutive model for plastic adaptation. For this purpose, the concept of cycle is discussed, many different categories of cyclic stress/strain paths are considered and some theoretical issues concerning both the flow and the strain‐hardening rules are tackled. Even though the paper focuses on soil behaviour, the conclusions can be extended to all materials exhibiting ratcheting due to volumetric behaviour.Copyright © 2013 John Wiley & Sons, Ltd.     

一种旋转塑性势面模型及非关联塑性流动法则

为了较好地描述软土塑性应变发展规律,提出了一种改进的塑性流动模型。该模型采用了与屈服函数形式相同,但具有一定倾角γ的塑性势函数。土体在变形过程中,塑性流动方向会依赖于塑性势面的旋转而变化,直至达到破坏状态。通过对常规三轴试验结果的分析可以发现:在剪切过程中,塑性势面旋转角的初值γ₀与终值γ_d较为稳定,不受围压变化影响。在此试验观察基础上,引入了归一化的旋转角参数γ’以及描述土体应力状态的参数ξ,在采用蛋形势函数的情况下二者具有良好的分段线性关系。利用该关系,建立了改进的塑性流动法则,只需要2个额外的模型参数。对所提出的塑性流动模型进行了验证,计算结果表明该模型能较好地反映塑性应变的变化趋势。     

Multiscale approach to geo‐composite cellular structures subjected to rock impacts

Geo‐composite cellular structures are an efficient technological solution for various applications in civil engineering. This type of structure is particularly well adapted to resisting rockfalls and can act as a defensive structure. However, the design of such structures is for the most part empirically based; this lack of research‐based design stagnates optimization and advanced development. In this paper, the mechanical behaviour of a geo‐composite cellular structure is investigated using a multi‐scale approach, from the individual cell made up of an assembly of rocky particles contained in a wire netting cage to the entire structure composed of a regular array of cells. Based on discrete modelling of both the cell and structure scales, a computational tool has been developed for design purposes. Copyright © 2007 John Wiley & Sons, Ltd.     

Efficient solution algorithms for multiaxial probabilistic elasto‐plastic constitutive simulations of soils

A Fokker‐Planck‐Kolmogorov (FPK) equation approach has recently been developed to probabilistically solve any elastic‐plastic constitutive equation with uncertain material parameters by transforming the nonlinear, stochastic constitutive rate equation into a linear, deterministic partial differential equation (PDE) and thereby simplifying the numerical solution process. For an uniaxial problem, conventional numerical techniques, such as the finite difference or finite element methods, may be used to solve the resulting univariate FPK PDE. However, for a multiaxial problem, an efficient algorithm is necessary for tractability of the numerical solution of the multivariate FPK PDE. In this paper, computationally efficient algorithms, based on a Fourier spectral approach, are presented for solving FPK PDEs in (stress) space and (pseudo) time, having space‐independent but time‐dependent coefficients and both space‐ and time‐dependent coefficients, that commonly arise in probabilistic elasto‐plasticity. The algorithms are illustrated by probabilistically simulating 2 common laboratory constitutive experiments in geotechnical engineering, namely, the unconfined compression test and the unconsolidated undrained triaxial compression test.