有限单元法分析九畹溪水库大坝渗流稳定.pdfVIP

第16卷第2期 水利科技与经济 V01．16 No．2 2010年 2月 WaterConservancyScienceandTechnologyandEconomy Feb．，2010 有 限 单 元 法 分 析 九 畹 溪 水 库 大 坝 渗 流 稳 定 谭远春 (秭归县水利水电勘察设计院，湖北 秭归 443600) [摘 要]对九畹溪水库大坝坝体稳定、渗流稳定进行了分析，其中渗流场的分析采用有限单 元法计算，运用加权残余法对模型进行离散，渗透坡降采用达西定律计算得到，采用强度储备 法研究了大坝的失稳破坏机理 ，并根据提 出的应力校核法和强度储备法的安全评价准则对大 坝的安全储备进行了评价，成果对实际工程应用具有一定的价值，也为其它类似工程提供 了参 考。 、 [关键词]有限单元法；渗流稳定；坝坡稳定 [中图分类号】 TV135 [文献标识码】 B [文章编号】 1006—7175(2010)02—0152—03 1 工程概况 3 大坝渗流分析 九畹溪流域位于长江三峡南岸，流域位置在 El1O。 3．1 渗流计算理论及计算程序 15～EllO。45，N30。43一N31。00之间，属长江 的一级支 渗流分析领域中的计算方法有很多，其中应用较为 流。九畹溪流域发源于秭归县云台荒南部千家屋场 ，海 广泛的3种方法：解析法、数值法和 电拟法。 拔高程2041m，至下游秭归县九畹溪镇聚集坊汇入长江。 解析法…可以直接求解渗流的问题 ，计算比较精确， 河流全长44．6km，控制 流域面积590km ，总落差 可以计算渗流场中任意一点处的渗流要素，但是这种方 l840m，平均总坡降20．78％，九畹溪流域流向大致由北 法计算比较复杂，目前只能对几种比较简单的渗流边界 向南。九畹溪水库位于九畹溪流域中下游的花桥河与主 有解答。 河道交汇处，坝址集水面积410km ，占总流域面积的 电模拟 是利用渗流与电流在数学表达式方面的相 69．49％。九畹溪流域已建电站 l4处，总装机24000kW， 似性，建立用电流现象模拟渗流现象的试验，但是这种方 其中l1座电站均在九畹溪水库坝址 以上，3座 电站建在 法模拟原型的复杂边界和 内部结构制作模型费时费力。 九畹溪水库坝址以下。九畹溪水库的建设除能保证其下 本文用的是数值方法 J，数值方法 目前主要有有 限 游九畹溪漂流产品的安全性和稳定性，同时对下游 电站 差分法、边界元法和有限单元法。 的发电效益和整个流域的防洪 、防汛及为三峡库区拦沙 有限差分法是从微分方程出发，将研究区域经过离 起到积极的作用。 散处理后 ，近似地用差分 、差商来代替微分 、微商。这样一 来，基本微分方程和边界条件的求解可归结为解一个线 2 坝址区工程地质条件 性方程组，所得结果为数值解。有限差分法的优点是其 本区位于鄂西南褶皱山地，属构造剥蚀型中高山地 原理易懂，算式简单，有较成熟的理论基础。其缺点是往 貌区。区域地质构造环境稳定，地震活动微弱。据 《中国 往局限于规则的差分网格，对曲线边界和渗透介质的各 地震动参数区划图》(GB18306～2001)，工程所在区域地 向异性模拟 比较困难。 震动峰值加速度为0．05g，地震动反应谱周期特征值为 边界元法 是通过把求解域边界剖分为若干个单 0．35s，相应地震基本烈度为Ⅵ度 。 元