09第5章并行多极展开GMRES的研究及其在边界元中应用.docVIP

第5章 多极展开GMRES算多极展开GMRES算法在边界元中的应用5.1.1 弹性边界元法简介 弹性问题边界元法，由Rizzo. F. J[74]在1967年给出，将积分方程——Somigliana公式[75]适用于边界点上，采用常单元离散边界而形成矩阵方程。弹性体基本解函数——Kelvin解，用Galerkin矢量求出。 借助加权残量法的思想，利用格林公式，静力学基本方程可以化为以下的边界积分方程 (5-1) 式中，和分别为求解问题的位移和面力，和为三维弹性问题的基本解，由所在的位置决定。 为了对边界积分方程求数值解，将整个边界离散划分为n个单元，采用一定的插值函数将每个边界元内的位移矢量和面力矢量表示为离散形式 (5-2) (5-3) 代入原积分方程后将其转化为代数方程 再将该代数方程中的未知边界量及其相应系数归到方程左边，已知边界量及相应系数移至右边，可以得到标准形式的线性代数方程组 求解方程组得到未知边界量进而求出边界及域内应力应变值。 5.1.2 多极展开GMRES在边界元中应用的基本思想 多极展开GMRESGMRES算法应用于边界元法中可大大简化计算，设边界上有N个节点M个单元，传统的边界元法要完成计算，需每个节点对所有单元积分，所需计算量与M×N成正比，如图5-1所示；若采用多极展开法，将边界划分为K个区域，不妨假设每个区域内的单元数同为M/K，则每个区域内所有单元向区域中心的规约计算量与M/K成正比，故K个区域的总计算量与M成正比，考虑到相邻区域需直接计算，总的计算量和(M+N(K-3+3M/K))成正比，如图5-2所示。 图5-1 单元间直接计算 Fig. 5-1 Direct Computing of Elements 图5-2 采用多极展开法计算 Fig. 5-2 Computing with FMM 5.1.3 边界元基本解核函数的等价变形 考虑有限体弹性问题的边界积分方程，并不影响扩展到无限体、含有杂质、裂纹等问题上。认为有限体域为，其表面边界为，面力已知表面为，位移已知表面为，，若不考虑体积力的影响，可得边界位移和面力的边界积分方程为 (5-4) 式中，为边界上的任意一点，为源点，为边界形状系数。和是三维弹性问题的位移和面力基本解核函数[75]，通常表示为 (5-5) (5-6) 式中，为Kronecker 函数，是二阶不变张量，为观测点和源点间的距离，为边界的外法矢。 边界积分方程(5-4)通过离散变成代数方程。边界上单元坐标表示为 (5-7) 式中，是单元的局部坐标，是单元的节点的坐标，是节点的形函数。 在点对(5-4)式数值积分，为单元的积分点，为 处的积分权函数。(5-4)式改写为 (5-8) 在(5-8)式中,积分在主单元区域进行，是与(5-7)式对应的雅可比行列式。在点用(5-8)式数值积分得 (5-9) 通过(5-9)可以看出，对每一个单元区域上的运算， 在所有边界点对其积分时是不变的。如果采用迭代的方法求解方程，和在迭代前赋值，那么在每一迭代步中，乘积项对各个单元也是定值，同理和相关的项具有同样的性质。对于(5-9)中的两个结构相同的求和项，先讨论其中第一个，为简化推理过程可令 (5-10) 对于基本解项，为适用于多极展开近似求和公式，需首先将其等价变形成项的函数，即如下形式 (5-11) 式中，，将其代入(5-9)式中的第一求和项，再提出不变项，可变形为 (5-12) 同理，(5-9)可全部改写为 (5-13) 显然，(5-13)式中的求和项可代入多极展开公式快速计算，而无需形成系数矩阵。因为，，所以第一项求和需要调用多极展开9次，第二项为3次，第三项为3次，第四项为1次，总共为16次。 5.1.4 多极展开GMRES多极展开GMRES多极展开GMRES时，矩阵乘法和使用多极展开法计算，使用多极展开法时，首先计算每个单元的多极系数，每个单元有16个系数。当和变化时单元的多极系数需要重新计算。如果近似解满足精度要求，迭代终止，否则以得到的近似解为初解继续迭代。计算流程如图5-3所示。 图5-3 边界元中多展开GMRES法的计算流程 Fig. 5-3 The Flow Chart of GMRES based on FMM in BEM 边界元中