一种RCM有限元带宽优化改进算法.pdfVIP

第27卷第4期 计 算 力 学 学 报 ol. 27, No. 4 2010 年8 月 Chinese Journal of Computational Mechanics August 2010 :100724708(2010)0420694204 RCM 1 * 1 2 1 1 杜宪亭, 夏 禾 , 龙佩恒, 余 竹 , 王少钦 ( 1.北京交通大学 土木建筑工程学院, 北京 100044;2.北京建筑工程学院土木工程系, 北京 100044) :应用RCM(Reverse Cuthil2l Mckee) 算法进行带宽优化时存在优化结果不稳定的问题, 通过对算法进行系 统 析发现,正序排列过程中非完全依靠节点之间的拓扑关系是问题的关键。本文在考虑层、联结度判据基础 上,通过新增列高和判据进行节点正序排列, 解决了RCM 算法存在的问题, 通过实际结构算例验证了改进后的 RCM 算法的稳定性, 并获得了列高和更小的优化方案, 实现了节省计算机内存和提高运算效率的目的。 :有限元; 带宽优化; 改进 RCM 算法;列高和; 拓扑关系 :O242.21;T115 :A 度矩阵, 需要存储上三角部分元素值个数为 1 (NDF+ 1) # NDF TN= # NNODE+ 在应用有限元方法分析土木工程结构时, 结构 2 总体刚度矩阵呈对称、稀疏带状分布的特点。为了 NDF # NDF# SUM- NNODE ( 1) 减小存贮空间与提高计算效率,常采用一维等带宽 式中 NDF 为节点自由度数, NNODE 为节点数, 或者变带宽方法存储结构的总体刚度矩阵。而在 SUM 为列高和。列高和越小意味着带宽越小、需 大型结构的精细化分析与非线性分析过程中, 网格 要存储的数据越少。例如, 由N 个节点组成的结 划分后的节点与单元数量数以万计,特别是在结构 构,其相应的节点号排序方案有N! 种;但通过将所 阶段施工分析时, 结构形状与施工步骤的变化使结 有方案全部展开,分别计算相应的带宽与列高和,然 构单元与节点数量不断变化, 方程组需重新编号。 后比较各个方案的优劣得到优化方案的方法计算量 如何高效解决此类问题一直是国内外学者研究的 太大。目前,解决这一问题的经典优化方法有CM [126] (Cuthill2Mckee) 法和RCM(Reverse Cuthill2Mckee) 热点 。 [7] 本文提出一种可以降低带宽和列高和的RCM 法 。CM 法是从图论的角度出发, 将节点之间的 带宽优化改进算法, 能够节省计算机内存, 缩短求 拓扑关系,按照一定方法展成一个以某节点为根的 [10] 解时间,在解决大型复杂结构问题时, 效果特别明 树形结构 。结构中有刚度联系的两点相互为邻 显。 点, 一个节点的邻点数作为该节点的联结度, 从联