多波前子结构并行有限元计算方法及应用.pdfVIP

第31卷增刊 防灾减灾工程学报 V01．31Suppl． 2011年9月 ofDisasterPreventionand Sep．2011 Journal MitigationEngineering 多波前子结构并行有限元计算方法及应用’ 张惠军，李俊杰 (大连理工大学建设工程学部，辽宁大连116024) 摘要：针对面板堆石坝结构规则、断面连续的特点，将总体结构划分为多个子结构，利用并行多波前处理技术实现 弹塑性有限元并行计算，将其应用在水布垭面板堆石坝施工加裁三雏有限元并行计算中，体现了施工分期填筑的 逐级加裁及堆石材料非线性变形特点。通过对计算结果分析和比较，表明该算法有效可行． 关键词：多波前法；并行计算；面板堆石坝；逐级加载 中图分类号：TU352．1文献标识码：A 文章编号：1672—2132(2011)增刊·0052—03 的优化策略，把大型稀疏矩阵的整个分解过程转化 引 言 为许多小的稠密矩阵的分解，从而使数值分解的运 算量减少到几乎最小。 在高面板堆石坝实际工程施工中，考虑到施工 并行多波前法(文献[23)E‘]在时间上不分先后， 期度汛与施工进度的需要以及堆石材料自身的特 采用的是一种“分而治之”的策略。基本思想是将总 性，一般采用分期填筑的方式，将坝体填筑和面板的 体结构按照某种剖分原则分解成若干个互不重叠的 浇筑分为不同的阶段。面板堆石坝有限元分析计算 子结构，利用多个波前独立并行地对每个子结构进 的逐级加载方法很大程度上贴近实际情况，能够很 行静凝聚。各子结构中波前过程集方程的组集和消 好地模拟坝体结构本身随施工过程的变化，反映出 元于一体，局部波前不断向前推移，直至达到子结构 材料刚度在施工蓄水分阶段变化过程的情况，更好 剖分时划定的界面，此时已完成了波前静凝聚过程。 地体现了堆石材料非线性变形的特点。然而，面板堆 在子结构凝聚过程中，每组装一次单元矩阵，对 石坝逐级加载三维有限元分析多级荷载步分析时间 已完全组集过的节点变量进行消元运算，消去第s 过长，精细化模拟受限于有限的内存空间和海量的 个变量对应的刚度矩阵和对应的右端荷载向量的运 数据处理，对其实现并行化处理具有重要的现实意 算过程如下： 义。 子结构法作为任务级并行的一种区域分解方 法，在有限元并行计算中，结合多波前法来凝聚子结 构对界面方程的贡献信息并求解，国内有多人做过 ∽H刖一∑譬 (2) 这方面的研究。胡宁等[1]在共享内存机上利用并行 式中 [K0]和(聍}为局部波前组装到第m个单元 多波前对结构两个子区凝聚。在网络机群环境下，余 时，这m个单元贡献的刚度矩阵和荷载向量的组集 天堂等[2]利用并行多波前法对2个子结构静凝聚， 值； 采用预条件共轭梯度法(PCG)并行求解界面方程； [K：]和[K0]为刚度矩阵中第s个变量对应 胡峥嵘等[3]将多波前方法用于面板堆石坝三维非线 的行和列； 性有限元并