Mesoscopic mechanics and microstructural evolution of triaxially compressed polydisperse pebbles: insights from DEM modeling

DK Pawar and RK Annabattula and N Swaminathan, GRANULAR MATTER, 27, 74 (2025).

DOI: 10.1007/s10035-025-01550-7

A computational framework that utilizes the discrete element method (DEM) was developed to conduct triaxial compression tests on lithium- based pebbles. In addition, the process for determining optimized simulation parameters for performing stress-controlled numerical hydrostatic and triaxial compression tests using LIGGGHTS has been outlined. This framework was employed to study the influence of the polydispersity of pebbles on the Drucker-Prager (D-P) parameter, which is the friction angle. The study revealed that the effect of polydispersity, measured by the polydispersity index (\documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda $$\enddocument) was negligible, as the value of friction angle (\documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\beta $$\enddocument) remained constant at \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\approx 29<^>\circ $$\enddocument for \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda \le 0.5$$\enddocument. However, in highly polydisperse samples with \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda =0.8$$\enddocument, \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\beta $$\enddocument increased to \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\approx 31<^>\circ $$\enddocument. Additionally the dilatancy angle (\documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\psi $$\enddocument) decreases as \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda $$\enddocument increases. The difference between \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\beta $$\enddocument and \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\psi $$\enddocument increases with \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda $$\enddocument, and thus, the associative flow rule is not suitable for highly polydisperse samples. To conduct a more detailed analysis of the mesoscopic mechanics, three parameters based on the mean number of contacts, which governed the microstructural similarity of the sample and the extent of particle participation, were examined. Additionally, two local parameters based on Voronoi tessellation were investigated. These parameters highlight how changing \documentclass12ptminimal \usepackageamsmath \usepackagewasysym \usepackageamsfonts \usepackageamssymb \usepackageamsbsy \usepackagemathrsfs \usepackageupgreek \setlength\oddsidemargin-69pt \begindocument$$\lambda $$\enddocument influences the local deformation and characterizes the local structural variation in the granular samples. In particular it was found that, the participation of particles to the total deformation was higher in samples with high polydispersity.

Return to Publications page