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.     

Laboratory modelling of shear behaviour of soft joints under constant normal stiffness conditions

Shear behaviour of regular sawtooth rock joints produced from casting plaster are investigated under constant normal stiffness (CNS) conditions. Test results obtained in this investigation are also compared with the constant normal load (CNL) tests. It is observed that the peak shear stress obtained under CNL conditions always underestimates the peak shear stress corresponding to the CNS condition. Plots of shear stress against normal stress show that a nonlinear (curved) strength envelope is acceptable for soft rock joints subjected to a CNS condition, in comparison with the linear or bilinear envelopes often proposed for a CNL condition. Models proposed by Patton (1966) and Barton (1973) have also been considered for the predictions of peak shear stress of soft joints under CNS conditions. Although Patton's model is appropriate for low asperity angles, it overestimates the shear strength in the low to medium normal stress range at higher asperity angles. In contrast, while Barton's model is realistic for the CNL condition, it seems to be inappropriate for modelling the shear behaviour of soft joints under CNS conditions. The effect of infill material on the shear behaviour of the model joints is also investigated, and it is found that a small thickness of bentonite infill reduces the peak stress significantly. The peak shear stress almost approached that of the shear strength of infill when the infill thickness to asperity height ratio (t/a) reached 1.40. This paper also introduces an original, empirical shear strength envelope to account for the change in normal stress and surface degradation during CNS shearing. © Rapid Science Ltd. 1998     

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.     

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

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.     

Soil–structure cyclic direct shear tests: a new interpretation of the direct shear experiment and its application to a series of cyclic tests

This paper presents a series of cyclic 2D direct shear tests on sand–rough material interfaces under constant normal load (CNL) and constant normal stiffness (CNS) conditions. The aim of these tests is to describe the behavior of the soil–pile contact subjected to a large number of cycles due to environmental or anthropic loadings. These cycles (typically 10⁴ or less due to an early rupture) are small (10, 20 and 40 kPa in terms of shear stress). A new interpretation of the direct shear tests is proposed. The sample of soil is schematically composed of a sheared interface and of a buffer under oedometric load. The problem of sand leakage between the shear box and the rough plate, classical phenomenon in this type of test, is focused. The effect of initial density, position of “center of cycles” in stress plane (mean cyclic variables) and cyclic amplitude is investigated. The cycles are defined by the initial mean cyclic normal stress, the level of initial mean cyclic stress ratio and the normalized cyclic amplitude. Under CNL condition, either dilation or contraction is exhibited, in agreement with the characteristic state developed by Luong (International symposium on soils under cyclic and transient loading, Swansea, 7–11 January, pp 315–324, 1980). The influence of a prescribed normal stiffness is especially considered. It can be highlighted that CNS cyclic paths are always contractive. This contraction results in a drop of mean cyclic normal stress often called degradation of friction.     

An Elasto-Plastic Constitutive Model for Rock Joints Under Cyclic Loading and Constant Normal Stiffness Conditions

An elasto-plastic constitutive model is introduced for rock joints under cyclic loading, considering the additional shear resistance generated by the asperity damage in the first forward shear cycle and sliding mechanism for further shearing. A series of cyclic loading direct shear tests was conducted on artificial joints with triangular asperities and replicas of a real rock asperity surface under constant normal stiffness (CNS) conditions. The model was calibrated and then validated using selected data sets from the experimental results. Model simulations were found to be in good agreement with the rock joints behaviour under cyclic loading and CNS conditions both in stress prediction and dilation behaviour. In addition, dynamic stability analysis of an underground structure was carried out, using Universal Distinct Element Code and the proposed constitutive model.     

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.     

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

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

Study of rock joints under cyclic loading conditions

Summary A conceptual model for the behaviour of rock joints during cyclic shear and under constant normal stresses was proposed according to results from shear tests with 50 concrete replicas of rock joints. The shear strength and deformability of joint samples were found to be both anisotropic and stress dependent. Based on these experimental results, a two-dimensional constitutive model was developed for rock joints undergoing monotonic or cyclic loading sequences. The joint model was formulated in the framework of non-associated plasticity, coupled with empirical relations representing the surface roughness degradation, appearance of peak and residual shear stresses, different rates of dilatancy and contraction, variable normal stiffness with normal deformation, and dependence of shear strength and deformability on the normal stress. The second law of thermodynamics was represented by an inequality and used to restrict the values of some of the material parameters in the joint model. The new joint model was implemented into a two-dimensional Distinct Element Method Code, UDEC, and its predictions agreed well with some well-known test results.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part I—tests with monotonic loading and stress cycles

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such a model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. The purpose of these two companion papers was to provide such database collected for a fine sand. The database consists of numerous undrained cyclic triaxial tests with stress or strain cycles applied to samples consolidated isotropically or anisotropically. Monotonic triaxial tests with drained or undrained conditions have also been performed. Furthermore, drained triaxial, oedometric or isotropic compression tests with several un- and reloading cycles are presented. Part I concentrates on the triaxial tests with monotonic loading or stress cycles. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

An experimental database for the development,calibration and verification of constitutive models for sand with focus to cyclic loading: part II—tests with strain cycles and combined loading

For numerical studies of geotechnical structures under earthquake loading, aiming to examine a possible failure due to liquefaction, using a sophisticated constitutive model for the soil is indispensable. Such model must adequately describe the material response to a cyclic loading under constant volume (undrained) conditions, amongst others the relaxation of effective stress (pore pressure accumulation) or the effective stress loops repeatedly passed through after a sufficiently large number of cycles (cyclic mobility, stress attractors). The soil behaviour under undrained cyclic loading is manifold, depending on the initial conditions (e.g. density, fabric, effective mean pressure, stress ratio) and the load characteristics (e.g. amplitude of the cycles, application of stress or strain cycles). In order to develop, calibrate and verify a constitutive model with focus to undrained cyclic loading, the data from high-quality laboratory tests comprising a variety of initial conditions and load characteristics are necessary. It is the purpose of these two companion papers to provide such database collected for a fine sand. Part II concentrates on the undrained triaxial tests with strain cycles, where a large range of strain amplitudes has been studied. Furthermore, oedometric and isotropic compression tests as well as drained triaxial tests with un- and reloading cycles are discussed. A combined monotonic and cyclic loading has been also studied in undrained triaxial tests. All test data presented herein will be available from the homepage of the first author. As an example of the examination of an existing constitutive model, the experimental data are compared to element test simulations using hypoplasticity with intergranular strain.     

Hyperelastic modelling of small‐strain stiffness anisotropy of cyclically loaded sand

Experimental evidence shows that soil stiffness at very small strains is strongly anisotropic and depends on the stress level and void ratio. In particular, stiffness anisotropy varies considerably in sand when subjected to cyclic loading, following the stress cycles applied. To model this behaviour, an innovative hyperelastic formulation based on the elastoplastic coupling is incorporated in a new kinematic hardening elastoplastic model. The proposed hyperelastic–plastic model is the first to be capable of correctly simulating all aspects of the small‐strain behaviour of granular materials subjected to monotonic and cyclic loads. This hyperelastic formulation is generally applicable to any elastoplastic model. Copyright © 2009 John Wiley & Sons, Ltd.     

循环荷载作用下饱和软黏土的损伤模型

基于各向同性弹塑性损伤和Prevost模型的基本理论,把弹塑性等向硬化、运动硬化和各向同性损伤结合起来,推导了循环荷载作用下不排水饱和软黏土的弹塑性动力损伤本构模型。由循环累积塑性偏应变建立损伤演化方程,以描述循环加载对软黏土结构的破坏作用,并通过对循环三轴试验的模拟验证了模型的有效性。结果表明,提出的模型能较好地描述饱和软黏土在循环荷载作用下的变形、孔压变化及模量损伤过程。     

A Laboratory Shear Cell Used for Simulation of Shear Strength and Asperity Degradation of Rough Rock Fractures

Different failure modes during fracture shearing have been introduced including dilation, sliding, asperity cut-off and degradation. Several laboratory studies have reported the complexity of these failure modes during shear tests performed under either constant normal load (CNL) or constant normal stiffness (CNS) conditions. This paper is concerned with the mechanical behaviour of synthetic fractures during direct shear tests using a modified shear cell and related numerical simulation studies. The modifications made to an existing true triaxial stress cell (TTSC) in order to use it for performing shear tests under CNL conditions are presented. The large loading capacity and the use of accurate hydraulic pumps capable of applying a constant shear velocity are the main elements of this cell. Synthetic mortar specimens with different fracture surface geometries are tested to study the failure modes, including fracture sliding, asperity degradation, and to understand failure during shearing. A bonded particle model of the direct shear test with the PFC2D particle flow code is used to mimic the tests performed. The results of a number of tests are presented and compared with PFC2D simulations. The satisfactory results obtained both qualitatively and quantitatively are discussed.     

Constitutive description of interface behavior including cyclic loading and particle breakage within the framework of critical state soil mechanics

The cyclic behavior of soil–structure interface can be very important in dynamic problems. The cyclic behavior of soil–structure interface may be nonlinear, which includes hysteresis, hardening, degradation, and particle breakage. The breakage of granular soil particles during shearing of granular soil–structure interface is associated with cyclic degradation and can be critical to the dynamic behavior of soil–structure system. The critical state soil mechanics concept formerly used to simulate the nonlinear monotonic behavior of granular soil–structure interface was modified and extended to describe the cyclic behavior, especially soil‐particle breakage and degradation during cyclic shearing. Soil‐particle breakage was assumed to relate to the energy consumption during cyclic shearing and the critical state line of the soil–structure interface was assumed to translate with the consumption of shearing energy during cyclic shearing as the threshold value is attained. The model was formulated in the framework of generalized plasticity and is capable of describing the salient features of granular soil–structure interface under cyclic loading. Most of the model parameters have straightforward physical meanings and are calibrated using monotonic or cyclic interface test results. The proposed model was calibrated and validated against published test results. The dependency of interface behavior on stress path and cyclic degradation can be successfully described by the proposed model. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

A Constitutive Model for Shear Behavior of Rock Joints Based on Three-Dimensional Quantification of Joint Roughness

A new constitutive model to describe the shear behavior of rock joints under constant normal stiffness (CNS) and constant normal load (CNL) conditions is proposed. The model was developed using an empirical approach based on the results of a total of 362 direct shear tests on tensile fractured rock joints and replicas of tensile joints and on a new quantitative roughness parameter. This parameter, the active roughness coefficient C _r, is derived from the features of the effective roughness mobilized at the contact areas during shearing. The model involves a shear strength criterion and the relations between stresses and displacements in the normal and shear directions, where the effects of the boundary conditions and joint properties are considered by the shape indices C _d and C _f. The model can be used to predict the shear behavior under CNS as well as CNL conditions. The shear behavior obtained from the experimental results is generally in good agreement with that estimated by the proposed model, and the effects of joint roughness, initial normal stress, and normal stiffness are reasonably reflected in the model.     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.