We present a three-dimensional distinct element model (DEM) able to handle populations of spherocylinders. We report on granular crystallization occurring when vibrating mono-disperse assemblies of spherocylinders that faithfully reproduce the corresponding results of physical experiments from the literature.