added 3rd order taylor-eytension based integrator for DEM model, and a second model for normal force