非线性动力有限元重叠区域分裂的隐式并行算法.pdfVIP

第23卷第6期 　 计 算 力 学 学 报 　 Vol . 2 3 ,No . 6 2006 年 12 月 Chinese Journal of Computational Mechanics December 2006 文章编号 (2006) 非线性动力有限元重叠区域分裂的隐式并行算法 付朝江 1 ,2 , 　张 　武1 ( 1. 上海大学 计算机工程与科学学院 , 上海 200072 ;2 . 南昌大学 建筑工程学院 , 江西 南昌330029) 摘 　要 :针对大规模结构非线性瞬态动力分析非常耗时 ,提出了相应的并行算法 。该算法采用无条件稳定的 Ne β ( ) wmark 方法 平均加速技术 进行时间积分 ,并结合区域分裂技术进行分析 。它不同于已有的采用非重叠区域 的并行算法 ,而是采用重叠区域的并行算法 。对给定结构有限元分析的质量 、阻尼 、刚度矩阵进行分裂可推出重 叠区域分裂算法的计算公式 。为改善每一步的求解 ,采用预估和校正子方案 。编写了该算法的程序 ,在工作站机 群上实现了数值算例 ,验证了算法的性能 。计算结果表明该算法优于非重叠区域分裂算法 。 关键词 :并行计算 ;非线性动力 ;有限元方法 ;Newmark 算法 ;非重叠区域分裂 。 中图分类号 :O242 ; O354 　　　文献标识码 :A [ 5 ] [ 6 ] 性 。此外 ,也提出 一些专门方法如组隐式算法 1 　引 言 ( ) GI 和隐式处理单个子区域和显式处理界面的并 区域分裂法是将计算区域分割成若干独立的 行混合时间积分算法 。为避免对通信开销最大的 规模较小的子区域 ,使原问题的求解转化为各子区 界面进行并行求解 ,开发了这种 GI 算法和混合算 域上子问题的求解 ,一个或若干个子问题用一个处 法 。然而 , 由于实际的时间步大小不定 , 这些 GI 理机来处理 。当所有的子问题均衡地映射到多个 和混合时间积分算法的数值精度受到影响 。至今 处理机上时 ,通过对子问题并行求解 ,提高解题速 所见到的文献中 ,动力有限元分析大都采用非重叠 度 。根据子区域是否重叠 ,区域分裂算法分为重叠 区域分裂法[7 ] 。 区域分裂算法和非重叠区域分裂算法 。对于重叠 有限元法的隐式公式提供一个可行的方法来 区域分裂 ,重叠区域中的变量与多个子区域相联 求解承受动力和准静力结构的时程反应 。大规 系 ;而对于非重叠区域分裂 , 只有区域之间的界面 模结构动力 问题 ,往往 同时包含非线性 , 如材料 上的变量与多个子区域相联系 。重叠区域分裂迭 非线性 、几何非线性特性和接触非线性 , 用这种 代速度一般更快 ,但子区域 由于与其他区域有重 隐式动力有限元方法计算大规模结构非线性动 叠 ,因而相对较大 。 力反应的计算复杂度是非常高 ,使得计算的难度 近 10 多年来 ,提出了用区域分裂方法来并行 和时间显著增加 ,并行计算可以实现大规模结