Hyperélasticité

L'hyperélasticité est un formalisme mathématique utilisé en résistance des matériaux pour décrire la relation contrainte-déformation de certains matériaux grandement déformables (polymères thermoplastiques, polymères thermodurcissables, élastomères, tissus biologiques). Contrairement à l'élasticité linéaire définie explicitement par la loi de Hooke pour les petites déformations, en hyperélasticité, on postule l'existence d’une densité d’énergie de déformation notée W dont les dérivées par rapport à la déformation dans une direction donnée donnent l'état de contrainte au sein du matériau dans cette même direction. Physiquement, W représente la quantité d’énergie élastique que le matériau emmagasine en fonction de l’étirement imposé.

Tous les polymères ne présentent pas le même comportement mécanique en fonction de leur microstructure (longueur des macromolécule, nature des composants chimiques, degré de réticulation, possibilité de cristallisation sous-tension, température d'utilisation, etc.). Une loi de comportement hyperélastique est donc une expression de W permettant de décrire le comportement particulier d'un matériau. Il en existe une multitude pouvant être séparées en deux catégories : les modèles phénoménologiques et les modèles statistiques.

Formalisme mathématique

En mécanique des milieux continus, il existe plusieurs tenseurs de déformation et de contrainte permettant de décrire le comportement du matériau dans une configuration donnée :

• le tenseur des contraintes de Cauchy représente la contrainte dite « vraie » qui est réellement ressentie par le matériau dans la configuration actuelle, c'est-à-dire la force divisée par la section à l'instant t. $\sigma ={F \over S(t)}$ Cette contrainte est conjuguée à la déformation dans la configuration actuelle, c'est-à-dire au tenseur d'Euler-Almansi ${\mathbf {e}}={1 \over 2}({\mathbf {Id}}-{\mathbf {B^{-1}}})$ où ${\mathbf {Id}}$ est la matrice identité et ${\mathbf {B}}={\mathbf {F.F^{T}}}$ le tenseur des dilatations de Cauchy-Green droit ;
• le tenseur de Piola-Kirchhoff 2 (abrégé « PK2 ») représente la contrainte « ingénieure », c'est-à-dire la contrainte que ressentirait le matériau s'il était dans sa configuration initiale, soit la force divisée par la section initiale $\pi ={F \over S_{0}}$. Cette contrainte est conjuguée à la déformation dans la configuration initiale, c'est-à-dire au tenseur de Green-Lagrange ${\mathbf {E}}={1 \over 2}({\mathbf {C}}-{\mathbf {Id}})$ où ${\mathbf {Id}}$ est la matrice identité et ${\mathbf {C}}={\mathbf {F^{T}.F}}$, le tenseur des dilatations de Cauchy-Green gauche.

En hyperélasticité, on postule l'existence d’une densité d’énergie de déformation notée W dont les dérivés donnent l'état de contrainte au sein du matériau. Ainsi, la manière de dériver W donnera des contraintes différentes. La dérivée de W par rapport à ${\mathbf {B}}$ donnera la contrainte de Cauchy alors que sa dérivée par rapport à ${\mathbf {C}}$ donnera la contrainte de PK2. Mathématiquement, on écrira :

${\mathbf {\sigma }}=2{\partial W \over \partial {\mathbf {B}}}$

${\mathbf {\pi }}=2{\partial W \over \partial {\mathbf {C}}}$

La relation entre ces deux tenseurs est donnée par : $J{\mathbf {\sigma }}={\mathbf {F.\pi .F^{T}}}$

où $J=det({\mathbf {F}})$ est le Jacobien de la transformation, c'est-à-dire la variation de volume.

Cas des matériaux isotropes

Si le matériau est isotrope, il est possible de développer les lois de comportement grâce aux invariants des tenseurs de dilatation :

${\mathbf {\sigma }}=2{\Biggl (}{\Biggl (}{\partial W \over \partial I_{1}}+I_{1}{\partial W \over \partial I_{2}}{\Biggr )}{\mathbf {B}}-2{\partial W \over \partial I_{2}}{\mathbf {B}}^{2}-{\partial W \over \partial I_{3}}I_{3}{\mathbf {Id}}{\Biggr )}$

${\mathbf {\pi }}=2{\Biggl (}{\Biggl (}{\partial W \over \partial I_{1}}+I_{1}{\partial W \over \partial I_{2}}{\Biggr )}{\mathbf {Id}}-{\partial W \over \partial I_{2}}{\mathbf {C}}-{\partial W \over \partial I_{3}}I_{3}{\mathbf {C^{-1}}}{\Biggr )}$

Avec $I_{1}=tr({\mathbf {C}})$, le premier invariant

$I_{2}={1 \over 2}(tr({\mathbf {C}})^{2}-tr({\mathbf {C}}^{2}))$, le second invariant

$I_{3}=det({\mathbf {C}})$, le troisième invariant.

Cas des matériaux incompressibles

Si, de plus, le matériau est incompressible ($I_{3}=1$), il est possible de réécrire ces équations en introduisant le multiplicateur de Lagrange noté p, assimilable à la pression hydrostatique, tel que :

${\mathbf {\sigma }}=2{\Biggl (}{\Biggl (}{\partial W \over \partial I_{1}}+I_{1}{\partial W \over \partial I_{2}}{\Biggr )}{\mathbf {B}}-2{\partial W \over \partial I_{2}}{\mathbf {B}}^{2}-p\,{\mathbf {Id}}{\Biggr )}$

${\mathbf {\pi }}=2{\Biggl (}{\Biggl (}{\partial W \over \partial I_{1}}+I_{1}{\partial W \over \partial I_{2}}{\Biggr )}{\mathbf {Id}}-{\partial W \over \partial I_{2}}{\mathbf {C}}-p\,{\mathbf {C^{-1}}}{\Biggr )}$

Lois de comportements existantes

Dans cette section, les expressions mathématiques des modèles les plus couramment utilisés seront détaillées.

Modèles phénoménologiques

Dans l’approche dite « phénoménologique », on cherche à décrire le comportement macroscopique de matériau sans nécessairement en comprendre les fondements physiques. Mooney (1940) fut le premier à proposer un potentiel d’énergie de déformation à partir des invariants du tenseur des dilatations de Cauchy-Green. Rivlin (1948) en a donné une expression simplifiée donnant le modèle « Néo-hookéen », puis une forme généralisée grâce à la décomposition en série de Taylor. Par la suite, plusieurs autres auteurs ont apporté des modèles phénoménologiques de plus en plus sophistiqués : Ogden, Gent, Yeoh, Marlow pour ne citer qu'eux. Les expressions de ces modèles sont données à partir d'invariants ($I_{1},I_{2},I_{3}$) ou d'élongations principales ($\lambda _{1},\lambda _{2},\lambda _{3}$). Rappelons qu’il est nécessaire de soustraire 3 à chaque invariant (resp. à chaque somme des élongations principales) afin d’assurer que la densité d’énergie de déformation soit nulle lorsque le matériau n’est pas sollicité, c'est-à-dire lorsque les invariants (resp. la sommes des élongations principales) sont égaux à 3.

Modèle Néo-Hookéen (1948)

$W(I_{1})=C_{0}(I_{1}-3)$ où $C_{0}$ est le module d'élasticité du matériau dont la dimension est celle d'une contrainte

Modèle de Mooney-Rivlin (1940-1948)

Le modèle de Rivlin généralisé est donné par :

$W(I_{1},I_{2})=\sum _{i,j}^{N}C_{ij}(I_{1}-3)^{i}(I_{2}-3)^{j}$

Un cas particulier est donné par le modèle de Mooney-Rivlin dans le cas où $N=1$ et $i\neq j$

$W(I_{1},I_{2})=C_{10}(I_{1}-3)+C_{01}(I_{2}-3)$

où les $C_{ij}$ sont les modules d'élasticité du matériau

Modèle de Gent

Contrairement au modèle Néo-hookéen et de Mooney-Rivlin qui montre une certaine linéarité en grande déformation, le modèle de Gent (en) reproduit correctement le durcissement important des élastomères avant rupture dû à la limite d'extensibilité des chaînes macromoléculaires. Gent a imaginé une densité d'énergie de déformation présentant une singularité quand le premier invariant atteint une valeur limite notée $I_{m}$

$W(I_{1})=-\mu (I_{m}-3).ln{\Bigl (}1-{I_{1}-3 \over I_{m}-3}{\Bigr )}$

Les paramètres de ce modèle sont donc $\mu$, module d'élasticité, et $J_{m}=I_{m}-3$.

Représentation graphique de différents modèles hyperélastiques.

Modèle d'Ogden

Le modèle d'Ogden pose : $W(\lambda _{1},\lambda _{2},\lambda _{3})=\sum _{i=1}^{N}{\mu _{i} \over \alpha _{i}}(\lambda _{1}^{\alpha _{i}}+\lambda _{2}^{\alpha _{i}}+\lambda _{3}^{\alpha _{i}}-3)$

où $\alpha _{i}$ et $\mu _{i}$ sont les paramètres du modèle. Contrairement aux autres modèles, ces paramètres n'ont pas de signification physique.

Modèle de Yeoh (1993)

Le modèle de Yeoh pose : $W(I_{1},I_{2})=C_{10}(I_{1}-3)+C_{20}(I_{2}-3)^{2}+C_{30}(I_{1}-3)^{3}$

où les $C_{ij}$ sont les modules d'élasticité du matériau.

Modèles statistiques

Dans l’approche statistique, il est question de reproduire le comportement mécanique macroscopique par des considérations microstructurales physiquement fondées. On distingue les modèles gaussiens et non-gaussiens. Flory est en 1953 le premier à proposer un modèle statistique pour décrire le comportement des polymères en utilisant une répartition gaussienne des chaînes, c’est-à-dire que la longueur des chaînes est supposée suivre une loi normale. Il considère que l’entropie du réseau est la somme des entropies élémentaires des chaînes constitutives, que la longueur de chaque chaîne varie dans le même rapport que les dimensions du milieu macroscopique et que la déformation est isochore (sans variation de volume). En prenant en considération toutes ces hypothèses physiques, une densité d’énergie de déformation est définie en fonction des trois élongations principales :

$W={nkT \over 2}(\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}-3)$

où n est le nombre de chaînes par unité de volume, k la constante de Boltzmann et T la température absolue.

On remarque premièrement qu’en posant $C_{0}={nkT \over 2}$ et en remplaçant la somme des carrés des élongations principales par le premier invariant du tenseur des dilations, on retrouve le modèle néo-hookéen. Le modèle de Flory ne s’avère précis que pour les petites déformations au même titre que le modèle néo-hookéen. De plus, la répartition gaussienne suppose que les chaînes moléculaires puissent atteindre une longueur infinie, ce qui est physiquement impossible. C’est la raison pour laquelle la théorie statistique à répartitions non-gaussienne a vu le jour. Ces derniers améliorent les prédictions pour une large gamme de déformation. Le premier modèle non-gaussien a été introduit par Kuhn et Grun (1942) décrivant l’élasticité d’une chaîne isolée composée de N monomères, tous de longueur L . La longueur initiale d’une chaîne $r_{0}=L{\sqrt {N}}$ est déterminée statistiquement, ainsi que sa longueur maximale $r_{0}=LN$. L’expression de la densité d’énergie d’une chaîne isolée est formulée comme suit :

$W={nkT \over 2}{\biggl (}{\lambda \over {\sqrt {N}}}\beta +ln{\biggl (}{\beta \over sinh(\beta )}{\biggr )}{\biggr )}$

où $\beta =L^{-1}{\biggl (}{\lambda \over {\sqrt {N}}}{\biggr )}$ , $\lambda ={r \over r_{0}}$ et $L^{-1}$ désigne la bijection réciproque de la fonction de Langevin. Cette fonction de Langevin inverse permet de reproduire le durcissement lors des grandes déformations.

De manière générale, tous les modèles statistiques sont basés sur cette théorique statistique non-gaussienne.

Modèle d'Arruda-Boyce (2016)

Le plus célèbre des modèles statistiques est le modèle d'Arruda-Boyce. Ellen M. Arruda et Mary C. Boyce ont développé le modèle dit « à huit chaînes » en considérant un élément cubique dont les demi-diagonales définissent les huit directions privilégiées de déformation. Le solide complet étant composé d'une multitude d'éléments cubiques. Du fait de la symétrie du modèle par rapport aux trois directions principales du tenseur de déformation, les huit chaînes se déforment toutes de la même quantité $\lambda _{ch}={\sqrt {I_{1}/3}}$:

$W={nkT \over 3}{\sqrt {N}}{\biggl (}\lambda _{ch}\beta _{ch}+{\sqrt {N}}ln{\biggl (}{\beta _{ch} \over sinh(\beta _{ch})}{\biggr )}{\biggr )}$

avec $\beta _{ch}=L^{-1}{\biggl (}{\lambda _{ch} \over {\sqrt {N}}}{\biggr )}$.

Modèle de Gent

En 2003, Cornelius O. Horgan, chercheur au département de Génie Civil de l'Université de Virginie (USA), et Giuseppe Saccomandi, chercheur au département d'Ingénierie de l'Innovation à l'Université de Lecce en Italie, ont prouvé que le modèle de Gent, originellement phénoménologique, est analogue au modèle d'Arruda-Boyce à huit chaînes. Ainsi, le modèle de Gent est relié à la microstructure du matériau selon les relations suivantes :

$\mu =nkT$

$J_{m}=3(N-1)$.

avec $\mu$ le module d'élasticité du modèle de Gent, $J_{m}$ un scalaire ayant le rôle de limite d'extensibilité, $n$ la densité de chaînes et $N$ le nombre moyen de monomères par chaîne.