围岩变形弹塑性分析

1、§ 2.1 道围岩重分布应力的计算隧道开挖前,岩体中每个质点均受到天然应力的作用而处于相对平衡状态；隧洞开挖后,洞壁岩体因失去了原有岩体的支撑,破坏了原有的平衡状态, 从而产生向洞内空间的膨胀变形,其结果又改变了相邻质点的相对平衡关系,引起应力、应变和能量的重新调整,到达新 的平衡关系,形成新的应力状态.2.1.1 弹性围岩重分布应力对于那些坚硬致密的块状岩体,当天然应力大约等于或小于其单轴抗压强度的一般时,隧道开挖后的围岩将呈弹性变形状态. 这类围岩可近似视为各向同性、连续、均质的线弹性 体,其围岩应力重分布可用弹性力学的根本理论来分析,隧洞半径相对于洞长很小时,可按平面应变问题考

2、虑,围岩重分布应力可用柯西 Kirsh课题求解.图2-1是柯西课题的简化模型.设无限大弹性薄板,在边界上受沿X方向的外力P作用, 薄板中有一半径为 R的圆形小孔.取如图极坐标,薄板中任一点 M ,.的应力及方向如图所示,按平面问题考虑,不计体力,那么M点的各应力分量,即径向应力 二、环向应力9和 剪应力与应力函数间的关系,根据弹性理论可表示为:11 1r .r r n(2-1 )=rp -产1加 1冽-2rr - r 1-r - t1上式的边界条件为:9rL = p 十1cos 26b 6K- 5sin 2日b L Rg > (2-2)b = R0设满足该方程的应力函数 中是:A = A

3、ln r +Br2+(Cr2+Dr/+F )cos 26(2-3)带入上式并考虑边界条件,可求得应力函数为：4=pR0- In r r 1 r -Ry cos 2日(2-4 )2 _2 苗 2苗 2r2代入可得各应力分量:p %3R44%r =(1 -3)+(1 + - - -4)C0S 2e i21r r r p 苗 3弓(2-5)=p (1 +：) _(1 + 一)cos 29921 r2r4p3R42R2r = -京1-4)s>n 22 r r式中,crx, ctq,分别为m点的径向应力、环向应力和剪应力,以压应力为正,拉应力为负；6为M点的极角,自水平轴x轴起始,反时针方向为正；

4、r为径向半径.式2-5是柯西课题求解的无限薄板中央孔周边应力计算公式,我们把它应用到隧道围岩重分布应力计算中来.假定隧道开挖在天然应力比值系数九的岩体中,那么问题可简化为如下图的无重板岩体力学模型.假设水平和垂直天然应力都是主应力,那么隧道开挖前板内的天然应力为：式中,5, Oh为岩体中垂直和水平天然应力；xz ,Ezx为天然剪应力.取垂直坐标轴为Z,水平轴为x,那么隧道开挖后,铅直天然应力引起的围岩重分布应力也可由式2-5 确定.在式2-5 中,p用.V代替,日是径向半径OM与x轴的夹角日'来表示,那么JTcos=cos 28 卜(2-7 )sin2u =一 sin 2这样由二v引起

5、的重分布应力为:(1)-(13R44R2)cos 2u仃日；1(13年)cos 212 (2-8)(13R42R2)sin 2ir水平应力h产生的重分布应力,可由式2-5直接求得:(1鸟(1 r3R：尸)cos1129 1.； v"T(1鸟一 (1 r3R4-4)cos 2ur卜2-93K2年)sin2u将以上两式联立求和,即可得到隧道弹性围岩重分布的弹性计算方程:中13R44M4- - 2")cos 2 71r rr -1 二- vR 1 - (102)-(1r22驾)cos 2.r2 (2-10)驾)sinr2u由上式可知当r = R0时,即围岩洞壁上的应力集中最大.2

6、.1.2塑性围岩重分布应力由弹性围岩重分布应力特点可知,隧道开挖后洞壁的应力集中最大.当洞壁重分布应力超过围岩屈服极限时,洞壁围岩就由弹性状态转化为塑性状态,并在围岩中形成一个塑性松动圈.这种塑性圈随着距洞壁距离的增大,径向应力由零逐渐增大, 应力状态由洞壁的单向应力状态逐渐转化为双向应力状态.莫尔应力圆由与强度包络线相切的状态逐渐内移,变为与强度包络线不相切,围岩的强度条件得到改善.围岩也就由塑性状态逐渐转化为弹性状 态,围岩中出现塑性圈和弹性圈.塑性圈的出现,使圈内一定范围内的应力因得到释放而明显降低,最大应力集中由原弹性区以外那么是原岩应力区,也来的洞壁移至塑、弹圈交界处,使弹性区的应力

7、明显升高.即应力根本没有发生变化的天然应力区,各局部的应力变化如图 9-8所示.此时围岩重分布应力应采用弹塑性理论求解.假设在均质、各向同性、连续的岩体中开挖一半径为 R隧道；开挖后形成的塑性圈半径为R ,塑性圈内岩体强度服从莫尔直线强度条件,塑性圈以外围岩体仍处于弹性状态.如图9-9所示,在塑性圈内取一微小小单元体abcd, bd面上作用有径向应力crr , ac面上的径向应力为 +dor,ab和cd面上作用有切向应力 仃$入=1,剪应力加6 = 0.当单元体处于极限平衡状态下时,作用在单元体上的全部力在径向应力在半径r上的投影为零,那么单元体上径向应力的平衡方程为：,、,、 八 .,di

8、-3rd 日(% + d%)( r +dr)dH + 2仃.sin( ) = 0(2-11)d) d)当d6很小时,sin( )电 o那么上式可简化为:(仃e - 仃r)dr = rdr (2-12)由摩尔强度理论,岩体的塑性条件为:二1 Cm COt m _ 1 . s|n m二 r - Cm COt m - 1 - Sin m(2-13)联立以上两式得：d(. ： Cmcot m),1 sin m ,、dr 2 sin m drm = (m - 1) 一 = m(2-14)二rCm COt m 1 一 s|nm r 1 - s|n m r2 sin m积分整理后得径向应力c,为：% =(P

9、i Cm cot .m)( )b,m .Cmcot m同理可求得:2sin %=(P+CmC0t 峭 NK-CmSf所以,塑性圈内围岩重分布应力的计算公式为:二 r =(Pi二-(Pi2 sin .+ Cm COt 售)(J)m 一 Cm COt 爆RD,2sin <h,1 + sin中mr e布+ Cm COt t)；m()a-Cm COt 着 >(2-15)1 - s1n *m RqCm,m为塑性圈岩体的内聚力和内摩擦角；Pi为洞壁支护力.塑性圈与弹性圈交界面(r = B)上的重分布应力,该面上弹性应力与塑性应力相等可得:二 R1二 4 r peLIR1R sin 虱二0Cm

10、cot m)2a(2-16)=；=rpe = R1 -sin m) - K COS=二0(1Sinm) - Cm COS m=0二(R2 一 R2)=二(R _ r)2 _(R _ -1)22R R sinm(二 0Cm cotm)D =D =ZR0R B2q.R0rpe,口脚,7rpe为r = 6处的径向应力,环向应力和剪应力；.0为岩体天然应力.§ 2.2 岩的变形隧道开挖后,岩体中形成一个自由变形空间,使原来处于挤压状态的围岩,由于失去了原有的应力状态而发生向洞内松胀变形；如果这种变形超过了围岩本身所能承受的水平,那么围岩就要发生破坏. 隧道围岩位移是表征和限制隧道工作状态的重

11、要参数,研究说明：当围岩应力超过围岩的极限强度时,围岩将立即发生破坏. 当围岩应力的量级介于岩体的极限强度和长期强度之间时,围岩需经瞬时的弹性变形及较长时期蠕动变形的开展方能到达最终的 破；当围岩应力的量级介于岩体的长期强度及蠕变临界应力之间时,围岩除发生瞬时的弹性变形外,还要经过一段时间的蠕动变形才能到达最终的稳定；当围岩应力小于岩体的蠕变临界应力时,围岩将于瞬时的弹性变形后立即稳定下来.2.2.1 围岩弹性位移计算在天然应力不大的情况下,围岩常处于弹性状态.这时洞壁围岩的位移可用弹性理论进 行计算.根据弹性理论,平面应变与位移的关系为：u：ru 1人 -I 2 (2-17)yri平面应变与

12、应力物理方程为:1,2,E=1- M - N1 +(2-18)匚me72,?目=三1 + N 入r H 匚 m,r H匚me联立以上两式得:曳=/E -图阿-九.+九户日Cr Eme -：+3包=/一1_.-心小历2-19 r r CQEme J将式2-10代入上式并进行积分,可求得在平面应变条件下的围岩位移为:.二vEme3K 4R2«"r cos 2% (1 + % )R2仃 h 一 巴r rEmeJ 2r22i _m 二h R42R2 .m r0r - -0- sin 2iEme2r r3R0 cos2. r3(2-20)小(1 + m ) I % - 0vEme:2

13、3R4 r frW sin 291r当r = 6时,可得洞壁上的弹性位移为:(1 -2)R二一-二h WZ -:v)cos 2u Eme2(1 - YR,- E( - h - - v)s i n 2-me(2-21 )式中N , V分别为围岩内任一点的径向位移和环向位移；Eme,Nm为岩体的弹性模量和泊松比；其余符号同前.2.2.2塑性位移计算当围岩中由于应力重分布作用形成塑性圈以后,这时围岩的塑性位移采用弹塑性理论来分析.由于开挖卸载形成塑性圈后,弹、塑性圈交界面上的径向应力增量(Aar)r = R1 和环向应力增量区仃a r R1为:二rr=R1 = - 01 -rR2RR. -0r r二二 5 - ：7021rcR2=Ri= '' 012rR2一 OR -二0 产=g R1 -仃0)* (2-22)R2aR 2 一 仃01 r代入式2-19 并积分得交界面上的径向位移为:1rR (2-23)Em、Gm为塑性圈岩体的变形模量和剪切模量,<rR为塑性圈作用于弹性圈的径向应力,由式2-16得:二rpe =:.-