Etude numérique de la fissuration d'un milieu viscoélastique - Pastel

FIG. 6.3 – La longueur de fissure en fonction du temps
Le facteur de rupture f est défini comme suit :
f = l − (l3 − ∆l23)
∆l23
6.2. Algorithme de la propagation de fissure
avec l est la longueur actuelle de la fissure. Le nœud 3 correspondant à la position actuelle de la fissure
se décolle lorsque le facteur de rupture f atteint la valeur 1.0 avec une tolérance ftol donnée :
6.2.2 Calcul préliminaire
1 − ftol ≤ f ≤ 1 + ftol
Avant d’appliquer la technique maître-esclave à l’essai de Rulob caractérisé par une géométrie
complexe et un comportement adéquat identifié, nous avons besoins d’examiner la méthode par
rapport à la sensibilité du maillage, la cinématique d’ouverture de la fissure lors de son implantation
dans le code Abaqus. C’est la raison pour laquelle, cette partie s’emploie essentiellement à étudier la
technique de propagation de fissure, les aspects géométrique et rhéologique du matériau sont simplifiés.
Modélisation du barreau viscoélastique fissuré
Le barreau est de forme cylindrique dont les propriétés géométriques et mécaniques sont données
au tableau 6.2. Pendant le temps de chargement, nous appliquons un déplacement imposé sur
la surface supérieure du cylindre. En raison de sa géométrie symétrique et axisymétrique, on ne modélise
qu’un quart du cylindrique comme le montre la figure 6.4.
Nous adoptons pour ce modèle, un matériau viscoélastique linéaire rhéologiquement simple de
type Zener (c.f §3.3.3). Pour intégrer ces propriétés viscoélastiques de ce matériau dans Abaqus, nous
avons besoins les transformer en séries de Prony comme suit :
τf = η E1 + E∞
E1.E∞
(6.1)
(6.2)
(6.3)
113