Numerical simulation of density distribution during compaction of iron powders