基于静态松弛法求解的三维离散单元法.pdfVIP

( ) 第19 卷　第 4 期 岩石力学与工程学报 19 4 : 453～ 458 2000 年 7 月 Ch inese J ou rna l of R ock M echan ics and E ng ineering J u ly , 2000 基于静态松弛法求解的三维离散单元法 1 1, 2 焦玉勇 　葛修润 ( 1 中国科学院武汉岩土力学研究所　武汉　430071) (2 上海交通大学建工学院　上海　200030) 摘要　将基于静态松弛解法的离散单元法推广至三维。论述了三维静态松弛离散单元法的特点和求解过程, 建立了 块体的静力平衡方程, 并计算了两个算例。结果表明所提出的计算公式是有效的, 静态松弛法比动态松弛法具有更 高的效率。 关键词　静态松弛, 三维, 离散单元法, 平衡方程 分类号　 242. 26　　　　文献标识码　 　　　　文章编号(2000) 04 0453 06 O A 1　引　言 2　静态松弛法的特点 Cundall ( 1971) 提出的离散单元法是基于动态松 所谓松弛法, 指的是块体位移的求解方法, 静态 弛解法的, 它采用牛顿第二定律或的差分格式求解 松弛法与动态松弛法的根本区别就在于块体位移的 块体的位移, 并用时步积分累积求算块体的大位移, 求解方法不同。众所周知, 动态松弛法的出发点是牛 因而动态松弛法存在着动能的耗散和时间步长的选 顿第二定律, 由加速度→速度→位移的求解格式使 取问题, 一般认为, 时步和阻尼的确定是动态松弛解 得阻尼系数和时间因素的引入成为不可避免, 继而 法中最为困难的问题之一。针对这一弊端, Stew art 导致系统振荡控制、加快收敛速度、防止块体相互嵌 ( 1982) 提出了静态松弛离散单元法。静态松弛法根 入等一系列问题。与此相对照, 静态松弛法则采取了 据块体产生不平衡力后达到再平衡时的力位移关系 较为简单的处理方法: 根据产生不平衡力的块体达 建立平衡方程组, 通过求解方程组得到块体的位移。 到再平衡时的力位移关系建立平衡方程组, 由方程 由于采用隐式格式和静力解法, 静态松弛法完全避 组隐式求解块体的位移。静态松弛法避开了动态松 [ 1 ] 开了时步和阻尼 , 并且计算效率大大提高。之后, 弛法在求解位移过程中的诸多不便, 其代价仅是引 文[2 ]研究了静态松弛离散单元法的理论和算法; 专 入一个 6 ×6 的线性方程组。两种松弛法的计算循环 著[3 ] 中辟专门章节对静态松弛法进行了介绍; 文 如图 1, 2 所示。 [4 ]提出了静态同步松弛的离散单元法; 文[ 5