等几何分析的多重网格共轭梯度法.pdf

　 ， ３５ ６ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１４年６月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３５，Ｎｏ．６，Ｊｕｎ．１５，２０１４ 文章编号：１０００⁃０８８７（２０１４）０６⁃０６３０⁃１０ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 等几何分析的多重网格共轭梯度法∗ １ ２ １ ２ １ 刘　 石 陈德祥 冯永新 徐自力 郑李坤 ，　 ，　 ，　 ，　 （１． 广东电网公司 电力科学研究院，广州 ５１００８０； ２． 西安交通大学 航天航空学院；机械结构强度与振动国家重点实验室，西安 ７１００４９） 摘要：　 提高ＮＵＲＢＳ基函数阶数可以提高等几何分析的精度，同时也会降低多重网格迭代收敛速 度．将共轭梯度法与多重网格方法相结合，提出了一种提高收敛速度的方法，该方法用共轭梯度法 作为基础迭代算法，用多重网格进行预处理．对Ｐｏｉｓｓｏｎ（泊松）方程分别用多重网格方法和多重网 格共轭梯度法进行了求解，计算结果表明：等几何分析中采用高阶ＮＵＲＢＳ基函数处理三维问题 时，多重网格共轭梯度法比多重网格法的收敛速度更快． 关　 键　 词：　 等几何分析；　 多重网格；　 共轭梯度法；　 Ｐｏｉｓｓｏｎ方程；　 迭代算法；　 ＮＵＲＢＳ 中图分类号：　 Ｏ２４１．８２　 　 　 文献标志码：　 Ａ ｄｏｉ：１０．３８７９／ ｊ．ｉｓｓｎ．１０００⁃０８８７．２０１４．０６．００５ 引　 　 言 ［１］ 等几何分析是求解偏微分方程的一种新数值方法 ，它与有限元法一样都基于Ｇａｌｅｒｋｉｎ 变分，但是采用非均匀有理Ｂ样条（ＮＵＲＢＳ）作为基函数．因为ＮＵＲＢＳ基函数在单元边界上能 ｋ ［２⁃３］ 够实现Ｃ （ｋ ≥１）连续性，提高光滑性ｋ可以提高数值解的精度 ，所以等几何分析在精度 ［４］ 方面优于传统有限元方法 ． 当等几何分析产生的代数方程规模很大时，直接求解计算成本过高，一般采用迭代方法求 ［５⁃６］ ［７⁃８］ 解．多重网格迭代是一种高效的迭代算法，该方法在有限容积法 、有限差分法 和有限元 ［９］ ［１０］ 方法 中的研究较多，但在等几何分析中的研究较少，仅Ｇａｈａｌａｕｔ等 在２０１３年研究了用多 重网格算法求解等几何分析中的代数方程．文献［１０］的结果表明多重网格方法具有很快的收 敛速度，但是当基函数的次数或光滑性增加时，多重网格的收敛速度降低．因此，在提高基函数 光滑性获得更高精度的同时，多重网格的效率也随之降低． 多重网格的基本思想是在密网格上进行光顺，使高频误差衰减，在粗网格上进行误差修 ［７］ 正，使低频误差衰减 ．常用的光顺算法是Ｊａｃｏｂｉ迭代或Ｇａｕｓｓ⁃Ｓｅｉｄｅｌ迭代，随着基函数光滑性 升高，Ｊａｃｏｂｉ迭代或Ｇａｕｓｓ⁃Ｓｅｉｄｅｌ迭代的光顺效果降低，多重网格的收敛速度降低．共轭梯度法 是一种Ｋｒｙｌｏｖ子空间迭代方法，它适用于系数矩阵为对称正定的情况．理论上若代数方程的维 数为ｎ，那么进行ｎ次共轭梯度法迭代将收敛到精确解；实际计算中共轭梯度法的收敛速度快 ［１１］ 慢