第五章 2 围岩应力分析
- 格式:doc
- 大小:984.50 KB
- 文档页数:44
5 地下洞室的围岩应力与围岩压力5.1 地下洞室的围岩应力计算及应力分布5.1.1 概述在岩体中开挖地下洞室,必然会破坏原来岩体内相对平衡的应力状态,并在一定范围内引起岩体天然应力状态的重分布。
岩体的强度和变形特性是否适应重分布以后的应力状态,将直接影响地下建筑物的安全。
为了正确评价地下建筑的稳定性,除进行必要的地质分析外,对围岩应力分布特征的分析和计算,也是评价围岩稳定性所必须的环节。
洞室开挖后,周围的岩石在一般情况下(侧压力系数<3)必然会在半径方向上发生伸长变形,在切线方向上发生压缩变形,这就使原来径向上的压缩应力降低,切向上的压缩应力增高,而这种降低和增高的程度随着远离洞壁逐渐减弱,达到一定距离后基本无影响。
通常将应力的这种变化称为应力重分布(即原始的应力状态变化到新的平衡的应力状态的过程)。
把应力重分布影响范围内的岩体称为围岩。
围岩内的应力称为围岩应力或二次应力(相对与天然应力)。
理论研究和实际测量结果表明,围岩应力的分布规律与开挖前岩体的天然应力状态及洞型等有关。
地下工程在设计、施工和使用时,总是要研究其稳定性问题。
在地下工程(井巷、隧道、洞室等)工作期内,安全和所需最小断面得以保证,称为稳定。
稳定如果用公式来表示的话,就是:Uu S <<max max σ 其中,σmax 、u max ——地下工程岩体或支护体中最大、最危险的应力与位移;S 、U ——岩体或支护材料的强度极限与位移。
无论无支护或有支护,凡涉及这方面研究的问题,统称为稳定性问题。
地下工程稳定性可分为两类:(1)自稳——能长期自行稳定的情况,如天然石灰岩溶洞、某些金属采矿场等。
通常不需要进行支护。
(2)人工稳定——需要依靠支护才能达到稳定的情况,如煤矿中的软岩巷道、表土洞室等,由于次生应力场的作用形成破碎带。
地下工程自身影响范围达不到地面的,称为深埋,否则称为浅埋。
深埋地下工程存在如下力学特点:(1)可视为无限体中的孔洞问题,孔洞各方向的无穷远处仍为原岩体;(2)当埋深Z 达到巷道半径或宽高之半的20倍及以上时,巷道影响范围内的岩体自重可忽略不计;原岩水平应力可以简化为均匀分布,通常误差不大(在10%以下);(3)深埋的水平巷道长度较大时,可作为平面应变问题处理。
围岩应力应变计算公式引言。
在岩石工程中,围岩的应力应变计算是非常重要的一部分。
通过对围岩的应力应变进行准确的计算,可以帮助工程师更好地了解围岩的稳定性,从而有效地设计和施工工程项目。
本文将探讨围岩的应力应变计算公式,以及如何应用这些公式来进行工程实践。
围岩的应力应变计算公式。
围岩的应力应变计算公式可以通过岩石力学理论来推导和建立。
在岩石力学中,通常使用弹性模型来描述岩石的应力应变关系。
弹性模型假设岩石在受力作用下会产生弹性变形,即在去除外力后能够完全恢复到原来的形状。
基于弹性模型,可以得到围岩的应力应变计算公式如下:ε = σ / E。
其中,ε表示围岩的应变,σ表示围岩的应力,E表示围岩的弹性模量。
这个公式表明了围岩的应变与应力之间的线性关系,弹性模量E越大,围岩的应变就越小。
应力应变计算公式的应用。
围岩的应力应变计算公式可以应用于岩石工程的各个方面,包括岩石的稳定性分析、隧道和坑道的支护设计、岩石爆破工程等。
下面将分别介绍这些方面的应用。
1. 岩石的稳定性分析。
在岩石工程中,经常需要对围岩的稳定性进行分析,以确保工程项目的安全性。
围岩的应力应变计算公式可以帮助工程师计算出围岩在受力作用下的变形情况,从而评估围岩的稳定性。
通过对围岩的应力应变进行分析,工程师可以确定围岩的变形情况,为工程项目的设计和施工提供重要的参考依据。
2. 隧道和坑道的支护设计。
在隧道和坑道工程中,围岩的应力应变计算公式可以用于支护结构的设计。
通过计算围岩的应力应变,工程师可以确定支护结构的尺寸和材料,以确保支护结构能够有效地抵抗围岩的变形和破坏。
应力应变计算公式还可以用于评估支护结构的稳定性,为支护设计提供科学依据。
3. 岩石爆破工程。
在岩石爆破工程中,围岩的应力应变计算公式可以用于评估爆破后围岩的变形情况。
通过计算围岩的应力应变,工程师可以确定爆破参数,以减小围岩的应力集中,降低围岩的破坏程度,从而提高爆破效果和保证工程安全。
四、 围岩应力分析地下洞室的开挖,会产生应力重分布,围岩应力不仅与天然应力场有关,而且还与洞室的开挖有关。
基本假定:①岩体均质,各向同性,连续体。
②无构造应力作用,仅由自重应力形成天然应力场,其大小为z z ⋅=γσ取K 0=0,1/3,1三种情况下的应力场③忽略洞室高度上的应力场变化,即认为洞顶和洞底处的天然应力相同,则有H P v γ=(H 为洞中心的深度)④为平面应力问题K = 10P = P h vPv = r . Z1. 圆形洞室洞室开挖前的天然应力为H P v γ=v h P K P 0=围岩中的径向应力r σ,切向应力θσ以及剪应力θτr 可按下述公式计算:(Ⅰ)r0——洞室半径(m)r——自洞室中心算起的径向距离(m)θ——自水平轴算起的极坐标中的角度p v——垂直方向的压应力(MPa)(=γH)p h——水平方向的压应力(MPa)(= K0 p v)讨论:1).K 0=1; p h =p v (K 0=1静水压力式的天然应力场) 洞室周围处于等压状态,(Ⅰ)式变为⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎭⎫ ⎝⎛+=⎪⎪⎭⎫ ⎝⎛-=011220220θθτσσr v v r r r p r r p , 推出⎪⎪⎪⎩⎪⎪⎪⎨⎧=+=-=011220220θθτσσr vvr r r p r r p取v rp σ或 vp θσ为应力集中系数,绘制应力集中系数分布曲线。
① σr 分布:洞室开挖后,围岩中的径向应力σr 始终小于岩体初始应力γH 。
即σr <P 0=γH②切向应力θσ大于v p ; 在洞壁上最大θσ=2P 0=2γH ③ 当r=6r 0时00611P P r ≈⎪⎭⎫ ⎝⎛-=σ 00611P P ≈⎪⎭⎫ ⎝⎛+=θσ离洞中心三倍洞直径的地方;其应力基本上为岩体的天然应力。
所以,洞室开挖的影响范围是三倍洞直径。
2).K 0=0即0=h p ,仅有垂向应力(单向受压状态)从式中看,当r 一定时,r σ,θσ,θτr 是θ的函数, 讨论① 0=θ(水平方向)σr 分布⎥⎦⎤⎢⎣⎡-=44022023r r r rp v rσ 0r r =时,0=r σ 0r r >>时,0→r σ取04r r =,059.0=r σ05r r =,0384.0=r σ; 010r r =,0099.0=r σvr p σ对r 求导,得:02r r =,rH r 83=σ最大。
巷道掘进后,围岩的原岩应力平衡状态被破坏,围岩应力重分布进入塑性状态,一方面导致塑性区扩大和两帮位移增长,另一方面使一些强度低的岩石因应力达到强度极限而破坏,随之产生裂缝或剪切位移,造成破坏岩石的大量脱落,为了保持围岩,稳定以及硐室结构的安全,在围岩稳定性评价的基础上进行支护,以改变围岩的物理力学性质,改善围岩内部的应力及应力分布,降低支承压力区的承载能力,使支承压力向围岩深部转移,从而提高围岩的稳定性。
巷道掘进时,巷道围岩压力以压应力为主,在巷道两帮和拱脚处应力较大,容易出现应力集中,导致出现片帮和巷道拱脚岩石破裂。
顶板由于卸载作用,巷道出现最大切向拉应力,并分布在顶板一定范围内,极易出现顶板岩石离层现象。
施工巷道的围岩应力及来压分析施工巷道的围岩应力及来压是施工过程中需要重点分析和考虑的问题。
在进行施工工程之前,首先需要对巷道围岩进行调查和评估,确定巷道围岩的性质、厚度、倾角等基本参数。
然后,根据巷道的设计方案和施工工艺,结合现场的实际情况,进行应力及来压分析。
巷道的围岩应力分析是指对围岩的应力状态进行研究。
一般来说,巷道围岩的应力状态可以分为两种情况:一是围岩受到一定的外力作用而形成的应力状态;二是围岩失去平衡后自身产生的应力状态。
在施工巷道中,由于施工过程中的爆破、挖掘和支护等工作会对围岩产生应力作用,因此需要对这些外力作用进行分析和评估。
巷道围岩的来压分析是指对围岩受到力的压力大小进行分析。
在施工巷道中,围岩受到上部土层、建筑物及水压力等的压力作用,这些压力会对巷道的稳定性和安全性产生影响。
因此,需要对来压力进行分析和计算,确定巷道围岩的受力情况,为巷道设计和施工提供依据。
对于施工巷道的围岩应力及来压分析,可采用数值模拟方法进行研究。
通过建立巷道围岩的数学模型,输入相关参数和边界条件,使用计算机软件对围岩的应力和来压进行模拟计算。
根据计算结果,可以评估巷道围岩的稳定性,并提出合理的支护措施和安全预警机制。
在进行施工巷道围岩应力及来压分析时,应注意以下几个方面:一是要充分考虑巷道围岩的地质条件,包括岩性、结构、断裂、节理等因素;二是要准确测量巷道围岩的应力和来压,尽可能获取真实的数据;三是要合理选择适应的数值模拟软件和方法,确保计算结果的可靠性;四是要结合施工工艺和支护措施,进行综合分析和评估。
总之,施工巷道的围岩应力及来压分析是保证巷道工程顺利施工和安全运营的重要工作。
通过合理的分析和评估,可以为巷道设计和施工提供科学依据,保障巷道的稳定性和安全性。
四、 围岩应力分析地下洞室的开挖,会产生应力重分布,围岩应力不仅与天然应力场有关,而且还与洞室的开挖有关。
基本假定:①岩体均质,各向同性,连续体。
②无构造应力作用,仅由自重应力形成天然应力场,其大小为z z ⋅=γσ取K 0=0,1/3,1三种情况下的应力场③忽略洞室高度上的应力场变化,即认为洞顶和洞底处的天然应力相同,则有H P v γ=(H 为洞中心的深度)④为平面应力问题K = 10P = P h vPv = r . Z1. 圆形洞室洞室开挖前的天然应力为H P v γ=v h P K P 0=围岩中的径向应力r σ,切向应力θσ以及剪应力θτr 可按下述公式计算:(Ⅰ)r0——洞室半径(m)r——自洞室中心算起的径向距离(m)θ——自水平轴算起的极坐标中的角度p v——垂直方向的压应力(MPa)(=γH)p h——水平方向的压应力(MPa)(= K0 p v)讨论:1).K 0=1; p h =p v (K 0=1静水压力式的天然应力场) 洞室周围处于等压状态,(Ⅰ)式变为⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=⎪⎪⎭⎫ ⎝⎛+=⎪⎪⎭⎫ ⎝⎛-=011220220θθτσσr v v r r r p r r p , 推出⎪⎪⎪⎩⎪⎪⎪⎨⎧=+=-=011220220θθτσσr vvr r r p r r p取v rp σ或 vp θσ为应力集中系数,绘制应力集中系数分布曲线。
① σr 分布:洞室开挖后,围岩中的径向应力σr 始终小于岩体初始应力γH 。
即σr <P 0=γH②切向应力θσ大于v p ; 在洞壁上最大θσ=2P 0=2γH ③ 当r=6r 0时00611P P r ≈⎪⎭⎫ ⎝⎛-=σ 00611P P ≈⎪⎭⎫ ⎝⎛+=θσ离洞中心三倍洞直径的地方;其应力基本上为岩体的天然应力。
所以,洞室开挖的影响范围是三倍洞直径。
2).K 0=0即0=h p ,仅有垂向应力(单向受压状态)从式中看,当r 一定时,r σ,θσ,θτr 是θ的函数, 讨论① 0=θ(水平方向)σr 分布⎥⎦⎤⎢⎣⎡-=44022023r r r rp v rσ 0r r =时,0=r σ 0r r >>时,0→r σ取04r r =,059.0=r σ05r r =,0384.0=r σ; 010r r =,0099.0=r σvr p σ对r 求导,得:02r r =,rH r 83=σ最大。
θσ 分布4402202321rr r r p v ++=θσ 0r r =时,v p 3=θσ, 为最大切向应力0r r >>时,v p →θσ, 为岩体初始应力0=θτr ②θ=90°时r σ44022023251rr r r p v r+-=σ0r r =时,0=r σ 0r r >>时,v r p rH ==σ 06r r =时,v r p 9317.0=σ 010r r =时,v r p 97515.0=σ对v r p σ求导,得02.1r r =时,最小, 为拉应力,由于较小,一般未予考虑。
θσ 分布4402202321rr r r p v -=θσ 0r r =时,rH p v -=-=θσ出现拉应力。
03r r =时,0=θσ 03r r >时,0>θσ但很小如03r r =时,rH 271=θσ 05r r >时,0→θσ0=θτr③θ=45°r σ分布⎪⎪⎭⎫ ⎝⎛-=220121r r p v rσ 0r r =时,0=r σ 06r r >时,v r p 21≈σ θσ分布⎪⎪⎭⎫⎝⎛+=220121r r p v θσ 0r r =时,v p =θσ06r r >时,v p 21→θσ0=θτr⎪⎪⎭⎫ ⎝⎛-+=44022032121r r r r p v r θτ 0r r =时,0=θτr : 02r r =时,v vr p p 65625.0169==θτ 03r r =时,v v r p p 5926.05432==θτ: 04r r =时,v r p 5566.0=θτ 05r r =时,v r p 5376.0=θτ: 06r r =时,v r p 5266.0=θτ 010r r =时,v r p 5098.0=θτ: 0r r >>时,vr p 21=θτ 对r 求导,得时,v r p 6667.0=θτ为最大13)K0=3(Ⅰ-3)讨论① 0θ=σ分布r⎪⎪⎭⎫ ⎝⎛-+=402032131r r r r p v rσ 0r r =时, 0=r σ 02r r =时, v r p 4375.0=σ03r r =时, v r p 395.0=σ04r r =时, v r p 371.0=σ 010r r =时, v r p 3399.0=σ0r r >> 时,h p =θσθσ 分布440220321rr r r p v ++=θσ 0r r =时,02r r =时, v p 229.1=θσ03r r =时, v p 086.1=θσ 04r r =时, v p 046.1=θσ 0r r >>时,v p =θσ0=θτr② 90=θ r σ44022021rr r r p v r+-=σ 0r r =时,0=r σ 02r r =时,v r p 5625.0=σ03r r =时,v r p 79.0=σ 04r r =时,v r p 879.0=σ05r r =时,v r p 9216.0=σ010r r =时,v r p 98.0=σ; 0r r >>时,v r p →σ 最大θσ分布⎪⎪⎭⎫ ⎝⎛-+=44022032131r r r r p v θσ 0r r =时,0=θσ;02r r =时,v p 4375.0=θσ; 03r r =时,v p 395.0=θσ 04r r =时,v p 371.0=θσ; 05r r =时,v p 3584.0=θσ010r r =时,v p 3399.0=θσ; 0r r >>时,v v p p 31333.0==θσ对r 求导,得03r r =时,v v r p p 444.094==σ 最大。
0=θτr ③θ=45°r σ 分布⎪⎪⎭⎫ ⎝⎛-=220132r r p v rσ 0r r = 时,0=r σ; 02r r =时,v r p 21=σ 04r r = 时,v r p 85=σ;010r r =时,v r p 66.0=σ02r r = 时,v r p 21=σ; 02r r =时,v r p 21=σ0r r >>时,v v r p p 667.032==σθσ分布⎪⎪⎭⎫ ⎝⎛+=220132r r p v θσ 0r r =时,v p 34=θσ; 02r r =时,v p 65=θσ; 03r r =时,v p 74.0=θσ 04r r =时,v p 708.0=θσ; 010r r =时,v p 673.0=θσ 0r r >>时,v v p p 667.032==θσθτr 分布⎪⎪⎭⎫ ⎝⎛-+=44022032131r r r r p v r θτ 0r r =时,0=θτr02r r =时,v r p 4375.0=θτ; 03r r =时,v r p 395.0=θτ 04r r =时,v r p 371.0=θτ; 05r r =时,v r p 3584.0=θτ 010r r =时,v r p 3399.0=θτ0r r >>时,v r p 31=θτ03r r =时,v v r p p 4445.094==θτ为最大综合上述讨论,有以下结论: ①无论K 0何值,洞壁上的径向应力σr =0 ;水平方向上,趋于水平应力的大小,如K 0=1时,σr =P h ;K 0=0时,σr =0; K 0=1/3时,σr =P h / 3;垂直方向上,趋于垂直应力的大小,即σr =P v ;不论 K 0=0,=1/3,或1。
② 不论K 0为任何值,洞壁上的剪应力0=θτr ,且在水平方向和垂直方向上的剪应力0=θτr , 其它方向的剪应力不为零。
③关于洞壁上的切向应力θσ:a.在水平方向上(θ=0°),洞壁上的θσ随K 0值的增大而减小,b.在垂直方向上(θ=90°,洞顶)洞壁上的θσ随K 0值的增大而增大,且当310<k 时,出现拉应力。
c.当 K 0<1时,σθ随θ的增大(0→90而减小。
d.当K 0>1时,σθ随θ的增大而增大。
④围岩在洞壁上的各应力中,只有切向应力不为零,且最大切向应力发生在洞壁上。
⑤在围岩应力中,r σ,θτr 都小于v p ,而θσ可以大于v p稳定验算部位:侧墙和顶拱,验算的应力是切向应力θσ,压或拉。
如K 0=1时,验算θσ压。
K 0=0时,验算侧墙处压θσ>r σ和顶拱θσ>r σ。
出现拉应力的区域为: 120~90~60=θ2、椭圆形洞室椭圆形洞室两个半轴设为q 1、q 2,则参数方程为βcos 1a q x = βsin 2a q y =由于洞室围岩中的最大切向应力发生在洞壁上,则在稳定性分析需计算洞壁上的切向应力。
按弹性力学原理,开挖后所引起洞室边界上的切向应力为: (Ⅱ)21q q m =β——椭圆偏心角由于P v 、P h 为主应力,则0=xy τ 将m 代入(II )并整理,得222222121222221sin )(2]sin ))[((q q q p q q q q q p p vh v t +-+-+-=ββσ(Ⅱ-1)与圆形洞室讨论一样,可取不同的q 1/q 2宽高比,讨论K 0取不同值时的变化(取β=0、45°、90°)讨论(1).取宽高比25.0//21==q q h B A 点 (β=900):K 0=0时,v t p -=σ(最大拉应力) K 0=1/3时,v t p 2=σ,压应力 K 0=1时,v t p 8=σ,压应力最大 B 点 (β=00):v t p )2~0(=σ,K 0=0、1/3、1(2).宽高比B/h=0.5 A 点(β=900):00=k ,v t p -=σ拉应力 310=k ,v t p 21=σ10=k ,v t p 4=σB 点(β=00):v t p )2~1(=σK 0 = 0、1/3、1(3).宽高比B/h=2A 点(β=900):K 0=0,σt = - p v , K 0=1/3,v t p 31-=σ 拉应力 K 0=1,σt = p v B 点(β=00):v t p )5~8.3(=σ 均为压应力。