3D mode discrete element method with the elastoplastic model