有限单元法的波阵解法.pdfVIP

清华穴字字于民 JOURN AL OF TSINGHUA UNIVERSITY Vo l. 23 No.1 1983 第 23 卷第1 期 有限单元法的波阵解法 清华大学水利系周维垣 水利电力部成都勘测设计院戴国平 摘要 本文叙述了波阵法的理论分析及使用方法。波阵法对解大型线性方程组有其特殊的 优越性，水工结构及岩体力主问题常常存在很大的刚度矩阵， 以致于此类问题的正确解 答及精度受计算机容量的限制J 波阵法在合理使用机器内存.主主与外存交换方面提供了 一个方便的方法。此法主要点是在计算机内只要求存储与消去现役行有关的元素， 机器一边消元，一边形成刚度系数，因而内存要求不多，一般中小机器都可使用。 前 目 主 在水工结构及岩体力学问题的有限元计算肘，常常要解算大型线性方程组。方程组 的系数矩阵要求极大的计算机内存容量，因此往往影响解题的精度。使用合理有效的解 也是研究有限元的重要课题之一。波阵解法又叫战前解法，它是指方程组的一种解 法。通常方程组的解法有两种:一是直接解怯，这里主要指的是高斯、乔累斯基的系数 矩阵三角分解法;二是选代解法。迭代解法的优点是需要计算机的内存容量较小，结点 编号可以灵活。缺点是计算工作量较大，有的题目收敛速度较慢。直接解法的优缺点正 好相反，优点是计算工作量小，fr{;点是需用的计算机内存容量较大。目前国内一般的计 算机内存量较小，在遇到解大型题目时，就出现小机解大题的矛盾。另外，使用直接 解法时，对结点的编号要给予注意，即所谓优化带宽的问题。但实际上由于各种具 体原因，往往不可能做得很好，刷品对于大型的题目更是如此。以上问题在战阵解法 中都可得到较好的解决。下面介绍波阵法的原理与方法。 一、按非自然顺序作的高斯消元法 众所周知，高斯消元法是解线性方程组的基本方法。改造的高斯消元法一一主元素 站，不是按未知数的自然顺序作消元，而是按照主元绝对值大小的次序进行消元的。选 择某种消元顺序，可以在每一步消元时，所修改的方程系数的个数和自由项的个数达到 可能的最小数目。被修改的系数需要保存在内存中，不修改的系数可以不必占用内存容 本文于 1982 年 3 月收到 学 学 ~挡 大 报 渭 一- 量。当然随着消元过程的进行，被修改的保存在内存中的系数也不断有所变化，但这个 容量总比总刚的全部系数要少得多。这就是战阵解法的基本出发点。所以，波阵解法实 质上就是以存储量最小为准则的主元素法。 已知下列方程组z K (O) u =F CO) ( 0 ) 设有n 个未知数，则 K(O) 为 nXn 阶方阵。 现假设 n 个未知数的消元顺序为 p p z, , 用 Pn J 第