34岩石的流变性质解析.docx

3.4岩石的流变性质在上节中所讨论的岩石变形特性都是在加载后瞬时的变形特性，这些变形特性与时间 是无关但实际上，种类岩土工程的变形都不同程度上与时间有关例如，在中硬以下岩 石及软岩中开掘的隧道、矿山巷道等地下工程，经常出现顶板下沉、边墙挤进和底板隆起 等工程使用空间缩小现象这就是岩石流变性质的显现研究岩石流变性质，对解决岩土 工程的维护设计和长期稳定性问题有十分重要的意义岩石流变现象的组成定义如下：蠕变：应力为常量，应变随时间延长而增大的现象流变Y尸 弹性后效：加载（或卸载）后经一段时间应变才增加（或E 减小）到一定值的现象粘性流动：卸载后，部分应变永久不恢复的现象I松弛：应变为常量，应力随时间延长而减小的现象其中，蠕变现象是岩土工程中显现最明显， 对工程稳定性影响最大的流变现象， 是岩石流变理论研究中的常规内容3.4.1典型蠕变曲线特征以应变鼻为纵座标，时间t为横座标，作应变与时间的关系曲线（如图 3.24所示），该 曲线就是蠕变曲线它的形状和特性与岩石性质、加载水平等多种因素有关，各种蠕变曲线 的形状和特性不尽相同图 3.24是一条典型的蠕变曲线从曲线形态上看，可将该曲线分成三个阶段：I .AB阶段，称作为瞬态蠕变阶段（或称初始蠕变阶段） 。

加载：首先岩石特产生瞬时的弹性应变，这一应变是与时间无关的，如图中所示的 OA段当外荷载维持一定的时间后，岩石将产生一部分随时间而增大的应变， 此时的应变速率将随时间的增长远渐减小，曲线呈下凹型，并向直线状态过渡卸载：岩石的瞬时弹性应变最先恢复，如图中的 PQ段之后，随着时间的增加，其剩 余应变亦能逐渐地恢复，如图中的 QR段QR段曲线的存在，说明岩石具有随时间的增长应变逐渐恢复的特性，这一特性被称作为弹性后效n .bc阶段，被称作为稳定蠕变阶段（或称等速蠕变阶段） 加载：在这一阶段最明显的特点是应变与时间的关系近似地呈直线变化， 应变速率为一常数，该应变率与作用的外荷载的大小和介质的粘滞系数 有关卸载：出现与第一阶段卸载时一样的特性， 弹性后效仍然存在，但是这时的应变已无法全部恢复，存在着部分不能恢复的永久变形m .C点以后阶段，为非稳态蠕变（或称加速蠕变阶段） 加载：当应变达到 C点后，岩石将进入非稳态蠕变阶段这时岩石的应变速率剧烈增 加，整个曲线至上凹型， 经过短暂的时间后试件特发生破坏 C点往往被称作为蠕变极限应力，其意义类似于屈服应力I—o —图3.24岩石典型的蠕变曲线3.4.2岩石蠕变的影响因素(1) 岩性岩石本身性质是影响其蠕变性质的内在因素。

图 3.25为花岗岩等三种性质不同的岩石在室温和10MPa压应力下的蠕变曲线，由图可知，像花岗岩一类坚硬岩石，其蠕变相对很小，加荷后在很短时间内变形就趋于稳定，这种蠕变常常可忽略不计； 而像页岩、泥岩一类软弱岩石，其蠕变就很明显，变形以常速率持续增长直至破坏此外，岩石的结构构造、孔 隙串及含水性等对岩石蠕变性质也有明显的影响图3.25 常温下几种岩石的典型蠕变曲线(2) 应力对同一种岩石来说，应力大小不同，蠕变曲线的形状及各阶段的持续时间也不同图3.26为雪花石膏在不同应力下的蠕变曲线，由图可知：在低应力 (小于12. 5MPa)下，曲线不出现加速蠕变阶段；在高应力 (大于25MPa)下，则几乎不出现等速蠕变阶段，由瞬时变形 很快过渡到加速蠕变阶段，直至破坏；而在中等应力条件下，曲线呈反“ s”型，蠕变可明显分为三个阶段，但其等速阶段所持续的时间随应力增大而缩短0 20 40 60 80 IOO图3.26雪花石膏在不同压应力下的蠕变曲线(3) 温度、湿度温度和湿度对岩石蠕变有较大的影响图 3.27为人造盐岩在围压 B3=102Mpa和不同温度下的蠕变曲线由图可见，随着温度的提高，岩石的总应变与等速阶段的应变速率都明 显增加了。

另外，试验研究表明岩性不同， 岩石的总应变及蠕变速率随温度增加的幅度也不 相同湿度对岩块蠕变也有类似的影响，如 Griggs(1940)将雪花石霄浸到不同溶液中进行单轴蠕变试验，发现其总应变及蠕变速率比干燥的大，且随溶液性质不同而不同图3.27 人造盐岩在围压 o3=102Mpa和不同温度下的蠕变曲线3.4.3 岩石的流变力学模型要研究岩石的流变现象，必须首先建立其考虑时间因素的应力-应变关系称之为流变力学方程，常用的流变力学方程中最重要的是蠕变方程 它是在试验的基础上建立的， 通常包括微分模型、经验模型和积分模型等三类其中微分力学模型是由多个能够表征不同力学性 质的元件组成的物理模型推导出来的，不仅直观、 易懂，而且各种元件可以任意组合， 以满足不同特性的岩石流变性质的表征需要， 因此，在国内外应用广泛 经验模型来自于对实际工程的长期连续测试数据的统计分析下面对这两类模型进行分别介绍3.4.3.1.微分流变模型这里将岩石的变形性质分解成刚性、 弹性、塑性和粘性的四种基本形式，并用相应的基 本元件来表征若干个不同元件组合成的模型又可用于表征复杂的岩石流变性质 这里先基 本建立元件的应力一应变关系，再推出组合元件的应力一应变关系。

前三种性质和时间无关， 在以往的由粘性体流动所代表的材料变形性质和变形 （运动） 的速率有关根据试验资料，常可直观地大致判断所涉及的基本元件及组合形式 这是模型 方法的基本特点正是这种直观性和可组合性， 现今的流变模型理论在国内外获得广泛应用一、基本元件上述提及的刚性、弹性、 塑性和粘性的四种基本元件，其中前 3种都与时间无关，曾经学习过下面仅介绍粘性元件，并将四种基本元件的力学特性列于表 3.8粘性通常用牛顿体来表征，如图 3.28所示，粘缸内装牛顿液体，设有一带孔活塞该元件不论何种加载方式， 都没有瞬时应变但由于活塞中的小孔沟通两侧的液体， 在恒定力作用下，随后活塞将发生缓慢移动，其应力与应变率成正比，即（牛顿体的本构关系） ：L L d ;以 z =—— （3.12 ）d式中：一应变率；叮一粘性系数（粘性模数），单位为泊（poise ）; 1 poise=0.1N - s/ mmo式（3.12 ）为线性关系，故牛顿液体又称为线粘性液体 与其相应，理想液体的Y] = 0,而非牛顿液体的 Y]丰const o二、二元件模型常用的二元件模型有，马克斯韦尔模型、 开尔文模型、宾厄姆模型等三个， 其中前两个属于粘弹模型，后一个属于粘塑模型。

1、马克斯韦尔模型( Maxwell,1868 )马克斯韦尔模型由虎克体 (弹簧)和牛顿体(阻尼器)串联而成，简称M体，如图3.29(a)所示用等式表示为： M = H - N其中，H表示虎克体，N表示顿液体，“一"表示串联，“I "表示并联(下同)a)物理模型（1） M体的本构关系静力平衡条件：a =1 =阪 （1）变形协调条件：6 = &1十布 （2）式中，E分别表示模型的总应力和总应变； 〔，各分别表示元件1的应力和应变；2 , d分别表示元件2的应力和应变已知元件1的本构关系：，对t求导得：L - L 二-1 二1 =—=— K K又知元件2的本构关系：」d 2 二 2 ^2 = =一 dt式（2）对t求导得：* = 4 +《，并将式（3）和式（4）代入得M体本构关系:(3)(4)（2） M体蠕变方程UU a crz =— +—K(3.13 )根据蠕变的定义设： 0 = const ；由M体的原模型可知其初始条件：当 t = 0时,芝（弹簧有瞬时应变） const代入本构万程（3.22 ）得：8=— + -=一，即K(5)代入初始条件确定积分常数 C =；,得M体蠕变方程crn crn = -0t • -0K(3.14 )由式(3.14 )作蠕变曲线，如图 3.29(b)所示。

可见 M体有瞬时应变，属于线性蠕变 模型3) M体松弛方程根据松弛的定义设： 耳=so =const ；由M体的原模型可知其初始条件：当 t = 0时，设 C- = C-0将尊=So=const代入本构方程： $ =柘=；+|^=0,为线性齐次方程分离变量，—=—gdt，积分得：ln _Kt +C由初始条件确定积分常数： C=lnso,代入上式得M体松弛方程：上t Kb0e 门=boexp(-*t) (3.15)由式(3.15 )作M体松弛曲线(图 3.29 C),可见，当tT *,qt 0通常把应力下降为初始应力的e分之一(1 &0.37 )所经历的时间称为松弛时间，记为： 丁「，由松弛方程eKK 」 ―一 r得：^=Q0exp(—Ki?), e =e ，由此得M体的松弛时间：［「=一 e K(4) M体弹性后效与粘性流动弹性后效与粘性流动是固体卸载过程中变形特征， 因此，设加载至t1时刻卸载，即，当0 = const,经历t1时段，在t =t〔时刻，由蠕变方程(3.14 )得，其应变为：t^—;在此基础上卸载，即， 0,则弹簧应变瞬时恢复 生0 = 0 ,而阻 K K尼器变形则不可恢复，把不可恢复的变形称为粘性流动，其值为：e = -^-t| = const (3.16 )式(3.16)就是M体的卸载方程，卸载后应变瞬时变成常量，如图 3.29 (d),可见M M体 无弹后，但有粘流。

5) M体的流变特征总结流变特征瞬变蠕变松弛弹后粘流M体有有有无有2、开尔文模型(Kelvin体或沃伊特(Voigt)体，1890)开尔文模型由弹簧阻和尼器并联而成，简称 K体，如图3.30 (a)所示用等式表示为：M = H /N1) K体的本构关系变形协调条件:- ；1 = 2同上理，将元件1和元件2的本构关系代入上两式，得到 K的本构方程:=K ;1 K t2 = K ; |(2)(c) K体松弛曲线 (d) K体弹性后和粘性流动图3.30开尔文模型(2) K体的蠕变方程设应力为常量： 0 = const ；无瞬时应变)0 = const代入本构方程(由K体物理模型得初始条件：t = 0,3.17)得：n涉 + K@ =8=0(阻尼器先求齐次方程解：r^+K&=0, 史=—Kdt, ln8=—^t+lnC z n nK.K齐次解为 ；=Ce =Cexp(_—t)可以看出，特解为一常数 ；=二K利用初始条件确定积分常数 C =-芝，得K体蠕变方程(3.18 )；=^[1—exp(—Kt)]由此方程作K体蠕变曲线如图3.30(b)由图3.30(b)可以看出t=0时，& =0,无瞬时应变；tT叫时，6 =竺；最终最大应K变仅等于弹性元件的瞬时应变，想当于推退弹性应变的出现。

故 K体又称为延退(退滞)模型当 t =曲=—时，e = ― [1 — exp(—— — )] = —— 。