La dégradation géométrique de la voie ballastée sous circulation commerciale nécessite des opérations de maintenance fréquentes et onéreuses. La caractérisation du comportement des pro-cédés de maintenance comme le bourrage, la stabilisation dynamique, est nécessaire pour proposer des améliorations en terme de méthode, paramétrage pour augmenter la pérennité des travaux. La simulation numérique d'une portion de voie soumise à un bourrage ou une stabilisation dynamique permet de comprendre les phénomènes physiques mis en jeu dans le ballast. Toutefois, la complexité numérique de ce problème concernant l'étude de systèmes à très grand nombre de grains et en temps de sollicitation long, demande donc une attention particulière pour une résolution à moindre coût. L'objectif de cette thèse est de développer un outil de calcul numérique performant qui permet de réaliser des calculs dédiés à ce grand problème granulaire moins consommateur en temps. La méthodologie utilisée ici se base sur l'approche Non Smooth Contact Dynamic s(NSCD) avec une discrétisation par Éléments Discrets (DEM). Dans ce cadre, une méthode de dé-composition de domaine (DDM) alliée à une parallélisation adaptée en environnement à mémoire partagée utilisant OpenMP sont appliquées pour améliorer l'efficacité de la simulation numérique. / The track deterioration rate is strongly influenced by the ballast behaviour under commercial traffic. In order to restore the initial track geometry, different maintenance processes are performed, like tamping, dynamic stabilisation. A better understanding of the ballast behaviour under these operations on a portion of railway track is a key to optimize the process, to limit degradationand to propose some concept for a better homogeneous compaction. The numerical simulation isdeveloped here to investigate the mechanical behaviour of ballast. However, the main difficultiesof this research action concerns the size of the granular system simulation increasing both in termof number of grains and of process duration. The purpose of this thesis is to develop an efficient numerical tool allows to realize faster computations devoted to large-scale granular samples. In this framework, the Non-Smooth Contact Dynamics (NSCD) of three-dimensional Discrete ElementMethod (DEM) simulations, improved by Domain Decomposition Method (DDM) and processedwith the Shared Memory parallel technique (using OpenMP) has been applied to study the ballast media mechanics.