This is an old revision of the document!

Commit 2014-11-08

Modifications

Réorganisation de la libraire mtContact

J'ai regroupé toutes les sources dans le nouveau sous-dossier src et j'ai ajouté les nouveaux cas tests qui testent l'interfaçage dans le nouveau sous-dossier tests.

Modification Interfaçage du Lagrangien Augmenté

Par rapport à mon commit précédent, j'ai raccourci la fonction :

            
alm.addContactInteraction(ci) devient alm.add(ci)

Ajout d'une fonctionnalité pour la méthode Area In Contact

Pour des cas-tests d'Yves, j'ai interfacé une option qui permet de calculer une aire de contact nodale lorsque le nœud de contact se trouve à l'intérieur du solide.

Pour se faire, il suffit d'utiliser les commandes suivantes :

prp2 = ElementProperties(Contact2D|3DElement)   
prp2.put(AREAINCONTACTBOUNDARYSELECTION, AICBS_INTERNAL)

Par défaut, l'option AREAINCONTACTBOUNDARYSELECTION est fixée à AICBS_EXTERNAL. Dans ce cas là, un noeud de contact à l'intérieur du solide possède une aire de contact nodale nulle et dès lors une force de contact nulle, quelque soit la valeur du gap !

J'ai corrigé un bug dans mes sources précédentes qu'Yves a mis en évidence lors d'une tentative de compilation sur le cluster :

La fonction valeur absolue “fabs” (pour les doubles ou float) nécessite d'inclure dans le header #include <math.h>.

J'ai ajouté la conversion des fichiers .txt dans l'opération dos2unix dans comp.py sur les stations ! Veuillez mettre à jour votre bin sur toutes les stations !

Modification des Cas Tests

Area In Contact

J'ai nettoyé tous les cas-tests où l'option Area In Contact est utilisée sous les formes suivantes (avec un booléen ou un entier au lieu de l'enum):

prp2 = ElementProperties(Contact2D|3DElement)        
prp2.put(AREAINCONTACT, True) => prp2.put(AREAINCONTACT, AIC_ONCE)
prp2.put(AREAINCONTACT, False) => prp2.put(AREAINCONTACT, AIC_NO)
prp2.put(AREAINCONTACT, 1) => prp2.put(AREAINCONTACT, AIC_ONCE)
prp2.put(AREAINCONTACT, 0) => prp2.put(AREAINCONTACT, AIC_NO)

Perspective

Dans la suite du commit, il faut

ajouter des cas-tests (état plan déformation, axisymétrique et 3D) pour tester l'interfaçage et éventuellement statuer sur le sort des cas-tests actuellement dans la batterie (2ième étape);
ajouter les cas-tests de référence de l'article de Simo et Laursen 1990 (3ième étape);
faire la documentation;

Fichiers ajoutés/supprimés

[a]:
[r]:

Tests ajoutés/supprimés

[r]:
[a]:

— Gaëtan WAUTELET 2014/11/08

Metafor

Table of Contents