GeoFEA的直接和迭代线性求解器-CAE软件库.PDFVIP

中国土木工程学会第十届土力学及岩土工程学术会议论文集 GeoFEA 的直接和迭代线性求解器 1 2 2 2 陈曦 ，方思瀚 ，王效宁 ，罗章权 （1. 北京伟思德克科技有限责任公司，上地信息路甲 28 号科实大厦 B 座 10 层 D， 北京市海淀区，100085； 2. 吉奥软件公司，2 Kaki Bukit Place, #07-00 Tritech Building, 新加坡, 416180 ） 摘 要：有限元法离散所产生的大规模线性方程组的求解是有限元模拟过程中最耗计算时间和内存的部分，因此，需要使用 高效的线性求解器来快速精确地进行有限元模拟。GeoFEA 的线性求解器分为直接和Krylov 子空间迭代求解器两大类。直接 求解器的特点是可以用固定计算量精确地得到线性方程系统的解，但对于涉及十万甚至上百万个自由度的大规模岩土工程问 题，现有的直接求解器所需要的内存和计算时间是难以接受的；迭代求解器的特点是内存需求量小，当结合高效的预处理子 时，对三维问题所产生大规模线性方程系统的计算时间要显著少于直接求解器。 关 键 词：GeoFEA ；线性方程系统；直接求解器；迭代求解器 1 引言 并不能确保都收敛； （2 ）直接求解法能够用固定的浮点运算数得到 岩土工程问题的大规模有限元分析中，空间离 精确的、可靠的解。 散和时间插分所产生的线性方程系统经常含有上 基于上面考虑， GeoFEA 既提供了波前 万、数十万、甚至成百万个未知数。求解这样大规 （Frontal ）直接求解器，也提供了多种可供选择的 模的线性方程系统，常规的高斯消去基础上的直接 对称和非对称 Krylov 子空间迭代求解器。 解法需要消耗大量计算时间和内存，以至于基于个 人计算机的模拟是十分困难的。自从 Hestenes 和 [1] 2 GeoFEA 的波前直接求解器 Stiefel 发现共轭梯度法（CG–Conjugate Gradient ） 以来，求解线性方程的 Krylov 子空间迭代法有了长 一个线性方程系统可用下面的通用形式表示： N×N N 足的发展，尤其是近二十年的发展尤为显著。基于 Ax = b, 这里A= [aij] ∈ℜ , 并且x , b ∈ℜ (1) 个人计算机的 Krylov 子空间迭代计算法已经成功 求解线性方程系统 （1）的各种直接求解法都是传统 地应用于大规模二维，尤其是三维岩土工程有限元 高斯消去算法的变体。其中，较著名的有 Irons[6]提 分析中 （见Smith and Griffiths[2]，Phoon et al[3]， 出的用于求解对称有限元矩阵方程的波前直接求解 [4] [5] 法。波前直接求解法采用单元刚度矩阵组装与高斯 Chen ）。大型三维有限元分析软件 GeoFEA 是由 消去相互交织的求解方式，基本求解原理可以概括 新加坡吉奥软件公司开发的新一代岩土工程软件， 为 本文对 GeoFEA 的波前直接求解器、对称和非对称 ⎧ * =− −1 Krylov 子空间迭代求解器的内在算法及适用范围做 a