武立偉++張健飛++張倩

摘要： 基于光滑聚集代數(shù)多重網(wǎng)格法實現(xiàn)一種用于結(jié)構(gòu)有限元并行計算的預(yù)條件共軛梯度求解方法。對計算區(qū)域進行均勻劃分，將這些子區(qū)域分配給各個進程同時進行單元剛度矩陣的計算，并組合形成分布式存儲的整體平衡方程。采用光滑聚集代數(shù)多重網(wǎng)格預(yù)條件共軛梯度法對整體平衡方程進行并行求解，在天河二號超級計算機上進行數(shù)值試驗，分析代數(shù)多重網(wǎng)格的主要參數(shù)對算法性能的影響，測試程序的并行計算性能。試驗結(jié)果表明該方法具有較好的并行性能和可擴展性，適合于大規(guī)模實際應(yīng)用。

關(guān)鍵詞： 有限元法； 光滑聚集； 代數(shù)多重網(wǎng)格； 共軛梯度法； 可擴展性

中圖分類號： TB121文獻標(biāo)志碼： A

收稿日期： 2017[KG*9〗07[KG*9〗26修回日期： 2017[KG*9〗09[KG*9〗04

基金項目： 中央高?；究蒲袠I(yè)務(wù)費專項（2016B06414）；國家自然科學(xué)基金委員會廣東聯(lián)合基金超級計算科學(xué)應(yīng)用研究專項

作者簡介： 武立偉（1993—），男，山東淄博人，碩士研究生，研究方向為高性能計算，（Email）liweiwu@hhu.edu.cn；

張健飛（1977—），男，江蘇海門人，副教授，研究方向為高性能計算與計算力學(xué)，（Email）jianfei@hhu.edu.cn

Implementation of parallel finite element computation

based on smoothed aggregation algebraic multigrid

WU Liwei，ZHANG Jianfei，ZHANG Qian

（College of Mechanics and Materials， Hohai University， Nanjing 211100， China）

Abstract： Based on the smoothed aggregation algebraic multigrid， a preconditioned conjugate gradient method for parallel structural finite element computation is proposed. The structure is divided into substructures evenly， the substructures are assigned to processes to calculate the element stiffness matrix simultaneously， and the global equilibrium equations in the distributed storage are assembled. The parallel solution of the global equilibrium equations are improved by the smoothed algebraic multigrid preconditioned conjugate gradient method. The numerical test is carried out on the Tianhe 2 super computer. The influence of the main parameters of algebraic multigrid on the algorithm performance is analyzed， and the parallel computation performance of the process is assessed. The numerical results show that the proposed implementation is of good parallel performance and scalability， and it can be used in largescale applications.

Key words： finite element method； smoothed aggregation； algebraic multigrid； conjugate gradient method； scalability

0引言

有限元法是工程結(jié)構(gòu)分析的一種重要數(shù)值方法。隨著工程規(guī)模不斷擴大、工程復(fù)雜性不斷增加和計算精度要求不斷提高，傳統(tǒng)的串行有限元程序的計算規(guī)模和計算速度都已經(jīng)不能滿足需求，迫切需要發(fā)展能夠在超級計算機上高效運行的可擴展有限元并行算法和程序。目前，傳統(tǒng)有限元并行計算中常用的并行算法主要有子結(jié)構(gòu)并行算法[1]、多波前法[2]和預(yù)條件共軛梯度法[34]等。這些方法中的直接法計算量和存儲量大、并行程度不高，預(yù)條件共軛梯度法中的預(yù)條件為提高并行性而降低收斂性，也無法適應(yīng)大規(guī)模應(yīng)用。代數(shù)多重網(wǎng)格法[56]不需要幾何網(wǎng)格信息，僅從方程組代數(shù)結(jié)構(gòu)出發(fā)，形成虛擬的粗細網(wǎng)格，以達到加速收斂的目的。該方法具有存貯量小、收斂快和可擴展性好等優(yōu)點，將其用于結(jié)構(gòu)有限元大規(guī)模并行計算，可以進一步增大計算規(guī)模、提高計算速度，能滿足現(xiàn)代工程結(jié)構(gòu)分析與設(shè)計的要求。

本文基于光滑聚集型代數(shù)多重網(wǎng)格法[7]，實現(xiàn)一種用于結(jié)構(gòu)有限元大規(guī)模并行計算的預(yù)條件共軛梯度方法，并在天河二號超級計算機上對不同的聚集策略、光滑迭代算法和粗網(wǎng)格求解方法進行影響分析，測試和分析程序的并行性能。

1有限元計算

有限元法是通過連續(xù)體的離散與分片插值，將求解物理問題的控制微分方程轉(zhuǎn)變?yōu)榍蠼饩€性代數(shù)方程組的一種近似數(shù)值算法。首先，將問題的求解區(qū)域剖分成有限個單元的集合，在單元內(nèi)用分片插值表示待定函數(shù)的分布；然后，由變分原理獲得定解問題的線性代數(shù)方程組；最后，求解該方程組即可得到待定函數(shù)的數(shù)值解。endprint