有限元快速稀疏静力及特征值问题求解器研究.pdfVIP

!!!!兰垒里塑型!兰堑：堡查：塞皇童竺丝丝查兰!三翌墨皇旦兰查叁堡丝圣叁 些塞：!!!垒!坠!兰! 有限元快速稀疏静力及特征值问题求解器 袁明武陈璞孙树立孟昭龙 f北京大学力学与丁程科学系t北京100871) 摘要奉丈提出了基于细胞索引概念和双向展开的快速稀疏静力及特征值问题求解器，谊算法仅适用于有限元分析． 其中采用了一些软件技术使得运行时占用的内存和磁盘空间都大大减少，这样也就减少了I／O。文中对若干实际工程问题给 出了几个静力求解器在Pc机上计算的耗时厦占用磁盘空间的性能比较。在静力分析算例中，与传统的轮廓线求解器相比， 本方法的求解速度提高了7—57倍．而存储空问减小到原来的l，3。对于通常的特征值问题，效率将近提高8倍。试验表明 采用了本文提出的细胞稀疏存储方案和循环展开后，无论是对于静力问题还是通常的特征值问题，求解性能都有大幅度的 提高。 关键词 高性能计算，稀疏矩阵，有限元分析。特征求解 1．引言 在～十1廿纪七十年代早期，轮蛳线活动列算法被公认为是求解稀疏线性方程组KX=，的Cholesky分 解和LDL7的有效算法。这里置代表正定系数矩阵。在有限元分析L|_I是刚度矩阵，，和x分别表示载荷向 节和未知的位移向景。但是在轮廓线F通常几乎每‘列都还有很多零。所有这些零都r々用存储空间并参与 运算。如果考虑剑采删外存，这些零导致了内存和外存问额外的FO。这些不必要的零浪费了磁盘空问并 ，々JLf=i了YO时问的相当‘部分。当线性方程组规模很大时，轮廓线方法就会受到存储空间不足和效率低下 的蚓扰。 同传统轮廓线和半带宽存储方案相比，稀疏存储方案只记录刚度矩阵中的非零元素，整个矩阵以一种 巧妙的方法处理成紧凑行向最的结合。运算只对非零元素实行，这样在大型问题上稀疏求解器所需的磁盘 空问比传统的采蝴轮廓线或半带宽存储方案的直接求解器更少，同时速度也更快。当然，稀疏求解器的编 科也比轮廓线或半带宽求解器更复杂。 在电脑上组织得很好的数值算法能极大地提高效率，但也不是影响计算效率的唯一一因素。最近提出 的处理矩阵乘法的循环展开技术使得寄存器和随机存取存储器之间不必要的数据交换大大减小，它的使用 在串行和并行计算方面都促进了计算性能的提高。循环展开方法在矩阵计算时减少了循环的次数而增加了 循环的内容。在矩阵三角分解中，展开技术在不同的机器和编泽器上最多能使速度提高至2～7倍。其基 本概念在”“中就能找到。 ＼ 文献【91L}J在传统存储方案的稀疏矩阵计算中执行‘维展开时对于胃的因子工7引入了超方程的概念。 在本文中，我们把超方稃的概念扩展到三维有限元分析的总体刚度矩阵并引入细胞稀疏存储方案。细胞稀 LDL 疏存储方案为■维循环展开提供了极好的数据结构，便于我们对被称为JKI 7分解的三重嵌套循环进 行两次外部循环的展开。这个实质性的变化引起了计算效率和存储空间需求方面的极火改进。 2．超方程和细胞存储方案 众所周知，求解一个线性系统的很大部分时间是进行如F的刚度矩阵分解 置=上Ⅳ=LDLr 2002年全国搬动(诊断．模态、噪声与结构动力学)．T-程度应用学术会议论文集 北京，2002 10．16-19 ‘ ’ 阵o ． Minimum 在这部分，假定已经用如ModifiedDegree或Nested 算法得到了最初的系数矩阵K。为了便于讨论，我们假定得到的是如下6x6的总体刚度矩阵B ll a b c 11 (口b) c 22 d 22