Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2019:08_30

Commit - LPx - 30 Aout 2019

AEJ : Correction Matrice de raideur analytique (hypo)

AEJ (suite) :

AejElement_mec.hpp :

  • MRT Numerique : j'avais commentarisé les portions de code relatives aux EqualityDofsConstraintSet du calcul des matrices de raideur tangentes numériques des AEJ. Comme ca n'a pas d'effet sur le calcul de la MRT numérique, je n'ai réactivé (pour le moment où on voudra tester cette possibilité)

TriangleAej2DElement.cpp & TetraAej3DElement.cpp :

  • J'avais laissé quelques détections de “pointOnBorder” dans la fonction updateAejProperties(), alors que j'ai stoqué l'information dans l'élément (⇒ j'utilise la variable stoquée au lieu de redétecter ….)

AEJ 2D vs 3D :

  • Lorsque Vincianne a commencé sa thèse, le 2D et le 3D n'étaient pas unifiés comme maintenant (au niveau des Points de gauss, tenseurs matériaux,…) ⇒ les AEJ montraient des différences dans le formalisme et son implémentation. Or depuis, une unification du 2D et 3D a été faite (du calcul des jacobiens qui même en 2D sont des matr3 (avec le terme 33 calculé selon que l'on est en axisym, état plan de défo ou état plan de défo généralisé (2.5D)). Les particularisations de calcul 2D/3D au sein des éléments AEJ n'ayant plus de raison d'être, j'ai tout unifié sur le 3D
  • computeQFromH : correction du calcul de
    Qijkl = (1./3) * Hijmn dmn dkl - (2./3) sigij dkl 

    (le “-” était appliqué 2 fois)

AEJ :

  • a ce jour, il semble que ces corrections ont corrigé les MRT analytiques pour les matériaux Hypo (on a les même solutions et la même convergence sur les quelques tests pour lesquels on peut valider la convergence MRT Ana vs Num : qulques tests en plus).
  • Par contre ca améliore légèrement la convergence des pk2 (neoHookean), mais elle reste quand même très différente ⇒ to be continued ….

Membrane de Cooks 3D

  • ajout d'éléments sur l'épaisseur des tests de membrane de cooks3D avec les AEJ, sinon tous les noeuds sont considérés sur la frontière et aucun effet AEJ n'est alors pris en compte (nz = 3)
  • pour que les tests avec MRT numérique ne prennent pas 45', j'ai réduit le maillage nx=4 et ny=4 (au lieu de nx=8 et ny=8)

Divers :

Fichiers ajoutés/supprimés :

Added : 
Deleted : 
Moved : 

Tests ajoutés/supprimés

added : oo_meta/apps/aej/cooks2DepNum.py
added : oo_meta/apps/aej/cooks2DincNum.py
added : oo_meta/apps/aej/cooks2DneoHookeanNum.py
added : oo_meta/apps/aej/cooks2Dep_ep2Num.py
added : oo_meta/apps/cooksMembrane/cooks3DElastTetraAejNum.py
Deleted : 
Moved : 

Luc Papeleux 2019/08/29

commit/2019/08_30.txt · Last modified: 2019/08/29 17:51 by papeleux