基于非线性有限单元法的堆石坝施工期坝坡稳定性研究.docxVIP

基于非线性有限单元法的堆石坝施工期坝坡稳定性研究 1 期安全渡未来 板岩水库是安全经济水库。它对地形、地质和气候有很强的适应性，在近几十年来，它一直在迅速发展，总结了整个项目和施工经验。混凝土面板堆石坝施工期临时渡汛坝体的坝坡稳定问题是面板堆石坝施工期安全渡汛中的重要问题之一,它不仅事关大坝施工期安全,而且对大坝后期的运行安全也有重大影响。随着面板堆石坝高度的增加,坝体填筑方量大,填筑周期长,施工期中往往由堆石体在没有混凝土面板的保护下来临时抵挡施工期洪水,此时,在施工期渡汛过程中坝坡的稳定性对于工期与施工安全以及后期坝体的工作性态都有重要影响。本文以重庆鲤鱼塘面板堆石坝为实例,在应用非线性有限元方法研究渡汛期坝体稳定渗流场的基础上,应用瑞典条分法分析面板堆石坝施工期临时渡汛断面坝坡的稳定性,给出了坝坡稳定性评价的结果,为合理地评价面板堆石坝施工期临时渡汛断面安全渡汛提供了参考。 2 上、下游坝坡稳定性 面板堆石坝施工期临时渡汛断面坝坡稳定性分析,主要是研究临时渡汛断面在汛期洪水及坝体自重作用下,上、下游坝坡的稳定性。鉴于堆石材料的压缩性远小于一般的土体,其单位贮存量较小,故在汛期洪水作用下坝体内形成的渗流场可近似视为稳定渗流场,渗流作用以渗透体积力的形式作用于坝体结构,计算中的这种处理是基于简化计算的目的提出的一种近似计算方法,会得出偏于安全的结论。 2.1 初始流速法求解渗流场的基于自由面位置 设水头函数为?=z+pr?=z+pr,水头?在求解区域内满足的基本方程为: ??x[kxx???x+kxy???y+kxz???z]+??y[kyx???x+kyy???y+kyz???z]+??z[kzx???x+kzy???y+kzz???z]-Q=0(1)??x[kxx???x+kxy???y+kxz???z]+??y[kyx???x+kyy???y+kyz???z]+??z[kzx???x+kzy???y+kzz???z]?Q=0(1) 同时,水头函数?还必须满足边界条件: ①已知水头边界:?=?0(2) ②已知流量边界:lxvx+lyvy+lzvz=vn(3) 稳定渗流场的计算可通过求解如下方程得到: [Η]{?}=[F](4)[H]{?}=[F](4) 解出各结点?值后,可进一步求得各点流速和渗透坡降;在渗流场的自由面上,压力p=0,且无渗流量通过,因此,其必须同时满足以下两个条件: ?=z及vn=0 由于自由面位置和下游坝面渗流逸出点的位置事先未知,故需通过迭代的方法来求解。本文应用初始流速法来进行迭代求解自由面位置,所得到的 式中R1为坝体浸润线以上区域。由于本方法令自由面以上区域内的流速为零,在计算{Fn+1}时只需判断各高斯点是否满足条件?nz即可,因此处理十分简单,且收敛速度较快。 2.2 法线与竖直线的交互作用 本文采用瑞典条分法、来分析坝坡稳定性,如图1,AC是假定的滑动面,其滑动圆心为O,半径为R;将滑动体ABC分成若干土条,其中任意第ⅰ条的受力分析为:①土条自重Wi,方向竖直向下,其值为Wi=γbihi,将Wi分解为通过滑弧圆心的法向力Ni和滑弧切向的剪切力Ti;若以θi表示该土条底部中心点的外法线与竖直线的交角,则有Ni=Wicosθi及Ti=Wisinθi;②作用在土条底面上的法线反力ˉΝi,与Ni大小相等,方向相反;③作用在土条底面上的剪力ˉΤi,其可能发挥的最大值等于土条底面上土的抗剪强度与滑弧长度li的乘积,方向与滑动方向相反。当坡体处于稳定状态并假定各土条底部滑动面上的滑动稳定安全系数均等于整个滑动面上的滑动稳定安全系数时,其实际发挥的抗剪力为: ˉΤi=τfiliFs=cli+ˉΝitgφFs 将整个滑动体内各土条对圆心取力矩平衡,∑ΤiR=∑ˉΤiRe 当用有效应力法进行分析,瑞典条分法的稳定安全系数为: Fs=∑(C′li+(Wicosθi-μili)tgφ′)∑Wisinθi 当坡体内存在稳定渗流场时,瑞典条分法的稳定安全系数为: Fs=∑(c′li+[Wicosθi-(μ2-γwzi)li]tgφ′)∑Wisinθi 式中,Wi—土条重量,浸润线以上用湿容重,浸润线与坡外水位之间用饱和容量,坡外水位以下用浮容重计算,Zi—坡外水位高出土条底面中点的距离。 3 坝体坝体概况 鲤鱼塘水库工程为105m高的混凝土面板堆石坝,位于重庆市开县境内,地处长江北岸二级支流桃溪河上游,坝址距开县县城47km。坝体上游坝坡1∶1.4,下游坝坡在高程400m以上为1∶1.4,400m以下为1∶1.5,最大底宽为318m;坝顶宽10m,坝顶长237m;坝体填筑料为灰岩、砂岩、粉砂岩及泥质岩和开挖料。大坝填筑分两期进行,施工期需临时断面渡汛一次,其填筑分区如图2所示。 3.1 降低+x轴正向回复突变型