N Billon

Approche visco hyperélastique basée sur une théorie de réseau; Application au PMMA 1

Keywords: Modèle de comportement, Anélasticité, Visco élasticité, Grandes déformations, Polymères, PMMA

Une extension d'une théorie de réseau modélise le comportement visco hyperélastique des polymères au dessus de Tg. Les effets anélastiques découlent de l'évolution de paramètres représentant le dé enchevêtrement ou le durcissement du matériau. Il s'agit d'une approche mésoscopique conservant un sens physique sans pour autant être une approche micromécanique. Ecrit dans un cadre thermodynamique rigoureux cette approche est couplée à une utilisation des principes d'équivalences temps température qui autorise la modélisation des effets combinés de la température et de la vitesse de chargement tout en réduisant le nombre de paramètres. Une application complète est décrite dans le cas d'un PMMA.

Introduction

La non linéarité du comportement des polymères, la combinaison d'effets viscoélastiques et/ou viscoplastiques, ou encore, la coexistence de processus d'endommagement et de durcissement structural ainsi que, la sensibilité à la température, régie par des relaxations discrètes, rendent la modélisation mécanique dans les polymères non triviale et interdisent, à ce jour, de proposer un modèle de comportement générique unique.

De nombreux modèles ont été proposés dans la littérature. Le plus souvent ils décomposent, a priori, les déformations globales en des composantes élémentaires réversibles ou irréversibles, dépendantes du temps ou non. La combinaison se fait sous hypothèse d'additivité soit des déformations soit des contraintes. Le formalisme mathématique associé est, quant à lui, construit soit à partir de modèles analogiques 1D généralisés soit dans un cadre thermodynamique où les contraintes sont régies par des potentiels [1,2]. En terme d'efficacité, ces deux routes sont équivalentes lorsqu'elles sont menées rigoureusement.

Dans les deux cas les formalismes impliquent des grandeurs, associés à chacune des composantes (modules, viscosité, seuils, écoulement plastiques etc.), qui peuvent dépendre beaucoup de la température. Il en résulte souvent un très grand nombre de paramètres qui n'ont pas toujours un plein sens physique et qui nécessitent parfois la mise en place de stratégie d'identification lourde. La présente étude [3] vise à mettre en place une voie de modélisation qui introduirait les effets irréversibles et les effets du temps autrement que par cette partition a priori. En parallèle, dans l'objectif de réduire le nombre de paramètres du modèle, nous proposons d'utiliser le principe d'équivalence « tempstempérature » pour exprimer, de manière combinée, les dépendances en vitesse d déformation et les dépendances en température.

Le premier objectif est atteint de manière rigoureuse dans un cadre thermodynamique. Le deuxième s'appuie sur une validation expérimentale et une modélisation phénoménologique.

Le présent article résume les hypothèses de base du modèle ainsi qu'une validation de sa limite 1D sur la base d'essais de traction réalisés sur un PMMA.

Fondement théorique

Comportement dépendant du temps

On admettra que le coeur du comportement de tout polymère est assimilable à celui d'un réseau virtuel de chaînes non Gaussiennes enchevêtrées. Les contraintes sont donc régies par un potentiel, W, qui dépend de la déformation élastique (i .e. des extensions principales, élastiques), de la température, T, et de variables d'états (ou de paramètres) dépendant de la microstructure, m i=1..n . Si nous limitons, dans une première approche, notre analyse aux matériaux incompressibles et initialement isotropes nous pouvons donc écrire que la contrainte σ i (i=1, 2, 3), obéit à l'équation (1).

( ) ( )      = - λ ∂ … ∂ λ = σ      = λ λ λ - λ ∂ … λ λ λ ∂ λ = σ 1 I P m , ,
       ∂ λ ∂ λ - = ∂ λ ∂ ⇔ = ∂ λ ∂ λ + ∂ λ ∂ λ = ∂ λ ∂ ≥ β ∂ λ ∂ λ ∂ ∂ β - = ∂ ∂ ∂ ∂ ∑ ∑ = = t ln t 0 t t t ) 1 ' avec ( t W ' 1 t m m W a i e i re i re i a i a i e i i 3 1 i re i e i n 1 j j j (2)
L'énergie additionnelle, source du paramètre β', peut être, quant à elle, dissipée en chaleur ou stockée sous forme de microstructure.

Nous supposons ici implicitement que les extensions principales élastiques, anélastiques et globales sont colinéaires. Cette hypothèse de confort pourra être levée si nécessaire mais demandera, bien sûr, une modélisation complémentaire. En l'état, il suffit, maintenant, de donner une forme mathématique aux termes t m j ∂ ∂ pour compléter le modèle. Pour cela il faut donner une forme à W et un sens physique à ses paramètres.

Choix du potentiel de déformation

Nous choisissons (Eq. ( 3)) le modèle d'Edwards et Vilgis [4] qui représente les polymères comme un réseau de chaînes non Gaussiennes comportant une densité de points de réticulation, N c , auxquels s'ajoutent une densité d'enchevêtrements, N s . La liberté conformationnelle apportée par les enchevêtrements est représentée par le concept d'anneaux glissants, introduit par Ball [5]. Dans l'expression finale du potentiel apparaît, de ce fait, un paramètre η (i.e. longueur de glissement) qui est tel qu'une valeur nulle correspond à un point de réticulation permanent. Une valeur positive correspond, quant à elle, à un enchevêtrement d'autant plus 'libre' que η est élevé. Les chaînes sont supposées avoir une extensibilité limitée, ce qui se réduit à un paramètre α positif inférieur à 1/3. Il en découle que (Eq. ( 3)) : Ce potentiel a déjà démontré sont aptitude à reproduire partiellement le comportement des polymères quand il est associé à des éléments visqueux ou plastiques dans plusieurs modèles construits « classiquement » (cf. plus haut) [7][8][9][10][11].

( ) ( ) ( ) ( ) ( )( ) 2
A ; 3 I 2 I Z ; I 1 X ; I I 1 Y X ln I X 1 2 N X ln Y ln Z Y X A 2 N W α - η + = η + η + = α - = η + η + η + =         + α - +         + + =

Cinétique d'évolution des paramètres

Dans notre cas, les effets de temps sont introduits par une évolution des paramètres du modèle, c'est-à-dire de

α η ou et N N s c / , , ~.
Pour illustrer l'idée, on peut imaginer qu'un processus de dé-enchevêtrement conduirait à une augmentation de η. Une cristallisation induite par la déformation pourrait conduire à une augmentation de α ou/et une augmentation de c N au détriment de s N .

Si on se limite, pour la lisibilité, à une formulation 1D (Eq. ( 2)), la vitesse de déformation anélastique devient au terme de ce modèle (Eq. ( 4)) :

( )

1 ' ' 2 2 ~≥         ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ = ∂ ∂ β α α η η σ β ε avec t W t N N W t N N W t W t s s c c a (3)
Le dernier point à traiter est de postuler une forme mathématique pour les cinétiques d'évolution :

t et t N t N t s c ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ 2 , , α η .
En première approximation des formes phénoménologiques, bâties à partir des observations expérimentales, peuvent être proposées. Nous pouvons faire l'hypothèse, a priori, que ces cinétiques dépendront de l'énergie élastique stockée dans le matériau (Eq. ( 4)) :

( ) ( ) ( ) ( ) ( ) ( ) 2 2 2 s 3 1 ln 1 ln 3 1 1 3 X ln Y ln Z Y X A f α - +       η + + α - α - -         + + = αλ ηλ ηλ ηλ αλ αη (4)
Plus précisément, on peut supposer que plus cette énergie sera grande dans les enchevêtrements plus le processus de dé-enchevêtrement sera rapide. Les potentialités de cette démarche sont décrites par ailleurs [3]. Nous ne reproduirons ici qu'un exemple simple où l'on admet qu'un seul processus de dissipation (seul 0 t ≠ ∂ η ∂

) existe, contrôlé par la forme la plus élémentaire de cinétique imaginable (Eq. ( 5)) :

s f t ξ = ∂ η ∂ (5)
Ainsi, cette démarche autorise de reproduire des effets viscoélastiques (boucles d'hystérésis et effets partiels de vitesse) sans décomposition, a priori, des déformations et au prix de seulement de 4 paramètres au plus. Ce nombre se réduit à 3 dans un réseau de chaînes Gaussiennes (α=0) (Fig. 1). 5)) et une vitesse de 0,01 s -1 (2) ou une vitesse de 0,001 s -1 (3). .

4

1 ' ; 0 ]; 11 [ 2343 , 0 initial ; N 1 , 0 N s c = β = α = η =

Effets de température et de vitesse

La prise en compte d'une cinétique d'évolution microstructurale suffit donc à introduire une certaine sensibilité à la vitesse de déformation. Cette méthode ne s'avère pourtant pas toujours suffisante pour reproduire toutes les sensibilités observées. Il faut donc admettre que tout ou partie de paramètres dépend de la température et de la vitesse de sollicitation. Partant d'observations expérimentales qui tendent à démontrer que le comportement des polymères peut parfois être paramétré, non pas en terme de couple -vitesse de sollicitation / température-mais en terme de vitesse de déformation à une température de référence, une écriture simplifiée peut être proposée. Cette équivalence est illustrée ici dans le cas d'un PMMA mais des conclusions équivalentes ont pu être dressées pour de PAs, des PETs, des PUs et un PS.

Dans ce cas (où il est possible de mettre en évidence une équivalence) nous suggérons d'exprimer les paramètres du modèle en fonction de la vitesse équivalente à une température de référence pour réunir la double sensibilité à la température et à la vitesse en une seule (Eq. ( 6)). Les facteurs de décalage, De ces observations nous pouvons écrire, de manière phénoménologique, que les paramètres de notre modèles obéissent à une loi donnée par l'Eq. ( 7) :

2 0 1 0 m 2 . T / T 2 ' 2 2 ~2 c ~2 s m 2 . T / T 1 ' 1 1 ~1 c ~1 s 0 ~0 c ~0 s c s a 1 1 N N 2 a 1 1 N N 2 1 N N ' N N - -       ε τ +               β ξ +       ε τ +               β ξ +               ξ =               β ξ (7)
Il est ainsi possible de reproduire des essais de traction non monotones réalisés entre 120 et 150°C et de 10 -4 à 10 -2 s -1 . Il est aussi possible de modéliser, avec les mêmes coefficients, des essais de traction interrompus par des relaxations (Figure 3).

Conclusion

Dans cette étude un nouveau modèle de comportement est proposé qui autorise de rendre compte de la visco élasticité et de la viscoplasticité des polymères sur une large gamme de températures et de vitesses jusqu'à de grandes déformations.

Ce modèle est construit de manière rigoureuse dans un formalisme thermodynamique 3D. Il est écrit en grandes transformations. Son originalité réside en ce que, à la différence des approches classiques, il est bâti à partir de concept de dynamique de chaîne, réalisant un mariage entre la mécanique et la physico chimie des polymères.

Il autorise la prise en compte des couplages thermomécaniques forts et a prouvé son efficacité lors de chargements non monotones, cyclés ou interrompus par des relaxations, au prix d'un faible nombre de paramètres tous identifiables. Il a été validé dans le cas d'un PMMA, d'un PS, d'un PET et celui d'un PA 6-6.

Des travaux futurs vont s'attacher à introduire explicitement une description des microstructures cristallines et à étende l'écriture aux cas de sollicitations de type fatigue.

  la constante de Boltzmann).

3 FIG. 1 -

 31 FIG. 1 -Courbes de traction (charge -décharge -recharge) calculées à vitesse de déformation constante. Comparaison entre : Le modèle d'Edwards Vilgis initial (1) et le modèle proposé pour un paramètre 01 , 0 = ξ (Eq.(5)) et une vitesse de 0,01 s -1 (2) ou une vitesse de 0,001 s -1 (3).

6 )

 6 estimés de manière classique sur la base d'essais DMA. Puis, les essais mécaniques, réalisés à la température T et la vitesse . ε , sont dépouillés en fonction de la vitesse équivalente : La démarche complète est ici illustrée à l'aide d'un PMMA ALTUGLAS V825 T commercialisé par la société ARKEMA et de masses molaires M n de 40 434 g/mol et M w de 93 000 g/mol.Des essais DMA autorisent la construction d'une courbe maîtresse dans la zone de visco élasticité linéaire. A 130 °C les facteurs de décalage obéissent à une loi de WLF (Eq. (6)) :La figure 2 illustre que les résultats expérimentaux à plus grandes déformations, i .e. hors de la zone linéaire, obéissent encore à cette équivalence, c'est-à-dire que la vitesse équivalente à 130 °C peut suffire à paramétrer les conditions expérimentales multiples.

FIG. 2 -FIG. 3 -

 23 FIG.2-Courbes de traction (charge -décharge) expérimentales et calculées à vitesse de déformation constantes pour un PMMA. Comparaison entre différentes conditions expérimentales qui conduisent à des vitesses équivalentes à la température de 130°C proches.

Ce travail a été réalisé dans le cadre du projet STREP Apt-Pack (STREP 505204-1 financé par la commission Européenne.