Kinematic hardening of a porous limestone

Summary A concept for a kinematic hardening yield surface in stress space for Cordova Cream limestone (Austin Chalk) developed by Celle and Cheatham (1981) has been improved using Ziegler's modification of Prager's hardening rule (Ziegler, 1959). Data to date agree with the formulated concepts. It is shown how kinematic hardening can be used to approximate the yield surface for a wide range of stress states past the initial yield surface. The particular difficulty of identifying the yield surface under conditions of unloading or extension is noted. A yield condition and hardening rule which account for the strain induced anisotropy in Cordova Cream Limestone were developed. Although the actual yield surface appears to involve some change of size and shape, it is concluded that true kinematic hardening provides a basis for engineering calculations.     

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.     

Rotational kinematic hardening model for three‐dimensional stress reversals in sand

A kinematic hardening mechanism has previously been proposed to capture the behavior of soil during large stress reversals in the triaxial plane. This mechanism is now extended to the principal stress space. It incorporates rotation and intersection of yield surfaces to achieve a consistent and physically rational fit with experimentally observed soil behavior during large three‐dimensional stress reversals. An existing elasto‐plastic model with isotropic hardening is used as the basic framework to which the rotational kinematic hardening mechanism has been added. The new combined model preserves the behavior of the isotropic hardening model under monotonic loading conditions, and the extension from isotropic to rotational kinematic hardening under three‐dimensional conditions is accomplished without introducing new material parameters. The framework of the model is presented here with some comparisons between theoretical and experimental directions of strain increment vectors to indicate the potential of the model. Copyright © 2008 John Wiley & Sons, Ltd.     

基于旋转运动硬化理论的修正剑桥模型的改进

修正剑桥模型是最早建立和得到广泛承认的经典土体弹塑性模型之一,但不能模拟应力路径转折时土体的应力-应变特性以及应力引起的各向异性。将旋转运动硬化理论引入到剑桥模型中,给出了椭圆屈服面的旋转运动硬化机制,在不增加任何模型参数的情况下,把等向硬化的修正剑桥模型扩展为旋转运动硬化模型。扩展的新模型既保留了单调加载时的等向硬化的特性,又能反映应力路径转折时土体的本构特性与应力诱发的各向异性,初步验证了模型的有效性。     

Effects of principal stress rotation on wave-induced soil response in a poro-elastoplastic sandy seabed

In this study, a constitutive model is developed in order to investigate wave–seabed interactions. This model takes into account the impact of principal stress rotation (PSR) and is based on the generalized plasticity theory, in which plastic strain generated by PSR is considered an additional item in the constitutive relationship of soil. The normalized loading direction and plastic flow direction were determined based on the stress tensor invariant. Comparisons between the present model and previous Hollow Cylinder Apparatus tests and geotechnical centrifugal wave tests show good agreement. Numerical results show the effects of PSR on predictions of liquefaction potential due to: (a) the cumulative impact of plastic strain in the seafloor and (b) the buildup of pore pressure. Parametric study shows that the model parameters, including the wave and seabed parameters, have significant effects on the wave-induced soil liquefaction.

     

非饱和土的本构模型研究

引用平均土骨架应力的概念,研究推导出非饱和土的刚度参数随吸力变化而变化的关系式,进而推导得到用平均土骨架应力表述的非饱和土LC屈服面函数以及硬化规律。从土力学原理推导,得到土样由于在净应力和吸力作用下产生体积变形引起土样饱和度变化的关系式。由平均土骨架应力推广,得到三轴应力状态的椭圆屈服函数,这一非饱和土本构模型的优点在于考虑了应力作用后土样饱和度的变化,通过对已有试验数据的初步验证,表明提出的非饱和土本构模型的合理性和适用性。     

Kinematic hardening of soft clay in simple shear

In a separate paper, the authors have proposed a normalized, non-degrading form of the shear stress–shear strain relationship for undrained, cyclic simple shear of soft clay. This relationship is described in the present paper, and it is seen to include a single fatigue parameter—the mean effective stress. Application of the relationship therefore requires knowledge of the history of the mean effective stress during any loading history. The present paper proposes an effective stress path model which may be used for prediction of this history. The model is developed within the framework of bounding surface kinematic and isotropic hardening plasticity. It incorporates an isotropic hardening bounding surface, and a kinematic hardening yield surface, in which the elastic region vanishes, and so the yield surface reduces to the stress point. The normalized shear stress–shear strain relationship, developed on the basis of Iwan's model, is used to establish the shape of the cap of the bounding surface. A new translation rule is also incorporated in the model, allowing improved prediction of stress path development within the bounding surface during regular or irregular cyclic loading. Use of the proposed model to simulate the behaviour of soft clay in laboratory undrained cyclic simple shear tests shows excellent qualitative agreement, with most of the major features of the actual behaviour being predicted.     

An elastoplastic model based on the shakedown concept for flexible pavements unbound granular materials

Nowadays, the problem of rutting of flexible pavements linked to permanent deformations occurring in the unbound layers is taken into account only by mechanistic empirical formulas. Finite element modelling of realistic boundary value problems with incremental rheological models will lead to unrealistic calculation time for large cycle numbers. The objective of the authors is to present a simplified model which can be used to model the flexible pavements rutting with the finite elements framework. This method is based on the shakedown theory developed by Zarka which is usually associated to materials like steels. It has been adapted for granular materials by introducing a yield surface taking into account the mean stress influence on the mechanical behaviour and a dependency of the hardening modulus with the stress state. The Drucker–Prager yield surface has been used with a non‐associated flow rule. Comparisons with repeated load triaxial tests carried out on a subgrade soil have been done. These comparisons underline the capabilities of the model to take into account the cyclic behaviour of unbound materials for roads. Finally, a discussion, dealing with the use of the simplified method within a finite element modelling of a full‐scale experiment, is presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulating the effects of induced anisotropy on liquefaction potential using a new constitutive model

The effects of induced anisotropy on the undrained behaviour of very loose and saturated sands have been a subject of intensive investigation, both experimentally and theoretically, by several authors in the past few years. This paper proposes an original constitutive model well‐adapted to simulate the behaviour of sands subject to complex stress histories, in particular, the preloading cycle along the classical drained stress path in compression. The developed model belongs to the family of critical state models. Its construction is based on a few theoretical concepts taken from the theory of ‘Bounding Surface Plasticity’ developed among others by Y. Dafalias and Popov (1975), the ‘Clay And Sand Model’ (CASM) of H. Yu (2006), the CJS model (B. Cambou and K. Jafari (1988)) and the hyperelastic isotropic model of P. Lade (1987). To accurately simulate volume changes, which represent a key element in soil behaviour, a state‐dependent dilatancy rule is proposed, which can account for the influences of stress and void ratio. The current void ratio depends implicitly on the irreversible strains already accumulated hence the strain history. A kinematic hardening is combined with an isotropic hardening, involving rotation and distortion of the bounding surface, in order to capture correctly the experimental observations. Comparisons of experimental results to numerical simulations show that the model is able to simulate with a good precision the major trends of undrained responses of loose and presheared sands. It predicts correctly rapid static liquefaction at small or null drained preloading, as well as the progressive transition to a completely stable behaviour typical of dense sands, while the sample is loose in reality. At intermediate to large amplitudes of preloadings, the model also predicts correctly the temporary stage of instability when the deviatoric stress decreases slightly before rising up again. Copyright © 2013 John Wiley & Sons, Ltd.     

Predictions of large stress reversals in true triaxial tests on cross‐anisotropic sand

A rotational kinematic hardening constitutive model with the capability of predicting the behavior of soil during three‐dimensional stress reversals has been developed. An existing elasto‐plastic constitutive model, the Single Hardening Model, utilizing isotropic hardening serves as the basic framework in these formulations. The framework of the kinematic hardening model was discussed in a companion paper. The previously proposed cross‐anisotropic Single Hardening Model is added to the present kinematic hardening mechanism to capture inherent anisotropy of sands in addition to the stress reversals. This model involves 13 parameters, which can be determined from simple laboratory experiments, such as isotropic compression, drained triaxial compression and triaxial extension tests. The results from a series of true triaxial tests with large three‐dimensional stress reversals performed on medium dense cross‐anisotropic Santa Monica Beach sand are employed for comparison with predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

基于硬化土模型的小应变本构模型研究

大量的工程现场监测结果表明,城市中强支护隧道施工过程中其周围大部分土体仍处于小应变状态,小应变情况下土体具有显著的高模量和非线性特性。土体小应变刚度特征和应力路径相关性是准确分析土与隧道相互作用的重要因素。虽然采用双刚度的硬化土模型(hardening soil model)能够基本反映应力路径的影响,但其卸载再加载模量并没有与应力-应变水平相关,不能模拟小应变情况下土体模量的高度非线性。为此,结合小应变刚度理论对硬化土模型的卸载再加载模量进行了改进,使之与应力-应变水平相关,并且考虑了土体卸载抗剪强度指标的变化以及侧向卸荷应力路径下不同的模量。通过与土体的应力路径试验结果进行对比,证明了模型的合理性。     

敏感环境下基坑数值分析中土体本构模型的选择

数值分析已成为敏感环境下基坑工程分析的最重要手段,其关键是选择合适的土体本构模型和计算参数。在分析了岩土数值分析中常用土体本构模型特点的基础上,通过算例较系统地对比了各类模型在基坑开挖数值分析中的适用性。敏感环境下的基坑工程需重点关注墙后土体的变形,从满足工程需要和方便实用的角度出发,建议采用能考虑黏土的塑性和应变硬化特征、能区分加荷和卸荷且刚度依赖于应力水平的硬化类弹塑性模型,如MCC模型和HS模型进行分析。具体工程实例的分析,表明了硬化类弹塑性模型在敏感环境下基坑开挖数值分析中的适用性。     

Constitutive modeling of dense gravelly soils subjected to cyclic loading

A constitutive model for dense gravelly soils was developed to reproduce their responses under cyclic loadings. Its application aims at nonlinear dynamic analyses of earth structures involving gravelly soils, such as rockfill dams and railroad ballasts. The framework of generalized plasticity was modified to incorporate the concept of stress distance for better simulation of unloading and reloading responses. It was then combined with the theory of critical state soil mechanics to develop the constitutive model. The model has the following important features: unified simulation of particle breakage through translating critical state line, smooth transition from unloading to reloading in the stress space, and proper modeling of cyclic hysteresis, cyclic densification, and cyclic hardening of dense gravelly soils. Most of the model parameters can be obtained through simple calculation using conventional triaxial test results, and their calibration process was discussed. The model was used to simulate the cyclic responses of three gravelly soils with satisfactory accuracy. Copyright © 2014 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.     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

结构性土UH模型

在描述重塑土的统一硬化模型（UH模型）基础上,以动态地移动正常固结线（MNCL）代替静态的正常固结线（NCL）作为参考线来确定参考应力,将土的结构性衰减体现在MNCL的演化中,从而把UH模型扩展为可考虑土结构性的结构性土UH模型。结构性土UH模型可以合理地并连续光滑地描述结构性土的等向压缩规律、应变硬化、应变软化、剪切体积收缩、剪切体积膨胀、不排水剪切减压软化以及结构性和密度耦合影响等力学特性。相比于UH模型,结构性土UH模型新增3个参数,分别用于描述结构性的程度、衰减速度、以及塑性流动规律。通过3种天然土的试验结果与模型预测对比验证了结构性土UH 模型的合理性。     

锁儿头滑坡滑带土不同含水率大剪试验研究

锁儿头滑坡属特大型断层破碎带滑坡。该滑坡滑带土主要由强风化的中、上志留统千枚岩、板岩等构成,滑带土中粒径大于2 mm的颗粒含量约为30%,在滑带土颗粒组成中占有重要地位,其强度特征对滑坡体的稳定性分析具有重要影 响。为此,利用室内大型直剪仪对不同含水率下锁儿头滑坡滑带土的剪切特性进行试验研究,结果表明,在控制干密度为1.7 g/cm3情况下,不同含水率滑带土试样的剪应力-应变关系均呈现出弱应变硬化特征,而且法向压力越大,应变硬化特征越明显;在剪应变ε为10%时,相同含水率下滑带土剪切模量G0.1随法向压力的增大而增大,同一法向压力下的剪切模量G0.1与含水率间呈负线性相关关系;利用邓肯-张双曲线模型对锁儿头滑坡滑带土不同含水率下大剪试验剪应力-应变关系进行数据拟合,拟合效果良好,但该模型不能反映土的应变软化性质。     

循环荷载下考虑主应力旋转的软土力学响应研究

目前用于轨道路堤分析计算的参数多是通过三轴试验获得,而列车荷载经过时土单元应力路径与循环三轴试验加载路径的显著差异可能导致预测失真而引发工程问题。针对某城市轨道进行了三维计算,并分析了移动列车荷载下土单元应力幅值、循环周数、主应力轴旋转的分布规律。在该基础上,采用空心圆柱循环扭剪试验对该复杂应力路径进行了模拟,研究了饱和软黏土的孔压及变形累积特性。结果表明,列车荷载经过时地基土单元大主应力将在(－90°, 90°)内发生旋转;该主应力轴旋转将显著促进软黏土孔压和应变的累积,竖向应力幅值为15 kPa时,循环扭剪试验产生的孔压值比循环三轴试验高77%,累积应变增大了近50%;随着循环应力水平提高,二者累积孔压及累积应变的差值还会进一步增大,甚至出现循环三轴下仅产生较小应变而循环扭剪下已破坏的本质性差异。     

An elastoplastic model with combined isotropic–kinematic hardening to predict the cyclic behavior of stiff clays

This paper presents a kinematic hardening model for describing some important features of natural stiff clays under cyclic loading conditions, such as closed hysteretic loops, smooth transition from the elastic behavior to the elastoplastic one and changes of the compression slope with loading/unloading loops. The model includes two yield surfaces, an inner surface and a bounding surface. A non-associated flow rule and a kinematic hardening law are proposed for the inner surface. The adopted hardening law enables the plastic modulus to vary smoothly when the kinematic yield surface approaches the bounding surface and ensures at the same time the non-intersection of the two yield surfaces. Furthermore, the first loading, unloading, and reloading stages are treated differently by applying distinct hardening parameters. The main feature of the model is that its constitutive equations can be simply formulated based on the consistency condition for the inner yield surface based on the proposed kinematic hardening law; thereby, this model can be easily implemented in a finite element code using a classic stress integration scheme as for the modified Cam Clay model. The simulation results on the Boom Clay, natural stiff clay, have revealed the relevance of the model: a good agreement has been obtained between simulations and the experimental results from the tests with different stress paths under cyclic loading conditions. In particular, the model can satisfactorily describe the complex case of oedometric conditions where the deviator stress is positive upon loading (compression) but can become negative upon unloading (extension).