岩土力学 ›› 2021, Vol. 42 ›› Issue (12): 3428-3439.doi: 10.16285/j.rsm.2021.1319

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

准三维有限元模型在多层地基区域渗流中 的扩展与适宜性分析

梅锦玲1, 2,曹洪1, 2,骆冠勇1, 2,潘泓1, 2   

  1. 1. 华南理工大学 土木与交通学院,广东 广州 510640;2. 华南理工大学 亚热带建筑科学国家重点实验室,广东 广州 510640
  • 收稿日期:2021-08-11 修回日期:2021-09-25 出版日期:2021-12-13 发布日期:2021-12-14
  • 通讯作者: 潘泓,男,1967年生,博士,教授,主要从事岩土方法的教学与研究工作。E-mail: hpan@scut.edu.cn E-mail: jinlingmei@163.com
  • 作者简介:梅锦玲,女,1989年生,博士研究生,主要从事渗流及岩土工程方面的研究工作。
  • 基金资助:
    国家自然科学基金(No.51978282);广东省自然科学基金(No.2018A030313383,No.2020A1515010583)

Extension and applicability analysis of a quasi-three-dimensional finite element model to simulate regional groundwater flow problems in multi-layered subsoils

MEI Jin-ling1, 2, CAO Hong1, 2, LUO Guan-yong1, 2, PAN Hong1, 2   

  1. 1. School of Civil Engineering and Transportation, South China University of Technology, Guangzhou, Guangdong 510640, China; 2. The State Key Laboratory of Subtropical Building Science, South China University of Technology, Guangzhou, Guangdong 510640, China
  • Received:2021-08-11 Revised:2021-09-25 Online:2021-12-13 Published:2021-12-14
  • Supported by:
    This work was supported by the National Natural Science Foundation of China(51978282) and the Natural Science Foundation of Guangdong Province (2018A030313383, 2020A1515010583).

摘要: 在多层地基区域地下水流动的数值模拟中,假设地下水在强透水层内为水平流动、在弱透水层内为垂向流动的传统准三维方法计算效率高,但多适用于强弱透水层隔层分布的多层情况,不利于推广。对准三维计算理论进行扩展以适用于一般的多层情况,推导了各层土层内以平均水头为控制变量、且考虑上下界面入渗的平面二维流动控制方程,利用界面处的流速连续条件耦合各层控制方程,并在各土层内对平均水头与上下界面水头进行二次插值获取水头分布的垂向变化。基于扩展的计算理论建立相应的准三维有限元模型,各土层内划分相同的水平面二维网格,并在各节点处建立流量平衡关系。经算例验证,扩展后的准三维方法打破了传统方法的局限,允许各土层的水平和垂向流动同时得到合理的模拟,具有更高的精度和通用性,并且保持了计算效率的优势。由于方便求解自由面问题和处理地下结构奇异问题,且可以获取较高精度的三维水头分布结果,该方法在分析多层地基区域地下水流动问题中可以成为一个高效和有竞争力的方案。

关键词: 渗流, 数值模拟, 多层地基, 区域地下水流动, 准三维模型

Abstract: The conventional quasi-three-dimensional technique has a high efficiency in the simulation of regional groundwater flow in multilayered subsoils, with the assumption that groundwater flow is essentially horizontal in the aquifers and vertical in the aquitards. This assumption requires a two orders of magnitude permeability contrast between the adjacent soils, which is not conducive to a wide application of the conventional technique. In this study, calculation theory is extended to be applied to regional groundwater flow in general multi-layered subsoils. Firstly, the partial differential equation governing two-dimensional horizontal flow in each soil is derived with the average hydraulic head as the main variable and considering the infiltration at the upper and lower interfaces. The velocity continuity conditions at the soil interfaces provide the coupling condition among the soils. Then, the vertical variation of the hydraulic head in individual soil is obtained by quadratic interpolation between the average head and the head at the upper and lower interfaces. A quasi-three-dimensional finite element model is presented based on the extended calculation theory. Each soil is discretized with the same two-dimensional horizontal grids, and the fluid mass balance is established at each node. The extended quasi-three-dimensional technique has a higher accuracy and greater applicability, as it can overcome the limitations of the conventional technique, and simultaneously get a reasonable simulation of the horizontal and vertical flow in each layer without sacrificing the calculation efficiency. The extended technique provides an efficient and competitive scheme to simulate regional groundwater flow in multilayered subsoils, because it is convenient to solve free surface problems and deal with singular problems of underground structures, and compute a high-precision three-dimensional representation of the hydraulic head distribution.

Key words: seepage, numerical simulation, multilayered subsoils, regional groundwater flow, quasi-three-dimensional model

中图分类号: 

  • TU 46
[1] 贺勇, 胡广, 张召, 娄伟, 邹艳红, 李星, 张可能, . 污染场地六价铬迁移转化机制与数值模拟研究[J]. 岩土力学, 2022, 43(2): 528-538.
[2] 魏天宇, 王旭宏, 吕涛, 胡大伟, 周辉, 洪雯. 湿化膨胀与掺砂率对混合型缓冲材料 THM耦合过程的影响分析[J]. 岩土力学, 2022, 43(2): 549-562.
[3] 侯晓萍, 樊恒辉. 基于COMSOL Multiphysics的非饱和 裂隙土降雨入渗特性研究[J]. 岩土力学, 2022, 43(2): 563-572.
[4] 马成昊, 朱长歧, 刘海峰, 崔翔, 王天民, 姜开放, 易明星, . 土的颗粒形貌研究现状及展望[J]. 岩土力学, 2021, 42(8): 2041-2058.
[5] 崔溦, 王利新, 江志安, 王超, 王枭华, 张社荣, . 基于修正立方定律的岩体粗糙裂隙网络 注浆过程模拟研究[J]. 岩土力学, 2021, 42(8): 2250-2258.
[6] 夏才初, 徐英俊, 王辰霖, 赵海鸥, 薛小代, . 基于非稳态渗流过程的压气储能 洞室空气渗漏率计算[J]. 岩土力学, 2021, 42(7): 1765-1773.
[7] 朱淳, 何满潮, 张晓虎, 陶志刚, 尹乾, 李利峰, . 恒阻大变形锚杆非线性力学模型 及恒阻行为影响参数分析[J]. 岩土力学, 2021, 42(7): 1911-1924.
[8] 赵海鹏, 黎学优, 万建宏, 郑翔之, 刘思威, . 基于高性能有限单元法的多层地基水平 受荷桩受力分析[J]. 岩土力学, 2021, 42(7): 1995-2003.
[9] 黄娜, 蒋宇静, 程远方, 刘日成, . 基于3D打印技术的复杂三维粗糙裂隙网络 渗流特性试验及数值模拟研究[J]. 岩土力学, 2021, 42(6): 1659-1668.
[10] 彭守建, 王瑞芳, 许江, 甘青青, 蔡果良, . 二次炭化温度对热压型煤力学性质 及微观结构影响的试验研究[J]. 岩土力学, 2021, 42(5): 1221-1229.
[11] 李跃, 徐卫亚, 易魁, 谢伟超, 张强, 孟庆祥, . 堆积体滑带土非饱和-饱和渗透特性试验研究[J]. 岩土力学, 2021, 42(5): 1355-1362.
[12] 王兆耀, 刘红军, 杨奇, 赵真, 胡瑞庚, . 波流作用下大直径单桩的局部冲刷特征分析[J]. 岩土力学, 2021, 42(4): 1178-1185.
[13] 雷华阳, 许英刚, 缪姜燕, 刘旭. 动渗耦合作用下软黏土动力特性试验研究[J]. 岩土力学, 2021, 42(3): 601-610.
[14] 陈猛, 崔秀文, 颜鑫, 王浩, 王二磊. 岩石−钢纤维混凝土复合层抗压强度预测模型[J]. 岩土力学, 2021, 42(3): 638-646.
[15] 石峰, 卢坤林, 尹志凯. 平移模式下刚性挡土墙三维被动滑裂面的确定与土压力计算方法研究[J]. 岩土力学, 2021, 42(3): 735-745.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 石玉玲,门玉明,彭建兵,黄强兵,刘洪佳. 地裂缝对不同结构形式桥梁桥面的破坏试验研究[J]. , 2009, 30(10): 2917 -2922 .
[5] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[6] 徐速超,冯夏庭,陈炳瑞. 矽卡岩单轴循环加卸载试验及声发射特性研究[J]. , 2009, 30(10): 2929 -2934 .
[7] 张力霆,齐清兰,魏静,霍倩,周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[8] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[9] 易 俊,姜永东,鲜学福,罗 云,张 瑜. 声场促进煤层气渗流的应力-温度-渗流压力场的流固动态耦合模型[J]. , 2009, 30(10): 2945 -2949 .
[10] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .