Anomalie #210
Ajouté par Gérard Rio il y a plus de 5 ans.
Mis à jour il y a plus de 5 ans.
Description
Bonjour Gérard,
J'observe un changement de comportement étrange entre la version 6.894 et les suivantes, un effet de bord assez incompréhensible aux vues des évolutions apportées au logiciel.
Dans mes calculs BSO je réalise un dernier incrément de calcul dans pour lequel tous les noeuds sont encastrés. Cela me permet de récupérer proprement la force ascensionnelle globale. Cela fonctionne depuis 4 ans !
Après la version 6.894 la force de réaction (N_tout) est nulle. J'ai testé 6.899 et non la 6.898.
Autre chose qui n'a je pense rien à voir, avec la version 6.903 le calcul plante sérieusement lors du passage au 2e incrément.
Fichiers
Effectivement une modification consécutive au contact, a laissée une petite trace ! dans le cas de l'algo de relaxation dynamique. Du coup les torseurs de réaction n'étaient pas correctement calculés.
C'est modifié et ce sera opérationnel dans les versions à partir de 6.907
J'ai également changé la prise en compte du critère mixte de convergence:
ARRET_A_EQUILIBRE_STATIQUE_ 2
qui ne fonctionnait pas en fait, correctement ici: le delta X mini n'était pas bien vérifié: cf. le paramètre max_deltaX_pourRelaxDyn_.
Du coup il faut faire attention à la valeur du paramètre : max_deltaX_pourRelaxDyn_
Dans l'exemple fourni on avait
max_deltaX_pourRelaxDyn_ 0
C'est difficile à obtenir !! et on observe que cela ne converge pas
Dans mon cas j'ai mis :
max_deltaX_pourRelaxDyn_ 0.1
là pas de pb de convergence
- Statut changé de Nouveau à Résolu
- % réalisé changé de 0 à 100
suite de la demande concernant la convergence:
1) si on s'arrête à t= 1s c-a-d incrément 2:
on constate que la version 6.894 et la version 6.907 conduisent aux mêmes résultats et que la version 6.907 est un peu plus lente (environ 1% de perte). Je pense que c'est dû au fait qu'il y a un ajout de vecteurs globaux (consécutifs au nouveau type de contact). Mais la différence est quand même à la marge.
2) si on continue jusqu'à t= 4s, avec les paramètres actuelles, effectivement on observe des relaxations successives avec la version 6.907 et on n'arrive jamais à la convergence. Du coup j'ai doublé arbitrairement le lambda final (de 2 -> 4). Là on observe une convergence avec la version 6.907 et surtout on s'aperçoit que le temps de calcul cette fois est beaucoup plus faible qu'avec la version 6.894 !! alors que le nombre d'itération est sensiblement le même !!!
6.894 avec lambda final= 2 -> cpu: 02:16.68
6.907 avec lambda final= 4 -> cpu: 00:59.08, c-a-d plus de 2 fois plus rapide
Je pense que cela provient d'une modif que j'ai fait lors de la résolution d'un système linéaire que j'ai constaté en fait non symétrique, alors qu'il était supposé symétrique (cf. wiki). Du coup les algos de CP CP2 et plis convergent plus rapidement et cela doit expliquer le gain obtenu
J'ai ensuite diminué le lambda final à 3 (2 -> 3) et j'obtiens:
avec la 6.894 avec lambda final= 3 -> cpu = 1:41.04 donc là aussi on voit qu'il est préférable d'augmenter un peu le lambda sinon on a des amortissements intempestifs
avec la 6.907 avec lambda final= 3 : cela converge également correctement, avec moins d'itération qu'avec la 6.894 (580 au lieu de 739 pour l'incrément 3) et on obtient un cpu = 00:37.33 donc quasiment un facteur 4 par rapport à la version 6.894 initiale !!
Donc à confirmer avec la version à venir : 6.907
Formats disponibles : Atom
PDF