›› 2018, Vol. 39 ›› Issue (7): 2639-2650.doi: 10.16285/j.rsm.2017.1702

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

三维距离势离散单元法

刘勋楠1,赵兰浩1,毛 佳1,许 栋2   

  1. 1. 河海大学 水利水电学院,江苏 南京210098;2. 天津大学 水利工程仿真与安全国家重点实验室,天津 300072
  • 收稿日期:2017-08-16 出版日期:2018-07-10 发布日期:2018-08-05
  • 通讯作者: 赵兰浩,男,1980年生,工学博士,教授,主要从事流固耦合和水工结构抗震的工作。E-mail: zhaolanhao@hhu.edu.cn E-mail:liuxunnanhhu@163.com
  • 作者简介:刘勋楠,男,1990年生,博士研究生,主要从事离散元方面的研究工作。
  • 基金资助:

    国家自然科学基金(No. 51279050);霍英东教育基金会第十五届高等院校青年教师基金(No. 151073);水利部公益性行业科研专项经费项目(No. 201501034);国家重点研发计划(No. 2016YFC0401601)。

Discrete element method using three dimensional distance potential

LIU Xun-nan1, ZHAO Lan-hao1, MAO Jia1, XU Dong2   

  1. 1. College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing, Jiangsu 210098, China; 2. State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China
  • Received:2017-08-16 Online:2018-07-10 Published:2018-08-05
  • Supported by:

    This work was supported by the National Natural Science Foundation of China (51279050), the 15th Fok Ying-Tong Education Foundation for Young Teachers in the Higher Education Institutions of China (151073), the Non-profit Industry Financial Program of MWR of China (201501034) and the National Key R & D Program of China (2016YFC0401601).

摘要: 通过定义距离势函数,提出一种适用于空间任意凸多面体单元的三维距离势函数离散单元法。该方法采用归一化的计算方式,将势函数表征为接触体间的距离函数,并基于此建立接触力计算方程,明确了势函数的物理意义,使接触力计算更加合理,无需对各种可能的接触形式进行差异化处理。新方法克服了原有势函数物理意义不明、接触力计算受单元形式影响等重要缺陷,并突破了四面体单元的限制,可采用空间任意凸多面体单元。通过若干算例说明新方法的正确性和有效性。数值模拟结果表明新方法能够很好处理空间任意多面体单元复杂接触变换过程,准确处理复杂非连续介质的运动过程。

关键词: 距离势函数, 任意凸多面体块体单元, 接触力, 距离势离散单元法

Abstract: A 3-D discrete element method is presented on the basis of new definition of distance potential function, which is suitable for the arbitrary polyhedral elements. In this approach, the potential function is characterized as a distance function between the contacts using a normalized computational algorithm. Meanwhile, the calculation formulations for contact interaction are deduced with the new function. Thus this new distance potential function implies sound physical meaning, and is reasonable and concise to avoid different contact situations. This method overcomes the problems of the original potential method in which the physical meaning of potential function is not clear and the calculation of contact force is affected by the element form. Furthermore, it gets rid of the restrain of tetrahedral shape and can be applied to arbitrary polyhedrons. Several examples are analyzed to demonstrate the accuracy and efficiency of this new method. This method shows good capability of handling the problem of contact transformation for arbitrary polyhedrons and the motion process of the complex discontinuous media.

Key words: distance potential function, arbitrary convex polyhedral blocks, contact interaction, distance potential discrete element method

中图分类号: 

  • O 158

[1] 杨石扣, 张继勋, 任旭华, . 基于改进数值流形法的接触裂纹问题研究[J]. 岩土力学, 2019, 40(5): 2016-2021.
[2] 顾晓强, 杨朔成, . 基于离散元数值方法的砂土小应变弹性特性探讨[J]. 岩土力学, 2019, 40(2): 785-791.
[3] 严成增,郑 宏,葛修润. 基于统一标定的势接触力计算[J]. , 2015, 36(1): 249-256.
[4] 朱小明 ,李海波 ,刘 博,. 循环剪切荷载作用下含二阶起伏体模拟岩石节理力学特性研究[J]. , 2014, 299(2): 371-379.
[5] 徐栋栋,孙冠华,郑 宏. 接触处理Munjiza方法在数值流形法中的应用[J]. , 2013, 34(2): 526-534.
[6] 蒋明镜,孙渝刚,李立青. 胶结颗粒接触力学特性测试装置研制[J]. , 2011, 32(1): 309-315.
[7] 陈万祥,郭志昆,黄 育. 岩土颗粒介质在应力波作用下的接触力与变形关系[J]. , 2010, 31(8): 2580-2586.
[8] 何思明,吴 永,李新坡. 滚石冲击碰撞恢复系数研究[J]. , 2009, 30(3): 623-627.
[9] 李建红,沈珠江. 结构性土的微观破损机理研究[J]. , 2007, 28(8): 1525-1532.
[10] 江 巍,郑 宏. 非连续变形分析方法中人为参数的影响[J]. , 2007, 28(12): 2603-2606.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!