Differences

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

--- commit:2015:10_29 [2015/10/29 11:22] – wautelet
+++ commit:2015:10_29 [2016/03/30 15:23] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
 ====== Commit 2015-10-29 ======
-==== Contact Bug Fix ====
+===== Time Step Manager =====
-=== Matrice de raideur de contact ===
+==== Nettoyage ====
-== Wire ==
+J'ai éliminé la variable stepIsRecomputed.
-Le calcul de la normale à un wire dépend de son orientation global et non de l'orientation de chacune des courbes le constituant. Dès lors, il fallait adapter la fonction du calcul de la normale !
+Pour information, un pas de temps est recalculé si il y a une erreur (Jacobien Négatif ou non convergence d'une loi matérielle (retour radial, etc ...) durant le calcul des forces (C'est le flag "ErrorDuringStep" situé dans Metafor) ou il est refusé (C'est le flag "StepIsRefused" situé dans TimeStepManager) sur base d'un critère d'évaluation de la qualité de l'intégration temporelle entre les deux instants discrets to et t1.
-== Augmented Lagrangien Method ==
+J'ai éliminé la variable oldTimeStepNo.
-J'ai corrigé un bug dans l'expression de la matrice de raideur de contact dans le cas collant ou avec frottement (Partie géométrique).
+Après investigations, cette variable était utilisée pour indiquer le step courant avec un pas de temps limité par la présence d'un archivage. Pour éviter de perturber l'intégration temporelle, on recalcule le pas de temps sur base du fait qu'il a été tronqué précédemment !
-=== Area in Contact ===
+A la place de cette variable, j'ai utilisé le flag suivant "timeStepIsLimitedByFac", plus explicite !
-J'ai adapté le calcul de la mise à jour des composantes de la force de contact lors de la phase ALE pour faire face au cas particulier de l'area in contact (mise à jour de l'aire nodale à toutes les itérations mécaniques).
+==== Restart ====
-==== Contact Augmented Lagrangien Method ====
-En suivant la procédure proposée par Romain pour la mise à jour des grandeurs de contact lors de la phase endALE(), j'ai proposé une méthode (**peu rigoureuse !**) pour mettre à jour les grandeurs liées au contact dans le cas du transfert des données d'un maillage à l'autre.
-Dans le cas du Lagrangien Augmenté, il faut transférer le lagrangien normal (scalaire) et le lagrangien tangentiel (vecteur) ainsi que le point de collement.
-Pour se faire, j'ai tout simplement utilisé la même méthode qu'en pénalité, c'est à dire j'ai adapté la valeur du lagrangien normal et du lagrangien tangentiel aux nouvelles valeurs de la pression de contact et du cisaillement de contact, si ceux-là étaient non nuls dans la configuration lagrangienne.
-<note important>La méthode fonctionne mais elle manque certainement de robustesse ! </note>
+J'ai corrigé un bug lors d'un restart d'un cas test. En fait, lors d'un restart, nous devons calculer le pas de temps du nouveau step. Il s'avère que si le pas de temps fut limité par la présence d'un archivage, le pas de temps était mal recalculé ! Ainsi, nous pouvons observer sur le cas test apps.complex.cont5_1 que le nombre de pas de temps pour l'intégration complète est de 76. Si nous recommençons l'intégration temporelle à la moitié, nous avons toujours 76 pas de temps au lieu de 62 pas de temps comme auparavant !
+Il y a toujours un problème pour la méthode de calcul du pas de temps basée sur l'erreur d'intégration.
-==== Nettoyage des cas tests ====
+===== Nettoyage des cas tests =====
 Pour tous les cas tests qui utilisent la méthode d'erreur d'intégration pour mettre à jour le pas de temps, j'ai indiqué explicitement dans le jeux de données la commande suivante :