MSC.Software2010中国用户大会论文*本文主体曾发表于我单位内部刊物《高性能计算发展与应用》2010年第2期,此次投稿略加改动.作者简介:王惠,工学硕士,HPC应用技术工程师,E-mail:hwang@ssc.net.cn基于并行计算的钢板加强钢筋混凝土三维侵彻仿真研究*王惠,丁峻宏,李根国(上海超级计算中心,上海201203)摘要:钢筋混凝土抗侵彻性能是建筑防护研究的重要课题。本文建立钢板加强钢筋混凝土有限元模型,研究了不同钢板厚度下钢筋混凝土的抗侵彻能力,通过对比说明在钢筋混凝土表面覆以钢板并通过增加钢板厚度的方式能有效提高钢板混凝土防护能力。计算中结合上海超级计算中心高性能计算平台,对计算模型采用基于区域分解法的并行计算方法,探讨了采用不同处理器个数对加速比和并行效率的影响,对如何提高计算效率提出了合理建议。关键字:钢筋混凝土钢板并行计算侵彻Three-DimensionalPenetrationSimulationResearchforSteelPlateReinforcedConcreteBasedonParallelComputingWangHuiDingJunhongLiGenguo(ShanghaiSupercomputerCenter,Shanghai201203,China)Abstract:Itisofmuchimportancetoinvestigatetheanti-penetrationperformanceofsteelplatereinforcedconcretetoprotectsomespecifiedarchitectures.Thispaperestablishedacorrespondingfiniteelementmodel,andevaluatedtheanti-penetrationcapability.Bycomparison,itshowsthattheprotectioncapabilitycouldbeeffectivelystrengthenedbyaddingsteelplateontheoutsideofreinforcedconcrete,andincreasingthethicknessofsteelplate.CalculationwascompletedbasedontheHPCplatformofSSC,andthedomaindecompositionmethodwasutilizedforparallelcomputing.TheinfluenceofCPUnumbersonscalabilityandcomputingefficiencywasstudied,andsomeproperadviceswerefinallyarrivedat.Keywords:ReinforcedConcrete;SteelPlate;ParallelComputing;Penetration1研究背景钢筋混凝土是目前工业与民用建筑中应用最广泛的结构工程材料之一,已被广泛应用于构建高层建筑物、长跨桥、大坝水电站、隧道、码头等。此外,在常规战争中,许多战略价值高的目标(如核设施、控制通讯掩体、民防掩蔽工程、桥墩、飞机库及机场跑道等)大都是由钢筋混凝土构建而成。随着现代高新技术武器的发展,精确制导武器命中率极高,对重点防护工程构成了严重的威胁,当这些建筑物被作为攻击的目标时.如何更行之有效地发挥防护工程的防护功能广受科研工作者关注。研究强冲击载荷下混凝土的侵彻特性,对于提高防护工程的防护能力以及探索研究新型防护结构具有重要意义。防护结构主要是钢筋混凝土材料构筑,通过增加结构层的厚度和提高混凝土材料强度成为提高防护能力的首要手段,但研究表明只靠增加混凝土防护层厚度并不能成为提高防护能力的首要措施,而在钢筋混凝土的表面覆以钢板可以在几乎同样厚度的条件下大幅度提高结构的防护能力,这种措施已经在土建结构中广泛应用并取得良好的效果,在国外核反应堆安全壳以及导弹发射井中均有应用。目前对钢筋混凝土的防护能力研究主要集中在钢筋混凝土的抗冲击、抗爆炸动力学特性研究,其中,用弹体侵彻钢筋混凝土来研究钢筋混凝土的防护能力一直是研究的热点问题。弹体对钢筋混凝土靶板的侵彻问题是一个过程十分短暂的高压力、大应变的非线性动态问MSC.Software2010中国用户大会论文2题。数值模拟可以非常逼真地再现侵彻全过程的所有细节,给出与侵彻相关的所有物理量,如应力场、变形场、微裂纹的演化过程等。目前,能够数值模拟强冲击载荷下混凝土侵彻特性的分析软件很多,在众多的分析软件中,LS-DYNA多年来一直是最具通用性的非线性动力分析软件之一,尤其在侵彻过程与爆炸成坑模拟分析、军用设备和结构设施受碰撞和爆炸冲击加载的结构动力分析等领域被广泛认可。MSC公司的Dytran求解器包含了LS-DYNA求解器,用户可在MSC.Dytran界面直接调用LS-DYNA求解器,LS-DYNA求解功能和并行计算功能均可方便实现,其优秀的分布式并行计算技术可以十分方便地将一个规模较大的题目分阶段进行运算。本文基于LS-DYNA软件模拟了强冲击载荷下混凝土侵彻特性,探讨了不同钢板厚度下钢筋混凝土的抗侵彻性能,通过对比说明,在钢筋混凝土表面覆以钢板并通过增加钢板厚度的方式是一种提高钢板混凝土防护能力行之有效的方法。强冲击载荷下混凝土侵彻是钢筋混凝土在极短的时间内在剧烈碰撞动态载荷作用下发生的复杂非线性动态响应过程。在侵彻过程中,钢筋混凝土结构经历复杂的变形,具有材料非线性、几何非线性和接触非线性的特点,因而使得计算过程复杂,计算量大,随着实际工程问题的规模日益增大,对求解精度也提出了更高的要求,因此需要大量的内存、硬盘资源和快速的求解能力。传统有限元软件均在单CPU上运行,这就为精确的分析带来了种种困难。随着计算机软硬件的飞速发展,商品化有限元软件结合并行求解技术,开始利用多个CPU并行处理求解大规模动态非线性复杂问题。本文基于区域分解的并行计算方法,探讨采用不同数目处理器对钢筋混凝土侵彻仿真计算效率的影响。2钢筋混凝土侵彻仿真的并行原理和方法2.1基于区域分解的并行有限元法钢筋混凝土侵彻过程是有着内部接触高度非线性的动力学问题,隐式求解方法往往无法保证收敛。本文采用显式中心差分方法来进行时间积分。对整个钢筋混凝土侵彻模型进行有限元单元划分后,引入应变位移关系、应力应变关系,t时刻系统动力学方程可描述为:cttttFFQM(1)其中:M是质量矩阵,t是加速度,tQ是外载荷,tF是内力,ctF是接触力。对于式(1),采用中心差分法,其关系表达式为:tttttt)2/()2/((2))2/(tttttt(3)联立(2)和(3)式,可以得到:2/)2(ttttttt(4)根据(1)和(4),可建立时间递推公式:]2)([21tttctttttMMFFQtM(5)其中t是时间步长,tt,t,tt为不同时刻的位移。采用集总的方法可以将质量矩阵M对角化,式(5)解耦成5n个方程(n为有限元模型节点总数),这些方程式不需要集合总体刚度矩阵,并采用中心单点积分,从而可以大大节省存储空间和求解机时。显示有限元并MSC.Software2010中国用户大会论文3行化就是对式(5)的并行计算。目前最常采用的方法是区域分解法。该方法就是把一个复杂系统或区域,按照一定的原则(如物理特性、几何形状、离散方式等)分解为若干个子系统或子区域,然后将各个子区域分配给不同的处理器分别并行计算,处理器之间通过交互机制(例如MPI方式)进行数据交换。区域分解的作用不仅在于将复杂的系统分解为简单的子问题,而且非常有利于并行计算,各个子域可以在各处理机上单独运行,然后将各个子域归并就可以得到问题的完整解。本文对钢筋混凝土侵彻仿真研究就是基于区域分解的并行计算方法。其基本过程:首先建立整个钢筋混凝土侵彻行为有限元模型、定义初始条件、边界条件、材料模型、接触条件和分析类型、求解控制、结果输出内容和方式,然后通过程序提供的分区方法将整个有限元模型分成若干子区域并将各个子区域分配给不同的处理器分别计算。在每一个时间步长内,各个处理器独立计算对应子区域中单元的位移、速度、加速度和应力等,相关处理器之间的通信仅需传送子区域边界单元在前一时间步长的位移计算结果,最后将每个子域的解进行合并即得到整个计算域解。2.2区域分解法的分区策略和原则区域分解是并行有限元分析的一个很重要的前处理部分,考虑到显式有限元绝大多数计算是针对于每个单元进行的,因此可以方便地分割单元区域。其基本思路如图1所示,当数据文件读入完毕,采用区域分割方法得到各个子区域的分区定义,并且为每一个子区域指定处理器负责计算,保证各个子区域在不同的处理器并行计算。由于有限元模型往往是大规模的、不规则的,区域分解的好坏是一个载荷平衡的问题,直接关系到并行求解的效率,所以一种效果好的区域分解算法是很重要的,在分区过程中应尽量使各CPU处理相同的计算和通讯量,才能实现高效率和高加速比的并行计算。区域分解的判断标准可以归纳为3点:(1)各子域的自由度数相同;(2)各子区域交界面上的节点数最小;(3)所分割的子区域有较好的长宽比。分区最重要的原则是分区包含的单元数应接近,以保证各个分区内的计算量大致相同。图1显式有限元并行计算流程图Fig.1FlowchartofexplicitFEMparallelcomputing有限元计算分区方法对并行计算的效率高低有着决定性的影响,目前主要的区域分解算法有模拟退火法、贪婪方法、递归谱二分法、递归坐标二分法(RecursiveCoordinateBisection,MSC.Software2010中国用户大会论文4简称RCB)等。本文钢筋混凝土侵彻有限元计算采用是非线性动力分析软件LS-DYNA,其并行版本默认的区域分割方式是坐标递归对分方法。该方法属于几何模型分区方式,该方法的基本原理如下:1)判断模型在三个坐标方向的长度,垂直于最长坐标方向将模型一个为二,区域边界的节点被复制,分布到对应的区域中;2)判断新模型在三个坐标方向的长度,仍然垂直于最长坐标方向将新模型一个为二,并复制边界节点,将其分布到对应区域;3)按照上面的策略,依次将模型分割直到满足分割收敛条件,算法就停止分割;4)输出模型分区结果。3典型应用实例本文钢筋混凝土侵彻数值计算是在上海超级计算中心的“曙光5000A”高性能计算机(又名魔方)上完成的。3.1钢板加强钢筋混凝土侵彻数值模型钢板加强钢筋混凝土主要由混凝土、钢筋、钢板、螺钉等组成,见图2。按照实际尺寸,分别建立混凝土、钢筋、钢板、螺钉的几何模型,其主要尺寸如下:混凝土长宽高分别为1800*1800*120mm;前后钢板厚度为2.4mm,钢筋直径为10mm,间距为85mm,钢筋距前钢板约为10mm,螺钉直径为20mm,长度为20mm,螺栓间距为40mm,螺钉用于把钢板固定在混凝土表面,有限元模型见图3。有的文献中利用分离式Hopkinson压杆(SHPB)对高刚度混凝土进行了动态冲击试验,本文综合相关资料后设计选用的冲击杆为圆柱体形状,长度为0.5m,杆直径为100mm,取名为SRC-2.4。图2钢板加强钢筋混凝土模型Fig.2CADmodelofsteelplatereinforcedconcreteofcertaintype冲击杆侵彻钢板加强钢筋混凝土有限元模型见图,整个计算模型单元数超过20万,模型详细情况见表1。计算中,混凝土材料模型采用MAT_CSCM_CONCRETE模型,其密度为24503/mKg,抗压强度为38.1MPa,失效控制系数为1.05。钢板采用45号钢,假设为线弹性材料模型,密度为78503/mKg,杨氏模量为206GPa,泊松比为0.3,屈服强度为356MPa,失效应变为0.26。冲