无横向水源补给非稳定渗流作用下边坡稳定性分析.docVIP

无横向水源补给非稳定渗流作用下边坡稳定性分析 某某 摘要：无横向水源补给的非稳定渗流是库区普遍存在的一种渗流工况，针对这种具体工况，结合边坡分析软件GEO-SLOPE和PLAXIS的计算特点，分析库水位下降过程中，各种土质边坡的渗流场、变形场与边坡的安全系数随水位变动的规律。得出变形场受渗流场影响的规律，以及土体塑性变形对边坡的安全储备有贡献，初步探讨了渗流、变形与稳定的关系。 关键词：无横向水源补给 非稳定渗流 边坡稳定性 三峡工程建成后库区每年将有30m的水位涨落，且三峡库区属多雨地区，库水位大幅度、周期性涨落，加之强降雨，将使滑坡区渗流场产生剧烈变化，从而威胁到库区两岸大量边坡的稳定[1]。目前有很多关于渗流对边坡稳定性影响的研究和不少边坡分析软件，但针对某种具体工况的详细分析还比较少，没有形成一种规范的方法，而实际工程中广泛采用的还是一种简化的粗略的计算方法。本文针对库区一种比较普遍的渗流情况即：无横向水源补给的非稳定渗流，结合大型边坡分析软件GEO-SLOPE和PLAXIS的特点进行边坡稳定性分析，探讨渗流场的规律以及对边坡稳定性的影响，为设计提供依据。 1．研究现状及发展趋势 无横向水源补给的非稳定渗流是库区一种比较普遍的渗流情况，其特点是边坡后沿的地下水水头不是定值，它是随坡外水位变化而变化的[2]。渗流场计算主要是确定浸润线的位置，目前有图解法、模拟实验法和数值解法等[3]，数值计算比较准确，但边界条件难以确定，目前趋向于使用现成的计算软件如：GEO-SLOPE和FLAC3D等。这方面不少学者进行了探讨如：时卫民,郑颖人根据Boussinesq非稳定渗流微分方程，得到了库水位下降时坡体内浸润线的简化计算公式[4]；陈昌禄等提出了确定稳定、非稳定渗流自由面的变单元渗透系数法，可应用于稳定和非稳定渗流自由面的计算[5]。但是针对某种具体工况的规律的分析还比较少。 在考虑渗流对边坡稳定性的影响时，主要采用“代替容重法”、“静水压力法”等基于传统极限平衡理论的方法[6]。也有不少研究成果如：时卫民，郑颖人通过算例分析得出库水位的下降过程中，坡体存在一个最危险的水位，在这个水位坡体的稳定系数最小，该位置一般在滑体的下1/3处[7]；刘新喜将暂态孔隙水压用于滑坡的极限平衡分析，结果表明渗流对滑坡稳定的影响受控于滑坡土的入渗能力和滑坡结构形态[8]。但变形与稳定是密切相关不可分割的，传统的极限平衡法的处理是不合理的，应该是将变形计算与稳定分析融合在一起，这也是边坡稳定性分析的一个趋向。目前有一些这方面的研究但还不够深入不够具体如：黄茂松、贾仓琴比较了非饱和非稳定渗流的极限平衡法和位移有限元法，但只用位移法分析了稳定渗流的一种简单情况，非稳定渗流情况的位移法还需要进一步研究[9]。 2. GEO-SLOPE与PLAXIS的理论基础和分析特点 GEO-SLOPE是一个多模块的大型边坡分析软件，其中SEEP是针对渗流计算的模块，是以非饱和土渗流理论为基础的有限元计算程序，渗流基本微分方程为： （1） 式中为渗透系数；为水的密度；重力加速度；为比水容量，定义为体积含水量对基质吸力的偏导数的负值；是基质吸力的函数。SEEP适于计算水位涨落、降雨等工况下的边坡渗流场。考虑渗流对边坡稳定影响的分析过程是，在SEEP中计算渗流场，然后将计算得到的浸润线位置传到另一个模块SLOPE中采用条分法分析边坡的稳定性。 PLAXIS是一个大型有限元分析软件，浸润线的位置可手动输入也可调用它的渗流计算模块，但只能计算稳定渗流的情况。它是基于位移变形理论的有限元程序，在渗流场计算的基础上采用PLAXIS，可以分析水位涨落工程中边坡的位移、应力、应变、塑性区等，还可以直接采用有限元强度折减法计算边坡的安全系数。强度折减即将土的抗剪强度参数和按下面的式子等比例减少，即 ， （2） 利用和通过有限元计算边坡的应力应变，考察塑性区的发展，直至破坏，则破坏时的折减系数即为安全系数，此时的安全系数具有强度储备的物理意义[9]。 PLAXIS在计算渗流时是用稳定渗流近似代替非稳定渗流，但对于无横向水源补给的情况，由于边坡后沿水头不是定值，PLAXIS是无法计算的。即使用其他方法得出了边坡后沿水头，它在某一时刻的分析结果与SEEP也有较大差别。图1图2分别是PLAXIS与SEEP对水位从175m高程下降到145m高程时的分析结果。由SEEP的结果可见，溢出点高于坡外水位线，有一段自由渗出高度是符合实际情况的[10]；而PLAXIS的计算结果可见溢出点于