等效裂隙岩体THM耦合并行有限元程序开发.pdfVIP

第26卷第10期 长 江 科 学 院 院 报 Vo1_26 No．10 2009年 10月 JournalofYangtzeRiverScientificResearchInstitute 0ct．20 0 9 文章编号：1001—5485(2009)10—0058—04 等效裂隙岩体 THM耦合并行有限元程序开发 张强林 。王 媛 ，曹国利 (1．黄河勘测规划设计有限公司，郑州 450003；2．河海大学 岩土所，南京 210098) 摘要：为了解决计算精度与计算时间之问的矛盾，针对三场耦合问题 ，根据区域分解策略，采用预处理共轭梯度法 进行方程组的并行求解，利用mpi+f77编制主从式并行有限元计算程序，并在河海大学高性能计算机群上进行调 试及计算，通过算例计算，显示并行计算可大幅节省计算时间，未来必会得到更广泛的应用。 关 键 词：并行有限元；THM耦合；区域分解；预条件共扼梯度法 中图分类号：TV641；TB115 文献标识码：A 有限元数值计算是研究岩土工程问题的一种重 区域分解策略，另一种则是方程组解法并行策略。 要方法，由于岩土工程领域研究对象所处的环境越 目前应用较多的则是 2种策略方法相结合，即先进 来越复杂，模型考虑越来越多，再加上原本的岩土介 行区域分解，而后以各种预条件子的共轭梯度为解 质材料非线性、各向异性等计算需要，使得数值计算 法求解。本文开发的三场耦合并行求解程序也是如 工作量越来越大，而计算时间也越来越长，如果采用 此 。 简化模型或粗网格等手段来减少计算量则计算精度 1．1 区域分解算法 会受到影响。计算精度与计算时间之间的矛盾，对 区域分解算法 的基本思想是通过将一个区 于现有的有限元模拟来说是很大的挑战，因此有限 域划分为若干个独立的、规模较小的子区域，对一个 元并行计算模拟即成为人们解决此 问题的首选方 或若干个子问题用一个处理机来处理，使原问题的 案，如张友 良等 J、茹忠亮等 针对水布垭地下厂 求解转化为各子区域上子问题的求解，通过并行求 房的研究，刘耀儒等 ．4针对二滩、锦屏等工程研究 解子问题而获得整个区域上大问题的解。区域分解 都是针对单场进行并行计算，Wang等 、江春波 算法通过对子问题的同步并行求解来缩小问题规 等 也对渗流问题的并行计算进行了研究。 模，降低求解复杂度，提高计算效率，体现的是 “分 在考虑多场耦合情况下，由于研究问题更复杂， 而治之”的思想。相对于传统算法，区域分解算法 考虑因素更多，多场耦合并行计算尚处于起步阶段， 的优越性体现在 以下几个方面： 耦合问题的并行计算并不多见，而对于 目前研究的 (1)它把大问题化为若干个小问题，缩小了计 热点之一——三场耦合 (温度一渗流一应力耦合)问 算规模 ； 题更是如此 ，纵观国内外针对三场问题的研究，仅仅 (2)子区域形状如果规则，可以使用熟知的快 有 Yu-Shuwu等 ¨ 利用TOUGH2并行版，对 Yuc- 速算法，如快速Fourier交换、谱方法等； ca山中的三场试验进行了研究，并取得了不错的效 (3)允许使用局部拟一致网格，无需使用整体 果，但其他的研究成果 目前 尚未发现。由此本文即 拟一致网格 ； 针对三场耦合问题，根据区域分解策略，开发三场耦 (4)允许不同子区域选用不同的数学模型，以 便更适合于工程物理实际情况； 合并行程序进行数值模拟，并在河海大学高性能计 算机群上利用算例验证了程序的可靠性 。