On double shearing in frictional materials

This paper evaluates the mechanical behaviour of yielding frictional geomaterials. The general Double Shearing model describes this behaviour. Non‐coaxiality of stress and plastic strain increments for plane strain conditions forms an important part of this model. The model is based on a micro‐mechanical and macro‐mechanical formulation. The stress–dilatancy theory in the model combines the mechanical behaviour on both scales. It is shown that the general Double Shearing formulation comprises other Double Shearing models. These models differ in the relation between the mobilized friction and dilatancy and in non‐coaxiality. In order to describe reversible and irreversible deformations the general Double Shearing model is extended with elasticity. The failure of soil masses is controlled by shear mechanisms. These shear mechanisms are determined by the conditions along the shear band. The shear stress ratio of a shear band depends on the orientation of the stress in the shear band. There is a difference between the peak strength and the residual strength in the shear band. While peak stress depends on strength properties only, the residual strength depends upon the yield conditions and the plastic deformation mechanisms and is generally considerably lower than the maximum strength. It is shown that non‐coaxial models give non‐unique solutions for the shear stress ratio on the shear band. The Double Shearing model is applied to various failure problems of soils such as the direct simple shear test, the biaxial test, infinite slopes, interfaces and for the calculation of the undrained shear strength. Copyright © 2006 John Wiley & Sons, Ltd.     

基于广义位势理论的土的非共轴特性研究

传统的塑性位势理论隐含了应力主方向和塑性应变增量主方向共轴的假定,无法客观地描述主应力轴旋转过程中的非共轴现象。基于广义位势理论提出的拟弹性弹塑性本构模型,把总的塑性应变分解为满足弹性分解准则的拟弹性部分和符合传统塑性理论假设的纯塑性部分,分解后建立的模型更为合理和简便,同时又可以解决土的非共轴问题。通过单剪试验结果的验证表明,基于广义位势理论的拟弹性弹塑性模型的模拟效果较好,传统的弹塑性模型（共轴模型）模拟得到的主应力方向和塑性主应变增量方向保持共轴,而拟弹性弹塑性模型（非共轴模型）的模拟结果则能够合理地描述主应力轴旋转过程中的非共轴特性,结果更符合实际,从而为解决土的非共轴特性问题提供了一种有效的方法。     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulations of simple shear with non‐coaxial soil models

This paper investigates the effects of a non‐coaxial model on simulated stress–strain behaviour of granular materials subject to simple shearing under various initial conditions. In most cases, a significant difference of predictions between coaxial and non‐coaxial modelling is found during the early stage in shearing. With the increase in shearing, non‐coaxial simulations approach and tend to coincide with coaxial simulations. It is also found that the roles of non‐coaxial modelling in simulating simple shear behaviour are considerably influenced by hardening rules, flow rules, initial static lateral pressure coefficients. In some cases, the non‐coaxial modelling gives a similar simulation as the coaxial modelling. In other cases, the non‐coaxial modelling decreases the hardening response or softening response of materials, compared with the coaxial modelling. Under certain conditions, the predicted peak strength of materials with non‐coaxial modelling is larger than that for coaxial modelling. Some of these observations can be attributed to the amount of principal stress rotation in various cases analysed. Others can be attributed to the difference between the directions of the non‐coaxial plastic flow and those for coaxial plastic flow. Copyright © 2005 John Wiley & Sons, Ltd.     

The role of non‐coaxiality in the simulation of strain localization based on classical and Cosserat continua

It is normally accepted that materials inside the shear band undergo severe rotation of the principal stress direction, which causes non‐coaxiality between the principal stress and principal plastic strain rate. However, classical plasticity flow theory implicitly assumes that the principal stress and the principal plastic strain rate are coaxial; thus, it may not correctly predict the onset of the shear band. In addition, classical continuum does not contain any internal length scales; as a result, it cannot provide a reasonable shear band thickness. In this study, the original vertex non‐coaxial plastic model based on the classical continuum is extended to the Cosserat continuum. The corresponding codes are implemented via the interface of the user defined element subroutine in ABAQUS. Through a simple shear test, the effectiveness of the user's codes is verified. Through a uniaxial compression test, the influence of non‐coaxiality on the onset, the orientation, and the thickness of the shear band is investigated. Results show that the onset of the shear localization is delayed, and the thickness of the shear band is widened when the non‐coaxial degree increases, while the orientation of the shear band is little affected by the non‐coaxial degree. In addition, it is found that the non‐coaxiality can weaken the micro‐polar effect to some extent; nonetheless, the Cosserat non‐coaxial model still has its advantage over the classical non‐coaxial model in capturing the pre‐bifurcation as well as the post‐bifurcation behaviors of strain localization. Copyright © 2016 John Wiley & Sons, Ltd.     

Comment on the postulate of three plane strain mechanisms

This note discusses the inconsistencies that are inherent in the postulate of three plane strain mechanisms. It is shown that this postulate violates the principle of invariance and one obtains different results depending on the choice of the reference axes. If formulated in the principal stress space, this postulate requires that the principal stress and principal plastic strain increment directions be coaxial. Constitutive models based on this postulate cannot be used for general loading situations involving principal stress rotation where significant non-coaxiality is obtained.     

A non‐coaxial constitutive model for sand deformation under rotation of principal stress axes

A constitutive model for the simulation of non‐coaxiality, an aspect of anisotropic behavior of sand subjected to the rotation of the principal stress axes, is presented in this paper. Experimental studies have shown that non‐coaxiality or non‐coincidence of principal plastic strain increments with principal stress axes under loadings involving the rotation of principal stress axes may be considerable. Besides, the rotation of the principal stress axes results in dramatic effects on stiffness and dilatant behavior of sand. Therefore, the consequences of principal stress axes rotation on deformational behavior, dilatancy and soil stiffness must be taken into account in theoretical and practical problems. To this aim, the following steps are taken: (1) A general relationship for flow direction with respect to possibility of non‐coaxial flow is developed. Moreover, special circumstances linking non‐coaxiality to instantaneous interaction between loading and soil fabric are proposed. (2) Proposing novel expressions for plastic modulus and dilatancy function, the model is enforced to provide realistic simulations when sand is subjected to the rotation of the principal stress axes. Finally, with numerous examples and comparisons, the model capabilities are shown under various stress paths and drainage conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

Simple shear of isotropic elasto–plastic soil

The implications of assuming isotropic elasto–plasticity to model the behaviour of soil under simple shear conditions are considered. For small strains, use of such a model implies the following three consequences: (1) strains and strain increments at any stage of shearing may be expressed as the sum of elastic and plastic components; (2) principal directions of stress and of plastic strain increment are collinear; (3) principal directions of stress increment and of elastic strain increment are collinear. These consequences are used in order to establish relationships between the stresses, stress increments and strains which develop in a simple shear test. No additional assumptions with regards the form of the yield function, the flow rule or the hardening function are required for this development. By defining the ratio of the plastic to the total shear strain increment on the horizontal plane (the plane of zero extension) as λ, it is possible to define the horizontal normal stress σ_x in terms of λ and other stresses and strains which are normally known during simple shear loading. As a result, all components of the stress tensor in the simple shear plane may be defined. Results of some direct simple shear tests on soft clay have been interpreted using the model and found to be generally consistent with some of the observations reported in the literature from tests in which boundary stresses were measured.     

非共轴本构模型在地基承载力数值计算中若干影响因素的探讨

土体在剪切变形过程中产生主应力方向的旋转时,主应变增量方向与主应力方向之间存在着非共轴现象,然而传统的弹塑性本构模型未能考虑该现象的影响。通过在屈服面的切线方向增加一项非共轴塑性应变增量,即可实现对非共轴现象的反映。采用显式积分算法和自动分步方法,将非共轴本构模型运用到桶形基础地基承载力问题的有限元计算中,并讨论了流动法则、内摩擦角、膨胀角等因素与非共轴模型的联系。计算结果表明：采用有限元程序默认容许误差时,该本构模型可达到理想的收敛精度,并且,该模型对关联、非关联流动法则均适用。采用共轴模型进行数值计算时,不同流动法则对计算结果的影响可以忽略;采用非共轴模型时,不同流动法则的计算结果之间存在差异。非共轴现象对地基承载力-位移曲线具有软化作用,并且,该软化作用在采用非关联流动法则时变得更加明显     

A non‐coaxial and unbalanced plastic flow rule for geomaterials and its application to the shear band orientation prediction

In this paper, the non‐coaxial relation between the principal plastic strain increments and the principal stresses, which results from the internal friction in geomaterials, is analyzed, and the phenomenon of the unbalanced development of plastic flow in two conjugate directions is discussed. A non‐coaxial, unbalanced plastic flow model for Coulomb frictional materials is developed and used to determine the orientation of shear band in geomaterials. It is shown that the unbalanced index r of plastic flow has important effect on the orientation of the shear band, and the orientation determined by the conventional plastic flow theory is only a special case of the proposed model when r=0. This result soundly explains the reason that the geomaterials with the same internal friction angle and dilatancy angle can have very different shear band orientations. In addition, the difference between the intrinsic and apparent dilatancy angles is analyzed, and it is emphasized that the dilatancy angle commonly used in practice is indeed the apparent dilatancy angle. Copyright © 2010 John Wiley & Sons, Ltd.     

Response of a circular opening in a friable low-permeability medium to temperature and pore pressure changes

A general analysis using an incremental elastic, perfectly plastic constitutive stress–strain relationship for poroelastoplastic materials is presented to simulate an opening in a low-permeability friable porous medium under non-isothermal conditions. Analytical solutions are obtained for the stresses and strains around a 2-D plane strain circular borehole. An expansion potential is introduced by combining the strains induced by temperature and pore pressure changes. Steady-state pressures and temperatures are considered, and a non-associated plastic flow rule is applied to calculate plastic strains. Focusing on stress distribution near a circular opening, the classic solutions for those stresses under dry and isothermal conditions are used to compare with the newly derived solution. The general poroelastoplastic effect and the thermal effect on sand production and borehole stability are addressed. We suggest that the knowledge of stress history is critical to achieve adequate solutions for displacement and stress in friable media such as clays, shales and oil sands.     

A kinematic hardening constitutive model for sands: the multiaxial formulation

This paper explores the possibility of using well-accepted concepts—Mohr-Coulomb-like strength criterion, critical state, existence of a small strain elastic region, hyperbolic relationship for representing global plastic stress–strain behaviour, dependence of strength on state parameter and flow rules derived from the Cam-Clay Model—to represent the general multiaxial stress–strain behaviour of granular materials over the full range of void ratios and stress level (neglecting grain crushing). The result is a simple model based on bounding surface and kinematic hardening plasticity, which is based on a single set of constitutive parameters, namely two for the elastic behaviour plus eight for the plastic behaviour, which all have a clear and easily understandable physical meaning. In order to assist the convenience of the numerical implementation, the model is defined in a ‘normalized’ stress space in which the stress–strain behaviour does not undergo any strain softening and so certain potential numerical difficulties are avoided. In the first part the multiaxial formulation of the model is described in detail, using appropriate mixed invariants, which rationally combine stress history and stress. The model simulations are compared with some experimental results for tests on granular soils along stress paths lying outside the triaxial plane over a wide range of densities and mean stresses, using constitutive parameters calibrated using triaxial tests. Furthermore, the study is extended to the analysis of the effects induced by the different shapes of the yield and bounding surfaces, revealing the different role played by the size and the curvature of the bounding surface on the simulated behaviour of completely stress- and partly strain-driven tests. Copyright © 1999 John Wiley & Sons, Ltd.     

基于临界状态模型的砂土非共轴本构模拟

传统的砂土本构理论隐含了应力和塑性应变率的共轴条件,无法客观描述主应力轴旋转试验中的非共轴现象,并且当密度和围压变化较大时也不适用。基于材料状态相关砂土临界状态概念,将Pietruszczak和Stolle所提出的砂土本构模型进行了改进,并在模型中引入非共轴塑性流动理论来描述非共轴现象。通过对单剪试验和空心圆柱试验进行数值模拟,表明基于临界状态理论的非共轴模型能够合理描述主应力轴旋转过程中砂土的非共轴变形特性     

Application of a non-coaxial soil model in shallow foundations

The influence of a non-coaxial model for granular soils on shallow foundation analyses is investigated. The non-coaxial plasticity theory proposed by Rudnicki and Rice (J. Mech. Phys. Solids 1975, 23, 371–394) is integrated into a Drucker–Prager model with both perfect plasticity and strain hardening. This non-coaxial model is numerically implemented into the finite-element program ABAQUS using a substepping scheme with automatic error control. The influence of the non-coaxial model on footing settlement and bearing capacity is investigated under various loading and boundary conditions. Compared with the predictions using conventional coaxial models, the non-coaxial prediction results indicate that the settlement of a footing increases significantly when the non-coaxial component of plastic strain rate is taken into consideration, although ultimate footing bearing capacities are not affected significantly. The non-coaxial model has a different effect on footing settlements under different loading and boundary conditions. In general, the discrepancies between coaxial and non-coaxial predictions increase with increasing rotation of principal stresses of the soil mass beneath a footing. It can be concluded that if the non-coaxial component of plastic strain rate is neglected in shallow foundation problems using the finite-element method, the results tend to be non-conservative when designs are dominated by settlement of footings.     

Effect of Temperature on the Strength Parameters at the Plastic Domain for Unconsolidated Sandstones

Heavy and extra heavy oil reservoirs are subjected to high temperatures during thermal recovery processes like SAGD and CSS. In-situ temperature state changes throughout these operations may generate variations in the mechanical properties or strength parameters of the rock such as Young and Bulk’s modulus, internal friction angle, cohesion and dilation angle. This paper presents results about the variation of the elastic modulus and mechanical properties, in both elastic and plastic region due to changes in temperature and effective confining stress for reconstituted samples of Colombian unconsolidated sands. In order to study these changes, several drained triaxial tests are performed in a range of temperature between 50 and 230 °C and a range of effective confining stress between 0.4 and 8.2 MPa, these tests are carried out on reconstituted samples of oil sands recovered in an outcrop in the Magdalena middle valley of Colombia. The experimental results are analyzed with an analytical model based on the constitutive model of Mohr–Coulomb. The analysis is divided in two ways: First, the elastic region is analyzed calculating the variation of the elastic modulus (Young and bulk’s modulus) and the mechanical parameters (cohesion and internal friction angle) with temperature and effective confining stress, these calculus are done in the linear zone of the stress–strain curve. Then, the plastic region is analyzed using Mohr–Coulomb constitutive model with a non-associative flow rule in order to represent in a better way the dilatant/contractive behavior of the material, from this analysis is possible to obtain the behavior of the mechanical properties during plastic strains. The stress–strain curves obtained from the triaxial compression tests show that for the highest and medium effective confining stress (8.2 and 4.0 MPa), there is a hardening strain behavior after the yielding point that implies a loose sand behavior under confining conditions. Meanwhile, for lowest effective confining stress (0.4 MPa) the samples show a softening strain, a typical behavior for dense sands under confining conditions.     

A new closed-form solution for circular openings modeled by the Unified Strength Theory and radius-dependent Young’s modulus

A new closed-form solution is presented for the stress and displacement distribution surrounding circular openings with finite external radii that are subject to uniform internal and external pressures under plane strain conditions. The specific solution for a deep circular tunnel in an infinite rock mass is also provided. It is assumed that the rock mass is elastic–brittle–plastic and governed by the Unified Strength Theory (UST). In the plastic zone, the radius-dependent Young’s modulus (RDM) model and a non-associated linear flow rule were adopted to establish the radial displacement solution. The new closed-form solution obtained in this paper is a series of results rather than one specific solution; hence, it is suitable for a wide range of rock masses and engineering structures. The traditional solutions, which are based on the Mohr–Coulomb failure criterion and the Generalized Twin Shear Stress yield criterion, can be categorized as special cases of this proposed solution. This new solution agrees reasonably well with the results of a borehole collapse test, a secondary development numerical simulation and an additional closed-form solution using the generalized non-linear Hoek–Brown failure criterion. Parametric studies were conducted to investigate the effects of intermediate principal stress, RDM and dilatancy on the results. It is shown herein that the effects of intermediate principal stress and dilatancy are significant; the RDM model is recommended as the optimum approach for calculating radial displacement and support pressure.     

A phenomenological constitutive model for rocks with shear failure mode

In this paper, a phenomenological constitutive model is proposed to simulate the stress–strain behaviours of intact rocks with shear failure mode. The model captures a wide range of behaviours of rock material such as elastic, plastic, strain softening, ‘Class II’, strain localization, elastic modulus degradation, etc. The sensitivity of the stress–strain relation on the parameters is also investigated. Typical results obtained by testing a number of granite and marble specimens are used to validate the proposed model. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical aspects of non-coaxial model implementations

This paper ascertains the reasons for the numerical problems when the yield vertex non-coaxial model is implemented in the finite element analysis to predict the behaviour of complicated geotechnical engineering problems. The numerical problem, reflected in the failure of convergence in the non-linear solutions in ABAQUS, is likely to happen when a smaller non-coaxial plastic modulus is used. It is found that a large non-coaxial influence at the start of elastoplastic loading, which causes the predominance of non-coaxial plastic strain rate over the coaxial plastic strain rate in stress–strain responses, is the reason for the numerical problem. The original yield vertex non-coaxial model is modified to overcome the numerical problem. Instead of a constant non-coaxial plastic modulus in the original yield vertex non-coaxial model, the non-coaxial modulus is made a function of cumulative deviatoric plastic strain. It shows that the modified non-coaxial model retains the functions of the original non-coaxial model. Meanwhile, it overcomes the non-convergence problem when a smaller non-coaxial plastic modulus, representing a larger non-coaxial influence, is used.     

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.

     

从各向异性的角度解释和模拟土的非共轴特性

从发挥面的角度出发,分析论证各向异性是引起岩土材料出现非共轴现象的根本原因,得到与材料力学一致的结论。当共轭的两发挥面与沉积面的夹角不相等时,主应力面上将出现塑性应变增量的切向分量,所以塑性应变增量的主方向与应力的主方向非共轴。按照这一结论,对非共轴的数值模拟,也应当根据各向异性本构模型进行。为考虑各向异性影响新近提出的各向异性变换应力法,改变了各应力分量的相对大小,得到的各向异性变换应力张量与真实应力张量的主方向不一致,因此也能反映非共轴。利用各向异性变换应力法,能够在现有的弹塑性本构模型的框架下,描述土的非共轴现象。以各向异性UH模型为例,预测各种加载条件下的非共轴变形,验证了该方法的有效性。