Cette étude, inscrite dans le Programme National de Recherche NOREV (Nouveaux REfractaires Verriers), porte sur la simulation numérique par éléments finis du refroidissement post-coulée de réfractaires électrofondus à très haute teneur en zircone (THTZ). Un modèle numérique thermique du refroidissement d'une dalle a été développé, à partir de mesures thermiques réalisées dans un moule lors d'une coulée en usine. Ce calcul a permis de situer dans le temps, pour les différents points de la dalle, les phénomènes intervenant au cours du refroidissement : solidification, fluage et transformation de la zircone. Des essais mécaniques isothermes et anisothermes ont permis de définir les équations constitutives du comportement en fluage (écrouissage isotrope) et lors de la transformation de phase de la zircone (modèle de plasticité de transformation de Leblond). A partir des chargements thermiques, les contraintes générées au cours du refroidissement ont été calculées, ainsi que les contraintes résiduelles après retour à température ambiante. La transformation de phase de la zircone, qui induit un gonflement volumique de 4%, semble s'avérer être le phénomène physique le plus critique pour ces deux types de contraintes. Les études paramétriques réalisées ont prouvé l'intérêt du modèle numérique en tant qu'outil d'aide à l'amélioration du procédé en jouant sur les propriétés de divers paramètres : propriétés du matériau coulé, géométrie du moule, propriétés de ses constituants,...