基于细胞稀疏存贮方案有限元刚度矩阵组装.pdfVIP

基于细胞稀疏存贮方案的有限元刚度矩阵组装1 1 1 陈璞 ，陈斌 1 北京大学力学与工程科学系，北京 100871 E-mail：chenpu@ 摘 要：在现代有限元分析的直接解法和迭代解法中都需要组装刚度矩阵。本文在[1,2]的基 础上，叙述了两种基于细胞稀疏存贮方案有限元的符号刚度矩阵的组装方法。在与传统的方 法进行了比较的基础上，提出了改进效率的方法。 关键词：有限元分析 稀疏矩阵 高性能计算 中图分类号：O241.82 1．引言 近年来，稀疏解法逐渐替代了带宽解法与变带宽解法，成为了有限元分析中首选的直接 解法[1,2,3] 。不同于带宽解法或变带宽解存贮方案，稀疏解法要求建立符号矩阵 symbolic matrix ，它是总体刚度矩阵对应的Bool值矩阵。刚度矩阵中的非零元在Bool值矩阵为1，其 余则为0 。在用图讨论矩阵运算时符号矩阵又称为邻接图 adjacent graph 。在稀疏解法中符号 矩阵涉及几乎所有的步骤，如：符号组装；填充元优化；符号分解；（数值）刚度矩阵的组 装；数值分解；消元与回带等。因此，减小符号矩阵的大小在稀疏解法中具有十分重要的意 义。 传统的稀疏符号矩阵的组装方法直接采用了基于方程的数据结构[5]，符号矩阵占用的内 存为数值刚度矩阵占用内存的1/2 （4字节整数，8字节实数）。对于给定的单元类型，符号矩 阵的大小一般正比于节点数。如果分析的问题不太大，符号矩阵可以整个在内存中处理。但 模型较大时有限元的符号矩阵的可能超出内存允许的范围。于是就出现了各种符号矩阵的压 缩存贮方法，其中最著名的是文献[6]的方法。在工程有限元分析中，因为求解位移向量场 的需要，总体刚度矩阵是由一系列小的子矩阵组成。利用这一性质，文[1,2]给出了一个称为 [6] 细胞稀疏存贮方案的符号矩阵压缩方法，其压缩比大于Sherman的方法 。 本文从两个方面改进符号矩阵的组装方法：1）利用文[1,2]的超方程和细胞符号矩阵概 念直接组装细胞符号矩阵。在工程有限元分析中，它的大小仅是传统符号矩阵的1/36到1/9； 2 ）给出了两种不同的组装方案。它们分别适用于不同的内存数据安排。数值试验表明，本 文提出的方法可以在仅有28M 内存的情形下在内存快速地组装100万阶3维8节点实体有限元 方程组的符号矩阵。 2 ．稀疏矩阵的存贮方案 为了更清楚地讨论问题，取一个 6 阶的总体刚度矩阵如下: 1 高等学校博士学科点专项科研基金资助项目 座机电话号码112 - 1 - 11 a b c ? 11 ? ? ? a b c ? ? ? ? ? 22 d ? ? 22 d ? ? 33 e f ? ? 33 e f ? ? ? ? ?