预条件GMRESm法迭代求解大规模边界元弹性问题研究.pdfVIP

第27卷专辑 固体力学学报 H 2006年12月 ACTAMECANlCASOLIDASINlCA 预条件GMRES(m)法迭代求解 大规模边界元弹性问题。 陈泽军 肖 宏 (燕山大学机械工程学院轧机研究所．秦皇岛，066004) 摘要 提出了一种基于单元节点的块雅可比预条件方法，扩大了边界元法的计算规模．使之可用于大规模 工程同题的求解．数值实验说明了这种顶条件技术的有效性．表明预条件GMRESCm)算法具有较好的收敛特性。 适合于求解大规模问题边界元弹性问题所形成的稠密非对称线性方程组． 关t词 广义极小残值算法(GMRES)，预条件．边界元法．弹性问题．迭代法 0 引言 之真正能够进行大规模工程问题的计算． 边界元法是将描述问题的偏微分方程转化为边 1 弹性问题边界元法 界积分方程，并吸收了有限元法的离散化技术而发 在不计体力的情况下，三维弹性问题的位移边 展起来的．边界元法具有高精度和降维等优点。已经 界积分方程可以写作 成为一种精确、高效的工程数值分析方法，被大量地 勺(2)坳(z)+IT。(z，y)uj(y)dr(y)= 用于科学与工程问题的计算分析，被公认为有限元 法一个最重要的补充．边界元法的主要运算量包括 I砜(z，y)白(y)dr(y)(1) 两个部分，影响系数矩阵的边界积分计算和形成的 线性方程组的求解．由于边界元法形成的线性方程 这里q(z)=岛+口。(z)．‘。o)是边界点工几何形 组的系数矩阵一般是非对称的满秩矩阵，若采用直 状的函数．如果工是光滑边界上的点，那么边界的法 接法求解(例如高斯消去法)，所需的运算量为 矢量在z点是连续的，则它可以写成c。(z)=如／Z． 0(Ⅳ)，这里N为问题的自由度数．因此，边界元法 难以求解大规模问题，这限制了边界元法的发展．为 它们被称为弹性问题的基本解． 了减少求解方程组运算量【1]，更多的是采用迭代法 三维弹性问题的基本解砜(z，y)和L(z，，) 进行求解．使用迭代法求解，其计算量约为 可以表示为[6】 O(rdV2)，咒为迭代次数．可见，当计算大规模问题 时，迭代法的优势将更加明显． 近几年，国内外许多学者在运用迭代法求解边 (2) 界元形成线性方程组方面进行了许多研究[2~5]，广 MinimalRESidual) 义极小残值算法(Generalized 是一种适合于求解稠密非对称线性方程组的有效算 法．随着问题规模的扩大，边界元法形成的线性方程 组逐渐呈现病态，导致迭代的收敛速度变慢，甚至发 (3) 生中断．这样对于某一个问题，达到一定的计算精度 其中，R是源点z和积分点Y之间的距离． 时所需的迭代次数花也随之增大，使得迭代法的优 边界离散后，(1)式可以表示为 势有所减弱．因此，在运用迭代法求解的过程中必须 白(2)吻(z)+ 对线性方程组进行有效的预条件处理．这样才能真 正发挥迭代法的优势，扩大边界元法的计算规模。使 ■一l、·一1“ J -国家自然科学基金资助． 专 辑 陈泽军等： 预条件GMRES(m)法迭代求解大规模边界元弹性问题 ·51·