v7.22.120 HSNV120 - Traction hyperélastique d’un barreau sous chargement thermique#

Résumé:

Ce test thermomécanique quasi-statique consiste à chauffer uniformément un barreau parallélépipédique, le soumettre à une traction importante pour finalement le laisser revenir dans un état déchargé. On valide ainsi la cinématique des grandes déformations hyperélastiques (commande STAT_NON_LINE, mot-clé COMPORTEMENT) pour une relation dec comportement élastique non-linéaire (ELAS_VMIS_LINE et ELAS_VMIS_TRAC) avec chargement thermique.

Le barreau est modélisé par un élément volumique (HEXA20, modélisationA).

Les résultats obtenus par Aster ne diffèrent pas de la solution théorique.

Solution de référence#

Méthode de calcul utilisée pour la solution de référence#

On cherche le champ de déplacement \(U\) sous la forme:

../../../../_images/Object_6124.svg

Le gradient de la transformation, la déformation et sa part mécanique sont alors:

../../../../_images/Object_7116.svg

avec:

../../../../_images/Object_8119.svg

Remarque:

\({({E}^{m})}_{\text{eq}}=\mid a-b\mid =a-b\) (on suppose que \(a>b\) )

La relation de comportement s’écrit:

../../../../_images/Object_10108.svg

avec:

../../../../_images/Object_11119.svg

Pour déterminer

../../../../_images/Object_12100.svg

en tenant compte de l’écrouissage linéaire, on introduit:

  • le module de cisaillement:

    ../../../../_images/Object_1385.svg
  • le module d’écrouissage:

    ../../../../_images/Object_1468.svg

,

La « pseudo variable interne »

../../../../_images/Object_1558.svg

vaut alors:

../../../../_images/Object_1667.svg

Finalement,

../../../../_images/Object_1749.svg

s’écrit:

../../../../_images/Object_1847.svg

En tenant compte des conditions aux limites:

../../../../_images/Object_1949.svg

Le système à résoudre s’écrit:

../../../../_images/Object_2040.svg

Il s’écrit aussi:

../../../../_images/Object_2179.svg

A

../../../../_images/Object_2238.svg

fixé, il s’agit donc d’un système non linéaire en

../../../../_images/Object_2344.svg

et

../../../../_images/Object_2440.svg

, puisque

../../../../_images/Object_2537.svg

est quadratique en

../../../../_images/Object_2636.svg

et

../../../../_images/Object_2733.svg

quadratique en

../../../../_images/Object_2834.svg

.

Néanmoins, on peut choisir de fixer

../../../../_images/Object_2930.svg

(donc

../../../../_images/Object_3030.svg

) et résoudre un système linéaire en

../../../../_images/Object_3169.svg

et

../../../../_images/Object_3225.svg

(duquel on déduit

../../../../_images/Object_3325.svg

et

../../../../_images/Object_3425.svg

):

Il reste alors à exprimer la contrainte de Cauchy:

../../../../_images/Object_3918.svg

Soit ici:

../../../../_images/Object_4013.svg

Quant à la force exercée sur la face [3,4], du fait de l’hypothèse de charges mortes, elle s’écrit simplement:

../../../../_images/Object_4127.svg

Résultats de référence#

On adoptera comme résultats de référence les déplacements, la contrainte de Cauchy et la force exercée sur la face \([3,4]\) (en \(\mathrm{3D}\) seulement):

Au temps \(t=\mathrm{2 }s\) ( \(\Delta T=100°C\) , traction

../../../../_images/Object_4319.svg

)

En fait, on cherche

../../../../_images/Object_4418.svg

tel que l’allongement:

../../../../_images/Object_4517.svg
  • \(K=166\mathrm{666 }\mathit{MPa}\) \(\mu =\mathrm{76 923 }\mathit{MPa}\) \(R'=\mathrm{2 020 }\mathit{MPa}\)

  • \(a=0.095\)

  • ../../../../_images/Object_5010.svg
  • ../../../../_images/Object_5129.svg
  • ../../../../_images/Object_5216.svg
  • ../../../../_images/Object_5316.svg

Au temps \(t=\mathrm{3 }s\)

../../../../_images/Object_5417.svg

Le barreau est revenu dans son état initial:

../../../../_images/Object_5517.svg

Incertitude sur la solution#

La solution est analytique. Aux erreurs d’arrondis près, on peut la considérer exacte.

Références bibliographiques#

On pourra se référer à :

    1. LORENTZ : Une relation de comportement hyperélastique non linéaire - Note interne EDF DER HI-74/95/011/0

Modélisation A#

Caractéristiques de la modélisation#

Modélisation volumique:

1 maille HEXA20 1 maille QUAD8

../../../../_images/10001B9600001C02000015309D009CE50DA8878C.svg

Conditions aux limites:

\(\mathit{N2}\) :

\(\mathit{N1}\) : \(\mathit{N6}\) :

../../../../_images/Object_5713.svg ../../../../_images/Object_5815.svg

\(\mathit{N9}\) , \(\mathit{N13}\) , \(\mathit{N14}\) , \(\mathit{N5}\) , \(\mathit{N17}\) :

Charge: Traction sur la face \([348711161915]\)

Caractéristiques du maillage#

Nombre de nœuds: 20

Nombre de mailles: 2

1 HEXA20

1 QUAD8

Grandeurs testées et résultats#

Identification

Référence

\(t=2\) Déplacement \(\mathit{DX}\) (\(\mathit{N8}\) )

100

\(t=2\) Déplacement \(\mathit{DY}\) (\(\mathit{N8}\) )

–37

\(t=2\) Déplacement \(\mathit{DZ}\) (\(\mathit{N8}\) )

–37

\(t=2\) Contraintes \(\mathit{SIGXX}\) (\(\mathit{PG1}\) )

1399.66

\(t=2\) Contraintes \(\mathit{SIGYY}\) (\(\mathit{PG1}\) )

11013.986

(4915)#\[ \begin{align}\begin{aligned}t=2\\\textrm{Contraintes :math:`\mathit{SIGZZ}` (:math:`\mathit{PG1}` )}\end{aligned}\end{align} \]

\(t=2\) Contraintes \(\mathit{SIGXY}\) (\(\mathit{PG1}\) )

0

(4915)#\[ \begin{align}\begin{aligned}t=2\\\textrm{Contraintes :math:`\mathit{SIGXZ}` (:math:`\mathit{PG1}` )}\end{aligned}\end{align} \]

\(t=2\) Contraintes \(\mathit{SIGYZ}\) (\(\mathit{PG1}\) )

0

\(t=2\) Variable

../../../../_images/Object_6012.svg
width:

361

height:

273

\(\mathit{VARI}\) (\(\mathit{PG1}\) )

8.9110–2

(4915)#\[ \begin{align}\begin{aligned}t=3\\\textrm{Déplacement :math:`\mathit{DX}` (:math:`\mathit{N8}` )}\end{aligned}\end{align} \]

\(t=3\) Déplacement \(\mathit{DY}\) (\(\mathit{N8}\) )

0

(4915)#\[ \begin{align}\begin{aligned}t=3\\\textrm{Déplacement :math:`\mathit{DZ}` (:math:`\mathit{N8}` )}\end{aligned}\end{align} \]

\(t=3\) Contraintes \(\mathit{SIGXX}\) (\(\mathit{PG1}\) )

0

(4915)#\[ \begin{align}\begin{aligned}t=3\\\textrm{Contraintes :math:`\mathit{SIGYY}` (:math:`\mathit{PG1}` )}\end{aligned}\end{align} \]

\(t=3\) Contraintes \(\mathit{SIGZZ}\) (\(\mathit{PG1}\) )

0

(4915)#\[ \begin{align}\begin{aligned}t=3\\\textrm{Contraintes :math:`\mathit{SIGXY}` (:math:`\mathit{PG1}` )}\end{aligned}\end{align} \]

\(t=3\) Contraintes \(\mathit{SIGXZ}\) (\(\mathit{PG1}\) )

0

(4915)#\[ \begin{align}\begin{aligned}t=3\\\textrm{Contraintes :math:`\mathit{SIGYZ}` (:math:`\mathit{PG1}` )}\end{aligned}\end{align} \]

\(t=3\) Variable

../../../../_images/Object_6125.svg
width:

361

height:

273

\(\mathit{VARI}\) (\(\mathit{PG1}\) )

0

\(t=2\) Force nodale \(\mathit{DX}\) (\(\mathit{N8}\) )

–1.0817108

(4915)#\[ \begin{align}\begin{aligned}t=2\\\textrm{Force nodale :math:`\mathit{DY}` (:math:`\mathit{N8}` )}\end{aligned}\end{align} \]

\(t=2\) Force nodale \(\mathit{DZ}\) (\(\mathit{N8}\) )

0

Remarques#

Calcul de la force nodale:

La force appliquée sur la face \([3,4]\) ,

../../../../_images/Object_6215.svg

, se répartit entre les différents nœuds suivant la pondération suivante:

  • nœuds sommets:

    ../../../../_images/Object_6315.svg
  • nœuds milieux:

    ../../../../_images/Object_6415.svg
../../../../_images/1000045A00000F6600000CD1D051E52336168495.svg

Synthèse des résultats#

Les résultats numériques et analytiques coïncident remarquablement.