›› 2015, Vol. 36 ›› Issue (7): 2081-2087.doi: 10.16285/j.rsm.2015.07.033

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

基于距离和局部Delaunay三角化控制的颗粒离散元模型填充方法研究

王秀菊1,石 崇2,李德杰2,梁邦炎3,严晨宇3   

  1. 1. 南京交通职业技术学院 建筑工程系,江苏 南京 211188;2. 河海大学 岩土工程科学研究所,江苏 南京 210098; 3. 中交第四航务工程局有限公司,广东 广州 510000
  • 收稿日期:2014-02-17 出版日期:2015-07-11 发布日期:2018-06-13
  • 通讯作者: 石崇,男,1978年生,博士,副教授,主要从事岩石力学与工程方面的教学与研究工作。E-mail: scvictory@hhu.edu.cn E-mail:wangxiujuhf@126.com
  • 作者简介:王秀菊,女,1981年生,硕士,讲师,主要从事岩土数值仿真与工程应用方面的研究与教学工作。
  • 基金资助:
    国家自然科学基金青年项目(No. 51309089,No. 11202063);江苏省基础研究计划(自然科学基金)(No. BK20130846)。

A model construction method for particle DEM with distance control and local Delaunay triangulation

WANG Xiu-ju1,SHI Chong2,LI De-jie2,LIANG Bang-yan3,YAN Chen-yu3   

  1. 1. Department of Architectural Engineering, Nanjing Communication Institute of Technology, Nanjing, Jiangsu 211188, China; 2. Geotechnical Research Institute, Hohai University, Nanjing, Jiangsu 210098, China; 3. CCCC Fourth Harbor Engineering Co., Ltd., Guangzhou, Guangdong 510000, China
  • Received:2014-02-17 Online:2015-07-11 Published:2018-06-13

摘要: 根据颗粒离散元Kelvin 接触力计算模型,分析了圆形颗粒体模拟材料力学特性应具备的条件,在此基础上提出了一种新颗粒模型构建方法。该方法首先在复杂模型域内随机生成种子,然后利用相切条件逐步扩展填充整个区域。填充过程中借助局部Delaunay三角化网格控制新颗粒的生成,采用复杂几何体距离控制颗粒与模型边界的相对位置,对靠近模型边界的颗粒进行容忍性优化填充,从而增加模型颗粒与边界的耦合性。同时对模型孔隙进行再填充,保证每个填充颗粒至少与3个颗粒相切,提高了模型内颗粒间的耦合性和模型的密度。最后采用任意多边形控制材料边界,将模型材料的设置简化为判断点是否在多边形内,简化了复杂模型材料属性的设置过程。结果表明:与膨胀颗粒生成法相比,该方法生成模型重叠量小、颗粒间及颗粒-边界相互耦合、填充率高。因此,颗粒黏结力破坏后不会造成飞溢现象,可适用于任意连通域模型的生成,能更好地实现复杂岩土细观介质变形破坏机制的模拟与研究。

关键词: 离散颗粒, 复杂边界, 局部Delaunay三角化, 重填, 边界耦合

Abstract: Based on Kelvin’s contact force calculation model, the conditions required for simulating the mechanical behavior of the porous materials consisting of spherical particles are analyzed, from which a new granular model construction method is proposed. In this model a seed position is stochastically generated in the complex solution domain; and then it is gradually expanded to fill the entire region. In the filling process a local Delaunay triangulation mesh control is used to generate new particles; and particle distance control in a complex geometry boundary is used to decide the relative position of a new generated particle. For particles near the model boundary, an optimized tolerance position and radius are calculated and filled, which enhances the coupling characteristic between boundary and the particles. Meanwhile, the pore position can be refilled to ensure that each particle is tangent with at least three particles in order to improve the coupling characteristic between model particles. Finally, arbitrary polygons are used to divide several materials. The model region will be simplified for judging whether the polygon point is in the polygon or not, which can simplify the setup process of complex model materials. The results show that the proposed method has small overlap between particles, high filling rate between particles, particle and boundary. The proposed method can decrease flying overflow phenomenon when the adhesive force between particles disappears, and can be used to model construction of random connected domain.

Key words: discrete particles, complicate boundary, local Delaunay triangulation, refilling, boundary coupling

中图分类号: 

  • O 241
[1] 李 毅 ,伍 嘉 ,李 坤. 基于FLAC3D的饱和-非饱和渗流分析[J]. , 2012, 33(2): 617-622.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 肖云华,王 清,陈剑平. 基于优化技术的权重计算方法在岩体质量评价中的应用[J]. , 2009, 30(9): 2686 -2690 .
[2] 范庆来,栾茂田,刘占阁. 软土中T型触探仪贯入阻力的数值模拟[J]. , 2009, 30(9): 2850 -2854 .
[3] 张安康,陈士海,杜荣强,魏海霞. 岩石类材料的能量基率相关弹塑性损伤模型[J]. , 2010, 31(S1): 207 -210 .
[4] 王小军,屈耀辉,魏永梁,杨印海,达益正. 郑西客运专线湿陷性黄土区试验路堤的沉降观测与预测研究[J]. , 2010, 31(S1): 220 -231 .
[5] 陈 瑜,曹 平,蒲成志,刘业科,李 娜. 水-岩作用对岩石表面微观形貌影响的试验研究[J]. , 2010, 31(11): 3452 -3458 .
[6] 赵延喜,徐卫亚. 基于AHP和模糊综合评判的TBM施工风险评估[J]. , 2009, 30(3): 793 -798 .
[7] 张其一,栾茂田. 复合加载情况下非均质地基上条形基础的极限承载力研究[J]. , 2009, 30(5): 1281 -1286 .
[8] 王俊卿,李 靖,李 琦,陈 立. 黄土高边坡稳定性影响因素分析 ——以宝鸡峡引水工程为例[J]. , 2009, 30(7): 2114 -2118 .
[9] 常林越,王金昌,朱向荣. 多级线性荷载下饱和软黏土一维大应变固结解析解[J]. , 2009, 30(8): 2343 -2347 .
[10] 龚彦峰,张俊儒. 隧道单层衬砌设计方法研究及应用[J]. , 2011, 32(4): 1062 -1068 .