广义塑性力学的加卸载准则与土的本构模型——广义塑性力学讲座(3)

首先介绍了广义塑性力学的加卸载准则,该准则能准确判断各应变分量的加卸载状态,可以方便地应用于数值分析。然后提出了基于广义塑性力学的土本构模型。通过试验,给出了硬化压缩土和硬化剪胀土的屈服条件与计算参数,并由算例说明了其合理性。     

双剪双参数强度准则硬化函数的探讨

目前双剪强度理论在工程实际中得到了广泛的应用。基于双剪理论和相关流动法则,采用各向同性强化法则,选用等效塑性应变作为强化参数,从一致性条件出发,推导出双剪双参数强度准则本构模型硬化函数的显式表达式,并利用单轴加载试验数据得到硬化函数与硬化参数之间的微分关系,推导出考虑硬化的双剪双参数强度准则本构模型的弹塑性刚度矩阵,给出了求解双剪强度理论硬化函数的另一思路。通过工程算例验证,得出一些结论对实际计算具有一定的指导意义。     

岩土材料应力路径无关硬化参量的构成方法

在分析砂土试验结果的基础上,揭示了基本硬化内参量（塑性体积应变、塑性剪应变）变化的应力路径相关性,提出了应力路径无关硬化参量的一般表达式。并指出了塑性功硬化参量及其提出的黏土和砂土统一硬化参量均是所提一般表达式框架下的特例。还将该统一硬化参量直接和修正剑桥模型结合,并与试验结果对比,表现出描述剪胀特性的简单性和有效性。     

一个考虑循环荷载作用的简化模型

基于塑性硬化模量场理论和多重屈服面模型,结合各向同性硬化准则和移动硬化准则,在新的应力空间建立了一个新型不排水循环荷载作用下的多屈服面模型,并推导了一个适合三轴试验的简化的多屈服面模型。在此基础上,结合一个循环荷载作用下的动孔压模型,进行了饱和软黏土的动三轴模拟试验。结果表明,文中建立的多屈服面模型能够较好地模拟循环三轴试验、直剪试验和平面应变条件下的试验。     

高低压下不同力学特性的砂土统一模型

为构建能够反映砂土高低压下不同力学特性的统一模型,分析了砂土在较大的压力范围内的力学试验结果并获取其强度、等向压缩以及临界状态特性参数。通过引入应力路径相关因子来修正塑性应变增量中与应力路径相关的部分,从而使得模型硬化参量能够反映密实砂土在常压下的剪胀特性。同时,基于砂土的临界状态特性提出能够与砂土内部状态相对应的潜在状态面概念,由屈服面与潜在状态面间的动态关系确定加载过程中的动态密实参数及潜在强度,进而使得硬化参量也能够反映密实砂土在常压下的软化特性及高压下的剪缩、硬化特性。分析模型屈服面及潜在状态面间的演化规律并对不同压力等级下的砂土受荷力学行为进行模拟预测,证实了该模型能够反映密实砂土常压下剪胀软化及高压下剪缩硬化的特性。     

基于热力学的临界状态模型的返回映射算法及实现

相对于其他临界状态模型,基于热力学的临界状态（TCS）模型不需要引入塑性势假设,能够自动满足热力学定律。通过对TCS的修正,使其能够模拟初始K0固结的影响,利用了返回映射算法进行了TCS模型的ABAQUS本构二次开发,通过与ABAQUS内嵌的修正剑桥模型（MCC）计算结果的对比,证实了该程序的可靠性。对模型中的2个控制屈服面形状的参数进行了讨论,分析了它们对应力-应变关系与剪胀关系的影响,修改TCS模型参数可以实现非椭圆的屈服面,进而拓展了模型的适用范围,不同的参数对于屈服面形状和大小的影响也有所不同。同时还比较了采用考虑K0固结及旋转硬化的TCS模型与不考虑K0固结及旋转硬化的MCC模型所得到的应力-应变关系与剪胀关系的差异。对于真实土体,K0固结及旋转硬化是土体的基本力学特性。最后证明了TCS模型较MCC模型（模型中没有考虑旋转硬化与K0固结对硬化规律的影响）可以很好地模拟该特性,所以更加有效。     

剪胀性砂土本构模型的研究

对于中密砂和密砂,剪胀现象是一个比较显著的特征,相变线是描述剪胀现象的一条特征线,它比较容易确定。因此,以相变线作为状态参考线,提出并定义了一种状态参量,在已有弹塑性模型基础上,将状态参量引入剪胀方程和塑性硬化模量表达式中,建立了一个新的考虑砂土剪胀性的本构模型。该模型尤其适合于中密砂和密砂,对于松砂则退化为原弹塑性模型。模型形式较简单,包含9个材料参数,且比较容易确定。数值模拟结果与试验数据吻合较好。     

考虑硬化软化和剪胀特性的绿泥石片岩力学模型

绿泥石片岩是一种典型的软岩,开挖过程中出现的围岩大变形和塌方对工程安全危害极大。为了深入研究绿泥石片岩的力学特征并建立其合理的力学模型,首先进行了不同围压下的三轴压缩试验,并考虑应力状态对塑性演化的影响,定义了新的内变量;分析了凝聚力和内摩擦角随内变量的演化规律,结果表明,凝聚力先近似呈线性减小,而后近似呈抛物线递减至残余破坏;而内摩擦角一直近似呈抛物线递增;研究了绿泥石片岩的剪胀特性,分析了剪胀角随内变量演化规律,结果表明,剪胀角与内摩擦角呈相反的演化趋势。考虑硬化、软化和剪胀特性,建立了绿泥石片岩的力学模型,对室内三轴压缩试验结果的模拟分析表明,该力学模型可以较好地描述绿泥石片岩的硬化、软化规律和剪胀等性质,为工程安全性分析提供了基本的力学模型,对于类似软岩力学性质的研究也具有重要的参考意义。     

基于相似性原理超固结土不排水扩张半解析解

在实际工程中,土体往往因卸载、再加载等复杂应力路径而处于超固结状态,而现有的圆孔扩张问题的计算模型往往不能反映超固结土中剪胀、软化等一些特殊性质。为了解决这一问题,基于相似性原理和统一硬化（UH）模型,结合相关联的流动法则和大变形理论,采用相似求解技术求解了超固结土不排水扩张问题的半解析解答。通过理想化算例分析了圆孔扩张挤土产生的应力和孔压响应,并通过分析不同超固结比OCR的土体应力路径的变化规律,讨论了UH模型的适用性。结果表明：对于轻超固结土,空腔周围土体孔压在塑性区沿径向单调递减,随着OCR增大,塑性区内孔压分布呈现出“S”形的趋势,孔壁附近的孔压逐渐减小,孔壁周围甚至出现负孔压。随着OCR增大,压力?扩张曲线收敛变慢。在扩孔过程中正常固结土一直处于剪缩硬化阶段。而对于超固结土,土体则经历了临界状态→剪胀硬化阶段→临界（特征）状态→剪缩硬化阶段。该研究成果不仅丰富了相似求解技术的应用,而且为超固结土中桩基承载力、隧道围岩变形预测和原位测试参数等岩土工程问题的计算提供了理论依据。     

考虑剪胀效应的深埋圆形隧道围岩应变软化弹塑性解

隧道围岩非线性体积膨胀影响施工安全,如何正确评价隧道围岩塑性区域内扩容机制非常重要。多数研究从理论上推导深埋圆形隧洞应变软化围岩应力-应变场时仅采用恒定或线性变化剪胀模型,针对此不足,提出了一种基于有限差分法的分析方法,能够合理地考虑围岩非线性剪胀效应及其应变软化特性。利用已有的研究成果,验证了该方法的合理性,并进一步探讨了不同质量石灰岩与支护压力下剪胀系数在围岩塑性区域内的影响因素,比较了恒定与非线性变化剪胀模型下围岩变形的不同。结果表明:对于地质强度指标(GSI)较小、质量较差的岩体,塑性区域内主要由围压控制剪胀效应程度;恒定与非线性变化剪胀模型下围岩洞壁变形差别显著。     

车辆荷载作用下软土地基塑性变形的计算

对于低路堤公路,合理地预测车辆荷载引起的地基塑性累积变形对于低路堤公路的设计有着重要的意义。从机制上讲,可将地基的塑性应变分为塑性体积应变和塑性剪切应变。其中对于塑性体积应变,提出采用不排水条件下的残余孔压模型与一维固结理论相耦合进行计算;对于塑性剪切应变,考虑了塑性体积应变引起的土体硬化,提出了修正的塑性剪切应变计算模型。然后基于所提模型,对萧山黏土的不排水与排水条件下的动三轴试验进行模拟,初步验证了所提模型的有效性。对于实际工程问题,首先运用数值的方法求得车辆荷载在地基中引起的附加应力,然后运用所提模型计算各土层的塑性应变,再结合分层总和的方法求得地基的累积塑性变形。以日本佐贺机场1号测点为例,说明了模型在工程边值问题中的有效性。     

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

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

弹塑脆性岩土材料中柱形孔扩张大应变解

在应力跌落模型基础上引进软化阈值,建立弹塑脆性模型,模拟岩土材料的脆性软化特性,采用双剪统一强度理论和不相关联流动法则,考虑软化区的大应变,对柱形孔扩张问题进行求解,得出了扩张率-扩张压力的解析表达式。算例分析表明,软化阈值对扩张变形和塑性区范围影响显著,软化阈值越小,扩张变形和塑性区范围越大,且软化区占整个塑性区的比重越大;是否考虑中间主应力对屈服的影响,对解答也有明显的影响。     

一种旋转塑性势面模型及非关联塑性流动法则

为了较好地描述软土塑性应变发展规律,提出了一种改进的塑性流动模型。该模型采用了与屈服函数形式相同,但具有一定倾角 的塑性势函数。土体在变形过程中,塑性流动方向会依赖于塑性势面的旋转而变化,直至达到破坏状态。通过对常规三轴试验结果的分析可以发现：在剪切过程中,塑性势面旋转角的初值 与终值 较为稳定,不受围压变化影响。在此试验观察基础上,引入了归一化的旋转角参数 以及描述土体应力状态的参数 ,在采用蛋形势函数的情况下二者具有良好的分段线性关系。利用该关系,建立了改进的塑性流动法则,只需要2个额外的模型参数。对所提出的塑性流动模型进行了验证,计算结果表明该模型能较好地反映塑性应变的变化趋势。     

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.     

Analysis of double shearing in frictional materials

The concept that the flow of granular materials is governed by shear on certain critical planes is used to formulate an elastic–plastic model. When the elastic strains are neglected, the Double Shearing model becomes identical to the rigid-plastic model of de Josselin de Jong which he named the Double Sliding model. After a discussion of the model, the small strain formulation is used to describe general boundary-value problems. A comparison is then made between the coaxial Mohr–Coulomb yielding and the Double Shearing model on the basis of several numerical simulations. The non-coaxiality of the axes of principal plastic strain rate and principal stress in the Double Shearing model leads to essentially different behaviour for stress rotations in comparison to the coaxial model. The Double Shearing model predicts in general lower limit loads because, for a given state of stress, it allows for several possible directions of plastic flow rather than a unique direction which derives from a plastic potential.     

Poroplastic modelling of the excavation and ventilation of a deep cavity

A plastic deviatoric model with hardening is developed on the basis of geomechanical tests performed in the saturated case on low permeable porous material such as argillite. This model is a generalized Mohr–Coulomb plastic criterion combined with a Drucker–Prager plastic potential and the hardening parameter is the plastic distortion. Three different hardening functions have been introduced on the basis of triaxial tests: an increase of friction angle, a decrease of cohesion after a threshold and a contractancy to dilatancy transition for volumetric plastic strain. This plastic model has been adapted to the partially saturated case. The effective stress is expressed thanks to the equivalent interstitial pressure π. Numerical results are presented for the excavation and monotonous ventilation of a deep cylindrical cavity. A first plastification due to excavation is followed by a second one due to desaturation. The extent of the non-saturated zone provokes an extent of a plastic zone in the rock mass. Analysis shows that the origin of the plastification can be found in the deviatoric stresses because mean effective stresses are compressive during drying. Copyright © 2004 John Wiley & Sons, Ltd.     

一个基于热力学的土体本构模型

修正剑桥模型假设体积塑性功的一半转化为耗散能,一半转化为由塑性变形约束的自由能,但通常体积耗散能和塑性体积自由能并不相等。基于热力学理论,通过引进间隔应力比来表征体积塑性功和塑性体积自由能的比例关系,构建了一个土的体积硬化模型。由于岩土材料耗散函数与所处的应力状态有关,因此,需要应用非关联流动法则,修正剑桥模型是适合应用相关联流动法则的一个特例,也是所提出模型的一个特例。     

A discrete approach for anisotropic plasticity and damage in semi-brittle rocks

In this paper, we propose an anisotropic plastic damage model for semi-brittle geomaterials based on a discrete thermodynamic approach. The macroscopic plastic deformation is generated by frictional sliding of weakness planes. The evolution of damage is related to growth of such weakness planes. The local frictional sliding in each family of weakness planes is described by a non-associated plastic model taking into account material softening and volumetric dilatancy. The damage evolution is coupled with plastic deformation and modelled by an isotropic damage criterion. The proposed model is applied to modelling mechanical responses of typical sandstone under different loading paths. There is good agreement between numerical predictions and experimental data. Further, the anisotropic distributions of plastic deformation and induced damage are analysed and discussed.     

Non‐normality and induced plastic anisotropy under fractional plastic flow rule: a numerical study

In this paper, an implementation of fractional plastic flow rule in the framework of implicit and explicit procedures is under consideration. The fractional plastic flow rule is obtained from a generalisation of the classical plastic flow rule utilising fractional calculus. The key feature of this new concept is that in general, the non‐associative flow is obtained without necessity of additional potential assumption. If needed, the model can cover the anisotropy induced by plastic deformation. Illustrative examples showing the unusual flexibility of this model are also presented. Copyright © 2015 John Wiley & Sons, Ltd.