基于CUDA-GPU架构的超二次曲面离散单元并行算法

作     者：王嗣强 季顺迎 WANG Siqiang;JI Shunying

作者机构：大连理工大学工业装备结构分析国家重点实验室 

出 版 物：《应用数学和力学》 (Applied Mathematics and Mechanics)

年 卷 期：2019年第40卷第7期

页      面：751-767页

核心收录：

学科分类：080103[工学-流体力学] 08[工学] 080102[工学-固体力学] 0801[工学-力学（可授工学、理学学位）] 

基　　金：国家重点研发计划(2016YCF1401505) 国家自然科学基金(11572067 11772085) 

主　　题：CUDA-GPU架构 超二次曲面单元 非球形颗粒 离散单元方法 并行算法 

摘      要：大规模离散元的并行计算通常基于理想的球体单元,然而自然界或工业生产中普遍存在的是由非球形颗粒组成的复杂体系,其在不同空间尺度下的动力学行为及力学性质与球形颗粒具有显著差异.基于连续函数包络的超二次曲面单元能有效地构造非球形颗粒的几何形态,并通过非线性Newton迭代算法准确计算单元间的作用力.针对非球形颗粒间接触判断的复杂性及其大规模离散元计算的需求,该文发展了基于CUDA-GPU构架下超二次曲面单元并行算法.该方法在球形颗粒并行计算的基础上,通过核函数建立单元包围盒的粗判断列表及Newton迭代的细判断列表,并优化了并行算法和内存访问模式以提高算法的计算效率.为检验超二次曲面并行算法的可靠性,对非球形颗粒的流动过程进行离散元模拟,并与试验结果进行对比验证.在此基础上,进一步分析了颗粒单元不同长宽比和表面尖锐度对颗粒材料流动特性的影响,为非球形颗粒材料的大规模离散元模拟提供了一种有效的数值方法.