结构分析中的gpu并行快速多极边界元法分析-gpu parallel fast multipole bem analysis in structural analysis.docxVIP

摘要设计是企业产品创新的源头，是制造业核心竞争力的关键所在。产品创新设计不仅 要满足产品结构形状的需求，而且还要满足结构性能要求。当前，大量企业产品设计应 用三维 CAD 软件进行几何建模，采用基于限元法的 CAE 软件进行结构性能分析，但 CAD 模型与 CAE 模型之间的模型转换、模型简化、网格剖分等前处理过程需耗费大量 的时间，影响产品设计效率。为此，本文重点对边界元法和快速多极边界元法进行深入 研究，提出一种 GPU 并行的快速多极边界元方法，并应用于产品结构性能分析。该方 法能有效简化传统有限元分析前处理过程，提高产品设计效率，同时也为新一代产品设 计 CAD/CAE 软件提供一种可行的集成化方法。论文主要研究工作如下：(1)对三维弹性力学问题的边界元法进行研究，包括边界积分方程建立、单元积分 方法、角点问题处理、边界面应力计算及 GMRES 迭代求解算法等。针对角点处面力不 连续问题，提出了一种边界条件相关的混合单元法，并利用三维模型 BREP 表达中角点 拓扑关系实现了混合单元的自动生成。与现有混合单元角点处理算法相比，该方法仅在 位移约束角点处采用非连续单元，有效减少非连续单元引起的附加自由度，降低结构分 析问题求解规模。(2)对快速多极边界元法的算法原理进行研究，提出了一种节点单元双重信息自适 应结构树的构建方法，实现了高阶边界单元积分的快速计算。基于该方法建立的双重信 息快速多极边界元法可将边界元法的时间和空间复杂度由 O(N2)降到 O(N)，且单元积分 计算量仅为采用全局节点法和节点分片法的快速多极边界元法的三分之一。此外，将快 速多极边界元法与给定边界条件结合，提出了一种适用于快速多极边界元法的刚体位移 特解法，解决了 1/r2 奇异积分和自由项系数的求解问题。(3)针对快速多极边界元法中多极展开向局部展开系数传递(M2L)计算过程存在效率低的问题，本文对基于指数展开的新型快速多极边界元法进行探索，研究表明该方法 在展开阶次较大时才能达到高的计算精度及显著的加速效果，且需额外增加存储。为此， 本文进一步研究子层结点向父层结点越层传递的 M2L 优化改进方法，实验结果显示， 越层传递 M2L 方法不需增加额外内存，且加速效果与展开阶次无关，有利于结构性能 分析的快速计算。(4)充分利用边界单元及自适应结构树结点的固有并行特征，提出了一种基于CUDA架构的自适应快速多极边界元GPU并行算法，对快速多极边界元法中多极展开、多极展开系数传递、多极展开向局部展开系数传递、局部展开系数传递以及近场节点单 元积分计算进行加速。实验结果表明，该算法不仅具有显著的加速效果，而且对不同形 状的三维模型均具有良好的适应性，有效提高产品结构性能分析效率。最后，以上述理论研究为基础，对集成化CAD/CAE产品设计软件技术及系统架构 进行研究，结合现有自主知识产权的三维参数化特征建模软件InterSolid，采用Visual C++ 集成开发环境，研制开发了集成化CAD/CAE设计分析原型系统软件。并以此为基础， 针对不同形状、不同复杂程度的典型三维产品实例进行结构分析工程计算验证，实验结 果表明，本文所提出的理论及算法具有计算效率高、求解规模大、适应性强等优势，具 有良好的工程应用前景。关键词：CAD/CAE 集成，结构性能分析，GPU 并行计算，弹性力学问题，边界元法， 快速多极AbstractDesign is the source of product innovation in enterprises and plays a key role in the core competitiveness of manufacturing industry. The innovative design of products not only satisfies the demands of structural shape, but also meets the requirements of structural performance. Nowadays, CAD and finite element method based CAE softwares are used in geometrical modeling and structural performance analysis respectively in most product designs of enterprises, but it takes much time in the pre-processings between CAD models and CAE models such as model conversion, model sim