北京工业大学学报
2018年3月JOURNAL OF BEIJING UNIVERSITY OF TECHNOLOGY
Vol.44 No.3 Mar. 2018
变强度参数蠕变模型研究
何利军1,2,3，张涛\\胡颖萍1
(1.南昌航空大学土木建筑学院，南昌330063; 2.广西防灾减灾与工程安全重点实验室，南
3.广西大学工程防灾与结构安全教育部重点实验室，南宁530004)
宁
530004;
摘要：为了解决非线性蠕变模型的参数取值问题，需要把握参数变化规律，提出了 FLAC3D的Cvisc模型中关于 强度参数随蠕变时间变化的FISH语言编写形式，探讨了黏聚力与内摩擦角不同的变化情形对加速蠕变变形的影 响规律，得出了湛江黏土在髙偏应力条件下的蠕变变形过程中，强度参数以黏聚力减小或先有所增大后减小为主 的初步结论.
关键词：加速蠕变；变强度参数；湛江黏土中图分类号：TU433
文献标志码：
A 文章编号：0254 -0037(2018)03 -0327 -09
doi： 10.11936/bjutxb2017070026
Study on Creep Model of Variable Strength Parameters
HE Lijun1，2,3，ZHANG Tao1，HU Yingping1
(1. College of Civil Engineering and Architecture，Nanchang Hangkong University，Nanchang 330063，China；
2. Key Laboratory of Disaster Prevention and Engineering Safety of Guangxi，Nanning 530004，China;
3. Key Laboratory of Disaster Prevention and Structural Safety of China Ministry of Education ，Guangxi University，
Nanning 530004，China)
Abstract ： In order to solve the problem of parameter selection of nonlinear creep model，it is necessary
to grasp the change of parameters. The FISH language of that the variation of intensity parameters with creep time in Cvisc model of FLAC3D was proposed. The influence of different cohesion and internal friction angle on thedeformation of accelerated creep was discussed， and a preliminary conclusion was drawn that the strength parameters of Zhanjiang clay were reduced or first increased and then decreased under creep of the condition of high deviatoric stress. It shows that the research contributes to exploring the law of accelerating creep，and has important theoretical significance and practical reference value for the correct simulation of creep phenomenon.
Key words： accelerated creep; variable strength parameters; Zhanjiang clay
对蠕变模型的研究由来已久，近5年来,赵延林 等[1]提出岩石黏弹塑性应变分离的蠕变试验方法 和数据处理技巧.在FLAC3D中实现BNMC蠕变损 伤本构模型的二次开发，提出该模型参数的辩识方 法.王军保等[2]通过将表征应力水平对岩石非线性
收稿日期：2017-07-17
基金项目：南昌航空大学博士启动资助项目（EA2012111000);广西防灾减灾与工程安全重点实验室开放课题资助项目
蠕变特性影响的函数引入到常规分数阶黏滞体的本 构关系中，提出一种改进的分数阶黏滞体，能够描述 岩石加速蠕变变形.杨振伟等[3]基于三维颗粒流程 序（PFC3D)，米用Burgers流变模型,分析了伯格斯 模型中的参数对瞬时强度特性和流变特性的影响.
(2016ZDK013);广西大学工程防灾与结构安全教育部重点实验室开放课题资助项目（2016ZDK013)作者简介：何利军（1977—），男，讲师，主要从事土力学基本理论与数值模拟的研究，E-mail:70457@nchu.edu.cn
328北京工业大学学报2018 年
杨红伟等[4]以细粒砂岩为研究对象,对三轴压缩条 件下岩石孔隙水压力分级加载蠕变试验进行了蠕变 特性及模型研究.王军保等[5]基于损伤力学理论, 通过引入Kachanov蠕变损伤演化规律，构建了可反 映岩石加载瞬时变形和加速蠕变的考虑时效损伤的 弹性体元件，建立了一个新的2元件非线性蠕变损 伤模型及其蠕变方程.康永刚等[6]提出用非牛顿黏 壶代替Kelvin模型中的牛顿黏壶，串联用一种非定 常、非牛顿黏壶和塑性体并联为黏塑性体，黏塑性体 关闭时,所得蠕变函数可以描述衰减蠕变和稳态蠕 与本文有直接联系的论文成果有:杨圣奇等[15] 基于在岩石剪切流变仪上得到的泥板岩剪切流变试 验结果,采用库仑强度准则获得了岩石长期强度参 数（长期黏聚力和长期内摩擦系数），与瞬时强度参 数进行对比，结果表明,岩石黏聚力对时间的敏感程 度高于内摩擦系数对时间的敏感程度，提出了一个 考虑黏聚力和内摩擦系数随时间变化的塑性元件, 将其与线性黏性元件并联起来，得到一个新的考虑 黏聚力和内摩擦系数的非线性黏塑性体（CF-
NVPB),同时，将CF-NVPB模型与三元件线性黏弹
变.黏塑性体开启时,蠕变函数可以描述加速蠕变 变形.蒋海飞等[7]采用幂函数、对数函数混合方程 对加速蠕变段进行拟合，再通过类比的方法提出一 个新的非线性黏性元件,并将其与塑性体并联，得到 一个可以反映岩石加速蠕变特性的非线性黏塑性模 型，将该模型与Burgers模型串联，构建一个新的六 元件非线性黏弹塑性蠕变模型.此外，还提出过一 个能反映不同孔隙水压作用下加速蠕变特征的二元 件黏塑性蠕变模型[8]，将该模型与修正的广义 Kelvin模型串联组成一个新的改进西原模型.范翔 宇等[9]利用饱依丁-汤姆逊体与宾汉姆体,建立储 气层煤岩蠕变模型.刘浪等[10]分析深部岩石在饱 水情况下的流变规律，并基于统一流变力学模型，对 深部饱水岩石的流变模型进行辨识,确定流变模型 的相关参数.乔丽苹等[11]提出采用不可恢复应变 作为内变量来描述复杂条件下岩石蠕变性质，并据 此提出了岩石的内变量蠕变模型.宋勇军等[12]将 含分数阶导数的软体元件与虎克体串联，引入能反 映应力水平和时间影响的损伤变量,提出一种四元 件非线性蠕变损伤模型.并借鉴经典元件模型的建 模思路[13]，将岩石的初始屈服强度作为蠕变硬化的 应力阈值，岩石的长期强度作为损伤软化的应力阈 值，引入能反映岩石硬化效应的硬化函数和损伤效 应的损伤变量,建立能够全面反映蠕变机制的岩石 非线性蠕变模型.高赛红等[14]基于蠕变过程中的 硬化-损伤机制，对高应力条件下衰减蠕变阶段的 黏滞系数(硬化黏滞系数）进行修正,结合定常蠕变 阶段的损伤蠕变机制，提出高应力条件下岩石损伤 效应的累积-扩散机理，并据此引入以累积损伤蠕 变量为判定准则的加速蠕变触发模型（非线性蠕变 体)描述岩石的不稳定蠕变特征.在此基础上，引入 瞬时塑性元件,与虎克体和黏弹塑性体串联,建立了 一个能够完整描述岩石蠕变全过程的非线性黏弹塑 性蠕变模型.
性流变模型串联起来,建立了一个新的考虑黏聚力 和内摩擦系数的岩石非线性黏弹塑性流变模型，并 推导了该非线性流变模型的蠕变方程与松弛方程.
佘成学等[16]在西原模型的基础上，研究建立了一种 新的岩石非线性蠕变模型.首先,在分析岩石宏、细 观蠕变过程的基础上，提出了蠕变参数变化率概念; 以此为基础，统一建立了等速蠕变和加速蠕变阶段 的岩石蠕变参数变化公式.然后，利用该非线性蠕 变参数以及根据岩石的屈服软化过程，以西原模型 为基础，建立了岩石的非线性弹-黏弹-黏塑性蠕变 模型.最后，将单轴压缩蠕变试验结果与数值计算 结果对比，验证了模型对于描述岩石全过程蠕变变 形破坏的正确性.
通常认为岩石强度伴随着结构破损，即塑性剪 切应变的增加而降低，与塑性剪应变或时间有关，相 对于岩石而言,特殊土的力学性质更为复杂、强度敏 感性更强，强度伴随时间的变化规律更复杂,即使没 有塑性剪切应变发生，原状土受到扰动，强度也会发 生先减后增的变化，探讨特殊土的关于变参数蠕变 模型的具体形式及应用方式的论文成果有：
孔令伟等[17]基于湛江黏土的蠕变变形在低偏 应力下较小，其蠕变等时曲线没有明显的屈服特征 而近似为直线，表现为线性蠕变性状;当偏应力达到 某一临界值后,在短时间内发展到加速蠕变阶段并 发生破坏.从湛江黏土蠕变各阶段的变形机制出 发,通过引入服从Mohr-Coulomb塑性流动规律的变 参数塑性元件,将该塑性元件与广义Kelvin模型串 联，提出描述蠕变全过程的黏弹塑性模型.同时将 该模型在FLAC3D中进行二次开发，与FLAC3D中 自带的Cvisc模型进行对比验证，结果表明，基于变 参数塑性元件的黏弹塑性模型能较好地描述湛江黏 土各个蠕变阶段的变形性状.
为了尽可能准确描述湛江黏土蠕变全过程，并 试图从土工材料流变学角度能够给出合理解释[18],
第3期何利军，等：变强度参数蠕变模型研究329
需要引入塑性元件，塑性元件服从的强度准则有很
多种，这里还是采用3描述土体常用的
的棘手问题，尤其在地基处理、边坡施工和围岩开挖 与支护等工程中，岩土长时间在恒定外荷载作用下， 产生的持续变形会破坏其上建筑（或构筑物）的稳 定性和安全性;尤其是加速蠕变的出现，给正在建设 或已经完工的工程带来巨大的破坏和损失.而由摩 尔-库伦强度准则，在外荷载不变的情况下，强度参 数的大小决定土体内部是否存在达到强度屈服准则 的区域，即塑性区，而一旦产生塑形区便意味着该部 分区域的变形量将可能不断增加.对于具有蠕变特 性的土体而言，即会出现加速蠕变现象.所以强度 Mohr-Coulomb塑性元件.传统塑性元件又称为摩阻
FLACD
元件，其是一个线性元件.当施加的应力大于该元 件的摩阻力时，元件将随时间产生持续增加的变形， 一旦外力卸去，已发生的变形也不会恢复，成为永久 应变.然而，由莫尔-库伦强度理论，土体的塑性发
生通常与土体的黏聚力c和内摩擦角渍两个强度参 数有关.本文提出变强度参数蠕变模型应用于湛江 黏土，探讨强度参数c和渍的变化对土体加速蠕变 阶段的影响，试图得到一些基本规律为实际工程建 设提供一些建议和参考.
本文的研究对于完善非线性蠕变模型理论，促 进与推动关于能描述加速蠕变阶段的蠕变全过程理 论模型的构建，具有重要的意义和价值.
1 Cvisc模型
Cvisc模型是由传统经典Burger模型和一■个塑
性元件串联而成，见图1，使得该本构模型能够描述 材料的弹性、黏性、黏弹及塑性变形.当施加的应力 滓 伦强度准则. 图1 Cvisc模型 Fig. 1 Cvisc model Cvisc模型一'维本构关系为 着(^)= ’E滓 滓 滓 -%、^ M + 浊M 1 + EK (1-e \"K)， 滓 < 滓 s滓 + 滓浊 1 + 滓、 E M E K (1 - e- K^Kt ) + 着v( t^)，滓〉滓s、乜 M(1) 式中sv(t滓）表示塑性体变形量是关于时间和应力 的函数，其他参数同流变学文献一致[19]，这里不再 赘述. 2 变强度参数 2.1变强度参数的引入 在工程建设中，岩土的蠕变问题是个需要重视 参数对土体加速蠕变阶段出现与否存在一定的 影响. 另外大量的土工试验及研究表明[20-21]: 土体在 受荷产生变形直至破坏的过程中，其强度参数并非 一成不变，而是随着位移的发生，强度有增有减，另 外土体的破坏过程区别于岩石，其具有一个明显的 渐变过程，属于塑性破坏，而土体受荷产生的位移量 是关于时间的函数.这种受荷呈现出强度随时间变 化的特征，表现出土体作为一种三相材料体，即具有 颗粒滑移的一'般特性又具有一'定的自身结构和黏结 能力的基本特征. 本文根据对湛江黏土三轴试验获得的瞬时强度 指标和长期强度指标来看，蠕变引起的强度衰减主 要表现为黏聚力c值的明显降低，而摩擦角渍变化 不大[22]. 基于上述2个指标在时间变化上表现出的差 异，引入变强度参数模型，考虑采用FLAC3D所自带 的 FISH语言功能定义一个函数，实现土体强度参 数c和渍随蠕变时间发生变化，试图通过该模型来 描述160 kPa偏应力下湛江黏土的加速蠕变阶段. 2.2变强度参数的建立 文中变强度参数模型是在前述 Cvisc模型的基 础上建立起来，通过FLAC3D内置FISH语言定义函 数，实现强度参数c和渍随着蠕变时间不断发生变 化，此外，实际的土体在受荷后，即使强度产生衰减 导致发生破坏，其依然具有一定的残余强度，也即其 强度参数c和渍不会为零，所以在建立强度参数随 时间变化的过程中，要么给定2个强度指标降低的 一个最低限值;要么确定蠕变时间，通过定义它们的 降低速率来控制最低值FISH . 下面是通过语言定义一个函数，在函数 内部定义2个强度参数值随蠕变时间不断减小，具 体实现代码如下defcvscs_gb ： _ 330北京工业大学学报2018 年 whilestepping p_z = zone_head loop while p_z # null p_coh = y_coh - j_xl * crtime p_fri = z_frw - f_js * crtime z_prop ( p_z, coh') = p_coh z_prop ( p_z, fri') = p_fri p_z = z_next ( p_z ) endloop end sety_coh = 38e3 ;初始黏聚力为 38kPa set j_xl = 2e4 setz_frw = 25 ;初始内摩擦角为25。 set f_js = 5 cvs_cs_gb 上述代码2个强度参数随蠕变时间是不断降低 的，用户可以根据实际情况添加if. . else..条件语句 来控制它们降低的最低限值.然后将代码插入到三 轴蠕变命令流文本中，通过solveagexx来设置2个 不同的求解时间点，分别在软件界面通过点击 Plootitems-Add-Block-property 来选择查看黏聚力和 内摩擦角是否发生变化，上述代码经过验证能够实 现2个强度参数随蠕变时间发生改变，见图2. 3 c和渍对加速蠕变的影响 3.1 c、渍变化的9种情形 采用的模型为如图2所示的圆柱体，高为1个 单位，直径为0. 5个单位，主要分析强度参数对湛江 黏土加速蠕变阶段出现的影响，虽然影响土体加速 蠕变阶段出现的因素很多，如荷载水平、湿度、颗粒 级配和结构等，但在这里暂不考虑.为了通过数值 模拟，准确分析2个强度参数对加速蠕变出现的影 响程度，需要将强度参数变化分成以下9种情形进 行简要分析： ) 黏聚力和内摩擦角都不变. ) 黏聚力减小，内摩擦角不变.) 黏聚力增大，内摩擦角不变.) 内摩擦角减小，黏聚力不变.) 内摩擦角增大，黏聚力不变.) 黏聚力先增后减，内摩擦角不变.) 黏聚力先减后增，内摩擦角不变.) 内摩擦角先增后减，黏聚力不变.) 内摩擦角先减后增，黏聚力不变. 显然上述9种情形不能完全包含2个强度参数 指标可能出现的所有变化范围.此外，还有黏聚力 和摩擦角可能呈拋物线形式同时变化，发生先增后 减、先减后增或呈指数函数形式变化等等，这些变化 形式应当结合工程实际进行具体分析，现仅讨论以 上9种情形所产生的影响，对于其他变化情况可参 照同样的分析方法.3.2数值计算参数 为了便于分析，将数值结果与试验数据进行比 较，选用湛江黏土在160 kPa偏应力下初始蠕变阶 段和稳定蠕变阶段试验数据，采用Burger三维本构 关系拟合相关参数，拟合方法与前述相同，强度参数 已由试验测定得到，初始黏聚力可选取为38 kPa，内 摩擦角为25。[17]，拟合出相关参数见表1. 1 234567 图2 2000时步强度参数c和渍 Fig. 2 Intensity of the parameters c and 渍 with 2 000 times 表1 160 kPa下Burger/Cvisc模型相关参数 Table 1 Relevant parameters of the Burger/Cvisc model under 160 kPa 滓砖 /kPa K/ kPa 39 018 GM/kPa 18 008 ^M/(kPa*h-1) 1 858 GK/kPa 12 963 ”K/(kPa.h-1) 227 R0.997 2 160 3.3 9种情形分析比较前2阶段试验数据拟合出的参数，在强度参数不变 的情况下，整个蠕变曲线是趋向于稳定，没有发生加 速蠕变现象，且2条曲线在初始蠕变阶段和稳定蠕 变阶段基本一致. 1)强度参数不发生变化，单纯Cvisc模型进行 数值计算，其结果与三轴蠕变试验数据对比见图3. 从图3可以看出：当依据在偏应力160 kPa下 第3期何利军，等：变强度参数蠕变模型研究331 10 8 6 4 2 C变阶段向加速蠕变阶段转化的明显转折点，那么从 图中可以看出随着黏聚力降低速率进一步增大，该 “突变点”不断左移，表明黏聚力降低，保持其他因 素不变，土体在相同条件下更容易出现加速蠕变达 到破坏.最后，从图中显而易见，在降低速率为70 X 103时，该数值曲线与三轴试验曲线基本一致.②黏聚力随蠕变时间不断升高，内摩擦角不 变；由于c = 38 kPa时不会出现加速蠕变，此时继续 增加黏聚力，显得毫无意义，所以为了表明黏聚力的 增加对加速蠕变的影响，这里将黏聚力初始值设为 %/_保挡叵捃0 0.1 0.2 0.3 t/h 0.4 0.5 0.6 图3 Cvisc模型数值计算结果与三轴蠕变试验数据对比 Fig. 3 Comparison of numerical results of Cvisc model with test data of triaxial creep 2)强度参数按情形2、3变化，采用一次函数 y = + b进行描述，a为黏聚力变化率，kPa/h,b为 初始黏聚力38 kPa;t为蠕变时间，蠕变步长设为5 X 10-5;虽然蠕变是一个历时很长的过程，但试验状态 下的加速蠕变历时很短，湛江黏土的加速蠕变时间 只能记录到0. 55 h，与此对应，计算模拟的总蠕变时 间定为0. 5h，内摩擦角25毅.情形2、3分别是黏聚 力下降和增加，而为了更直观表明黏聚力变化对加 速蠕变阶段出现时间的影响，将上述2种情况依据 黏聚力的变化速率快慢又各分为4小类，具体如下: ①黏聚力随蠕变时间不断降低，摩擦角不变; 下降情况见图4,所得数值计算结果见图5. 12 %180/_ 傲6挡 4叵 捃 r2 v-0 0.1 0.2 03 t/h 0.4 0.5 0.6 图5 黏聚力递减蠕变曲线 Fig. 5 Creep curves with decreasing cohesion 由图5可以看出，在变强度参数模型中，随着给 定黏聚力降低幅度的增大，试样会逐渐出现加速蠕 变现象，且加速蠕变量随着黏聚力的降低而不断增 大.另外，如果将稳定蠕变阶段的延长线与加速蠕 变阶段切线的交点称之为“突变点冶，意即由稳定蠕 8 kPa.升高情况如图6,所得数值计算结果见图7. : c =8xl03c'=8xl03+30xl03f 1 6c+4〇xl〇3f ，8xl03+50xl03f ±」% /_12傲± 挡-8叵据 -4 0 一 0.1 0.2 0.3 t/h 0.4 0.5 0.6 图7 黏聚力递增蠕变曲线 Fig. 7 Creep curves with increased cohesion 从图7可知：初始黏滞系数为8 kPa，4条模拟 曲线都是“S”型曲线，而不同于“反倒S”型的三轴 试验数据的，此时蠕变曲线只包含初始蠕变和稳定 蠕变2个阶段；曲线的竖直位置次序表明随着黏聚 力增长速率的提高，轴向应变量逐渐减小，稳定蠕变 的“稳定点”，也就是曲线由斜向上到水平的拐点逐 渐左移. 3)强度参数按情形4、5变化，同样分为2种情 况，采用一次函数y = ct + d;d为初始内摩擦角，c为 降低或升高速率，t为蠕变时间. ①内摩擦角随蠕变时间不断降低，给定初始值 为25毅，黏聚力不变;变化情况见图8，所得数值计算 结果见图9. 由图9显然可知：内摩擦角对土体的加速蠕变 332北京工业大学学报2018 年 图9 内摩擦角递减蠕变曲线 Fig. 9 Creep curves with decreasing internal friction angle 阶段存在影响.影响规律为随着内摩擦角降幅的提 高，加速蠕变阶段的出现越发明显，随着内摩擦角降 幅的增加，蠕变曲线的“突变点”不断左移，意味着 加速蠕变阶段的提早出现. ②内摩擦角随蠕变时间不断增加，黏聚力不 变;同样为了表明内摩擦角增加对加速蠕变阶段的 影响，给定内摩擦角初始值为5毅.具体变化情况见 图10,所得数值计算结果见图11. vs. time curves 图11表明：在内摩擦角增幅最低时，c = 10,得 到最早加速蠕变曲线；当增幅达到c =20时，得到稳 定蠕变曲线，未出现加速蠕变阶段，所以随着内摩擦 角增幅逐渐增大，蠕变曲线会逐渐变缓，导致对应的 变形量明显降低，加速蠕变阶段不再显现.4)黏聚力按情形6变化，内摩擦角不变.初始 黏聚力设为5 kPa，按照先增后减变化规律，分为4 图1 1内摩擦角递增蠕变曲线 Fig.11 Creep curves with increasesd internal friction angle 种变化幅度，采用二次函数（拋物线形式）来描述黏 聚力的变化情况，具体见图12,依据变化情况所得 数值计算曲线见图13. 图1 2黏聚力先增后减曲线 Fig. 12 First increaseing and then decreasing cohesion vs. time curves+/5=5+他 5 o 4o3 o 2 o o 1 0 0.1 0.2 0.3 t/h 0.4 0,5 0.6 图1 3黏聚力先增后减蠕变曲线 Fig. 13 Creep curves with first increasing and then decreasing cohesion 由图13可以看出：在黏聚力先增后减变化幅度 最小时，出现最大加速蠕变，随着黏聚力变化幅度的 增大，蠕变曲线的“突变点”逐渐右移，加速蠕变阶 段出现时间推后且对应的蠕变量明显降低. 5)黏聚力按情形7变化，内摩擦角不变.根据 黏聚力先减后增的变化幅度不同，分为4种情况，同 样采用二次函数(拋物线形式）来进行描述，具体见 图14,依据变化情况所得数值计算曲线如图15. 在图15中，曲线为“S”型曲线，黏聚力先减后 第3期何利军，等：变强度参数蠕变模型研究 333 图1 4黏聚力先减后增曲线 Fig. 14 First decreasing and then increasing cohesion vs. time curves 图1 5黏聚力先减后增蠕变曲线 Fig. 15 Creep curves with first decreasing and then increasing cohesion 增变化幅度较大时，轴向应变量较大，随着幅度的降 低，“稳定点”逐渐左移. 6)内摩擦角按情形8变化，黏聚力不变.内摩 擦角根据先增后减的变化幅度不同，分为4种情况， 也是采用二次函数(拋物线形式）来进行描述，具体 见图16,依据变化情况所得数值计算曲线见图17. 图1 6内摩擦角先增后减曲线 Fig. 16 First increasing and then decreasing internal friciton angles vs. time curves 在图17中，4条模拟曲线都是“反倒S”型，“突 变点”的水平位置在前3条上并没有明显变化，轴 向应变最大的第4条更多的是受内摩擦角偏小的 30 # ^^=-320^+160^+5 /25—^=-240^+120f+5 _ z3=-160^+80^+5 20傲z4=-8〇P+4(^+5 挡15+ 三轴试验数据^ 叵1 0捃 1 0 0.1 0.2 03 t/h 0.4 0.5 0.6 图1 7内摩擦角先增后减蠕变曲线 Fig. 17 Creep curves with first decreasing and then increasing internal friction angle 影响. 7)内摩擦角按情形9变化，黏聚力不变.内 摩擦角根据先减后增的变化幅度不同，分为4种 情况，同样也是采用二次函数（拋物线形式）来进 行描述，具体见图18,依据变化情况所得数值计算 曲线如图19: c/2 0«1 5 «1 -»E: 0 5 0 0.1 0.2 t/h 03 0.4 0.5 图1 8内摩擦角先减后增曲线 Fig. 18 First decreasing and then increasing internal friction angle vs. time curves 4 o —zz6=160f2-80?+25 ：s=m2-40t+25 3 o z=24Ot2-U0t+25o z=320t2-l60t+25 2+三轴试验数据 o 1 0 0.1 0.2 03 t/h 0.4 0.5 0.6 图1 9内摩擦角先减后增蠕变曲线 Fig. 19 Creep curves with first decreasing and then increasing internal friction angle 4条模拟曲线都是“S”型，与试验曲线的“反倒 S”型不同.同样的，随先减后增幅度的增大，轴向应 变量增大. 综合对上述5种情形的分析，可以得出如下 结论： 1) 黏聚力和内摩擦角对土体的加速蠕变阶 存在显著的影响. 2) 黏聚力或内摩擦角增加，会使得土体难以 入加速蠕变阶段；反之，土体会提前进入加速蠕变 334北京工业大学学报2018 年 240-248. YANG Z W，JING A B，ZHOU Y，et al. Parametric 3) 从模拟曲线与三轴试验曲线的变化趋势的 analysis of Burgers model and creep properties of rock with 对比来看，黏聚力和内摩擦角按线性减小、黏聚力按 particle flow code [J]. Rock and Soil Mechanics，2015, 抛物线先增后减与内摩擦角先增后减时，模拟曲线 36(1)： 240-248. (in Chinese) 与试验曲线变化规律能够保持一致. [4] 杨红伟，许江，聂闻，等.渗流水压力分级加载岩石蠕 4) 模拟曲线与三轴试验曲线的变化趋势最不 变模型[J].宜宾学院学报，2015, 15(6)： 1-5. 一致的是:黏聚力和内摩擦角按线性减增加、黏聚力 YANG H W，XU J，NIE W，et al. Study on the creep 按抛物线先减后增与内摩擦角先减后增时,模拟曲 model of rocks under step loading of seepage pressure[J]. 线以“S”型曲线为主,而试验曲线是“反倒S”型.Journal of Yibin University， 2015，15 (6)： 1-5. (in Chinese) 4结论 [5] 王军保，刘新荣，邵珠山，等.岩石非线性蠕变损伤模 阶段. 1) 本文从采用Cvisc模型无法描述湛江黏土加 速蠕变阶段，且依据对湛江黏土的长期和瞬时强度 指标的差异来看,引人变强度参数模型，讨论强度参 数变化对土体加速蠕变阶段的影响.然后将强度参 数变化分为9种情形，采用FISH语言定义变强度参 数模型进行数值计算,将模拟曲线与三轴试验曲线 的变化趋势进行对比，结合试验研究得出的湛江黏 土内摩擦角变化不大的结论，可以得出湛江黏土在 能导致加速蠕变的高偏应力条件下的强度参数变化 规律：以黏聚力逐渐减小或先有所增大后减小为主 要变化趋势. 2) 综合而言，本文是通过数值模拟探讨强度参 数的变化规律对于高偏应力条件下蠕变变形的影 响，所做研究工作对于探索加速蠕变的内部机理，具 有较为重要的意义和价值. 参考文献：[1]赵 延林，唐劲舟，付成成，等.岩石黏弹塑性应变分离 的流变试验与蠕变损伤模型[J].岩石力学与工程学 报，2016, 35(7): 1297-1308.ZHAO Y L，TANG J Z，FU C C，et al. Rheological test of separation between viscoelastic-plastic strains and creep damage model [ J]. Chinese Journal of Rock Mechanics and Engineering， 2016， 35 ( 7 )： 1297-1308. (inChinese)[2] 王军保，刘新荣，王铁行.基于改进分数阶黏滞体的岩 石非线性蠕变模型[J].中南大学学报（自然科学版）， 2015, 46(4) ： 1461-1467.WANG J B，LIU X R，WANG T X. A nonlinear creep model for rocks based on modified fractional viscous body [J]. Journal of Central South University ( Science and Technology)，2015，46(4) ： 1461-1467. (in Chinese)[3]杨 振伟，金爱兵，周喻，等.伯格斯模型参数调试与岩 石蠕变特性颗粒流分析[J].岩土力学，2015,36(1): 型研究[J].现代隧道技术， 2014, 51(3)： 79-84. WANG J B，LIU X R，SHAO Z S，et al. A study of the nonlinear creep damage model for rocks [ J ] . Modern Tunneling Technology，2014，51 ( 3 )： 79-84. ( in Chinese) [6] 康永刚，张秀娥.一种改进的岩石蠕变本构模型[J] 岩土力学，2014,35(4)：1049-1055.KANG Y G， ZHANG X E. An improved constitutive model for rock creep [ J]. Rock and Soil Mechanics， 2014, 35(4)： 1049-1055. (in Chinese) [7] 蒋海飞，刘东燕，赵宝云，等.高围压高水压条件下岩 石非线性蠕变本构模型[J].采矿与安全工程学报， 2014, 31(2) ： 284-291.JIANG HF，LIUDY，ZHAO BY，et al. Nonlinear creep constitutive model of rock under high confining pressure and high water pressure [J]. Journal of Mining & Safety Engineering，2014，31(2) ： 284-291. (in Chinese) [8] 蒋海飞，刘东燕，黄伟，等.高围压下不同孔隙水压作 用时岩石蠕变特性及改进西原模型[J].岩土工程学 报，2014, 36(3) ： 443^51. JIANG H F，LIU D Y，HUANG W，et al. Creepproperties of rock under high confining pressure and different pore water pressures and a modified Nishihara model [J]. Chinese Journal of Geotechnical Engineering， 2014, 36(3) ： 443^51. (in Chinese) [9] 范翔宇，张千贵，艾巍，等.煤岩储气层岩石蠕变特性 与本构模型研究[J].岩石力学与工程学报，2013,32 (增刊 2)：3732-3739. FAN X Y，ZHANG Q G，AI W，et al. Research on creep property and constitutive model of rock in coal bed gas reservoir [ J]. Chinese Journal of Rock Mechanics and Engineering， 2013，32 ( Suppl 2 ) ： 3732-3739. (in Chinese)[10] LIU L，WANG G M，CHEN J H，et al. Creep experiment and rheological model of deep saturated rock [ J ] . Transactions of Nonferrous Metals Society of第3期何利军，等：变强度参数蠕变模型研究 报（工学版） 335 China，2013, 23(2) ： 478-483.[11] 乔丽苹，王者超，李术才，等.岩石内变量蠕变模型 研究 [J].岩土力学，2012, 33 ( 12): 3529-3537, 3603. QIAO L P，WANG Z C，LI S C，et al. An internal- variable creep model for rocks [ J]. Rock and Soil Mechanics，2012，33(12)： 3529-3537，3603. (in Chinese) [12] 宋勇军，雷胜友.基于分数阶微积分的岩石非线性蠕 变损伤力学模型[J].地下空间与工程学报，2013,9 ，2009 , 42(1)： 25-28. SHE C X， CUI X. Research on nonlinear creep model of rock [ J ] . Engineering Journal of Wuhan University， 2009, 42(1)：25-28. (in Chinese) [17] 孔令伟，何利军，张先伟.湛江黏土的蠕变模型与变 参数塑性元件[J].岩土力学，2012,33(8)： 22412246. KONG L W， HE L J， ZHANG X W. Creep model of Zhanjiang clay and plastic components with variable parameters [ J]. Rock and Soil Mechanics， 2012， 33 (1)： 91-95, 112. SONG Y J，LEI S Y. Mechanical model of rock nonlinear creep damage based on fractional calculus [J]. Chinese Journal of Underground Space and Engineering，2013，9 (1)：91-95, 112. (in Chinese) [13] 宋勇军，雷胜友，刘向科.基于硬化和损伤效应的岩 石非线性蠕变模型[J].煤炭学报，2012,37(增刊 2) ： 287-292. SONG Y J，LEI S Y，LIU X K. Non-linear rock creep model based on hardening and damage effect[ J] . Journal of China Coal Society，2012, 37 (Suppl 2)： 287-292.( in Chinese) [14] 高赛红，曹平，汪胜莲，等.改进的岩石非线性黏弹 塑性蠕变模型及其硬化黏滞系数的修正[J].煤炭学 报，2012,37(6)：936-943. GAO S H，CAO P，WANG S L，et al. Improved nonlinear viscoelasto-plastic rheological model of rock and its correction of hardening coefficient of viscosity [ J] . Journal of China Coal Society，2012，37(6) ： 936-943.( in Chinese) [15] 杨圣奇，朱运华，于世海.考虑黏聚力与内摩擦系数 的岩石黏弹塑性流变模型[J].河海大学学报（自然 科学版），2007,35(3)：291-297.YANG S Q， ZHU Y H， YU S H. Nonlinear visio-elasto- plastic rheological model of rock with consideration of cohesion and internal friction coefficient [J]. Journal of Hohai University (Natural Sciences )， 2007， 35 (3)： 291-297. (in Chinese)[16] 佘成学，崔旋.岩石非线性蠕变模型[J].武汉大学学 (8)：2241-2246. (in Chinese)[18]黄 明奎.公路工程材料流变学[M].成都：西南交通 大学出版社，2010. [19] 孙钧.岩土材料流变及其工程应用[M].北京：中国 建筑工业出版社，1999. [20] 刘志伟.土体强度参数对边坡隧道影响数值模拟分析 [J].水利与建筑工程学报，2014, 12(4)： 176-180. LIU Z W. Numerical stimulation analysis for the effects of soil mass strength parameters on the tunnels with side slopes [J]. Journal of Water Resources and Architectural Engineering，2014，12(4)： 176-180. (in Chinese) [21] 周春梅，赵子鹏，鲁阳.含水量对滑带土强度变形参 数及滑坡稳定性的影响[J].防灾减灾工程学报， 2016, 36(2) ： 213-219. ZHOU C M，ZHAO Z P，LU Y. The influence of water content on strength and deformation parameters of sliding zone and slope stability [ J ]. Journal of Disaster Prevention and Mitigation Engineering， 2016，36 (2)： 213-219. (in Chinese) [22] 孔令伟，臧濛，郭爱国，等.湛江强结构性黏土强度 特性的应力路径效应[J].岩土力学，2015,36(A1)： 19-24. KONG L W，ZANG M，GUO AG，et al. Effect of stress path on strength properties of Zhanjiang strong structuredclay[J]. Rock and Soil Mechanics，2015，36(A1)： 1924. (in Chinese) (责任编辑郑筱梅）
因篇幅问题不能全部显示，请点此查看更多更全内容