第八章线性方程组求解.pdfVIP

计算力学 第八章 线性方程组求解 引言 • 解法分类： • 1直接解法；以高斯消元法为基础，低自 由度时效率较高。 • 2 迭代解法：CG为代表的迭代解法，适 用于高自由度，精度要求高的计算模型。 特别适合实体单元的分析。 直接解,波前法, 内存限制 直接解,迭代法难 以收敛的模型,计 算速度较快。 迭代解法,单自由 结点模型,计算速 度要求高的模型。 迭代解法,速度 快，易迭代,范 围广 迭代解法,速度快 （实体模型）， 范围广. 类似PCG，8 个处理器 网络并行计算 ANSYS提供的方程组求解方法 8.1 矩阵的存储方案 ——对称假设 1. 二维等带宽存储 D 半带宽：D * * (i,j ) (i ,j ) 二维数 组的界 n ×D i*=i 二维等带宽存储在局部带宽不大时 * 效率较高，程序易于编制，但局部 j =j-i+1 带宽大时的效率不高。 2. 一维变带宽存储（轮廓线法） 一维数组A [m]：存储轮廓线内 k11 k12  元素    k22 k23 k24  列指示数组M [n+1] ：存储各列  k k  主对角元在A数组中的位置  33 34  k  44   kn 2,n  每一列的列高：N =M [i+1]-M [i]   i  kn 1,n   knn  kij A [N] N=M [j ]+j -i 一维变带宽存储的与二维等带宽存 储相反，程序编制困难，但是在局 部带宽较大时的计算效率很高。 起始行号 m =i-N +1 i i 3. 索引存储 k11 k12  一维数组AH [m]：按列存放k 中   所有非零元素 k k