留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

考虑颗粒级配和形态的颗粒柱坍塌特性离散元模拟

崔溦 魏杰 王超 王枭华 张社荣

崔溦, 魏杰, 王超, 王枭华, 张社荣. 考虑颗粒级配和形态的颗粒柱坍塌特性离散元模拟[J]. 机械工程学报, 2021, 43(12): 2230-2239. doi: 10.11779/CJGE202112009
引用本文: 崔溦, 魏杰, 王超, 王枭华, 张社荣. 考虑颗粒级配和形态的颗粒柱坍塌特性离散元模拟[J]. 机械工程学报, 2021, 43(12): 2230-2239. doi: 10.11779/CJGE202112009
CUI Wei, WEI Jie, WANG Chao, WANG Xiao-hua, ZHANG She-rong. Discrete element simulation of collapse characteristics of particle column considering gradation and shape[J]. JOURNAL OF MECHANICAL ENGINEERING, 2021, 43(12): 2230-2239. doi: 10.11779/CJGE202112009
Citation: CUI Wei, WEI Jie, WANG Chao, WANG Xiao-hua, ZHANG She-rong. Discrete element simulation of collapse characteristics of particle column considering gradation and shape[J]. JOURNAL OF MECHANICAL ENGINEERING, 2021, 43(12): 2230-2239. doi: 10.11779/CJGE202112009

考虑颗粒级配和形态的颗粒柱坍塌特性离散元模拟

doi: 10.11779/CJGE202112009
基金项目: 

国家自然科学基金项目 U1765106

国家自然科学基金项目 52079092

详细信息
    作者简介:

    崔溦(1977— ),男,博士,教授,主要从事岩土工程和水工结构物静动力分析研究。E-mail:cuiwei@tju.edu.cn

  • 中图分类号: TU43

Discrete element simulation of collapse characteristics of particle column considering gradation and shape

  • 摘要: 颗粒形态和级配情况等是影响碎屑颗粒流(如滑坡、泥石流、岩崩等)运动的重要因素。基于沃洛诺伊镶嵌原理的随机生成方法创建了不同长细比不同级配的多面体颗粒,引入势粒子算法用于考虑颗粒间的接触作用,根据室内试验确定了离散元接触模型的各项参数,对考虑级配和形态的颗粒柱坍塌特性开展数值试验,研究结果表明:①颗粒柱的归一化堆积高度随颗粒的长细比和中值粒径d50的减小而减小,归一化跑出距离则随其减小而增加。②堆积过程中不同工况的相对静止角α为61.49°~64.99°,且变化规律与归一化堆积高度变化一致。③不同工况的归一化能量耗散为27.1%~35.5%,且转动动能仅占平动动能的8.20%~9.05%。④归一化动能和颗粒配位数呈现负相关的关系,归一化动能达到峰值时颗粒配位数也达到最小值。⑤崩塌过程中强力链一般分布在滑动堆积体的中下部区域,形成力链传递的“拱效应”。中值粒径d50和长细比增大会减少强力链的数量,接触力传递的路径少而集中,从而限制颗粒在堆积过程中的运动。

     

  • 图  随机形态多面生成方法[13]

    Figure  1.  Random shape multi-faceted generation method

    图  势粒子算法

    Figure  2.  Algorithm of potential particles

    图  室内试验设备布置

    Figure  3.  Layout of indoor test equipments

    图  接触摩擦系数测量装置

    Figure  4.  Measuring devices for friction coefficient contact

    图  试验模型

    Figure  5.  Experimental model

    图  物理试验和数值试验对比图

    Figure  6.  Comparison between physical and numerical tests

    图  滑动体积占比和相对静止角

    Figure  7.  Sliding volume ratios and relative angles of repose

    图  颗粒粒径曲线(PSD)

    Figure  8.  Curves of particle-size distribution(PSD)

    图  不同长细比不同粒径的颗粒建模

    Figure  9.  Modeling of particles with different aspect ratios and particle sizes

    图  10  不同工况颗粒柱坍塌过程

    Figure  10.  Collapse process of particle column under different working conditions

    图  11  不同长细比不同级配颗粒柱最终堆积形态对比

    Figure  11.  Comparison of final packing morphologies of particle columns with different aspect ratios and gradations

    图  12  相对静止角

    Figure  12.  Relative rest angles

    图  13  不同工况的相对静止角

    Figure  13.  Angles of repose under different working conditions

    图  14  PSD3时不同长细比下能量变化图

    Figure  14.  Variation of energy of PSD3 with different aspect ratios

    图  15  末时刻不同工况归一化能量耗散图

    Figure  15.  Diagram of normalized energy dissipation under different working conditions at the end

    图  16  PSD3时不同长细比下平动动能和旋转动能随标准时间的变化

    Figure  16.  Variation of translational kinetic energy and rotational kinetic energy with standard time under different aspect ratios in PSD3

    图  17  PSD3时不同长细比颗粒柱配位数、归一化动能曲线

    Figure  17.  Curves of coordination number and normalized kinetic energy of particle column with different aspect ratios at PSD3

    图  18  不同工况末时刻配位数变化曲线

    Figure  18.  Curves of coordination number under different working conditions at the end

    图  19  不同工况强力链占比

    Figure  19.  Proportions of strong chains under different working conditions

    图  20  末时刻不同工况强力链占比

    Figure  20.  Proportions of strong chains under different working conditions at the end

    图  21  PSD3、AR=1时颗粒柱的力链网络图

    Figure  21.  Diagram of force chain network of granular column with different aspect ratios at PSD3 and AR=1

    图  22  末时刻不同工况平均接触力

    Figure  22.  Average contact forces under different conditions at the end

    表  1  离散元参数表

    Table  1.   Parameters of discrete elements

    类别密度/(kg·m³)摩擦角/(°)法向刚度Kn /(108 N·m-1)切向刚度Ks /(108N·m-1)局部阻尼系数
    十二面体250030110.4
    长方体槽350035110.4
    下载: 导出CSV

    表  2  摩擦系数表

    Table  2.   Friction coefficients

    类别摩擦系数
    颗粒-挡板0.24
    颗粒-底板0.25
    下载: 导出CSV

    表  3  室内试验和数值模拟结果对比

    Table  3.   Comparison between indoor tests and numerical simulations

    校准参数最大跑出距离Lf/mm最终堆积高度Hf/mm相对静止角α/(°)滑动体积占比w/%
    室内试验78827447.6256.97
    数值试验77127047.2053.11
    相对误差/%2.161.460.886.7
    下载: 导出CSV

    表  4  不同级配颗粒的中值粒径

    Table  4.   Median particle sizes of different graded particles

    级配PSD1PSD2PSD3
    中值粒径d50/mm27.7519.5215.49
    下载: 导出CSV
  • [1] 刘广煜, 徐文杰, 佟彬, 等. 基于块体离散元的高速远程滑坡灾害动力学研究[J]. 岩石力学与工程学报, 2019, 38(8): 1557-1566.

    LIU Guang-yu, XU Wen-jie, TONG Bin, et al. Study on dynamics of high-speed and long Run-out landslide hazards based on block discrete element method[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(8): 1557-1566. (in Chinese)
    [2] MARKS B, ROGNON P, EINAV I. Grainsize dynamics of polydisperse granular segregation down inclined planes[J]. Journal of Fluid Mechanics, 2012, 690: 499-511.
    [3] 张雪, 盛岱超. 一种模拟土体流动的连续体数值方法[J]. 岩土工程学报, 2016, 38(3): 562-569.

    ZHANG Xue, SHENG Dai-chao. Continuum approach for modelling soil flow in geotechnical engineering[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(3): 562-569. (in Chinese)
    [4] OREFICE L, KHINAST J G. Deformable and breakable DEM particle clusters for modelling compression of plastic and brittle porous materials—Model and structure properties[J]. Powder Technology, 2020, 368: 90-104.
    [5] ZENIT R. Computer simulations of the collapse of a granular column[J]. Physics of Fluids, 2005, 17(3): 31703.
    [6] LUBE G, HUPPERT H E, SPARKS R S J, et al. Collapses of two-dimensional granular columns[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2005, 72(4 Pt 1): 041301.
    [7] UTILI S, ZHAO T, HOULSBY G T. 3D DEM investigation of granular column collapse: Evaluation of debris motion and its destructive power[J]. Engineering Geology, 2015, 186: 3-16.
    [8] PHILLIPS J C, HOGG A J, KERSWELL R R, et al. Enhanced mobility of granular mixtures of fine and coarse particles[J]. Earth and Planetary Science Letters, 2006, 246(3/4): 466-480.
    [9] 张成功, 尹振宇, 吴则祥, 等. 颗粒形状对粒状材料圆柱塌落影响的三维离散元模拟[J]. 岩土力学, 2019, 40(3): 1197-1203.

    ZHANG Cheng-gong, YIN Zhen-yu, WU Ze-xiang, et al. Three-dimensional discrete element simulation of influence of particle shape on granular column collapse[J]. Rock and Soil Mechanics, 2019, 40(3): 1197-1203. (in Chinese)
    [10] ŠMILAUER V. Yade Documentation[M]. 2nd ed. 2015. The Yade Project. doi: 10.5281/zenodo.34073 (http://yade-dem.org/doc/)
    [11] CUNDALL P A, STRACK O D L. Discussion: a discrete numerical model for granular assemblies[J]. Géotechnique, 1980, 30(3): 33-336.
    [12] LANDAUER J, KUHN M, NASATO D S, et al. Particle shape matters - Using 3D printed particles to investigate fundamental particle and packing properties[J]. Powder Technology, 2020, 361: 711-718.
    [13] ELIÁŠ J. Simulation of railway ballast using crushable polyhedral particles[J]. Powder Technology, 2014, 264: 458-465.
    [14] BOON C W, HOULSBY G T, UTILI S. A new algorithm for contact detection between convex polygonal and polyhedral particles in the discrete element method[J]. Computers and Geotechnics, 2012, 44: 73-82.
    [15] BOYD S, VANDENBERGHE L. Convex Optimization[M]. Cambridge: Cambridge University Press, 2004.
    [16] COETZEE C J. Review: Calibration of the discrete element method[J]. Powder Technology, 2017, 310: 104-142.
    [17] ZHAO S W, ZHOU X W, LIU W H. Discrete element simulations of direct shear tests with particle angularity effect[J]. Granular Matter, 2015, 17(6): 793-806.
    [18] MINDLIN R D. Compliance of elastic bodies in contact[J]. Journal of Applied Mechanics, 1949, 16(3): 259-268.
    [19] LI Y J, XU Y, THORNTON C. A comparison of discrete element simulations and experiments for ‘sandpiles’ composed of spherical particles[J]. Powder Technology, 2005, 160(3): 219-228.
    [20] 王玉峰, 程谦恭, 朱圻. 汶川地震触发高速远程滑坡-碎屑流堆积反粒序特征及机制分析[J]. 岩石力学与工程学报, 2012, 31(6): 1089-1106.

    WANG Yu-feng, CHENG Qian-gong, ZHU Qi. Inverse grading analysis of deposit from rock avalanches triggered by Wenchuan earthquake[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(6): 1089-1106. (in Chinese)
    [21] 边学成, 李伟, 李公羽, 等. 基于颗粒真实几何形状的铁路道砟剪切过程三维离散元分析[J]. 工程力学, 2015, 32(5): 64-75, 83.

    BIAN Xue-cheng, LI Wei, LI Gong-yu, et al. Three-dimensional discrete element analysis of railway ballast's shear process based on particles' real geometry[J]. Engineering Mechanics, 2015, 32(5): 64-75, 83. (in Chinese)
    [22] 杨舒涵, 周伟, 马刚, 等. 粒间摩擦对岩土颗粒材料三维力学行为的影响机制[J]. 岩土工程学报, 2020, 42(10): 1885-1893.

    YANG Shu-han, ZHOU Wei, MA Gang, et al. Mechanism of inter-particle friction effect on 3D mechanical response of granular materials[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1885-1893. (in Chinese)
    [23] ZHAO X L, EVANS T M. Numerical analysis of critical state behaviors of granular soils under different loading conditions[J]. Granular Matter, 2011, 13(6): 751-764.
    [24] 孙其诚, 王光谦. 静态堆积颗粒中的力链分布[J]. 物理学报, 2008, 57(8): 4667-4674.

    SUN Qi-cheng, WANG Guang-qian. Force distribution in static granular matter in two dimensions[J]. Acta Physica Sinica, 2008, 57(8): 4667-4674. (in Chinese)
    [25] 戴北冰, 杨峻, 刘锋涛, 等. 散粒土自然堆积的宏细观特征与形成机制[J]. 岩土工程学报, 2019, 41(增刊2): 57-60.

    DAI Bei-bing, YANG Jun, LIU Feng-tao, et al. Macro-and micro-properties and formation mechanisms of granular piles[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S2): 57-60. (in Chinese)
  • 加载中
图(22) / 表(4)
计量
  • 文章访问数:  93
  • HTML全文浏览量:  91
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-23
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-12-01

目录

    /

    返回文章
    返回