Identification des mécanismes de fissuration dans un alliage d ...

112 Etude de la fissuration sous chargement de fretting et de fatigue
Pour cela, un modèle simple d’éprouvette trouée 2D développé par ailleurs
sous Abaqus © pour calculer les contraintes en déformations planes, est modifié
conformément à la figure 3.39. La méthode consiste à réserver deux espaces carrés
en bord de trou pour y placer les fissures (fig. 3.39a). Ensuite, un programme
Matlab 8 permet de remplir automatiquement chaque espace avec une fissure en
générant un maillage adapté pour le calcul du facteur d’intensité de contrainte
(fig. 3.39b) [95]. Enfin le chargement est appliqué et le résultat du calcul fournit
la valeur du facteur K I par évaluation de l’intégrale de contour (fig. 3.39c).
a) b)
pointe de
fissure
bord
du trou
c)
Fig. 3.39: Méthode de calcul du facteur d’intensité de contrainte dans une
éprouvette trouée 2D en déformations planes; a) réservation de deux espaces
pour les fissures dans le maillage de base, b) génération du maillage avec
la fissure et les contours autour de la pointe (agrandissement de la partie
encadrée en a), c) visualisation de la contrainte équivalente de Von Mises
en pointe de fissure.
Trois calculs sont effectués pour 3 longueurs de fissures différentes : 200, 500 et
1000 µm. La figure 3.40 regroupe le tracé du facteur d’intensité de contrainte en
fonction de la longueur de fissure rapportée à la largeur de l’éprouvette 2a/W
(voir fig. 1.3 pour la définition des notations), pour les différentes méthodes de
calcul. Sur la figure 3.40, on constante que pour des valeurs faibles de a, les
deux descriptions suivent une même tendance. Ceci n’est plus du tout vrai dès
que a > 500 µm. Il apparait aussi une bonne corrélation entre les calculs par
8 développé par C. Meunier dans le cadre de son DEA au laboratoire LTDS