基于重叠区域分解思想的有限元程序并行化研究.pdfVIP

基于重叠区域分解思想的有限元程序并行化研究 张磊张国新 中国水利水电科学研究院结构材料所，100038 *i*n目}自青mt#序*”女m一*女￡$m”t±■∞Ⅱ域》*i*(f#日 FET[EBE)，toi**女m普《≈m较*堆L作量≈{*8”目Ⅱ单一是*№ ￡．奉t*女一竹基十AZTEC十parMctls库∞|叠＆B№e*}*，*艘日戒I井行dB 自^n有mt*有步#*女n*#m．#有目#展H＆章*i式是镕高《等特{，井 Ⅱ§t女战，#枷*#*#％#∞I#t 芙链目：*”有限iE域e*AZTEC并n目M》 引言 有限元计算主要有前处理阶段、单元分析、单元矩阵和向量装配成总体矩阵和向量、 约束处理、方程组求解、后处理选六部分，实际上单元分析、组装和方程组求解是有限 元分析的关键所在，特别是方程组掌解．一般至少能占到整个计算量的70％．但是在问 题的规模非常大的情况下．或者受存储的限制，前后处理的计算时间也不窖忽视，比如 流体力学计算等每个计算步骤要变化罔格的情况下．前处理耗时也很大，因此对于大规 模复杂结构的有限元并行计算．星好能对H上六部分同时进行井行实现，传统而容易实 现的并{T肇略是将求解方程组部分并行计算．但这种方式可扩展性差．有时候还需要用 创外存。区域分解方法具有井行度高、可扩展性好的优点．已经成为当前并行计算巾最 活跃的研究与应用{噩域之，也垦有阻元并行计算的最有趣的方法之一．比如子结构法、 FETI法及其它子结构法变种咀及并行EBE方法等，这些并行空现JJ法普遍都比较繁琐， 特别是用这些方法将串行有限元程序改成并行，工作量很大，且求解器部分针对性单一． 灵活性差。奉文提出一种基十区域分解思想的全局有限元并行实现策略。 二、并行实现思想 2：；：娑i嚣 子匡域．各个灶理帆保存每个子区域的信息(单元、节点，约柬)，以平面三角形单元，。嚣亍；ii要i：篓嚣翠；：蒸糯聋墨警鬻躲j 时糖…烈州皑g] m“’ 似卜惦罐】P嘛卜瞄荨j 如果将网格进行不重叠分区．对于各于医域边界 点形成总刚时，难免有些[捌元素在其他的于区 域，固此每个子区域计算时需蔓在边界上进行通 讯．这样必彝!；会增加程序编写的复杂性。而采用 重叠的分区方法，也就是在每个子区域的边界再 加一层单元，多计算一层单元矩阵，子区域总刚 组集时避免通信，而疗程组求解时通信量最小。 港 这样仅仅增加一点存储量和计算晕，一个测试”I 表明时于一十810000个单元的六面体单元嗣格， 增加的实际计算量的单元数也就5000十左右．增 目1重盎型E*ⅫH 加点¨尊量来减少大规模的通信．并Ⅱ减少程序编写的复杂度，显然是报划掉的。 这种思想的并行实现主要肓两个模块重叠l蔓域分解模块和并行求解模块。区域分 解模块土要是两方面内容：并行阿格生成和并行分区。并行阿格生成目前还没有好的商 业软件和开源软件．并行嗣格生成目前也是研究的鼎点和难点，对于水工结构计算通常 不用变换阿格(静态图)，孵此没必要并行生成网格．井行分区目前主要是基于并行匣 剖分方法．许多开源库可以实现，比如J∞tle、PⅡMeth等．本文采用PⅡMetis。如图 2l所示．采Ⅲ基于节点划分”啪区域分解法．区域划分后相当于把总体矩阵按行分配到 每个处理机．每个子区域内的单元分析和总体矩阵的形成都是完全井行的：在对整个结 构分区豹过程中，同样可以对边界进行分区旌加，因此形成右端向量以及约束处理也可 以完哞=并行虫现。这样，系数矩阵和右端项都分布存储在各个处理机。这种数据的独立 性．还可通过井行I／0．进一步提高程序性能。井行求解横块通过对AZTEC库进行并 行编程宴现。AZTEC提供了多种迭代法和顶条件子．将其为我所用增加了程序解法的 灵活性。由于数据分布存储，后处理阶段也可以并行实现。 三、程序实现侧试 基于上述有限元并行化思想，作者开发了区域分解 模块和求解罂模块，并且肘一个三维弹性力学问题有限 元半行程序进行了并行化。区域分解模块基于ParMetis