基于初流量法的子结构技术在渗流分析中的应用-中国水利水电科学.pdfVIP

第11卷 第4期 中国水利水电科学研究院学报 Vol.11 No.4 2013年12月 Journal of China Institute of Water Resources and Hydropower Research December，2013 文章编号：1672-3031（2013）04-0241-08 基于初流量法的子结构技术在渗流分析中的应用 1，2 1 1 李新强 ，王思洪 ，王京宇 （1. 水利部建设管理与质量安全中心，北京 100038；2. 中国水利水电科学研究院 岩土研究所，北京 100048） 摘要：水利水电工程中地下厂房或洞室往往采用密集排水孔的排水措施，以确保其长期稳定和安全运行。对其作 渗流计算时，由于渗流自由面的求解是一个非线性问题，而密集排水孔的模拟条件较为复杂，这两个问题往往交 织在一起，计算难度很大。本文在引入连续调整系数对初流量法进行改进的基础上，采用子结构技术，并将子结 构作为整体网格结构的一部分进行处理，统筹兼顾整体和局部（密集排水孔）渗流场，较好的解决了上述问题。并 将该方法应用于泵站地下主厂房的渗流场分析中，取得较好的效果。 关键词：岩石力学；渗流；子结构；初流量法 中图分类号：TU43 文献标识码：A 1 研究背景 有自由面的渗流数值计算，由于在计算中自由面边界是未知的，使得其计算变得复杂和困难。 最初采用移动网格法，后来国内外学者提出了有自由面渗流计算的新方法——固定网格法。固定网 ［1-2］ ［3］ ［4］ 格法可分为两类：一是以单元渗透矩阵调整法 、剩余流量法 、改进剩余流量法、初流量法 、变 ［5］ ［6］ ［7］ 分不等式法 、节点虚流量法 、虚单元法和截止负压法 为代表的饱和渗流模型；二是建立地下水 饱和非饱和联合渗流计算模型。其研究核心是采用扩大渗流区域和固定边界来求解，以达到求解过 程中单元网格不变的目的。初流量法具有只需一次形成和分解总体传导矩阵，并不需在每次迭代步 中确定自由面的近似位置和判别自由面与单元相交的实际情形等诸多优点，因此是目前不变网格求 解有自由面渗流问题的一种较为有效的方法。该方法得到普遍认可，但在实际应用中出现了自由面 ［8］ 较大变动的现象。为此，王媛 通过分析初流量法的基本原理，提出区域识别函数的新概念，并利用 ［9］ 连续的区域识别函数进行线性微调，使之具有较好的收敛性。李新强 对初流量法进行了仔细分析， 发现单元区域和有初流量域的积分误差是造成自由面较大变动的根本原因，提出了连续调整系数的 概念，较为完善地解决了计算收敛性问题。水利水电工程中地下厂房或洞室往往处于高水头作用 下，降低地下厂房或洞室所受的外水压力，是保证其长期稳定和安全运行最为可靠的手段，目前最 为经济有效的方法是采用密集排水孔排水措施。在用有限元法对排水系统模拟时，由于排水孔直径 小、密集程度高，全局网格剖分使得原始数据准备工作量大而繁琐，且难以考虑多种计算方案。有 限元子结构法具有灵活性强，可以在全局网格的基础上进行二次网格划分从而模拟细部结构，因 此，密集排水孔排水的模拟采用子结构技术是明智的选择。 在具有较高地下水位或较强入渗的岩石工程中，往往同时存在具有自由面求解和模拟密集排水 孔的渗流计算问题。本文采用改进初流量法求解有自由面的渗流场，并用排水子结构法模拟地下厂 房周围布设的密集排水孔，对某泵站工程地下厂房围岩渗流场进行有限元分析，研究了地下厂房围 岩不同渗控方案下的渗控效果。 收稿日期：2012-07-19 基金项目：中国水利水电科学研究院科研专项（岩集1027）；水利部公益性行业科研专项（201301021） 作者简介：李新强（1972-），男，山东滕州人，博士，高级工程师，主要从事水工渗流研究工作。E-mail：lixq@