论文部分内容阅读
EFG(Element-Free Galerkin)法不仅摆脱了网格的束缚,而且具有计算稳定性好、收敛速度快以及精度高等优点,另EFG法拓扑优化可有效解决有限元法拓扑优化中出现的重构网格等问题,并可有效消除棋盘格现象,所以EFG法及其拓扑优化吸引了越来越多学者们的关注。不过由于EFG法具有计算效率低的缺点,严重束缚了它们在大规模问题中的研究及应用。与此同时,随着GPU(Graphic Processing Unit)并行计算的迅速发展,其现已被广泛应用于众多领域。因此,本文针对EFG法计算效率低和耗时长等缺点,开展EFG法及其拓扑优化的GPU并行计算研究,主要内容有:(1)三维EFG结构分析方法的研究及应用。引入三角形单元(2D)或四面体单元(3D)代替EFG法传统算法中的四边形单元(2D)或六面体单元(3D)作为背景积分网格,并采用Hammer积分取代传统算法中的Gauss积分实现EFG中的积分计算,推导出相关公式且通过数值算例验证了该方案的可行性。(2)EFG法的GPU并行计算研究。首先基于交叉节点对思想,提出了一种总体刚度矩阵、总体惩罚刚度矩阵和预处理矩阵的联合组装方案及其GPU并行算法;其次提出了一种耦合GPU和PCG(Preconditioned Conjugate Gradient)法的EFG法加速算法,并给出了该加速算法的相关公式推导过程及流程图;而后通过数值算例,研究了PCG法求解EFG法总体离散控制方程的收敛速度及计算耗时,并探讨了求解域离散节点数及GPU线程块大小对加速比的影响。(3)EFG法拓扑优化的GPU并行计算研究。首先基于交叉节点对思想,提出了一种目标函数灵敏度的并行算法,并结合CUDA架构(Compute Unified Device Architecture)设计了目标函数灵敏度和OC(Optimization Criterion)法的GPU并行算法;然后提出了一种EFG法拓扑优化的GPU并行算法并给出了该算法流程图。通过数值算例验证了上述算法的可行性及有效性,获得了24倍的加速比,且探讨了设计域离散节点数对加速比及优化结果的影响。本文利用GPU并行计算显著缩短了EFG法及其拓扑优化的计算时间,提高了其计算效率,该研究对EFG法及其拓扑优化应用于工程问题具有重要的理论参考价值和工程实际意义。