岩土力学 ›› 2022, Vol. 43 ›› Issue (7): 1961-1968.doi: 10.16285/j.rsm.2021.1754

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

可变形圆化多边形离散单元法

周洋诗琦,赵兰浩,邵琳玉,毛佳   

  1. 河海大学 水利水电学院,江苏 南京 210098
  • 收稿日期:2021-10-18 修回日期:2022-05-10 出版日期:2022-07-26 发布日期:2022-08-05
  • 通讯作者: 赵兰浩,男,1980年生,博士,教授,主要从事流固耦合和水工结构抗震方面的研究工作。E-mail: zhaolanhao@hhu.edu.cn E-mail:zysq@hhu.edu.cn
  • 作者简介:周洋诗琦,女,1998年生,硕士研究生,主要从事离散元方面的研究工作。
  • 基金资助:
    国家自然科学基金(No. 52009034);霍英东教育基金会第十五届高等院校青年教师基金(No. 151073);中央高校基本科研业务费专项基金(No. B210201036)。

A deformable spheropolygon-based discrete element method

ZHOU Yang-shi-qi, ZHAO Lan-hao, SHAO Lin-yu, MAO Jia   

  1. College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing, Jiangsu 210098, China
  • Received:2021-10-18 Revised:2022-05-10 Online:2022-07-26 Published:2022-08-05
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (52009034), the 15th Fok Ying-Tong Education Foundation for Young Teachers in the Higher Education Institutions of China (151073) and the Fundamental Research Funds for the Central Universities (B210201036).

摘要: 将基于圆化多边形离散单元法与有限元方法结合,提出一种可变形圆化多边形离散单元法。此法对块体离散元进行圆化处理,可较好地表征不规则块体外形,又保留了颗粒离散元计算高效的优势。在求解接触力时,消除了角点处法向奇异等问题,同时增强计算的稳定性和简化接触判断。同时对切向接触力计算模型进行修正,使得接触力计算效率得到提高。此法突破了圆化多边形刚体假设的限制,可以精确计算任意形状不规则离散单元之间的相互作用,对单元的运动和变形进行模拟。通过超静定梁冲击试验、不规则块体单轴压缩试验和料斗流动“卡阻”试验3个数值模拟算例,论证此法可以有效地捕捉单元的碰撞、分离和变形等空间运动和自身特性以及其细观力学表征。

关键词: 圆化多边形离散单元法, 可变形单元, 接触力

Abstract: A deformable spheropolygon-based discrete element method (DSEM) which combines the spheropolygon-based discrete element method (DEM) and the finite element method (FEM) is proposed. In this method, with the rounded discrete elements, the shape of the irregular blocks can be well characterized while the high efficiency advantage of particle discrete element method is also kept. The proposed method can not only eliminate the normal singularity at the corner but also improve the stability and simplify the contact judgment during contact force calculation. Meanwhile, the calculation model of the tangential contact force is also modified to improve the efficiency of contact force calculation. The DSEM breaks through the limitation of the rigid body hypothesis of the spheropolygon. Thus, the interaction between discrete elements with arbitrary shape can be accurately calculated, and the movement and deformation of the elements can also be simulated. Finally, three numerical examples including the dynamic response of statically indeterminate beam, the uniaxial compression test of irregular block, and hopper flow test is simulated using the DSEM to prove that the proposed method can effectively capture both the spatial motion as well as characteristics of elements (such as collision, separation and deformation) and their microscopic mechanical characterization.

Key words: spheropolygon-based DEM, deformable particles, contact force

中图分类号: 

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


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 石玉玲,门玉明,彭建兵,黄强兵,刘洪佳. 地裂缝对不同结构形式桥梁桥面的破坏试验研究[J]. , 2009, 30(10): 2917 -2922 .
[5] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[6] 徐速超,冯夏庭,陈炳瑞. 矽卡岩单轴循环加卸载试验及声发射特性研究[J]. , 2009, 30(10): 2929 -2934 .
[7] 张力霆,齐清兰,魏静,霍倩,周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[8] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[9] 易 俊,姜永东,鲜学福,罗 云,张 瑜. 声场促进煤层气渗流的应力-温度-渗流压力场的流固动态耦合模型[J]. , 2009, 30(10): 2945 -2949 .
[10] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .