GPU加速的三维弹性问题分析求解算法.doc

GPU加速的三维弹性问题分析求解算法 王 钢，王启富?收稿日期：2011 ?收稿日期：2011-**-**；修回日期：2011-**-**. 基金项目:国家自然科学基金；教育部新教师基金(200804871128)。王 钢（1987），男，硕士研究生，主要研究方向为CAD/CAE，数值计算，GPU并行计算；王英俊(1984)，男，博士研究生，主要研究方向为计算机图形学、边界元计算; 王启富（1961）， HYPERLINK mailto:wangqf@hust.edu.cn wangqf@hust.edu.cn通讯作者，男，博士，教授，博士生导师，主要研究方向为CAD/CAE/CAM、计算机图形学、计算机仿真；黄运保（1975），男，博士，副教授，硕士生导师，主要研究方向为计算机图形学，反求工程；魏一雄（1988），男，博士研究生，主要研究方向为计算机图形学，计算机仿真。 (华中科技大学国家CAD支撑软件工程技术研究中心 武汉 430074 ) (HYPERLINK mailto:wangqf@hust.edu.cn%20wangqf@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]是求解边界元