Differences

This shows you the differences between two versions of the page.

--- commit:2014:11_08 [2014/11/08 17:53] – created wautelet
+++ commit:2014:11_08 [2016/03/30 15:23] (current) – external edit 127.0.0.1
@@ Line 3: / Line 3: @@
 ===== Modifications =====
-==== Réorganisation de la libraire mtContact ====
+==== Matériau de contact  ====
-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.
+J'ai modifié la manière dont sont calculées les forces de contact tangentielles et la matrice de raideur matérielle dans le cas du Lagrangien Augmenté avec frottement afin d'accroître la robustesse du modèle et être capable de traiter le cas où la pression de contact vaut zéro et la prédiction du cisaillement de contact est différente de zéro (on se trouve au sommet du Cône de Coulomb d'où la difficulté de définir la direction de glissement !).
-==== Modification Interfaçage du Lagrangien Augmenté  ====
+==== Modification Interfaçage du matériau de contact de Tresca  ====
-Par rapport à mon commit précédent, j'ai raccourci la fonction :
+J'ai modifié l'interfaçage python des paramètres matériaux afin qu'il correspond mieux à ce qui est implémenté dans le code source. J'ai adapté tous les cas-tests et j'ai mis à jour la documentation !
-<code>
-alm.addContactInteraction(ci) devient alm.add(ci)
-</code>
-==== Ajout d'une fonctionnalité pour la méthode Area In Contact ====
+==== Nettoyage de l'ancien interfaçage ====
-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 :
-<code>
-prp2 = ElementProperties(Contact2D|3DElement)
-prp2.put(AREAINCONTACTBOUNDARYSELECTION, AICBS_INTERNAL)
-</code>
-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 !
-<note important>
-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>.
-</note>
-<note warning>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 ! </note>
-==== 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):
-<code>
-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)
-</code>
+J'ai nettoyé les anciennes fonctions appelées par l'algorithme du Lagrangien Augmenté dans les librairies mtElements, mtFEM, mtFEMBase. Il reste toutefois quelques fonctions à nettoyer dans les méthodes de contact et dans les matériaux de contact, mais c'est localisé dans la libraire mtContact !
 ==== Perspective ====