The discrete element method (DEM), based on a soft-sphere approach, is commonly used to simulate powder compaction. With these simulations a new macroscopic constitutive relation can be formulated. It is able to de-scribe accurately the constitutive material of powders during the cold compaction process. However, the force-law used in the classical DEM formulation does not reproduce correctly the stress evolution during the high density compaction of powder. To overcome this limitation at a relative density of about 0.85. the high density model is used. This contact model can reproduce incompressibility effects in granular media by implementing the local solid fraction into the DEM software, using Voronoi cells. The first DEM simulations using the open-source YADE software show a fairly good agreement with the multi-particle finite element simulations and experimental results. (C) 2010 Elsevier B.V. All rights reserved.