基于CUDA的并行计算在瞬态分析中的应用.pdfVIP

江　苏　航　空 2010增刊　 基 于 CUDA的 并 行 计 算 在 瞬 态 分 析 中 的 应 用 李　杨　陈国平　邓　军 (南京航空航天大学航空宇航学院 , 南京 , 2 10016) 摘要: 针对 瞬态分 析耗时长的问题 ,研 究了基 于 CU DA 的并行 计算方 法 。 CU DA是 NV IDI A 公司推 出 的基 于 GPU 的并行编程模型 ,采用 C语言编程 。它 利用 GP U 多处理器的特点来实现大规模数据 的并行计算 , 从而显著 提高 问题的求解速度 。本文将 CU DA 技术应用到瞬态分析 问题的显式求解上 ,较大的提高了求解效率 。 关键词: 中心差分法 ; CU DA; 并行计算 引　　言 1　瞬态分析问题的显式积分过程 瞬态分析 问题主要研究结构在外载荷作用下 系统的运动微分方程 ¨ 的时间历程响应 ,计算结构在振动时的动位移和动 ext Mut + Cut + Kut = F t ( 1) 应力的大小以及变化规律 ,是结构动力学的一个重 其 中 、、 分别为系统的质量矩阵、阻尼矩阵和 M C K ¨ · 要研究方向。求解瞬态分析 问题通常有显式积分和 刚度矩阵 ; t、 t、t 是系统的加速度 向量、速度 向 u u u 隐式积分两种求解格式 。由于求解此类问题 的时间 量 、位移向量 ;F t 是外载荷 向量 。 步长一般较小 ,实现时间序列上 的求解需要耗费较 不考虑系统的阻尼 ,则上式可以改写为 ¨ 多的时间 ,因此瞬态分析问题的并行算法研究得到 ext int Mut = F t - F t ( 2) [1 ] 了很多学者 的重视。自 N ASA的 N oor A K 第一 int 其 中F t 是系统的节点内力向量 次在向量处理机上实现 了显式积分 的并行处理 以 int F t = Kut ( 3) 来 , 国内外大量学者在瞬态分析 的并行计算方面开 　　 中心差分法是求解瞬态 问题 的一种有效的方 展 了卓有成效的研究。 法 。由于在每一步迭代运算时 ,