局域网MPI环境实现并行有限元求解.pdfVIP

512 土木工程与高新技术 局域网MR环境实现并行有限元求解 表 勇 曹 骥 (同济大学建筑工程系 上海 200092) [提 要1本文从软硬件著手讨论了在局域网MPI环境下实现有限元计算的 并行求解。在分析了有限元的基本过程以后，本文看重讨论了并行总体刚度矩阵 组装以及结构平衡方程的并行求解，并给出了一个分析多高层结构的、可以在局 城网内用徽机实现的并行有限元程序以及算例。文末是结论和展望。 关【锐词】有限元 并行计算 局城网 MPI x 众所周知，作为现代结构分析的核心手段，有限元越来越受到人们的重视。由于常规的 有限元数值分析是一个耗时的过程;因此，计算效率能否提高，已经成为有限元分析手段能 否得以推广、普及并全面反映实际结构的真实情况的一道屏障。在过去40年发展过程中， 计算机硬件结构的冯·诺依受体系在CPU、总线等方面都制约着计算效率的提高。高性能计 算中的并行计算有望打破这一坚冰。随着计算机硬件条件的迅猛发展，目前的微机和工作 站的计算能力已经达到或者是超过上个世纪90年代小型机的水平。同时，集群技术(Clus- tering)使得局域网内的多个微机或工作站可以组建成具有专业并行计算机系统计算能力的 并行计算硬件环境。而并行支撑环境MPI(MessagePassingInterface)的出现为并行计算的 局域网化体提供了软件支持环境。因此，如何充分利用已有的软件、硬件条件，挖掘局域 网内各个闲散设备的能力，实现有限元并行局域网计算将是本文讨论的重点。 自上个世纪九十年代开始，国内外对有限元的并行计算就展开了研究工作。文献 [1 一3〕对94,95年以前的发展情况作了综述。95年以后虽然有多篇文献涉及有限元方法的 并行求解，但是在总体仍然采用以前的研究方式，主要领域有: (1)使用子结构法将问题求解域划分为若干个可并行求解区域比4，司。 (2)使用直接法或者迭代法求解结构平衡方程组阁。 (3)求解结构动力相应或者是弹塑性分析的并行算法[[7] 值得注意的是，大多数文献中实现的算法都需要专用的并行计算机，例如文献 【1, 6]中的Transputer计算机等。对局域网的软硬件条件在并行计算中的应用还未涉及太多。 一、局域网并行计算的软硬件条件 ，.硬件条件 为实现在局域网中进行有限元并行计算，同济大学HPC实验室 (高性能计算和高性 能混凝土研究实验室)组建了如下8节点的集群高性能计算系统: 第九篇 计算机应用及信息技术 513 CPU: 奔腾3代，主频667 内存: 256M 硬盘: 4GG 网卡: loom网卡 集线器: loom集线器 图11IPC实验室8节点高性能计算机配置 以上设备为市场常见设备，无须专门采购。另外，上海大学自行组装的 “自强2 128节点高性能计算系统也采用的是类似的配置。 2.软件条件 IIPC实验室使用流行的Wlndows2000作为并行计算的操作系统，并使用VC6作为编程 环境。在并行支持环境方面，HP(:实验室使用美国ArgonneNationalLab发布的MPICH.NT (版本号:1.2.3)。所有的有限元软件全部自行开发，编程语言是C十十。 二、有限元基本过程 以弹性体的静力分析为例。该过程是整个结构有限元分析的基础，主要由以下步骤完成: (1)结构的离散化 这一步是把要分析的结构划分成有限个单元体，并在单元制定位置设置节点，把相邻 单元在节点处连接起来组成单元的集合体，以代替原来的结构。 (2)选择位移函数 首先假定单元内任意一点的位移是坐标的某种简单函数，称之为位移函数。也即: