Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2009:12_10

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

A Hill48PlasticCriterion2.h/cpp
R

Tests ajoutés/supprimés

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

commit/2009/12_10.txt · Last modified: 2016/03/30 15:23 (external edit)