Le comportement mécanique des matériaux granulaires, et en particulier des poudres métalliques ductiles, peut être appréhendé par l'intermédiaire de simulations discrètes dans lesquelles les grains sont modélisés. Cette méthode a l'avantage de permettre l'exploration de chemins de chargement inaccessibles à une approche expérimentale. La méthode développée dans cette thèse consiste à modéliser des assemblages de sphères maillées dans le logiciel d'éléments finis ABAQUS pour effectuer des simulations en appliquant divers chemins de chargement. La première partie de ce travail consiste à utiliser cette méthode pour enrichir la méthode des éléments discrets en l'étendant au domaine de la compression à haute densité relative. La méthode de simulations discrètes par éléments finis donne des informations précieuses sur les déformations non-linéaires des grains aux contacts. Ces informations permettent de proposer un modèle de contact qui gère les interactions complexes entre les diérents contacts à haute densité relative. Par la suite, ce modèle est introduit dans le code d'éléments discrets YADE, et les résultats des deux méthodes peuvent être comparés. Dans la seconde partie de cette thèse, un échantillon numérique est soumis à des sollicitations dans toutes les directions de l'espace des contraintes, pour obtenir les surfaces de charge macroscopiques. L'accent est mis sur l'évolution des surfaces de charge lors de chargements complexes, avec des changements de direction, afin d'obtenir des indications sur la mémoire du matériau. Enfin, l'observation de la microstructure de l'échantillon permet de comprendre les phénomènes micromécaniques qui sont à l'origine de l'évolution macroscopique des surfaces de charge.