利用GPU加速Abaqus求解结构问题.pdfVIP

利用GPU 加速Abaqus 求解结构问题 陈小龙 郭培卿 上海超级计算中心 上海 201203 摘要：SIMULIA 最新发布的Abaqus 6.11 中加入了GPU 加速Abaqus/Standard 求解模块的功能， 从而使用户在面临规模日益增长的计算模型时，Abaqus 能够在并行计算的基础上，发挥GPU 浮 点计算的特长，进一步提高问题的求解效率。本文在上海超级计算中心“魔方”平台上通过Abaqus 6.11 求解典型结构力学问题，针对Abaqus/Standard 求解模块的GPU 加速性能进行了测试分析。 结果表明，在CPU 并行规模小于16 的情况下开启GPU 加速后，多种类型结构问题的求解效率 均有不同程度的提高，并且结构静力学与非滑移接触问题的加速效果最为明显，个别单元类型问 题GPU 加速无法得以体现。随着CPU 与GPU 异构协同处理技术的不断发展，GPU 加速在Abaqus 软件中必将发挥越来越重要的作用。 关键词：Abaqus/Standard ，结构力学，GPU 加速 一、Abaqus 与GPU 通用计算简介 Abaqus是DS SIMULIA旗下一套功能强大的工程模拟有限元分析软件，可以分析复杂的固体 力学结构力学系统，特别是能够驾驭非常庞大复杂的问题和模拟高度非线性问题[1] ，在大量的高 科技产品研究中发挥着巨大的作用[2-4] 。Abaqus 中包含Abaqus/Standard 和Abaqus/Explicit两个主 求解器模块，Abaqus/Standard提供并行稀疏矩阵求解器，对各种大规模计算问题都能十分可靠地 快速求解。在解决实际问题时，计算规模往往非常庞大，除了强大的并行功能外，Abaqus/Standard 还包含许多新颖的求解技巧来提高求解速度。Dassault Systèmes于2011 年5 月发布了最新版本 Abaqus 6.11，在最新版本中增加了GPU加速求解功能，能够利用英伟达[5] 的Quadro或Tesla系列 GPU，实现由GPU加速的CAE计算。 目前，通用科学与工程计算正由CPU 中央处理向CPU 与GPU 协同处理的方向发展。在 Abaqus/Standard 6.11 中，CPU负责控制作业的启动和终止，GPU负责求解繁重的计算任务，并将 结果通过PCIE ×16 接口返回给CPU 。根据官方数据[6] ，对于Abaqus/Standard求解器时间 （CPU Time ）占主导地位的相关问题，GPU加速能够明显提高Abaqus 的计算效率。 二、测试工况 本文选取结构静力学、材料非线性以及状态非线性（接触）三类典型结构问题，其中静力学 问题采用两种不同单元类型，设计四个测试工况，工况设置如表1 所示。 表 1 测试工况表 Case 问题描述 力学问题分类 单元类型 DOFs a 悬臂梁 静力学 C3D8R 713,691 b 法兰盘受拉 静力学 CAX4R 577,059 c 斜板受压 材料非线性 S8R5 1,027,950 d 实体碰撞 接触 C3D8R 606,600 典型的计算模型示意图如图1 所示。 1 图1 测试模型 （Case a ） 在静力学问题求解过程中增量步与迭代步数目均为 1，非线性问题求解时的最大增量步长为 0.2，在20 个增量步之内完成。对每个模型首先在不同的并行规模下进行求解（CPU核心数分别 为 2 、4 、8、16），记录下CPU求解所需时间TCPU ，然后在Abaqus/Standard求解器中开启GPU加 速选项，再针对同一问题进行求解，记录下带有GPU加速下的计算耗时TCPU+GP