fixed some errors in force calculation for DEM model