有限元第12讲.pdfVIP

结构有限元分析 结构有限元分析 北京交通大学土木建筑工程学院 北京交通大学土木建筑工程学院 2008 年6 月 2008 6 年 月 第12讲（2学时） 第12讲（2学时） 第三章 线性方程组的 第三章 线性方程组的 若干直接解法 若干直接解法 T 3-4 实对称正定矩阵的LDL T分解 3-4 实对称正定矩阵的LDL 分解 T 3-5 等带宽存储下的LDL T分解算法 3-5 等带宽存储下的LDL 分解算法 T 3-6 变带宽存储下的LDL T分解算法 3-6 变带宽存储下的LDL 分解算法 T 3-4 实对称正定矩阵的 T LDL 分解 3-4 实对称正定矩阵的LDL 分解 高斯消元法的优缺点： 高斯消元法的优缺点： 程序简单，有较好的计算精度，使用广泛； 程序简单，有较好的计算精度，使用广泛； 消元过程结束之后，系数矩阵A 中的内容已经改 消元过程结束之后，系数矩阵A 中的内容已经改 变。这对于需要多次求解方程组A X = B 是不利的。 变。这对于需要多次求解方程组A X = B 是不利的。 i i i i 特别是，在有些问题中，B 需在算出X 后方能确 特别是，在有些问题中，B 需在算出X 后方能确 i +1 i i +1 i 定，这就需要对同一个矩阵A 进行多次同样的消元 定，这就需要对同一个矩阵A 进行多次同样的消元 运算，从而使得算法效率显著下降；并且还需设置 运算，从而使得算法效率显著下降；并且还需设置 一个数组用以存储A 的原始内容，这将导致计算机 一个数组用以存储A 的原始内容，这将导致计算机 内存容量不足的困难。 内存容量不足的困难。 3-4-1 对称正定矩阵的Cholesky分解 3-4-1 对称正定矩阵的Cholesky分解 定义： 定义： 任何自然数都可以分解成有限个质因数的乘积。 任何自然数都可以分解成有限个质因数的乘积。 与此类似，在一定条件下,将某个矩阵分解成两个 与此类似，在一定条件下,将某个矩阵分解成两个 或多个矩阵的乘积，这叫做矩阵的因子分解。 或多个矩阵的乘积，这叫做矩阵的因子分解。 4 2 12 2 0 0 2 1 ⎤⎡6 ⎤ ⎡ ⎤ ⎡ 1 2 1 0 1 2 ⎡ ⎤ ⎡