La fatigue de contact est un des modes de défaillance prédominants des composants tels que les engrenages ou les roulements. Les mécanismes d’initiation de fissures associés à ce mode de défaillance sont fortement liés à la microstructure du matériau. Cependant, la plupart des modèles utilisés pour prédire la durée de vie se situent à l’échelle macroscopique. Un modèle basé sur une représentation de type Voronoi des grains (échelle mésoscopique) est développé afin d’analyser les mécanismes d’initiation. Le concept d’endommagement est appliqué aux joints de grain modélisés par la méthode des zones cohésives. L’objectif de ce modèle est (i) de contribuer à une meilleure compréhension de l’influence de paramètres tels que ceux liés aux conditions de contact (rugosité, lubrification) ou aux matériaux (présence d’inclusions ou gradients de propriétés et contraintes résiduelles générés par les traitements de surface…) sur les mécanismes d’initiation et (ii) de fournir une estimation de la durée de vie jusqu’à cette initiation. Un premier modèle 2D isotrope a permis de mettre en place l’approche proposée et d’analyser le comportement numérique des éléments cohésifs : influence de la valeur des raideurs cohésives et apparition de singularités aux jonctions triples. Cette singularité semble inévitable, mais l’approche consistant à considérer le joint de grain comme une unique entité, et donc à utiliser des valeurs moyennes le long du joint de grain permet de s’affranchir de cette singularité. La représentativité du modèle a ensuite été améliorée par la modélisation de l’anisotropie cristalline. Un modèle de type élasticité cubique a été utilisé pour modéliser le comportement des grains. Enfin, une analyse approfondie de l’application du concept d’endommagement aux joints de grains a permis de proposer une nouvelle formulation entraînant une influence plus réaliste de cet endommagement sur le cisaillement intergranulaire et conduisant à une durée de vie estimée (apparition des premières micro-fissures) d’un ordre de grandeur comparable à celles données par l’expérience. / Contact fatigue is the predominant mode of failure of components subjected to a repeated contact pressure, like rolling element bearings or gears. This phenomenon is known as rolling contact fatigue (RCF). A large number of models have been developed to predict RCF, but there is today no complete predictive life model, and understanding RCF failure mechanism remains a significant challenge. RCF failure mechanisms are known to be very sensitive to a large number of parameters linked to contact conditions (roughness, lubrication) or materials (inclusions, gradients properties, residual stresses…). To improve knowledge about the influence of these parameters on failure mechanisms and life, a numerical model is developed to simulate the progressive damage of a component subject to rolling contact fatigue. Mechanisms associated with the initiation stage of failure process are located at a scale lower than the macroscopic scale. The proposed approach is to develop a grain level model (mesoscopic scale) in order to focus on initiation mechanisms. A Voronoi tessellation is used to represent the material microstructure. The progressive deterioration is simulated by applying the concept of damage mechanics at grain boundaries represented by cohesive elements. This approach has been first applied to a 2D isotropic model. The numerical behaviour of cohesive elements has been investigated: the influence of cohesive stiffness has been analysed and singularities at the triple junctions has been highlighted. The representativeness of the original model was improved by modelling crystal anisotropy. A cubic elasticity model was used to represent the behaviour of grains. Finally, a thorough analysis of the application of the damage concept at grain boundaries highlighted that the initial formulation results in a very low influence of the damage on the intergranular shear stress. A new formulation leading to a direct influence of the damage on the intergranular shear stress has been proposed. This new formulation has resulted in (i) a change in the distribution of micro-cracks, with coalescence between the different micro-cracks, and (ii) a large increase in the RCF life estimated by the model. The order of magnitude of the number of cycles corresponding to the first micro-cracks is comparable to that given by experiments.