《第九章有限元线性方程组的解法》-公开课件.pptVIP

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