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.     

A cyclic elastoplastic-viscoplastic constitutive model for soils

An elastoplastic-viscoplastic constitutive model for soils is presented in this study, based on an original approach concerning viscous modelling. In this approach, the viscous behaviour is defined by internal viscous variables and a viscous yield surface. The model has been developed from a basic elastoplastic model (CJS model) by considering an additional viscous mechanism. The evolution of the viscous yield surface is governed by a particular hardening called ‘viscous hardening’. This model is able to explain the time-dependent behaviour of soils such as creep (primary, secondary and un-drained creep rupture), stress relaxation and strain rate effects in static and cyclic loadings. The existing problems in the classical elasto-viscoplastic models related to the plasticity failure, the rapid loading and the cyclic loading are solved in the proposed model. The physical meanings and the identification strategy of model parameters are clearly given. The validation on certain triaxial test results and the simulation of cyclic triaxial test indicate the capacity of this model in prediction of time-dependent behaviour of clayey soils.     

A constitutive model for crushed salt

A constitutive model for crushed salt is presented in this paper. A creep constitutive model is developed first and compared with test results. The constitutive model presented here concentrates on creep deformation because saline media behave basically in a ductile and time‐dependent way. An idealized geometry is used as a common framework to obtain stress–strain macroscopic laws based on two deformation mechanisms: fluid‐assisted diffusional transfer creep and dislocation creep. The model is able to predict strain rates that compare well with results from laboratory tests under isotropic and oedometric conditions. Macroscopic laws are written using a non‐linear viscous approach, which incorporates also a viscoplastic component, based on critical state theory. The viscoplastic term is intended for non‐creep deformation mechanisms such as grain reorganization and crushing. Copyright © 2002 John Wiley & Sons, Ltd.     

高岭土一维增湿的湿陷-流变分析方法研究

塑性单元体与黏性单元体并联之后再与弹性单元体串联即构成三要素模型,其各单元体都具有非线性本构关系。通过考虑增湿过程中饱和度对三要素模型中的黏性单元体和塑性单元体的影响,建立了弹-塑-黏性与湿陷性耦合的湿陷流变本构模型。湿陷对塑性单元体的影响表现为应变不变条件下屈服应力的降低;而对黏性单元体的影响表现为速率敏感系数的升高。两种影响可以通过试验确定。以此模型为基础,分析了高岭土的单轴压缩固结试验中的湿陷-流变变形,提出了一种湿陷-流变分析方法。试验中,在干燥条件下观测了蠕变;在增湿过程中观测了湿陷。试验中对含水量的测量精度高。理论所得含水率-应变-时间关系与实测结果接近,并区分了增湿过程中的湿陷和干燥条件下的蠕变,证明了本方法的有效性。     

Delayed plastic model for time‐dependent behaviour of materials

A delayed plastic model, based on the theory of plasticity, is proposed to represent the time‐dependent behaviour of materials. It is assumed in this model that the stress can lie outside the yield surface and the conjugate stress called static stress is defined on the yield surface. The stress–strain relation is calculated based on the plastic theory embedding the static stress. Thus, the stress–strain relation of the model practically corresponds to that of the inviscid elastoplastic model under fairly low rate deformation. The delayed plastic model is coupled with the Cam‐clay model for normally consolidated clays. The performance of the model is then examined by comparing the model predictions with reported time‐dependent behaviour of clays under undrained triaxial conditions. It is shown that the model is capable of predicting the effect of strain rate during undrained shear and the undrained creep behaviour including creep rupture. Copyright © 2001 John Wiley & Sons, Ltd.     

A visco‐elasto‐plastic model for granular materials under simple shear conditions

The numerical simulation of rapid landslides is quite complex mainly because constitutive models capable of simulating the mechanical behaviour of granular materials in the pre‐collapse and post‐collapse regimes are still missing. The goal of this paper is to introduce a constitutive model capable of capturing the response of dry granular flows from quasi‐static to dynamic conditions, in particular when the material experiences a sort of solid‐to‐fluid phase transition. An ideal assembly of identical spheres under simple shear conditions is considered. In the constitutive model, void ratio and granular temperature have been chosen as state variables, and both shear and normal stresses are computed as the sum of two contributions: the quasi‐static one and the collisional one. The former is determined by using a perfect elasto‐plastic model including the critical state concept, while the latter is derived from the kinetic theory of granular gases. The evolution of the granular temperature, fundamentally governing the material phase transition, is obtained by imposing the kinetic fluctuating energy balance. The constitutive relationship has been integrated, under both constant pressure and constant volume conditions, and the influence of shear strain rate, initial void ratio and normal pressure on the mechanical response has been investigated. Copyright © 2015 John Wiley & Sons, Ltd.     

A constitutive law for rock salt based on creep and relaxation tests

Summary The paper deals with uniaxial relaxation tests on rock salt which are the basis for a constitutive equation. Since so far no regard was paid to relaxation behaviour, corresponding test results are not available and hence a special device for performing of uniaxial relaxation tests had to be constructed. Some interesting test results are discussed in the paper. These results are useful for the verification of constitutive equations. A critical analysis of conventional constitutive laws and their fundamentals shows that these material laws have considerable imperfections and are rested on substantial restrictions. Particularly, they are not able to describe relaxation behaviour sufficiently. Consequently, a constitutive law is proposed consisting of a strain hardening approach with separate creep and relaxation functions. By post-calculation of different laboratory tests it could be shown that in comparison to conventional steady-state creep equations this material law describes the viscous behaviour of rock salt more realistically.     

Modeling the viscous behavior of frozen soil with hypoplasticity

This paper presents a hypoplastic constitutive model for the viscous behavior of frozen soil. The model is composed of a ‘solid’ part and a ‘fluid’ part. The solid part is based on the extended hypoplastic model, and the fluid part is dependent on the second time derivative of strain. The performance of the model is demonstrated by simulating some uniaxial compression tests at different strain rates. Moreover, the model can describe in a unified way the three stages of typical creep tests, that is, primary, secondary, and tertiary stage. Copyright © 2016 John Wiley & Sons, Ltd.     

基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

A viscoplastic model including anisotropic damage for the time dependent behaviour of rock

This paper presents a new constitutive model for the time dependent mechanical behaviour of rock which takes into account both viscoplastic behaviour and evolution of damage with respect to time. This model is built by associating a viscoplastic constitutive law to the damage theory. The main characteristics of this model are the account of a viscoplastic volumetric strain (i.e. contractancy and dilatancy) as well as the anisotropy of damage. The latter is described by a second rank tensor. Using this model, it is possible to predict delayed rupture by determining time to failure, in creep tests for example. The identification of the model parameters is based on experiments such as creep tests, relaxation tests and quasi‐static tests. The physical meaning of these parameters is discussed and comparisons with lab tests are presented. The ability of the model to reproduce the delayed failure observed in tertiary creep is demonstrated as well as the sensitivity of the mechanical response to the rate of loading. The model could be used to simulate the evolution of the excavated damage zone around underground openings. Copyright © 2005 John Wiley & Sons, Ltd.     

Solid–fluid transition modelling in geomaterials and application to a mudflow interacting with an obstacle

Given the contrasting behaviour observed for geomaterials, for example, during landslides of the flow type, this contribution proposes an original constitutive model, which associates both an elasto‐plastic relation and a Bingham viscous law linked by a mechanical transition criterion. This last is defined as the second‐order work sign for each material point, which is a general criterion for divergence instabilities. Finite element method with Lagrangian integration points is chosen as a framework for implementing the new model because of its well‐known ability to deal with both solid and fluid behaviours in large deformation processes. A first boundary model considering a sample of initially stable soil, a slope and an obstacle is performed. The results show the power of the constitutive model because the consistent evolution of initiation, propagation and arrest of the mudflow is described. A parametric study is led on various plastic and viscous parameters to determine their influence on the flow development and arrest. Finally, forces against the obstacle are compared with good agreement with those of other authors for the same geometry and a pure viscous behaviour. Copyright © 2014 John Wiley & Sons, Ltd.     

Time‐dependent modelling for soils and its application in tunnelling

Monitoring of the progressive convergence of a tunnel shows that deformations occurring in the soil surrounding a tunnel exhibit a strong evolution with time. This time‐dependent behaviour can be linked to three essential factors: the distance from the point of interest to the working face over time, the distance of unsupported tunnel to the working face and the viscous properties of the soil. The objective of this paper is to propose a constitutive model of the time‐dependent behaviour of soil which has been developed within the framework of elastoplasticity–viscoplasticity and critical state soil mechanics. The consideration of viscoplastic characteristic sets the current model apart from the CJS (Cambou, Jafari and Sidoroff) model as the basic elastoplastic model, and introduces an additional viscous mechanism. The evolution of the viscous yield surface is governed by a particular hardening called ‘viscous hardening’ with a bounding surface. The proposed constitutive model has been applied in the analysis of tunnelling. Two kinds of numerical calculations have been used in the analysis, axisymmetric analysis and plane strain analysis. Monitoring of the progressive convergence of a tunnel conducted in the railway tunnel of Tartaiguille (France), has been used to describe the calculation procedure proposed and the capability of the model. The finite difference software, fast Lagrangian analysis of continua (FLAC), has been used for the numerical simulation of the problems. The comparison of results shows that the observed deformations could have been reasonably predicted by using the constitutive model and calculation strategy proposed. Copyright © 2004 John Wiley & Sons, Ltd.     

On incremental non‐linearity in granular media: phenomenological and multi‐scale views (Part I)

On the basis of fundamental constitutive laws such as elasticity, perfect plasticity, and pure viscosity, many elasto‐viscoplastic constitutive relations have been developed since the 1970s through phenomenological approaches. In addition, a few more recent micro‐mechanical models based on multi‐scale approaches are now able to describe the main macroscopic features of the mechanical behaviour of granular media. The purpose of this paper is to compare a phenomenological constitutive relation and a micro‐mechanical model with respect to a basic issue regularly raised about granular assemblies: the incrementally non‐linear character of their behaviour. It is shown that both phenomenological and micro‐mechanical models exhibit an incremental non‐linearity. In addition, the multi‐scale approach reveals that the macroscopic incremental non‐linearity could stem from the change in the regime of local contacts between particles (from plastic regime to elastic regime) in terms of the incremental macroscopic loading direction. Copyright © 2005 John Wiley & Sons, Ltd.     

2D viscoplastic finite element modelling of slow landslides: the Portalet case study (Spain)

This paper proposes a hydro-geomechanical finite element model to reproduce the kinematic behaviour of large slow landslides. The interaction between solid skeleton and pore fluids is modelled with a time dependent u–p _w formulation and a groundwater model that takes into account recorded daily rainfall intensity. A viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behaviour and the delayed creep deformation. The proposed model is applied to Portalet landslide (Central Spanish Pyrenees). This is an active paleo-landslide that has been reactivated by the construction of a parking area at the toe of the slope. The stability analysis reveals that, after the constructive solutions were undertaken, the slope is in a limit equilibrium situation. Nevertheless, time-dependent analysis reproduces the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. Overall, the model reproduces a 2-m displacement in the past 8  years that coincides with in situ monitoring data. The proposed model is useful for short- and mid-term predictions of secondary creep. However, long-time predictions remain uncertain, stability depends strongly on the position of the water table depth and new failures during tertiary creep due to soil temporal microstructural degradation are difficult to calibrate.     

盐岩的分数阶非线性蠕变本构模型

在非牛顿流体黏滞阻尼元件的基础上提出了分数阶非线性黏壶元件。基于分数阶导数的定义,得到了该元件的理论表达式,其表达式与Abel黏壶的表达式相一致,进一步比较发现该元件能反映蠕变曲线的非线性加速特征。通过引入该元件建立了新的非线性蠕变本构模型。基于蠕变试验数据,对该模型进行了拟合分析。结果表明,非线性蠕变本构模型与试验数据吻合的较好,并且能很好地反映全过程蠕变曲线的特征,尤其是非线性加速蠕变阶段。     

一种基于修正西原模型的冻土蠕变本构关系

冻土在高应力水平下通常会表现出加速蠕变特性。经典西原（Nishihara）模型可以反映衰减蠕变和稳态蠕变,但不适用于描述加速蠕变。通过对Nishihara模型中黏性元件进行修正,将黏弹性元件中的定常的黏滞系数修正为时间的函数,并在黏塑性元件中采用非定常、非牛顿黏性元件,得到的蠕变柔量可以考虑应力水平和时间的影响,使其能够反映其加速蠕变的过程。用提出的模型分别对一维和三维应力、不同应力水平下的冻土蠕变试验数据进行了拟合,拟合结果能够反映冻土从衰减蠕变和稳定蠕变状态进入加速蠕变的过程,说明了该模型对于描述冻土在不同应力下不同蠕变状态的适用性。     

基于断裂及高温损伤的岩石蠕变模型研究

为了反映热-力耦合作用下岩石蠕变变形的全过程,依据断裂力学原理提出了岩石裂纹扩展的临界损伤应力和一个新的可描述岩石在稳态蠕变阶段与临界损伤应力相关的非线性黏性分量,在传统西原模型和Burgers模型的基础上,将指数形式的损伤变量、临界损伤应力以及与其有关的非线性黏性分量引入到流变微分方程,通过叠加原理推导了考虑温度效应的单轴和三轴压缩条件下岩石的流变本构关系,建立了岩石的热-力耦合损伤蠕变本构模型。利用不同温度、不同应力条件下花岗岩的三轴蠕变试验曲线和本文蠕变模型的计算曲线进行比较,结果表明本文蠕变模型能够较好地模拟岩石在初始瞬态、稳态和加速蠕变阶段全过程的变形规律,验证了所建模型的有效性和合理性。该模型为分析高温、高应力环境下岩石工程的长期变形和稳定情况提供了理论依据。     

On the time-dependent effects in advancing tunnels

A finite element study is presented of the short-term effects that develop when a tunnel is driven in a ground showing viscous behaviour associated with the deviatoric deformations. Axisymmetric conditions around the tunnel centreline are assumed and the process of excavation is simulated by means of a step-by-step time incremental technique. The Kelvin's model is used to approximate the medium creep behaviour and a simple procedure is presented that allows the determination of the creep model constants from the laboratory test data. Various rates of advancing, the presence of an unsupported zone close to the tunnel face and the temporary interruption of the excavation process are considered. For the viscous model and the analytical technique used, simple, approximate relationships are presented for the displacements around the tunnel and the overexcavation volume as functions of the rate of excavation advancing.