›› 2017, Vol. 38 ›› Issue (4): 1211-1217.doi: 10.16285/j.rsm.2017.04.036

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

基于加密物理片的数值流形法中局部网格加密

刘治军1,郑 宏1,董 苇2,葛修润1,孙冠华1   

  1. 1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071;2. 湖北省水利水电科学研究院,湖北 武汉 430070
  • 收稿日期:2015-06-02 出版日期:2017-04-11 发布日期:2018-06-05
  • 作者简介:刘治军,男,1988年生,博士,主要从事计算力学和岩土力学方面的研究工作
  • 基金资助:

    国家自然科学基金项目(No. 11572009,No. 51538001)

Local mesh refinement in numerical manifold method based on refined physical patches

LIU Zhi-jun1, ZHENG Hong1, DONG Wei2, GE Xiu-run1, SUN Guan-hua1   

  1. 1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China; 2. Hubei Water Resources Research Institute, Wuhan, Hubei 430070, China
  • Received:2015-06-02 Online:2017-04-11 Published:2018-06-05
  • Supported by:

    This work was supported by the National Natural Science Foundation of China(11572009, 51538001).

摘要: 数值流形法的数学网格不需要适应求解域内各种不连续界面和边界,因此,总可以用规则的结构化网格建立数学覆盖。但对于大多数问题,在整个求解域上布置统一密度的网格显得浪费。因此,需要研究在结构化网格上实施局部加密,建立了加密物理片的方法用以解决这一问题。具体的实施过程中,确定了需要加密的区域后,先在这些区域布置规则的精细网格,然后找到这些区域中包含的原始网格中的物理片,用精细网格上建立的插值代替被加密物理片上的局部近似,从而提高了局部近似的阶次。数值算例结果表明,该方法收敛性良好。另外,如果所有物理片上的局部近似都采用0阶多项式(常数),那么将会得到正定的刚度矩阵。

关键词: 数值流形法, 加密物理片, 结构化网格, 局部加密

Abstract: The mathematical mesh of numerical manifold method (NMM) does not have to accommodate to various boundaries of physical domains, and thus the mathematical coverage is always built by the regular structured mesh. However, for most problems, it is wasteful to use the mesh with the uniform density on the the whole physical region. Therefore, it is necessary to study the implementation of the local refinement on structured mesh, and a method of refining physical patches is proposed to solve the problem. For a practical problem, firstly we determine each region in which the mesh needs to be refined, and it is found that the physical patches entirely is contained by the refined mesh. Then, an interpolation on the refined mesh is constructed inside each physical patch, and the original local approximation of the physical patch is replaced by the new interpolation. Thus, the order of the local approximation is improved. Numerical results show the proposed method has good convergence. In addition, for two-dimensional analysis, the stiffness matrix obtained by the proposed method is positive definite if local approximations on all physical patches are constant.

Key words: numerical manifold method, refined physical patches, structured mesh, local refinement

中图分类号: 

  • O 242

[1] 聂治豹, 郑宏, 万涛, 林姗. 弹性静力学问题的边界积分形式的数值流形法[J]. 岩土力学, 2020, 41(4): 1429-1436.
[2] 杨石扣, 张继勋, 任旭华, . 基于改进数值流形法的接触裂纹问题研究[J]. 岩土力学, 2019, 40(5): 2016-2021.
[3] 刘登学, 张友良, 丁秀丽, 黄书岭, 裴启涛, . 数值流形法中基于适合分析T样条的 局部网格加密算法[J]. 岩土力学, 2019, 40(4): 1584-1595.
[4] 杨石扣,张继勋,任旭华,. 基于数值流形法的三维裂纹扩展研究[J]. , 2018, 39(S1): 488-494.
[5] 杨石扣,任旭华,张继勋,. 基于数值流形法的重力坝水力劈裂研究[J]. , 2018, 39(8): 3055-3060.
[6] 杨石扣,任旭华,张继勋,. 改进的数值流形法在水力劈裂中的应用[J]. , 2018, 39(10): 3875-3881.
[7] 陈 立,张 朋,郑 宏,. 土石混合体二维细观结构模型的建立与数值流形法模拟[J]. , 2017, 38(8): 2402-2410.
[8] 张友良,刘登学,刘高敏. 数值流形法的T样条局部加密算法[J]. , 2016, 37(8): 2404-2410.
[9] 徐栋栋,杨永涛,郑 宏,邬爱清,. 基于数值流形法的重力坝多裂纹扩展研究[J]. , 2016, 37(12): 3598-3607.
[10] 徐栋栋 ,杨永涛 ,郑 宏 ,邬爱清,. 自由度具有物理含义的线性无关高阶数值流形法[J]. , 2016, 37(10): 2984-2992.
[11] 杨石扣,张继勋,任旭华,张道法. 三维数值流形法在裂纹扩展中的应用研究[J]. , 2016, 37(10): 3017-3025.
[12] 徐栋栋,郑 宏. 数值流形法在处理强奇异性问题时的网格无关性[J]. , 2014, 35(8): 2385-2393.
[13] 刘 建,陈 佺. 岩体黏弹性蠕变计算的高阶数值流形法研究[J]. , 2012, 33(7): 2174-2180.
[14] 林兴超 汪小刚 王玉杰 李 旭 韩鑫晔. 数值流形法中“质量守恒”的探讨[J]. , 2011, 32(10): 3065-3070.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!