Un modèle numérique mis au point pour étudier l'hystérésis des propriétés élastiques des céramiques poreuses en fonction de la température mesure la distribution de l’orientation cristallographique expérimentalement et les valeurs du coefficient de dilatation. Des grains poly-cristallins générés sont assimilés à du vide pour simuler la porosité du matériau et des éléments cohésifs sont insérés le long des joints de grains pour modéliser la fissuration inter-cristalline. Le modèle est appliqué à une céramique de cordiérite poreuse. Les propriétés de fissure obtenues sont utilisées ensuite pour simuler les courbes contrainte-déformation macroscopiques non linéaires générales de la cordiérite, validant ainsi le modèle.