Ce travail de thèse a porté sur la simulation numérique du comportement thermique et mécanique d'un combustible nucléaire à particules. Il s'agit d'un composite réfractaire constitué d'une matrice de graphite comportant 45 % en fraction volumique de particules sphériquesd'UO2 revêtues de deux couches de pyrocarbone. L'objectif était de développer une modélisationmulti-échelles de ce composite afin d'estimer son comportement moyen, ainsi que les hétérogé-néités des champs mécaniques au sein des constituants. Nous avons modélisé la microstructuredu combustible et généré des échantillons numériques en 3D. Pour cela, des outils de générationde distributions aléatoires de sphères, de maillage et de caractérisation microstructurale, tellela covariance, ont été développés dans le code de calcul Cast3M. Une centaine d'échantillonsnumériques de différentes tailles ont été réalisés. Le comportement thermo-élastique du combustiblea été caractérisé à partir de ces échantillons, à l'aide de calculs de microstructures paréléments finis. Nous avons étudié l'influence de divers paramètres de la modélisation, dont lesconditions aux limites. Nous proposons une méthode pour s'affranchir des effets des conditionsaux limites sur les résultats, appelée méthode d'érosion. Elle s'appuie sur l'analyse des résultatssur un érodé du volume élémentaire. Nous avons alors déterminé les propriétés effectives ducomposite (modules d'élasticité, conductivité thermique, dilatation thermique), ainsi que lesdistributions des champs mécaniques locaux au sein de la matrice. Enfin, nous avons proposéun modèle de changement d'échelles permettant d'obtenir, non seulement les valeurs moyennesdes variables mécaniques dans chaque phase, mais également leurs variances et covariances pourtout chargement macroscopique imposé. Cette approche statistique de changement d'échellespermet ainsi d'estimer la distribution des grandeurs mécaniques au sein de chaque phase ducomposite. / The aim of this work was to perform numerical simulations of the thermal and mechanical behavior of a particle-based nuclear fuel. This is a refractory composite material made of UO2spherical particles which are coated with two layers of pyrocarbon and embedded in a graphitematrix at a high volume fraction (45 %). The objective was to develop a multi-scale modelingof this composite material which can estimate its mean behavior as well as the heterogeneity ofthe local mechanical variables. The first part of this work was dedicated to the modeling of themicrostructure in 3D. To do this, we developed tools to generate random distributions of spheres,meshes and to characterize the morphology of the microstructure towards the finite elementcode Cast3M. A hundred of numerical samples of the composite were created. The secondpart was devoted to the characterization of the thermo-elastic behavior by the finite elementmodeling of the samples. We studied the influence of different modeling parameters, one of themis the boundary conditions. We proposed a method to vanish the boundary conditions effectsfrom the computed solution by analyzing it on an internal sub-volume of the sample obtained byerosion. Then, we determined the effective properties (elastic moduli, thermal conductivity andthermal expansion) and the stress distribution within the matrix. Finally, in the third part weproposed a multi-scale modeling to determine the mean values and the variance and covarianceof the local mechanical variables for any macroscopic load. This statistical approach have beenused to estimate the intra-phase distribution of these variables in the composite material.