逆矩阵求解算法-维普网.pdfVIP

维普资讯 第23年牵第12智月 ACTAME固CH体AN力ICA学SO学LI报DASINICA DVecoe1m．2b3eNro20．402 2002 一 种改进 的 2n因子乘积形式的 逆矩阵求解算法 吕毅宁 柳玉起 胡 平 (吉林大学南岭校区汽车覆盖件成形技术研究所，长春 ，130025) 摘 要 复杂结构的有限元分析需要耗费巨大的内存和计算机时间．在此过程中，有限元 线性方程组的求解时间占很大比例，因此，发展高效的线性方程组求解算法是提高有限元分析 效率的关键．针对矩阵求逆问题 ，该文将有限元方程组的系数矩阵视为等分块矩阵，并基于等分 块矩阵的概念推广了传统的逆矩阵的 n因子和2n因子乘积形式．推广后得到的这组乘积形式 在适当条件下又分别可以退化到传统的逆矩阵的 n因子和2n因子乘积形式 ．应用基于推广后 得到的这组乘积形式的求逆算法来求解板材冲压成形有限元数值模拟 中的大型有限元线性方 程组 ，结果表明，该算法可以显著地提高大型有限元线性方程组的求解效率． 关键词 分块矩阵，逆矩阵，有限元 ，板材冲压成形 1 引言 复杂结构的有限元分析需要耗费巨大的内存和计算机时间．对采用 UpdatingLagrange有 限元方法 ¨对板材冲压成形过程进行有限元数值模拟的过程进行分析，发现求解有限元线 性方程组的时间占整个数值模拟过程的比例超过 90％．因此，提高有限元线性方程组的求 解效率对提高复杂结构的有限元分析效率至关重要 ． 求解线性方程组的数值方法有很多种．总体上可以分为直接法和迭代法两类．由于复杂 结构的变形过程往往是非线性，甚至是多重非线性的，对这样 的变形过程进行有限元分析时 遇到的有限元线性方程组的系数矩阵的性态会有很大差异．用迭代法对该线性方程组进行 求解时常常会出现不收敛或收敛速度太慢的情况．虽然直接解法需要的内存量 比迭代法要 多，但它能够精确地对这些线性方程组进行求解 ．可 以预见，随着计算机技术的迅速发展，计 算机性能不断提高，通过直接法可以求解 的方程组的规模将会逐步提高．因此，直接法还是 首选的方法 ． 特别地，对于大型有限元方程组来说，由于系数矩阵(刚度矩阵)具有稀疏、正定、对称等 特性，因此可以在高斯法的基础上发展更为有效的直接解法 ．常见的这类解法，如二维等带 宽存储的高斯消去法和一维变带宽存储的LD 法等．实际上通过对有限元离散网格模型的 节点进行适当的排序，用逆矩阵的2n因子乘积形式来求解有限元刚度方程，可以更好地利 用刚度矩阵的稀疏特性 ，合理建立存储刚度矩阵的数据结构，进一步减少 内存使用，并缩短 求解时间 ，，． * 国家 自然科学基金重点项 目、教育部重点科技攻关项 目(99O34)和 国家杰 出青年科学基金项 目 联合资助． 2001-03-19收到第 l稿 ，2001—11-05收到修改稿 ． 维普资讯 第 4期 吕毅宁等 ： 一种改进的2n因子乘积形式的逆矩阵求解算法 ‘447 ’ 在用有限元方法分析问题时，考虑到每个节点一般有多个 自由度 ，则所形成的刚度矩阵 可以视为一个特殊的分块矩阵．基于这一点，作者将传统的逆矩阵的n因子和2n因子乘积 形式l3 加以推广，得到基于等分块矩阵的逆矩阵的n因子和2n因子的乘积形式及基于此 乘积形式的线性方程组求解算法 ． 2 等分块矩阵的定义 为了下文叙述方便 ，做如下定义： 定义：将形如 A=(a) ，以子矩阵a ∈R～ ，(i，=1，2，…，n)为元素的分块矩阵 称为等分块矩阵．这里，～ 为实rn阶方阵的集合 ．