[2018年最新整理]FE-Ch-6线性代数方程组解.pptVIP

弹性（材料线性）小变形（几何线性）问题 不随 的变化而变化。 是一组线性代数方程组。 线性代数方程组的解法 直接解法：对一给定的方程组，事先按一定的算法步骤计算出它所需要的算术运算操作数，直接给出最后的结果。 适用于小于100000阶 的方程组。 以高斯消去法为基础，求解效率高。 一. Gauss循序消去法___直接解法的基础 对于n阶线性代数方程组的求解： 第一步:消元 n阶线性代数方程组，共需进行 n-1次消元。 第m次消元： 对i 行 m列( i m )的消元公式，（将 m列从m+1个元素开始变为0） 消元过程可写为： 2 . 三角分解递推公式： 所以 任一的元素： 采用上三角储存 所以 原矩阵 中的 记为： 即二维存储中新的编号 为： 在m次消元过程中， 二维等带宽存储 (n?ND) 二维等带宽存储 存储方法： ND 相关结点：所有与结点i共单元的结点， 称为结点i的相关结点。 如果结点j是结点i的相关结点， 则 ， 不是结点i的相关结点， 则 改变矩阵中元素的编号 2. 一维变列高存储 M：[1, 2, 4, 6, 10, 12, 16, 18, 22 ] 一维存储必须按列存储，辅助数组M：[n+1] 主对角元素的位置 ND 列高：m(j+1)-m(j) 在一维存储中的位置： m(j) + j - i j列最上面非零元素的行号 ： r0 = j - m(j+1)+m(j) + 1 3. 两种存储的比较 二维等带宽 内存 cpu 较多 程序编制：简单 一维变带宽 内存 cpu 较少 程序编制：较复杂 变列高 找元素 存储、消元 所以，在计算机内存允许时，一般常用二维 等带宽存储方法。 二维等带宽存储的Gauss循序消去法 一. 工作三角形 由于系数矩阵呈带状 每次消元只涉及包括 主元在内的一个三角 形内的元素，称为工 作三角形。 每一次Gauss消元,最多 修正元素的个数为： 二维等带宽存储工作区为三角形 一般Gauss消去法回代公式： 二维等带宽存储回代公式： 二. 有限元中的三角分解公式： 除对称外，还有带状、稀疏， 1. 二维等带宽存储 (n?ND) 采用按行分解 I=i , J=j-i+1 新的循环界： r=max(j-ND+1, i-1) 2. 一维变列高存储 采用按列分解，j 列高：m(j+1)-m(j) j列最上面元素的行号 ： r0 = j - m(j+1)-m(j) + 1 在一维存储中的位置： m(j) + j - i 分解递推公式： 2. 一维变列高存储 三.Gauss消去法与三角分解法比较： 相同点：1). 二种解法都是基于Gauss消去法的原理 2). 工作区的运动方向一致 不同点：3). 但具体的计算 步骤（次序） 不同，因此引 起的工作区不 同。三角分解 占机时少。 §6.5 分块分解（主要解决计算机内存容量不够时 的一种方法） 分析与思考： 在有限单元法中，一般解题过程是: 先集成，后消元。 Gauss消去法的特点： 1. 第m次消元过程中，所涉及的元素仅在三角形工作区内 消元行元素 (i)第m行元素必须集成完毕。 （主消元行） 我们注意到： (ii)待消元行可以没有集成完毕。 集成与消元可以同时进行。 但必须在作为消元行前集成完毕。 2. 三角分解递推公式： 或： 在整个消元过程中， 工作区自上向下运动。 为分块解法奠定基础。 前面的元素不再参加消元 后面的元素尚未参加消元 由前面内容知： 对第 i 行元素 即： 二.分块解法： 1.边集成边消元 2.假设允许内存NA： 公共块 半带宽 2.先分块集成，后分块消元 相关块数： 每块允许元素行数NQ: 此时采用在内存放两小块： K1块、 K2块 若带宽ND比较大，上述方法不能进行