Discrete element modelling of the quasi-static uniaxial compression of individual infant formula agglomerates

► A 3D DEM was developed for individual agglomerates of infant formula. ► The Taguchi method was used to calibrate the DEM parameters with quasi-static loading data. ► Weibull moduli for the simulation results and the experimental data were similar. ► The energy dissipation rates were highest where bond breakage events occurred.