岩土力学 ›› 2019, Vol. 40 ›› Issue (11): 4515-4522.doi: 10.16285/j.rsm.2018.1515

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

基于显式时间积分的球颗粒DDA计算方法

赵强1, 2,焦玉勇3,张秀丽1, 2,谢壁婷1, 2,王龙4,黄刚海5   

  1. 1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071;2. 中国科学院大学,北京 100049; 3. 中国地质大学(武汉)工程学院,湖北 武汉 430074;4. 国家电网山东省电力公司经济技术研究院,山东 济南 250021; 5. 湖南科技大学 土木工程学院,湖南 湘潭 411201
  • 收稿日期:2018-08-20 出版日期:2019-11-11 发布日期:2019-12-01
  • 通讯作者: 焦玉勇,男,1968年生,博士,教授,博士生导师,主要从事岩土力学领域中数值计算方法的研究。E-mail: ddarf@qq.com E-mail:zhaoqiang_dda@qq.com
  • 作者简介:赵强,男,1990年生,博士研究生,主要从事岩土力学领域中数值计算方法的研究
  • 基金资助:
    国家重点基金项(No. 41731284);国家自然科学基金面上资助项目(No. 11672360)。

Explicit time integration based spherical DDA calculation method

ZHAO Qiang1, 2, JIAO Yu-yong3, ZHANG Xiu-li1, 2, XIE Bi-ting1, 2, WANG Long4, HUANG Gang-hai5   

  1. 1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Engineering School, China University of Geosciences, Wuhan, Hubei 430074, China; 4. State Grid Electric Power Company Economic and Technological Research Institute of Shandong, Jinan, Shandong 250021, China; 5. School of Civil Engineering, Hunan University of Science and Technology, Xiangtan, Hunan 411201, China
  • Received:2018-08-20 Online:2019-11-11 Published:2019-12-01
  • Supported by:
    This work was supported by the Key Program for National Natural Science Foundation of China (41731284) and the National Natural Science Foundation of China General Program (11672360).

摘要: 非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。

关键词: 非连续变形分析, 显式时间积分, 质量矩阵, 最大位移准则, 岩土工程

Abstract: Discontinuous deformation analysis (DDA) is a new numerical method parallel to the finite element method. Based on the principle of minimum potential energy, the method unifies the deformation, motion and contact of each discrete block into the equilibrium equation for implicit solution. However, the traditional DDA method requires assembling the whole stiffness matrix to solve the equations in parallel, which occupies a large amount of memory, takes a long time and has a very low computational efficiency in the three-dimensional numerical simulation of large geotechnical engineering problems. Therefore, a three-dimensional spherical DDA method based on explicit time integration is proposed. This method does not require assembling the whole stiffness matrix. Because of the mass matrix is diagonal matrix, it can be stored as a one-dimensional vector with less memory, and can be solved block by block with higher efficiency in the acceleration solving process. The maximum displacement criterion is used to simplify the contact algorithm in the contact judgment, and the smaller time step ensures the accuracy of calculation. Finally, the accuracy and efficiency of the method are verified by several classical examples.

Key words: discontinuous deformation, explicit time integration, mass matrix, maximum displacement criterion, geotechnical engineering

中图分类号: 

  • O 242
[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] 范文亮,王余乐,魏奇科,杨朋超,李正良, . 岩土工程可靠度分析的改进四阶矩方法[J]. , 2018, 39(4): 1463-1468.
[6] 黄明华,赵明华,陈昌富. 锚固长度对锚杆受力影响分析及其临界值计算[J]. , 2018, 39(11): 4033-4041.
[7] 熊自明,卢 浩,王明洋,钱七虎,戎晓力,. 我国大型岩土工程施工安全风险管理研究进展[J]. , 2018, 39(10): 3703-3716.
[8] 刘泉声,蒋亚龙,何 军. 非连续变形分析的精度改进方法及研究趋势[J]. , 2017, 38(6): 1746-171.
[9] 方砚兵,苏永华,肖 旺,梁 斌. 基于子区间法的隐式功能函数非概率可靠性方法研究[J]. , 2017, 38(4): 1171-1178.
[10] 黄小福,张迎宾,赵兴权,余鹏程,邢 昊,张 珏,陈岩岩,. 地震条件下危岩崩塌运动特性的初步探讨[J]. , 2017, 38(2): 583-592.
[11] 闫澍旺,林 澍,霍知亮,楚 剑,郭 伟,. 桶形基础液压下沉过程的耦合欧拉-拉格朗日有限元法分析[J]. , 2017, 38(1): 247-252.
[12] 付晓东,盛 谦,张勇慧,冷先伦, . 非连续变形分析(DDA)线性方程组的高效求解算法[J]. , 2016, 37(4): 1171-1178.
[13] 张友良,谭 飞,张礼仁,施明明. 岩土工程亿级单元有限元模型可扩展并行计算[J]. , 2016, 37(11): 3309-3316.
[14] 张 蕾 ,唐小松 ,李典庆 ,曹子君 , . 基于Copula函数的岩土结构物系统可靠度分析[J]. , 2016, 37(1): 193-202.
[15] 付晓东,盛 谦,张勇慧. 非连续变形分析方法中的阻尼问题研究[J]. , 2015, 36(7): 2057-2062.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!