基于刚体极限平衡的土石坝稳定性研究.docxVIP

基于刚体极限平衡的土石坝稳定性研究 摘要：服役条件下坝坡稳定性验算是土石坝设计与运行中的一项重要内容,其稳定性直接关系到土石坝工程的安全与稳定。服役条件下库水位升降、坝体材料劣化及渗透能力的改变均对坝坡稳定性产生重要影响,为研究各因素对土石坝稳定性的影响,基于刚体极限平衡原理,结合Fredlund非饱和土抗剪强度理论,采用简化Bishop法分析了库水位升降、材料劣化及坝体渗透能力变化等条件下某长期服役土石坝坝区土体的基质吸力及坝坡稳定性。结果表明,坝区土体基质吸力随土体渗透系数及体积含水率的增大而减小;由于材料劣化导致的抗剪参数降低会促使坝坡的稳定性降低;在水位较低时,由材料劣化导致的土石坝稳定性降低占主导作用;在水位较高时,由水位升高导致的土石坝稳定性降低占主导作用。该方法适用于长期服役条件下土石坝稳定性的研究,可为类似工程提供参考。 1 引言土石坝是一种重要的挡水建筑物，其稳定性直接关系着运行期的安全。土石坝坝坡稳定性受多种因素共同影响，其中材料劣化导致的抗剪强度变化、水位升降引起的土体孔隙水压力与渗透系数改变等对土石坝的稳定性影响较大2 刚体极限平衡法原理2.1 土石坝渗流计算库水位升降会引起非稳定的渗流场，在非稳定渗流场作用下，在进行坝坡稳定性研究时，传统的Mohr Coulomb屈服准则将不再适用，对此FREDLUND D G等2.2 坡外水的处理坡外水的处理采用“等效置换法”，坡外水压力处理示意图见图1。计算公式为：式中，u2.3 简化Bishop法基于刚体极限平衡理论，采用简化Bishop法计算坝坡稳定性时，假定滑动面上的土体条分，并假设滑裂面，根据条块受力特征推导出坝坡安全系数的表达式。其中由竖直方向受力平衡可得：式中，W由滑裂面极限平衡可得：式中，c由式（3）、（4）可得：式中，b采用Bishop法计算时，滑裂面OAB的滑动力矩为：抗滑力矩为：式中，N为边坡条分数量；R为滑裂面半径；E由力矩平衡可得坝坡稳定系数的表达式为：3 算例分析3.1 水位升降过程中土石坝稳定性某土石坝为均质折线坝，轴线呈“L”形，折角约106°，坝高31.40 m，坝顶长343 m，坝顶宽5m，坝顶高程1 756.50m。折角左边大坝上游坡比为1∶2.6，下游坡比从高至低依次为1∶2.3、1∶2.5、1∶3.0。该土石坝死水位、正常蓄水位、设计洪水位、校核洪水位分别为1 736.00、1 751.60、1753.50、1 754.00m。选取折角左边大坝为研究对象，其剖面见图2。该土石坝1954年修筑完成，现已服役66年。勘察后发现坝体在长期的服役过程中，土石坝填料的材料抗剪强度参数在水的长期作用下逐渐劣化，该土石坝在设计初期和服役66年后的土体材料物理力学参数统计见表1。该土石坝由坝体区、排水区、坝基区3部分组成。基于Fredlund理论，根据表1土石坝材料物理力学参数的设计值及当前值，采用商用软件Geostudio中的SEEP/W模块计算得到坝体区土体基质吸力与渗透系数、体积含水率的关系曲线见文献坝体区设计值曲线始终位于当前值曲线下侧，且在坝体区土体基质吸力逐渐增大的过程中，土体的渗透系数和体积含水率均逐渐减小。这是基质吸力作为土体基质对水吸持能力的体现，其通过吸附作用和毛细作用产生影响，基质吸力增大，使得土体中部分水被附着于土体颗粒上，土粒孔隙间距降低，渗透系数减小，此外土体中部分水被附着导致孔隙水含量降低，土体体积含水率随之降低。将水位上升过程中的初始水位设定为死水位1 736.00 m，上升的最大高度为校核洪水位1 754.00 m，速度为1 m/d。水位下降过程中的初始水位定为校核洪水位1 754.00m，下降的最小高度为死水位1 736.00m，速度为1m/d。图3给出了由Geostudio中SEEP/W模块计算得到的库水位升降过程中坝体内浸润线的变化过程。为研究土石坝稳定性与水位升降之间的关系，共设计了4种工况：(1)工况1。水位上升设计值；(2)工况2。水位上升当前值；(3)工况3。水位下降设计值；(4)工况4。水位下降当前值。对4种工况下的土石坝进行稳定性研究，得到各工况下土石坝稳定性与库水位升降之间的关系，见图4。由图4可知，水位上升时，土石坝的稳定性逐渐降低；水位下降时，土石坝的稳定性逐渐增加；此外，无论是水位上升还是下降，设计值的曲线始终位于当前值曲线的上方，且水位较低时，两曲线的间距较大，说明水位较低时，由抗剪强度参数所导致的土石坝稳定性降低占主导，而水位较高时，由水位所致的土石坝稳定性降低占主导。3.2 抗剪参数对土石坝稳定性影响在长期服役过程中，土石坝的抗剪强度参数会逐渐劣化。图5给出了水位