第八章地下洞室围岩稳定性分析第一节概述1.地下洞室(underground cavity):指人工开挖或天然存在于岩土体中作为各种用途的构筑物。
2.我国古代的采矿巷道,埋深60m,距今约3000年左右(西周)。
目前,地下洞室的最大埋深已达2500m,跨度已过50m,同时还出现有群洞。
3.分类:按作用分类:交通隧洞(道)、水工隧洞、矿山巷道、地下厂房仓库、地铁等等;按内壁有无水压力:有压洞室和无压洞室;按断面形状为:圆形、矩形或门洞形和马蹄形洞室等;按洞轴线与水平面间的关系分为:水平洞室、竖井和倾斜洞室三类;按介质,土洞和岩洞。
4.地下洞室→引发的岩体力学问题过程:地下开挖→天然应力失衡,应力重分布→洞室围岩变形和破坏→洞室的稳定性问题→初砌支护:围岩压力、围岩抗力(有内压时)(洞室的稳定性问题主要研究围岩重分布应力与围岩强度间的相对关系)第二节围岩重分布应力计算1.围岩:指由于人工开挖使岩体的应力状态发生了变化,而这部分被改变了应力状态的岩体。
2.地下洞室围岩应力计算问题可归纳的三个方面:①开挖前岩体天然应力状态(一次应力、初始应力和地应力)的确定;②开挖后围岩重分布应力(二次应力)的计算;③支护衬砌后围岩应力状态的改善。
3.围岩的重分布应力状态(二次应力状态):指经开挖后岩体在无支护条件下,岩体经应力调整后的应力状态。
一、无压洞室围岩重分布应力计算1.弹性围岩重分布应力坚硬致密的块状岩体,当天然应力()c v h σσσ21≤、,地下洞室开挖后围岩将呈弹性变形状态。
这类围岩可近似视为各向同性、连续、均质的线弹性体,其围岩重分布应力可用弹性力学方法计算。
重点讨论圆形洞室。
(1)圆形洞室深埋于弹性岩体中的水平圆形洞室,可以用柯西求解,看作平面应变问题处理。
无限大弹性薄板,沿X 方向的外力为P ,半径为R 0的小圆孔,如图8.1所示。
任取一点M (r ,θ)按平面问题处理,不计体力。
则:……………………①式中Φ为应力函数,它是x 和y 的函数,也是r 和θ的函数。
边界条件:()()()()()⎪⎪⎪⎩⎪⎪⎪⎨⎧===>>-=⎪⎭⎫ ⎝⎛--=>>+=-++=====003103131R b 0)(2sin 22sin 2)(2cos 222cos 22b r r b r r br r b r r R b p R b p p θθτσθθσστθθσσσσσ ………………② 设满足方程①的应力函数φ为:()θ2cos ln 222F Dr cr Br r A ++++=Φ-………………………………③由③代入①,并由②可得:2R F ,4-D ,4-c ,4B ,2204020p pR pp pR A ====-= ⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂Φ∂-∂Φ∂=∂Φ∂=∂Φ∂+∂Φ∂=θθτσθσθθr r r r r r r r r 22222221111图 8.1柯西课题分析示意图故应力函数φ为:⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛-----=Φθ2cos 2212ln 222020220220r R R r R r r pR …………………………④ 由④代入①可得各应力分量:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-+-=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛+-⎪⎪⎭⎫ ⎝⎛+=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛+-+⎪⎪⎭⎫ ⎝⎛-=θτθσθσθθ2cos 32122cos 31122cos 34112440220440220440220220r R r R p r R r R p r R r R r R p r r ………………………………⑤ 岩体天然应力比值系数为λ,假定岩体为无重板的力学模型如图8.2所示。
若水平和铅直的天然应力均为主应力,则开挖前板内的天然应力为:⎪⎩⎪⎨⎧=====0zx xzv h v vz ττλσσσσσ (1)由铅直天然应力v σ引起(产生)的重分布应力:θπθθσ→='→2,-p v代入⑤得:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-+=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛+--⎪⎪⎭⎫ ⎝⎛-=θστθσσθσσθθ2sin r 3r 122cos r 31r 122cos r 3r 41r 12440220r 440220*********r R R R R R R R v v v ……………………………⑥图8.2 圆形洞室围岩应力分析模型(2)由水平天然应力h σ产生的重分布应力:p v h →=λσσ代入⑤得:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-+-=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++⎪⎪⎭⎫ ⎝⎛+=⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫⎝⎛+-+⎪⎪⎭⎫ ⎝⎛-=θστθσσθσσθθ2sin r 3r 2122cos r 31r 122cos r 3r 41r 12440220r 440220*********R R R R R R R h h h v ………………………………⑦ ⑥+⑦得由v σ和h σ同时作用时引起圆形洞室围岩重分布应力的计算公式:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛-+--=⎪⎪⎭⎫⎝⎛+--⎪⎪⎭⎫ ⎝⎛++=⎪⎪⎭⎫ ⎝⎛+--+⎪⎪⎭⎫ ⎝⎛-+=θσστθσσσσσθσσσσσθθ2sin 32122cos 312122cos 34121244022********40220220r R r R r R r R rR r R r R v h r v h v h v h v h v …………………⑧由⑧式可知:当v σ、h σ和R 0恒定时,重分布应力是研究点位置()θ,r 的函数。
当0R r =时,洞壁上的重分布应力:()⎪⎩⎪⎨⎧=--+==02cos 20θθτθσσσσσσr v h v h r ………………………………………………⑨ ※重分布应力的影响范围分析:图 8.3 V σσθ/随λ的变化曲线设0v h : ,1σσσλ===即,由⑧式可得:⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎭⎫ ⎝⎛+=⎪⎪⎭⎫ ⎝⎛-=0112202200θθτσσσσr r r R r R …………………………………………………………………⑩图8.4 r σ,θσ与r 之间的变化曲线(1)当0r 02 ,0 ,σσσθ===时R r ,洞壁上应力差最大,且处于单向受力状态,最易发生破坏;(2)00r , ,σσσσθ↓→↑→∞→r 那么? , ?0r 呢即时θθσσσσσ=→=r r 计算得出%)8.2( ,60θσσ==r R r 时因此,一般认为,地下洞室开挖引起的围岩分布应力范围为6R 0。
在此范围之外,不受开挖影响。
(2)其他形状洞室 非圆形由前面可知,重分布应力的最大值在洞壁上,且仅有θσ,只要在θσ作用下洞壁围岩不发生破坏,那么洞室围岩一般就是稳定的。
引入“应力集中系数”→为了研究各种洞形洞壁上的重分布应力及其变化情况应力集中系数:地下洞室开挖后洞壁上一点的应力与开挖前洞壁处该点天然应力的比值。
如圆形洞室洞壁处的应力:()()()vh v h v h σθσθθσσσσσθ2cos 212cos 21 2cos 2++-=--+=令θβθα2cos 21,2cos 21+=-=,则:v h βσασσθ+=其中α、β为应力集中系数,其大小仅与点的位置有关。
同理,根据光弹实验或弹性力学方法可求得不同形状洞室的应力集中系数α和β,见教材P 144。
(3)软弱结构面对围岩重分布应力的影响假定岩体中结构面无抗拉能力,且f τ很小,在剪切过程中,结构面无剪胀作用。
(1)软弱结构面v σ⊥,沿水平直线方向0 ,0==θτθr 时θσ、r σ均为主应力,结构面上无剪切应力,无影响。
不会沿结构面产生滑动,结构面对围岩重分布应力的弹性分析无影响。
图 8.5 沿圆形洞水平轴向方向发育结构面的情况及应力分析示意图(2)软弱结构面平行v σ,沿铅直方向0 ,90=︒=θτθr 时,不会影响围岩弹性应力分布。
但是,当31<λ时,顶底板产生拉应力。
λλ2310-=∆R h结构面被拉开,形成应力降低区,有影响。
2.塑性围岩重分布应力岩体受结构面切割使其整体性丧失,强度降低,在重分布应力作用下,很容易产生塑性变形。
s σσθ≥时,洞壁围岩由弹性→塑性状态→形成一个塑性松动圈。
r →↑,围岩中出现三个区:图8.6 软弱结构面对重分布应力的影响示意图图 8.7 围岩中出现塑性圈时的应力状态塑性圈(裂隙增多,C 、φ和E ↓)→弹性圈(θσ↑)→天然应力区(r=6R 0) (“单向应力状态→双向应力状态”) 此时,就必须采用弹塑性理论求解。
如何求解塑性圈内的重分布应力? 假设①岩体均质、各向同性,连续;②开挖洞室半径为R 0,塑性圈半径为R 1,岩体中的天然应力o v h σσσ==,(11=-=μμλ时情况);③圈内岩体强度服从Mohr 直线强度条件即(Coulomb -Naiver 准则)。
分析如图8.8所示:图 8.8 塑性圈围岩应力分析图1=λ, 0==∴r r θθττ由0=∑r F (取向外为正,向内为负)得:()()02sin 2=⎪⎭⎫⎝⎛+++-θσθσσθσθd dr d dr r d rd r r r 22sin ,θθd d dr =⎪⎭⎫ ⎝⎛很小 上式整理变为:()r rr rd rd dr σσσσθ==-…………………………………………………………①塑性圈内的θσ和r σ是主应力,设岩体满足如下塑性条件(Mohr 准则):mmm m r m m ctg C ctg C φφφσφσθsin 1sin 1-+=++………………………………………………………….②由①式可得:r rrd rd σσσθ+=代入②积分后得: ()A r ctg C mmm m r +-=+ln sin 1sin 2ln φφφσ…………………………………………….③A —常数边界条件:i r P R r ==σ ,0(i P 为洞室内壁上的支护力)代入③可得:()0ln sin 1sin 2ln R ctg C P A mmm m i φφφ--+=……………………………④④代入③可得:()m m mmm m i r ctg C R r ctg C P φφφφσ--⎪⎪⎭⎫⎝⎛+=sin 1sin 20 同理求得环向应力θσ:()m m mmm m i ctg C R r ctg C P mm φφφφσφφθ-⎪⎪⎭⎫ ⎝⎛-++=-sin 1sin 20sin 1sin 1塑性圈内围岩重分布应力的计算公式:()()⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧=-⎪⎪⎭⎫ ⎝⎛-++=-⎪⎪⎭⎫ ⎝⎛+=--0sin 1sin 1sin 1sin 20sin 1sin 20θφφθφφτφφφφσφφσr m m m m m m i m m m m i rctg C R r ctg C P ctg C R r ctg C P m mmm…………………………⑤ 式中:m C 、m φ为塑性圈岩体的内聚力和摩擦角;r 为向径; i P 为洞壁支护力;0R 为洞半径。