渗流转角奇异问题的数值流形单元法

doi:10.13722/j.cnki.jrme.2017.0057

摘要
图/表
参考文献(36)
相关文章 (15)

全文: PDF (386 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要

在水利工程中的渗流计算中，大量存在着转角奇异性问题。一般有限元方法需要构建特殊的奇异形函数和单元属性来处理，存在不便利性，而且广泛使用的等参单元又忽视了奇异性问题。运用数值流形方法，提出了采用构造局部特殊片的方法来解决该奇异性问题。首先，将转角奇异的解析解构建为一个基函数，与常数项基函数一起形成奇异片覆盖函数，在奇异点影响到的星点上采用奇异片覆盖函数，其他星点采用常数片覆盖函数进行求解，计算结果逼近真解，并能很好地反映转角中的强奇异性。其次，给出这类问题的流形元离散列式，将经典正三角形流形元改进为正四边形，并介绍程序的编制步骤，和该特殊物理片的施加方法。再次，3个工程算例分别为坝基不透水帷幕有压渗流解析解对比、经典流速跌落试验对比以及堰下板桩有压渗流与已有计算流网对比，算例表明了该方法的有效性和正确性，为进一步分析奇异性对结构和复杂工程的影响打下基础。该方法也适用于其他基于单位分解插值的数值分析方法，也为其他场变量连续而偏导数不连续问题提供了思路。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	冯强1
	2
	郑宏1
	2
	3

关键词 ：数值分析, 数值流形法, 渗流, 转角奇异性

Abstract：

In seepage calculation of hydraulic engineering，there are a large number of singularity problems at turning point. It is inconvenience that the general finite element method needs to construct special singularity shape function and element property to deal with the singularity problem，which used to be ignored when the isoparametric element was widely used. In this paper，a new method by constructing a local special patch is proposed to solve this singularity problem by using numerical manifold method. Firstly，the singularity covering function is formed together by the constant term basis function and a new basis function constructed by the analytic solution at turning point. The singularity covering function is used on the star points affected by the singularity point，at the same time，the constant patch covering function is used on the other star points，and the calculation result approximates the real solution and can well reflect the strong singularity in the corner. Secondly，the discrete formula of the numerical manifold method about this kind of problem is given，and the classical regular triangle of manifold element is improved to square，as well as，the procedure of the program and the application method of the special physical patch are introduced. Thirdly，three engineering examples are

respectively confined seepage of dam foundation with impervious curtain contrast with the analystic solution，the comparison of the classical flow velocity drop test，and the confined seepage of weir with sheet pile compared with flow net which had been calculated，those show the effectiveness and correctness of the method，and provide a basis for further analysis of the effects of singularity on the structure and complex engineering. Finally，The proposed method is also applicable to other numerical analysis methods based on the partition of unity method of interpolation，and also provides an approach for those problems with continuous field and discontinuous partial derivatives.

Key words： numerical analysis numerical manifold method seepage turning point singularity

引用本文:

冯强1，2，郑宏1，2，3. 渗流转角奇异问题的数值流形单元法[J]. 岩石力学与工程学报, 2017, 36(S2): 3892-3901.
FENG Qiang1，2，ZHENG Hong1，2，3. Turning point singularity problem in seepage by numerical manifold method. , 2017, 36(S2): 3892-3901.

链接本文:

https://rockmech.whrsm.ac.cn/CN/10.13722/j.cnki.jrme.2017.0057 或 https://rockmech.whrsm.ac.cn/CN/Y2017/V36/IS2/3892

[1] MELENK J M. Finite element methods with harmonic shape functions for solving Laplace's equation[M. S. Thesis][D]. Maryland：The University of Maryland，1992：1–89.

[2] MELENK J M. On generalized finite element methods[Ph. D. Thesis][D]. Maryland：The University of Maryland，1995：1–187.

[3] BABUŠKA I，OSBORN J E. Special finite element methods for a class of second order elliptic problems with rough coefficients[J]. Siam Journal on Numerical Analysis，1994，31(4)：945–981.

[4] BABUŠKA I，MELENK J. The partition of unity method[J]. International Journal for Numerical Methods in Engineering，1997，40(4)：727–758.

[5] DUARTE C A，ODEN J T. Hp clouds-an hp meshless method[J]. Numerical Methods for Partial Differential Equations，1996，12(6)： 673–706.

[6] SUKUMAR N，BELYTSCHKO T. Arbitrary branched and intersecting cracks with the extended finite element method[J]. International Journal for Numerical Methods in Engineering，2000，48：1 741–1 760.

[7] BELYTSCHKO T，LU Y Y，GU L. Element‐free Galerkin methods[J]. International Journal for Numerical Methods in Engineering，1994，37(2)：229–256.

[8] SHI G H. Manifold method of material analysis[C]// Transactions of the 9th Army conference on Applied Mathematics and Computing. Minneapolis：U. S. Army Research Office，1991：57–76.

[9] 李晓春，陈剑平，石丙飞，等. 交叉裂隙渗流无网格法初探[J]. 岩土力学，2007，28(增1)：371–374.(LI Xiaochun，CHEN Jianping，SHI Binfei，et al. A study on the meshless method on seepage of intersected fractures[J]. Rock and Soil Mechanics，2007，28(Supp.1)：371–374.(in Chinese))

[10] JIE Y X，LIU Y. Simulated annealing based algorithm for node generation in seepage analysis with meshless method[J]. Mechanics Research Communications，2012，43：96–100.

[11] ZHENG H，LIU F，LI C. Primal mixed solution to unconfined seepage flow in porous media with numerical manifold method[J]. Applied Mathematical Modelling，2015，39(2)：794–808.

[12] HUO Z Y，QIAN G X，ZHENG D J. Seepage analysis of the structure with cracks based on XFEM[J]. Open Civil Engineering Journal，2015，9(1)：90–97.

[13] 冯强. 数值流形法及其在无压渗流中的应用[硕士学位论文][D]. 宜昌：三峡大学，2005：1–74.(FENG Qiang. Numerical manifold method and its application in unconfined seepage[M. S. Thesis][D]. Yichang：China Three Gorges University，2005：1–74.(in Chinese)) .

[14] 姜清辉，邓书申，周创兵. 有自由面渗流分析的三维数值流形方法[J]. 岩土力学，2011，32(3)：879–884.(JIANG Qinghui，DENG Shushen，ZHOU Chuangbing. Three-dimensional numerical manifold method for seepage problems with free surfaces[J]. Rock and Soil Mechanics，2011，32(3)：879–884.(in Chinese))

[15] WANG Y，HU M，ZHOU Q，et al. Energy-work-based numerical manifold seepage analysis with an efficient scheme to locate the phreatic surface[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2014，38(15)：1 633–1 650.

[16] WANG Y，HU M，RUTQVIST J. Confined-unconfined seepage analysis using numerical manifold method with an energy-work-based model[C]// Frontiers of Discontinuous Numerical Methods and Practical Simulations in Engineering and Disaster Prevention，New York：CRC Press，2013：349–355.

[17] WANG Y，GONG J K. Simulation of seepage in porous medium by Numerical Manifold Method[C]// 10th International Conference on Analysis of Discontinuous Deformation- Advances in Discontinuous Numerical Methods and Applications in Geomechanics and Geoengineering. [S. l.]：[s. n.]，2012：275–280.

[18] 蔡永昌，廖林灿，张湘伟. 高精度四节点四边形流形单元[J]. 应用力学学报，2001，18(2)：75–80.(CAI Yongchang，Liao Lincan，ZHANG Xiangwei. A high-order precision quadrilateral manifold method[J]. Chinese Journal of Applied Mechanics，2001，18(2)：75–80.(in Chinese))

[19] 蔡永昌，张湘伟. 使用高阶覆盖位移函数的数值流形方法及其应力精度的改善[J]. 机械工程学报，2000，36(9)：20–25.(CAI Yongchang，ZHANG Xiangwei. Expansion to high-order cover function and improvement of the stress accuracy in numerical manifold method[J]. Chinese Journal of Mechanical Engineering，2000，36(9)：20–25.(in Chinese))

[20] 苏海东，谢小玲，陈琴. 高阶数值流形方法在结构静力分析中的应用研究[J]. 长江科学院院报，2005，22(5)：74–77.(SU Haidong，XIE Xiaoling，CHEN Qin. Application of high-order numerical manifold method in static analysis[J]. Journal of Yangtze River Scientific Research Institute，2005，22(5)：74–77.(in Chinese))

[21] 林绍忠，祁勇峰，苏海东. 基于矩阵特殊运算的高阶流形单元分析[J]. 长江科学院院报，2006，23(3)：36–39.(LIN Shaozhong，QI Yongfeng，SU Haidong. Element analysis of high-order numerical manifold method based on special matrix operations[J]. Journal of Yangtze River Scientific Research Institute，2006，23(3)：36–39.(in Chinese))

[22] 田荣，栾茂田，杨庆，等. 高阶流形方法及其应用[J]. 工程力学，2001，18(2)：21–26.(TIAN Rong，LUAN Maotian，YANG Qing，et al. Fundamentals and applications of high-order manifold method[J]. Engineering Mechanics，2001，18(2)：21–26.(in Chinses))

[23] 张国新，彭静. 二阶流形元法与结构变形分析[J]. 力学学报，2002，34(2)：261–269.(ZHANG Guoxin，PENG Jing. Second- order manifold method in structure failure analysis[J]. Acta Mechanica Sinica，2002，34(2)：261–269.(in Chinese))

[24] 邓安福，朱爱军，曾祥勇. 高低阶覆盖函数混合的数值流形方法[J]. 土木工程学报，2006，39(1)：75–78.(DENG Anfu，ZHU Aijun，ZENG Xiangyong. The unmerical manifold method using low-order and high-order cover functions[J]. China Civil Engineering Journal，2006，39(1)：75–78.(in Chinese))

[25] 彭自强，葛修润. 数值流形方法中覆盖函数选用的建议[J]. 岩土力学，2004，25(4)：624–627.(PENG Ziqiang，GE Xiurun. Suggestion to the choice of physical cover functions in numerical manifold method[J]. Rock and Siol Mechanics，2004，25(4)：624–627.(in Chinese))

[26] 林绍忠，祁勇峰，苏海东. 数值流形法中覆盖函数的改进形式及其应用[J]. 长江科学院院报，2006，23(6)：55–58.(LIN Shaozhong，QI Yongfeng，SU Haidong. Improved local cover function of numerical manifold method and its application[J]. Journal of Yangtze River Scientific Research Institute，2006，23(6)：55–58.(in Chinese))

[27] 周小义，邓安福，黄润英. 数值流形方法中覆盖位移函数的改进[J]. 岩土工程技术，2007，21(1)：1–3.(ZHOU Xiaoyi，DENG Anfu，HUANG Runying. Improvement of cover displacement functions in the numerical manifold method[J]. Geotechnical Engineering Technique，2007，21(1)：1–3.(in Chinese))

[28] 郭朝旭，郑宏. 高阶数值流形方法中的线性相关问题研究[J]. 工程力学，2012，29(12)：228–232.(GUO Chaoxu，ZHENG Hong. Study on linear dependence problem in high-order numerical manifold method[J]. Engineering Mechanics，2012，29(12)：228–232.(in Chinese))

[29] 苏海东. 固定网格的数值流形方法研究[J]. 力学学报，2011，43(1)：169–178.(SU Haidong. Study on numerical manifold method with fixed meshes[J]. Chinese Journal of Theoretical and Applied Mechanics，2011，43(1)：169–178.(in Chinese))

[30] 蔡永昌，朱合华，夏才初. 流形方法覆盖系统自动生成算法[J]. 同济大学学报：自然科学版，2004，32(5)：585–590.(CAI Yongchang，ZHU Hehua，XIA Caichu. Automatically forming of cover system in numerical manifold method[J]. Journal of Tongji University：Natural Science，2004，32(5)：585–590.(in Chinese))

[31] 祁勇峰，苏海东，崔建华. 部分重叠覆盖的数值流形方法初步研究[J]. 长江科学院院报，2013，30(1)：65–70.(QI Yongfeng，SU Haidong，CUI Jianhua. Preliminary study on numerical manifold method with partially overlapping covers[J]. Journal of Yangtze River Scientific Research Institut，2013，30(1)：65–70.(in Chinese))

[32] AALTO J. Singularity elements for seepage analysis[J]. International Journal for Numerical and Analytical Methods in Geomechanics，1985，9(2)：185–196.

[33] HARR M E. Groundwater and seepage[M]. [S. l.]：McGraw-Hill，1962：101–136.

[34] 苑莲菊. 工程渗流力学及应用[M]. 北京：中国建材工业出版社，2001：33–45.(YUAN Lianju. Engineering seepage mechanics and its application[M]. Beijing：China Building Materials Industry Press，2001：33–45.(in Chinese))

[35] 石根华. 数值流形方法与非连续变形分析[M]. 裴觉民，译. 北京：淸华大学出版社，1997：11–12.(SHI Genhua. Numerical manifold method and the discontinuous deformation analysis[M]. Translated by PEI Juemin. Beijing：Tsinghua University Press，1997：11–12(in Chinese))

[36] 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.