数值分析

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

展开
  • 1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071;2. 湖北省水利水电科学研究院,湖北 武汉 430070
刘治军,男,1988年生,博士,主要从事计算力学和岩土力学方面的研究工作

收稿日期: 2015-06-02

  网络出版日期: 2018-06-05

基金资助

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

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

Expand
  • 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 date: 2015-06-02

  Online published: 2018-06-05

Supported by

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

摘要

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

本文引用格式

刘治军,郑 宏,董 苇,葛修润,孙冠华, . 基于加密物理片的数值流形法中局部网格加密[J]. 岩土力学, 2017 , 38(4) : 1211 -1217 . DOI: 10.16285/j.rsm.2017.04.036

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.
文章导航

/