Metafor

ULiege - Aerospace & Mechanical Engineering

User Tools

Site Tools


commit:2013:06_27

Commit 2013-06-27

Modifications

Initialisation des éléments de contact

Lors de l'initialisation des éléments de contact, on calcule la projection de chaque nœud esclave sur l'entité maître pour déterminer si le nœud est initialement en contact selon la convention utilisée (Bilatéral, Unilatéral positif ou négatif). Dans le cas échéant, on sauvegarde le GPstate et CSstate de contact dans celui qui correspond à la prise de contact (GP00 et CS00). Ainsi, lors du premier pas de temps ou d'une phase initiale d'équilibrage quasi-statique, nous avons une force tangentielle de contact, puisque nous connaissons le point de collement.

Désormais, l'initialisation des éléments de contact est indépendante de la méthode de contact utilisée.

Ainsi, j'ai réussi à éliminer cette portion de code dans les méthodes conservatives.

void
ConsistentTimeIntegration::elementsInitialisation()
{
    //IntegrationType ndynOld=(IntegrationType)metafor.getIntegerData(MDE_NDYN);
    //chge de meth pour l'initialisation des contact
    //metafor.setIntegerData(MDE_NDYN, QUASI_STATIQUE);
    metafor.elementsInitialisation();
    //metafor.setIntegerData(MDE_NDYN, ndynOld);
}

Cas-Test

J'ai ajouté dans la battery des cas tests 2D et 3D pour tester le contact collant bilatéral utilisée pour “coller” deux corps déformables avec une phase initiale d'équilibrage quasi-statique pour la mise en rotation, selon différents schémas d'intégration temporelle.

Cas Test : initialBalancingAube2to3Moyeu2D Cas Test : initialBalancingAube2to5Moyeu2D Cas Test : initialBalancingAube2to3Moyeu3D Cas Test : initialBalancingAube2to5Moyeu3D

Perspectives

Méthode du lagrangien augmenté :

  • Débugger la méthode
  • Nettoyage du code
  • Ajout de nouveaux critères d'arrêt
  • Implémentation des cas test de référence de l'article de J. C. Simo and T. A. Laursen, 1990
  • Mettre à jour l'aire de contact nodale au cours des augmentations
  • Nettoyer et implémenter des nouveaux extracteurs de contact
  • Regarder si il est possible de faire des GPState en fonction du matériau et de la méthode de contact utilisée. (Gestion de mémoire + clarté du code)

Fichiers ajoutés/supprimés

[r]:
[a]:

Tests ajoutés/supprimés

[r]:
[a]: oo_meta\apps\contactMec\initialBalancingAube2to3Moyeu2D.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to3Moyeu2DConsDiss.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to5Moyeu2D.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to5Moyeu2DConsDiss.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to3Moyeu3D.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to3Moyeu3DConsDiss.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to5Moyeu3D.py
[a]: oo_meta\apps\contactMec\initialBalancingAube2to5Moyeu3DConsDiss.py

Gaëtan WAUTELET 2013/06/27

commit/2013/06_27.txt · Last modified: 2016/03/30 15:23 (external edit)