Simulation of yielding and stress–stain behavior of shanghai soft clay

In this paper, a simple bounding surface plasticity model is used to reproduce the yielding and stress–strain behavior of the structured soft clay found at Shanghai of China. A series of undrained triaxial tests and drained stress probe tests under isotropic and anisotropic consolidation modes were performed on undisturbed samples of Shanghai soft clay to study the yielding characteristics. The degradation of the clay structure is modeled with an internal variable that allows the size of the bounding surface to decay with accumulated plastic strain. An anisotropic tensor and rotational hardening law are introduced to reflect the initial anisotropy and the evolution of anisotropy. Combined with the isotropic hardening rule, the rotational hardening rule and the degradation law are incorporated into the bounding surface formulation with an associated flow rule. Validity of the model is verified by the undrained isotropic and anisotropic triaxial test and drained stress probe test results for Shanghai soft clay. The effects of stress anisotropy and loss of structure are well captured by the model.     

Constitutive modelling of granular materials using a contact normal-based fabric tensor

This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials.

     

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.     

软黏土的各向异性临界状态模型

剑桥模型只适用于正常固结软黏土,不能描述不等向固结土的应力-应变行为的各向异性特性。基于剑桥模型,在其椭圆屈服面中引入各向异性张量和一个形状参数,建立了一个各向异性屈服面,提出了一个适用于等向和不等向固结软黏土的本构模型。各向异性张量的初始值由初始固结应力状态确定,其演化过程由一个与塑性剪应变和塑性体应变都有关的硬化法则描述。形状参数的引入保证了各向异性屈服面的灵活性和适应性。通过对Boston Blue黏土、高岭土和Otaniemi黏土的三轴试验结果的模拟,验证了模型的模拟能力。     

基于热力学的临界状态模型的返回映射算法及实现

相对于其他临界状态模型,基于热力学的临界状态（TCS）模型不需要引入塑性势假设,能够自动满足热力学定律。通过对TCS的修正,使其能够模拟初始K0固结的影响,利用了返回映射算法进行了TCS模型的ABAQUS本构二次开发,通过与ABAQUS内嵌的修正剑桥模型（MCC）计算结果的对比,证实了该程序的可靠性。对模型中的2个控制屈服面形状的参数进行了讨论,分析了它们对应力-应变关系与剪胀关系的影响,修改TCS模型参数可以实现非椭圆的屈服面,进而拓展了模型的适用范围,不同的参数对于屈服面形状和大小的影响也有所不同。同时还比较了采用考虑K0固结及旋转硬化的TCS模型与不考虑K0固结及旋转硬化的MCC模型所得到的应力-应变关系与剪胀关系的差异。对于真实土体,K0固结及旋转硬化是土体的基本力学特性。最后证明了TCS模型较MCC模型（模型中没有考虑旋转硬化与K0固结对硬化规律的影响）可以很好地模拟该特性,所以更加有效。     

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).     

Implicit and explicit integration schemes in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay

Two integration algorithms, namely the implicit return mapping and explicit sub-stepping schemes, are adopted in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay and are implemented into finite element code. The model is a representative of a series of bounding surface models that have typical characteristics, including isotropic and kinematic hardening rules and a rotational bounding surface to capture complex but important cyclic behaviours of soils, such as cyclic shakedown and degradation. However, there is no explicit current yield surface in the model to which the conventional implicit algorithm returns the stress state back or the sub-stepping integration corrects the drift of the stress state. Hence, necessary modifications have been made for both of the integration schemes. First, the image stress point is mapped or corrected to the bounding surface instead of mapping back or correcting the stress state to the yield surface. Second, the unloading–loading criterion is checked to determine the image stress point rather than checking the yield criterion after giving the trial stress state in a conventional way. Comparative studies on the accuracy, stability and efficiency of the two integration schemes are conducted not only at the element level but also in solving boundary value problems of monotonic and cyclic bearing behaviours of rigid footings on saturated clay. For smaller strain increments, there is no significant difference in the accuracy between the two integration schemes, but the explicit integration shows a higher efficiency and accuracy. For relatively larger increments, the implicit return mapping algorithm presents good accuracy and more robustness, while the sub-stepping algorithm shows deteriorating accuracy and suffers the convergence problem. With the tolerance used in the present model, the bearing capacity of the rigid footing predicted by the return mapping algorithm is closer to the available analytical and numerical solutions, while the bearing capacity predicted by the sub-stepping algorithm shows a marginal increase.     

A critical state model for overconsolidated structured clays

This paper presents a generalised critical state model with the bounding surface theory for simulating the stress–strain behaviour of overconsolidated structured clays. The model is formulated based on the framework of the Structured Cam Clay (SCC) model and is designated as the Modified Structured Cam Clay with Bounding Surface Theory (MSCC-B) model. The hardening and destructuring processes for structured clays in the overconsolidated state can be described by the proposed model. The image stress point defined by the radial mapping rule is used to determine the plastic hardening modulus, which varies along loading paths. A new proposed parameter h, which depends on the material characteristics, is introduced into the plastic hardening modulus equation to take the soil behaviour into account in the overconsolidated state. The MSCC-B model is finally evaluated in light of the model performance by comparisons with the measured data of both naturally and artificially structured clays under compression and shearing tests. From the comparisons, it is found that the MSCC-B model gives reasonable good simulations of mechanical response of structured clays in both drained and undrained conditions. With its simplicity and performance, the MSCC-B model is regarded as a practical geotechnical model for implementation in numerical analysis.     

考虑饱和黏土初始各向异性的循环边界面塑性模型

天然土体的初始各向异性通常可对其后继循环特性产生显著影响。现有考虑循环载荷作用的土体弹塑性模型,往往采用类似修正剑桥模型的椭圆形屈服面,已有研究表明,该椭圆形屈服面因其拉伸弹性区域偏大,针对天然K0固结状态的土体,其计算精度较差。基于新近提出的广义各向同性硬化准则,在边界面方程中引入初始各向异性张量,并采用空间滑动面破坏准则（SMP）的变换应力法,建立了能考虑饱和黏土初始各向异性的循环边界面塑性模型。分别针对等压和偏压固结的饱和黏土静、动三轴试验进行模拟,结果表明,该模型能合理反映土体的初始各向异性及其后继循环动力特性。     

循环荷载下砂土的剪切硬化边界面本构模型

基于临界状态土力学框架,建立了一个适用于往返循环荷载作用的砂土边界面本构模型。采用无纯弹性域假设,认为受到反向荷载的瞬时土体就产生塑性变形,砂土的弹性区域退化为一个点。屈服面为倒子弹头型,由于砂土孔隙比与压力之间不存在惟一对应的关系,使得屈服面大小无法与体积应变直接耦合,故采用塑性偏应变而不是剑桥模型那种塑性体应变作为硬化参数。流动法则采用加入状态参数的修正的Rowe应力剪胀关系,体现了依赖状态的剪胀思想。屈服面大小的比值 反映了塑性模量的演化,并推导了 的表达式。只用1套参数,该模型就能合理地模拟砂土在不同密度和固结压力下循环荷载的应力-应变关系曲线。     

考虑循环载荷下饱和黏土软化的损伤边界面模型研究

研究表明,循环载荷作用下饱和黏土将发生软化,其机制主要有两个：一是孔压的积累;二是土体原有结构的不断损伤和新结构的不断重塑。针对上述机制,基于广义各向同性硬化准则建立了考虑饱和黏土循环软化的损伤单面模型。该模型在有效应力空间中引入损伤变量,表征土体结构的损伤和重塑程度,在连续的循环加载下,损伤不断累积,边界面则随着损伤的累积不断收缩,以模拟饱和黏土刚度和强度的软化;以应力反向点作为边界面的广义各向同性硬化中心和映射法则的映射中心,灵活地选择塑性模量的插值公式以模拟塑性变形和孔压的累积以及应力-应变的滞回特性。应用该模型对不排水循环三轴试验进行模拟,并且考查了循环周次、循环应力水平和固结历史对饱和黏土循环软化特性的影响,并与相关试验比较,验证了模型的有效性。     

A fully coupled simple model for unsaturated soils

Different phenomena influence the strength and volumetric behavior of unsaturated soils. Among the most important are suction hardening, hydraulic hysteresis, and the influence of volumetric strain on the soil-water retention curves. Fully coupled hydro-mechanical models require including all three phenomena in their constitutive relationships. Among these phenomena, suction hardening is the most influencing as it determines the apparent preconsolidation stress, the position of the loading-collapse yield surface, and the shift of both the isotropic consolidation and the critical state lines. In this paper, a fully coupled hydro mechanical model is presented. It is based on the modified Cam-Clay model but includes a yield surface with anisotropic hardening that takes account of the shift of the critical state line with suction. For highly overconsolidated materials, the sub-loading surface concept has been included in order to increase the precision of the model for these materials. The shift of the retention curves produced by volumetric strains is simulated using a hydraulic model based on the grain and current pore size distribution of the soil.     

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

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

适用于砂土循环加载分析的边界面塑性模型

基于临界状态土力学框架,建立了一个适用于砂土排水循环加载的边界面塑性模型。采用了考虑虚拟峰值应力比的偏应变硬化准则,初始加载阶段应力点位于边界面上,反向加载阶段以历史最大屈服面作为边界面,同时实现了对密砂软化现象的模拟和对历史所受最大应力的记忆。边界面采用修正的椭圆形,引入考虑密度与应力水平的状态相关剪胀函数,采用非相关联流动法则和以应力反向点作为映射中心的径向映射准则。模型仅有10个参数,通过常规三轴试验即可确定,并且使用一套参数可以模拟不同围压、密度的单调和循环加载情况。分别对饱和砂土的单调、循环排水三轴试验进行模拟,结果表明,该模型能够合理地反映饱和砂土排水条件下的应力-应变特性。     

A two-surface plasticity model for stiff clay

This paper presents a constitutive model for describing some important features of the behavior of natural stiff clay evidenced experimentally such as the limited elastic zone, the presence of strain hardening and softening, and the smooth transition from elastic behavior to a plastic one. The model, namely ACC-2, is an adapted Modified Cam Clay model with two yield surfaces: similarly to bounding surface plasticity theory, an additional yield surface—namely Inner yield surface—was adopted to account for the plastic behavior inside the conventional yield surface. A progressive plastic hardening mechanism was introduced with a combined volumetric-deviatoric hardening law associated with the Inner yield surface, enabling the plastic modulus to vary smoothly during loading paths. The main feature of the proposed model is that its constitutive equations can be simply formulated based on the consistency condition for the Inner yield surface, so that it can be efficiently implemented in a finite element code using a stress integration scheme similar to that of the Modified Cam Clay model. Furthermore, it is proved to be an appropriate model for natural stiff clay: the simulations of a set of tests along different mechanical loading paths on natural Boom Clay show good agreement with the experimental results.