各向异性弹黏塑性模型在ABAQUS中的研发

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明：(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。     

A viscoplastic subloading soil model for rate‐dependent cyclic anisotropic structured behaviour

This paper presents a new purely viscoplastic soil model based on the subloading surface concept with a mobile centre of homothety, enabling the occurrence of viscoplastic strains inside the yield surface and avoiding the abrupt change in stiffness of the traditional overstress viscoplastic models. This is required for overconsolidated soils. The model is formulated to reproduce the soil rate‐dependent behaviour under cyclic loading (changes in loading direction) and incorporates both initial and induced anisotropy, as well as destructuring. The model shows good qualitative response to some imposed three‐dimensional stress paths under quasi‐inviscid (elastoplastic) behaviour. Some of the main time‐dependent aspects of soil behaviour that the model is capable of reproducing were also illustrated. The capability of the model to adequately reproduce the results from an undrained triaxial test performed on stiff overconsolidated clays from the Lisbon region (Formação de Benfica), with an unloading–reloading deviatoric stress cycle at constant mean stress, that incorporates a series of staggered fast loading and creep stages, was evaluated. The model was able to reproduce well the main observed aspects of the time‐dependent stress–strain response and pore pressure evolution of a stiff overconsolidated clay under complex loading. The revised and generalised viscoplastic subloading surface concept is viable and can be applied to a consistent extension to viscoplasticity, including in the interior of the yield surface, of existing elastoplastic models formulated for soils and other materials. Copyright © 2016 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.     

An anisotropic hardening model for soils and its application to cyclic loading

An anisotropic hardening model for soils is proposed by applying the concept of a field of hardening moduli developed previously for metals. Besides the yield surface, a set of nesting surfaces in the stress-space specifies the variation of hardening moduli during the deformation process. Both drained and undrained soil behaviour can be treated and distortional as well as volumetric strain cycles can be considered. The model can be applied in studying soil behaviour under cyclic loading and in particular to describe densification or liquefaction phenomena.     

Integration of creep into a modified hardening soil model for time-dependent analysis of a high rockfill dam

This paper proposes an approach to extend a modified hardening soil model into the time-dependent analysis of high rockfill dams. The key feature of the extended model is the integration of a creep model and the double yield surface elasto-plastic model, via modifying the hardening functions. The integrated model is validated by a laboratory multistage creep test. Results from three-dimensional analyses including the without creep and with creep considerations of the Nam Ngum 2 dam have been compared with the in-situ measurements. The prediction of dam deformation can be significantly enhanced by the proposed model.     

Mathematical modelling of monotonic and cyclic undrained clay behaviour

The proposed general analytical model describes the anisotropic, elasto-plastic, path-dependent, stress-strain-strength properties of inviscid saturated clays under undrained loading conditions. The model combines properties of isotropic and kinematic plasticity by introducing the concept of a field of plastic moduli which is defined in stress space by the relative configuration of yield surfaces. For any loading (or unloading) history, the instantaneous configuration is determined by calculating the translation and contraction (or expansion) of each yield surface. The stress-strain behaviour of clays can thus be determined for complex loading paths and in particular for cyclic loadings. The stress-strain relationships are provided for use in finite element analyses. The model parameters required to characterize the behaviour of any given clay can be derived entirely from conventional triaxial or simple shear soil test results. The model's extreme versatility is demonstrated by using it to formulate the behaviour of the Drammen clay under both monotonic and cyclic loading conditions. The parameters are determined by using solely the results from monotonic and cyclic strain-controlled simple shear experimental tests, and the model's accuracy is evaluated by applying it to predict the results of other tests such as (1) cyclic stress-controlled simple shear tests, (2) monotonic triaxial loading compression and unloading extension tests, and (3) cyclic stress- and strain-controlled triaxial tests on, this same clay. The theoretical predictions are found to agree extremely well with the experimental test results.     

Fractional order model for granular soils under drained cyclic loading

The theory of fractional calculus has been successfully applied to model the triaxial behaviour of soils under static loading conditions. However, limited work has been carried out in using the fractional calculus to describe the cyclic behaviour of granular soils. In this paper, a fractional order constitutive model for granular soils under drained cyclic loading is proposed by incorporating the concept of fractional rate for strain accumulation. The fractional rate for strain accumulation is obtained from the analysis of the experimental data by utilizing the fractional calculus. Comparison between the test results and model predictions is presented. The key feature of the proposed model is that it can reasonably characterise the cyclic deformation of granular soils under both low and high loading cycles. Copyright © 2016 John Wiley & Sons, Ltd.     

基于分数阶微积分的粗粒料静动力边界面本构模型

分数阶微分理论在土体静力黏弹性本构模型中得到了广泛应用,然而,其在动力弹塑性模型中的应用尚不多见。为此,基于分数阶微积分理论分析了粗粒料在循环荷载下的变形特性,提出了粗粒料在循环荷载下的分数阶应变率;并以此为基础,进一步建立了粗粒料受静动力荷载作用下的边界面塑性力学本构模型。所提出模型包含10个参数,均可以运用常规三轴试验获得。为了验证所提出模型,选取了几种已有不同文献中的不同粗粒料试验数据进行了模拟,发现,所提出的模型可以较好地模拟粗粒料在静动力加载下的应力-应变行为,对于循环荷载下的长期变形也能较好地预测。     

A hierarchical approach for constitutive modelling of geologic materials

A hierarchical concept is proposed for the development of constitutive models to account for various factors that influence behaviour of (geologic) materials. It permits evolution of models of progressively higher grades from the basic model representing isotropic hardening with associative behaviour. Factors such as non-associativeness and induced anisotropy due to friction and cyclic loading, and softening are introduced as corrections or perturbations to the basic model. The influence of these factors is captured through non-associativeness manifested by deviation from normality of the plastic strain increments to the yield surface, F. Details of four models: isotropic hardening with associative behaviour, isotropic hardening with non-associative behavioural anisotropic hardening and strain-softening with a damage variable are presented. They are verified with respect to laboratory multiaxial test data under various paths of loading, unloading and reloading for typical soils, rock and concrete. The proposed concept is general, yet sufficiently simplified in terms of physical understanding, number of constants and their physical meanings, determination of the constants and implementation.     

FINITE ELEMENT MODELLING OF A DEEP SEA CLAY IN LONG-TERM LABORATORY CREEP TESTS

A new finite element program is introduced and its predictive capabilities are compared to results from two long-term, drained laboratory creep tests on a deep sea clay. The constitutive behaviour is based on Cam clay critical state plasticity theory with creep and time-dependent hardening. Creep is computed using either Singh–Mitchell's three-dimensional equation or Taylor's secondary compression relationship. The experimental creep data include a triaxial specimen subjected to two deviatoric stress increments and a one-dimensional consolidation specimen subjected to three vertical stress increments. In addition, the pore pressure behaviour following an increase in stress is examined in the triaxial sample. Predictions compare favourably to test data, which provide confidence for applying the chosen constitutive model and numerical formulation to solve seabed-related problems on the continental slope that are of interest to geologists, the oil industry and the navy, among others.     

Predictive potential of Perzyna viscoplastic modelling for granular geomaterials

This paper reappraises Perzyna-type viscoplasticity for the constitutive modelling of granular geomaterials, with emphasis on the simulation of rate/time effects of different magnitude. An existing elasto-plastic model for sands is first recast into a Perzyna viscoplastic formulation and then calibrated/validated against laboratory test results on Hostun sand from the literature. Notable model features include (1) enhanced definition of the viscous nucleus function and (2) void ratio dependence of stiffness and viscous parameters, to model the pycnotropic behaviour of granular materials with a single set of parameters, uniquely identified against standard creep and triaxial test results. The comparison between experimental data and numerical simulations points out the predicative capability of the developed model and the complexity of defining a unique viscous nucleus function to capture sand behaviour under different loading/initial/boundary and drainage conditions. It is concluded that the unified viscoplastic simulation of both drained and undrained response is particularly challenging within Perzyna's framework and opens to future research in the area. The discussion presented is relevant, for instance, to the simulation of multiphase strain localisation phenomena, such as those associated to slope stability problems in variably saturated soils.     

New Viscoplastic Model for Design Analysis of Tunnels in Squeezing Conditions

Summary  This paper is intended to describe the SHELVIP (Stress Hardening ELastic VIscous Plastic) model, a new viscoplastic constitutive law which has been developed to incorporate the most important features of behaviour observed in tunnels excavated in severe to very severe squeezing conditions. This model couples the elastoplastic and time-dependent behaviour by using a plastic yield surface, as frequently adopted in tunnel design analysis, and the definition of a state of overstress referred to a viscoplastic yield surface. The model is formulated in all its detailed aspects. The related analytical closed-form solution for representing triaxial creep deformations is developed. Also developed is an incremental numerical solution for describing the triaxial stress–strain behaviour under constant strain rate conditions. The model is shown to fit very satisfactorily the results of creep tests on clay shales and relaxation tests on coal specimens, as recently performed for design analysis of tunnels in squeezing conditions. Correspondence: D. Debernardi, Department of Structural and Geotechnical Engineering, Politecnico di Torino, Torino, Italy     

Anisotropic bounding-surface plasticity model for the cyclic shakedown and degradation of saturated clay

The cyclic behaviours of embedded offshore structures under different cyclic loading levels are related to the cyclic shakedown and degradation of the surrounding soils. In the present study, a damage-dependent bounding-surface model based on a newly proposed hardening rule was developed to predict the cyclic shakedown and degradation of saturated clay and the effect of the initial anisotropic stress state. By extending the Masing’s rule to the bounding-surface plasticity theory, the stress reversal point is taken as the generalised homological centre of the bounding surface. With movement of the generalised homological centre, at lower stress amplitudes, the cyclic process ends at a steady state, and cyclic shakedown is reached. At higher stress amplitudes, a damage parameter related to the accumulated deviatoric plastic strain is incorporated into the form of the bounding surface, which is hence able to contract to model degradations in stiffness and strength. To take into account the effects of initial anisotropic conditions on the cyclic behaviour of soils, an initial anisotropic tensor is introduced in the bounding surface. The developed model is validated through undrained isotropic and anisotropic cyclic triaxial tests in normally consolidated and overconsolidated saturated clay under both one-way and two-way loadings. Both cyclic shakedown and degradation are well reproduced by the model, as is the anisotropy effect induced by the initial anisotropic consolidation process.     

Constitutive model for cyclic behaviour of cohesionless sands

This paper presents a constitutive model for describing the stress-strain response of sands under cyclic loading. The model, formulated using the critical state theory within the bounding surface plasticity framework, is an upgraded version of an existing model developed for monotonic behaviour of cohesionless sands. With modification of the hardening law, plastic volumetric strain increment and unloading plastic modulus, the original model was modified to simulate cyclic loading. The proposed model was validated against triaxial cyclic loading tests for Fuji River sand, Toyoura sand and Nigata sand. Comparison between the measured and predicted results suggests that the proposed modified model can capture the main features of cohesionless sands under drained and undrained cyclic loading.     

Nonstationary flow surface theory for modeling the viscoplastic behaviors of soils

This paper presents a three-dimensional elastic viscoplastic model that can describe the time-dependent behaviors of soft clays. The constitutive model is formulated based on the nonstationary flow surface theory and incorporates new developments, including (i) an improved definition of the nonstationary flow surface that is capable of capturing the stress–strain behaviors under different loading paths, (ii) a unique stress–strain—viscoplastic-strain-rate equation that is able to explicitly describe the nonstationary flow surface, and (iii) a final stable state concept that identifies the final equilibrium state at the end of creep and stress relaxation, which is also used to simplify the loading criteria. The consistency condition is validated for the proposed model, and the viscoplastic multipliers are calculated by solving the consistency equations. The model performance is investigated and validated via simulation of both oedometer and triaxial tests. The numerical results demonstrate that the proposed model is able to reproduce the main viscoplastic behaviors of soils, including creep, undrained creep rupture, stress relaxation, rate effect and accumulated effect.     

Non‐isothermal plasticity model for cyclic behaviour of soils

On the one hand, it has been observed that liquefaction‐induced shear deformation of soils accumulates in a cycle‐by‐cycle pattern. On the other hand, it is known that heating could induce plastic hardening. This study deals with the constitutive modelling of the effect that heat may have on the cyclic mechanical properties of cohesive soils, a relatively new area of interest in soil mechanics. In this paper, after a presentation of the thermo‐mechanical framework, a non‐isothermal plasticity cyclic model formulation is presented and discussed. The model calibration is described based on data from laboratory sample tests. It includes numerical simulations of triaxial shear tests at various constant temperatures. Then, the model predictions are compared with experimental results and discussed in the final section. Both drained and undrained loading conditions are considered. The proposed constitutive model shows good ability to capture the characteristic features of behaviour. Copyright © 2007 John Wiley & Sons, Ltd.     

机车荷载作用下青藏铁路多年冻土区普通路基的蠕变分析

基于青藏铁路北麓河试验段的机车现场实时振动监测,在三轴流变试验的基础上,利用ABAQUS有限元软件中的时间硬化与Druker－Prager屈服破坏准则耦合的蠕变模型,对青藏铁路多年冻土区的普通素土路基在等效机车荷载作用下的蠕变效应进行分析。研究表明,ABAQUS软件中时间硬化与蠕变模型能较好地模拟路基运营过程中的时间效应,素土路基的竖向位移由路基顶部中心位置向路基内部和两侧逐渐减小,一年中的最大位移达14 mm,与该段的年沉降变形观测值一致。路基中心的蠕变应变值大于路基顶面和天然地表的应变值,路基中心受机车荷载的振动影响最大     

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

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