在线咨询
中国工业与应用数学学会会刊
主管:中华人民共和国教育部
主办:西安交通大学
ISSN 1005-3085  CN 61-1269/O1

工程数学学报 ›› 2018, Vol. 35 ›› Issue (5): 534-544.doi: 10.3969/j.issn.1005-3085.2018.05.005

• • 上一篇    下一篇

组合杂交四边形元的多重网格预处理共轭梯度方法

王惠玲,   聂玉峰,   张   玲   

  1. 西北工业大学计算科学研究中心,西安  710129
  • 收稿日期:2016-06-29 接受日期:2017-12-08 出版日期:2018-10-15 发布日期:2018-12-15
  • 通讯作者: 聂玉峰 E-mail: yfnie@nwpu.edu.cn
  • 基金资助:
    国家自然科学基金(11471262; 11501450).

The Multigrid Preconditioned Conjugate Gradient Method for Combined Hybrid Quadrilateral Element

WANG Hui-ling,   NIE Yu-feng,   ZHANG Ling   

  1. Research Center for Computational Science, Northwestern Polytechnical University, Xi'an 710129
  • Received:2016-06-29 Accepted:2017-12-08 Online:2018-10-15 Published:2018-12-15
  • Contact: Y. Nie. E-mail address: yfnie@nwpu.edu.cn
  • Supported by:
    The National Natural Science Foundation of China (11471262; 11501450).

摘要: 组合杂交元方法是一种求解弹性力学问题的稳定化有限元方法.为了快速求解组合杂交元离散得到的大型、稀疏、对称正定系统,本文研究了多重网格预处理共轭梯度方法.首先,通过选用合适的网格转移算子和光滑策略,得到了有效的多重网格预处理器.其次,通过分析数值试验结果证明所得到的多重网格预处理共轭梯度方法是有效可行的,利用该预处理方法大大降低了系数矩阵的条件数,提高了计算效率.此外,对于一类高性能的组合杂交元,多重网格预处理共轭梯度方法在网格畸变时依然收敛.

关键词: 组合杂交元, Wilson插值, 网格转移算子, 多重网格方法, 预处理共轭梯度方法

Abstract: Combined hybrid finite elements method, applied to linear elasticity problem, is a stabilized finite element method. It is a large sparse symmetric positive definite system arising from combined hybrid quadrilateral element discretization. In order to solve the system quickly, the multigrid preconditioned conjugate gradient method (MGCGM) is introduced in this paper. Firstly, by choosing the appropriate intergrid transfer operators and smoothing strategy, an effective multigrid preconditioner is obtained. Then, numerical results show that MGCGM proposed is efficient, that is the condition number of the stiffness matrix is greatly decreased by the multigrid preconditioner. Moreover, the method still has a good convergence for the combined hybrid element with high performance in the case of the mesh distortion.

Key words: combined hybrid element, Wilson interpolation, intergrid transfer operator, multigrid method, preconditioned conjugate gradient method

中图分类号: