汽车碰撞过程的计算负载均衡并行有限元技术.docVIP

汽车碰撞过程的计算负载均衡并行有限元技术 　　摘要： 为了在有限的计算资源下挖掘更高的汽车碰撞仿真计算效率，基于当前普遍采用的并行处理区域分解算法，结合一个典型碰撞案例分析可知，汽车碰撞模型计算过程中的主要时间开销与结构大变形相关.提出基于坐标递归对分（Recursive Coordinate Bisection，RCB）的计算负载均衡区域分解算法，并将此法应用于某整车碰撞仿真模型的3类典型工况.计算结果表明此法能将并行计算效率提高10%～24%. 　　关键词： 　　汽车碰撞仿真； 并行计算； 区域分解； 计算负载均衡； 坐标递归对分； 计算效率 　　中图分类号： U462 　　文献标志码： B 　　0 引 言 　　在汽车被动安全研究领域中，计算机仿真是经济、有效的方法.汽车整车碰撞过程仿真是典型的大变形条件下的大规模冲击大变形计算，涉及高度的几何非线性、材料非线性和边界条件非线性等问题.为确保仿真计算结果有足够的精度，汽车整车有限元模型的单元数量一般以百万计，甚至达到千万级别.[1]得益于计算机软硬件的飞速发展，人们普遍采用大规模并行系统（Massively Parallel Processors，MPP）求解此类大规模动态非线性问题.因此，适合于MPP并行计算机的非线性有限元结构分析的并行求解技术引起高度重视，其中，基于区域分解算法的并行有限元法是当前并行计算最活跃的研究和应用领域之一.[29] 　　本文简要介绍基于区域分解算法的并行有限元法，分析并行计算区域分解算法的主要时间开销，提出基于坐标递归对分（Recursive Coordinate Bisection，RCB）的负载均衡区域分解算法[1011]，最后将此法应用于某整车碰撞仿真模型，探讨前部正面刚性墙碰撞、正面40%偏置碰撞以及可变形移动壁障侧面碰撞等3种工况下提高并行计算效率的可能性. 　　1 基于区域分解算法的并行计算技术 　　1.1 基于区域分解的显式并行有限元法 　　汽车碰撞过程的计算机仿真一般采用显式有限元求解，其基本流程[12]如下. 　　在MPP平台上实现以上计算流程.对于汽车碰撞这类大规模的仿真问题，目前最常采用的方法是区域分解算法. 　　区域分解算法是将一个复杂大系统，按照一定的规则，分解成与参加计算处理器的数目相等的子系统，各处理器承担各自子系统内的计算负荷，并通过交互机制（例如MPI方式）进行数据通信和数据交换.理论分析和实际计算[912]均已表明区域分解算法能够为大规模问题提供高度并行的、可扩展的鲁棒算法. 　　1.2 并行计算区域分解算法的主要时间开销 　　对于某个计算模型，采用并行计算时，其执行时间tP的计算公式[13]为 　　式中：tPa为并行部分的执行时间，通常为各个处理单元执行并行部分所耗时间的最大值；ti/o为输入/输出所耗费的时间，输入和输出部分通常都由单机完成[9]；tc为通信所耗费的时间，即通信开销. 　　加速比SP和并行效率η的计算公式分别定义为 　　式中：ts为单个处理单元计算该模型所耗费的时间；N为并行计算时所用的处理单元总数. 　　在采用区域分解的并行计算中，需要预先给每个处理单元分配相应的计算任务.分析式（6）和（7）可知，要提高加速比和并行效率，就要设法减小并行部分的执行时间tPa和通信开销tc， 　　式中：m为消息长度；n为处理器数目；t0（n）为启动时间，取决于具体的硬件系统；tc（n）为每字节通信时间，依赖于具体的计算平台. 　　因此，区域分解有以下3个原则[3，9]： 　　（1）各个子系统的计算负载大致平衡，以减少并行部分的执行时间tPa. 　　（2）最大相邻的处理单元数尽可能小，以减少通信次数. 　　（3）分子系统之间的边界节点数量尽可能小，以减少单次通信量. 　　区域分解的核心是静态负载均衡问题，系统的分解优劣关系到并行求解的效率.模拟过程中的动态负载（如接触力）均衡则更是不可或缺，?M量使各处理单元有平衡的负载，才能实现高加速比和高并行效率的并行计算. 　　RCB算法作为一种静态负载均衡算法[810]，可以使各个子系统的单元数量大致相等，其默认子系统分解见图1.由于该方法只考虑有限元模型的静态几何信息而未考虑模型在整个碰撞时域内实际载荷的动态分布，对于汽车碰撞这类大规模冲击大变形问题，默认子系统分解方式不能有效保证各处理单元计算负载均衡，所以并不是最理想的区域分解方案. 　　汽车碰撞模型计算过程中的时间开销分布情况见表1.整个模型时耗的90%以上都花费在单元处理和接触计算上，即第1.1节描述的第（2）和第（5）步.耗时多的计算都与结构大变形有关. 　　为此，本文探讨在保证各处理单元间单元数量平衡的基础上，通过均匀分布大变形的区域分解方案实