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

作     者：赵强 焦玉勇 张秀丽 谢壁婷 王龙 黄刚海 ZHAO Qiang;JIAO Yu-yong;ZHANG Xiu-li;XIE Bi-ting;WANG Long;HUANG Gang-hai

作者机构：中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室湖北武汉430071 中国科学院大学北京100049 中国地质大学(武汉)工程学院湖北武汉430074 国家电网山东省电力公司经济技术研究院山东济南250021 湖南科技大学土木工程学院湖南湘潭411201 

出 版 物：《岩土力学》 (Rock and Soil Mechanics)

年 卷 期：2019年第40卷第11期

页      面：4515-4522页

核心收录：

学科分类：07[理学] 070102[理学-计算数学] 0701[理学-数学] 

基　　金：国家重点基金项(No.41731284) 国家自然科学基金面上资助项目(No.11672360)~~ 

主　　题：非连续变形分析 显式时间积分 质量矩阵 最大位移准则 岩土工程 

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