离散单元法的计算参数和求解方法选择Ξ.pdfVIP

第 16 卷第 1 期 计 算 力 学 学 报 V o l. 16 N o. 1 1999 年 2 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 1999 离散单元法的计算参数和求解方法选择 邢纪波　　俞良群　　张瑞丰　　　　　王泳嘉　　　　 ( ) ( ) 　　　　 烟台大学土木工程系, 烟台, 264005 　　　　　 东北大学资源及土木工程学院, 沈阳, 110006 摘　要　本文对离散单元法的一些基本问题进行了研究, 给出了刚度计算公式, 提出了 两种新的用于求解静态问题的自适应阻尼, 并以自由落体运动为例, 探讨了中心差分 法、 法和 法对离散单元法解的影响。 W ilson N ewm ark 关键词　离散单元法; 计算参数选择; 直接积分法 分类号　TB 115;O 242 2 1　前　　言 离散单元法是一种适用于解决非连续介质力学问题的数值计算方法, 已在结构工程、岩土 工程、散体力学以及爆炸力学等领域得到了成功应用, 取得了一些令人瞩目的研究成果[ 1, 2, 3 ]。 离散单元法假定单元之间的作用力与它们之间的相对位移量成正比, 然后应用牛顿力学 建立单元的运动方程, 最后采用中心差分法进行显示迭代求解。离散单元法的计算原理虽然简 单, 但是在采用该方法求解具体问题时, 人们感到最棘手的难题是如何选择计算参数, 如接触 刚度和粘性阻尼系数等。 在离散单元法计算中, 单元之间的接触刚度一般由实验方法测定, 不仅费时费力, 而且有 时测得的值又不能满足计算要求。本文根据离散单元法进行动态松弛求解时的限制条件, 即在 一个计算时步内, 任何单元所受的扰动只能传递到与它相连接的单元, 推导出接触刚度的理论 计算公式。 以往用离散单元法求解块体系统的静平衡问题时, 往往采用与速度成正比的粘性阻尼吸 收块体的动能。此时所遇到的困难是求系统的最佳阻尼比和振动圆频率, 由于这两个参数不能 采用求连续介质振动的特征值方法计算, 而往往采用试算方式确定, 因此应用起来显得极为不 便。此外, 一旦确定了这些阻尼参数, 则阻尼将被等同地施加于所有块体上, 实际上这是完全不 必要的。因为所研究的问题常常表现为一部分块体处于平衡状态, 而另一部分块体则处于非稳 定流动状态。因此, 较为合理的作法是对不同状态的块体取不同的阻尼值。正是基于以上分析, 本文提出了两种自适应阻尼以吸收块体的动能。 在求解方法方面, 以往研究者多采用动态松弛法求解离散单元的运动方程。本文以自由落 体运动为例, 探讨了用中心差分法、 法和 法求解离散单元运动方程时, 如 W ilson N ewm ark 何选取计算时步、迭代次数和阻尼系数等问题, 并比较了它们各自的优缺点。 山东省自然科学基金资助项目 　收稿日期: 修改稿收到日期: 　邢纪波: 男, 1962 年生, 博士, 副教授 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 48 计 算 力 学 学 报 　　16 卷 2　接触刚度 2 1　法向刚度 当一个单元处于平衡状态时, 该单元所受的各接触力是一定的。由于单元之间的接触力与 它们之间的叠合值及刚度成正比, 当接触刚度大时, 单元之间的叠合量则小; 反之, 当接触刚度 小时, 单元之间的叠合量则大。因此, 接触刚度无须精确计算。本文根据离