Numerical Simulation of Mechanical Properties and Force Chain Evolution of Granular System during Densification Process

doi:10.12382/bgxb.2025.0690

Abstract

Abstract:

To explore the mechanical properties and microscopic mechanisms of particle system densification process，a measurement ball generation method with automatic position adjustment is proposed and a uniaxial compression discrete element numerical model of particle system is established.The evolution laws of displacement field and velocity field of particle system during densification process are analyzed，and the influence laws of particle size and loading speed on porosity and coordination number of particle system are studied.The evolution characteristics of force chain network of particle system during densification process are analyzed，and the influence laws and microscopic mechanisms of mechanical properties of particle system during densification process are revealed.The research results indicate that the densification process of particle systems can be divided into three stages：rapid densification，slow densification，and steady densification.A wider particle size distribution is conducive to improving the filling density of particles and accelerating the densification process.During the densification process，the displacement field and velocity field exhibit a centralization phenomena.The loading speed has a significant impact on porosity The porosity decreases from 35.23% to 27.50%，26.2%，24.4%，and 22.3%，respectively，at speeds of 0.05m/s，0.5m/s，1.5m/s，and 3m/s.Increasing the loading speed is beneficial for the compaction of the particle system.The larger the particle size is，the greater the particle deformation and contact strength are during the loading process.The combination of 0.8-1.5mm mixed particle sizes significantly affects the formation of force chains and the uniformity of stress，leading to a stable state of short chain dominance and long chain shrinkage during the densification process of particle system.

Key words: particle system, mechanical property, force chain evolution, discrete element

YU Yang, CHEN Fanxiu, SHEN Jun, LIU Rui, ZHOU Jinqiang, WU Chengcheng. Numerical Simulation of Mechanical Properties and Force Chain Evolution of Granular System during Densification Process[J]. Acta Armamentarii, 2025, 46(S1): 250690-.

Figures/Tables 15

References 23

[1]	OUYANG Q M, KANG G, ZHUANG X Y, et al. Development of crack surface tracking algorithm for explosive fracture simulation with three-dimensional numerical manifold method[J]. Engineering Fracture Mechanics, 2025, 313:110645.
[2]	REDDY N S C, HE H, SENETAKIS K. Multi-scale mechanics of polydisperse granular materials:From micro-scale and wave propagation experiments to DEM analysis[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2024, 48(14):3515-3536.
[3]	NIE J Y, CUI Y F, WANG G D, et al. A comprehensive numerical investigation of multi-scale particle shape effects on small-strain stiffness of sands[J]. Géotechnique, 2024, 75(3):323-336.
[4]	孙其诚, 刘晓星, 张国华, 等. 密集颗粒物质的介观结构[J]. 力学进展, 2017, 47(1):263-308.
	SUN Q C, LIU X X, ZHANG G H, et al. Mesoscopic structure of dense granular matter[J]. 2017, 47(1):263-308. (in Chinese)
[5]	YU Y, CHEN F X, GU Y J, et al. Optimization of 3D reconstruction of granular systems based on refractive index matching scanning[J]. Optics & Laser Technology, 2025, 186:112662.
[6]	张催, 张才鑫, 杨亚飞, 等. 造型颗粒统计特征参量CT量化表征方法[J]. 含能材料, 2023, 31(12):1262-1268.
	ZHANG C, ZHANG C X, YANG Y F, et al. CT quantitative characterization method of statistical characteristic parameters of modeling particles[J]. Energetic Materials, 2023, 31(12):1262-1268. (in Chinese)
[7]	ROGERS B M, HARDY A J, COFFEE K R, et al. Segmentation of RDX and TNT in X-ray computed tomography reconstructions of melt-cast explosives[J]. Propellants,Explosives,Pyrotechnics, 2025, 50(2):e202400065.
[8]	GU Y J, CHEN F X, YU Y, et al. Structural properties and mechanical behavior of three-dimensional cylindrical particle-like systems under in situ loading[J]. Powder Technology, 2025, 449:120365.
[9]	MOR Y, SAKAI M. Advanced DEM simulation on powder mixing for ellipsoidal particles in an industrial mixer[J]. Chemical Engineering Journal, 2022, 429:132415.
[10]	GONG J, PANG X W, Tang Y, et al. Effects of particle shape,physical properties and particle size distribution on the small-strain stiffness of granular materials:a DEM study[J]. Computers and Geotechnics, 2024, 165:105903.
[11]	KETTERHAGEN W, CARL W. A perspective on calibration and application of DEM models for simulation of industrial bulk powder processes[J]. Powder Technology, 2022, 402:117301.
[12]	CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1):47-65.
[13]	ZHANG M, KONIETZKY H, SONG Z Y. Evolution of force network,contact network,and tensile force chains in rock-like bonded granular materials under unconfined and confined compression:a DEM study[J]. Environmental Earth Sciences, 2025, 84(12):320.
[14]	HORABIK J, MOLENDA M. Parameters and contact models for DEM simulations of agricultural granular materials:a review[J]. Biosystems Engineering, 2016, 147(2):206-225.
[15]	THAN V D, KHAMSEH S, TANG A M, et al. Basic mechanical properties of wet granular materials:a DEM study[J]. Journal of Engineering Mechanics, 2017, 143(1):C4016001.
[16]	TAHA H, NGUYEN N S, MAROT D, et al. Micro-scale investigation of the role of finer grains in the behavior of bidisperse granular materials[J]. Granular Matter, 2019, 21(2):28.
[17]	HU J, WU H W, GU X Q, et al. Particle shape effects on dynamic properties of granular soils:a DEM study[J]. Computers and Geotechnics, 2023, 161:105578.
[18]	ADESINA P, O’SULLIVAN C, WANG T. DEM study on the effect of particle shape on the shear behaviour of granular materials[J]. Computational Particle Mechanics, 2024, 11(1):447-466.
[19]	ZHOU W, YANG L F, MA G, et al. DEM analysis of the size effects on the behavior of crushable granular materials[J]. Granular Matter, 2016, 18:1-11.
[20]	李旭, 刘彦, 安丰江, 等. 行星式球磨颗粒流场分布与形貌变化规律[J]. 兵工学报, 2022, 43(4):876-891. doi: 10.12382/bgxb.2021.0252
	LI X, LIU Y, AN F J, et al. Flow field distribution and morphology variation of particles in planetary ball milling[J]. Acta Armamentarii, 2022, 43(4):876-891. (in Chinese) doi: 10.12382/bgxb.2021.0252
[21]	XIAO Y C, GONG T Y, ZHANG X W, et al. Multiscale modeling for dynamic compressive behavior of polymer bonded explosives[J]. International Journal of Mechanical Sciences, 2023, 242:108007.
[22]	SUN H T, FAN X G, ZHAN M, et al. Mechanism of uneven densification in PBX compression molding[J]. International Journal of Mechanical Sciences, 2024, 282:109683.
[23]	MA H Q, ZHOU L Y, LIU Z H, et al. A review of recent development for the CFD-DEM investigations of non-spherical particles[J]. Powder Technology, 2022, 412:117972.

参数	数值
墙体法向刚度k_n/（N·m^-1）	2×10¹²
墙体切向刚度k_s /（N·m^-1）	2×10¹²
墙体摩擦系数μ₁	0.25
阻尼系数	0.7
颗粒密度ρ/（kg·m^-3）	1190
颗粒摩擦系数μ₂	0.25
泊松比ν	0.3
弹性模量E/MPa	109
颗粒粒径R/mm	0.8~1.5
颗粒数量	2×10⁴~6×10⁴

参数	数值
墙体法向刚度k_n/（N·m^-1）	2×10¹²
墙体切向刚度k_s /（N·m^-1）	2×10¹²
墙体摩擦系数μ₁	0.25
阻尼系数	0.7
颗粒密度ρ/（kg·m^-3）	1190
颗粒摩擦系数μ₂	0.25
泊松比ν	0.3
弹性模量E/MPa	109
颗粒粒径R/mm	0.8~1.5
颗粒数量	2×10⁴~6×10⁴

加载速度/ （m·s^-1）	粒径范围/mm
加载速度/ （m·s^-1）	0.8~1.0	1.0~1.2	1.2~1.5	颗粒数量
0.05	100	0	0	66477
0.5	0	100	0	43264
1.5	0	0	100	25796
3	50	25	25	55153

加载速度/ （m·s^-1）	粒径范围/mm
加载速度/ （m·s^-1）	0.8~1.0	1.0~1.2	1.2~1.5	颗粒数量
0.05	100	0	0	66477
0.5	0	100	0	43264
1.5	0	0	100	25796
3	50	25	25	55153

粒径范围/mm	加载速度/ （m·s^-1）	初始状态	快速压密阶段	缓慢压密阶段	平稳阶段
0.8~1.0	0.05	35.23	27.50	24.52	23.16
	0.5	35.23	26.22	22.31	21.48
	1.5	35.23	24.41	20.82	18.82
	3	35.23	22.35	18.33	17.64
1.0~1.2	0.05	35.33	27.50	24.22	23.16
	0.5	35.33	26.82	23.31	21.95
	1.5	35.33	25.41	21.80	19.88
	3	35.33	22.35	19.33	18.59
1.2~1.5	0.05	35.54	31.51	28.01	25.94
	0.5	35.54	29.52	26.06	24.20
	1.5	35.54	27.02	22.80	21.70
	3	35.54	23.35	20.86	20.32
0.8~1.5	0.05	36.51	27.90	24.74	24.1
	0.5	36.51	26.42	23.55	22.52
	1.5	36.51	24.80	20.08	18.82
	3	36.51	22.02	18.02	16.41