Les combustibles à particules sont constitués de quelques milliers de billes d'un millimètre de diamètre composées d'oxyde d'uranium enrobé de couches de confinement qui sont noyées dans une matrice graphite pour former un élément combustible. L'objectif de ce travail est de développer un outil de simulation du comportement thermomécanique de ces combustibles sous irradiation permettant une estimation fine des chargements locaux sur ces particules. Le choix s'est porté vers la méthode des éléments-finis au carré, où interviennent deux échelles distinctes de discrétisation : une structure '' macroscopique'' homogène dont les propriétés en chaque point d'intégration sont calculées sur une seconde structure ''microscopique'' hétérogène, le Volume Élémentaire Représentatif (VER). La première partie du travail a porté sur la définition de ce VER. Un indicateur morphologique basé sur la distribution des distances minimales centre-à-centres a été proposé pour sélectionner des tirages aléatoires de microstructure. La réponse macroscopique élastique des VER, calculée par éléments finis, a été comparée à un modèle analytique. Des indicateurs de représentativité thermique et mécanique du chargement local ont été construits à partir les modes de rupture de la particule. Une étude statistique de ces critères sur une centaine de VER a démontré l'importance de sélectionner une microstructure représentative. Il a dans cette optique été développé un modèle empirique reliant l'indicateur morphologique à l'indicateur mécanique. La seconde partie du travail traite de deux méthodes de changement d'échelle qui sont basées sur l'homogénéisation des milieux périodiques. Considérant un problème de thermique linéaire avec terme source en régime permanent, il a été montré que l'hétérogénité de la source de chaleur implique l'utilisation d'une méthode au second ordre pour relocaliser correctement le champ de température. Le problème mécanique non-linéaire a lui été traité en utilisant l'algorithme itératif de cast3M, en substituant à l'intégration de la loi de comportement un calcul élément finis sur le VER. Cet algorithme a été validé, puis couplé à la résolution thermique afin de simuler un chargement d'irradiation. Un calcul sur un élément combustible complet a mis en évidence une forte interaction entre les deux échelles, ce qui confirme l'intérêt d'un tel modèle pour simuler le comportement de ces combustibles.