›› 2012, Vol. 33 ›› Issue (6): 1896-1903.

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

裂隙网络无压渗流分析的初流量法

姚 池1,姜清辉1, 2,叶祖洋1,周创兵1, 2   

  1. 1. 武汉大学 水资源与水电工程科学国家重点实验室,武汉 430072;2. 武汉大学 土木建筑工程学院,武汉 430072
  • 收稿日期:2011-10-31 出版日期:2012-06-11 发布日期:2012-06-14
  • 通讯作者: 姜清辉,男,1972年生,博士,教授,主要从事岩土力学数值方法与边坡稳定分析方面的教学与科研工作 E-mail:122354748@99.com
  • 作者简介:姚池,男,1986年生,博士研究生,主要从事岩土力学数值方法、高坝复杂地基与岩石高边坡设计的研究工作
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2011CB013506);国家自然科学基金项目(No. 51079110)。

Initial flow method for unconfined seepage problems of fracture networks

YAO Chi1,JIANG Qing-hui1, 2,YE Zu-yang1,ZHOU Chuang-bing1, 2   

  1. 1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China; 2. School of Civil and Architectural Engineering, Wuhan University, Wuhan 430072, China
  • Received:2011-10-31 Online:2012-06-11 Published:2012-06-14

摘要: 借鉴初流量法的思想,引入初流速来抵消在裂隙网络干区实际不存在的流速,将Darcy定理扩展到整个区域,定义了在整个区域上的非线性边值问题,并将潜在溢出面边界条件归纳为Signorini型边界条件,建律了等价的变分不等式提法。通过结合连续型的Heaviside函数,给出了基于变分不等式的初流量有限元迭代算法,研制开发了相应的计算程序,通过两个典型算例验证了本文算法在求解复杂裂隙网络渗流自由面的有效性和鲁棒性。

关键词: 裂隙网络, 自由面, 渗流, 罚函数, 变分不等式

Abstract: Based on the idea of the initial flow method, initial velocity is introduced to eliminate the velocity that does not exist in dry zone. Darcy formula is extended to the entire region and a nonlinear boundary value problem of the entire region is defined. Through categorizing the potential seepage face as a boundary condition of Signorini’s type, an equivalent variational inequality is established. By combining the continuous Heaviside function, a finite element iterative algorithm of the variational inequality based on initial flow method is given. With this algorithm, a program is developed. Two typical examples are employed to demonstrate the efficiency and robustness of the algorithm proposed here in locating the free surface of complex fracture networks subject to underground water.

Key words: fracture network, free surface, seepage, penalty function, variational inequality

中图分类号: 

  • O 357.3
[1] 苏永华, 李诚诚. 强降雨下基于Green-Ampt模型的边坡稳定性分析[J]. 岩土力学, 2020, 41(2): 389-398.
[2] 刘成禹, 陈博文, 罗洪林, 阮家椿, . 满流条件下管道破损诱发渗流侵蚀的试验研究[J]. 岩土力学, 2020, 41(1): 1-10.
[3] 刘忠玉, 夏洋洋, 张家超, 朱新牧. 考虑Hansbo渗流的饱和黏土 一维弹黏塑性固结分析[J]. 岩土力学, 2020, 41(1): 11-22.
[4] 夏才初, 喻强锋, 钱 鑫, 桂 洋, 庄小清. 常法向刚度条件下岩石节理剪切−渗 流特性试验研究[J]. 岩土力学, 2020, 41(1): 57-66.
[5] 支永艳, 邓华锋, 肖瑶, 段玲玲, 蔡佳, 李建林. 微生物灌浆加固裂隙岩体的渗流特性分析[J]. 岩土力学, 2019, 40(S1): 237-244.
[6] 许江, 邬君宇, 刘义鑫, 雷娇, . 不同充填度下岩体剪切−渗流耦合试验研究[J]. 岩土力学, 2019, 40(9): 3416-3424.
[7] 刘伟俊, 张晋勋, 单仁亮, 杨昊, 梁辰, . 渗流作用下北京砂卵石地层多排管局部 水平冻结体温度场试验[J]. 岩土力学, 2019, 40(9): 3425-3434.
[8] 李博, 黄嘉伦, 钟振, 邹良超, . 三维交叉裂隙渗流传质特性数值模拟[J]. 岩土力学, 2019, 40(9): 3670-3768.
[9] 张强, 李小春, 周英博, 石露, 白冰, . 高压孔隙CO2/水作用下完整四川三叠系 砂岩剪切特性的试验研究[J]. 岩土力学, 2019, 40(8): 3028-3036.
[10] 王鹏飞, 谭文辉, 马学文, 李子建, 刘景军, 武洋帆, . 不同粗糙度和隙宽贯通充填裂隙 渗流特性试验研究[J]. 岩土力学, 2019, 40(8): 3062-3070.
[11] 张天军, 庞明坤, 蒋兴科, 彭文清, 纪翔, . 负压对抽采钻孔孔周煤体瓦斯渗流特性的影响[J]. 岩土力学, 2019, 40(7): 2517-2524.
[12] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[13] 刘忠玉, 崔鹏陆, 郑占垒, 夏洋洋, 张家超. 基于非牛顿指数渗流和分数阶Merchant模型的 一维流变固结分析[J]. 岩土力学, 2019, 40(6): 2029-2038.
[14] 刘 健, 陈 亮, 王春萍, 马利科, 王 驹. 一种非稳态气体渗流条件下岩石渗透特性 参数计算方法及应用[J]. 岩土力学, 2019, 40(5): 1721-1730.
[15] 徐 强, 肖 明, 陈俊涛, 倪少虎, . 渗流监测数据缺失处理与渗透稳定判断[J]. 岩土力学, 2019, 40(4): 1526-1534.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘 晓,唐辉明,刘 瑜. 基于集对分析和模糊马尔可夫链的滑坡变形预测新方法研究[J]. , 2009, 30(11): 3399 -3405 .
[2] 胡大伟,周 辉,谢守益,张 凯,邵建富,冯夏庭. 大理岩破坏阶段Biot系数研究[J]. , 2009, 30(12): 3727 -3732 .
[3] 师旭超,韩 阳. 卸荷作用下软黏土回弹吸水试验研究[J]. , 2010, 31(3): 732 -736 .
[4] 朱建明,彭新坡,姚仰平,徐金海. SMP准则在计算煤柱极限强度中的应用[J]. , 2010, 31(9): 2987 -2990 .
[5] 原喜忠,李 宁,赵秀云,杨银涛. 东北多年冻土地区地基承载力对气候变化敏感性分析[J]. , 2010, 31(10): 3265 -3272 .
[6] 白 冰,李小春,石 露,唐礼忠. 弹塑性应力-应变曲线的斜率恒等式及其验证和应用[J]. , 2010, 31(12): 3789 -3792 .
[7] 唐利民. 地基沉降预测模型的正则化算法[J]. , 2010, 31(12): 3945 -3948 .
[8] 李占海,朱万成,冯夏庭,李绍军,周 辉,陈炳瑞. 侧压力系数对马蹄形隧道损伤破坏的影响研究[J]. , 2010, 31(S2): 434 -441 .
[9] 蔡辉腾,危福泉,蔡宗文. 重庆主城区粉质黏土动力特性研究[J]. , 2009, 30(S2): 224 -228 .
[10] 宋 玲 ,刘奉银 ,李 宁 . 旋压入土式静力触探机制研究[J]. , 2011, 32(S1): 787 -0792 .