Ce mémoire de thèse synthétise trois années de recherches dédiées à l'étude numérique et théorique de la compression à haute densité de poudres métalliques. Des diérentes phases qu'intègrent la métallurgie des poudres, la phase de compression à froid de la poudre est l'une des phases les plus sensibles de ce procédé de fabrication, car elle influence les propriétés mécaniques de la pièce finale. Il est donc nécessaire de mettre en place une approche numérique qui permet de contrôler et d'optimiser la compression de poudre jusqu'à de fortes valeurs de compacité (compacité supérieure à 0:9). Pour cela, nous proposons de reproduire par la méthode des éléments discrets le comportement de la poudre observé expérimentalement sous diérents types de chargement. A ce jour, les simulations via cette méthode sont limitées à une valeur de compacité ne dépassant pas 0:85. Pour dépasser ces limitations, nous présentons un modèle de contact implémenté dans un code éléments discrets libre (Yade). Ce nouveau modèle de contact est développé sur la base de la loi de contact normal qui intègre le terme de densité locale des particules dans son expression, afin de prendre en compte l'incompressibilité des grains se produisant à des valeurs de compacité supérieures à 0:85. Dans le but de procéder à des simulations plus réalistes, un nouvel algorithme géométrique de génération d'empilements de sphères polydisperses est développé. Ce nouvel outil numérique est capable de générer très rapidement de grands assemblages de sphères en contact tout en contrôlant diérents paramètres comme la distribution de la compacité, la taille minimale et maximale des sphères. Avec le modèle de contact capable de reproduire l'interaction entre les grains et la création d'un algorithme pouvant générer des assemblages de sphères similaires à un tas de poudres, nous procédons à des simulations de compression isostatique et en matrice pour diérents types de poudres (cuivre, aluminium, fer). Les résultats obtenus sont directement comparés à ceux issus des simulations éléments finis multi-particules et de l'expérience. Ces comparaisons permettent ainsi de valider et de tester la robustesse du modèle de contact développé. Pour finir, nous investiguons sur la base de nos divers développements validés, l'évolution d'une poudre d'aluminium avec un gradient de compacité au cours d'une compression en matrice.