岩土力学 ›› 2019, Vol. 40 ›› Issue (2): 785-791.doi: 10.16285/j.rsm.2017.1517

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

基于离散元数值方法的砂土小应变弹性特性探讨

顾晓强1, 2,杨朔成1, 2   

  1. 1. 同济大学 地下建筑与工程系,上海 200092;2. 同济大学 岩土与地下工程教育部重点实验室,上海 200092
  • 收稿日期:2017-07-18 出版日期:2019-02-11 发布日期:2019-02-19
  • 作者简介:顾晓强,男,1981年生,博士,副研究员,主要从事土动力学和宏微观土力学等方面的研究。
  • 基金资助:
    国家重点研发计划(No. 2016YFC0800200);国家自然科学基金(No. 41772283;No. 51822809);中央高校基本科研业务费。

Numerical investigation on the elastic properties of granular soils by discrete element method

GU Xiao-qiang1, 2, YANG Shuo-cheng1, 2   

  1. 1. Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China; 2. Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China
  • Received:2017-07-18 Online:2019-02-11 Published:2019-02-19
  • Supported by:
    This work was supported by the National Key R&D Program of China (2016YFC0800200), the National Natural Science Foundation of China (41772283, 51822809) and the Fundamental Research Funds for the Central Universities.

摘要: 土体在非常小应变(小于10?6)时表现为线弹性。以往大量试验研究表明,土体的弹性特性主要受到土体有效围压、孔隙比、颗粒级配和试样制备方法等影响。通过离散元方法模拟Duffy和Mindlin中的物理试验,建立球形颗粒规则排列的三维模型,继而进行小变形条件下的数值加载试验来确定弹性特性,并从细观层面揭示决定土体弹性特性的内在因素。模拟结果表明:粒状土的弹性模量取决于土体的颗粒接触配位数、法向接触力的大小和分布,跟宏观土体孔隙比、有效围压作用相对应。弹性模量应力指数n高于Hertz-Mindlin接触法则对应的1/3,主要是由于应力增加过程中配位数的增加和接触力趋于更均匀导致。同时,泊松比也随着配位数的增加而减小,并非保持不变,跟试验结果趋势一致。

关键词: 离散元, 弹性模量, 泊松比, 配位数, 接触力

Abstract: It’s well recognized that the soil response is generally linear and elastic at very small strain level (i.e. below 10-6). Previous numerous experiments showed that the soil elastic properties mainly depended on the effective confining pressure, void ratio, soil grading and sample preparation method. In this paper, the physical experiments on regular packing of steel spheres by Duffy and Mindlin were simulated by three dimensional discrete element method in order to reveal the fundamental mechanism related to soil elastic properties at meso-scale. The simulations showed that the elastic stiffness of the soil depended on the coordination number, the normal contact force and its distribution, which correspond to the macroscopic parameters of void ratio and effective confining pressure. It was found that the stress exponent n was higher than 1/3 as predicted by the Hertz-Mindlin contact model and it was mainly due to the increase of coordination number and uniformity of contact force distribution during the increase of confining pressure. Meanwhile, Poisson's ratio decreased with increasing coordination number instead of a constant as commonly assumed, which was consistent with experimental results.

Key words: discrete element method, elastic modulus, Poisson’s ratio, coordination number, contact force

中图分类号: 

  • TU 443
[1] 薛亚东, 周杰, 赵丰, 李兴. 基于MatDEM的TBM滚刀破岩机理研究[J]. 岩土力学, 2020, 41(S1): 337-346.
[2] 鲍宁, 魏静, 陈建峰. 桩承式路堤土拱效应三维离散元分析[J]. 岩土力学, 2020, 41(S1): 347-354.
[3] 周跃峰, 杨哲, 饶锡保, 肖国强, 周黎明. 不同含水率下黏土弹性参数的 弯曲−伸展元试验研究[J]. 岩土力学, 2020, 41(S1): 387-393.
[4] 吴再海, 纪洪广, 姜海强, 齐兆军, 寇云鹏, . 尾砂胶结含盐冻结充填体力学特性研究[J]. 岩土力学, 2020, 41(6): 1874-1880.
[5] 赵军, 郭广涛, 徐鼎平, 黄翔, 胡偲, 夏跃林, 张頔. 三轴及循环加卸载应力路径下深埋 硬岩变形破坏特征试验研究[J]. 岩土力学, 2020, 41(5): 1521-1530.
[6] 徐东升, 黄明, 黄佛光, 陈成. 不同级配珊瑚砂水泥胶结体的破坏行为分析[J]. 岩土力学, 2020, 41(5): 1531-1539.
[7] 吴祁新, 杨仲轩. 基于应变响应包络的颗粒材料增量力学行为研究[J]. 岩土力学, 2020, 41(3): 915-922.
[8] 孙博, 杨怀德, 谷玲, 刘跃, 唐碧华, 赵桂连, 张连明. 基于UDEC颗粒模型的不确定性分析[J]. 岩土力学, 2019, 40(9): 3679-3688.
[9] 马秋峰, 秦跃平, 周天白, 杨小彬. 多孔隙岩石加卸载力学特性及本构模型研究[J]. 岩土力学, 2019, 40(7): 2673-2685.
[10] 王蕴嘉, 宋二祥. 堆石料颗粒形状对堆积密度及强度影响的 离散元分析[J]. 岩土力学, 2019, 40(6): 2416-2426.
[11] 赵丁凤, 梁 珂, 陈国兴, 熊 浩, 周正龙, . 剪切-体积应变耦合的孔压增量模型试验[J]. 岩土力学, 2019, 40(5): 1832-1840.
[12] 陈 峥, 何 平, 颜杜民, 高红杰, . 考虑土拱效应的管棚合理间距计算方法[J]. 岩土力学, 2019, 40(5): 1993-2000.
[13] 杨石扣, 张继勋, 任旭华, . 基于改进数值流形法的接触裂纹问题研究[J]. 岩土力学, 2019, 40(5): 2016-2021.
[14] 赵兰浩, 芮开天, 刘勋楠. 非均匀颗粒快速线性接触检测算法[J]. 岩土力学, 2019, 40(3): 1187-1196.
[15] 李静, 孔祥超, 宋明水, 汪勇, 王昊, 刘旭亮, . 储层岩石微观孔隙结构对岩石力学特性 及裂缝扩展影响研究[J]. 岩土力学, 2019, 40(11): 4149-4156.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!