二维弹性力学快速多极边界元法的并行计算.pdf

·305· 二维弹性力学快速多极边界元法的并行计算 雷霆，姚振汉，王海涛 (清华大学工程力学系， 北京100084) Mt本文给出了二维快速多极边界元法弹性静力学问题的并行算法，包括针对分布式存储环境下的自适应树结构的划分方 案、负载均衡的实现、以及迭代算法的并行处理。在机群系统上通过采用消息传递函数库MPI实现了该算法的计算程序，数 值结果表明了该算法的较高的并行效率和可扩展性。 关键词：边界元：快速多极：分布式并行计算：机群 n姒LLELCOMPUTATIONOF2DELASTICSOLID USINGFAST MUIJIPoLEBoUNDARYELEMENTMETHOD LEI Haitao Zhenhan，WANG Ting，YAO (DepartmentofEngineeringMechanics，TsingIIuaUniversity，Beijing100084，China) this Abstract：In of2Delasticsolid fast Element paper,theparallelcomputation using multipoleBoundary Methodis the schemeof ofiterativemethodandload presented，including adaptive partition tree，parallelprocess inadistributedenvironment．Wethe balancing implement storage parallelprogramusingMPI(Message Passing acluster resultshows and Interface)insystem．Numerical scalability． highparallelefficiency Keywords：boundaryelement；fastmultipole；distributedparallel computation；cluster 1引言 边界元法作为一种广泛应用的数值方法，它只在构形的边界上进行离散，具有降低问题维数，求解精 度高等特点。但它形成的矩阵是非对称满阵，采用传统的直接法进行求解，求解算法时间和空间复杂度分 别为0(Ⅳj)和D(∥)，其中Jv是自由度数，对求解速度和规模都有很大的制约。 近年来提出的快速多极算法【1】通过引入多极展开和局部展开来逼近距离源点较远处的单元对于积分方 程的贡献，并采用不平衡的四叉树数据结构存储展开系数，从而隐式地完成了矩阵向量的相乘，结合 GMRES算法，使计算量级和存储量级均降为线性。 常规的边界元并行算法采用的是基于分而治之的分域解法忙J，把结构分为若干个子域(子结构)，通过 交界面上的位移连续条件和面力平衡条件，形成系统矩阵进行并行求解。 基于快速多极算法的边界