Voronoi杂交渗流通量有限元方法研究

罗超; 张蕊; 郭然; 盖文海

doi:10.21656/1000-0887.460148

Voronoi杂交渗流通量有限元方法研究

doi: 10.21656/1000-0887.460148

罗超^1,,
张蕊^1, ,,
郭然¹,
盖文海²

1.
昆明理工大学建筑工程学院，昆明 650500
2.
云南交通职业技术学院公路与建筑学院，昆明 650500

(本刊编委郭然来稿)

基金项目:

国家自然科学基金 12062007

云南省“兴滇英才”支持计划基金 12072135

详细信息

作者简介:
罗超(1997—)，男，硕士生(E-mail: 20232110004@stu.kust.edu.cn)

通讯作者:
张蕊(1982—)，女，教授，博士生导师(通信作者. E-mail: zrui0310@163.com)

中图分类号: O343
计量
- 文章访问数: 107
- HTML全文浏览量: 39
- PDF下载量: 23
- 被引次数: 0
出版历程
- 收稿日期: 2025-08-19
- 修回日期: 2025-09-12
- 刊出日期: 2026-06-01

Investigation of a Voronoi-Based Hybrid Seepage Flux Finite Element Method

LUO Chao^1
,,
ZHANG Rui^{1
, ,},
GUO Ran¹,
GAI Wenhai²

1.
Faculty of Civil Engineering and Mechanics, Kunming University of Science and Technology, Kunming 650500, P.R. China
2.
School of Highway and Construction Engineering, Yunnan Communications Vocational and Technical College, Kunming 650500, P.R. China

(Contributed by GUO Ran, Member of the Editorial Board of AMM)

摘要

摘要: 针对土石混合体和各向异性特征的天然土质坝基，将其视为一种复合多孔介质土体模型，假定渗流为二维Darcy流，水头在土石混合体内连续. 构建了考虑土石界面的非均质土石混合体渗流场和各向异性均质坝基渗流场的杂交渗流通量有限元(HS-FEM)求解模型，该模型在单元域内假设独立的高阶渗流通量变量，同时在单元边界上给定水头值，仅需少数单元即可模拟考虑土石界面的承压渗流工况，避免了传统有限元法在界面处进行密集网格细化的弊端，亦可解决二维正交各向异性稳态线性渗流问题. 验证算例表明：该方法在保持稀疏离散化的前提下，获得了与传统密集网格有限元相当的计算精度.
- 杂交渗流通量有限元法(HS-FEM) /
- 渗流分析 /
- 复合材料
Abstract: The soil-rock mixtures and naturally anisotropic dam foundations were treated as composite porous media, and the seepage flow was assumed to follow 2D Darcy's law with continuous hydraulic head distributions within the soil-rock mixture. A hybrid seepage flux finite element method (HS-FEM) model was developed for analyzing seepage fields in soil-rock mixtures with interfaces and anisotropic homogeneous dam foundations. For this model, independently assumed higher-order seepage flux variables was adopted within element domains with hydraulic head values prescribed on element boundaries. Only a limited number of elements were required to effectively simulate confined seepage conditions involving soil-rock interfaces, thereby to overcome the drawback of traditional FEMs necessitating dense mesh refinement at material interfaces. Additionally, the method is capable of solving 2D orthotropic steady-state linear seepage problems. Numerical examples demonstrate that the proposed HS-FEM achieves comparable accuracy to traditional dense-mesh FEMs while maintaining computational efficiency through sparse discretization.
- hybrid seepage flux finite element method (HS-FEM) /
- seepage analysis /
- composite material
edited-by

edited-by
注释:

1) (本刊编委郭然来稿)

HTML全文

图 1 单元域划分

Figure 1. Unit domain divisions

下载: 全尺寸图片幻灯片

图 2 土石混合体模型对比

注为了解释图中的颜色，读者可以参考本文的电子网页版本，后同.

Figure 2. Comparison of soil-rock mixture models

下载: 全尺寸图片幻灯片

图 3 水头及渗流通量分布云图对比

Figure 3. Comparative contour plots of hydraulic head and seepage flux

下载: 全尺寸图片幻灯片

图 4 工况条件下的各向异性均质坝基模型

Figure 4. The anisotropic homogeneous dam foundation model under operational conditions

下载: 全尺寸图片幻灯片

图 5 坝基模型对比

Figure 5. Comparison of dam foundation models

下载: 全尺寸图片幻灯片

图 6 水头等值线对比

Figure 6. Comparison of hydraulic head contours

下载: 全尺寸图片幻灯片

图 7 渗流通量云图对比

Figure 7. Comparison of seepage flux contours

下载: 全尺寸图片幻灯片

图 8 y=58 m渗流通量路径比较

Figure 8. Comparison of seepage fluxes paths at y=58 m

下载: 全尺寸图片幻灯片

图 9 基函数矩阵P项数对渗流通量q_x的影响

Figure 9. Influences of the basis function matrix P term number on seepage components q_x

下载: 全尺寸图片幻灯片

参考文献(31)

[1]	NEUMAN S P, WITHERSPOON P A. Finite element method of analyzing steady seepage with a free surface[J]. Water Resources Research, 1970, 6(3): 889-897. doi: 10.1029/WR006i003p00889
[2]	ZHANG W, DAI B, LIU Z, et al. Unconfined seepage analysis using moving Kriging mesh-free method with Monte Carlo integration[J]. Transport in Porous Media, 2017, 116(1): 163-180. doi: 10.1007/s11242-016-0769-9
[3]	DAI Q, LEI Y, ZHANG B, et al. A practical adaptive moving-mesh algorithm for solving unconfined seepage problem with Galerkin finite element method[J]. Scientific Reports, 2019, 9: 6988. doi: 10.1038/s41598-019-43391-4
[4]	SHARMA V, FUJISAWA K, MURAKAMI A. Space-time finite element method for transient and unconfined seepage flow analysis[J]. Finite Elements in Analysis and Design, 2021, 197: 103632. doi: 10.1016/j.finel.2021.103632
[5]	CHUGH A K, FALVEY H T. Seepage analysis in a zoned anisotropic medium by the boundary element method[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1984, 8(4): 399-407. doi: 10.1002/nag.1610080407
[6]	柴军瑞, 仵彦卿. 采用边界元法确定渗流自由面的一种改进方法[J]. 陕西水力发电, 2000(1): 11-12. CHAI Junrui, WU Yanqing. A modified iterative approach for solving seepage free surface by BEM[J]. Journal of Shaanxi Water Power, 2000(1): 11-12. (in Chinese)
[7]	葛锦宏, 李广信, 介玉新. 无单元法在有自由面渗流计算中的应用[J]. 计算力学学报, 2003, 20(2): 241-243. GE Jinhong, LI Guangxin, JIE Yuxin. Application of element-free method to the seepage calculation with free surface[J]. Chinese Journal of Computational Mechanics, 2003, 20(2): 241-243. (in Chinese)
[8]	张石峰, 李茜, 高佩玲. 天然差分方法在地下水渗流问题中的应用[J]. 计算物理, 2007, 24(3): 307-312. ZHANG Shifeng, LI Qian, GAO Peiling. Natural difference method for groundwater seepage[J]. Chinese Journal of Computational Physics, 2007, 24(3): 307-312. (in Chinese)
[9]	BARDET J P, TOBITA T. A practical method for solving free-surface seepage problems[J]. Computers and Geotechnics, 2002, 29(6): 451-475. doi: 10.1016/S0266-352X(02)00003-4
[10]	李凤志. 渗流自由面分析的比例边界有限元法[J]. 计算物理, 2009, 26(5): 665-670. LI Fengzhi. Scaled boundary finite-element method for seepage free surfaces analysis[J]. Chinese Journal of Computational Physics, 2009, 26(5): 665-670. (in Chinese)
[11]	李珏, 金吾根, 薛馨, 等. Trefftz间接法解自由面渗流问题[J]. 计算力学学报, 2005, 22(3): 295-298. LI Jue, JIN Wugen, XUE Xin, et al. Solution of free surface seepage problems by Trefftz indirect method[J]. Chinese Journal of Computational Mechanics, 2005, 22(3): 295-298. (in Chinese)
[12]	薛馨, 金吾根, 李珏, 等. Trefftz直接法解渗流自由面问题[J]. 岩石力学与工程学报, 2005, 24(13): 2322-2326. XUE Xin, JIN Wugen, LI Jue, et al. Solution of free surface seepage problems by Trefftz direct method[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(13): 2322-2326. (in Chinese)
[13]	徐文杰, 王永刚. 土石混合体细观结构渗流数值试验研究[J]. 岩土工程学报, 2010, 32(4): 542-550. XU Wenjie, WANG Yonggang. Meso-structural permeability of S-RM based on numerical tests[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(4): 542-550. (in Chinese)
[14]	YAN L, MENG Q X, XU W Y, et al. A numerical method for analyzing the permeability of heterogeneous geomaterials based on digital image processing[J]. Journal of Zhejiang University (Science A): Applied Physics & Eingineering, 2017, 18(2): 124-137.
[15]	WEI W, JIANG Q, YE Z, et al. Equivalent fracture network model for steady seepage problems with free surfaces[J]. Journal of Hydrology, 2021, 603: 127156. doi: 10.1016/j.jhydrol.2021.127156
[16]	陈远强, 郑宏, 屈新. 基于数值流形法的降雨入渗与坡面径流耦合算法研究[J]. 应用数学和力学, 2023, 44(12): 1499-1511. CHEN Yuanqiang, ZHENG Hong, QU Xin. A coupling analysis of rainfall infiltration and slope surface runoff based on the numerical manifold method[J]. Applied Mathematics and Mechanics, 2023, 44(12): 1499-1511. (in Chinese)
[17]	JIA Z, ZHENG H. A new procedure for locating free surfaces of complex unconfined seepage problems using fixed meshes[J]. Computers and Geotechnics, 2024, 166: 106032. doi: 10.1016/j.compgeo.2023.106032
[18]	ZHAO L, ZHANG H, MAO J, et al. Numerical simulation of the bubbling phenomenon ahead of the dam using a two-phase seepage free surface model[J]. Computers and Geotechnics, 2023, 155: 105239. doi: 10.1016/j.compgeo.2022.105239
[19]	ZHAO L, ZHANG H, MAO J, et al. An ICLS-based method for solving two-phase seepage free surface considering compressible gas in porous media[J]. Computers and Geotechnics, 2022, 141: 104528. doi: 10.1016/j.compgeo.2021.104528
[20]	李树光, 曲凯. 多孔介质中单相气体局部流动的均质化建模[J]. 应用数学和力学, 2024, 45(2): 175-183. doi: 10.21656/1000-0887.440246 LI Shuguang, QU Kai. Homogenization modeling of single-phase gas local flow in porous media[J]. Applied Mathematics and Mechanics, 2024, 45(2): 175-183. (in Chinese) doi: 10.21656/1000-0887.440246
[21]	PIAN T H H. Derivation of element stiffness matrices by assumed stress distributions[J]. AIAA Journal, 1964, 2(7): 1333-1336. doi: 10.2514/3.2546
[22]	PIN T, PIAN T H H. A variational principle and the convergence of a finite-element method based on assumed stress distribution[J]. International Journal of Solids and Structures, 1969, 5(5): 463-472. doi: 10.1016/0020-7683(69)90036-5
[23]	PIAN T H H, TONG P. Basis of finite element methods for solid continua[J]. International Journal for Numerical Methods in Engineering, 1969, 1(1): 3-28. doi: 10.1002/nme.1620010103
[24]	GHOSH S, MUKHOPADHYAY S N. A material based finite element analysis of heterogeneous media involving Dirichlet tessellations[J]. Computer Methods in Applied Mechanics and Engineering, 1993, 104(2): 211-247. doi: 10.1016/0045-7825(93)90198-7
[25]	GHOSH S, MALLETT R L. Voronoi cell finite elements[J]. Computers & Structures, 1994, 50(1): 33-46.
[26]	GUO R, SHI H J, YAO Z H. Modeling of interfacial debonding crack in particle reinforced composites using Voronoi cell finite element method[J]. Computational Mechanics, 2003, 32(1): 52-59.
[27]	ZHANG R, GUO R. Voronoi cell finite element method for fluid-filled materials[J]. Transport in Porous Media, 2017, 120(1): 23-35. doi: 10.1007/s11242-017-0898-9
[28]	ZHANG R, GUO R. Determination of crack tip stress intensity factors by singular Voronoi cell finite element model[J]. Engineering Fracture Mechanics, 2018, 197: 206-216. doi: 10.1016/j.engfracmech.2018.04.042
[29]	KARCHEWSKI B. Multi-scale modelling of geomechanical behaviour using the Voronoi cell finite element method (VCFEM) and finite-discrete element method (VCFEM-DEM)[D]. Hamilton, Ontario: McMaster University, 2015.
[30]	FELIPPA C A. Parametrized multifield variational principles in elasticity, Ⅰ: mixed functionals[J]. Communications in Applied Numerical Methods, 1989, 5(2): 79-88. doi: 10.1002/cnm.1630050204
[31]	FELIPPA C A. A survey of parametrized variational principles and applications to computational mechanics[J]. Computer Methods in Applied Mechanics and Engineering, 1994, 113(1/2): 109-139.