Herez++_test

Concernant le test de validation en attente: Test_R_contact_TYPE_MATRICE
mis en place par Julien.
J'ai mis des commentaires (*** GR ->) qui je propose de conserver d'où ce ticket à moi même

------------------------------------------------------
But du test
------------------------------------------------------
test purement "informatique" de l utilisation conjointe des fonctionnalités :
- contact
- para_syteme_lineaire $\rightarrow$ TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK
- para_syteme_lineaire $\rightarrow$ MATRICE_S_SECONDAIRE_S_
- para_syteme_lineaire $\rightarrow$ TYPE_RESOLUTION_S_SECONDAIRE_S_

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
--- à effacer du README ---

Sans contact, le calcul tourne (TYPE_MATRICE + MATRICE_S_SECONDAIRE_S_= OK)

Sans matrice secondaire, le calcul tourne (contact + TYPE_MATRICE = OK).

La conjonction des 3 ne tourne pas (contact + TYPE_MATRICE + MATRICE_S_SECONDAIRE_S_ = ECHEC)

La conjonction des 3 tourne avec la version 6.768

• GR -> il y avait eu une évolution non complètement terminée, c'est réparé dans la version 6.775

IMPORTANT : le .CVisu actuel fait une sortie .maple pour tous les noeuds des 2 maillages de diverses grandeurs
liées au contact et aux réactions. La plupart sont nulles. Peut-etre que certaines d entre-elles
ne devrait pas etre nulles ?
=> sorties à vérfier avant de générer un .maple.ref

• GR > il y a un message de warning concernant VECT_REAC_N
  a priori le message est normal.
  lorsque l'on fait une sortie en gmsh, ce vecteur doit être calculé pour tous les points donc il est
  définit en tout point et il n'y a plus de message de warning
  - lorsque l'on fait une sortie maple (sans gmsh) ce vecteur n'existe "que" pour les noeuds en contact
  du coup, s'il n'y a pas de contact, herezh prévient l'utilisateur qu'il n'y a pas de réaction, mais ce n'est
  qu'un warning

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

------------------------------------------------------
Description du calcul
------------------------------------------------------
compression uniaxiale suivant Y d un cube unitaire. La compression est appliquée par contact
avec un second maillage (également un maillage à un seul élément).
Les déplacements de ce second maillage sont entièrement gouvernés par les conditions limites
(cube infiniment rigide). Le cube déformable a des conditions de symétrie usuelles pour une compression uniaxiale.