›› 2011, Vol. 32 ›› Issue (3): 879-884.

• 数值分析 • 上一篇    下一篇

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

姜清辉1, 2,邓书申2,周创兵2   

  1. 1. 武汉大学 土木建筑工程学院武汉 430072;2. 武汉大学 水资源与水电工程科学国家重点实验室,武汉430072
  • 收稿日期:2009-12-27 出版日期:2011-03-10 发布日期:2011-03-21
  • 作者简介:姜清辉,男,1972年生,博士,教授,主要从事岩土力学数值计算方法与岩土工程稳定分析方面的教学与研究工作。
  • 基金资助:

    国家自然科学基金项目(No. 50839004,No.51079110 );教育部新世纪优秀人才支持计划项目(No. NCET-07-0632)。

Three-dimensional numerical manifold method for seepage problems with free surfaces

JIANG Qing-hui 1,2 , DENG Shu-shen2 , ZHOU Chuang-bing2   

  1. 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
  • Received:2009-12-27 Online:2011-03-10 Published:2011-03-21

摘要:

提出了求解有自由面渗流问题的三维数值流形方法,通过构造任意形状流形单元的水头函数,推导了流形单元的渗透矩阵和无压渗流分析的总体控制方程,并给出了自由面的迭代求解策略和渗透体积力的计算方法。典型算例的数值分析表明,该方法采用数学网格覆盖整个材料区域,在自由面的迭代求解过程中数学网格保持不变,只考虑自由面以下渗流区的介质,只对自由面以下的流形单元形成总体渗透矩阵,具有精度高、收敛速度快、编程简单等优点,而且能够通过单纯形积分精确计算被自由面穿越单元的渗透作用力,因此,特别适用于有自由面渗流问题的模拟。

关键词: 三维渗流, 数值流形方法, 自由面, 数学网格, 流形单元

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 strictly realize the seepage analysis of the saturated domain on the condition of mathematical meshes keeping unchanged. Furthermore, the seepage force acting on the transitional elements cut by the free surface can be accurately calculated by the manifold method. Therefore, the proposed method is featured in high accuracy, fast convergence rate and simple programming, especially applicable to simulate the unconfined seepage problem with free surface.

Key words: three-dimensional seepage, numerical manifold method, free surface, mathematical mesh, manifold element

中图分类号: 

  • O342
[1] 刘登学, 张友良, 丁秀丽, 黄书岭, 裴启涛, . 数值流形法中基于适合分析T样条的 局部网格加密算法[J]. 岩土力学, 2019, 40(4): 1584-1595.
[2] 韩智铭,乔春生,朱 举. 含两组交叉贯通节理岩体的强度及破坏特征分析[J]. , 2018, 39(7): 2451-2460.
[3] 杨石扣,任旭华,张继勋. 基于布尔交运算的三维流形单元生成研究[J]. , 2016, 37(9): 2706-2711.
[4] 徐杨青,刘国锋,吴西臣. 武汉长江Ⅰ级阶地含水层水文地质参数研究及工程应用[J]. , 2016, 37(10): 2915-2920.
[5] 刘登学,张友良,谭 飞,张礼仁. 六面体覆盖的高阶数值流形方法的探讨[J]. , 2014, 35(S2): 662-666.
[6] 刘泉声,刘学伟. 多场耦合作用下岩体裂隙扩展演化关键问题研究[J]. , 2014, 299(2): 305-321.
[7] 李树忱 ,王兆清 ,袁 超 . 岩土体渗流自由面问题的重心插值无网格方法[J]. , 2013, 34(7): 1867-1873.
[8] 叶祖洋 ,姜清辉 ,姚 池 ,周创兵 . 岩体裂隙网络非稳定渗流分析与数值模拟[J]. , 2013, 34(4): 1171-1178.
[9] 姚 池 ,姜清辉 ,叶祖洋 ,周创兵 . 裂隙网络无压渗流分析的初流量法[J]. , 2012, 33(6): 1896-1903.
[10] 刘 建,陈 佺. 一种模拟岩石蠕变的数值流形方法[J]. , 2012, 33(4): 1203-1209.
[11] 李 毅 ,伍 嘉 ,李 坤. 基于FLAC3D的饱和-非饱和渗流分析[J]. , 2012, 33(2): 617-622.
[12] 郝亚飞,李海波,刘恺德,杨风威,夏 祥,邹 飞,朱小明. 单自由面爆破振动特征的炮孔堵塞长度效应[J]. , 2011, 32(10): 3105-3110.
[13] 谭 飞,邓 帮,张东明,王元汉. 杂交边界点法求解带自由面的渗流问题[J]. , 2010, 31(S2): 75-78.
[14] 焦 健,乔春生,徐干成. 开挖模拟在数值流形方法中的实现[J]. , 2010, 31(9): 2951-2957.
[15] 周小义,邓安福. 六面体有限覆盖的三维数值流形方法的非线性分析[J]. , 2010, 31(7): 2276-2282.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 肖云华,王 清,陈剑平. 基于优化技术的权重计算方法在岩体质量评价中的应用[J]. , 2009, 30(9): 2686 -2690 .
[2] 张鸿飞,程效军,高 攀,周鑫鑫. 隧道衬砌空洞探地雷达图谱正演模拟研究[J]. , 2009, 30(9): 2810 -2814 .
[3] 张安康,陈士海,杜荣强,魏海霞. 岩石类材料的能量基率相关弹塑性损伤模型[J]. , 2010, 31(S1): 207 -210 .
[4] 王小军,屈耀辉,魏永梁,杨印海,达益正. 郑西客运专线湿陷性黄土区试验路堤的沉降观测与预测研究[J]. , 2010, 31(S1): 220 -231 .
[5] 陈 瑜,曹 平,蒲成志,刘业科,李 娜. 水-岩作用对岩石表面微观形貌影响的试验研究[J]. , 2010, 31(11): 3452 -3458 .
[6] 李明超,王忠耀,刘 杰. 水库库岸滑坡体稳定性与三维可视化模拟分析系统研究[J]. , 2009, 30(1): 270 -274 .
[7] 赵延喜,徐卫亚. 基于AHP和模糊综合评判的TBM施工风险评估[J]. , 2009, 30(3): 793 -798 .
[8] 张其一,栾茂田. 复合加载情况下非均质地基上条形基础的极限承载力研究[J]. , 2009, 30(5): 1281 -1286 .
[9] 王俊卿,李 靖,李 琦,陈 立. 黄土高边坡稳定性影响因素分析 ——以宝鸡峡引水工程为例[J]. , 2009, 30(7): 2114 -2118 .
[10] 龚彦峰,张俊儒. 隧道单层衬砌设计方法研究及应用[J]. , 2011, 32(4): 1062 -1068 .