降雨入渗对岩质边坡的稳定影响分析.pdfVIP

降雨入渗对岩质边坡的稳定影响分析 张卓1 练继建1 杨晓慧2 1天津大学建筑工程学院，天津，300072 z秦皇岛市保神交通监理有限公司，河北秦皇岛，066000 摘要 雨水入渗改变了岩质边坡非饱和区渗流场的分布，是引发边坡滑坡的主要因素之一。本文以 非饱和土抗剪强度理论为基础，引入了渗透力的概念，利用有限单元法进行了降雨条件下岩质边坡稳 定性的分析。最后，利用编制的跚rAI狐，n、Y边坡稳定程序进行了算例分析。结果表明，岩体边坡的稳 定安全系数随着雨水的不断渗入在逐渐地降低，降雨入渗对边坡的稳定性有很大的影响。 关键词 稳定 非饱和 有限单元法渗透力 1前言 很多地区的边坡失稳与降雨或雾化雨的人渗有着密切的关系。雨水入渗会使边坡内非饱 和带的吸力锐减，暂态水荷载增加，地下水位升高，这不仅增加了边坡的滑动力，同时也减 小了边坡的有效凝聚力和摩擦角，使抗剪强度降低。随着实际工程中雾化雨引发边坡滑坡事 故的不断发生，人们对降雨条件下边坡稳定性方面的研究在不断地深入，并取得了一定的成 果。以刚性极限平衡为原理的条分法由于便于计算，故在边坡稳定分析中一直得到人们的青 睐，但该方法在结合渗流场计算时存在着一些问题。目前应用的条分法主要采用条块周边的 孔隙水压力来考虑水的渗透作用，这种算法只考虑了条块侧边水压力大小而忽略了作用点所 发生的力矩影响，因此，在渗流场变化很大的情况下，条分法计算的精度就会受到很大的影 响。由于渗透力等价予浮力与作用在条块周边各水压力的总和[1]，故本文采用把岩质边坡划 分成有限单元，将岩体单元各边上的水压力转换成一个渗透力来考虑雨水渗流的影响。该法 能很好地利用有限元法计算的饱和一非饱和渗流场分析边坡的稳定性。 2非饱和土抗剪强度理论 雨水的人渗将改变边坡内非饱和区渗流场中压力水头的分布，尤其是强降雨或暴雨使这 种分布改变更大，有时甚至形成暂态饱和区，使得非饱和区的基质吸力大大降低，抗剪强度 明显下降。在降雨条件下，由于非饱和区基质吸力的变化对边坡稳定性影响很大，所以，本 文采用非饱和土的抗剪强度理论进行边坡稳定分析。 非饱和土的抗剪强度理论公式如下[2]： r，一C’+o一“。)tan97+(“。一‰)tan矿 (1) 式中：c7为有效凝聚力；(盯～‰)为破坏时滑动面上的净法向应力；口为滑动面上的法向总应 力；Ua为滑动面上的孔隙气压力；97为与净法向应力(口--Ua)有关的内摩擦角；(‰一‰)为破 坏时滑动面上的基质吸力；“。为滑动面上的孔隙水压力；在非饱和区，矿为与基质吸力 (‰一“。)有关的内摩擦角，在饱和区，矿=97。 在计算中，设孔隙气压力为大气压，即Ua=O，则非饱和土的抗剪强度公式可写为 rf—f7-Jr-atanPt一‰ta矽 (2) 3有限单元法‘3]I[4][5] 下面以圆弧形滑动面为代表，来说明用于边坡稳定分析中的有限单元法。 作用在单元上的力如图1所示。作用在三角形单元巧忌上的渗透体积力为 图1单元受力图和滑动面 F，一rJA F一(￡}一挫{竞)一一尬{塞h} fa ] L ay J G—yl△ (3) (4) (5) 式(4)和式(5)中：’，为水容重；h为测压管水头；J：和J，分别为单元在X、y方向的渗 透坡降；)，，为岩体容重(浸润线以下为饱和区，取浮容重；浸润线以上为非饱和区，取自然 容重)；A为三角形单元面积。 单元的有效自重和渗透力均作用在单元形心上，将它们分解为沿圆弧滑动的切向力T与 半径方向的法向力N。滑动力T作用于形心，而法向力N只有作用在滑动面才起阻滑作用。 将抗滑力的所有分力移置滑动面，则得 fN=[(yl+Ⅳy)cosal一。sinal]△ IT一[(M 4-玎，)sina+玎，co趿]△ (6) 由图l可知，抗滑力的垂直分力向下延伸到滑动面上，力矩效果不变；水平分力E由单 元形心移置到滑动面上时就被增大了一个力矩M1=F，(Rcosa，--rcosa)。利用非饱和土抗剪 强度理论，建立围绕圆心的力矩平衡方程式，得 ·310· 7∑(n—M1)一R(∑cl+∑Ntan驴7一∑“。Ita矽)一0 (7) 由式(7)得边坡稳定安全系数77，即 R(∑f7z+∑Ntanp7一∑甜。ltang，b)旷———1而i面厂一 (8) 式中：f7、∥分别为滑动面上的有效粘聚力和内摩擦角；R为滑弧半径；r为单元形心与滑弧 圆心的距离；a为单元重心偏离铅垂线的角度；a7为单元重心与下方滑弧交点偏离铅垂线的角 度。 根据上述计算公式，编制了雨水入渗条件下的岩坡稳定计算程序STABILITY。值得指 出的是，‰在非饱和区为基质吸力，用负值表示；在饱和区为压力水头，用正值表示。在有 限单元法中，边坡被划分成许多单元，一个单元的孔隙水压力用单元各节点的压力水头平均 值来表示。渗透坡降J的