• Hill 48 version 2

• 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

A Hill48PlasticCriterion2.h/cpp
R

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

— Luc PAPELEUX 2009/12/10