一种求解三维弹性问题有限元方程的并行DDM预条件子.doc

II Abstract Elasticity problem has wide applications in many ?elds, and Finite element method is one of the most commonly used methods for solving it. However, the condition number of sti?ness matrix strongly depends on mesh size, which leads to great necessity to study the corresponding fast algorithms. In this paper, we ?rst propose a new preconditioner for the linear ?nite element discrete system of 3D elastic problem by combining DDM and AMG. It is constructed based on simple coarse spaces and only needs to solve three classes of subsystems. Then corresponding fast algorithms are discussed respectively according to their speci?c characteristics. In particular, a so-called AMG-T solver is designed for the third class of subsystems by modifying the coarsening and interpolation strategies in the classical AMG(C-AMG). When the size of the subsystem is su?ciently large, experimental results show that AMG-T is better than C-AMG and CG both in the number of iteration and computing time. We also accomplish the parallel implement of the new preconditioner by using proper data structure. Numerical results show that PCG based on our new preconditioner is robust and e?cient, and that it also has good scalability. Keywords: three dimensional elasticity; preconditioner; non-overlapping DDM; AMG; parallel computing III IV I Abstract III 1 3 §2.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 §2.1.1 . . . . . . . . . . . 3 §2.1.2 Green PCG . . . . . . . . . . . . . . . . . . . 5 §2.1.3 AMG . . . . . . . . . . . . . . . . . . . 6 §2.2 . . . . . . . . . . . . . . . . . . . 7 DDM 11 §3.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 §3.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 §3.3 AMG . . . . . . . . . . . . . . . . . . . . . . . . . 13 15 §4.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 §4.2 B amg . . . . . . . . . . . . . . . . . . . . . 18 ddm §4.3 . . . . . . . . . . . . . . . . . . . . . 20 §4.4 .