›› 2016, Vol. 37 ›› Issue (4): 1171-1178.doi: 10.16285/j.rsm.2016.04.032
付晓东1,盛 谦1,张勇慧1,冷先伦1, 2
FU Xiao-dong1,SHENG Qian1,ZHANG Yong-hui1,LENG Xian-lun1, 2
摘要: 非连续变形分析(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
[1] | 何静斌, 冯忠居, 董芸秀, 胡海波, 刘 闯, 郭穗柱, 张聪, 武敏, 王振, . 强震区桩−土−断层耦合作用下桩基动力响应[J]. 岩土力学, 2020, 41(7): 2389-2400. |
[2] | 师旭超, 孙运德. 线性卸荷作用下软土超孔隙水压力 变化规律分析[J]. 岩土力学, 2020, 41(4): 1333-1338. |
[3] | 田威, 王震, 张丽, 余宸, . 高温作用后3D打印岩体试样力学性能初探[J]. 岩土力学, 2020, 41(3): 961-969. |
[4] | 徐栋栋, 邬爱清, 李 聪, 汪 斌, 蒋昱州, 曾 平, 杨永涛, . 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学, 2019, 40(3): 1169-1178. |
[5] | 赵强, 焦玉勇, 张秀丽, 谢壁婷, 王龙, 黄刚海, . 基于显式时间积分的球颗粒DDA计算方法[J]. 岩土力学, 2019, 40(11): 4515-4522. |
[6] | 范文亮,王余乐,魏奇科,杨朋超,李正良, . 岩土工程可靠度分析的改进四阶矩方法[J]. , 2018, 39(4): 1463-1468. |
[7] | 黄明华,赵明华,陈昌富. 锚固长度对锚杆受力影响分析及其临界值计算[J]. , 2018, 39(11): 4033-4041. |
[8] | 熊自明,卢 浩,王明洋,钱七虎,戎晓力,. 我国大型岩土工程施工安全风险管理研究进展[J]. , 2018, 39(10): 3703-3716. |
[9] | 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171. |
[10] | 方砚兵,苏永华,肖 旺,梁 斌. 基于子区间法的隐式功能函数非概率可靠性方法研究[J]. , 2017, 38(4): 1171-1178. |
[11] | 黄小福,张迎宾,赵兴权,余鹏程,邢 昊,张 珏,陈岩岩,. 地震条件下危岩崩塌运动特性的初步探讨[J]. , 2017, 38(2): 583-592. |
[12] | 闫澍旺,林 澍,霍知亮,楚 剑,郭 伟,. 桶形基础液压下沉过程的耦合欧拉-拉格朗日有限元法分析[J]. , 2017, 38(1): 247-252. |
[13] | 张友良,谭 飞,张礼仁,施明明. 岩土工程亿级单元有限元模型可扩展并行计算[J]. , 2016, 37(11): 3309-3316. |
[14] | 张 蕾 ,唐小松 ,李典庆 ,曹子君 , . 基于Copula函数的岩土结构物系统可靠度分析[J]. , 2016, 37(1): 193-202. |
[15] | 付晓东,盛 谦,张勇慧. 非连续变形分析方法中的阻尼问题研究[J]. , 2015, 36(7): 2057-2062. |
|