数值流形法渗流分析中处理弱不连续界面的新方法

贾真; 杨冬梅; 郑宏

doi:10.11988/ckyyb.20230742

PDF(4490 KB)

长江科学院院报 ›› 2024, Vol. 41 ›› Issue (12) : 133-137. DOI: 10.11988/ckyyb.20230742

岩土工程

数值流形法渗流分析中处理弱不连续界面的新方法

贾真 ¹ ,
杨冬梅 ² ,
郑宏 ¹

作者信息 +

A New Method to Deal with Weak Discontinuous Interface in Seepage Analysis Based on Numerical Manifold Method

Author information +

文章历史 +

摘要

数值流形法采用数学覆盖和物理覆盖的双覆盖系统,具有灵活处理边界问题、高效进行网格划分以及便捷提高近似精度等优点,是一种很有前景的数值方法。不同于传统数值流形法根据界面来切割数学覆盖形成物理覆盖的做法,数值流形法基于弱不连续物理覆盖,在流行单元中利用折射定律构造出一种新的权函数,以此建立局部近似,并将其应用在稳定渗流问题中。通过对典型算例的计算分析,结果表明该方法在解决不连续界面问题中具有准确性和便利性。

Abstract

The numerical manifold method is a promising numerical method which uses a dual coverage system consisting of both mathematical and physical coverages. It has the advantages of flexible handling of boundary problems, efficient meshing and convenient improvement of approximation accuracy. Different from the traditional numerical manifold method which cuts the mathematical coverage to form the physical coverage according to the interface, a new weight function is constructed by using the refraction law in the manifold element to establish the local approximation based on the weak discontinuous physical coverage, and it is applied to the steady seepage problem. The accuracy and convenience of the method in solving discontinuous interface problems is demonstrated through the analysis of the calculation of typical seepage flow problems.

导出引用

贾真, 杨冬梅, 郑宏. 数值流形法渗流分析中处理弱不连续界面的新方法[J]. 长江科学院院报. 2024, 41(12): 133-137 https://doi.org/10.11988/ckyyb.20230742

JIA Zhen, YANG Dong-mei, ZHENG Hong. A New Method to Deal with Weak Discontinuous Interface in Seepage Analysis Based on Numerical Manifold Method[J]. Journal of Yangtze River Scientific Research Institute. 2024, 41(12): 133-137 https://doi.org/10.11988/ckyyb.20230742

中图分类号： TU43 (土力学)

参考文献

列表( 原文顺序 | 文献年度倒序 | 文中引用次数倒序 ) 可视化分析

[1]

崔皓东, 朱岳明, 张家发, 等. 深埋洞室群围岩渗流场分析及渗控效果初步评价[J]. 长江科学院院报, 2009, 26(10): 71-75.

摘要

采用三维有限元精细模拟技术，对某抽水蓄能电站工程复杂渗流场进行计算，重点研究了深埋高压输水隧洞和岔洞钢筋混凝土衬砌开裂状态下洞室群围岩的渗流场分布特征。研究表明，在高压岔洞下游侧高压支管区设置防渗帷幕和排水孔幕组合防渗体系，对衬砌开裂状态下的高压渗水有显著控制作用，使该区域渗控效果基本满足要求，为下游厂房长期安全运行提供了保障。

(CUI

Hao-dong

, ZHU

Yue-ming

, ZHANG

Jia-fa

, et al. Analysis on Seepage Field and Preliminary Evaluation on Seepage Control Effect of Rockmass around Deep Buried Caverns[J]. Journal of Yangtze River Scientific Research Institute, 2009, 26(10): 71-75.) (in Chinese)

本文引用 [1] 摘要

The complicated seepage flow field of a pumped storage power station project was simulated by 3D finite element refined simulation technology, the distribution characteristics of surrounding rock seepage field was studied in detail under lining cracks of the high pressure tunnel. The results show that the comprehensive measures for water tight curtain and drainage holes in downstream side of tunnel is effective under cracking hydraulic tunnel, so the seepage field can be controlled to reach up to the requirements, and the safety of the powerhouse is ensured in operating period.

[2]

张俊霞, 李莉, 朱登峰. 西霞院水库电站厂房坝段三维渗流计算分析[J]. 长江科学院院报, 2009, 26(10):67-70.

摘要

根据西霞院反调节水库电站厂房坝段工程的地质情况，该段工程的地基处理采用水平放置的“［”形混凝土防渗墙防渗。为了解厂房坝段地下水运动规律、渗流控制措施作用，对电站厂房坝段进行了三维渗流计算分析，结果认为：电站厂房段渗流控制措施的选择和布置，从整体上看比较合理；从扬压力等值线分布图来看，海漫段底板所承受的渗透压力较大。因此建议，对于海漫的砼底板厚度进行渐变加厚处理，亦即海漫段上游端至下游端，底板厚度由1.0m 渐变至0.60m ，从而确保电站厂房坝段的安全运行。

(ZHANG

Jun-xia

, LI

, ZHU

Deng-feng

. Three-dimensional Seepage Analysis of Powerhouse Dam Section of Xixia Institute Anti-regulation Reservoir[J]. Journal of Yangtze River Scientific Research Institute, 2009, 26(10):67-70.) (in Chinese)

本文引用 [1] 摘要

The Xixia Institute Anti-regulation Reservoir is a coordinated project of Yellow River Xiaolangdi Water Resources Project. According to the project geological situation of powerhouse dam section of the reservoir, this section project foundation treatment uses horizontal middle-brackets-like “［” shape concrete cut-off wall as anti seepage. For understanding the motion law of underground water and the control effect of seepage-proof measures, three dimensional seepage calculation analysis for the power plant workshop dam section is conducted. The analytic result shows as follows：The selection and arrangement of seepage control measures in the powerhouse dam section are reasonable overall. The bottom plate at the apron section withstands a rather great pressure, seen from the contour-line distribution diagram of uplifts. Therefore, it is suggested that the thickness of concrete bottom plate of the apron should be taken a gradual thickness from 1.0-0.6 meter to ensure the powerhouse dam section safe operation.

[3]

张伟, 许季军, 陈劲松. 江口水电站坝址区三维渗流计算分析[J]. 长江科学院院报, 2002, 19(2):31-33.

摘要

针对江口水电站坝址区地质条件，在水平和垂直方向同时进行了渗透性分区；采用三维有限元法，并考虑了岩体各向异性的特点和降雨入渗补给边界，对坝址区渗流场进行了计算分析。该计算模型能真实反映实际情况。通过计算，论证了现有防渗帷幕和排水设计方案的明显渗控效果，同时为设计提供了可靠的依据。

(ZHANG

Wei

, XU

Ji-jun

, CHEN

Jin-song

. Analysis on Three-dimensional FEM Seepage of Jiangkou Hydroelectric Station’s Dam Area[J]. Journal of Yangtze River Scientific Research Institute, 2002, 19(2):31-33.) (in Chinese)

本文引用 [1] 摘要

In order to simulate truly the underground water seeppage of Jiangkou Hydroelectric Station's dam area, the threedimensional seepage model was founded by 3D FEM under the conditions of the features of terrane anisotropism and supplementary boundaries of rainfall infiltration. The dam area is separated into several zones according to penetrability in horizontal and vertical directions. A series of models is considered. The results of these models show that the seepage control of proof curtain and draining system is efficient, and the data may be used as the basis for design.

[4]	BELYTSCHKO T, GRACIE R, VENTURA G. AReview of Extended/Generalized Finite Element Methods for Material Modeling[J]. Modelling and Simulation in Materials Science and Engineering, 2009, 17(4): 043001. 本文引用 [1]

[5]	FRIES T P, BELYTSCHKO T. The Extended/Generalized Finite Element Method: an Overview of the Method and Its Applications[J]. International Journal for Numerical Methods in Engineering, 2010, 84(3): 253-304. 本文引用 [1]

[6]

董玉文, 任青文, 苏琴. 接触摩擦问题的扩展有限元数值模拟方法[J]. 长江科学院院报, 2009, 26(5): 45-49.

摘要

研究了采用扩展有限元法分析接触摩擦问题的方法。该方法基于单位分解思想，在位移插值形式中引入不连续函数和裂尖渐近位移场函数分别反映接触面的不连续性和接触面端点的奇异性；根据虚功原理建立有限元支配方程，接触面不连续函数积分采用分区高斯积分的方法。扩展有限元法在分析接触问题时不需要引入接触面单元，接触面可以位于单元内任意位置。因此该方法在分析接触摩擦问题以及压剪裂纹的开裂扩展时有显著的优势。

(DONG

Yu-wen

, REN

Qing-wen

, SU

Qin

. Extended Finite Element Method of Frictional Contact Problem[J]. Journal of Yangtze River Scientific Research Institute, 2009, 26(5): 45-49.) (in Chinese)

本文引用 [1] 摘要

A new method for modeling frictional contact problem by extended finite element method (XFEM) is researched. On the basis of the idea of partition of unity, a discontinuous function and the two-dimensional asymptotic crack-tip displacement fields are added to the displacement interpolation formula to account for discontinuity of the contact surface and stress singularity near the contact surface tip respectively. Governing equations are deduced by virtual work principle and the integral of the discontinuous functions are calculated by Gauss integral through sub-elements cut by discontinuous surface. Thus it is unnecessary for using contact element in contact analysis by XFEM compared with FEM, and the element can be meshed without considering the location of the contact surface. So XFEM is promising in modeling frictional contact problem and the crack growth of compression-shear mixed mode crack.

[7]	STROUBOULIS T, BABUŠKA I, COPPS K. The Design and Analysis of the Generalized Finite Element Method[J]. Computer Methods in Applied Mechanics and Engineering, 2000, 181(1/2/3):43-69. 本文引用 [1]

[8]	STROUBOULIS T, COPPS K, BABUŠKA I. The Generalized Finite Element Method[J]. Computer Methods in Applied Mechanics and Engineering, 2001, 190(32/33): 4081-4193. 本文引用 [1]

[9]	李录贤, 刘书静, 张慧华, 等. 广义有限元方法研究进展[J]. 应用力学学报, 2009, 26(1):96-108,214. (LI Lu-xian, LIU Shu-jing, ZHANG Hui-hua, et al. Researching Progress of Generalized Finite Element Method[J]. Chinese Journal of Applied Mechanics, 2009, 26(1):96-108,214.) (in Chinese) 本文引用 [1]

[10]	MELENK J M, BABUŠKA I. Approximation with Harmonic and Generalized Harmonic Polynomials in the Partition of Unity Method[J]. Computer Assisted Mechanics and Engineering Sciences, 1997, 4(3/4): 607-632. 本文引用 [1]

[11]	郑宏. 数值流形法[M]. 北京: 科学出版社, 2022. (ZHENG Hong. Numerical Manifold Method[M]. Beijing: Science Press, 2022.) (in Chinese) 本文引用 [1]

[12]	AN X, MA G, CAI Y, et al. A New Way to Treat Material Discontinuities in the Numerical Manifold Method[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(47/48): 3296-3308. 本文引用 [1]

[13]	AN X, MA G, CAI Y, et al. Arbitrary Discontinuities in the Numerical Manifold Method[J]. International Journal of Computational Methods, 2011, 8(2): 315-347. 本文引用 [1]

[14]	李伟, 郑宏. 基于数值流形法的渗流问题边界处理新方法[J]. 岩土工程学报, 2017, 39(10): 1867-1873. (LI Wei, ZHENG Hong. New Boundary Treatment for Seepage Flow Problem Based on Numerical Manifold Method[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(10): 1867-1873.) (in Chinese) 本文引用 [1]

[15]	HU M, WANG Y, RUTQVIST J. Development of a Discontinuous Approach for Modeling Fluid Flow in Heterogeneous Media Using the Numerical Manifold Method[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2015, 39(17): 1932-1952. 本文引用 [1]

[16]	HU M, WANG Y, RUTQVIST J. An Effective Approach for Modeling Fluid Flow in Heterogeneous Media Using Numerical Manifold Method[J]. International Journal for Numerical Methods in Fluids, 2015, 77(8): 459-476. 本文引用 [2]

[17]	AALTO J. Finite Element Seepage Flow Nets[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1984, 8(3): 297-303. 本文引用 [1]

[18]	TRACY F T, RADHAKRISHNAN N. Automatic Generation of Seepage Flow Nets by Finite Element Method[J]. Journal of Computing in Civil Engineering, 1989, 3(3): 268-284. 本文引用 [2]