光照碾压混凝土坝坝基三维渗流场的有限元精细求解.DOCVIP

光照碾压混凝土坝坝基三维渗流场的有限元分析 PAGE 6 PAGE 5 光照碾压混凝土坝坝基三维渗流场的有限元分析 作者简介：吴世勇，1982年出生，男，汉族，安徽省肥东县人，硕士研究生，主要研究高坝及地下结构的渗控计算与分析。 光照碾压混凝土坝坝基三维渗流场的有限元分析 吴世勇，朱岳明 （河海大学水利水电工程学院，南京，210098） 摘要：本文利用固定网格结点虚流量法，通过对坝基岩体内复杂排水系统的精细模拟，详细分析了光照坝基以及两岸岩体内的渗流状态和渗控效果，对光照碾压混凝土坝基础的渗控措施作出评价。 关键词：光照水电站；坝基；渗流场；渗控措施；有限元分析 中国分类号：TV 139.1 文献标志码：A 1 工程概况 北盘江光照水电站位于贵州省关岭县与晴隆县交界的北盘江中游，是北盘江干流梯级的龙头水电站，电站以发电为主，其次为航运，兼顾灌溉、供水及其他综合效益。正常蓄水位745m时，相应库容31.35亿m3，调节库容20.37亿m3，为不完全多年调节水库。电站总装机容量1040MW，保证出力180.2MW，年发电量27.54亿kW.h。枢纽工程主要由碾压混凝土重力坝、坝身泄洪系统、右岸引水系统、地面厂房和左岸预留远景通航建筑物等组成。 坝址区岸坡平直，地形完整，两岸地形较为对称，为典型的“V”型横向河谷，山体雄厚，河谷狭窄，岸坡上缓下陡，坡角约45°。坝区断层不发育，仅有F1和F2两条逆断层，其中F1规模稍大，并且穿过坝基，其影响范围约20~30m宽。坝区地下水靠大气降水分散式补给或洼地落水洞集中式补给，两岸透水岩组地下水位均高于或接近同期河水位，地下水以岩溶管道水、岩溶泉等形式排向北盘江。 光照水电站碾压混凝土重力坝最大坝高达200.50m，最大底宽159.05m，顶宽12m，坝顶长410m，大坝由左、右岸挡水坝段和河床溢流坝段组成。大坝建成并在水库蓄水后，在库水位以及各种渗控措施作用下，坝基岩体的渗流状态将发生了很大的变化，为了了解大坝基础的渗流状态和各种渗控措施的作用，笔者进行了光照水电站大坝基础岩体的三维渗流场有限元分析。 计算理论 根据广义达西定理，各向异性连续介质中的三维稳定渗流的控制方程和边界条件为 （1） 其中 为坐标，；为达西渗透系数张量；，为总水头，为位置水头，为压力水头；为源汇项。 边界条件有 （2） （3） （4） （5） （6） 其中 为已知水头函数；为边界外法线方向余弦，；分别为第一类渗流边界、第二类渗流边界、渗流自由面和渗流逸出面；为法向流量，流出为正；为第个溢出型排水孔；为第个溢出型排水孔的内部边界；为溢流型排水孔编号，＝1，2，……，M，其中M为溢流型排水孔总数；为第个溢流型排水孔的所有单元临空面之和；为第个溢流型排水孔的孔顶高程。各种类型排水孔见图1。详细的计算方法请参照文献[1][2][3][4]和[5]。 计算及成果分析 渗流计算建模及单元网格剖分 本次计算采用工程类比法选定了坝体材料的渗透系数；不同岩基的渗透系数则参照业主提供的压水试验结果选取；坝基断层F1与F2的影响带都取10m宽，渗透系数参照压水试验的结果（见表1） 由于资料有限，所以三维整体渗流场计算域的上、下游侧面、两岸山体侧面所取的范围达不到坝高的2倍，至于计算域底面到河床坝段的建基面的距离有最大坝高的2.5倍左右。底面及上、下游两个侧面视为不透水面，由于坝体的渗透系数比岩基的小两个数量级左右，即岩基于坝体的渗透交换水量非常小，可以忽略不计，所以坝体建基面也视为不透水面。计算模型坐标系的X轴与坝纵轴线重合，指向左岸；Y轴与坝横轴0+0.000重合，指向上游；Z轴垂直于XOY平面，正向竖直向上，原点落在工程零高程面上。用空间八结点六面体单元并辅以少量的五面体单元对计算域进行剖分，图2为渗流计算模型与单元网格图，图3为坝基防渗帷幕的单元剖分图，图4为河床及两岸坝基排水孔网格剖分图，其中逸流型和溢流型排水孔分别有120和508个。计算网格共有结点和单元数分别是63584和51348个，解题规模相当的大。本次计算上、下游面取正常水位，分别是745.00m和583.50m；左、右岸地下水位分别取636.00m和753.40m。 图2 渗流计算模型与单元网格剖分 图3 坝基防渗帷幕的单元