有限元网格节点编号优化方法.pdfVIP

第24卷 第2期 太 原 重 型 机 械 学 院 学 报 Vo1．24 No．2 2003年6月 JOURNALOFTAIYUAN HEAVY MACHINERY INSTITUTE Jun．20o3 文章编号：1000—159X(2003)02—0124—04 有限元网格节点编号优化方法 郭晓霞，刘建生，陈慧琴 (太原重型机械学院，太原 030024) 摘 要：本文采用的有限元网格节点编号优化方法是一种归纳、演绎方法，通过对节 点编号最优的典型网格的各种拓扑特性进行分析，从中归纳出三条规律，然后再利用该规 律去反复调整节点编号未优化的网格，直至网格中的带宽值不变为止。算例的结果表明， 对于任意划分的有限元网格，利用该方法对节点重新编号，可得到较小的半带宽值，节省 了计算机的内存空间，提高了计算效率。 关键词：有限元网格；节点编号；带宽减小 中图分类号 ：TG302 文献标识码 ：A - - 4--+ 一—·卜-—·卜-+ -—·卜-—·卜-+ 一+ -—·卜 · D，可以表示为 0 引言 D = (n 一r,／i)P 式中：n一和n 分别为单元e内节点的最大最小整 在有限元分析中，各单元刚度方程求出后，要 体号码，P为节点的自由度。由上式可知，总刚度矩 按整体 自由度集合的原则，把它们集合起来形成整 阵的带宽和所有单元的节点编号的最大差值有关。 体总刚度方程，该方程是一个线性方程组，其系数 对于同一个有限元网格，不同的节点编号方式所对 矩阵即总刚度矩阵是对称的，因此在程序设计时可 应的总刚度矩阵的带宽值可能不一样，只有当节点 只存储一个上三角 (或一个下三角)矩阵。但 由于 的编号适当时，才可得到较优的带宽值，使刚度矩 总刚度矩阵的稀疏性，所 占的内存空间还是比较 阵中的非零项尽可能排列在仅靠对角线的两边，从 大。因此，根据刚度矩阵中非零元素呈带状分布的 而节省计算机的内存空间和计算时间，提高计算的 特点，即上三角矩阵中的全部非零元素都包含在 以 速度和效率 。 对角线为一边的等带宽中，存储时常采用二维等带 本文提出了一种简单而有效的节点号码重编 宽存储 ，存储所 占的内存空间与带宽成正比。对于 的算法 j，它可以适用于任何类型单元的二维或 n阶矩阵，通常取最大半带宽为带宽D，即矩阵中每 一 三维网格，下面以四边形单元来说明算法的基本原 行从主对角线元素到最后一个非零元素个数中 坪 的最大值，其具体确定方法是先求出任一单元的最 大半带宽D，，然后再计算出所有单元中的最大值， 收稿 日期：2002—07—05 作者简介：郭晓霞 (1976一)，女，太原重型机械学院研究生，研究方向为金属塑性成形数值模