某心墙堆石坝二维平面有限元渗流分析.pdfVIP

科技信息 。建筑与工程o SCIENCETECHNOLOGYINFORMATION 2012年 第 35期 某心墙堆石坝二维平面有限元渗流分析 蒋德勇 杨培章z 1．四川华电金川水电开发有限公司 四川 成都 610041； 2_四川华电木里河水电开发有限公司 EIJII成都 610016 摘【 要】渗流问题是影响土石坝安全最为重要的因素，因而渗流问题的研究具有重要的理论意义和工程应用价值。本文结合土石坝相关 理论和计算机语言，编制了土石坝二维平面渗流有限元计算程序，并以某高心墙堆石坝为例，运用实际监测值验证计算的渗流参数正确性。 关【键词】渗流；土石坝；监测；程序 0 引言 表 1 堆石坝各结构的材料渗透系数 渗透系数 土石坝是我国历史最为悠久的一种坝型 ．该坝型具有施工简单、 造价低、速度快、对 自然条件适应性强等优点，因而备受坝工界青睐。 Kx cm／s Ky cm／s 堆石坝包括混凝土面板堆石坝和土质心墙堆石坝．其中前者主要是运 上游堆石 5．0E—02 5oE-02 用混凝土面板控制渗流．而后者是通过弱透水心墙以及防渗墙等措施 下游主堆石 5．0E-02 5．0E-02 来解决渗流问题 渗流问题是影响土石坝工程安全的重要因素．尤其 下游次堆石 1．0E一01 1．0E-01 是高土石坝 我国百米以上高度的土石坝已超过20座 ．如双江I I 高 反滤层 1．0E一02 1．OE-02 314m 、古水 高 305m 、两河口 高292m 、糯扎渡 高261．5m 、水布 垭 高232m 等。 过渡层 1．OE-02 1．OE-02 心墙 1．0E-07 1．0E-07 1 渗流分析原理 混凝土防渗墙 5．OE-07 5．0E-07 根据达西定律和水流连续性方程．不考虑土体压缩性时．其土体 覆盖层 1 1．OE一05 1．0E-05 各向异性连续介质稳定流域n内的水头函数 H满足如下偏微分方 覆盖层 2 1．OE-03 1．OE—o3 程 ： 覆盖层 3 1．0E-05 2．OE-06 軎 +ok， +0k 0 1 覆盖层 4 1．0E—O5 2．OE-06 式中：H为总水头；为x方向渗透系数；为Y方向渗透系数；： 4 计算成果分析 为z方向渗透系数 水头边界条件为： 本文根据土石坝渗流有限元分析原理和土体的相关特性 ．运用 日IFi H ，Y，Z，t 2 FORTRAN语言编制了土石坝渗流平面有限元计算程序 根据堆石坝 流量边界为： 的实际情况选取两种水位工况．进行该土石坝典型剖面的平面渗流有 限元分析。坝体内部埋设渗压计，取典型剖面上的5个典型测点，进行 一 -- k lF--q 3 2 验证程序计算值的准确性 根据程序分析堆石坝典型剖面渗流．得出 式中：r为上下游及渗出面边界之和： 为不透水边界、潜流边 相应节点单元各个工况下的压力水头值。 界、补给边界、自由面边界等之和；n为r的外法线方向； 为n方向 分析各工况下典型测点水头的计算值与监测差值和相应水头的 的渗透系数 变化规律可以看出．两种水位工况下典型剖面测压计的监测值与计算 对于无压渗流的自由面边界除应满足上式外，还需满足：H z。 值差值最大为O．97m．占坝高的O．78％ 坝高 125m 。其结果表明程序 的计算值能较好的与监测值吻合．且各测点水头的变化规律基本一 2 工程概况 致。另外由于该堆石坝的地质条件十分复杂，地下水丰富，有限元计算 模型的无法真实的模拟其结构．另外本文是分析堆石坝平面渗流．未 某心墙堆石坝为一个水库工程的挡水建筑物．最大坝高 125m，上 考虑平面以外的土体的影响．故存在一定的误差．符合堆石坝的有限 游坝坡 1：2．0，下游坝坡 1：1．8。该工程地质条件较复杂，左岸为石英闪 元渗流的一般规律 长岩．山体雄厚．基岩面自左岸向右岸向下深切 ．河床及右岸由中、上 更新统卵砾石层、粉质壤土和块碎石土夹硬质土层组成．大坝主要持 5 结语 力层为第三岩组[Q32—1 Ⅲ ]卵砾石层与粉质壤土互层，河床部份为第 本文首先介绍了土石坝渗流的相关理论 ．并结合计算机语言 二岩组fQ31 II 1块碎石土夹硬质土层，硬质土层约占该岩组总厚度的 FORTRAN编制了土石坝平面有限元渗流计算程序 以某心墙堆石坝 20～40％．此层为坝基深部相对隔水层 为了解决该堆石坝的渗流问 为例．运用程序分析不同工况下典型剖面的渗流问题．得出堆石坝相 题 ．采用心墙和混凝土防渗墙以及帷幕灌浆。 应的压力水头值．绘制其压力等值线图。另外取5个典型的渗流监测 3 计算模型及参数 点．分析其各个工况下压力水头的计算值与监测值的差异，以及相应 的变化规律 从分析的结果看．其程序计算的压力水头值基本正确