基于光滑聚集代数多重网格的有限元并行计算实现方法.docx

基于光滑聚集代数多重网格的有限元并行计算实现方法武立伟++张健飞++张倩摘要：基于光滑聚集代数多重网格法实现一种用于结构有限元并行计算的预条件共轭梯度求解方法。对计算区域进行均匀划分，将这些子区域分配给各个进程同时进行单元刚度矩阵的计算，并组合形成分布式存储的整体平衡方程。采用光滑聚集代数多重网格预条件共轭梯度法对整体平衡方程进行并行求解，在天河二号超级计算机上进行数值试验，分析代数多重网格的主要参数对算法性能的影响，测试程序的并行计算性能。试验结果表明该方法具有较好的并行性能和可扩展性，适合于大规模实际应用。关键词：有限元法；光滑聚集；代数多重网格；共轭梯度法；可扩展性中图分类号：TB121文献标志码：A收稿日期：2017[KG*9〗07[KG*9〗26修回日期：2017[KG*9〗09[KG*9〗04基金项目：中央高校基本科研业务费专项（2016B06414）；国家自然科学基金委员会广东联合基金超级计算科学应用研究专项作者简介：武立伟（1993—），男，山东淄博人，硕士研究生，研究方向为高性能计算，（Email）[e-mail]；张健飞（1977—），男，江苏海门人，副教授，研究方向为高性能计算与计算力学，（Email）[e-mail]elfiniteelementcomputationbasedonsmoothedaggregationalgebraicmultigridWULiwei，ZHANGJianfei，ZHANGQian（CollegeofMechanicsandMaterials，HohaiUniversity，Nanjing211100，China）Abstract：basedonthesmoothedaggregationalgebraicmultigrid，apreconditionedconjugategradientmethodforparallelstructuralfiniteelementcomputationisproposed.Thestructureisdividedintosubstructuresevenly，thesubstructuresareassignedtoprocessestocalculatetheelementstiffnessmatrixsimultaneously，andtheglobalequilibriumequationsinthedistributedstorageareassembled.Theparallelsolutionoftheglobalequilibriumequationsareimprovedbythesmoothedalgebraicmultigridpreconditionedconjugategradientmethod.ThenumericaltestiscarriedoutontheTianhe2supercomputer.Theinfluenceofthemainparametersofalgebraicmultigridonthealgorithmperformanceisanalyzed，andtheparallelcomputationperformanceoftheprocessisassessed.Thenumericalresultsshowthattheproposedimplementationisofgoodparallelperformanceandscalability，anditcanbeusedinlargescaleapplications.Keywords：finiteelementmethod；smoothedaggregation；algebraicmultigrid；conjugategradientmethod；scalability0引言有限元法是工程结构分析的一種重要数值方法。随着工程规模不断扩大、工程复杂性不断增加和计算精度要求不断提高，传统的串行有限元程序的计算规模和计算速度都已经不能满足需求，迫切需要发展能够在超级计算机上高效运行的可扩展有限元并行算法和程序。目前，传统有限元并行计算中常用的并行算法主要有子结构并行算法[1]、多波前法[2]和预条件共轭梯度法[34]等。这些方法中的直接法计算量和存储量大、并行程度不高，预条件共轭梯度法中的预条件为提高并行性而降低收敛性，也无法适应大规模应用。代数多重网格法[56]不需要几何网格信息，仅从方程组代数结构出发，形成虚拟的粗细网格，以达到加速收敛的目的。该方法具有存贮量小、收敛快和可扩展性好等优点，将其用于结构有限元大规模并行计算，可以进一步增大计算规模