一种有限元网格节点编号的优化算法.pdf

一种有限元网格节点编号的优化算法1 张媛媛1，侯华2，程军2 1.中北大学电子与计算机科学技术学院，太原（030051） 2. 中北大学材料科学与工程学院，太原 (030051） E-mail：yuan_yuanzhang@163.com 摘 要：为提高有限元的计算效率，减少计算数据所需在存储量，本文提出的节点编号优化 算法是一种归纳、演绎方法。通过对在网格半带宽最小时的节点排列的一些特征进行分析， 从其中归纳出三条规律，通过大量的算例来验证该规律的正确性和可靠性，结果表明，对于 任意划分的有限元网格，利用该方法对节点重新编号，可得到较小的半带宽值，节省了计算 机的内存空间，提高了计算效率。 关键词：有限元；网格；节点编号；半带宽；优化 1 引言 有限元分析时经常会碰到具有稀疏矩阵刚度的方程式。减少稀疏矩阵的带宽可以减少有 限元的计算时间，和减少占用计算机的储存空间。在有限元的计算过程中.可通过节点重编 号来减少带宽。对于节点重编号的研究较多，最典型的有：Cuthill和Mckee[1]基于图论的节 点分层编号方法，可在一定程度上减小刚度矩阵的带宽，此方法后来称为CM方法。Gibbs[2] 提出了更加有效的基于图论的减少稀疏矩阵带宽和外形的方法，此方法后来称为GPS方法。 Akhras[3]通过计算节点商方法来缩小带宽,该方法编程方便，但他的方法过分依赖初始编号， 算法不稳定。此方法后来称为AD方法。欧阳兴提出了前沿法与矩形法，该法特别适用于平 面单元和壳单元[4]。 2 优化算法分析 在对物体进行有限元分析时，首先要进行网格划分，通常初始划分网格和确定节点序号 是按照某种规律进行，这样有利于计算机自动划分网格。一般来说，编排出来的网格和节点 所形成的有限元总刚度矩阵，有时半带宽值很大，需要对网格节点进行重新排序，以达到半 带宽值最小的目的。本文采用的有限元网格节点优化排序方法是一种归纳、演绎方法。即对 一种最佳的网格节点编号图进行分析，根据其单元上节点的排列，达到构成刚度矩阵半带宽 值最小时的一些特征进行分析，从其中归纳出若干条规律，通过大量的算例来验证该规律的 正确性、可靠性和优越性。该算法的特点，可通过分析图 1 所示具有较优节点编号图来予以 说明。 1本课题得到山西省攻关项目铸造工艺设计优化集成系统软件（基金编号：51128）的资助。 -1- 3 6 9 12 15 ⑤ ⑥ ⑦ ⑧ 2 5 8 11 14 ① ② ③ ④ 1 4 7 10 13 图 1 具有较优节点编号图