岩土材料的各向异性强度准则

岩土材料的强度往往表现出很强的各向异性,而已有的各向同性强度准则不能够描述这一特性。提出一个岩土材料的各向异性强度准则。为了描述材料的各向异性,引入了一个由应力张量和组构张量的联合不变量表达的各向异性参数。该参数可以描述加载方向和材料组构方向的夹角。强度准则是基于材料在子午面和偏平面上的破坏特性而建立的,这为描述广义的材料强度各向异性提供了方便。与原各向同性强度准则相比,各向异性强度准则只引入了两个新的模型参数,而且所有的模型参数都可以通过常规的室内试验结果确定。该准则的预测结果与砂土、黏土、天然黏土和岩石的试验结果比较表明,它能够很好地描述岩土材料强度的各向异性     

基于修正应力法的横观各向同性摩尔-库仑准则及被动土压力公式

采用修正应力法考虑各向异性对土的抗剪强度的影响,该方法引入组构张量调整不同方向应力分量的相对大小,使得各向异性土在修正应力空间中等效成各向同性土。用修正应力张量代替真实应力张量,就能将摩尔-库仑强度准则发展至横观各向同性,公式的形式不发生改变,强度参数仍为与加载方向无关的常量。通过对比根据3种不同的修正应力公式得到的内摩擦角与加载方向的关系曲线,分析了组构值对内摩擦角变化规律的影响,在偏平面上绘出了连续规则的强度包线,并预测了不同各向异性岩土类材料的三轴压缩和真三轴试验结果。最后,根据横观各向同性摩尔-库仑准则计算水平沉积地层中挡土墙上的被动土压力,得到了简洁的显式表达式,验证了修正应力法的实用性。     

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

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

宏细观结合考虑主应力轴旋转的砂土破坏特性

采用宏细观结合各向异性破坏准则对主应力轴旋转条件下砂土的破坏特性进行分析。该准则是加载应力、组构各向异性程度和应力与组构几何关系3个因素的函数,可描述细观特性对任意应力旋转角度条件下破坏特性的影响。根据空心圆柱扭剪试验的特点推导一般正交坐标系下主应力轴旋转条件下的破坏关系式,考虑应力与砂土细观组构的几何关系,推导的关系式即可分析该条件下破坏特性。材料为各向异性时,主应力轴旋转造成破坏特性发生变化,细观各向异性程度越小变化越小;材料为各向同性时,则不会造成砂土破坏特性的变化。该式表明主应力轴旋转条件下不同破坏特性存在的根本原因是砂土各向异性的存在。采用空心圆柱试验结果进行验证,结果表明建立的关系式能较好描述不同应力加载角度条件下砂土的破坏特性。初步验证了由于砂土各向异性的存在使得主应力轴旋转造成了不同的破坏规律。     

地震动下考虑各向异性土体-盾构隧道数值模拟

针对层状土体-盾构隧道地震反应分析,引入了横观各向同性弹塑性模型,建立了横观各向同性介质的双渐近-多向透射边界条件。针对盾构隧道抗震设计的特点,基于横观各向同性弹塑性模型,研制了考虑层状土体各向异性和施工开挖效应,适合于地下结构动力计算的弹塑性动力有限元程序。在程序中对不同的材料采用了不同的本构关系和单元形式,并采用了关联流动法则和多种屈服准则,可同时进行横观各向同性土体与地下结构的平面应力、平面应变和轴对称问题的静动力有限元分析。最后,利用所研制的程序对上海地铁二号线石门一路站附近区间隧道在不同超越概率地震波输入下的地震反应进行了计算。结果表明,在层状土体-地铁盾构隧道的抗震分析中考虑土体各向异性的影响是必要的,所提出的计算模型是合理易行的。     

横观各向同性土的简明破坏机制解释

由于沉积等作用使得天然土材料内部颗粒排列表现为各向异性。土材料的各向异性使得材料内各平面滑动的内摩擦角不同。因此,在分析各向异性土的破坏问题中,除了考虑来自外部的应力分布以外,还需考虑土材料内部的强度分布。以横观各向同性土为例,阐述了同时考虑材料外部应力分布与材料内部强度分布下各向异性土的破坏机制。为说明,文中假定横观各向同性土内部各平面的内摩擦角正切值随空间方向而线性分布。外部加载进行时,横观各向同性土材料内部各处的应力不断变化,当其内部某处首次出现的应力达到该处的强度时,土材料发生破坏。进一步通过类比得出：区别于金属,各向同性土不是在最大剪应力面而是在最大剪压比面发生破坏;区别于各向同性土,各向异性土不是在最大剪压比面而是在最大剪压强比面发生破坏。     

各向异性初始应力状态下圆柱孔扩张理论弹塑性分析

传统的圆柱孔扩张理论假定初始应力为各向同性的。在城市地下管道的铺设、隧道工程等的水平掘进施工过程中,圆柱孔受到的初始竖向应力不等于初始水平应力,并不满足各向同性初始应力条件。考虑土体的各向异性初始应力条件,假定土体满足Mohr-Coulomb屈服准则,推导了各向异性初始应力状态下的圆柱孔扩张问题弹塑性解答。算例计算结果表明：土性参数一定时,在相同内压力作用下,各向异性初始应力状态下的圆柱孔扩张形成的塑性区半径大于各向同性初始应力状态下的圆柱孔扩张的塑性区半径,但前者的极限扩孔压力小于后者。     

复杂应力条件下土体各向异性及其建模思路

由于加荷方式不同,土体在复杂应力状态下在各主应力方向上应力-应变关系表现出显著应力各向异性,在常规三轴试验基础上,采用经典弹塑性理论各向同性土体模型对此不能合理描述。通过真三轴试验,总结应力各向异性柔度矩阵规律,结合试验规律进行相应理论研究,用非线性各向异性弹性矩阵代替弹塑性模型的弹性矩阵,用具有各向异性屈服准则的弹塑性模型描述塑性部分,建立非线性各向异性弹性-塑性模型,可以改善柔度矩阵矩阵形态,反映复杂应力状态下土体应力各向异性特征。     

一个与弹性应变能相协调的脆性横观各向同性材料破坏准则

岩土材料广泛存在于自然界中,在重力沉积作用下表现出明显的横观各向同性。为了更好地研究横观各向同性材料的破坏机制,基于Christensen提出的破坏准则应与弹性应变能函数相协调的思想,推导了横观各向同性材料的弹性应变能,通过应变能的分解和整理,分析了应变能中各项所表示的物理意义,确定了级数展开的基,从而使得应变能函数和材料破坏准则在形式上具有一致性,最终提出了一个适用于脆性横观各向同性材料的破坏准则。该准则由5项能够与弹性应变能相协调的部分组成,仅包含5个独立且物理力学意义清晰的参数。基于常规试验给出了5个独立参数的确定方案。通过与Tsai-Wu准则及已发表的多组试验数据的对比,验证了准则的正确性。     

土体强度各向异性基坑的抗隆起稳定性分析

天然土体通常为非均质、各向异性土体,目前常用的基坑抗隆起计算方法多是针对均质各向同性土体的。采用非均质非线性各向异性岩石破坏准则,考虑土体原生各向异性,假定Prandtl破坏模式,采用地基承载力模式进行基坑的抗隆起稳定性分析。根据滑移线理论计算承载力,利用有限差分法结合边界条件编制出相关计算程序,并定义基坑抗隆起稳定性系数。研究表明,各向异性比的提高会降低安全系数;内摩擦角和插入比的增大可有效提高基坑的稳定性。     

Anisotropic UH model for soils based on a simple transformed stress method

A simple method called anisotropic transformed stress (ATS) method is proposed to develop failure criteria and constitutive models for anisotropic soils. In this method, stress components in different directions are modified differently in order to reflect the effect of anisotropy. It includes two steps of mapping of stress. First, a modified stress tensor is introduced, which is a symmetric multiplication of stress tensor and fabric tensor. In the modified stress space, anisotropic soils can be treated to be isotropic. Second, a TS tensor is derived from the modified stress tensor for the convenience of developing anisotropic constitutive models to account for the effect of intermediate principal stress. By replacing the ordinary stress tensor with the TS tensor directly, the unified hardening model is extended to model the anisotropic deformation of soils. Anisotropic Lade's criterion is adopted for shear yield and shear failure in the model. The form of the original model formulations remains unchanged, and the model parameters are independent of the loading direction. Good agreement between the experimental results and predictions of the anisotropic unified hardening model is observed. Copyright © 2016 John Wiley & Sons, Ltd.     

A yield criterion for isotropic and cross‐anisotropic cohesive‐frictional materials

This paper proposes a yield and failure criterion for cohesive and frictional materials. The function is given by the combination of a Lode dependence for the behaviour in the deviatoric plane and a meridian function for the pressure‐dependent behaviour. A variety of shapes can be achieved with the proposed criterion including Lode dependences which are able to reproduce the behaviour of isotropic and cross‐anisotropic materials in the deviatoric plane. The criterion is validated through the comparison with experimental data based on multiaxial experimental tests on clays, sands, rocks and concrete. Finally, the convexity of the criterion is analysed and discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

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.

     

Modelling the elastic behaviour of granular materials

A review of the literature indicates that the elastic behaviour of granular materials is isotropic and that Poissony's ratio is constant, whereas Young's Modulus, the bulk modulus and the shear modulus vary with the mean normal stress and the deviatoric stress. A nonlinear, isotropic model for the elastic behaviour is developed on the basis of theoretical considerations involving the principle of conservation of energy. Energy is therefore neither generated not dissipated in closed-loop stress paths or in closed-loop strain paths. The framework for the model consists of Hooke's law, in which Poission's ratio is constant and Young's modulus is expressed as a power function invlving the first invariat of the stress tensor and the second invariant of the deviatoric stress tensor. The characteristics of the model are described, and the accuracy is evaluated by comparison with experimental results from triaxial tests and three-dimensional cubical triaxial tests with a variety of stress paths. Parameter determination from unloading–reloading cycles in conventional triaxial compression tests is demonstrated, typical parameter values are given for granular materials and extension of the model to soils with effective cohesion is described.     

A micromechanics-based model for estimating localized failure with effects of fabric anisotropy

This paper presents a micromechanics-based approach to investigate the effects of fabric anisotropy on the behavior of localized failure in granular materials. Based on a micromechanical analysis, the origin of deviatoric stress is decomposed into two components: contact force anisotropy and fabric anisotropy. Using a micro–macro approach, the back stress is interpreted as an contribution to the change of the fabric’s principal direction. The evolution of the back stress is deduced from the stress–fabric relationship and determined with reference to the deviation of the principal directions between the rate of the reduced stress tensor and the actual reduced stress tensor. With this micro–macro framework, a mixed (isotropic–kinematic) hardening model is developed based on the classical isotropic hardening theory. A laboratory simple shear test is first analyzed to validate the proposed model and illustrate the kinematic-hardening mechanism which is usually displayed under non-proportional loading. The analysis further focuses on the anisotropic aspect of localized failure. It has been discovered that the fabric anisotropy can play an important role in the occurrence of shear banding. An increasing degree of fabric anisotropy tends to delay the initiation of the strain localization and result in higher failure strength. The effects of fabric anisotropy have also been illustrated by comparing the theoretical predictions and measured results on the shear band inclination angle, shear strain level and dilatancy at bifurcation.     

软黏土的各向异性临界状态模型

剑桥模型只适用于正常固结软黏土,不能描述不等向固结土的应力-应变行为的各向异性特性。基于剑桥模型,在其椭圆屈服面中引入各向异性张量和一个形状参数,建立了一个各向异性屈服面,提出了一个适用于等向和不等向固结软黏土的本构模型。各向异性张量的初始值由初始固结应力状态确定,其演化过程由一个与塑性剪应变和塑性体应变都有关的硬化法则描述。形状参数的引入保证了各向异性屈服面的灵活性和适应性。通过对Boston Blue黏土、高岭土和Otaniemi黏土的三轴试验结果的模拟,验证了模型的模拟能力。     

控制应力路径散体材料真三轴试验强度特征的离散元模拟与细观机制分析

通过一系列真三轴离散元数值试验,模拟了不同应力路径下的等b试验中散体材料的强度特征。根据模拟结果详细地分析了三维应力条件下中主应力和应力路径对散体材料峰值强度的影响,研究了峰值摩擦角、峰值应力比的变化规律,并根据真应力的概念和组构张量的演化结果分析了散体材料的强度成因。研究表明,在不同类型的数值试验中峰值偏应力随b参数的变化规律不同,但采用初始围压归一化后的应力-应变曲线规律一致。峰值强度线的斜率只与b值有关而与应力路径无关,且随着b值的增加,峰值应力比qf /pf逐渐减小,数值模拟结果与室内试验结果吻合较好;随着应变的发展,数值试样的组构也随之发生变化,产生了明显的应力诱发各向异性;散体的强度为颗粒摩擦及材料各向异性共同作用的结果;理论上,组构比-应力比坐标系中破坏点位置仅取决于颗粒摩擦角 ,而数值模拟结果与理论值的差异源于颗粒间咬合和滚动摩擦的影响,其影响与颗粒表面摩擦系数有关,也受空间应力状态的影响。     

A critical state sand plasticity model accounting for fabric evolution

Fabric and its evolution need to be fully considered for effective modeling of the anisotropic behavior of cohesionless granular sand. In this study, a three‐dimensional anisotropic model for granular material is proposed based on the anisotropic critical state theory recently proposed by Li & Dafalias [2012], in which the role of fabric evolution is highlighted. An explicit expression for the yield function is proposed in terms of the invariants and joint invariants of the normalized deviatoric stress ratio tensor and the deviatoric fabric tensor. A void‐based fabric tensor that characterizes the average void size and its orientation of a granular assembly is employed in the model. Upon plastic loading, the material fabric is assumed to evolve continuously with its principal direction tending steadily towards the loading direction. A fabric evolution law is proposed to describe this behavior. With these considerations, a non‐coaxial flow rule is naturally obtained. The model is shown to be capable of characterizing the complex anisotropic behavior of granular materials under monotonic loading conditions and meanwhile retains a relatively simple formulation for numerical implementation. The model predictions of typical behavior of both Toyoura sand and Fraser River sand compare well with experimental data. Copyright © 2013 John Wiley & Sons, Ltd.     

A generalized nonlinear failure criterion for frictional materials

A generalized nonlinear failure criterion formulated in terms of stress invariants is proposed for describing the failure characteristics of different frictional materials. This failure criterion combines a power function and a versatile function in the meridian and deviatoric plane, respectively, which is a generalization of several classic criteria, including the Tresca, Drucker–Prager, Mohr–Coulomb, Lade–Duncan and Matsuoka–Nakai failure criterion. The procedure for determination of the strength parameters was demonstrated in detail. Comparisons between the failure criterion and experimental results were presented for uncemented/cemented Monterey sand, normally consolidated Fujinomori clay, rockfill, concrete, Mu-San sandstone and granite, which reveal that the proposed failure criterion captures experimental trend quite well.