岩土力学 ›› 2022, Vol. 43 ›› Issue (S1): 580-590.doi: 10.16285/j.rsm.2020.1194

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

离散单元法软件EDEM中接触半径对岩石 力学特性的影响及其应用

申浩翰,张海,范俊锴,徐瑞阳,张小明   

  1. 河南理工大学 机械与动力工程学院,河南 焦作 454002
  • 收稿日期:2020-08-12 修回日期:2021-03-05 出版日期:2022-06-30 发布日期:2022-07-15
  • 作者简介:申浩翰,男,1996年生,硕士研究生,主要从事岩石冲击破碎方面的研究。
  • 基金资助:
    国家自然科学基金项目(No.51975189)。

Influence of contact radius on rock mechanical property and its application in discrete element method software EDEM

SHEN Hao-han, ZHANG Hai, FAN Jun-kai, XU Rui-yang, ZHANG Xiao-ming   

  1. School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo, Henan 454002, China
  • Received:2020-08-12 Revised:2021-03-05 Online:2022-06-30 Published:2022-07-15
  • Supported by:
    This work was supported by the National Natural Science Foundation of China(51975189).

摘要: 以EDEM软件的Hertz-Mindlin接触模型为对象,对接触半径与黏结键强度的关系及其对岩石抗拉和抗压强度的影响进行研究,提出了一种建立高压拉强度比岩石离散元模型的方法。根据接触模型的强度准则和几何准则,提出了黏结键受拉状态下的“强键”、“弱键”以及临界接触半径系数的概念,推导出了接触颗粒间法向运动和切向运动时的临界接触半径表达式。以接触半径系数a 为自变量进行了岩石单轴压缩和单轴拉伸仿真试验,结果表明:抗拉强度随a 变化呈三阶段分段线性关系,1.000 0≤a<1.000 8区间内,曲线斜率很大;1.000 8≤ a ≤1.002 5区间内,曲线斜率较小;a>1.002 5后,抗拉强度几乎不变化。抗压强度随a 变化呈两阶段分段线性关系,在1.000 0≤a ≤1.000 2区间内,曲线斜率很大;当a>1.000 2后,抗压强度几乎不变化。在1.000 1≤a ≤1.000 3区间内,抗压强度和抗拉强度之比约为6.5~14.0,可以较好地拟合真实岩石的抗拉和抗压特性,从而构建出具有高压拉强度比的岩石离散元模型。

关键词: EDEM软件, 离散元法, 岩石, 接触半径

Abstract: Based on the Hertz-Mindlin contact model in EDEM, the relationship between contact radius and bond strength, and the influence of contact radius on tensile and compressive strengths of rock model were investigated. A method for establishing discrete element model of rock with high compressive-tensile strength ratio was suggested. According to the strength and geometric criteria of the contact model in EDEM, the concepts of “strong bond”, “weak bond” and critical contact radius coefficient were proposed. The formulae of critical contact radius when two particles in contact move along normal or tangential directions were derived. The rock uniaxial compression and tension simulation experiments were carried out with the contact radius coefficient α as an independent variable, and the results show that: the tensile strength shows a three-stage piecewise linear relationship with the change of α. At the interval of 1.000 0≤<1.000 8, the slope of the curve is very large. At the interval of 1.000 8≤≤1.002 5, the slope of the curve is small. When >1.002 5, the tensile strength is almost constant. The compressive strength shows a two-stage piecewise linear relationship with the change of α. At the interval of 1.000 0≤ ≤1.000 2, the slope of the curve is very large, and after >1.000 2, the compressive strength is also almost constant. In the interval of 1.000 1≤≤1.000 3, the compressive-tensile strength ratio is about 6.5 to 14.0, and it can be used to reproduce the tensile and compressive properties of real rock, so as to construct the discrete element model of rock with high compressive-tensile strength ratio.

Key words: EDEM software, discrete element method, rock, contact radius

中图分类号: 

  • TU452
[1] 岑夺丰, 刘畅, 黄达. 灰岩层面拉剪力学特性及层面起伏效应研究[J]. 岩土力学, 2022, 43(S1): 77-87.
[2] 朱星, 刘汉香, 胡桔维, 范杰, . 砂岩破坏声发射临界慢化前兆特征试验研究[J]. 岩土力学, 2022, 43(S1): 164-172.
[3] 匡智浩, 李邵军, 杜灿勋, 邱士利, 吝曼卿, 杜三林, . 考虑应力变化速率的岩石脆性评价指标[J]. 岩土力学, 2022, 43(S1): 293-300.
[4] 刘学伟, 刘泉声, 汪志强, 刘滨, 康永水, 王传兵, . 基于格栅拱架的软岩巷道分步联合控制技术研究[J]. 岩土力学, 2022, 43(S1): 469-478.
[5] 杨磊, 涂冬媚, 朱启银, 吴则祥, 余闯, . 考虑变温幅值影响的颗粒循环热 固结离散元法试验研究[J]. 岩土力学, 2022, 43(S1): 591-600.
[6] 李然, 王圣涛, 张顶立, 陈平, 潘红桂, 李奥, . 小净距隧道中夹岩对拉锚杆控制机制与工程应用[J]. 岩土力学, 2022, 43(7): 1865-1876.
[7] 王刚, 宋磊博, 刘夕奇, 包春燕, 吝曼卿, 刘广建, . 非贯通节理花岗岩剪切断裂力学特性及 声发射特征研究[J]. 岩土力学, 2022, 43(6): 1533-1545.
[8] 崔国建, 张传庆, 周辉, 卢景景, 高阳, 胡明明, 胡大伟, . 动力扰动作用下多功能岩石结构面剪切 试验装置研制与应用研究[J]. 岩土力学, 2022, 43(6): 1727-1737.
[9] 唐旭海, 许婧璟, 张怡恒, 何琦, 王正直, 张国平, 刘泉声, . 基于微观岩石力学试验和NWA13618陨石的 小行星岩石力学参数分析[J]. 岩土力学, 2022, 43(5): 1157-1163.
[10] 姜玥, 周辉, 卢景景, 高阳, . 空心圆柱砂岩真三轴试验研究[J]. 岩土力学, 2022, 43(4): 932-944.
[11] 程坦, 郭保华, 孙杰豪, 田世轩, 孙崇轩, 陈岩, . 非规则砂岩节理剪切变形本构关系试验研究[J]. 岩土力学, 2022, 43(1): 51-64.
[12] 张超, 杨楚卿, 白允. 岩石类脆性材料损伤演化分析及其模型方法研究[J]. 岩土力学, 2021, 42(9): 2344-2354.
[13] 杜崧, 肖明, 陈俊涛, . 洞室块体危险性分析的突变级数评价法研究[J]. 岩土力学, 2021, 42(9): 2578-2588.
[14] 申浩翰, 张海, 范俊锴, 徐瑞阳, 张小明. 基于多参数拟合的EDEM岩石建模方法[J]. 岩土力学, 2021, 42(8): 2298-2310.
[15] 李博, 崔逍峰, 莫洋洋, 邹良超, 伍法权, . 法向应力作用下砂岩错位裂隙变形行为研究[J]. 岩土力学, 2021, 42(7): 1850-1860.
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 .