有自由面渗流分析三维数值流形方法.pdfVIP

第32 卷第3 期 岩 土 力 学 Vol.32 No.3 2011 年3 月 Rock and Soil Mechanics Mar. 2011 文章编号：1000－7598 (2011) 03 －879－06 有自由面渗流分析的三维数值流形方法 1, 2 2 2 姜清辉 ，邓书申 ，周创兵 （1. 武汉大学 土木建筑工程学院武汉 430072 ；2. 武汉大学 水资源与水电工程科学国家重点实验室，武汉430072 ） 摘 要：提出了求解有自由面渗流问题的三维数值流形方法，通过构造任意形状流形单元的水头函数，推导了流形单元的渗 透矩阵和无压渗流分析的总体控制方程，并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表 明，该方法采用数学网格覆盖整个材料区域，在自由面的迭代求解过程中数学网格保持不变，只考虑自由面以下渗流区的介 质，只对自由面以下的流形单元形成总体渗透矩阵，具有精度高、收敛速度快、编程简单等优点，而且能够通过单纯形积分 精确计算被自由面穿越单元的渗透作用力，因此，特别适用于有自由面渗流问题的模拟。 关 键 词：三维渗流；数值流形方法；自由面；数学网格；流形单元 中图分类号：O342 文献标识码：A Three-dimensional numerical manifold method for seepage problems with free surfaces JIANG Qing-hui1,2 2 2 , DENG Shu-shen , ZHOU Chuang-bing (1.School of Civil and Architectural Engineering, Wuhan University, Wuhan, 430072, China; 2.State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan, 430072, China) Abstract: A three-dimensional numerical manifold method for seepage problems with free surfaces is proposed. The hydraulic potential functions for arbitrarily shaped manifold element are constructed and the element conductivity matrix is derived in detail. The global governing equations for unconfined seepage analysis are established by minimizing the flow dissipation energy. The proposed method employs the tetrahedral mathematical meshes to cover the whole material volume. In the process of iterative solving for locating the free surface, the numerical manifold method can stric