有限元若干问题并行处理.pdf

西北丁业大学颀卜论文 第一章绪论：有限元法的并行处理技术 1．1引言 随着科学技术和工程技术领域的不断丌拓和发展，出现了各种各样大型和超 大型的计算的复杂结构。这些结构不但是具有很大自由度的组合结构，还包含非 线性本构关系、随机载荷和复杂的边界条件等多种因素。对这些结构的分析必须 借助于高阶的数值分析模型和大规模数学计算，以保证数值解的精度。而传统的 串行有限元分析方法受求解此类问题所需的计算时间长和内存储量大的限制，难 以实现这些结构的计算，使其应用受到束缚。因而需借助于已出现的各种各样的 并行机进行结构有限元并行处理。 目前己提出的有限元法并行化的集中途径包括：①硬件方面，在并行化时利 用通用的硬件或设计专用的硬件，如各种向量机，有限元机…；②算法方面，由 于有限元结构分析中最耗时的计算阶段是重复求解大型线性代数方程系统。因 此，人们一直在研究并行算法，这些算法包括并行方程求解器及基于区域分裂的 并行算法、运算子分裂技术或单元接单元策略，以及在有限元分析过程中的各个 层次探索提高并行度的各种策略和技术【2】。 1．2有限元分析的并行处理技术 并行有限元分析理论与并行计算机的发展是密切相关的。早期的并行有限元 大部分都使用向量流水处理机、阵列处理机等专业并行机。自从美国国家宇航局 限元并行计算的研究已经历了近30年的发展历史，并取得了众多令人瞩目的研究 成果，包括有限元并行算法研究和并行程序设计等成果。程建钢【31，李晓梅‘4】曾 对1995年以f；{『的并行有限元发展作了综述，主要讨论了并行有限元在向量机、阵 列处理机、共享存储多处理机等上的应用情况。 20世纪90年代，网络技术的快速增长对并行计算机的发展产生重大的影响， 例如，用专门网络紧密耦合多个计算结点而构成的大规模并行处理机系统大规模 并行机(massivelyprocessor，MPP)，以及由多台高性能计算结点经计算 parallel 机网络连接组成的集群处理环境机群并行机(cIuster 典型的分布存储计算机系统逐渐成为高性能计算机领域的两大重要研究方向。尤 第一帝绪论：自限几的并}。+处世救术 其是cOw具有投资风险小、结构灵活、可扩展性强、通用性好、异构能力强等 较多优点，已成为高性能计算领域的一个新的发展热点。 1．2．1分布存储环境并行有限元的研究进展 在分布存储并行环境下，目前的有限元的并行策略主要有： ①基于线性方程组并行数值求解器。对传统有限元最耗时的过程——系统方 程组的求解进行并行化以加快整个分析过程的求解时间： ②并行子结构法。根掘传统区域分裂算法的理论划分子域，在边界结点上建 立整体结构的平衡方程，降低问题的自由度： element)法。在单元级上实施有限元整体运算的并行 ⑧并行EBE(elementby 化迭代算法，避免组集总刚矩阵； elementa11d ④其他方法，如FETI(finitetearinginterconnecting)法。 1．2．2基于系统方程求解的并行有限元 并行有限元最直接和最仞的实现是利用数值并行求解器求解系统方程组，其 原理如图l，1所示，即单元分析和系统方程的形成是串行过程(在一台处理机完 幽1．1 成)，而系统方程组的求解是并行过程。采用这种方案的最简单做法是利用一些 pq北T业大学硕卜沦文 L 现成的线性方程组软件包。MaIluel 一个钢筋混凝土空问框架结构非线性分析的并行程序；李根国等【4】将商业著名有 限元软件NAsnuN在“神威I”超级计算机系统进行了移植和并行计算功能的 二次开发。他们的丌发思路是将有限元的求解过程(线性方稃组的求解或矩阵的 特征值问题)进行并行化，并采用了PETsc软件包，前后的处理工作由NAsTRAN 软件完成，模块l’日J通过文件进行数掘交换。 另一种做法是将刚度矩阵按行进行“分割”，各处理机存储与每块范围相应 的有关数据(单元、结点、荷载等)，这样单元分析和系统方程组装过程可并行进 行，不需任何数据通信(当然是以少