Is there any computing module which can compute a heap of disk elements in DEM(discrete element method)?

Is there any computing module which can compute a heap of disk elements in DEM(discrete element method)?it'd be better if it is written in matlab or C/C++. I'm doing research on the combined DEM-FEM. I need such module so as to save my time.

Thank you