工程学院科技论文报告会参会论文基金项目：国家自然科学基金项目(40872175)作者简介：李男，女，1989年生，中国地质大学（武汉）工程学院工程地质专业052073班本科生。

指导老师：胡新丽，中国地质大学（武汉）工程学院教授。

紫红色泥岩蠕变损伤本构关系研究李男，徐辉，胡新丽，马越野(中国地质大学(武汉)工程学院,湖北 武汉430074)摘 要：岩体的蠕变性是影响边坡变形与长期稳定性的重要因素之一。

以湖北巴东地区的五里堆2号滑坡为工程背景，对滑坡内典型的软弱夹层紫红色泥岩，开展室内岩石剪切蠕变试验，研究其蠕变变形特性，进行了蠕变本构模型的辨识，最终选定西原模型并通过计算得到其蠕变参数，结果显示西原模型能够较好地模拟泥岩的初始和稳定蠕变阶段，但无法模拟加速蠕变阶段。

于是，在西原模型的基础上，引入损伤演化方程，从而建立了可以描述该泥岩全过程蠕变变形的损伤模型，并具体介绍了该模型参数的求解方法，拟合结果证明了新模型的正确和合理性。

关键词：岩体力学；泥岩；蠕变；损伤本构模型；参数求解方法STUDY ON CREEP DAMAGE CONSTITUTIVE RELATION OFFUCHSIA MUDSTONELI Nan, XU Hui, HU Xin-li(College of Engineering, China University of Geosciences, Wuhan, Hubei, 430074, China)Abstract: The time-dependent behavior of rock mass is one of the important factors that influence the deformation and long term stability of landslide. As an engineering example, Wu-li-dui landslide is studied ，which is located in Ba-dong city, Hubei Province, with typical soft interlayer- mudstone. The creep deformation characteristics of rock are studied through laboratory shear creep tests. Then, the Xiyuan creep constitutive model is chosen from many ones, and the parameters are identified accordingly, the results show that Xiyuan model can simulate initialized and stable creep deformation, but cannot simulate accelerated creep deformation. Therefore, a new damage developing function is tumbled in Xiyuan creep model to found the damage creep model which can describe whole-process creep deformation characteristics, and the methods to calculate the model parameters are also given in detail. The results indicate that the damage creep model is very correct and reasonable.Key words: rock mechanics; mudstone; creep; damage constitutive model; parameter calculating method1.引言工程实践与研究表明，岩石流变是岩石地下工程、构筑物基础、边坡及滑坡产生大变形乃至失稳的重要原因之一。

目前，国内外因对岩土体流变特性研究不够而导致的延误施工甚至工程失败的先例不胜枚举，如1963年意大利瓦依昂(Vajont)库岸的蠕滑破坏，巴拿马运河滑坡区的片状粘土在运河边坡完工几年后突然崩塌等等。

在第1届国际岩石力学会议上，Zischinsky 用流变学模型描述了高边坡的变形，并指出岩土体的蠕变在高边坡变形中起重要作用，自此，国内外对边坡岩体流变特性的研究越来越多，成果也十分显著[1]~[12]。

五里堆滑坡区位于湖北省恩施市巴东县官渡口镇五里堆村，处于巫峡与西陵峡之间的长江北岸，与巴东新城区隔江相望，滑坡区总面积达0.476平方公里，经勘查资料证实其为多次滑移后形成的一个滑坡群。

其下伏基岩为紫红色粉砂岩夹泥岩，软硬相间，单层厚度大于0.3m 。

由于该滑坡滑床为软硬互层且软弱结构面发育，属于典型的“红层岩体”。

相关研究表明，软弱岩层是红层岩体的薄弱部位，且具有明显的蠕变特性，其发育特征和空间形态对红层岩体的稳定性有重要的控制作用[13~14]，根据对滑带土矿物成分的研究分析，可以认为该滑坡基本上是沿软弱泥岩夹层发生蠕滑变形的，因此很有必要对该软弱岩层进行蠕变试验并开展蠕变特性及其本构关系的研究，以为五里堆2号滑坡的形成机制和蠕滑机理的探索以工程学院科技论文报告会参会论文及治理工作的开展提供重要的理论和数据支持。

2.泥岩蠕变试验及其结果分析2.1蠕变试验概况和蠕变特性分析本蠕变试验仪器采用中国科学院岩土研究所研制的JQ系列剪切流变仪，精度0.001mm。

试验操作严格按照规程[15]中的要求进行，加载方法采用分级连续加载，当连续24h总变形量小于0.001mm时,开始加下一级荷载增量。

试样取于五里堆2号滑坡的典型部位，岩性为紫红色泥岩，致密块状结构，共分为5组，其垂向应力分别为表盘上读数2，4，6，8，10MPa。

取其中一组试验数据进行分析研究，如图1所示：图1 0.617MPa法向应力下蠕变试验结果Fig.1 Results of creep test for 0.617MPa norm stress从图1可以看出，每一级应力施加瞬间，泥岩试件都会出现瞬时弹性应变，其量值随着剪应力水平的增大而增大。

当剪应力较小时，试样主要以瞬时弹性应变为主，一般在1~2d内就能够达到稳定蠕变阶段。

随着剪应力的增加，初始蠕变阶段所经历的时间逐渐加大，蠕变特征逐渐明显，蠕变总应变量也逐渐增大。

当剪应力等于3.580MPa时，试样反常地经历了大约15个小时的非衰减蠕变过程，然后趋于稳定，其原因可能是泥岩试件发生了微小结构性破坏，随后经过重新调整后恢复稳定状态。

当剪应力等于7.486MPa时，泥岩试件经历了初始蠕变阶段、稳定蠕变阶段和加速蠕变阶段后最终破坏。

对于泥岩试件，其加速蠕变阶段经历的时间较为短暂，破坏前出现明显的征兆，表现为延性破坏形式。

2.2 长期剪切强度的确定岩石的蠕变性与其所承受的应力水平密切相关，当应力水平低于某一界限值时，蠕变过程中蠕变速率持续衰减直至为零，而当应力水平达到或高于该应力界限值时，就会出现稳态蠕变(0constε=≠)或加速蠕变现象，这一应力水平界限值称为岩石的长期强度[16]。

目前，确定长期强度主要有以下两种方法[9]：一种是将梯级加载的试验结果在对数坐标下绘制应力-变形等时曲线，将对数坐标中曲线上明显拐点对应的应力作为长期强度，这是因为在非衰减蠕变荷载阶段变形值会急剧增长；另一种方法是根据流动速度来判断，绘制应力与应变速率的关系曲线，曲线和应力轴的交点即为长期强度。

这里采用第一种方法来确定泥岩的长期剪切强度值，绘制的等时曲线如图2所示,由此确定0.617MPa垂直压力作用下泥岩的长期剪切强度值τs 为6.30MPa。

第工程学院科技论文报告会参会论文图2 应变-应力等时曲线Fig.2 Isochronal curves of strain-stress3 泥岩蠕变本构模型3.1衰减蠕变阶段本构模型的辨识与拟合在黏弹性体假设的前提下，利用Boltzmann迭加原理对实验数据进行处理，可以得到各剪切力水平下的剪应变-时间蠕变曲线，如图3所示。

图3不同剪应力下蠕变曲线及拟合结果Fig.3 Creep curves of different shear stresses and their fitting results在实际应用中，需要根据岩石或岩体的真实流变特征及变形性态来选择模型进行实际工程问题的分析。

故根据室内或现场试验所得的数据进行蠕变模型的辨识就显得尤为重要。

目前，蠕变模型的辨识主要有直接筛选法、后验排除法或者两种方法的综合运用[16]。

本文将采用直接筛选法，即根据应变-时间曲线特征并参照表1直接进行模型的辨识。

表1 常用蠕变模型的蠕变特性模型弹性应变应变速率黏性流动Maxwell模型有不变有Kelvin模型无递减无工程学院科技论文报告会参会论文从图2可以看出：在任何一级剪切荷载作用下，都存在着瞬时弹性应变现象；当剪应力小于长期强度时，随着时间的增长，应变速率逐渐降低，近似呈负指数的形式趋于某一渐近线；当剪应力τ0=6.835MPa>τs 时，曲线中出现了略微的黏性流动现象。

根据上述分析，并参照表1，只有Burgers 模型和西原模型符合泥岩的蠕变特性，又由于应力小于长期强度时Burgers 模型有4个参数，而西原模型退化为三参量H-K 模型，参数的增多会增加计算分析和实际应用的难度，因此，本文选用西原模型来模拟泥岩的蠕变过程，其拟合结果见图2。

从图2中可以看出，西原模型能够较好地模拟泥岩的初始蠕变阶段和衰减或非衰减的稳定蠕变阶段，但不能模拟复杂的非线性加速蠕变阶段，那是因为剪应力大于长期强度时西原模型只是增加了一个简单的线性黏塑性元件，故只能模拟线性黏性流动现象。

3.2全过程损伤模型 3.2.1 损伤模型的建立为了更好地模拟泥岩蠕变变形的全过程性态，本文将提出一种全过程蠕变损伤模型。

基于黏弹塑性介质模型理论，我们可以认为黏弹性部分与黏塑性部分为串联的关系，由此得到如下的关系式：pe σσσ== (1)pe εεε+= (2)式中，σ为总应力，σe 为黏弹性部分的应力，σp为黏塑性部分的应力；ε为总应变，εe 为黏弹性部分的应变，εp 为黏塑性部分的应变。

将式(2)对时间求导后，可得：epεεε=+ (3) 对于黏弹性流变变形部分，其算子形式的通用表达式为：e P(D)Q(D)εσ=(4)式中，P(D)，Q(D)均为黏弹性模型算子函数，e s P(D)Q(D)εσσσ=< (5)对于具有应变软化特性的岩石的黏塑性部分，则有s s 21()(1())p w t εσσσση=-≥- (6)式中，w(t)为损伤演化函数，用以表征岩石损伤的演变规律，考虑到应变量的大小往往能够表征岩体的变形特征、微裂隙的发展以及判断岩石破坏与否，因此本文由文献[17]将损伤演化方程引申为：1()(1())()cw t k w t εϕε=- (7) 其中2()1()01k ccccεεεεεϕεεε⎧≥⎪⎪=⎨⎪<⎪⎩(8) 式中，εc 为稳定蠕变与加速蠕变阶段的分界点所对应的应变量。