------------------------------------------------------
Auteur
------------------------------------------------------
Julien Troufflard (troufflard@univ-ubs.fr)
Gérard Rio (gerard.rio@univ-ubs.fr)

------------------------------------------------------
Mots-cles
------------------------------------------------------
axisymetrique
contact
ISOELAS
HYPO_ELAS3D
comparaison Abaqus

------------------------------------------------------
But du test
------------------------------------------------------
test du contact dans le cas d un maillage axisymetrique (sur la base des 
travaux de these de Emilie Vieville 2014) :
""
  TYPE_DE_CALCUL
    non_dynamique

  para_contact
    CONTACT_TYPE   2
""

------------------------------------------------------
Description du calcul
------------------------------------------------------
simulation de la compression confinee d un joint dans une chemise
  - 2 maillages (joint:Bague14_QC.her, chemise:chem_fine_QC.her)
  - loi ISOELAS pour les 2 maillages
  - maillage axisymetrique quadratique complet (9 points d integration)
  - bien noter que le maillage chem_fine_QC.her a subi une rotation de -90 degres (en particulier pour la position du point d integration 8 de l element 16)
  - utilisation de reference de point d integration (G_) dans les sorties maple
  - precision de calcul egale a celle par defaut dans Abaqus :
         => PRECISION  5e-3
  - pour une comparaison avec Abaqus : voir rubrique "resultats Abaqus"
  
  NB: Dans le cas d'Abaqus, il s'agit d'élément à 8 noeuds donc d'une interpolation quadratique incomplète.

-------------------------------------------------------------
Grandeurs de comparaison
-------------------------------------------------------------
joint :
  pour les noeuds 27 et 77 :
    - position  : X1 X2
    - reactions : R_X1 R_X2
  pour le point d integration 6 de l element 14 :
    - deformation : def_duale_mises
    - contrainte : contrainte_mises

chemise :
  pour le noeud 75 :
    - position  : X1 X2
    - reactions : R_X1 R_X2
  pour le point d integration 8 de l element 16 :
    - deformation : def_duale_mises
    - contrainte : contrainte_mises

---------------------------------------------------
Informations sur les fichiers facultatifs
---------------------------------------------------
- Compression_confinee.inp : mise en donnees Abaqus v6.10

---------------------------------------------------
Comparaison avec des codes de calcul
---------------------------------------------------
** comparaison avec Abaqus (fichier Compression_confinee.inp)

(ci-dessous quelques resultats suite a un calcul fait le 24/02/2015 sur lo-lg2m-001 avec Abaqus v6.10)

joint :
  noeud 83 (equivalent du noeud HZ++ 27) : X1= 16.99199867  X2= 21.51000023  R_X1= 0  R_X2= -7505.22998047
  noeud 49 (equivalent du noeud HZ++ 77) : X1= 22.99478912  X2= 17.75500679
  point 6 element 14 : contrainte_mises= 66.91905975

chemise :
  noeud 801 (equivalent du noeud HZ++ 75) : X1= 22.98790932  X2= 15.74831772
  point 4 element 171 (equivalent du point 8 element 16 HZ++) : contrainte_mises= 160.24816895
  
GR: de manière à approcher au mieux la simulation dans Abaqus, au final une loi hypo-élastique est utilisée d'où par intégration on retrouve une approximation d'une loi élastique avec une mesure logarithmique cumulée. 
NB: Lorsque le pas de temps est limité à 0.02s la différence avec Abaqus se réduit ce qui montre l'importance de la mesure de déformation, même pour des déformations faibles.