地下硐室围岩应力计算及稳定性分析说课讲解.ppt

地下硐室围岩应力计算及稳定性分析,本章内容,3-1概述3-2弹性理论计算坑道围岩与衬砌应力3-3坑道围岩应力分布的弹塑性力学分析法3-4坑道围岩位移3-5围岩压力计算3-6坑道支护3-7竖井围岩应力计算及稳定性分析,本章的重点难点：,1、圆形坑道围岩应力弹塑性理论分析方法；2、围岩与支护相互作用原理；3、弹塑性理论计算围岩压力4、块体平衡理论计算围岩压力；5、压力拱理论计算围岩压力；6、太沙基理论计算围岩压力；7、喷锚支护的力学作用；8、圆形竖井围岩应力分布与稳定性评价。,要求,1、掌握本课程重点难点内容；2、掌握圆形坑道围岩应力分布规律；3、了解椭圆形、矩形坑道周边应力分布；4、掌握有内压圆形坑道围岩与衬砌的应力计算5、了解塑性区半径、松弛区半径及围岩位移的计算公式；6、了解岩体构造对井壁稳定性的影响；7、掌握井壁压力的平面挡土墙计算方法；8、了解井壁压力空心圆柱体挡土墙计算方法。,3-1概述,一、地下硐室的分类地下硐室(undergroundcavity)是指人工开挖或天然存在于岩土体中作为各种用途的构筑物。按用途：矿山巷道（井）、交通隧道、水工隧道、地下厂房（仓库）、地下军事工程按硐壁受压情况：有压硐室、无压硐室按断面形状：圆形、矩形、城门洞形、椭圆形按与水平面关系：水平硐室、斜硐、垂直硐室（井）按介质类型：岩石硐室、土硐按应力情况：单式硐室、群硐,二、地下硐室围岩应力分析方法,块状结构岩体：块体平衡理论分析碎裂和松散结构岩体：松散体力学分析,各向同性岩体,各向异性岩体,3-1概述,完整结构的岩体：弹塑性力学分析,普氏压力拱理论,太沙基理论,根据围岩的结构不同，可采用不同的分析方法。,3-2弹性理论计算坑道围岩与衬砌应力,基本假定：岩体为均质、连续和各向同性的介质。,将巷道和围岩视为无重量的有孔平板的平面应变问题，平板所受到的外力即原岩应力。巷道上部和下部的初始应力不相等，但当巷道埋深大于其高度的20倍时，这种应力差即可略去。于是，当p=q,即1，可视为二向等压下有孔平板平面应变问题，当pq时，即1，则视为二向不等压的有孔平板平面应变问题。计算结果表明，采用这种计算误差不超过1。,研究围岩二次应力状态的方法：,一、无内压坑道围岩应力分布1、圆形坑道围岩应力分布,设原岩垂直应力为p，水平应力为q，作用在围岩边界，忽略围岩自重的影响，按弹性理论中的基尔希公式计算围岩中任一点M（r,）的应力：,（1）当r时，,(3-1),(3-2),上式即为极坐标中的原岩应力。,（2）当ra时，即坑道周边的应力为：,(3-3),或：,式中：=q/p为侧压力系数。,(3-4),(3-1),由：,可见，与和密切相关。当0，时，,当3/2，/2时，,由于岩体的抗拉强度很小，认为岩体不抗拉，因此，坑道周边不能出现拉应力的条件为：,解得：,当0，时，,当3/2，/2时，,不同的下，坑道周边切向应力的分布：,不同的下，坑道周边切向应力的分布：,不同的下，坑道周边切向应力的分布：,不同的下，坑道周边切向应力的分布：,（3）当p=q,即=1时，,(3-1),(3-5),可见，、r与无关，=1（轴对称）时对圆形坑道围岩应力分布最有利。,当r=a,坑道周边应力为：,(3-6),圆形坑道开挖应力扰动范围为坑道半径的35倍。,当r时，坑道原岩应力为：,(3-7),几何模型,物理模型,圆形巷道周围sigmaYY等值线图,圆形巷道右侧sigmaYY变化曲线,圆形巷道周围竖直方向的位移等值线图,圆形巷道右侧竖直方向位移变化曲线,2、椭圆形坑道周边应力分布,在单向应力p0作用下，椭圆形坑道周边任一点的径向应力r、切向应力、剪应力r，根据弹性力学计算公式为：,式中：my轴上的半轴b与x轴上的半轴a的比值，即mb/a；洞壁上任意一点M与椭圆形中心的连线与x轴的夹角；荷载p0作用线与x轴的夹角；p0外荷载。,若0，p0p,则：,若900，p0p,则：,在原岩应力p、p作用下,则由（1）（2）得：,(1),(2),(3-8),上式也可表示为：,坑道周边两帮中点处（0,)切向应力为：,若（a)=(b),即12，则可得：,坑道周边顶底板中点处（3/2,/2)切向应力为：,(3-8),(a),(b),(c),由（c)可得：,可见，在原岩应力(p,p）一定的条件下，随轴比m而变化。为了获得合理的应力分布，可通过调整轴比m来实现。,(c),满足上式的轴比叫等应力轴比。在等应力轴比的条件下，椭圆形坑道顶底板中点和两帮中点的切向应力相等，周边应力分布比较均匀。,(3-9),例：1/4条件下，不同轴比m对应的顶底板和两帮中点处的：,(1)当m1，顶底板中的出现拉应力，故在1/4条件下，应选m1.(2)当m4时，巷道两帮中点和顶底板中点的应力为1.25p，出现切向应力相等的应力状态，即等应力轴比状态。,在等应力轴比状态下，即,将上式代入（3-8）：,在等应力轴比条件下，与无关，周边切向应力为均匀分布。,可见，椭圆形长轴与原岩最大主应力方向一致时，坑道周边不出现切向拉应力，应力分布较合理，等应力轴比时最好。,3、矩形坑道围岩应力分布,由实验和理论分析可知，矩形巷道围岩应力的大小与矩形形状（高宽比）和原岩应力（）有关。,高宽比1/3,1,矩形坑道围岩应力分布特征：（1）顶底板中点水平应力在坑道周边出现拉应力，越往围岩内部，应力逐渐由拉应力转化为压应力，并趋于原岩应力q；（2）顶底板中点垂直应力在坑道周边为0，越往围岩内部，应力越大，并趋于原岩应力p；（3）两帮中点水平应力在坑道周边为0，越往围岩内部，应力越大，并趋于原岩应力q.,（4）两帮中点垂直应力在坑道周边最大，越往围岩内部，应力逐渐减小，并趋于原岩应力p；,高宽比1/3,1,（5）巷道四角处应力集中最大，其大小与曲率半径有关。曲率半径越小，应力集中越大，在角隅处可达68。,例：不同和不同轴比m下，矩形坑道周边顶底板和两帮中点处的：,矩形坑道断面长轴与原岩最大主应力方向一致时，围岩应力分布较合理，等应力轴比时最好。,4、坑道围岩分布的共同特点：,（1）无论坑道断面形状如何，周边附近应力集中系数最大，远离周边，应力集中程度逐渐减小，在距巷道中心为35倍坑道半径处，围岩应力趋近于与原岩应力相等。（2）坑道围岩应力受侧应力系数、坑道断面轴比的影响，一般说来，坑道断面长轴平行于原岩最大主应力方向时，能获得较好的围岩应力分布；而当坑道断面长轴与短轴之比等于长轴方向原岩最大主应力与短轴方向原岩应力之比时，坑道围岩应力分布最理想。这时在巷道顶底板中点和两帮中点处切向应力相等，并且不出现拉应力。,（3）坑道断面形状影响围岩应力分布的均匀性。通常平直边容易出现拉应力，转角处产生较大剪应力集中，都不利于坑道的稳定。（4）坑道影响区随坑道半径的增大而增大，相应地应力集中区也随坑道半径增大而增大。如果应力很高，在周边附近应力超过岩体承载能力而产生的破裂区半径也将较大。（5）上述特征都是在假定坑道周边围岩完整的情况下才具备的。在采用爆破方法开挖的坑道中，由于爆破的松动和破坏作用，坑道周边往往不是应力集中区，而是应力降低区，此区域又叫爆破松动区。该区域的范围一般在05m左右。,4、坑道围岩分布的共同特点：,二、有内压坑道围岩与衬砌的应力计算,1、内压引起的围岩附加应力（1）厚壁筒应力公式设一弹性厚壁筒，内径为ri,外径为R,内压为pi，外压为pa,由弹性理论拉密解答，在距中心为r处的径向应力和切向应力为：,(3-9)厚壁筒应力公式,（2）水工隧道中内压引起的围岩附加应力,将隧道围岩看成厚壁筒，内径为ri=a,外径为R=,隧道充水后所产生的内压为pi，外压为pa0,由弹性理论拉密解答：,得出在距中心为r处的径向应力和切向应力为：,(3-10),(3-9)厚壁筒应力公式,在r=a（洞周边）：,在距中心为r处的径向应力和切向应力为：,(3-10),上式即是内压pi引起的附加应力。,（3）原岩应力为p（=1)、水工隧道中内压为pi时的围岩应力：,2、有内压坑道围岩与衬砌的应力计算,（1）无裂隙围岩A、刚度系数法求衬砌的应力a衬砌外周边的径向位移,设混凝土衬砌坑道的内径为ri,外径为a,围岩对衬砌的压力pa，内压为pi,混凝土的弹性模量和泊松比分别为Ec和c,混凝土衬砌内距坑道中心为r处的径向位移为u,由弹性理论有：,将拉密公式代入得：,当r=a时，即得衬砌外周边的位移,式中：t=a/ri,2、有内压坑道围岩与衬砌的应力计算,当r=a时，衬砌外周边的径向位移：,b坑道周边围岩的变形设刚度系数为k，坑道周边围岩在压力pa作用下发生的变形：,c根据变形协调条件，坑道周边围岩变形与衬砌变形相等，即式（3-12)=(3-11),则有：,(3-11),(3-12),即：,(3-13),令pa/pi=k1,则pa=k1pi,将pa、pi代入厚壁筒公式得到混凝土衬砌内距坑道中心为r处的应力为,由于是平面应变问题，故轴向应力为：,(3-14),(3-15),(3-16),(3-13),B、内压分配法求围岩应力,设内压pi通过衬砌传递到围岩上的压力为pa，papi，为内压分配系数。假设衬砌与围岩紧密接触。,设围岩的弹性模量为E,泊松比，由弹性力学得围岩内半径为r处的径向应变为：,在r=a处，即坑道壁面：r=pa,=-pa,(3-17),对u积分，并令r=a得坑道壁面围岩位移：,由（3-11）衬砌外周边处径向位移：,(3-11),坑道壁面围岩位移：,(3-17),式（3-11）式（3-17），于是：,（3-18),B、内压分配法求围岩应力,（3-18),求出后，即可按（3-10）求出围岩任一点由内压引起的附加应力，按厚壁筒公式（3-9）求出衬砌内任一点的应力。,(3-10),例题：P187,（2）有裂隙围岩,设围岩有径向裂隙，其深度为d，沿岩石表面的径向压力可假定为：,(3-19),(3-20),在裂隙岩体任一深度处（rd):,(3-21),在裂隙岩体外边界处（rd)，压力为:,(3-22),（2）有裂隙围岩,(3-23),在围岩内任一点(d2），而坑道顶部岩体不稳定，会发生冒落而形成自然平衡拱。,模型2：假定坑道两帮岩体也不稳定（f2），发生剪切破坏，导致平衡拱的跨度扩大。,（1）模型1及相应的计算方法,A、平衡拱形状平衡拱形状如图所示，为一条抛物线，其方程为：,B、平衡拱跨度平衡拱跨度等于坑道跨度，即2a。C、平衡拱高度b:,E、单位长度平衡拱内岩体重量W：,F、单位长度坑道上作用在坑道支架上的顶压pv:,D、平衡拱面积：,（2）模型2及相应的计算方法,当fz，即,则井筒的最大主应力为切向应力，井筒破坏条件：,竖井井筒周边（r=a)：,式中c为岩体的单轴抗压强度，hcr为临界深度。,当q2，则在井筒两壁产生应力集中，井筒两壁可能发生剪切破坏，如图所示。,如果q2/q13,井壁B、D点则产生拉应力，可能产生拉破坏,2、岩体构造对井筒围岩稳定性的影响,（1）井筒附近断层对井筒围岩稳定性的影响,（2）斜截断层对井筒围岩稳定性的影响,（3）节理岩体中竖井的稳定性问题,在节理化岩体中，采空区围岩移动方向和移动范围，往往受主节理控制，移动角接近于主节理倾角。,（三）竖井井壁压力计算,竖井井壁压力计算是设计井筒衬砌或支护的依据。1、表土层中井壁压力计算的依据表土层中井壁压力显现的形式：A、土体产生剪切破坏，形成空心圆锥形滑移体，作用在衬砌上一个水平分布力；B、粘土矿物含量高，遇水膨胀或冻胀，对井壁产生膨胀压力；C、流沙等流动性土体对井壁产生流动性压力。,（1）平面挡土墙计算法,土体产生剪切破坏，形成空心圆锥形滑移体。将井壁视为平面挡土墙，将土体或破碎岩体视为无粘结力的松散体，作用在井壁衬砌上的压力按主动土压力计算。,式中：pns,pnx为第n层上界点和下界点处井壁压力，hi为第i层土厚度，n为第n层土的内摩擦角。,（2）空心圆柱体挡土墙计算法,将井筒按空间结构物考虑，将下滑土体看成是空心圆锥体，采用空间轴对称极限平衡方程求解，其计算公式为：,式中：qn计算土层以上的土层传来的均布压力,Nn计算简化系数，Nn2tgntg(450+n/2);RnRn=a+htg(450-n/2),h计算土层上界点至计算点的高度；a井筒半径；Cn,n,n第n层岩层的粘结力、重度、内摩擦角。,（3）悬浮体计算法,将地下水位线以下的土体视为悬浮体，分别计算悬浮体和地下水对井筒衬砌的压力，从而对平面挡土墙和空心圆柱体挡土墙计算作了2点共同的修正。A、用悬浮体重度代替土体天然重度,式中：Gs土颗粒比重；V土体体积，V=VsVv;Vs，Vv分别为土体中土粒的体积和孔隙的体积。,B、地下水对井壁的压力,地下水对井壁的压力单独计算，等于计算点处的静水压力。如果井壁衬砌渗漏，可考虑一个折减系数(=0.80.9），即,式中：m地下水位以上土层数;w水的重度，w9.8kN/m3；n由地表至计算土层的土层数。,平面挡土墙公式：,根据以上2点的修正，平面挡土墙公式变为：,