J'ai ajouté une option globale de Debug accessible dans Metafor Options dans l'interface graphique :

L'option “save values to disk at each iteration” qui est appelée à chaque fois qu'il y a un calcul des forces.

Dès lors, tous vos extracteurs extraient la valeur à extraire à chaque pas de temps et à chaque itération !

Cette option est désactivée par défaut !

Dans mon cas de figure, je peux regarder l'évolution des forces de contact au cours des augmentations du Lagrangien Augmenté pour observer si l'algorithme se comporte correctement.

J'ai proposé une hiérarchie des matériaux de contact et des GPStates associés. J'en ai profité pour nettoyer toutes les fonctions membres des matériaux. C'est plus facile à gérer maintenant et il est plus facile d'ajouter une nouvelle loi matérielle de contact.

J'ai nettoyé la mise en forme des différents fichiers dans les sources. Il reste à corriger les fichiers suivants :

 '#include "gisoSimplifyMesh.h"' not the first file included in [F:\Offi1\oo_meta\geniso\src\gisoSimplifyMesh.cpp]
 no EOL at the end of the file in [F:\Offi1\oo_meta\mtElements\boundaries\Contact2DElement.cpp]
 no EOL at the end of the file in [F:\Offi1\oo_meta\mtElements\boundaries\Contact3DElement.cpp]
 RESOURCE_H not found 3x in [F:\Offi1\oo_meta\mtMain\resource.h]
 '#include "mtMain.h"' not found in [F:\Offi1\oo_meta\mtMain\resource.h]
 RESOURCE_H not found 3x in [F:\Offi1\oo_meta\stp2py\resource.h]
 UNISTD_H not found 3x in [F:\Offi1\oo_meta\stp2py\win32\unistd.h]

• Allocation de la méthode de traitement du contact comme dans les éléments finis volumiques avec les méthodes d'intégration
• Isoler le calcul de la force de contact dissipative et sa contribution à la matrice de raideur de l'élément de contact
• Ajouter la contribution géométrique à la matrice de raideur de contact dans le cas Rigide-Déformable et interfacer l'option via l'element properties
• Placer les variables stockées dans les GPStates de contact au bon endroit (pour la gestion de la mémoire)
• Simplifier l'implémentation de la méthode du Lagrangien augmenté
• Proposer un nouvel interface d'utilisation de la méthode du Lagrangien augmenté et également si possible des méthodes de contact (ce n'est pas intuitive du tout à l'heure actuelle ….)
• Ajoutez des cas-tests simples pour tester la méthode du Lagrangien augmenté 2D/3D
• Débugger la méthode du Lagrangien augmenté

[r]:mtFEM/ContactGpState.h => [a]:mtMaterials/ContactGpState.h
[r]:mtFEM/ContactGpState.cpp => [a]:mtMaterials/ContactGpState.cpp
[r]:mtFEM/ContactGpState.inl
[r]:mtFEM/ContactState.h => [a]:mtElements/ContactState.h 
[r]:mtFEM/ContactState.cpp => [a]:mtElements/ContactState.cpp 
[r]:mtFEM/ContactState.inl => [a]:mtElements/ContactState.inl 
[r]:mtFEM/ContactMaterial.h => [a]:mtMaterials/ContactMaterial.h 
[r]:mtFEM/ContactMaterial.cpp => [a]:mtMaterials/ContactMaterial.cpp 
[r]:mtFEM/ContactMaterial.inl 
[r]:mtFEM/TypeOfContact.h => [a]:mtMaterials/TypeOfContact.h 
[r]:mtFEM/TypeOfContact.cpp => [a]:mtMaterials/TypeOfContact.cpp 
[a]:mtXFEM/XFEMContactElShcuts.h
[a]:mtXFEM/XFEMContactElShcuts.cpp
[a]:mtMaterials/CoulombContactGpState.h
[a]:mtMaterials/CoulombContactGpState.cpp
[a]:mtMaterials/TrescaContactGpState.h
[a]:mtMaterials/TrescaContactGpState.cpp
[a]:mtMaterials/StickingContactGpState.h
[a]:mtMaterials/StickingContactGpState.cpp
[a]:mtMaterials/FrictionlessContactGpState.h
[a]:mtMaterials/FrictionlessContactGpState.cpp
[a]:mtMaterials/StickingVariablePenaltyContactGpState.h
[a]:mtMaterials/StickingVariablePenaltyContactGpState.cpp
[a]:mtMaterials/FrictionlessVariablePenaltyContactGpState.h
[a]:mtMaterials/FrictionlessVariablePenaltyContactGpState.cpp
[a]:mtMaterials/FEMContactGpState.h
[a]:mtMaterials/FEMContactGpState.cpp
[a]:mtMaterials/FrictionalMechanicalContactGpState.h
[a]:mtMaterials/FrictionalMechanicalContactGpState.cpp
[a]:mtMaterials/NormalMechanicalContactGpState.h
[a]:mtMaterials/NormalMechanicalContactGpState.cpp
[a]:mtMaterials/MechanicalContactGpState.h
[a]:mtMaterials/MechanicalContactGpState.cpp
[a]:mtMaterials/TangentialMechanicalContactGpState.h
[a]:mtMaterials/TangentialMechanicalContactGpState.cpp
[a]:mtMaterials/FEMContactMaterial.h
[a]:mtMaterials/FEMContactMaterial.cpp
[a]:mtMaterials/TangentialMechanicalMaterial.h
[a]:mtMaterials/TangentialMechanicalMaterial.cpp
[a]:mtMaterials/MechanicalMaterial.h
[a]:mtMaterials/MechanicalMaterial.cpp
[a]:mtMaterials/NormalMechanicalMaterial.h
[a]:mtMaterials/NormalMechanicalMaterial.cpp
[a]:mtMaterials/FrictionalMechanicalMaterial.h
[a]:mtMaterials/FrictionalMechanicalMaterial.cpp

[r]:
[a]:

— Gaëtan WAUTELET 2013/08/28