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.     

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.     

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.     

Non‐coaxial elasto‐plasticity model and bifurcation prediction of shear banding in sands

Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 John Wiley & Sons, Ltd.     

Investigation of influence of particle characteristics on the non‐coaxiality of anisotropic granular materials using DEM

As a result of deposition process and particle characteristics, granular materials can be inherently anisotropic. Many researchers have strongly suggested that the inherent anisotropy is the main reason for the deformation non‐coaxiality of granular materials. However, their relationships are not unanimous because of the limited understanding of the non‐coaxial micro‐mechanism. In this study, we investigated the influence of inherent anisotropy on the non‐coaxial angle using the discrete element method. Firstly, we developed a new discrete element method approach using rough elliptic particles and proposed a novel method to produce anisotropic specimens. Secondly, the effects of initial specimen density and particle characteristics, such as particle aspect ratio A _m, rolling resistance coefficient β , and bedding plane orientation δ , were examined by a series of biaxial tests and rotational principal axes tests. Findings from the numerical simulations are summarized as follows: (1) the peak internal friction angle ? _p and the non‐coaxial angle i both increase with the initial density, A _m and β , and they both increase initially and then decrease with δ in the range of 0–90°; (2) among the particle characteristics, the influence of A _m is the most significant; and (3) for anisotropic specimens, the non‐coaxial angle can be calculated using the double slip and rotation rate model. Then, an empirical formula was proposed based on the simulation results to depict the relationship between the non‐coaxial angle and the particle characteristics. Finally, the particle‐scale mechanism of non‐coaxiality for granular materials was discussed from the perspective of energy dissipation. Copyright © 2016 John Wiley & Sons, Ltd.     

基于统一强度理论深埋圆形隧道围岩的剪胀分析

基于连续介质理论中岩体的剪胀角与围压和塑性剪切应变密切相关,隧洞周边岩体的应力状态因开挖卸荷而发生应力重分布,迫使其围压由原地应力逐渐衰减,塑性剪切应变不断增加,引起剪胀效应呈非线性变化。首先,基于统一强度理论和非关联流动法则,将潜在塑性区围岩按等围压释放划分为若干同心圆,提出了考虑中间主应力和非线性剪胀性的有限差分法,计算应变软化围岩的力学问题,并以实例验证其正确性。其次,通过参数分析,研究塑性区内岩体的剪胀角受中间主应力、临界软化系数和支护力的影响规律。研究结果表明,中间主应力主要影响剪胀角的峰值,随着中间主应力效应增加,剪胀峰值增加;临界软化系数主要影响剪胀角的变化率,随着临界软化系数的增加,剪胀角变化缓慢;中间主应力和临界软化系数共同影响塑性区剪胀角的变化;随着支护力的增加,洞壁处的剪胀角增加;双剪强度理论计算的位移值较小,应谨慎采用,同时采用Mohr-Coulomb强度准则时可以适当考虑围岩的承载潜力。     

平面应变下原状黄土的剪切带研究

针对黄土地区工程建设中的平面应变问题,使用改进的西安理工大学真三轴仪进行了不同侧向固结压力、不同含水率下原状黄土的平面应变剪切试验,得出了剪切破坏时剪切带倾角试验值,根据应力-应变关系曲线求出了剪切带倾角的Mohr-Coulomb、Roscoe理论解,分析了相应控制因素下原状黄土强度、剪切带破坏模式与倾角的变化。结果表明：应力-应变曲线随含水率、侧向固结压力的增大从原生软化型向次生硬化型转变;剪切带破坏模式随侧向固结压力的增大依次呈现单一剪切带、锥型剪切带和X型剪切带;随着含水率的增大,黏聚力、内摩擦角减小导致Mohr-Coulomb剪切带倾角减小,剪胀作用增强导致Roscoe剪切带倾角增大;随侧向固结压力增大,动内摩擦角增大及剪胀作用增强,导致Mohr-Coulomb、Roscoe剪切带倾角增大。     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

Undrained shear band in water saturated granular media: A critical revisiting with numerical examples

By taking into account the continuity of stresses and pore pressure across the boundaries of a shear band, it is theoretically shown that incipient localization may take place in saturated loose contractive sand in the regime of decreasing deviator stress under undrained conditions. The undrained shear band orientation primarily depends on the Poisson's ratio and the dilatancy characteristics of the material. However, under strict isochoric constraint, localization is precluded in dense dilative sand with deviator stress increasing only. If any local volume change takes place with a certain mechanism, undrained shear band would be possible in dense sand, most likely in the regime of decreasing friction angle after the peak friction angle is mobilized. Numerical examples are given to demonstrate the variation of shear band inclinations with void ratios and initial consolidation pressures. Copyright © 2011 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.     

基于损伤控制方法的大理岩非线性剪胀特性试验研究

剪胀角是描述岩石体积膨胀扩容的常用参数,在非关联流动法则中,连续介质理论通常假设剪胀角为0;在关联流动法则中,其值恒定且等于内摩擦角。岩石三轴压缩全过程体应变曲线表明,其体积剪胀性依赖于围压和塑性参量,破坏过程中不仅其特征强度随围压和塑性参量呈非线性变化,而且剪胀特性也表现出非线性特征。基于塑性力学理论,针对锦屏大理岩损伤控制的全过程三轴加、卸载试验,采用双参数非线性函数拟合方法建立了能同时考虑围压效应和塑性参量的非线性剪胀角模型。结果表明,对于大理岩、中硬岩,在破坏过程中扩容行为强烈依赖围压和岩石塑性参量,均表现出先快速增加至峰值后,随着塑性变形增加逐渐减小的非线性演化规律。提出的双参数非线性剪胀角模型很好地描述了岩石破坏过程中的体积扩容特性,其结果对于研究地下工程围岩应力变化诱发的围岩剪胀破坏机制、体积扩容膨胀区范围预测和围岩支护的合理设计均具有一定的理论和工程应用价值。     

基于抗转模型的颗粒材料宏−细观关系研究

接触模型的宏?细观参数标定是成功使用离散元方法的关键。在离散元的接触模型中线性接触模型与抗转线性接触模型均可用于模拟砂性土的力学行为,其中抗转线性接触模型在模拟密砂的剪胀性方面具备优势。采用抗转线性接触模型对室内密实砂土三轴试验进行了离散元模拟,验证了抗转线性接触模型的可靠性;进而系统分析了颗粒间摩擦系数、刚度比和抗转动系数等细观参数与砂土峰值内摩擦角、残余内摩擦角、峰值剪胀角等宏观参数的相关关系并进行了验证;揭示了偏应力作用下,细观参数对密实砂土试样内部剪切带宽度与倾角变化的影响规律,提出了考虑剪胀角的剪切带倾角经验公式。通过研究建立了抗转线性接触模型宏?细观参数的量化关系并给出了标定参数的具体流程图,提出了快速标定宏观参数的方法并应用实例进行了验证,为采用抗转线性接触模型精准模拟密实砂土的力学特性提供依据。     

Numerical simulation of the shear behavior of rockfills

A strain hardening model is proposed for simulation of the shear behavior of rockfills, including the non-linear stress–strain relationship, stress-dependence of stiffness, non-linear strength envelope, and particularly the shearing contraction and dilatancy. This model utilizes the Mohr–Coulomb yield criterion, which is defined by the mobilized friction angle varying with plastic octahedral shear strain. The development of Rowe’s stress dilatancy theory is reviewed, with its limitation for rockfills highlighted. A modified solution is then proposed to derive the mobilized dilatancy angle of rockfills. Compared with large-scale triaxial testing results, it has been demonstrated that this simple model is capable of predicting the shear behavior of rockfills with fair accuracy.     

考虑受压侧岩体反力非线性作用的锚杆抗剪理论

针对锚杆受压侧岩体或砂浆体反力非线性作用及结构面的剪胀效应问题,基于经典梁理论推导锚杆轴力与轴向变形及横向剪切力与横向变形的理论公式,建立了锚杆抗剪力计算公式。通过加锚结构面直剪试验验证理论计算有效性,并分析结构面剪胀系数、围岩强度、锚杆安装角（倾角）对锚杆变形和抗剪力的影响。结果表明：锚杆抗剪理论计算与室内试验结果吻合较好;结构面剪胀系数越大,越能较快调动锚杆抗剪作用,相反锚杆塑性强化特征越不明显,改善加锚结构面的阻滑抗剪作用,主要依靠结构面固有抗剪强度;随着围岩强度降低,锚杆需经一定变形才能发挥较大抗剪作用,而随着围岩强度增大,锚杆将迅速达到屈服状态,并且锚杆由轴向张拉破坏逐渐转为拉剪破坏;锚杆最优安装角随结构面内摩擦角增大而增大,依据实际工程中结构面内摩擦角取值范围,可估算锚杆最优安装角为30°～68°。     

Discrete element method analysis of non‐coaxial flow under rotational shear

This study focuses on non‐coaxial flow behavior of cohesionless soil undergoing cyclic rotational shear, with a special interest in the effects of particle‐scale characteristics. To this end, we perform a series of 2D discrete element simulations with various particle shapes, inter‐particle coefficient of friction, initial density, and stress ratios. The validity and efficacy of the numerical model is established by systematically comparing numerical simulation results with existing laboratory testing results. Such comparison shows that the numerical simulations are capable of capturing mechanical behavior observed in laboratory testing under rotational shear. We further demonstrate and quantify a strong yet simple relationship between the deviatoric part of the normalized strain increment and the non‐coaxial angle, denoted by and ψ, respectively. This quantitative correlation between ψ and is independent of applied stress ratio, initial and current void ratio, and the number of cycles applied, but dependent on the principal stress orientation and particle‐scale characteristics. At the same , specimens with higher inter‐particle friction angle or smaller particle aspect ratio show greater non‐coaxial angles. A simple model is able to fit this ψ‐ relationship well, which provides a useful relationship that can be exploited in developing constitutive models for rotational shearing. Copyright © 2014 John Wiley & Sons, Ltd.     

应变软化扩容对含随机缺陷岩石的渐进变形及破坏前兆特征的影响

在单轴平面应变压缩条件下, 采用FLAC模拟了剪切扩容对含随机缺陷岩石破坏前兆及变形特征的影响。密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则, 破坏之后呈现应变软化-理想塑性行为。缺陷在破坏之后经历理想塑性行为。随着轴向应变的增加, 试样内部破坏的单元数目增加, 直到达到一个常数, 该常数随着扩容角的增加而增加。当扩容角较高时, 计算得到的泊松比在峰前就可以超过0.5;剪切扩容于峰前发生; 变形后试样的最终体积大于初始体积。剪切局部化(导致了毗邻块体之间的相对滑动)及剪切扩容(发生于剪切带内部)是非零扩容角试样峰后体积膨胀的原因。在峰前, 通过观察剪切应变增量、破坏的单元数目、侧向应变、计算得到的泊松比及体积应变可以发现, 扩容角越高, 试样破坏的前兆越明显。在低扩容角时, 由于弯曲的剪切带边界, 试样内部充分发展的剪切带的倾角比较分散, 剪切带的倾角更接近Arthur理论。      

土石混合料大型直剪试验的颗粒离散元细观力学模拟研究

土石混合料作为一种特殊的岩土介质越来越受到国内外众多研究者的重视。基于3维颗粒离散元PFC3D,建立了土石混合料直剪试验模型,进行了不同含石量、不同岩性的土石混合料直剪试验模拟研究。颗粒离散元模拟结果表明,土石混合料的石料岩性和含石量在很大程度上控制了土石混合料的抗剪强度特性。硬岩混合料的摩擦角普遍比软岩混合料大6°～ 7°,含石量为60%～80%时达到最大。土石混合料的剪切面不再是一个平面,其起伏度随含石量增加而增大。剪切过程中软岩混合料在低正应力下表现为剪胀,高正应力下表现为剪缩,并产生软化现象,硬岩混合料表现为剪胀和塑性;软岩土石混合料剪切过程中能量以应变能和动能为主,而硬岩土石混合料的能量以摩擦能和动能为主。     

Bearing capacity of non-associative coaxial granular materials by upper bound limit analysis and finite elements

In this paper, an upper bound estimate of the limit load on non-associative coaxial granular materials is presented. The kinematic approach of the upper bound limit analysis has been utilised. The failure mechanism is assumed to coincide with the direction of the shear bands at every point throughout the body. The shear band orientation in non-associative coaxial materials, i.e. those with the same major principal stress and major principal strain increment directions, can be found based on the angle of dilation and the major principal stress direction. Therefore, having known the stress field at limiting equilibrium, the orientation of the shear bands and hence, the failure mechanism can be obtained. In this study, the stress field is first determined by the method of stress characteristics. Then, the finite element interpolation technique is used to interpolate the stress field and to find the orientation of the shear bands at every point within the field. Once the failure mechanism and the stress state at every point along velocity discontinuities have been found, the upper bound limit analysis has been performed to estimate the limit load.     

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

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