›› 2014, Vol. 35 ›› Issue (8): 2401-2407.

• 数值分析 • 上一篇    下一篇

基于OpenMP的非连续变形分析并行计算方法

付晓东,盛 谦,张勇慧   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2012-11-27 出版日期:2014-08-12 发布日期:2014-08-14
  • 作者简介:付晓东,男,1986年生,博士,主要从事岩土力学领域中数值计算方面的研究工作。
  • 基金资助:

    国家自然科学基金资助项目(No.11272331)资助;国家自然科学基金重大研究计划集成项目(No. 91215301)。

Parallel computing method for discontinuous deformation analysis using OpenMP

FU Xiao-dong, SHENG Qian, ZHANG Yong-hui   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2012-11-27 Online:2014-08-12 Published:2014-08-14

摘要: 非连续变形分析(DDA)方法严格满足平衡要求和能量守恒,具有完全的运动学及数值可靠性,但对大规模岩土工程问题的数值模拟耗时太长,尤其是线性方程组求解,并行计算可以很好地解决该问题。首先基于DDA方法的基本理论,阐述了适用于DDA方法中的基于块的行压缩法和基于“试验-误差”迭代格式的非零位置记录;其次,引入块雅可比迭代法并行求解DDA方法的线性方程组,并改进了相应的非零存储方法;最后,基于OpenMP实现了DDA线性方程组求解并行计算,并将其应用于地下洞室群的破坏过程分析,以加速比为并行效率的指标评价,结果表明,该并行计算策略可以极大提高DDA的计算效率,而且适合各种规模的问题。

关键词: 岩土工程, 非连续变形分析, 并行计算, OpenMP, 块雅可比迭代法, 非零存储

Abstract: The discontinuous deformation analysis (DDA) strictly meets balance requirements and energy conservation, and has complete kinematics and numerical reliability. However, the computation time of simulating large-scale geotechnical engineering is too long, especially for solving linear equations; the parallel computing will be a good solution to the problem. Firstly, two contents of non-zero storage in DDA have been described. One is the block compressed sparse row method, and the other is the iterative scheme of non-zero position recording based on the trial-error. Furthermore, we have introduced the block Jacobi iterative method to parallel solving linear equations of DDA, and improved the non-zero storage method. Finally, benefit from the sub-matrices storing, parallel computing of the DDA has been achieved based on Opening Multi-Processing. The improved DDA program has been used to simulate the failure process of the underground caverns; the speedup is obtained to evaluate the parallel efficiency. The results demonstrate that the parallel computing strategy can greatly improve the computational efficiency of DDA for the problems with various sizes. The investigations of this paper will promote DDA to the engineering simulation.

Key words: geotechnical engineering, discontinuous deformation analysis, parallel computing, OpenMP, block Jacobi iterative method, non-zero element memory

中图分类号: 

  • TU 45
[1] 徐栋栋, 邬爱清, 李 聪, 汪 斌, 蒋昱州, 曾 平, 杨永涛, . 破裂全过程模拟的改进非连续变形分析方法[J]. 岩土力学, 2019, 40(3): 1169-1178.
[2] 赵强, 焦玉勇, 张秀丽, 谢壁婷, 王龙, 黄刚海, . 基于显式时间积分的球颗粒DDA计算方法[J]. 岩土力学, 2019, 40(11): 4515-4522.
[3] 范文亮,王余乐,魏奇科,杨朋超,李正良, . 岩土工程可靠度分析的改进四阶矩方法[J]. , 2018, 39(4): 1463-1468.
[4] 黄明华,赵明华,陈昌富. 锚固长度对锚杆受力影响分析及其临界值计算[J]. , 2018, 39(11): 4033-4041.
[5] 熊自明,卢 浩,王明洋,钱七虎,戎晓力,. 我国大型岩土工程施工安全风险管理研究进展[J]. , 2018, 39(10): 3703-3716.
[6] 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171.
[7] 方砚兵,苏永华,肖 旺,梁 斌. 基于子区间法的隐式功能函数非概率可靠性方法研究[J]. , 2017, 38(4): 1171-1178.
[8] 黄小福,张迎宾,赵兴权,余鹏程,邢 昊,张 珏,陈岩岩,. 地震条件下危岩崩塌运动特性的初步探讨[J]. , 2017, 38(2): 583-592.
[9] 闫澍旺,林 澍,霍知亮,楚 剑,郭 伟,. 桶形基础液压下沉过程的耦合欧拉-拉格朗日有限元法分析[J]. , 2017, 38(1): 247-252.
[10] 付晓东,盛 谦,张勇慧,冷先伦, . 非连续变形分析(DDA)线性方程组的高效求解算法[J]. , 2016, 37(4): 1171-1178.
[11] 王 翀 ,刘承论,. 基于3D-FSM的预紧力端锚数值模拟及其多核并行化研究[J]. , 2016, 37(3): 875-882.
[12] 张友良,谭 飞,张礼仁,施明明. 岩土工程亿级单元有限元模型可扩展并行计算[J]. , 2016, 37(11): 3309-3316.
[13] 张 蕾 ,唐小松 ,李典庆 ,曹子君 , . 基于Copula函数的岩土结构物系统可靠度分析[J]. , 2016, 37(1): 193-202.
[14] 付晓东,盛 谦,张勇慧. 非连续变形分析方法中的阻尼问题研究[J]. , 2015, 36(7): 2057-2062.
[15] 董威信,王翔南,王 远,于玉贞. 三维过渡等参单元在岩土工程有限元分析中的应用[J]. , 2015, 36(5): 1455-1462.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘 晓,唐辉明,刘 瑜. 基于集对分析和模糊马尔可夫链的滑坡变形预测新方法研究[J]. , 2009, 30(11): 3399 -3405 .
[2] 胡大伟,周 辉,谢守益,张 凯,邵建富,冯夏庭. 大理岩破坏阶段Biot系数研究[J]. , 2009, 30(12): 3727 -3732 .
[3] 师旭超,韩 阳. 卸荷作用下软黏土回弹吸水试验研究[J]. , 2010, 31(3): 732 -736 .
[4] 魏 纲,郭志威,魏新江,陈伟军. 软土隧道盾构出洞灾害的渗流应力耦合分析[J]. , 2010, 31(S1): 383 -387 .
[5] 朱建明,彭新坡,姚仰平,徐金海. SMP准则在计算煤柱极限强度中的应用[J]. , 2010, 31(9): 2987 -2990 .
[6] 白 冰,李小春,石 露,唐礼忠. 弹塑性应力-应变曲线的斜率恒等式及其验证和应用[J]. , 2010, 31(12): 3789 -3792 .
[7] 李占海,朱万成,冯夏庭,李绍军,周 辉,陈炳瑞. 侧压力系数对马蹄形隧道损伤破坏的影响研究[J]. , 2010, 31(S2): 434 -441 .
[8] 尹光志,王登科,张东明,魏作安. 基于内时理论的含瓦斯煤岩损伤本构模型研究[J]. , 2009, 30(4): 885 -889 .
[9] 侯公羽,牛晓松. 基于Levy-Mises本构关系及D-P屈服准则的轴对称圆巷理想弹塑性解[J]. , 2009, 30(6): 1555 -1562 .
[10] 蔡辉腾,危福泉,蔡宗文. 重庆主城区粉质黏土动力特性研究[J]. , 2009, 30(S2): 224 -228 .