
----------------------------------------
Description
----------------------------------------
comparaison entre un calcul Herezh et un calcul UMAT de type Abaqus-Herezh dans un cas irrotationnel
pour une loi additive :
      ISOHYPERBULK3
    + HART_SMITH dont le coefficient C1 est thermo-dependant
    + HYSTERESIS_3D avec une thermo-dependance via une loi des melanges

On choisit la temperature 20C (donc C1=1 et mu=1, Q0=1)

 1) dans le calcul Herezh, on impose completement la cinematique sur un cube unitaire :
    blocage de type traction aux 3 symetries et deplacement impose selon X, Y et Z sur les autres faces
      (UX different de UY different de UZ)
    => le calcul converge donc immediatement en un seul increment une seule iteration

    le repertoire calcul_HZ_reference/ contient ce calcul de reference Herezh


 2) dans le cas de l UMAT, on passe par le programme fortran UMAT_HZ.f faisant appel a Herezh (selon les etapes ci-dessous).
    cette fois, on impose directement un gradient de transfo F sur un point materiel (egalement sur une 
    seule iteration). Dans le fichier UMAT_HZ.f, on se place dans le meme cas que le calcul
    Herezh i.e le gradient de transfo suivant :
                    | 1    0    0 |
           F(t)  =  | 0    1    0 |
                    | 0    0    1 |

                   | 1+UX    0     0   |
         F(t+dt) = |  0     1+UY   0   |
                   |  0      0    1+UZ |

  

----------------------------------------
Etapes pour simuler l UMAT :
----------------------------------------

1) configurer le repertoire MAT_HS_HYST_THERMO/ qui contient l UMAT Herezh
   => commande : ./configuration.pl
      (rq : ce script compile appelc_1.cc et cree les pipes dans
            le repertoire MAT_HS_HYST_THERMO/)

2) compiler fichier fortran
   => commande : source src_compil
      (rq : creation de l executable UMAT_HZ)


3) lancer Umat
   => dans autre terminal :
        se placer dans le repertoire MAT_HS_HYST_THERMO/ et lancer :
          HZppfast_Vn-1 -f UMAT_MAT_HS_HYST_THERMO.info


4) lancer simulation d UMAT Abaqus : 
   => dans terminal d origine : ./UMAT_HZ

    (rq : cet executable affiche les contraintes obtenues pour SIG11, SIG22 et SIG33)



----------------------------------------
Resultat :
----------------------------------------
  le resultat donne par l UMAT est NON CONFORME au resultat Herezh de reference
  donne par le calcul calcul_ref_MAT_HS_HYST_THERMO.info

  resultat Herezh au pti 1 :
    => SIG11 =   0.5817952589717
       SIG22 =  -0.3120538436448
       SIG33 =  -0.2697414153449
  resultat UMAT :
    => SIG11 =   0.68322530954213967
       SIG22 =  -0.36561577905139864
       SIG33 =  -0.31760953049074186
