===== Commit 2009-12-10 =====

  * Hill 48 version 2

===== Modifs =====

  * ** Hill 48 version 2**
     * Ecriture d'un critère de Hill48 avec le tenseur normal normé (pour pouvoir faire du cinématique)
     * le tenseur normal étant très non linéaire aux contraintes, le calcul d'une nouvelle contrainte nécessite un NR avec au besoin un line-search (procédure totale chère!!!)
     * pas encore de matrice de raideur analytique (utiliser la semi-num voir ci dessous)
    
  * ** Matrice de raideur numérique hypomaterial **
     * Pour sa thèse & débugger les matériaux, PP avait programmé une matrice de raideur numérique au niveau du matériau evpIsoHHypoMaterial. 
     * Cette méthode a été nettoyée (fallait un ptit coup de Karcher) et généralisée à tous les HypoMaterials.
     * Accessible via paramètre matériau à condition que l'élément soit intégré analytiquement (je dois encore nommer le paramètre): 

    prp.put(STIFFMETHOD, STIFF_ANALYTIC)
    mater.put(MATERIALSTIFFMETHOD, 1)


  * ** Precision de l'intégration des matériaux (hypomaterial) **
     * CONSTITUTIVE_LAW_PREC : précédemment : 1.0e-8, maintenant 1.0e-12
     * ajout d'une précision sur l'endo (ancienne valeur): DAMAGE_PREC  1.E-8
     * modification dans la normation du critère de convergence de la loi constitutive : afin de prendre en compte:  
       * les matériaux "normaux"
       * les calculs dans des unités moins usuelles 
       * les matériaux à adoucissement 
       * les matériaux simulés avec une loi de hollomon en remplacement de l'écrouissage 
       * les matériaux visco-plastiques à limite élastique nulle
       * ...
     * Le paramètre de normation du résidu est pris égal à la moyenne de la limite élastique initiale et du cirère (J2) du prédicteur élastique.
     * GK.initialise ne remet plus le Gamma à zero. On recommence donc le calcul de la plasticité à partir de la solution de l'ité / du pas de temps précédent (si il est admissible). Ce qui accélère la convergence (intérêt important pour les matrices de raideur numérique)
     * On fait toujours au moins une ité NR pour la plastification (suite à la modif précédente, on ne différentiait pas les pertu de la matrice de raideur).

  * ** Divers :  **
     * correction stat.py (merci robo & pp : quand les grands esprits se rencontrent...)
     * CMVIM : Unification de l'accès au matériau dans ''fillMechanicalAnalyticalStiffness''
     * toolbox/tetgen.py : correction en provenance de Marlène
===== Fichiers ajoutés/supprimés =====
<code>
A Hill48PlasticCriterion2.h/cpp
R
</code>

===== Tests ajoutés/supprimés =====
<code>
A apps/monosMaterials/evpIsoH3dTracRotHill48v2.py
A apps/monosMeca/tracRotHill48AxesSemiNum.py
A apps/monosMeca/tracRotHill48MatSemiNum.py
A apps/monosMeca/tracRotHill48SemiNum.py
A apps/monosMeca/tracRotHill48v2AxesSemiNum .py
A apps/monosMeca/tracRotHill48v2MatSemiNum.py
A apps/monosMeca/tracRotHill48v2Num.py
A apps/monosMeca/tracRotHill48v2SemiNum.py
A apps/qs/cont2MatAna.py
A apps/qs/cont2MatNum.py
A apps/qs/godetNumisheetHill48v1.py
A apps/qs/godetNumisheetHill48v2.py
R
</code>



 --- //[[L.Papeleux@ulg.ac.be|Luc PAPELEUX]] 2009/12/10 //