Mechanics of binary and polydisperse spherical pebble assembly

The micromechanical behavior of an assembly of binary and polydisperse spherical pebbles is studied using discrete element method (DEM) accounting for microscopic interactions between individual pebbles. A in-house DEM code has been used to simulate the assemblies consisting of different pebble diameters and the results of the simulations are compared with that of mono-size pebble assemblies. The effect of relative radii and volume fraction of the pebbles on the macroscopic stress–strain response is discussed. Furthermore, the effect of packing factor and coefficient of friction on the overall stress–strain behavior of the system is studied in detail. The shear (tangential) stiffness between the particles is also another influencing parameter. For a very small shear stiffness the system shows a strong dependence on the packing factor while a pebble material dependent shear stiffness shows a rather moderate dependence on the packing factor. For a similar packing factor, the mono-size assembly shows a stiff behavior during loading compared to binary assembly. However, the simulations do not show a significant difference between the two behaviors in contrast to the observations made in the experiments. The discrepancy can be attributed to (i) probable difference in packing factors for mono-size and binary assemblies in the experiments, (ii) arbitrary friction coefficient in the current model and (iii) the tangential interaction (constant shear stiffness) implemented in the present model which needs further modification as a function of the load history on the pebbles. Evolution of other micromechanical characteristics such as coordination number, contact force distribution and stored elastic energy of individual pebbles as a function of external load and system parameters is presented which can be used to estimate important macroscopic properties such as overall thermal conductivity and crushing resistance of the pebble beds.