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

前面，我们主要讨论了静力平衡问题的有限元格式。在确定了离散所需要的单元形式后，需要进行单元特性矩阵的计算，最后由单元特性矩阵集合而成的有限元求解方程组; 在线性静力分析中，解代数方程组的时间在整个解题时间中占很大的比重。而在动力分析和非线性分析中这部分比重也是相当大的。若采用不适当的求解技术，不仅计算费用大量增加，更严重的是有可能导致求解过程的不稳定和求解的失败。 在有限单元法中，需求解的线性代数方程组（6-1）的系数矩阵K，即刚度矩阵，具有大型、对称、稀疏、带状分布以及正定、主元占优势的特点。在求解方程组时必须利用上述特点，以提高方程求解的效率，否则将无谓的提高计算费用。 线性联立方程组的解法可以分作两大类：直接解法和迭代法。 直接解法以高斯消元法为基础，求解效率高。在方程组的阶数不是??别高时（例如不超过10000阶），通常采用直接解法。当方程组的阶数过高时，由于计算机有效位数的限制，直接求解法中的舍入误差，消元中有效位数的损失等将会影响方程求解的精度，此时可用迭代解法。 本章主要讨论几种常用的比较有效的直接解法，LU分解法和波前法。; 有限元法中，线性代数方程组的系数矩阵是对称的，因此可只存储一个上三角（或下三角）矩阵。但是由于矩阵的稀疏性，仍然会发生零元素占绝大多数的情况。考虑到非零元素的分布呈带状特点，在计算机中系数矩阵