一种有限元并行计算前处理方法.pdfVIP

第 18 卷第 3 期 计 算 力 学 学 报 V o l. 18 N o. 3 2001 年 8 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS A ugust 2001 文章编号:(2001) 一种有限元并行计算前处理方法 1 2 余天堂 , 　仲济刚 ( 1 河海大学 土木工程学院, 江苏 南京 210098; 2. 河海大学 水利水电工程学院, 江苏 南京 210098) 摘　要: 本文研究一种有限元并行计算前处理方法。只需将各子结构的超单元信息发送给各处理机, 就可并行自动剖分各子结构和并行优化各子结构单元编号, 以达到极小最大波宽和波宽的均方根。 算例表明, 该方法可大大减少有限元并行计算前处理时间。 关键词: 有限元并行计算; 前处理; 子结构; 超单元; 优化; 单元编号 中图分类号: O 246　　　文献标识码:A 1　引　言 2　并行自动剖分方案及算法 ( ) 任务的划分 区域分裂或子结构划分 是并行计 2 1　方案 算的第一步。为了能获得较好的并行效果, 就需要区 设空间结构 为长方体结构, 如图 1。 域分裂算法能产生相同复杂性的子区域。区域分裂, ( 就是将计算区域分裂为与机器数相同 或为机器数 ) 的倍数 的子域, 分裂后的子域要满足各子域间的载 荷要均衡和相互间的通信量最小。关于区域分裂算 法的研究一直异常活跃, 目前已形成了很多有效的 算法, 这些算法可以分为两大类: 工程法[15 ] 和图论 法[68 ]。目前所有的区域分裂算法都需要有整个区域 ( 1) 将结构 划分为 n 个子结构 , 且 的网格单元信息和节点信息, 而要获得大型实际问 i n n 题的单元信息和节点信息又是一件非常繁琐的事。 = ∪ , ∩ = 。每个子结构由一个或多个超单元 i i i= 1 i= 1 有限元并行计算的一个主要途径是利用多波前 组成。 [9 ] 并行处理的子结构方法 。该方法既可获得较好的 (2) 为了保证整个结构上网格的连续性, 相邻 并行加速比, 又能有效地节省内存量。在进行有限元 ( 子结构界面上的节点必须保持一致 节点数及节点 并行分析之前, 首先要剖分各子结构, 形成计算数 ) 坐标相同 。确定各界面上节点数。 据; 然后优化单元编号以达