2. Solution de référence#

2.1. 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_6.svg

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

_images/Object_7.svg

avec:

_images/Object_8.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_10.svg

avec:

_images/Object_11.svg

Pour déterminer

_images/Object_12.svg

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

  • le module de cisaillement:

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

    _images/Object_14.svg

,

La « pseudo variable interne »

_images/Object_15.svg

vaut alors:

_images/Object_16.svg

Finalement,

_images/Object_17.svg

s’écrit:

_images/Object_18.svg

En tenant compte des conditions aux limites:

_images/Object_19.svg

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

_images/Object_20.svg

Il s’écrit aussi:

_images/Object_21.svg

A

_images/Object_22.svg

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

_images/Object_23.svg

et

_images/Object_24.svg

, puisque

_images/Object_25.svg

est quadratique en

_images/Object_26.svg

et

_images/Object_27.svg

quadratique en

_images/Object_28.svg

.

Néanmoins, on peut choisir de fixer

_images/Object_29.svg

(donc

_images/Object_30.svg

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

_images/Object_31.svg

et

_images/Object_32.svg

(duquel on déduit

_images/Object_33.svg

et

_images/Object_34.svg

):

  • _images/Object_35.svg
  • _images/Object_36.svg
  • _images/Object_37.svg
  • _images/Object_38.svg

Il reste alors à exprimer la contrainte de Cauchy:

_images/Object_39.svg

Soit ici:

_images/Object_40.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_41.svg

2.2. 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_43.svg

)

En fait, on cherche

_images/Object_44.svg

tel que l’allongement:

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

  • \(a=0.095\)

  • _images/Object_50.svg
  • _images/Object_51.svg
  • _images/Object_52.svg
  • _images/Object_53.svg

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

_images/Object_54.svg

Le barreau est revenu dans son état initial:

_images/Object_55.svg

2.3. Incertitude sur la solution#

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

2.4. 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