The subloading isotropic plasticity as a variable modulus model

The subloading concept is an extension of mathematical plasticity which defines an internal surface to the conventional yield surface. It is indeed a versatile approach, especially for the modelling of soils under quasi-static cycles with smooth transitions from pure elastic to elastoplastic behaviour. For the case of isotropic hardening models, this paper demonstrates that the subloading isotropic plasticity is equivalent to a variable modulus approach and therefore a simpler and equivalent methodology can be adopted instead. In addition to demonstrating this equivalence, an alternative formulation that was presented elsewhere and that uses only one surface is briefly discussed. The alternative formulation can then be easily applied to popular models for soils such as the Cam clay model. Finally, some numerical predictions are presented in order to illustrate the capabilities of the subloading isotropic plasticity and the corresponding variable modulus approach.     

加筋土坡临界高度的研究

假定破裂面为一过坡脚的斜平面,分别以传统塑性理论和广义塑性理论为基础,导出了加筋土坡临界高度的计算公式。通过与前人试验结果的比较发现,虽然基于广义塑性理论极限法求得的加筋土坡临界高度值要比基于传统塑性理论极限法求得的临界高度值略偏大,且更接近于试验实测值,但二者在工程中均是可靠的,都可用于加筋土坡设计时的参考。     

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.     

Implicit numerical integration for a kinematic hardening soil plasticity model

Soil models based on kinematic hardening together with elements of bounding surface plasticity, provide a means of introducing some memory of recent history and stiffness variation in the predicted response of soils. Such models provide an improvement on simple elasto‐plastic models in describing soil behaviour under non‐monotonic loading. Routine use of such models requires robust numerical integration schemes. Explicit integration of highly non‐linear models requires extremely small steps in order to guarantee convergence. Here, a fully implicit scheme is presented for a simple kinematic hardening extension of the Cam clay soil model. The algorithm is based on the operator split methodology and the implicit Euler backward integration scheme is proposed to integrate the rate form of the constitutive relations. This algorithm maintains a quadratic rate of asymptotic convergence when used with a Newton–Raphson iterative procedure. Various strain‐driven axisymmetric triaxial paths are simulated in order to demonstrate the efficiency and good performance of the proposed algorithm. Copyright © 2001 John Wiley & Sons, Ltd.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

An extended finite element framework for slow‐rate frictional faulting with bulk plasticity and variable friction

We present an extended finite element (FE) approach for the simulation of slow‐rate frictional faulting in geologic media incorporating bulk plasticity and variable friction. The method allows the fault to pass through the interior of FEs without remeshing. The extended FE algorithm for frictional faulting, advocated in two recent articles, emanates from a variational equation formulated in terms of the relative displacement on the fault. In the present paper we consider the combined effects of bulk plasticity and variable friction in a two‐dimensional plane strain setting. Bulk plasticity is localized to the fault tip and could potentially be used as a predictor for the initiation and propagation of new faults. We utilize a variable velocity‐ and state‐dependent friction, known as the Dieterich–Ruina or ‘slowness’ law, formulated in a slip‐weakening format. The slip‐weakening/variable friction model is then time‐integrated according to the generalized trapezoidal rule. We present numerical examples demonstrating the convergence properties of a global Newton‐based iterative scheme, as well as illustrate some interesting properties of the variable friction model. Copyright © 2009 John Wiley & Sons, Ltd.     

Kinematic hardening plasticity formulation of small strain behaviour of soils

An objective of this paper is to demonstrate that the small strain model developed by the authors can be incorporated into the conventional kinematic hardening plasticity framework to predict pre‐failure defor mations. The constitutive model described in this paper is constituted by three elliptical yield surfaces in triaxial stress space. Two inner surfaces are rotated ellipses of the same shape, representing the boundaries of the linear elastic and small strain regions, while the third surface is the modified Cam clay large‐scale yield surface. Within the linear elastic region, the soil behaviour is elastic with cross‐coupling between the shear and volumetric stress–strain components. Within the small strain region, the soil behaviour is elasto‐plastic, described by the kinematic hardening rule with an infinite number of loading surfaces defined by the incremental energy criterion. Within the large‐scale yield surface, the soil behaviour is elasto‐plastic, described by kinematic and isotropic hardening of the small strain region boundary. Since the yield surfaces have different shapes, the uniqueness of the plastic loading condition imposes a restriction on the ratio between their semi‐diameters. The model requires 12 parameters, which can be determined from a single consolidated undrained triaxial compression test. Copyright © 2000 John Wiley & Sons, Ltd.     

A METHOD TO CORRECT YIELD SURFACE DRIFT IN SOIL PLASTICITY UNDER MIXED CONTROL AND EXPLICIT INTEGRATION

When applying an explicit integration algorithm in e.g. soil plasticity, the predicted stress point at the end of an elastoplastic increment of loading might not be situated on the updated current yield surface. This so-called yield surface drift could generally be held under control by using small integration steps. Another possibility, when circumstances might demand larger steps, is to adopt a drift correction method. In this paper, a drift correction method for mixed control in soil plasticity, under drained as well as undrained conditions, is proposed. By simulating triaxial tests in a Constitutive Driver, the capability and efficiency of this correction method, under different choices of implementation, have been analysed. It was concluded that the proposed drift correction method, for quite marginal additional computational cost, was able to correct successfully for yield surface drift giving results in close agreement to those obtained with a very large number of integration steps. © 1997 by John Wiley & Sons, Ltd.     

ACMEG‐TS: A constitutive model for unsaturated soils under non‐isothermal conditions

This paper introduces an unconventional constitutive model for soils, which deals with a unified thermo‐mechanical modelling for unsaturated soils. The relevant temperature and suction effects are studied in light of elasto‐plasticity. A generalized effective stress framework is adopted, which includes a number of intrinsic thermo‐hydro‐mechanical connections, to represent the stress state in the soil. Two coupled constitutive aspects are used to fully describe the non‐isothermal behaviour. The mechanical constitutive part is built on the concepts of bounding surface theory and multi‐mechanism plasticity, whereas water retention characteristics are described using elasto‐plasticity to reproduce the hysteretic response and the effect of temperature and dry density on retention properties. The theoretical formulation is supported by comparisons with experimental results on two compacted clays. Copyright © 2008 John Wiley & Sons, Ltd.     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 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.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

3D subsidence analyses above gas reservoirs accounting for an unconventional plasticity model

The coupled hydro‐mechanical state in soils coming from consolidation/subsidence processes and undergoing plasticity phenomena is here evaluated by means of the subloading surface model. The most important feature of this theory is the abolition of the distinction between the elastic and plastic domains, as it happens in the conventional elastoplastic models. This means that plastic deformations are generated whenever there is a change in stress and a smoother elastoplastic transition is produced. The plasticity algorithm has been implemented in the PLASCON3D FE code, coupling hydro‐thermo‐mechanical fields within a saturated (locally partially saturated) porous medium subjected to external loads and water/gas withdrawals from deep layers (aquifers/reservoirs). The 3D model has been first calibrated and validated against examples taken from the literature, and then subsidence analyses at regional scales due to gas extractions have been developed to predict the evolution of settlements and pore pressure in soils for long‐term scenarios. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling the behavior of sandstone based on generalized plasticity concept

With the concept of generalized plasticity, a constitutive model for describing the deformation behavior of sandstone is proposed in this paper. This proposed model is characterized by the following features: (i) nonlinear elasticity under hydrostatic and shear loading; (ii) associated flow rule for pre‐peak simulation; (iii) substantial plastic deformation during shear loading; and (iv) significant shear dilation and distortion prior to the failure state. This model requires 10 material parameters, including three for elasticity and seven for plasticity. All of the parameters can be determined, in a straightforward manner, by the suggested procedures. The proposed model has been validated by comparing the triaxial test results of the Mushan sandstone under different hydrostatic stress, different stress paths, and cyclic loading condition. It is also versatile in simulating the deformation behaviors of two other sandstones. Upon slight modification of the model, the post‐peak behavior of sandstone can be reasonably predicted using proposed model. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

上限法在加筋土结构物稳定性分析中的应用

将极限分析上限法用于加筋土结构物稳定性分析的关键是滑动面上内力功如何计算的问题。在总结了基于传统塑性理论上限法内力功的计算方法的基础上,将基于广义塑性理论的上限法应用到杨俊杰的内力功计算方法中。对比基于传统和广义塑性理论的上限解及离心试验结果,认为基于传统塑性理论的上限法应用简单,计算结果基本可靠;基于广义塑性理论的上限法虽然理论上比较完善,但当滑动面为曲线时,应用繁琐。     

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

Validation of a three‐dimensional constitutive model for nonlinear site response and soil‐structure interaction analyses using centrifuge test data

The capability of a bounding surface plasticity model with a vanishing elastic region to capture the multiaxial dynamic hysteretic responses of soil deposits under broadband (eg, earthquake) excitations is explored by using data from centrifuge tests. The said model was proposed by Borja and Amies in 1994 (J. Geotech. Eng., 120, 6, 1051‐1070), which is theoretically capable of representing nonlinear soil behavior in a multiaxial setting. This is an important capability that is required for exploring and quantifying site topography, soil stratigraphy, and kinematic effects in ground motion and soil‐structure interaction analyses. Results obtained herein indicate that the model can accurately predict key response data recorded during centrifuge tests on embedded specimens—including soil pressures and bending strains for structural walls, structures' racking displacements, and surface settlements—under both low‐ and high‐amplitude seismic input motions, which was achieved after performing only a basic material parameter calibration procedure. Comparisons are also made with results obtained using equivalent linear models and a well‐known pressure‐dependent multisurface plasticity model, which suggested that the present model is generally more accurate. The numerical convergence behavior of the model in nonlinear equilibrium iterations is also explored for a variety of numerical implementation and model parameter options. To facilitate broader use by researchers and practicing engineers alike, the model is implemented as a “user material” in ABAQUS Standard for implicit time stepping.     

A bounding surface elasto‐viscoplastic constitutive model for non‐isothermal cyclic analysis of asphaltic materials

An elasto‐viscoplastic constitutive model for asphaltic materials is presented within the context of bounding surface plasticity theory, taking into account the effects of the stress state, void binder degree of saturation, temperature and strain rate on the material behaviour. A stress state dependent non‐linear elasticity model is introduced to represent time‐independent recoverable portion of the deformation. The consistent visco‐plasticity framework is utilised to capture the rate‐dependent, non‐recoverable strain components. The material parameters introduced in the model are identified, and their determination from conventional laboratory tests is discussed. The capability of the model to reproduce experimentally observed response of asphaltic materials is demonstrated through numerical simulations of several laboratory test data from the literature. Copyright © 2016 John Wiley & Sons, Ltd.