第2章 岩土流变力学

1、第二章 岩土流变力学2-1 2-1 岩石工程中的流变问题岩石工程中的流变问题 岩体的变形随时间增长而变化岩体的变形随时间增长而变化地下硐室的开挖地下硐室的开挖岩基地基岩基地基岩石边坡岩石边坡一、流变的概念一、流变的概念岩石的流变性岩石的流变性: :岩石应力应变关系随时间而变化的性质。岩石应力应变关系随时间而变化的性质。外部条件不变时，应力或变形随时间而缓慢变化。外部条件不变时，应力或变形随时间而缓慢变化。流变性（粘性）流变性（粘性）蠕变蠕变松弛松弛弹性后效弹性后效长期强度长期强度2-2 岩石流变力学属性岩石流变力学属性蠕变蠕变: :蠕变是当应力不变时，变形随时间增加而增长的现象蠕变是当应力不变

2、时，变形随时间增加而增长的现象松弛松弛: :松弛是当应变不变时，应力随时间增加而减小的现象松弛是当应变不变时，应力随时间增加而减小的现象弹性后效弹性后效: :弹性后效是加载或卸载时，弹性应变滞后于应力弹性后效是加载或卸载时，弹性应变滞后于应力的现象的现象长期强度长期强度: :在长期载荷持续作用下岩石的强度在长期载荷持续作用下岩石的强度 t t的强度的强度长期强度方法方法1 1方法方法2 2二、岩石的蠕变性能二、岩石的蠕变性能1 1、岩石的蠕变特性、岩石的蠕变特性 通常用蠕变曲线（通常用蠕变曲线（-t-t曲线）表示岩石的蠕变特性。曲线）表示岩石的蠕变特性。（1 1）稳定蠕变稳定蠕变：岩石在较小的

3、恒定力作用下，变形随时：岩石在较小的恒定力作用下，变形随时间增加到一定程度后就趋于稳定，不再随时间增加而变化，间增加到一定程度后就趋于稳定，不再随时间增加而变化，应变保持为一个常数。稳定蠕变一般不会导致岩体整体失稳。应变保持为一个常数。稳定蠕变一般不会导致岩体整体失稳。（2 2）非稳定蠕变非稳定蠕变：岩石承受的恒定荷载较大，当岩石应：岩石承受的恒定荷载较大，当岩石应力超过某一临界值时，变形随时间增加而增大，其变形速率力超过某一临界值时，变形随时间增加而增大，其变形速率逐渐增大，最终导致岩体整体失稳破坏。逐渐增大，最终导致岩体整体失稳破坏。2 2、岩石的典型蠕变曲线及其特征、岩石的典型蠕变曲线及

4、其特征典型的蠕变曲线可分为典型的蠕变曲线可分为4（3）个阶段：个阶段： (1)(1)瞬时弹性变形阶段（瞬时弹性变形阶段（OA）：）： E00 (2)(2)一次蠕变阶段（一次蠕变阶段（AB）：）： （瞬态蠕变段（瞬态蠕变段/ /第第一蠕变阶段一蠕变阶段/ /初始蠕变段初始蠕变段/ /减速蠕变阶段减速蠕变阶段） 此阶段卸载此阶段卸载 一部分应变瞬时恢复一部分应变瞬时恢复(PQ段段) 一部分应变随时间逐渐恢复（一部分应变随时间逐渐恢复（OROR段）段）022 tdd 粘弹性粘弹性 (4)(4)三次蠕变阶段（三次蠕变阶段（CD）：应变速率迅）：应变速率迅速增加，直到破坏速增加，直到破坏（第三（第三蠕变

5、阶段蠕变阶段/ /加速蠕变段）加速蠕变段） 022 tdd (3)(3)二次蠕变阶段（二次蠕变阶段（BC）：应变速率不变）：应变速率不变 （第二蠕变阶段第二蠕变阶段/ /等速或稳定蠕变段）等速或稳定蠕变段） 022 tdd 此阶段卸载此阶段卸载 一部分应变瞬时恢复一部分应变瞬时恢复 一部分应变随时间逐渐恢复一部分应变随时间逐渐恢复 一部分应变不能恢复（一部分应变不能恢复（v v）粘弹塑性粘弹塑性 当应力水平当应力水平 较低时，可能无此阶段较低时，可能无此阶段（稳定蠕变）（稳定蠕变）蠕变变形总量：蠕变变形总量：=0 0+ +1 1(t)+(t)+2 2(t)+(t)+3 3(t)(t)式中式中：

6、0 0为瞬时弹性应变；为瞬时弹性应变；1 1(t)(t)，2 2(t)(t)，3 3(t)(t)为与时间有关的一次蠕为与时间有关的一次蠕变、二次蠕变、三次蠕变。变、二次蠕变、三次蠕变。v v 为粘塑性应变，为粘塑性应变， Q Q 为粘弹性应变。为粘弹性应变。3 3、岩石的蠕变曲线类型、岩石的蠕变曲线类型类型类型1 1：稳定蠕变：稳定蠕变 。曲线包含瞬时弹性变形、瞬态蠕变和稳定蠕。曲线包含瞬时弹性变形、瞬态蠕变和稳定蠕变变3 3个阶段（压应力个阶段（压应力10MPa10MPa，12.5MPa12.5MPa），无第三阶段蠕变），无第三阶段蠕变类型类型2 2：典型：典型蠕变蠕变 。曲线包含。曲线包含

7、4 4个阶段（压应力个阶段（压应力15MPa15MPa，18.1MPa18.1MPa）类型类型3 3：加速：加速蠕变蠕变 。曲线几乎无稳定蠕变阶段，应变率很高（压。曲线几乎无稳定蠕变阶段，应变率很高（压应力应力20.5MPa20.5MPa，25MPa25MPa）变形近似直线状急剧发展，迅速破坏）变形近似直线状急剧发展，迅速破坏3-3 岩石的流变模型岩石的流变模型 岩石的流变本构模型岩石的流变本构模型 ：用于描述岩石应力应变关系随：用于描述岩石应力应变关系随时间变化的规律。它是通过试验理论应用证实而得到的。时间变化的规律。它是通过试验理论应用证实而得到的。 本构模型分类：本构模型分类： 1 1、

8、经验公式模型、经验公式模型：根据不同试验条件及不同岩石种类求得：根据不同试验条件及不同岩石种类求得的数学表达式，这种表达式通常采用幂函数、指数函数、的数学表达式，这种表达式通常采用幂函数、指数函数、对数函数的形式表达。对数函数的形式表达。2 2、积分模型、积分模型：是在考虑施加的应力不是一个常数时的更一：是在考虑施加的应力不是一个常数时的更一般的情况下，采用积分的形式表示应力应变时间关系般的情况下，采用积分的形式表示应力应变时间关系的本构方程。的本构方程。3 3、组合模型、组合模型：将岩石抽象成一系列简单元件（弹簧、阻尼：将岩石抽象成一系列简单元件（弹簧、阻尼器、摩擦块），将其组合来模拟岩石的

9、流变特性而建立的器、摩擦块），将其组合来模拟岩石的流变特性而建立的本构方程。本构方程。 （一）经验公式模型（一）经验公式模型 45690. 0101610. 1t1 1、幂函数型、幂函数型 ： DtBte lg)( )10()( nAttn 式中：式中：A A和和n n是经验常数，其值取决于应力水平、材料物是经验常数，其值取决于应力水平、材料物理特性及温度条件。理特性及温度条件。 2 2、对数型、对数型 ： 式中：式中： e 为瞬时弹性应变；为瞬时弹性应变；B B，D D取决于应力性质及水平的待定常数。取决于应力性质及水平的待定常数。大理岩试验大理岩试验（ 轴向、侧向）轴向、侧向）第一、第二阶

10、段第一、第二阶段轴向蠕变方程轴向蠕变方程 45044. 0104205. 0t第一、第二阶段第一、第二阶段侧向蠕变方程侧向蠕变方程 通常形式通常形式通常形式通常形式如：如：式中：式中： e 为瞬时弹性应变；为瞬时弹性应变；A A为蠕变系数。为蠕变系数。罗伯逊（罗伯逊（RoberstsonRoberstson）根据开尔文（）根据开尔文（KelvinKelvin）粘弹性模型通过试）粘弹性模型通过试验曲线修正后得到半经验公式：验曲线修正后得到半经验公式：如：如：tAteln)(cnEA)(单轴单轴cnGA)2(31三轴三轴n nc c 为蠕变指数为蠕变指数 )(exp1)(tfAt 3 3、指指数型

11、数型 ： 式中：式中： A A为试验常数，为试验常数，f(t)f(t)是时间是时间t t的函数。的函数。通常形式通常形式如：如：伊文思（伊文思（EvansEvans）对花岗岩、砂岩研究得到：）对花岗岩、砂岩研究得到：)1exp(1 )(4 . 0ctAt式中：式中： A A、c c为试验常数为试验常数4 4、幂函数、幂函数、指指数函数、对数函数数函数、对数函数 混合型混合型： 如：如： 干燥的钙质石灰岩：干燥的钙质石灰岩：6651. 010)48lg52822()(ttt 干燥的白云质石灰岩：干燥的白云质石灰岩：6489. 049. 010)567 . 0648()(tett 干燥的砂岩：干燥

12、的砂岩：6687. 001. 010)410581858()(tett经验法简单实用，经验法简单实用，对特定的岩石试验而得，难以推广到对特定的岩石试验而得，难以推广到所有情况所有情况（二）积分模型（二）积分模型 （一维）流变方程：流变方程：0),(tf应力为常数，蠕变方程应力为常数，蠕变方程0),(tf应变为常数，松弛方程应变为常数，松弛方程蠕变方程蠕变方程应力随时间应力随时间的变化规律的变化规律每时刻在给定应力下的应变每时刻在给定应力下的应变应力不为常数时应力不为常数时)(0tJ蠕变方程蠕变方程恒定应力t的函数0 0时刻：作用应力：时刻：作用应力：0 0-t-t时刻：作用应力：时刻：作用应力

13、： 0 0+0 0 t t时刻：应变：时刻：应变：)()(00tJtJ设应力增量设应力增量作用在作用在0 0时刻：时刻：时刻的应变为：时刻的应变为：)0()(000tJtJ t ti i时刻应力相对前一时间时刻应力相对前一时间步步t ti-1i-1增加增加i i，相应，相应应变应变增量为：增量为：)(iiitJ 引起的引起的总应变：总应变：tiiitJtJ00)()( 积分形式积分形式：)()()()()()(0000dtJtJdtJtJtt积分形式的流变方程积分形式的流变方程)()()(00dtJtJt 令：令：vduuvudv)(),(ddvtJu 利用分部积分：利用分部积分：dtJtJJ

14、tJtdJtJtJtttt)()()0()()()0()()()()|()()()(00000dtJtJJtJtt)()()0()()()0()()(00EJ1)0(式中：式中： 为为0 0时刻的应力应变关系，有：时刻的应力应变关系，有： )0(J)(1)(tKEtJ 令：令：材料蠕变核材料蠕变核dtKtEtt)()()(1)(0积分形式的流变方程常用通式积分形式的流变方程常用通式dtRtEtt)()()()(0积分形式的松弛方程常用通式积分形式的松弛方程常用通式同理同理材料松弛核材料松弛核（三）组合模型（三）组合模型1 1、流变模型元件、流变模型元件（1 1）弹性介质及弹性元件（虎克体）弹性

15、介质及弹性元件（虎克体） ： 满足虎克定律满足虎克定律 E 弹性介质性质：弹性介质性质： （1）具有瞬时变形性质；）具有瞬时变形性质；（2）常数，则常数，则保持不变，故无保持不变，故无应力松弛性质；应力松弛性质；（3）常数，则常数，则也保持不变，故也保持不变，故无蠕变性质；无蠕变性质；（4）0（卸载），则（卸载），则0，无，无弹性后效。弹性后效。 可见，可见，、与时间与时间t无关无关(无蠕变无蠕变)（2）粘性介质及粘性元件（牛顿体）满足牛顿定律 dtdct t 0 加载瞬间，无变形加载瞬间，无变形即当即当t=0t=0时时, ,=0 0,=0,=0,则则 c=0粘性介质性质：粘性介质性质： （1

16、）当）当0时，时， 说明在受应力说明在受应力 0作用，要产生相应的变形作用，要产生相应的变形必须经过时间必须经过时间t，表明无瞬时变形，表明无瞬时变形，粘性元件具有蠕变性质；粘性元件具有蠕变性质；t 00 （2）0（卸载），则（卸载），则常数，故无弹性后效，有永久变形。常数，故无弹性后效，有永久变形。（3）常数，则常数，则=dd/dt/dt0，粘性元件不受力，故无应力松弛，粘性元件不受力，故无应力松弛性质。性质。牛牛顿顿粘粘性性系系数数（3）塑性介质及塑性元件（圣维南体） 当当: :s s ，=0=0 s s , , 可模拟刚塑性体的变形性质。 牛顿体具有粘性流动的特点。塑性元件具有刚塑性体变

17、形牛顿体具有粘性流动的特点。塑性元件具有刚塑性体变形（塑性变形也称塑性流动）的特点。（塑性变形也称塑性流动）的特点。粘性流动粘性流动：只要有微小的力就会发生流动。：只要有微小的力就会发生流动。塑性流动塑性流动：只有当应力：只有当应力达到或超过屈服极限达到或超过屈服极限s才会产生才会产生流动。流动。粘弹性体粘弹性体：研究应力小于屈服极限时的应力、应变与时间：研究应力小于屈服极限时的应力、应变与时间的关系；的关系；粘弹塑性体粘弹塑性体：研究应力大于屈服极限时的应力、应变与时：研究应力大于屈服极限时的应力、应变与时间的关系；间的关系；2、岩石的组合流变模型 基本元件基本元件岩石一种性质岩石一种性质岩

18、石性质岩石性质弹性弹性塑性塑性粘性粘性组合组合组合方式组合方式串联串联并联并联应力：应力：n 21应变：应变：n 21应力：应力：n 21应变：应变：n 21（1 1）理想弹塑性介质模型）理想弹塑性介质模型 os 当当: :s s ，=s s , , 保持不变，保持不变， 持续增大，持续增大，。E 无蠕变、无松弛、无弹性后效无蠕变、无松弛、无弹性后效（2 2）理想粘塑性介质模型）理想粘塑性介质模型 1 1+2 2， 1 1 =2 2:0:11ss2阻尼器：阻尼器：摩擦片：摩擦片：ssss:0:本构方程本构方程sss:0:A A、蠕变曲线：、蠕变曲线：当当保持不变，保持不变， 即即 0 0常数常

19、数通解为：通解为：cts0初始条件：初始条件：加载瞬间加载瞬间00时，tts0ss:考虑：蠕变方程：蠕变方程：蠕变曲线蠕变曲线B B、卸载曲线：、卸载曲线：当当t=tt=t1 1时卸载，时卸载，10ts 由于：由于：1 1 =2 2 1 1 为塑性变形，故为塑性变形，故为为永久变形，永久变形， 无弹性后效无弹性后效C C、松弛曲线、松弛曲线：当当保持不变，保持不变，d/dtd/dt=0,=0, 代入代入ss 无松弛无松弛（3）马克斯威尔模型（Maxwell）（粘弹性体） 该模型由弹性元件和粘性元件串联而成，可模拟变形随该模型由弹性元件和粘性元件串联而成，可模拟变形随时间增长而无限增大的力学介质

20、。时间增长而无限增大的力学介质。1 1 2 2 设弹簧和粘性元件的应力、设弹簧和粘性元件的应力、应变分别为应变分别为1 1,1 1和和 2 2 ,2 2，组合模型的总应力为组合模型的总应力为和和。弹簧弹簧: :tttdddddd21 E11 tEtEtdd1dd1dd11 22ddt由由( (b):b):粘性粘性元件元件: :则则 1 12 2， (a) (a) 1 1 2 2 (b)(b) tE dd1马克斯威尔模马克斯威尔模型本构方程型本构方程马克斯威尔模型本构方程：A A、蠕变曲线：、蠕变曲线：当当保持不变，保持不变， 即即 0 0常数，常数，d/dtd/dt=0,=0, 代入上式得：代

21、入上式得： tEtdd1dd 0dd t通解为：通解为：ct 0初始条件：初始条件：加载瞬间加载瞬间Et000 时，时， 得：得： c = c = 0 0tEt 0000 蠕变方程：蠕变方程：0 1 1 2 2 B B、卸载曲线：、卸载曲线：当当t=tt=t1 1时卸载，弹时卸载，弹性变形性变形0 0立即恢复，则卸载曲线立即恢复，则卸载曲线为：为：10t 这是不可恢复的塑性变形。这是不可恢复的塑性变形。tEt 0000 1 1 2 2 0 C C、松弛曲线、松弛曲线：当：当保持不变，保持不变， 即即0 0常数，常数，d/dtd/dt=0,=0, 代入上式得：代入上式得： tEtdd1dd01

22、dtdE通解为：通解为：ctE ln初始条件：初始条件：00 时，时， t得：得：c = lnc = ln0 0tEe 0松弛方程：松弛方程：1 1 2 2 可见：马克斯威尔模型可见：马克斯威尔模型具有瞬时变形、蠕变和松弛的具有瞬时变形、蠕变和松弛的性质，可模拟变形随时间增长而无限增大的力学介质。性质，可模拟变形随时间增长而无限增大的力学介质。1 2 2 1 （4）开尔文(Kelvin)模型 设弹簧和阻尼元件的应力、应设弹簧和阻尼元件的应力、应变分别为变分别为1 1 、1 1和和2 2 、2 2，组组合模型的总应力为合模型的总应力为和和 。弹簧弹簧: : EE 11ttdddd22 由由( (

23、a):a):阻尼元件阻尼元件: :则则 1 1+2 2， (a) (a) 1 1 =2 2 (b)(b)开尔文开尔文模型本构方程模型本构方程tEdd (c)(c)(d)(d)开尔文模型本构方程：A A、蠕变曲线、蠕变曲线：当：当保持不变，保持不变， 即即 0 0常数，代入上式得：常数，代入上式得：通解为：通解为：tEAeE 0初始条件：初始条件：加载瞬间，粘性元件不加载瞬间，粘性元件不变形，即变形，即00时，时， t得：得：)1(0tEeE 蠕变方程蠕变方程：tEdd tEdd0 EA0 ( (c)c)1 2 2 1 ( (e)e)1(0tEeE 可见：可见：当当t=0t=0时时, , =0=

24、0，当当t t 时，时， 0 00 0/E ,/E ,即弹性变形即弹性变形趋于常数趋于常数 ( (e)e)蠕变方程1 2 2 1 凯尔文模型能模拟稳定蠕凯尔文模型能模拟稳定蠕变，不能模拟瞬时弹性变形。变，不能模拟瞬时弹性变形。若在若在t tt t1 1 时卸载，时卸载，0 0，由本构方程：由本构方程：0dd tE B、卸载曲线方程tEdd 得：得：通解为：通解为：tEAe 1 2 2 1 初始条件：初始条件：11ttt ，得卸载曲线：得卸载曲线：11tEteA 当当t t时，时，0,0,即卸载后，变形慢慢恢复到即卸载后，变形慢慢恢复到0 0（弹性（弹性后效）。后效）。得得: :)(11ttEt

25、e 1 2 2 1 弹性后效弹性后效C C、松弛曲线：、松弛曲线：当当保持不变，保持不变， 即即0 0常数，常数，d/dtd/dt=0,=0, 代入开尔文代入开尔文模型本构方程模型本构方程0 E 可见，应变恒定，应力恒定。即无应力松弛现象可见，应变恒定，应力恒定。即无应力松弛现象该模型反映了有弹性后效现象和稳定蠕变，无应力松弛的该模型反映了有弹性后效现象和稳定蠕变，无应力松弛的性质。性质。开尔文开尔文模型是一种粘弹性模型。模型是一种粘弹性模型。tEdd 1 2 2 1 （5）广义开尔文模型 一个弹簧一个弹簧 一个开尔文元件一个开尔文元件弹簧弹簧: :111E 1 12 2= = (a) (a)

26、 1 1+2 2 = = (b)(b)E E1 1, ,1 1E E2 2, ,2 2, ,2 2串联串联111E广义开尔文广义开尔文模型本构方程模型本构方程E E1 1, ,1 1E E2 2, ,2 2, ,2 2)()()()(112112EEEE开尔文开尔文体体: :222 E即即: :2121)1 (EEEE整理整理: :当当保持不变，即保持不变，即 0 0常数，常数，代入本构方程代入本构方程)1 (22010tEeEE蠕变方程蠕变方程：100EE E1 1, ,1 1E E2 2, ,2 2, ,2 2A A、蠕变曲线：、蠕变曲线：初始条件：初始条件：解微分方程解微分方程100Et

27、=0t=02010EEt t蠕变范围蠕变范围：10E2010EE广义开尔文蠕变曲线广义开尔文蠕变曲线10E10E20E若在若在t tt t1 1 时卸载，时卸载，0 0，由本构方程：由本构方程：B、卸载曲线方程10E10E20E通解通解初始条件：初始条件：特解特解C C、松弛曲线：、松弛曲线：当当保持不变，保持不变， 即即0 0常数，常数，d/dtd/dt=0,=0, 代入代入本构方程本构方程02122121)(211tEEeEEEEEEE（6）宾汉姆(Bingham)模型 一个弹簧一个弹簧 一个理想粘塑性元件一个理想粘塑性元件串联串联EE 11串联串联21,s串联串联22,1,E弹簧：弹簧：

28、理想粘塑性元件：理想粘塑性元件：sss:0:121本构方程本构方程21sssEEE:,:1当当保持不变，即保持不变，即 0 0常数常数A A、蠕变曲线：、蠕变曲线：sss:11变只有弹簧有变形，无蠕Ets00蠕变方程蠕变方程：解微分方程解微分方程E00初始条件：初始条件：微分方程微分方程B B、松弛方程：、松弛方程：当当保持不变保持不变即即0 0常数常数d/dtd/dt=0=00:101sssEE无松弛常数，本构方程本构方程sssEEE:,:10sEtEsse)(0松弛方程松弛方程：解微分方程解微分方程0初始条件：初始条件：微分方程微分方程松弛范围：松弛范围：s0（7）伯格斯（Burgers）模型串联串联马克斯威尔模型马克斯威尔模型开尔文模型开尔文模型串联串联E E1 1E E2 2开尔文开尔文体体: :1111 E马克斯威尔模型本构方程马克斯威尔模型本构方程：222E2121111121E)()(212211EE 12122121112212)(EEEEEEEE数学变换数学变换本构方程本构方程 12122121112212)(EEEEEEEE本构方程本构方程当当保持不变，即保持不变，即 0 0常数常数A A、蠕变方程：、蠕变方程：200E初始条件：初始条件：0210)11()1 (11102020t