岩土力学 ›› 2024, Vol. 45 ›› Issue (8): 2483-2491.doi: 10.16285/j.rsm.2023.1397

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

三维有限元中奇异点源的简化分析方法

潘泓1, 2,许嘉显1, 2,骆冠勇1, 2,彭斯格1, 2,曹洪1, 2   

  1. 1. 华南理工大学 土木与交通学院,广东 广州 510640;2. 华南理工大学 亚热带建筑科学国家重点实验室,广东 广州 510640
  • 收稿日期:2023-09-18 接受日期:2023-12-04 出版日期:2024-08-10 发布日期:2024-08-12
  • 通讯作者: 骆冠勇,男,1979年生,博士,副教授,主要从事岩土工程方面的教学与研究工作。E-mail: luogy@scut.edu.cn
  • 作者简介:潘泓,男,1967年生,博士,教授,主要从事岩土工程方面的教学与研究工作。E-mail: hpan@scut.edu.cn
  • 基金资助:
    国家自然科学基金(No. 52308353,No. 51978282);广东省自然科学基金(No. 2023A1515011571,No. 2023A1515011683);华南理工大学亚热带建筑科学国家重点实验室开放研究项目(No. 2022ZB21)

Simplified analysis method of singular point source in three-dimensional finite element calculation

PAN Hong1, 2, XU Jia-xian1, 2, LUO Guan-yong1, 2, PENG Si-ge1, 2, CAO Hong1, 2   

  1. 1. School of Civil Engineering and Transportation, South China University of Technology, Guangzhou, Guangdong 510640, China; 2. State Key Laboratory of Subtropical Building Science, South China University of Technology, Guangzhou, Guangdong 510640, China
  • Received:2023-09-18 Accepted:2023-12-04 Online:2024-08-10 Published:2024-08-12
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (52308353, 51978282), the Natural Science Foundation of Guangdong Province (2023A1515011571, 2023A1515011683) and the Open Fund of State Key Laboratory of Subtropical Building Science, South China University of Technology (2022ZB21).

摘要: 减压井是渗流控制的重要手段,减压浅井井底的入渗量不可忽略。在渗流场有限元有关井的计算中,含水层与井半径有巨大尺寸差异,只有将井周围的网格尺寸限制在厘米级或更少,才能达到合理的精度,但效率过低。为提高效率,分析时通常将井底简化为点源。由于点源附近的水头梯度极大,是奇异点,有限元中通常要进行等效处理才能保证全局的精度。提出了以零体积点元模拟点源的方法,该方法基于修正井水位法,利用点源解析解及四面体渗流单元,推导了点元单元的显式解。使用该点元单元模拟井底,可以在不失全局精度的条件下用大网格对点源进行模拟,避免了迭代计算,极大提高了利用三维有限元模拟点源的效率。详细介绍了该点元显式解的推导过程,验证其计算精度,并探究出简化计算方法的适用范围。最后介绍了该方法在实际排水减压抗浮工程中的应用。

关键词: 渗流, 有限元法, 井点, 奇异点单元

Abstract: Relief wells are an important means of seepage control, and the infiltration at the bottom of shallow relief wells cannot be ignored. When utilizing the finite element method for seepage field analysis in relation to wells, a significant disparity exists between the aquifer's dimensions and the well's radius. Achieving reasonable accuracy necessitates limiting the mesh size around the well to centimeters or less, albeit at the cost of reduced efficiency. To enhance efficiency, the well bottom is commonly approximated as a point source in analytical models. The head gradient near the point source is extremely large, which is a singular point. In the finite element method, equivalent treatment is usually required to ensure the global accuracy. This paper presents a method for simulating point sources with zero-volume point elements which is based on the correction well water level method, the method uses the point source analytical solution and the tetrahedral seepage element to deduce the explicit solution of the point source element. By using this point source element to simulate the bottom of well, the point source can be simulated by using a large grid without losing the global accuracy, which greatly improves the efficiency of simulating the point source by using the three-dimensional finite element. The paper thoroughly outlines the derivation process for the explicit solution of the point element, verifies its calculation accuracy, and investigates the applicability range of simplified calculation methods. Lastly, the practical application of this method in drainage decompression and anti-floating engineering is presented.

Key words: seepage, finite element method, well point, singular point element

中图分类号: TU 46+3
[1] 张宪成, 池宝涛, 于先泽, 郭前建, 袁伟, 张耀明, . 面向近场有限元法的非结构网格划分及其断裂损伤分析[J]. 岩土力学, 2025, 46(S1): 467-476.
[2] 江文豪, 王浩, 廖光志, 陈滨华, . 时变降雨场景下双层非饱和土中水分一维瞬态渗流解析解[J]. 岩土力学, 2025, 46(9): 2721-2737.
[3] 王昕琪, 冯子军, 陈正男, 高祺, 阴伟涛, 靳佩桦, 李玉彬, . 超临界水作用下花岗岩裂隙渗流特性演化规律研究[J]. 岩土力学, 2025, 46(9): 2847-2858.
[4] 张弛, 邓龙传, 庄欠伟, 李晓昭, 王秋平, 乔梁, . 基于试验和数值模拟的滚刀旋转破岩受力和效率研究[J]. 岩土力学, 2025, 46(9): 2995-3006.
[5] 靳贵晓, 林劭聪, 姜启武, 黄明, 李熹, . 基于Kozeny-Carman方程的酶诱导碳酸钙沉淀固化砂土的渗流数学模型[J]. 岩土力学, 2025, 46(8): 2376-2386.
[6] 吴兵, 盛建龙, 叶祖洋, 周新, . 单裂隙岩体非线性渗流-法向应力耦合模型研究[J]. 岩土力学, 2025, 46(8): 2495-2504.
[7] 陈灯红, 张心瀚, 刘云辉, 胡昊文, 刘云龙, 梁羽翔, . 基于八叉树比例边界有限元法的高拱坝-不规则地基-库水系统非线性地震响应分析[J]. 岩土力学, 2025, 46(8): 2586-2599.
[8] 张海燕, 胡新丽, 李亚博. 周期性渗流-应力耦合下滑带土的蠕变特性试验研究[J]. 岩土力学, 2025, 46(7): 2189-2198.
[9] 任旺, 苗君, 雷万钧, 王络, 陈益峰, . 混凝土坝坝基立体排水的数值模拟与优化设计[J]. 岩土力学, 2025, 46(7): 2281-2295.
[10] 刘杰, 钟伦伟, 钟振, 张宪尚, 隆清明, 应鹏, . 恒定法向刚度边界条件下受剪岩石裂隙非线性渗流特性的数值模拟研究[J]. 岩土力学, 2025, 46(6): 1919-1933.
[11] 邓汉鑫, 武东生, 王晓光, 王东坡. 非均质孔隙介质中循环驱替对两相渗流及毛细捕获的影响[J]. 岩土力学, 2025, 46(2): 505-514.
[12] 王浩, 侯泓冰, 江文豪, . 考虑根系对水文特性影响下非饱和植被土中水分一维瞬态渗流解析解[J]. 岩土力学, 2025, 46(11): 3329-3345.
[13] 张仁君, 肖碧, 杨志兵, 郑小康, 胡冉, 陈益峰. 裂隙介质中毛细黏聚作用影响颗粒运移-堵塞机制研究[J]. 岩土力学, 2025, 46(11): 3473-3484.
[14] 唐劲舟, 唐文豪, 杨科, 赵延林, 刘钦节, 段敏克, 谭哲, . 循环荷载作用下含倾斜单裂隙砂岩力学响应特征及渗流演化规律[J]. 岩土力学, 2025, 46(1): 199-212.
[15] 张瑨, 李书恒, 朱其志, 石玲玲, 邵建富, . 基于深度学习的岩石短长期本构模型及灰砂岩变形预测[J]. 岩土力学, 2025, 46(1): 289-302.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!