有限元快速稀疏静力及特征值问题求解器第十一届全国结构工程学术会议特邀报告研究.pdfVIP

!三堡垄茎!塑型兰垡兰巳——————————————————竺 有限元快速稀疏静力及特征值问题求解器 一第十一届全国结构工程学术会议特邀报告一 袁明武，陈璞，孙树立，孟昭龙 北京大学力学与工程科学系。北京100871 摘要： 本文提出了基于细胞索引概念和双向循环展开技术的快速稀琉静力及特征值问题求解器-该算法仅适用 于有限元分析．其中采用了一些软件技术使得运行时占用的内存和磁盘空间都大大减少·这样也就减少了I，o· 文中对若干实际工程问题给出了几个静力求解器在PC机上计算的时间及占用磁盘空间的性能比较-在静力分析 算例中，与传统的轮廓线求解器相比，本方法的求解速度提高了7～57倍，而存储空阃减小到原来的lf3．对于 通常的特征值问题，效率提高差不多8倍．试验表明采用了本文提出的细胞稀疏存储方案和循环展开后，无论是 对于静力问题还是通常的特征值问题，求解性能都有大幅度的提高。 关键词： 高性能计算，稀疏矩阵，有跟元分析，特征求解 1 引言 在二十世纪七十年代早期，轮廓线活动列算法被公认为是求解稀疏线性方程组j江=f的Cholesky分 解和LDL7的有效算法。这里置代表正定系数矩阵，在有限元分析中是刚度矩阵，，和z分别表示载荷向 量和未知的位移向量．但是在轮廓线活动列算法中，每一列从第一个非零元素开始存储，因而在轮廓线下 通常都还有很多零元素。所有这些零都占用存储空间并参与运算。如果考虑到采用外存，这些零导致了内 存和外存间额外的I／O。这些不必要的零浪费了磁盘空间并占用了I／o时间的相当大一部分．当线性方程 组规模很大时，轮廓线方法就会受到存储空间不足和效率低下的困扰。 同传统轮廓线和半带宽存储方案相比，稀疏存储方案只记录刚度矩阵中的非零元素，整个矩阵以一种 巧妙的方法处理成紧凑行向量的结合。运算只对非零元素实行，这样在大型问题上稀疏求解器所需的磁盘 空间比传统的采用轮廓线或半带宽存储方案的直接求解器更少，同时速度也更快。当然，稀疏求解器的编 程也比轮廓线或半带宽求解器更复杂。 在电脑上组织得很好的数值算法能极大地提高效率，但也不是影响计算效率的唯一因素。最近提出的 处理矩阵乘法的循环展开技术使得寄存器和随机提取存储器之间不必要的数据交换大大减小，它的使用在 串行和并行计算方面都促进了计算性能的提高。循环展开方法在矩阵计算时减少了循环的次数而增加了循 环的内容。在矩阵三角分解中，展开技术在不同的机器和编译器上最多能使速度提高至2～7倍。其基本 概念在㈣中就能找到。 文献91中在传统存储方案的稀疏矩阵计算中执行一维展开时对于K的因子L7；JlAT超方程的概念。在 本文中，我们把超方程的概念扩展到三维有限元分析的总体刚度矩阵并引入细胞稀疏存储方案．细胞稀疏 存储方案为二维循环展开提供了极好的数据结构，便于我们对被称为JKI厶D正7分解的三重嵌套循环进行 两次外部循环的展开。这个实质性的变化引起了计算效率和存储空间需求方面的极大改进。 2超方程和细胞存储方案 众所周知，求解一个线性系统的很大部分时间是进行如下的刚度矩阵分解： K=LU=LDLr (I) 阵。 8)，男-教授·从事计算力学、结构动力学和大型结构分析软件等方面舶研究 作者简介：塞明些(1939 1……… 妻柑登；：：篓：?：％I蛐～D～mm…chdat态兰三萝111～5：,．嗍!侉，从事计算力学、结构动力学和大登螭构分析软件等方面的研究 1 銎婪妾!}9681’，?·．苎塞盔主!蛄j列教授．从事计算力学、结构动力学和大型结构并薪软件再另蔷；油陵1 。一’～…。 孟昭龙(1977 2)，男．北京大学硕士生 ·106· 《工程力学》增H2002锌