基于CT扫描的不规则外形颗粒三维离散元建模

上海交通大学学报（自然版） ›› 2011, Vol. 45 ›› Issue (05): 711-715.

基于CT扫描的不规则外形颗粒三维离散元建模

杜欣，曾亚武，高睿，颜敬

(武汉大学土木建筑工程学院，湖北武汉 430072)

收稿日期:2010-10-19 出版日期:2011-05-30 发布日期:2011-05-30
基金资助:
中央高校基本科研业务费专项资金资助(20102100101000093），湖北省自然科学基金重点项目(2010CDA057)

3-D Modelling of Irregular Shape Particles for Discrete Element Method Based on X-ray Tomography

DU Xin, ZENG Ya-Wu, GAO Rui, YAN Jing

(School of Civil Engineering, Wuhan University, Wuhan 430072, China)

Received:2010-10-19 Online:2011-05-30 Published:2011-05-30

摘要/Abstract

摘要： 为实现不规则颗粒材料三维离散元建模，通过X射线断层扫描技术获取真实颗粒的切片图像，并对其进行二值化处理和边缘点探测提取计算，以重建三维颗粒模型而获取颗粒空间坐标信息.引入最小基本球体半径、有效空间系数和分类半径作为控制参数，采用VC++语言编写颗粒模型并生成程序，运用基本球体进行可重叠组合，自动生成不规则颗粒模型.结果表明，在相同参数条件下，该程序运行结果惟一，且与初始条件无关.同时，通过调控控制参数，能够获得不同精度的模型，并使颗粒模型精度与所需球体之间达到平衡关系.

关键词: X射线断层扫描, 不规则外形颗粒, 离散元方法, 数值建模

Abstract: In order to build 3-D models of irregular shape particles for discrete element method, the real ballast was scanned with X-ray tomography and its tomography images were obtained, a binary image processing was implemented for each image, and with the method of edgepoints detection, the three dimensional coordinates of the ballast were collected. Then, three governing parameters: minimum radius of basicspheres, efficient space factor, and splitradius were introduced and the generating program of particle was compiled based on VC++. The program combines the basicspheres automatically under an overlap-able algorithm to build the model of irregular shape particles. The research results indicate that with the same preset parameters, the process of modelling can be repeated, and the program is not related to the initial conditions. Meanwhile, the balance of model accuracy and required ball number can be achieved by adjusting governing parameters.

Key words: X-ray tomography (CT), irregular particles, discrete element method (DEM), numerical modelling

中图分类号:

TU441.93

杜欣, 曾亚武, 高睿, 颜敬. 基于CT扫描的不规则外形颗粒三维离散元建模[J]. 上海交通大学学报（自然版）, 2011, 45(05): 711-715.

DU Xin, ZENG Ya-Wu, GAO Rui, YAN Jing. 3-D Modelling of Irregular Shape Particles for Discrete Element Method Based on X-ray Tomography[J]. Journal of Shanghai Jiaotong University, 2011, 45(05): 711-715.

参考文献

［1］Coetzee C J, Els D N J. Calibration of granular material parameters for DEM modeling and numerical verification by bladegranular material interaction［J］. Journal of Terramechanics, 2009, 46(1): 1526.

［2］Taylor M A, Garboczi E J, Erdogan S T, et al. Some properties of irregular 3D particles［J］. Powder Technology, 2006, 162(1): 115.

［3］Cundall P A, Strack O D L. A discrete element model for granular assemblies［J］. Geotechnique, 1979, 29(1):4765.

［4］徐泳, 孙其诚, 张凌, 等. 颗粒离散元法研究进展［J］. 力学进展, 2003, 33(2): 251260.

XU Yong, SUN Qicheng, ZHANG Ling, et al. Advances in discrete element methods for particulate materials［J］. Advances in Mechanics, 2003, 33(2): 251260.

［5］王安麟, 邵萌, 梁波, 等. 三维流固两相流的颗粒群轨道柔性模型［J］. 上海交通大学学报, 2005, 39(5):790794.

WANG Anlin, SHAO Meng, LIANG Bo, et al. Particles trajectory flexible model of 3D fluidsolids twophase flows［J］. Journal of Shanghai Jiaotong University, 2005, 39(5): 790794.

［6］Lu M, McDowell G R. The importance of modelling ballast particle shape in the discrete element method［J］. Granular Matter, 2007, 9(12):6980.

［7］Wang L B, Park J Y, Fu Y. Representation of real particles for DEM simulation using Xray tomography［J］. Construction and Building Materials, 2007, 21(2): 338346.

［8］Wang L B, Frost J D, Lai J S. Threedimensional digital representation of granular material microstructure from Xray tomography imaging ［J］. Journal of Computing in Civil Engineering, 2004, 18(1): 2835.

［9］Garcia X, Latham J P, Xiang J, et al. A clustered overlapping sphere algorithm to represent real particles in discrete element modeling［J］. Geotechnique, 2009, 59(9): 779784.

［10］Ferellec J F, McDowell G R. A simple method to create complex particle shapes for DEM［J］. Geomechanics and Geoengineering, 2008, 3(3): 211216.

［11］巍柏, 彭刚, 戚永乐, 等. 基于CT图像的混凝土细观结构的有限元重建［J］. 混凝土，2008, 30(2): 6465.

WEI Bai, PENG Gang, QI Yongle, et al. Reconstructing finite element model of concrete mesolevel structure with CT scanning image［J］. Concrete, 2008, 30(2): 6465.