A bounding surface plasticity model for unsaturated soil at small strains

Most existing hydromechanical models for unsaturated soils are not able to fully capture the nonlinearity of stress–strain curves at small strains (less than 1%). They cannot therefore, for example, accurately predict ground movements and the performance of many earth structures under working conditions. To tackle this problem, a state‐dependent bounding surface plasticity model has been newly developed. Particularly, the degradation of shear modulus with strain at small strains ranging from 0.001% to 1% is focused. The proposed model is formulated in terms of mean average skeleton stress, deviator stress, suction, specific volume and degree of saturation. Void ratio‐dependent hydraulic hysteresis is coupled with the stress–strain behaviour. Different from other elastoplastic models for unsaturated soils, plastic strains are allowed inside bounding surfaces. In this paper, details of model formulations and calibration procedures of model parameters are presented. To evaluate the capability of the new model, it is applied to simulate a series of triaxial compression tests on compacted unsaturated silt at various suctions. Effects of suction, drying and wetting as well as net stress on unsaturated soil behaviour are well captured. The model shows good predictions of the degradation of shear modulus with strain over a wide range of strains from 0.001% to 1%. Copyright © 2015 John Wiley & Sons, Ltd.     

CASM: a unified state parameter model for clay and sand

The purpose of this paper is to present a simple, unified critical state constitutive model for both clay and sand. The model, called CASM (Clay And Sand Model), is formulated in terms of the state parameter that is defined as the vertical distance between current state (v, p′) and the critical state line in v–ln p′ space. The paper first shows that the standard Cam-clay models (i.e. the original and modified Cam-clay models) can be reformulated in terms of the state parameter. Although the standard Cam-clay models prove to be successful in modelling normally consolidated clays, it is well known that they cannot predict many important features of the behavior of sands and overconsolidated clays. By adopting a general stress ratio-state parameter relation to describe the state boundary surface of soils, it is shown that a simple, unified constitutive model (CASM) can be developed for both clay and sand. It is also demonstrated that the standard Cam-clay yield surfaces can be either recovered or approximated as special cases of the yield locus assumed in CASM. The main feature of the proposed model is that a single set of yield and plastic potential functions has been used to model the behaviour of clay and sand under both drained and undrained loading conditions. In addition, it is shown that the behaviour of overconsolidated clays can also be satisfactorily modelled. Simplicity is a major advantage of the present state parameter model, as only two new material constants need to be introduced when compared with the standard Cam-clay models. © 1998 John Wiley & Sons, Ltd.     

Multiphysics implementation of advanced soil mechanics models

Using multiphysics computer codes has become a useful tool to solve systems of partial differential equations. However, these codes do not always allow for the free introduction of implicitly defined state functions when automatic differentiation is used to compute the iteration matrix. This makes it considerably more difficult to solve geomechanical problems using non-linear constitutive models. This paper proposes a method for overcoming this difficulty based on multiphysics capabilities. The implementation of the well-known Barcelona basic model is described to illustrate the application of the method. For this purpose, without including formulation details addressed by other authors, the fundamentals of its implementation in a finite element code are described. Examples that demonstrate the scope of the proposed methodology are also presented.     

考虑物态变化的粗粒土亚塑性本构模型

建立了可描述粗粒土单调和循环力学特性的一个亚塑性本构模型。基于亚塑性理论的基本框架,引入临界状态参数,建立了一个新的粗粒土亚塑性模型,给出了数学公式及参数确定方法。采用该模型对粗粒土单调和循环加载试验进行了模拟和预测。该模型无需判断加卸载、参数较少、易于三维化,能够较全面地描述单调和循环荷载作用下粗粒土的主要力学特性,如强度与围压的非线性关系,胀缩规律与围压相关、卸载体缩、体变随加载过程累积等主要体变特性等。     

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.     

一类新型的散粒型土体本构理论

详细地介绍了国际上近30才发展起来的具有独特建模思想的一类新型的散粒型土体本构建模理论--亚塑性本构理论。将该理论与经典的弹塑性理论进行比较后发现,亚塑性理论的本构表达式自动隐含了弹塑性理论中相对应的一些基本概念和假定,无需额外引入,因此增加了模型的客观性。以Gudehus-Bauer亚塑性模型为例,分析了模型在临界状态时的本构特征。此外,通过三轴试验和循环剪切试验的亚塑性数值模拟表明,亚塑性本构模型能很好地反映无黏性散粒型土体的非弹性、非线性及剪胀（剪缩）性等主要应力-应变特性。     

土体小应变特性研究中的边界面模型

研究表明，土体小应变特性对准确预测土体的变形起着十分重要的作用，而这一特性很难用传统的本构模型对它加以描述。边界面模型能描述一些经典塑性理论所不能描述的土的真实特性，如何将之进一步加以改进，并有效地应用于土体小应变特性的模拟中是一个有待深入研究的课题。     

武汉黏土Duncan本构模型参数研究

通过常规三轴压缩排水试验,确定了武汉正常固结黏土的Duncan Et-?t模型、Et-Bt模型及常值?t模型参数;通过三轴试样的有限元变形分析,验证了模型参数的合理性与有效性,为武汉地区Duncan本构模型的应用提供了参考。针对Duncan Et-?t模型基于?1-?3曲线为双曲线关系假设,而实际试验曲线常常偏离理想的双曲线关系,造成?t参数拟合的困难,提出利用通过应力水平为70 %与95 %对应点的直线拟合?3 /?1与?3之间线性关系从而确定?t模型参数的改进方法,并将该模型参数用于三轴试样的体积变形预测,经与Duncan Et-Bt、Duncan Et和常值?t模型参数预测结果的对比,发现该方法与试验结果有更好的一致性,并具有简单方便的优点。     

基于细观分析的黏土本构模型

根据黏性土在固结排水条件下的加、卸载变形特征,分析了孔隙水,固体骨架在不同变形阶段的力学响应行为,按照细观力学分析中的自洽方法,建立了在固结排水条件下的黏土损伤本构模型。模型中考虑了孔隙水和固体骨架在加、卸载阶段的不同特性,认为在损伤阶段整体剪切模量的降低是由固体骨架颗粒接触面滑移而引起,与骨架中滑动相的体积百分比和滑动相的剪切模量有关,并给出了求解整体变形模量的解析方法,最后将模型预测与不同初始固结压力,不同应力路径的排水试验结果作了比较,证明该模型是合理的。     

On generalization of constitutive models from two dimensions to three dimensions

In this paper, a study is made of the generalization of constitutive models for geomaterials from two‐dimensional stress and strain states to three‐dimensional stress and strain states. Existing methods of model generalization are reviewed and their deficiencies are highlighted. A new method is proposed based on geometries of the model imprints on two normal planes. Using the proposed method, various three‐dimensional failure criterions suitable for geomaterials are implemented directly into a two‐dimensional model and the generalized model is identical to its original form for the axially symmetric condition. To demonstrate the application of the proposed method, the Modified Cam Clay model is extended using the Matsuoka–Nakai failure criterion. Simulations of soil behaviour for loading in the principal stress space are presented and analysed. Copyright © 2008 John Wiley & Sons, Ltd.     

Sensitive bounding surface constitutive model for structured clay

The main purpose of the paper is to present a relatively simple, yet realistic, constitutive model for simulations of structured sensitive clays. The proposed constitutive model can simulate 1‐D and isotropic consolidation, and drained and undrained shear response of sensitive structured clay. The proposed sensitive bounding surface model is based on concepts from the modified Cam clay model 8 and bounding surface plasticity 27 , with the addition of a simple degradation law. The key material parameters are M, λ, κ, and ν from the modified Cam clay framework, h from the bounding surface framework to model a smoothed elasto‐plastic transition, and ω_v, ω_q, and S_sr to model softening associated with destructuration. The model has separate parameters to model destructuration caused by volumetric strain and deviatoric strain. The model is capable of modeling unusual behavior of strain softening during 1‐D compression (i.e., a reduction of effective stress as void ratio decreases). A good match between test results and the model simulation is demonstrated. Copyright © 2016 John Wiley & Sons, Ltd.     

SANISAND-FN: An evolving fabric-based sand model accounting for stress principal axes rotation

SANISAND is the name of a family of bounding surface plasticity constitutive models for sand within the framework of critical state theory, which have been able to realistically simulate the sand behavior under conventional monotonic and cyclic loading paths. In order to incorporate the important role of evolving fabric anisotropy, one such model was modified within the framework of the new anisotropic critical state theory and named SANISAND-F model. Yet the response under continuous stress principal axes rotation requires further modification to account for the effect of ensuing noncoaxiality on the dilatancy and plastic modulus. This modification is simpler than what is often proposed in the literature, since it does not incorporate an additional plastic loading mechanism and/or multiple dilatancy and plastic modulus expressions. The new model named SANISAND-FN is presented herein and is validated against published data for loading that includes drained stress principal axes rotation on Toyoura sand.     

堆石料的临界状态与考虑颗粒破碎的本构模型

高应力水平时堆石料的颗粒破碎对其强度和变形机制有重要影响。临界状态土力学理论对重塑土的应力-应变关系的描述较为成功,但目前颗粒破碎对堆石料的临界状态的影响及其数学描述鲜有研究。对堆石料进行了固结应力从0.4 MPa到4 MPa的18组固结排水和固结不排水常规三轴压缩试验,以及6组等向压缩试验。试验结果表明：在排水条件和不排水条件下,不同的固结应力试样都趋于临界状态;堆石料的临界状态在q-p′平面和e-lgp′平面均为非线性变化。基于此试验结果,通过引入状态参数,在广义塑性力学的理论框架下,建立了考虑颗粒破碎的堆石料本构模型,并给出了模型参数的确定方法。与长河坝料的试验进行了对比,结果表明所建议的本构模型可以较好地模拟堆石料从低围压0.4 MPa到高围压 3.5 MPa下的应力-应变特性     

基于广义塑性力学的土体次加载面循环塑性模型(Ⅰ): 理论与模型

简要地介绍了次加载面理论的基本思想、假设及其物理解释。在广义塑性力学的框架内，引入次加载面的思想，把常规的椭圆-抛物线双屈服面模型，扩展为次加载面循环塑性模型，以反映循环荷载作用下土体的曼辛效应与棘轮效应。模型能考虑塑性应变增量对应力增量的相关性，既能反映土体的循环加载特性，又能反映正常固结土和超固结土的单调加载特性。     

Selection of sand models and identification of parameters using an enhanced genetic algorithm

Numerous constitutive models of granular soils have been developed during the last few decades. As a consequence, how to select an appropriate model with the necessary features based on conventional tests and with an easy way of identifying parameters for geotechnical applications has become a major issue. This paper aims to discuss the selection of sand models and parameters identification by using genetic algorithm. A real‐coded genetic algorithm is enhanced for the optimization with high efficiency. Models with gradually varying features (elastic‐perfectly plastic modelling, nonlinear stress–strain hardening, critical state concept and two‐surface concept) are selected from numerous sand models as examples for optimization. Conventional triaxial tests on Hostun sand are selected as the objectives in the optimization. Four key points are then discussed in turn: (i) which features are necessary to be accounted for in constitutive modelling of sand; (ii) which type of tests (drained and/or undrained) should be selected for an optimal identification of parameters; (iii) what is the minimum number of tests that should be selected for parameter identification; and (iv) what is the suitable and least strain level of objective tests to obtain reliable and reasonable parameters. Finally, a useful guide, based on all comparisons, is provided at the end of the discussion. Copyright © 2015 John Wiley & Sons, Ltd.     

A small‐strain overlay model

The paper begins with a discussion of the phenomenon of small‐strain stiffness and presents the small‐strain overlay model, a simple model that takes into account the non‐linear stiffness of soils at small strains. The new model can enhance already established elastoplastic formulations for non‐linear stiffness variation at small strains in a similar way that intergranular strain enhances the hypoplastic model. The overlay model is driven by the material's strain history and only two additional material constants, both with clear physical meaning. Therefore, the proposed model is a step towards the incorporation of small‐strain stiffness into routine design. In this paper, the new small‐strain overlay model is formulated. Its ability to take into account the influence of various strain histories on soil stiffness is illustrated in several examples. Copyright © 2008 John Wiley & Sons, Ltd.     

软粘土的各向异性和小应变条件下的本构模型(ASM)

以能量方程为基础推导出各向异性和小应变条件下弹塑性的本构方程,并将屈服面与修正桥模型和试验测得的屈服面进行了对比,表明新的屈服面比修正剑桥模型更接近试验结果,采用从初始应力到状态边界面的距离为参数的系数来修正硬化准则以模拟土体小应谱条件下的应力－应变特征。将新的本构模型编入有限元程序,对小应变试验进行了计算分析,计算结果表明,新的模型比修正剑桥模型能更好地反映土体的主要力学特性。     

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.