非连续变形分析（DDA）线性方程组的高效求解算法 - 岩土力学.PDF

第37 卷第4 期 岩 土 力 学 Vol.37 No. 4 2016 年4 月 Rock and Soil Mechanics Apr. 2016 DOI ：10.16285/j.rsm.2016.04.032 非连续变形分析 （DDA ）线性方程组的 高效求解算法 1 1 1 1, 2 付晓东 ，盛 谦 ，张勇慧 ，冷先伦 （1 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室，湖北 武汉 430071 ； 2. 中南勘测设计研究院 水能资源利用关键技术湖南省重点实验室，湖南 长沙 410000 ） 摘 要：非连续变形分析(DDA)方法对大规模工程问题的数值模拟耗时太长，其中线性方程组求解耗时可占总计算时间的 70% 以上，因此，高效的线性方程组解法是重要研究课题。首先，阐述了适用于DDA 方法的基于块的行压缩法和基于试验- 误差迭代格式的非0 位置记录；然后，针对DDA 的子矩阵技术，将块雅可比迭代法 (BJ)、预处理的块共轭梯度法 (PCG， 包括Jacobi-PCG、SSOR-PCG) 引入DDA 方法，重点研究了线性方程组求解过程中的关键运算；最后，通过两个洞室开挖 算例，分析了各线性方程组求解算法在DDA 中的计算效率。研究表明：与迭代法相比，直解法无法满足大规模工程计算需 要；BJ 迭代法与块超松弛迭代法(BSOR)的效率差别不大，但明显不如PCG 迭代法。因此，建议采用PCG 迭代法求解DDA 线性方程组，特别是 SSOR-PCG 值得推广；如果开展并行计算研究，Jacobi-PCG 是较好的选择，当刚度矩阵惯性优势明显 时，BJ 迭代法同样有效。 关 键 词：岩土工程；非连续变形分析；子矩阵技术；刚度矩阵存储；线性方程组求解；块雅可比迭代与预处理的块共轭 梯度法 中图分类号：TU 452；O 241 文献标识码：A 文章编号：1000－7598 (2016) 04 －1171－08 High efficient algorithms for solving linear equations in discontinuous deformation analysis