三维弹性静力学边界元系数积分的并行计算.docVIP

GPU加速的三维弹性问题分析求解算法 王 钢，王启富(，王英俊，黄运保，魏一雄 (华中科技大学国家CAD支撑软件工程技术研究中心 武汉 430074 ) wangqf@hust.edu.cn) 摘 要：目前三维弹性问题的边界元法分析已经受到广泛关注与研究，然而边界元法的系统方程存在非对称稠密特性，求解过程需要消耗大量计算资源和时间。为高效实现边界元法分析，本文提出了双重信息补偿的广义极小残值法(GMRES-DC)，同时结合GPU的并行计算能力，完成了算法的并行实现。实验结果表明，并行GMRES-DC算法能够高效的完成边界元法的求解过程。 关键词：三维弹性问题；边界元法；GMRES；双重信息补偿技术；GPU 中图法分类号：TP312 GPU-Accelerated Three-Dimensional Elasticity Analysis and Solution Algorithm Wang Gang, Wang Qifu, Wang Yingjun, Huang Yunbao, Wei Yixiong (National CAD Support Software Engineering Research Center, Huazhong University of Science and Technology, Wuhan 430074) Abstract: Boundary element method(BEM) to analyze three-dimensional elasticity problems has been widely studied. But the system equation of BEM is nonsymmetric and dense, the solving process costs considerable memory and time. In order to efficiently implement BEM analysis, Generalized Minimal Residual Method with Dual Compensation(GMRES-DC) algorithm is applied. Combined with the parallel computing power of GPU, we carry out the algorithm in parallel. The example results show that the parallel GMRES-DC algorithm can efficiently accomplish solving process of BEM. Key words: Three-dimensional elasticity; Boundary element method; GMRES; Dual compensation; GPU 三维弹性静力学问题是工程技术领域中应用最为广泛的一类问题，目前通常采用有限元法(Finite Element Method, FEM)进行分析求解，边界元法(Boundary Element Method, BEM)是继有限元法后一种重要的数值方法，可将问题的偏微分方程边界值问题转化为边界积分方程，相比于有限元法，具有降低问题维数，缩减未知数个数，减小网格剖分工作量，同时能够方便处理无限域问题等优点，引起了广泛的关注与研究[1]。然而离散后的边界元积分方程通常会产生大规模、非对称、稠密的线性系统方程，求解过程需要消耗大量计算资源和时间，限制了边界元法在实际工程中的应用。 大规模方程求解通常使用迭代法，最常用的一类迭代方法是Krylov子空间法[2]，其中广义极小残值法[3]是求解边界元系统方程较为有效的方法[4,5]。GMRES算法具有良好的收敛特性，在实际应用中通常结合重启技术（即迭代每进行m步重新开始一轮）来控制算法的存储计算量。但是，使用重启技术丢弃了重启前Krylov子空间中可能包含的有价值的信息，破坏了算法的收敛性。目前，针对这一问题，通常所采用的解决方法是信息补偿技术，主要分为两类。第一类是基于特征向量的补偿技术，将重启前求解过程中得到的特征谱信息添加到下轮循环中。R.B.Morgan[6]首先提出了扩张Krylov子空间的GMRES算法(AGMRES)，将重启前所得多个小近似特征值对应的特征向量直接添加到下轮迭代循环产生的子空间末尾。第二类是基于误差向量的补偿技术，将两轮迭代循环得到近似解间的差值信息添加到下轮循环中。A.H.Baker等人[7]提出了松弛的GMRES算法(LGMRES)，在每次重启时保留近似误差向量，并添加到下轮迭代循环产生的子空间末尾。 随着