Simulation of triaxial response of granular materials by modified DEM