第九章有限元线性方程组的解法.pptVIP

前面，我们主要讨论了静力平衡问题的有限元格式。在确定了 离散所需要的单元形式后，需要进行单元特性矩阵的计算，最后由 单元特性矩阵集合而成的有限元求解方程组 R K ? ? 是一组联立的线性代数方程组。这组方程在静力平衡问题中就 是以结点位移为基本未知量的系统结点平衡方程。有限元求解的效 率很大程度上取决于这组线性代数方程组的解法。有限元分析可以 通过细分单元的网格来提高解的精度。因此，当有限元分析采用越 来越多的单元离散模型来近似实际的结构时，线性联立方程组的阶 数越来越高。曾经有相当一部分的研究工作是围绕如何有效地求解 这组庞大的线性方程组的。 （ 9-1 ） 在线性静力分析中，解代数方程组的时间在整个解题时间中占 很大的比重。而在动力分析和非线性分析中这部分比重也是相当大 的。若采用不适当的求解技术，不仅计算费用大量增加，更严重的 是有可能导致求解过程的不稳定和求解的失败。 在有限单元法中，需求解的线性代数方程组（ 6-1 ）的系数矩 阵 K ，即刚度矩阵，具有 大型、对称、稀疏、带状分布以及正定、 主元占优势 的特点。在求解方程组时必须利用上述特点，以提高方 程求解的效率，否则将无谓的提高计算费用。 线性联立方程组的解法可以分作两大类：直接解法和迭代法。 直接解法以高斯消元法为基础，求解效