基于分块数据结构冲击问题并行计算.docVIP

基于分块数据结构冲击问题并行计算 　　摘 要：针对三维冲击问题，基于分块数据结构在共享内存并行机上实现OpenMP并行计算.分块数据结构不仅能有效利用计算机多层存储结构，而且增加OpenMP的并行粒度.数值实验表明：在使用分块数据结构后，串行程序的计算速度能提高3倍.通过柱体冲击平板数值模拟实验讨论并行程序的加速比和效率，表明并行程序能有效减少总计算时间. 　　关键词：OpenMP； 冲击； 并行计算； 分块数据结构 　　中图分类号：TP311 文献标志码：??A?? 　　 　　Parallel computation of impact problems based onblock data structure 　　XIAO Yonghao， HUANG Qingnan 　　(Computer Application Institute, China Academy of Engineering Physics, Mianyang 621900, Sichuan, China) 　　Abstract： As to 3D impact problem, OpenMP parallel computation is implemented in shared memory parallel computer based on block data structure. Block data structure is not only effective use of computer multi-storage structure, but also improves the OpenMP parallel granularity. The numerical experiment indicates that the serial program with block data structure is three times faster than original version. The speedup and efficiency of parallel program are discussed by anumerical simulation on impact of cylinder into a plate and the result shows that the parallel program can decrease the total computation time effectively. 　　Key words： OpenMP; impact; parallel computing; block data structure 　　 　　0 引 言 　　冲击问题大量存在于工程应用和现实生活中，如汽车碰撞和物体跌落等.冲击问题属非线性大变形问题，可用LS_DYNA和Abaqus等商业数值模拟软件对此类问题进行数值模拟计算.随着计算模型尺寸和复杂度的提高以及数值模拟精确性要求的增加，现有的单处理器性能已不能满足模拟需求，需进行冲击问题的并行算法和并行程序实现研究. 　　从计算机的体系结构来说，在21世纪初，单处理器的主频已达到顶峰，各处理器生产厂商开始生产多核处理器.2005年底，SUN公司发布的UltraSPARC T1处理器包含8个内核，每个内核能同时运行4个线程.IBM的Cyclops-64处理器在每块芯片上集成80个内核，能同时运行160个线程.未来的高端计算将支持多线程与消息传递相结合的计算模式.［1］ 　　OpenMP［2］作为科学计算领域中共享内存并行机主要并行化范式将成为主流，运行在共享内存节点上的OpenMP程序相对于MPI程序，在可扩展性、网络性能、内存需求和编程效率等各方面都具有其优点［1］. 　　1 冲击算法 　　大变形条件下冲击过程的数值求解是将空间进行有限元离散，某时刻的有限元方程以矩阵形式表示为?И?Mu¨+f???┆?int????=f???┆?ext?????И? 式中：??M=∫??V?[WTHX]N????T??N??d??V??为质量阵，N为形函数；u¨为加速度向量；??f???┆?int????=∫??VB????T???[WTBX]??d??V??为内力；f???┆?ext????[WT]为外力. 　　冲击计算可分为：内力计算，主要用于更新单元应力、应变，最终计算节点力；接触力计算，主要搜索接触对，进行接触计算. 　　在冲击计算中，接触通过边界条件施加到总体方程中.接触力计算分为接触面搜索、全局搜索、局部搜索和接触力计算等4个步骤.［3-4］由于接触关系的时变性质，判断接触面（接触面搜索）和搜索接触关系（全局搜索）是个全