基于岩体结构面网络的边坡稳定性分析.docxVIP

基于岩体结构面网络的边坡稳定性分析 在漫长的地质历史进化过程中，岩体形成并形成。其起源和结构复杂，岩石多样，包括结构面和结构体这两个基本结构单元。就一般岩石而言,边坡稳定性的最主要的影响因素并非岩石的性质,结构面对于岩体和边坡的失稳破坏起着主要的控制作用,并在岩体工程的设计与施工中占有重要的地位。 对于无明显滑面岩质边坡,基于岩体结构面网络,运用数值分析方法,定量研究无明显滑面岩质边坡稳定性,能为此类岩质边坡的加固设计与施工提供科学依据和指导。 1 初始基准结构面的搜索和确定 基于岩体结构面网络和边界条件,在靠近岩体结构边坡表面附近选择满足一定条件的一条结构面作为初始基准结构面,由此基准结构面端点确定一个搜索区域,一般为1/4的直角扇形区域,在搜索区域内搜索所有包含其内的岩体结构面,从中选择满足特定条件的一结构面作为新的基准结构面,继续搜索循环,直至边界,最后按照搜索循环步骤首尾连接所有这些搜索到的基准结构面,对于该初始基准结构面可搜索到一系列的滑移路径,形成一树枝形结构。按上述搜索基本思想,对靠近边坡每条满足条件的初始基准结构面进行循环搜索,从而可确定无明显滑面岩质边坡所有可能的滑移路径。 在计算区域内搜索滑移路径时,必须先确定初始的基准结构面,将以此作为初始步,初始基准结构面应满足3个假设条件:① 结构面必须位于计算区域内;② 结构面端点距边界的距离不大于平均迹长;③ 结构面的视倾角0≤θ≤90°。根据这些条件对计算区域内的所有结构面进行筛选,可确定所有满足上述条件的初始基准结构面,其可能有多条,如图1所示。 假设AB为搜索到的某一条基准结构面,在结构面端点B点,以平均迹长的一半为半径,即r=0.5Lˉˉˉr=0.5Lˉ,其中,r为搜索半径;LˉˉˉLˉ为平均迹长长度。作1/4的直角扇形区域,如图1中BCD区域所示,则此扇形即为所要确定的搜索区域。 在搜索区域内确定所有满足下述条件的结构面:① 结构面至少有一端点位于搜索区域内;② 结构面任一端点和B点的连线与水平方向的夹角满足0≤θ≤90°的条件;③ 结构面的视倾角满足0≤θ≤90°;④ 与基准结构面相交以及重叠不相交的特殊结构面,也得考虑将它作为一待选基准结构面。基于岩体结构的复杂性,在搜索区域内可能会出现多条待选基准结构面。若搜索区域内没有待选基准结构面,则判断是否已接近边界,判断依据为结构面端点到计算区域边界的距离小于搜索半径,若判别条件成立,即接近边界,则结束搜索;否则,加倍搜索半径,继续搜索,直至边坡边界。根据上述搜索原则,从计算区域内的每一条初始基准结构面出发,都可搜索出相应的一系列的滑移路径,形成树枝形结构形式,循环搜索从而可以确定计算区域内所有可能的滑移路径。 2 滑动岩体滑动岩体滑动稳定性分析 针对任意可能滑移路径,我们均可利用ANSYS计算出其对应的安全系数F的大小,令∑R为下部岩体为阻止其上部岩体的相对滑动在滑移路径上所产生的最大抗滑力,∑T表示上部滑动岩体相对于下部稳定基岩在滑移路径上所产生的实际剪力值,即下滑力,则边坡安全系数F为∑R与∑T之比。由岩体结构面网络可知,任意条滑移路径为一岩桥和结构面首尾相接的多线段组成的不规则的折线图形,因此抗滑力∑R和下滑力∑T分别由岩桥和结构面所产生,故边坡安全系数为 F=∑R∑T=∑Ri+∑Rj∑Ti+∑TjF=∑R∑Τ=∑Ri+∑Rj∑Τi+∑Τj 式中,∑Ri、∑Ti为滑移路径上所有岩桥所产生的抗滑力和下滑力;∑Rj、∑Tj为滑移路径上所有结构面所产生的抗滑力和下滑力。 2.1 岩桥黏聚力计算的计算 设α为剪切方向角,β为岩桥倾角,γ为产生拉裂纹的主应力方向与水平方向的夹角,θ为拉裂纹与剪切方向的夹角,τxy、σx、σy均源自于有限元计算得到的应力场数据,根据Lajtai和Einstein研究理论,当βγ+0.5π,岩桥将发生受剪破坏;当β≥γ+0.5π,则岩桥将发生受拉破坏。以下为岩桥在不同破坏形式下Ri、Ti的计算公式。 (1) 受剪破坏(βγ+π2)(βγ+π2)。 d=Lcos(β?α)σni=0.5(σx+σy)+0.5(σx?σy)cos(π?2α)?τxysin(π?2α)τsi=0.5(σx?σy)sin(π?2α)+τxycos(π?2α)τri=0.5c2i?2ciσni?????????√,Ri=∫d0τrids,Ti=∫d0τsidsd=Lcos(β-α)σni=0.5(σx+σy)+0.5(σx-σy)cos(π-2α)-τxysin(π-2α)τsi=0.5(σx-σy)sin(π-2α)+τxycos(π-2α)τri=0.5ci2-2ciσni,Ri=∫0dτrids,Τi=∫0dτsids 式中,ci为岩桥的黏聚力;σni为作用在岩桥上的法向应力;τri为岩桥的抗剪应力