电磁场数值方法FEM2012章节.pptVIP

处理强加边界条件小结 将对角线元素 置1 该行的右端项改为强加电位值 第m行与m列的其他元素全部置零 除第m行外，其他各行右端项为原右端项减去强加电位置 与对应第m列未变换前的系数的乘积 若强加边界条件节点共有N0个，则如法处理N0次 N阶方阵 与强加边界节点电位值及有关量组成的列阵 系数矩阵的存储 利用其对称性 、稀疏性 等带宽存储 变带宽存储 具体作法 将[K]矩阵的下三角部分带内元素按行的顺序， 从各行的第一个非零元素起， 至主对角元素（包括带内的零元素及对角元素）为止的元素 （共有MD+1），依次存入一维数组KL(n)中 ★ ★ ★ 说明KL(n)数组中各元素在原矩阵[K]中的行、列位置 指针数组L(n) ﹢ 存储原[K]阵中主对角元素 在KL(n)中的位置 需要计算的量 [K]的下三角阵各行的半带宽MD 第i行的半带宽MD=i-min(j) 主对角元素在KL(n)的地址JO JO=L(i)=L(i-1)+MD+1 非主对角元素在KL(n)的地址JO JO=L(i)-(i-j) 非主对角元素的地址为该行主对角元素的地址向前推移(i-j)个位置 所有以节点i为顶点的各单元上节点总体编号的最小值 前一行主对角元的地址 行号i MD 在KL(n)中 位置（即L(n)值） 在KL(n) 中位