有限元分析的并行程序设计－－基于网络和PVM的分布式并行环境研究.pdfVIP

.108 · 工《程力学》增刊 1999年 有限元分析的并行程序设计 一基于网络和PVM的分布式并行环境 余天堂 姜弘道 任青文 (河海大学士木工程学院，南京210098) 组 架 基于网络和PVM的分布式并行环境的有限元并行计算是有限元并行分析的一个重要方向，此种 并行计算方法具有投资少、灵活性强和得到的加速比高等优点。本文给出了基于网络和PVM的分布式并 行环境的有限元分析的并行程序设计方法，有限元并行分析采用多波前并行处理的子结构并行分析方法， 界面方程的求解采用预条件共辘梯度法 (PCG)。文中还给出了有限元并行前处理方法.此种设计方法易 于实现，且胎获得较高的加速比。 关.询 并行程序设计.网络，PVM.子结构，加速比 1 引言 并行化是有限元计算方法发展的趋势1[11121。有限元并行计算的一个主要途径是利用多 波前并行处理的子结构方法，即利用多个波前对各子结构并行地进行静凝聚，再并行求解 界面方程，然后并行回代求内点位移和计算应变、应力。该方法既可获得较高的加速比， 又能有效地节省内存量.并行程序的设计与有效实现强烈地依赖于并行机硬件的计算模 型。由于网络并行计算具有巨大的计算潜能、良好的性能价格比及可扩展性和灵活的体系 结构等优点，以及以PVM.MPI,EXPRESS1211+1等为代表的一批基于消息传递的并行程序设 计软件平台的出现，使得有限元分布式网络并行计算成了有限元并行计算的一个重要方 向。本文详细介绍了基于PVM的网络并行环境下有限元并行分析程序设计方法。此种设计 方法易于实现。算例表明，此种设计方法能大大地节省计算时间。 2井行环境分析 余天堂，男，1971.10出生。博士研究生 工《程力学》增刊 1999年 (09 · 本文采用网络分布式并行环境，网络并行编程环境为基于消息传递的PVM(Parallel VirtualMachine)系统。PVM具有简洁、高效及广泛的适用性、易用性等特点，它提供 了丰富的消息传递函数库。本文采用主从式 (Master-Slave)编程模式(a3，这种计算模式 有一个控制程序称为主进程 (Master)，负责进程的生成、初始化、收集并输出计算结果， 其余的从进程 ((Slave)执行实际计算，其负载或者是由主进程分配，或者是由进程本身 分配。一个问题的求解，由主进程和从进程共同完成。进程的生成、管理及相互间的通信 等并行过程通过调用PVM库函数实现。主从式编程模式适用于粗粒度并行计算，因此在设 计算法时并行粒度的划分不可太细，通讯与运算尽量重叠。 3 多波前子结构井行分析算法 子结构并行分析法，就是并行形成和组集各子结构的劲度矩阵和荷载列阵，再对各 子结构内点方程井行地进行静凝聚，然后并行求解界面方程，最后并行回代求得各子结构 内点位移和进行应力分析。为了节省内存，子结构内点方程的静凝聚一般采用波前法，即 一边组集一边消元。 3.1多波前井行处理技术 多波前并行处理技术的关键是在各子结构内独立进行静凝聚.首先，将整体结构分 成若干个子结构，定义各子结构界面为总体波前和子结构内的局部波前，然后利用波前存 贮法在各子结构内并行组集，并对内点变量进行消元，在子结构中任一点对应的劲度矩阵 和对应的右端项的运算过程为 (1) R{dqi。{、卜:KdiK]}_dR],dq (2) 这里，K[yd}和{R,0}是各单元贡献的组集值，K[0}是将要消去的点“变量组集后的主对 角元值。在完成上述过程后，局部波前到达子结构边界。这样，从劲度矩阵中划去内点对 应的行和列和从荷载列阵中划去内点对应的行，便得到各子结构边界凝聚后的劲度矩阵和 荷载列阵。 3.2PCO法井行求解界面方理 前面得到的劲度矩阵和荷载列阵为各子结构的，为了求解界面结点位移，可采用两种 方法，一个是组集整个结构的界面方程，然后进行求解;另一个是采用迭代法求解。PCG(6] 法是一种非常有效的迭代法，运用