面板堆石坝面板出现裂缝工况下渗流分析.pdfVIP

水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期 文章编号: 0559-9350(2005)04-0420-06 面板堆石坝面板出现裂缝工况下的渗流分析 1 1 2 张嘎,张建民,洪镝 (1.清华大学 水利水电工程系，北京 100084；2.黄河上游水电开发有限责任公司，青海 化隆 810902) 摘要：以公伯峡面板堆石坝为例，分析了面板出现裂缝的非常工况下的渗流场，以探讨裂缝对大坝渗流场和渗透 稳定性的影响。根据大坝设计剖面和地形剖面剖分网格，采用固定网格有限元全场渗流分析方法，计算了堆石坝 防渗结构发生水平或竖向裂缝非常工况下大坝渗流场和渗流量。计算结果表明，基岩上部弱透水层的渗透性对大 坝渗流场具有重要影响；面板附近垫层过渡料区中出现较大的渗透坡降，但随着离裂缝距离的增加较快衰减；裂 缝的方位和宽度对渗流量和其附近的垫层过渡料区的渗流场有较大影响。 关键词：公伯峡；渗流；有限元；面板堆石坝；裂缝；宽度；方位 中图分类号：TV41 文献标识码：A 近年来土石坝防渗结构的渗透特性及渗透稳定性研究受到广泛重视，在材料试验、模型试验、原型观 测以及数值模拟等方面均取得了较大进展[1～9] ，目前已发展了多种在防渗结构出现裂缝等非常工况下土石 坝渗流场的分析方法。采用有限元方法并将裂缝作为一个或若干个结点单独赋以边界条件是一种较常用的 方法，不过这种方法不能很好地反映裂缝的尺寸效应，可能会导致较大的计算误差，为了解决这一问题， 张有天[3] [4] 等发展了一种有限元边界元耦合分析方法，杜延龄等 提出了一种特殊的单元格式，但这些方法 编程较为复杂，通用性也受到限制。另外一种较常用的方法是将裂缝与周围的介质概化为各向异性介质进 行渗流分析，不过该方法只适用于裂缝较多情况，且难以精确地得到裂缝附近的渗流场[5] [6] 。许国安 等在 分析小浪底水库坝基防渗墙裂缝对渗流控制的影响时采用了两次计算的方法，即先不考虑裂缝得到整个土 石坝的渗流场，然后再对防渗体裂缝附近区域进行局部网格加密，通过第2次有限元计算得到裂缝附近的 渗流场，这种方法虽然计算量较小，稳定性好，但精度较低，且第2次计算时的边界条件处理较为复杂。 本文基于统一的饱和-非饱和土渗流控制方程，采用固定网格有限元分析方法，选用向量化的共轭梯 度法求解方程组，以提高计算效率。在此基础上，不再采用特殊的单元格式，而是采用全场有限元一次计 算的方法对当面板出现水平和竖向裂缝非常工况条件下的公伯峡面板堆石坝的渗流场进行三维分析，探讨 不同裂缝方位和宽度对大坝渗流场的影响，并初步评价大坝的渗透稳定性。 1 工程概况 公伯峡水电站位于青海省循化撒拉族自治县和化隆回族自治县交界处，是黄河上游龙羊峡至青铜峡河 段中的第4个大型梯级电站，总库容6.2亿m3 ，总装机容量1500MW，是一个以发电为主、兼顾灌溉及供水的 一等大(1)型水电枢纽工程。大坝为混凝土面板堆石坝，最大坝高133m，计划于2004年竣工蓄水。图1和图 2分别给出了公伯峡面板堆石坝的标准横剖面和沿坝轴线的纵剖面。 收稿日期：2004-06-24 作者简介：张嘎(1976-)，山东人，博士，讲师，研究方向为岩土工程。Email:zhangga@tsinghua.edu.cn 1 水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期