ID值索引技术在工程结构分析中的应用研究.pdf

·356· I D值索引技术在工程结构分析中的应用 张艳霞1，邹积麟2，张丽霞3，张艳4 (i北京建筑工程学院土术工程采．北京100044i2中盥建筑科学研究院结构所，北京10008 3．珲春市第三中学．吉林省珲春133300：4．北京建筑工程学院数学教研室．北京10004—1) 摘要：基于面向对象的程序设计思想．本文提出了基于ID值索引技术的线性代数方程组迭代解法，不仅节 省了存贮空间，而且显著提高了求解速度。 关键词：线性代数方程纽；结构剐度矩阵；迭代解法：ID值索引技术； 1 引言 在工程结构分析中，求解线性代数方程组是经常遇到的问题，也是影响结构分析质量和速度的重要问 题，因而越来越受到重视。在电子计算机出现以前，由于线性代数方程组的求解量太大，在这方面的研究 很少。计算机普及以后，尤其是在有限元出现以后，研究者在这方面作了大蠡的工作。但是现在工程界使 用的线性代数方程组的解法原理和过程自从电子计算机出现以来，几乎没有任何改变，而计算机的硬件和 软件技术都有了很大的进步，这就为数值计算技术的进一步提高提供了条件。本研究课题将传统的线性代 数方程组的解法同计算机软硬件技术的最新发展相结合，针对土木工程中的结构求解问题，对线性代数方 程组对求解方法作了研究，旨在将结构分析水平提高到与计算机软硬件发展水平相近的程度。 2线性代数方程组求解方法介绍 线性代数方程组的解法按照求解原理，大致可以分为碡大类：直接解法和迭代解法【IJ【“。直接解法是 指对方程组的系数矩阵直接进行处理，在满足非奇异的条件下，～步就可以求出满足方程组的解的方法。 送代法是指对方程缀的解进行有限次近似，最终得到满足预定精度的方法。 2．1直接解法 针对不同的线性代数方程组系数矩阵的形式，有各种不同的直接解法：可莱姆(Gram)法、高斯消去法、 LL分解法、平方根法、LDL。分解法、追赶法等等。 可莱姆法是直接解法的鼻祖，它是用矩阵行列式的形式给出方程组当解。在证明当系数矩阵不奇异时． 线性代数方程组当解存在唯一性方面发挥了重要的作用，并给出了这个唯一解。但在实际应用中，由于Gram 方法涉及到到计算耋太大，所以该方法很少采用，只用在进行理论证明时，才采用这一方法。 高斯消去法是针对满秩线性代数方程组系数矩阵的一般求解方法，它通过对方程组增广矩阵的一系列 列消去，最终得到方程组的解。高斯消去法由于计算量大，而且无法适应于结构分析中结构刚度矩阵大对 阵、稀疏特性，在结构分析中很少采用。 LU分解法是指对系数矩阵进行分解，采用回代的方法求出方程组的解的方法。任何一个非奇异的矩 阵．由矩阵运算原理，都可以分解成为A=工·U的形式。其中，￡为对角线元素为1的严格下三角矩阵， U为严格上三角矩阵。将矩阵进行Lu分解后，进行两步回代，就可以得出方程组的解。 平方根法是针对对称正定矩阵进行的方法。 LDL7分解法是针对对称矩阵进行的方法，不要求矩阵是正定的，但矩阵要求非奇异(这是对系数矩阵 基金礓目：北京建筑工程学院青年基金赉助项目 作者摘舟：张话菇(1970)．盘．讲师．工学碰士．土要从事工程结构分析方蕊的研究 譬积一(1971)．舅．耐爵究最．芏事博士，誊羹飙摹互繁精掩分孵和{博方萄的精究 张重矗(I靳7^女．庸擐教师。埘控长。主要飙事鼓掌方龌扣嗣F究 张绝‘1972)．盘．讲师．瑚士．圭薹械事教学方面的研究 ·3S7· 的一般要求)。 追赶法是针对特殊对称系数矩阵进行的方法。 2．2迭代解法 当结构刚度矩阵的自由度数很大时，由于大量运算舍入误差的影响，有时直接解法的精度不能满足要 求。尤其在进行结构弹塑性分析时，随着结构中破坏杆件数目的增加，结构刚度矩阵可能形成病态矩阵。 这种情况下，直接解法的精度远不能满足要求，甚至是无法接受的，随着电子计算机的出现和数值计算方 法的蓬勃发展，出现了大量的迭代解法。 共扼梯度法(CG法)从理论上讲是一种精确解法，因为在月步内可以得到方程的精确解。和其它的迭代 方法相比，cG法具有收敛速度快的特点。 大量的研究发现，如果在进行CG迭代求解以前，现对方程组的系数矩阵进行预处理．再对等价线性 方程组进行求解，收敛速度可以更快，这就是预处理的共扼