第7章 本构与强度理论
2、本构关系的复杂性
塑性阶段本构关系包括三组方程：
1）屈服条件：塑性状态的应力条件。
2）加-卸载准则：材料进入塑性状态后继续塑性变形或 回到弹性状态的准则，通式写成：
( ij , Ha ) 0
(i x, y, z， ; j x, y, z )
a。 1, 2
（7-7）
式中： 垂直于 轴的平面上平行于 轴的应力 j i ij 为某一函数关系， H 为与加载历史有关的参数， a 3）本构方程：
表7-1 几种岩石长期强度与瞬时强度比值
岩石名称
s / s0
粘土 0.74
石灰石 0.73
盐岩 0.70
砂岩 0.65
白垩 0.62
粘质页岩 0.50
§7.3 岩石强度理论与破坏判据
岩石强度理论：研究岩石在一定的假说条件下在各种
应力状态下的强度准则的理论。 强度准则：又称破坏判据，是表征岩石破坏条件的应 力状态与岩石强度参数间的函数关系，可用如下的方程 表示:
据广义虎克定理有：
x
1 ( ) x y z E 1 y ( ) y z x E 1 z ( ) z x y E 1 1 1 yx yz , zx zz , xy xy G G G
①弹性本构关系：线性弹性、非线性弹性本构关系。 ②弹塑性本构关系：各向同性、各向异性本构关系。 ③流变本构关系：岩石产生流变时的本构关系。流变
性是指如果外界条件不变，应变或应力随时间而变化的 性质。 岩石强度理论：研究岩石在一定的假说条件下在各种 应力状态下的强度准则的理论。
岩石的强度是指岩石抵抗破坏的能力。 岩石材料破坏的形式:断裂破坏、流动破坏（出现显 著的塑性变形或流动现象）。断裂破坏发生于应力达 到强度极限，流动破坏发生于应力达到屈服极限。 岩体的力学性质可分为变形性质和强度性质， 变形性质主要通过本构关系来反映，强度性质主要 通过强度准则来反映。 本章分别研究岩石、岩体的本构关系与强度理论。
中的E换成
E ， 2 v 1
换成
。
1
2. 空间问题弹性本构方程
1 x x ( y z ) E 1 y ( ) y z x E 1 z z ( x y ) E 2(1 ) 2(1 ) yz yz , zx zx E E 2(1 ) xy xy E
或均匀收缩（软化）。如果 是初始屈服面，后继屈服面为： f* 0
f f * ( ij ) H ( ) 0
②随动硬化模型：塑性变形发展时，屈服面的大小和形状保持不变，
仅是整体地在应力空间中做平动，后继屈服面为：
③混合硬化模型：介于等向硬化-软化和随动硬化之间的模型，后继
屈服面为：
f f * ( ij ijp ) H ( ) 0
在一系列的岩石流变试验基础上建立反映岩石流变性质 的流变方程，通常有二种方法：
1．经验法
岩石蠕变经验方程：
ε
D ε 3( t)
ε 2( t) ε 1( t)
A ε0
O
C
(t ) 0 1 (t ) 2 (t ) 3 (t )
B t
(t )为时间 t 的应变； 式中： 0 瞬时应变； (t ) 等速段应变； 1 (t )初始段应变； 3 (t )加速段应变。
p ij
F 其总应变增量表示为： d ij D d ij ij
1
由一致性条件可推出待定有限量 为：
1 F d ij A ij
对于系数A： 理想塑性材料： 硬化材料：
A0
F Q F Q A D ij p ij ij kl u
p 式中： ij 为总应力， ij 为塑性应力， 为标量的内变量，它可以代表塑性功，塑性体
积应变，或等效塑性应变。
屈服面：屈服条件在几何上可以看成是应力空间中的超曲面。
初始屈服面和后继屈服面。
分类：按塑性材料屈服面的大小和形状是否发生变化。 理想塑性材料（不变化）和硬化材料（变化）。
（ 1）硬化材料的屈服面模型 ①等向硬化-软化模型：塑性变形发展时，屈服面作均匀扩大（硬化）
2
图7-5 岩石的典型蠕变曲线
典型岩石蠕变方程：幂函数方程、指数方程、幂指数对数混合方程
2. 理论模型模拟法
将介质理想化，归纳成各种模型，模型可用理想化的具有基本性能（包括 弹性、塑性和粘性）的元件组合而成。
岩石的长期强度:由于流变作用，岩石强度随外载 作用时间的延长而降低，通常把作用时间 t 的强 度（最低值）称为长期强度。 对于大多数岩石，长期强度/瞬时强度（ ） 一般为 s /0.4 s0 ～0.8，软的和中等坚固岩石为0.4～0.6， 坚固岩石为0.7～0.8。表7-1中列出某些岩石瞬时 强度与长期强度的比值。
ij R( ij ) 或 d ij R(d ij )
式中：R为某一函数关系
1）岩石屈服条件和屈服面
初始屈服条件：从弹性状态开始第一次屈服的条件。
f ( ij ) 0
后继屈服条件：当产生了塑性变形，屈服条件的形式发
生了变化的屈服条件。
f ij , ijp , 0
2）增量理论：描述应力和应变增量间关系的理论，又称 流动理论。 当应力产生一无限小增量时，假设应变的变化可分成弹性 的及塑性的两部分：
d ij d d
e ij
p ij
弹性应力增量与弹性应变增量之间仍由常弹性矩 阵D 联系，塑性应变增量由塑性势理论给出，对弹 塑性介质存在塑性势函数Q，它是应力状态和塑性应 变的函数，使得： Q p
式中：E为物体的弹性模量； 为泊松比；G为剪切弹性模量，
E G 2(1 )
对于平面应变问题：因εz ＝γzx ＝γyz ＝ 0 ，故τyz
＝τzx＝0，可知：
z ( x y )
y x 1 1 2 y x y E 1 2(1 ) xy xy E 1 2 x E
对于平面应力问题：σz＝τ zx＝τ zy＝0
1 x x y E 1 y y x E 2(1 ) xy xy E
对比平面应力问题与平面应变的本构方程，
可以看出，只要将平面应力问题的本构关系式
力的平衡关系（平衡方程） 位移和应变的关系（几何方程）+ 边界条件 = 应力和应变的关系（物理方程 或本构方程）
应力场
位移场
依据适合的强度理论，判断岩体的破坏及其破坏形式。 岩体本构关系:指岩体在外力作用下应力或应力速率与其应变 或应变速率的关系。
岩石或岩体的变形性质：弹塑性或粘弹塑性。 本构关系：弹塑性或粘弹塑性本构关系。 本构关系分类：
f f * ( ij ijp ) 0
式中：α 是材料参数。
（2）塑性岩石力学最常用的屈服条件
库仑（Coulomb）、德鲁克-普拉格（Drucker-Prager）屈服条件。Leabharlann 2. 塑性状态的加-卸载准则
塑性加载：对材料施加应力增量后，材料从 一种塑性状态变化到另一种塑性状态，且有新的 塑性变形出现； 中性变载：对材料施加应力增量后，材料从一 种塑性状态变化到另一种塑性状态，但没有新的 塑性变形出现； 塑性卸载：对材料施加应力增量后，材料从塑 性状态退回到弹性状态。
3.本构方程
塑性状态时应力-应变关系是多值的，取决材料性质和加卸载历史。 1）全量理论：描述塑性变形中全量关系的理论，称形变理 论或小变形理论。 汉基（Hencky）、依留申等依据类似弹性理论的广义胡克定 律，提出如下公式 ：
xx m 2G xx m , xy G xy yy m 2G yy m , yz G yz zz m 2G zz m , zx G zx
三、 岩石塑性本构关系
塑性本构关系特点： 1、应力-应变关系的多值性 同一应力有多个应变值与σ 它相对应。 本构关系采用应力和应变 增量的关系表达。 塑性状态描述：除用应力、 应变，还需用塑性应变，塑 0 性功等内状态变量来刻画塑 性变形历史。
B A 加载 卸载 C
ε
图7-1 加－卸载应力－应变曲线
式中： u 为塑性功，这样加载时的本构方程为：
可以唯一地确定应变增量 d ij 。 应用增量理论求解塑性问题，能够反映应变历史对塑性 变形的影响，因而比较准确地描述了材料的塑性变形规律。
1 1 Q F d ij D d kl A ij kl p 对任何一个状态 kl , kl , u ，只要给出了应力增量，就
σ1= f (σ2 ,σ3 ,σC ,σt ,C ,Ф )
或处于极限平衡状态截面上的剪应力 和正应力 间 的关系方程：
f
§7.3 岩石强度理论与破坏判据
一、 库仑强度准则
二、 莫尔强度理论 三、 格里菲斯强度理论 四、 Griffith强度准则的三维推广 （Murrell强度准则） 五、 德鲁克一普拉格准则
四、岩石流变理论
流变：指材料的应力-应变关系与时间因素有关的性 质，材料变形过程中具有时间效应的现象称为流变现 象。 蠕变：当应力不变时，变形随时间增加而增长的现 象。 松弛：当应变不变时，应力随时间增加而减小的现 象。 弹性后效：加载或卸载时，弹性应变滞后于应力的 现象。
蠕变试验表明：
1)当岩石在某一较小的恒定荷载 持续作用下，其变形量虽然随时 间增长有所增加，但蠕变变形的 速率则随时间增长而减少，最后 变形趋于一个稳定的极限值，这 种蠕变称为稳定蠕变。
岩体本构关系与强度理论
Rockmass Mechanics
第七章 岩体本构关系与强度理论
概述 岩石的本构关系 岩石强度理论 岩体变形及本构关系 岩体破坏机制及破坏判据
概述
－－岩体的力学性质表现为弹性、塑性和粘性或三者之间的组合，如
粘弹性、弹粘性、弹塑性、弹塑粘性等。
如何求解岩体的基本力学问题呢？
3
(12c)1
1 3c
3(12)t
那么最大正应变强度判据可写成：
f C tg
1 1 3 c
当岩石应力条件满足以上判据时，岩石发生张破裂。取其中最小的
比较。 该强度理论适用于无围压和低围压及脆性岩石。
0
0
2. （二次）抛物线型 岩性较坚硬至较弱的岩石，
如泥灰岩、泥岩、砂岩、泥页 岩、页岩等岩石的强度包络线 近似于二次抛物线。
3．双曲线型 岩性坚硬、较坚硬的岩石，
如砂岩、灰岩、花岗岩等，其强 度包络线近似于双曲线。
二次抛物线型莫尔强度包络线 双曲线型莫尔强度包络线
0
说明：莫尔－库仑判据（强度理论）实质上是一种剪应力强度判据（理 论），既适用于塑性岩石，又适用于脆性岩石的剪切破坏，应用很广。
那么，岩石承受低于瞬时强度的荷载作用，是否会破坏？
只要收到长期荷载作用下，由于流变作用，岩石完全可能发生破坏。 而且岩石强度随外荷载作用时间的延长而降低，通常将时间t→∞的强度 （最低值）称为岩石的长期强度。
对大多数岩石，长期强度/瞬时强度（S∞/S0）一般为0.4～0.8，软的 和中等坚固岩石为0.4～0.6，坚固岩石为0.7～0.8。
2.莫尔－库仑强度理论
（Coulomb-Mohr Strength theory）
岩石的强度理论与本构关系
岩石的强度理论与本构关系朱浮声(东北大学土木系,沈阳110006)朱浮声,1948年6月生于黑龙江齐齐哈尔11976年毕业于东北大学,1983年获中国矿业大学工学硕士学位,1991年获东北大学博士学位11988年曾在美国南伊利诺大学作访问学者,1993年在瑞典皇家工学院任客座教授1现任东北大学土木工程系教授,辽宁省力学学会理事1主要研究方向为计算岩土力学和岩土加固技术1在国内外学术刊物上发表论文50余篇,出版5锚喷加固设计方法6等学术专著2部,译著1部1摘要本文简要介绍了岩石强度理论和本构关系的发展和现状,讨论了它们不同的特点与适用条件1关键词岩石,岩体,强度理论,本构关系1前言随着电子计算机的飞速发展和计算技术的逐步完善,对岩石强度理论和本构关系提出了更高要求,以便更真实描述岩石和岩体力学特征,求解复杂的工程岩石力学问题1由于岩石材料力学性质的某些相似性和其它历史原因,岩石强度理论和本构关系的早期研究曾大量引用了土力学成果,并提出了一些适用于岩土介质的强度理论和本构关系1随着岩石力学的发展,人们认识到,岩石和岩体的物理力学性质不仅有别于其它非摩擦工程材料,而且,与土或混凝土等摩擦材料也存在较明显差异1例如,岩石破坏包括脆性、延性及由脆性向延性转化等复杂类型;岩体的力学特性受控于岩块和不连续面的力学特性;岩石工程的稳定性通常受主要不连续面控制等1因此,近年来又提出了适用于岩石、不连续面和岩体的强度理论或本构方程式1本文旨在介绍这些理论研究的最新进展,并对已有岩土强度理论和本构关系的适用条件和局限性加以简要评价1限于篇幅,本文仅涉及与时间无关的各向同性和等向强化模型12岩土共用的强度理论和本构关系211弹性均质、各向同性或横观各向同性模型曾被广泛用于描述岩土力学特征,特别是峰值强度前的应力-应变关系,并得到了大量解析解和实用近似解1考虑到应力-应变曲线的明显非线性特性,曾将非线性弹性理论与计算机技术相结合,提出了一批数值算法,并在60~70年代的岩土力学分析中不断被引用1例如,以曲线各点的割线模量取代弹性常数,构成了各种超弹性模型[1],或以增量形式描述非线性弹性应力-应变关系,形成了亚弹性模型[2]等1但是,由于这些模型只考虑到岩土材料的弹性特征,并且,随着模型阶次增高,待定常数的数目往往过多,因而,限制了它们的广泛应用1212 理想塑性强度理论在计算岩土力学中,广泛采用了莫尔-库伦强度准则(Mohr -Coulomb)和德鲁克-普拉格准则(Drucker -Prager)1莫尔-库伦准则可以表述为R 1-B R 3=C (1)式中,B 和C 一般是常数1在主应力空间,式(1)表示一个以静水应力轴为中心轴,具有不规则六角形截面的角锥体表面(图1)1这个准则由于较好地表征了岩土介质在压缩条件下的某些弹塑性力学特征,因而得到了较广泛的应用1但是,由于忽略了中间主应力对破坏的作用,存在明显的缺陷1另外,由于屈服面在三维应力空间中存在/角隅0,给数值计算带来了诸多困难1为了解决上述问题,曾对莫尔-库伦准则进行修正,将米赛斯准则(M ises)加上一个静水应力因子,形成了著名的德鲁克-普拉格准则,不仅考虑到三个主应力对破坏的影响,并且消除了屈服面存在的角隅1这个准则可表述为A J 1+J 2=C (2)式中,J 1和J 2是主应力不变量,A 和C 是正常数1图1 莫尔-库伦准则和德鲁克-普拉格准则屈服面在主应力空间,式(2)是一个以静水应力轴为中心轴的圆锥体(图1)1它虽然克服了莫尔-库伦准则的上述缺点,但在破坏状态下,该准则给出了较大的材料体积膨胀,这与岩土介质的试验结果明显不符1我国学者俞茂由正交八面单元体的三个主应力出发,提出了双剪强度理论和适用于岩土体的广义双剪强度理论(包括屈服准则)[3],并得到了双剪统一强度理论[4]R 1-A 1+b (b R 2+R 3)=R t ,R 2F R 1+A R 31+A11+b (R 1+b R 2)-A R 3=R t ,R 2E R 1+A R 31+A (3)式中,A 和b 是常数,R t 是材料单轴抗拉强度1在主应力空间,式(3)表示一个以静水应力轴为中心轴,具有不等边十二边形截面的锥体表面1可以证明,广义双剪理论和莫尔-库伦准则在P 平面上的屈服曲线分别是各种岩土屈服准则的上限和下限1213应变硬化(软化)一般地,岩土体应力状态满足屈服准则时,将出现屈服应力随变形增大而不断增高(硬化)或降低(软化)现象1对于前者,屈服面在主应力空间是连续扩大的;对于后者,则表现为屈服面的不断收缩1当满足破坏条件时,将形成屈服面与破坏面(残余破坏面)相互重合,而屈服面与破坏面始终相一致的情况仅发生在完全塑性材料中1因此,为了建立岩土介质完整的本构关系,必须同时考虑屈服准则、流动法则和软(硬)化定律等三方面1其中,对材料硬(软)化特图2帽盖模型屈服面性的研究多借助于控制该材料硬化特性的屈服面,称之为硬化帽盖(图2)1根据不同的岩土介质和试验,提出了不同形状的帽盖[5],其一般表达式为f1(J1,J c2,k1)=0(4)式中,J c2为应力偏量第二不变量,k1为硬化参数1除了屈服帽盖,岩土帽盖模型还包括一个固定屈服面,例如,通常以初始德鲁克-普拉格破坏面与米赛斯屈服面光滑相接表示1一般地,固定屈服面取为强度理论限定的破坏面f2(J1,J c2)=0(5) 3岩石和岩体强度理论与本构关系如前所述,在一定条件下,可以使用相同强度理论分析岩土力学问题1但在一般情况下,岩石的破坏面具有如下特征:(1)在主应力空间,破坏曲面在原点附近的顶角是张开的;(2)岩石破坏包络线,即破坏面在伦杜列克面(Rendulic)上的子午线不是直线,而是曲线;(3)岩石有抗拉强度1通常,前述岩土体的屈服和强度准则都可以满足条件(1)和(3),为了满足条件(2),需要进行必要修正1311岩土强度理论的修正为了使强度理论满足上述条件(2),从而应用于岩石力学问题分析中,早期的工作多采用对破坏曲线近似逼近方法,例如,以双曲线或抛物线取代莫尔-库伦准则中的直线等1更一般的方法是直接采用莫尔强度理论,并通过对P平面上多边形屈服曲线角点的光滑化得到各种角隅模型[6]1典型范例是关于岩石的吉姆-拉德破坏准则(Kim-Lade)1拉德曾提出如下土体两参数破坏准则[7](J21/J3-27)(J1/P a)m=G1(6)式中,P a是该应力状态下大气压力,m和G1是常数,其中,破坏面在原点附近的张角随G1变化,而子午线曲率随m值变化(图3)1在主应力空间,拉德准则是一个以静水应力轴为中心轴,具有带圆角三角形截面的子弹头形曲面,该曲面顶点位于原点1图3拉德破坏准则破坏面[7]为了得到适用于岩石的强度准则,吉姆和拉德对式(6)进行了修正,即考虑岩石凝聚力和抗拉强度的作用,在式(6)的法向主应力分量叠加一个常应力项aP aR x=R x+aP aR y=R y+aP a(7)R z=R z+aP a研究式中,a是一个无因次常数,aP a的值反映了岩石的抗拉强度1不难看到,这个三参数强度准则较好地反映了岩石破坏面的上述3个特征,同时,原作者通过87组不同岩石的试验数据对模型进行了多次验证1312岩石的脆性破坏准则岩石三轴试验结果表明,在侧限压力较低时,岩石试件的破坏应力随变形增大而不断降低,在很小或完全未出现永久变形的情况下发生突然的脆性破坏1随着侧限压力增大,通常出现由脆性向延性破坏的转变,这种现象可以由塑性变形机制来解释(岩石破裂流动与颗粒滑移等)1岩石脆性破坏准则研究仍处于发展阶段1其中,格里菲斯理论(Griffith)是一个基于理想脆性假定的二维准则,由此理论预测的脆性材料单轴抗压与抗拉强度R c和R t的关系式R c=-8R t1由于岩石裂纹随围压增高将出现闭合,此时应考虑闭合裂纹表面间摩擦作用,因此,提出了关于岩石修正的格里菲斯理论[8]S=2R t+R n tg<(8)式中,S和R n为裂纹面上切应力与法应力,<为内摩擦角1显然,这种修正是将低应力条件下的格里菲斯理论转化为高应力条件下的莫尔理论1M urrell考虑到中间主应力的作用,提出了一个三维脆性破坏准则,这个准则预测R c/P a=12 |R t/P a|1在主应力空间,这个准则表示为一个以静水压力轴为中心轴的椭球面与一个处于拉应力区的四棱锥面相切得到的曲面(图4)1需要指出,这个准则虽计及中间主应力影响,并具有弯曲的子午线,但它的基本出发点却是基于单轴抗拉强度判据1图4M urrell三维脆性破坏准则破坏面[9]313岩石破坏的经验准则由于岩石和岩体力学特征的复杂性,针对不同岩石和荷载条件提出了大量实用经验准则1其中,霍克-布朗准则是应用最广的1这个准则依据格里菲斯和修正格里菲斯理论的基本概念,采用/试凑法0得到了分别适用于岩石和岩体的经验准则[10]R1=R3+m R c R3+S R c2(9)式中,m和S是表征岩石或岩体性质及其破坏范围的常数1在主应力空间,这是一个以静水应力轴为中心轴、具有6条抛物线围成的6边形截面的锥体表面(图5)1这个准则给出R c=-(7~25)R t,这与大量试验结果接近,因而,得到较广泛应用1图5霍克-布朗经验破坏准则破坏面4节理和节理岩体天然岩体由节理和岩块组合而成1对于起控制作用的节理,通常采用/节理单元0来模拟1早期的节理单元是一个非线性弹性模型,给出了节理面两侧力-位移的增量表达式1为了考虑节理延性和切向-法向作用的相互影响,普遍采用了遵循莫尔-库伦准则的弹塑性节理模型1但是,如果采用关联流动法则,这个模型将产生一个无法消除的剪胀率1因此,罗伯茨等(Roberts)建议用非关联流动法则,相应塑性势函数Q为Q=|R s|-R n tg W(10)式中,R s和R n是节理切向和法向应力,W是节理剪胀角,可由试验确定1由于试验水平和理论的限制,节理面的理论模型尚不成熟,在应用中最可靠、最广泛的是巴登(Barton)提出的经验准则[11]F=|R s|-R n tg<(JRC lg(JCS/R n)+<r)=0(11)式中,JR C是节理面粗糙度系数,JCS是节理面抗压强度,<和<r分别是节理面摩擦角和残余摩擦角1这些参数都可以由简单原位试验得到1国际岩石力学学会制订了规范用于获取和解释这些参数,从而,推动了这一准则的普遍应用1对于受多组节理切割的岩体,由于很难同时模拟这几组节理,通常需找到节理岩体的本构关系1目前,此项研究仍处于开始阶段,应予充分重视1对于等距排列的平行节理(未充填),若节理连续分布且尺寸远小于岩体或结构物尺寸,提出了/节理岩体层状模型0(Multilaminate model)[12]1这实质上是一种等效材料模型,在最终形成的弹塑性或粘弹塑性本构关系中,以不同力学模型分别描述各层岩石和节理面的力学特征,并同时考虑它们对岩体力学的影响,得到节理岩体总的粘塑性应变速率ÛE VP=C i3F i45Q i5R+E nJ=1C J3F J45Q J5R J5R J5R(12)式中,F i和F J分别是第i层岩石和第J层节理的屈服(破坏)函数,Q i和Q J为相应塑性势,R J 表示J组节理面上法向和切向应力,C i和C k是相应粘性参数1式(12)中等号右端第一项与岩石特性有关,第二项则涉及n个节理面的力学特征1如果不考虑岩体的流变特性,采用关联流动法则,可以给出弹塑性节理岩体的类似表达式15结束语(1)非线性弹性模型曾在岩土力学中应用1由于高次模型待定常数过多,且为区分加、卸载情况需给出复杂应力状态下加载条件,限制了它们的使用范围,在岩石力学中应用较少,并主要用于比例加载条件下1(2)莫尔-库伦准则和德鲁克-普拉格准则在岩土力学分析中得到广泛应用1前者的缺点是忽略了中间主应力的作用,并且,在三维主应力空间存在屈服面角隅,给计算带来了困难1后者虽然克服了上述问题,但在破坏状态下给出较大体积膨胀,这与岩土试验结果严重不符1广义双剪强度理论及其角隅模型展示了广阔应用前景,但需大量试验与工程验证1各种帽盖模型考虑到岩土介质的应变硬(软)化,计及剪胀或剪缩,但公式推导中加入种种补充假定,模型的可靠性需进一步验证1(3)岩石、岩体和土体的强度理论和本构关系相似又相区别1吉姆-拉德的三参数岩石破坏准则经过87组岩样检验1有较高可信度1各种岩石脆性破坏准则都源于单向抗拉强度判据,它们的可靠性有待检验1岩石和岩体经验破坏准则在应用中占重要地位,其中,霍克-布朗准则适用于R1>314R3条件下延性岩石(体),在无控制作用节理存在的岩体工程分析中得到普遍应用1(4)在主要节理面的模拟中采用了莫尔-库伦准则和关联/非关联流动法则1已提出的8参数节理模型可用于研究节理加载-卸载-再加载过程[12],但最可靠、应用最广的仍是巴登经验准则1对无控制性不连续面的节理岩体提出了等效模型和相应本构关系1节理岩体的强度理论和本构关系研究仍处于初始阶段,是目前主要研究方向之一1参考文献1Fung Y C1Foundations of Solid M echanics1Prentice-Hall,19652T ruesdell C1Hypoelas ticity1J Ration M ech A nal,1955,4:83~1333俞茂等1双剪应力强度理论及其推广1中国科学A辑,1985,28(11)4俞茂1统一强度理论及其应用1强度理论研究新进展1西安:西安交通大学出版社,1993133~445S chofield A N,Worth C P1Critical S tate Soil M echanics1M cGraw-Hill Book Company,19686Zienkiew icz O C,Pande G N1Som e useful forms of isotropic yield surfaces for soil and rock mechanics1Finite Elements in Geome-chanics1Gudehus G(eds)1John Wiley&Sons,19777Kim M K,Lade P V1M odelling rock strength in three-dimensions1Int J Roc k M ec h M in S ci&Geomech,1984,21:21~338M cCli ntok F A,Walsh J B1Fri ction on Griffith cracks under pressure1Proc4th US Nat Congr Appl M ech,196211015~10219M urrel l S A F1A cri terion for bri ttle fracture of rocks and concrete under triaxial stress and the effect of pore pres sure on the creter-i on1Proc5th US S ym p Rock M ech,Pergamon Press,19631563~57710Hoek E,Brown E T1岩石地下工程1连志升等译,北京:冶金工业出版社,198611Barton N R,Choubey V1The shear strength of rock joints in theory and practice1Rock M ech,1977,10:1~5412Pande G N,Beer G,Williams J R1Numerical M ethods in Rock M echanics1W i ley,199013Pande G N1A constitutive model of rock joints1Proc Int Symp Fund Rock Joints1Center Pub11985(1996年5月31日收到第1稿,1996年8月12日收到修改稿)(上接第7页)29Levy D,Powell K,van Leer B1An i m plementati on of a grid-i ndependent upwind scheme for the Euler equati ons1AIAA89-1931-CP 30Rumsey C L,van Leer B,Roe P L1A grid-independent approximate Riemann solver w ith applications to the Euler and Navier-Stokes equations1J Comput Physics,1993,105(2):306~32331Roe P L1Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics1J Comp ut Phys,1986,63: 458~47632Roe P L1Discontinuous sol utions to hyperbolic system s under operator splitting1Nu merical M ethods f or Par tial Diff erential Equa-tions,1991,7:277~29733Lacor C,Hirsch Ch1Genuinely upw ind algorithms for th e multidimensionalEuler equations1AIA A J,1992,30(1):56~6334Parpia I H,M ichalek D J1Grid-independent upw ind sch eme for multidimensi onal flow1AIA A J,1993,31(4)135Hartwich P M1Comparison of coordinate-i n variant and coordinate-aligned u pw ind for the Euler equations1A IAA J,1994,32(9): 1791~179936van Leer B1Advancing the accuracy and efficiency of explicit Euler solvers1AIAA90-001237Zhang X D,Trepanier J-Y,Reggio M,et al1Grid i nfluence on upw ind schemes for Euler and Navier-Stokes equations1A IAA J, 1996,34(4):717~72738Abarbanel S,Duth P,Gottlieb D1Splitting methods for low M ach number Euler and Navier-Stokes equations1Comp uters&Fluids, 1989,17(1):1~12(1996年11月29日收到第1稿,1997年3月30日收到修改稿)。
材料力学第20讲 Chapter7-4第七章 强度理论
低碳钢圆截面试件，实验表明： 在单向拉伸时会发生显著的屈服现象。
若在圆试件中部切出一个环形槽（如 图a所示）。 试 验表明：直到拉断都看不到显著的 屈服现象和塑性变形，而是在最弱部 位发生脆断。其断口平齐，与铸铁拉 伸断口相似（b）。 这是因为在最弱截面处，材料处于三向拉伸应力状态，斜截面 上的剪应力较小，不可能出现屈服现象，只可能发生脆断。
只要微元内的最大拉应力 1 达到了单向拉伸
的强度极限 b ，就发生断裂破坏。
脆性断裂的判据（或极限条件） 1 u
强度条件 1
19
《评价》
二向时：当 1 2 0 该理论与实验基本一致
三向时：当 1230同上
当主应力中有压应力时，只要 3 1 同上
当主应力中有压应力时，只要 3 1 误差较大
理论与实验基本符合 比第三理论更接近实际
29
二、相当应力（强度准则的统一形式）
r [ ] r —相当应力(equivalent stress)
r1 1
r21(23)
r3 13
r 4 1 2 [1 22 2 3 2 3 1 2 ]
[]1n{b,0.2,s}
30
强度理论应用于许用拉应力和许用切应力间的换算
m
在平均应力作用下，单元体的形
m
状不变, 仅发生是体积改变
m
7
按迭加原理（应力）
1
m
1-m
m
2
3
m
2-m 3-m
交互项
体积改变能密度
v v
1 2
3
v i
v i
i 1
3 2
mm
形状改变能密度 (畸变比能)
v d
1 2
岩体本构关系与强度理论
18
注意事项
通过以上教学步骤和注意事项 的安排，相信学生能够更好地 掌握表格数据的图表化技能， 为未来的学习和工作打下坚实
的基础
··· ···
-
谢谢您的 ·观看·
BUSINESS TRIP PROJECT PLAN
汇报人：XXX
汇报时间：XXXXX
了解图表的基本 类型及适用场景
课程目标
Stage 3 学会根据需：求 调整图表的样式 和布局
Stage 5
了解如何根据数 据变化更新图表
Stage 2
掌握在Ex：cel 中创建图表的方 法
Stage 4
掌握数据标签和 图例的添加方法
3 课程内容
1. 图表基础知 识
什么是图表？
图表的基本类型 及特点
岩体本构关系与强度 理论
1 课程简介 3 课程内容 5 教学步骤
-
2 课程目标 4 教学方法 6 注意事项
1 课程简介
课程简介
本课程旨在介绍如何使 用Excel软件进行表格
数据的图表化
通过学习这个课程，你 将掌握如何将表格数据 转化为直观的图表形式， 以更好地分析和理解数 据
2 课程目标
Stage 1ຫໍສະໝຸດ 5. 高级图表的 制作复合图表的制作
自定义图表的制作
数据透视表的创 建和应用
课程内容
4 教学方法
教学方法
5 教学步骤
作业布置：布置与课程内容 相关的作业，让学生进行实 际操作，加深对知识点的理
解和掌握
总结与回顾：对本节课的内 容进行总结与回顾，帮助学
生巩固所学知识
教学步骤
材料力学强度理论
x=24MPa
x
x=28MPa
解： 首先计算危险点处的主应力。 已知 x =28MPa、 y=0、 x= – 24MPa
由主应力计算式得
1 3
28 2
另一主应力
28 2
2
242
2 0
41.8
MPa
13.8
按莫尔强度条件，得
rM
1
t c
3
41.8
50 150
13.8
讨论
1、双剪应力强度理论与大多数金属材料的实验结果符 合得较好,对于铝合金在复杂应力状态下的实验结果，
较第四强度理论更为接近。 2、该理论也适用于岩石及土壤等材料，并与实验结果有
良好的符合。注意，其失效状态不再是屈服，而是 剪断或滑移。 3、该理论可看作是宏观固体力学中引用微观晶体滑移 理论而提出的一种进似，
A
B
C
0
从理论上讲，确定了 ABC 包络线，就可以确定各种应力状态 下的极限应力圆。
不同材料，包络线不同；
同一材料，包络线唯一。
3、莫尔强度理论的简化
以单向拉、压试验数据得两个极限应力圆，该两圆 的公切线代 替包络线，再除以安全系数。
强度条件：
压
拉
1
t c
3
t
0
相当应力：
3 c
1 t
rM
120 MPa r4
1 2
1
2 2
2
3 2
3
1
2
(a)
1 2
0
1202
120
1202
120
02
120MPa
（2）对于图 b 所示的单元体， 已知 1=14 0MPa，2=110MPa , 3=0
岩石力学与工程岩石本构关系与强度理论
其蠕变曲线和弹性后效曲线，如图3-15所示。
蠕变曲线
0 k2
0 k1
弹性后效
0
t1
t
图3-15 广义开尔文体蠕变曲线和卸载曲线
2020/11/3
34
3.4.4.5 饱依丁-汤姆逊体（PTh：H/M）
一、力学模型 k1,1
k2 , 2
图3-16 饱依丁-汤姆逊体力学模型
二、本构方程
本模型是由马克斯威尔体与虎克体并联而成，由并 联规则：
2020/11/3
19
3.4.4 组合流变模型
三种基本元件进行组合时应力、应变的计算规则：
1.串联组合体中各元件的应力相等；应变等于各元件应 变之和。 2.并联组合体中各元件的应变相等；应力等于各元件应 力之和。
5.4.4.1 圣维南体（St.V:H-C）
一、力学模型
图3-5 圣维南体力学模型
2020/11/3
y
x
21
xy E xy
2020/11/3
4
4.边界条件
（1）位移边界条件
us us，vs vs
（2）应力边界条件
l x m yx s f x s
m y
l xy
s
f y s
（3）混合边界条件
（在 su上）
（在 s 上）
2020/11/3
5
3.4 岩石流变理论
3.4.1概念
四、卸载方程
0
k
kt
1 e
在t t1 时卸载，即 0，代入本构方程：
k 0
2020/11/3
28
解上述微分方程可得：
kt
A1e
当t t1 时， 1 ，结合蠕变方程，可得卸载方程 ：
岩体本构关系与强度理论
对莫尔强度理论的评价：
优点：①适用于塑性岩石，也适用于脆性岩石的剪切破坏； ②较好解释了岩石抗拉强度远远低于抗压强度特征； ③解释了三向等拉时破坏，三向等压时不破坏现象； ④简单、方便:同时考虑拉、压、剪,可判断破坏方向. 不足：①忽视了σ2 的作用，误差：±10％； ②没有考虑结构面的影响； ③不适用于拉断破坏； ④不适用于膨胀、蠕变破坏。
判据的表达式 剪应力表达式:
4 ) t( t
2
主应力表达式:
2 ( ) 1 3 8 t 1 3 3 t
3 0 1 3 3 0 1 3
由格里菲斯判据得
1.当σ３＝０时，σ１＝σｃ＝８σｔ，即 c t =８，与库仑-纳维尔判据接近。 2. 适用脆性岩石拉破坏。 修正的格里菲斯判据（考虑摩擦效应f,C）
b.最大剪应力理论破坏面上剪应力最大； 而岩石破坏面上剪应力不是最大。
④歪形能理论
只与σ1 、σ2 和σ3三者之间的差的绝对值有关；
而与应力大小无关，这与岩石破坏现象不符。
1 2 2 2 [ ( ) ( ) ( )] 1 2 2 3 3 1 2
古典强度理论与岩石强度表现不符：
①最大拉应力理论没有考虑σ
2
和σ
3
的影响。
②最大伸长线应变理论虽考虑σ2 和σ3 的影响，但多向拉 比单向拉安全，与事实矛盾。
③最大剪应力理论与岩石试验结果不符
σ1－σ3≤[σ]
a.最大剪应力理论破坏面与σ1 的夹角为45°；
而岩石破坏面与σ1 的夹角为45°－φ/2。
第七章 岩体本构关系与强度理论
§ 7.1 概述 § 7.2 岩石的本构关系
第7章岩体本构关系与强度理论
σ σc
σ
利用图7-10中的关系，有：
σ3
1 2
(1 3)
1 2
(1
3)
ctg 2
sin 2
1.双向压7 缩应4力2圆，2.双向拉压应力圆，
3..双向拉伸应力圆 图7-10 二次抛物型强度包络线
其中：
n( t )
1 3 2
sin 2
(
1 3 )2 2
2
(
1
3
)
2
规定：
1、σ1为最大主应力 、σ2 为中间主应力、 σ3 为最小主应力 ；
2、压应力为正，拉应力为负，剪应力以逆时 针为正。位移与应变的规定也一样。
二、 岩石弹性本构关系 1．平面弹性本构关系
据广义虎克定理有：
成E/（1- μ 2） ，μ换成μ/（1- μ）。
2. 空间问题弹性本构方程
x
1 E
x
( y
z )
y
1 E
y
( z
x )
z
1 E
z
( x
y )
yz
2(1 E
) yz , zx
1
1 f f2
2
f
f
)
σ1
1 tan2 c
1 3tg 2 (45 / 2) 2ctg(45 / 2)
σc
arc( tan2 θ)
本构与强度理论教学课件
强度理论的应用
金属材料
适用于描述金属材料的强度行为 ，用于设计金属结构、评估金属
材料的承载能力和安全性。
复合材料
适用于描述复合材料的强度行为 ，用于设计复合结构、评估复合
材料的承载能力和安全性。
工程应用
强度理论广泛应用于土木工程、 航空航天、机械工程等领域，用 于评估结构的承载能力和安全性 ，优化结构设计，提高工程安全
02
强度理论是材料力学和断裂力学 中的重要概念，用于评估材料在 不同受力状态下的承载能力和安 全性。
强度理论的分类
01
02
03
最大剪应力理论
认为材料在最大剪应力达 到某一临界值时发生屈服 或断裂。
最大应变理论
认为材料在最大应变达到 某一临界值时发生屈服或 断裂。
米塞斯理论
认为材料在等效应力达到 某一临界值时发生屈服或 断裂。
实际应用的发展趋势
1 2
拓展应用领域
随着新材料和新技术的不断涌现，本构与强度理 论的应用领域将不断拓展，涉及航空航天、能源 、环保等领域。
提高工程安全性和可靠性
通过本构与强度理论研究，提高工程结构和设备 的承载能力和安全性，降低事故风险。
3
智能化和自动化技术的应用
结合智能化和自动化技术，实现本构与强度理论 的在线监测、实时分析和预警等功能。
本构与强度理论教学 课件
目录
• 本构理论概述 • 强度理论概述 • 本构理论与强度理论的联系与区别 • 本构与强度理论的实际应用 • 本构与强度理论的发展趋势与展望
01
本构理论概述
本构关系的定义
总结词
本构关系是描述材料在受力状态下性质如何变化的关系。
详细描述
本构关系是指材料在受到外力作用时，其内部性质如何随外力变化而变化的关 系。这种关系反映了材料的力学特性，是材料科学和工程领域中非常重要的概 念。
混凝土的本构关系
以主应力和主应变表示
则为:
式中切线弹性模量 和 ，泊松比 随应力状态和数值的变 化按下述方法确定。
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Darwin-Pecknold 本构模型
材料在双轴受压
应变为：
• 等效单轴应力-应变关系
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Darwin-Pecknold 本构模型
2、混凝土非线弹性本构模型____Ottosen本构模型
定义一非线性指标 ，表示当前应力状态
至混凝土
破坏（包络面）的距离，也即塑性变形发展的程度。假定
保持不变，压应力 增大至 时混凝土破坏，则
混凝土的多轴应力应变关系采用Sargin的单轴受压方程，即
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Ottosen本构模型
式中参数以多轴应力状态的相应值代替：
代入得一元二次方程，解之得到割线模量：
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Ottosen本构模型
混凝土的泊松比很难从试验中精确测定。Ottosen本构模型取割 线泊松比 随 的变化如图，计算式为：
式中可取：
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Ottosen本构模型
单轴受压应力-应变
多轴应力-应变
Ottosen本构模型
泊松比
§7.1.4 混凝土的本构关系
2、混凝土非线弹性本构模型____Ottosen本构模型 非线性指标
• 根据非线性指标 的定义， 值计算要通过破坏包络
面先求 ，在一般情况下需要经过多次迭代方能求出；