Modelling and numerical study of 3D effect on glassy polymer fracture

Le sujet concerne l’influence des effets tridimensionnelles sur les champs de déformation et de contrainte au voisinage d’une éprouvette chargée en mode I. La loi de comportement caractéristique des polymères amorphes avec un seuil d’écoulement viscoplastique suivi d’un adoucissement et d’un durcissement à mesure que la déformation augmente est prise en compte. La loi est implantée dans une UMAT abaqus. Les champs au voisinage de l’entaille sont analysés et les résultats 3D comparés à ceux correspondant à un calcul 2D sous l’hypothèse de déformation plane. L’influence de l’épaisseur de l’échantillon est étudiée et nous montrons qu’au-delà d’un rapport épaisseur t sur rayon d’entaille rt, t/rt>20, les champs de déformation plastique sont qualitativement similaires entre les calculs 3D et 2D. En revanche, nous montrons que la répartition des contraintes et notamment celle de la contrainte moyenne est sur-estimée avec un calcul 2D en comparaison à une simulation 3D. Nous prenons en compte par la suite la rupture par craquelage, modélisée avec un modèle cohésif. Une étude paramétrique est menée afin de définir une procédure d’identification des paramètres caractéristiques du modèle cohésif. Par ailleurs les simulations montrent qu’au-delà d’un rapport t/rt supérieur à 20, une ténacité minimum peut être estimée : ceci constitue un résultat important pour la détermination expérimentale de la ténacité des polymères ductiles.We investigate 3D effect of crack tip palsticity and the influence of the thickness on 3D glassy polymer fracture. The characteristuc constitutive law with a viscoplastic yield stress followed by softening and progressive hardening is accoutne for and implemented in a UMAt routine, in abaqus. The crack tip fields are investigated and 2D plasne strain versus 3D calculations compared. Qualitatively, the palstic distribution are comparableas soon as the ratio thickness over crack tip radius is larger than 20. However, the 2D calculations over estimate the stress distribution compared to the 3D cases. We have accounted for failure by crazing that is described with a cohesive models. A parametric study sheds light on the methodology to use for the calibration of the cohesive parameters. The simulations show that for a ratio thickness over craci tip radius larger than 20, a minimum tuoghness can be observed. This results has implication on the definition of a thickness larger enough experimentally