距离势离散单元法.PDF

第 36 卷 第 6 期 岩石力学与工程学报 Vol.36 No.6 2017 年 6 月 Chinese Journal of Rock Mechanics and Engineering June ，2017 距离势离散单元法 1 1 1 2 刘勋楠 ，赵兰浩 ，毛 佳 ，许 栋 (1. 河海大学 水利水电学院，江苏 南京 210098 ；2. 天津大学 水利工程仿真与安全国家重点实验室，天津 300072) 摘要：势函数离散单元法在接触力计算中存在势函数物理意义不明确、接触力计算与客观实际不符等重大缺陷。 在已有研究的基础上，提出一种新型任意形状块体距离势离散单元法。该方法将势函数定义为接触体间归一化的 距离函数，可以客观衡量单元间的嵌入程度，具有明确的物理意义，真实反映单元的接触刚度，使法向接触力计 算更为合理。此外，新方法重新定义了切向方向矢量，提出切向接触力计算方法。新方法保留原势函数法所有优 点，同时完善接触力计算受单元形式影响的不足，并突破三角形单元的限制，能够应用于任意凸多边形单元。通 过若干算例说明新方法的正确性和有效性，可用于精确描述复杂非连续介质的运动过程。 关键词：数值分析；离散单元法；任意多边形块体单元；距离势函数；接触力 中图分类号：O 241 文献标识码：A 文章编号：1000 –6915(2017)06 –1353 –10 Discrete element method based on distance potential 1 1 1 2 LIU Xunnan ，ZHAO Lanhao ，MAO Jia ，XU Dong (1. College of Water Conservancy and Hydropower，Hohai University ，Nanjing ，Jiangsu 210098 ，China；2. State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300072 ，China) Abstract ：The discrete element method based on the potential function has many deficiencies. The physical meaning of the potential function is not clear and the contact force calculation is inconsistent with the reality. A novel discrete element method based on the distance potential for the arbitral convex polygonal elements was thus proposed using the new definition of the distance potential function. This method describes the potential function as a normalized distance function between the