Anomalie #342
Pb RESTART Calcul RD+dynamique_implicite
Description
Gérard,
Je cherche à modéliser la montée dynamique du ballon après le lâcher. Mon ballon est modélisé en éléments axisymétrique (40 bielles).
J'utilise l'algo combiner avec l'enchainement suivant :
t de 0->1 RD : le noeud bas est bloqué
t de 1->3 dynamique_implicite : le noeud bas est libéré (temps_fin CL)
J'ai ajouté une masse additionnelle au noeud bas pour représenté la masse nacelle.
Le calcul inextenso se déroule très bien. Par contre lorsque je fais un RESTART à l'incrément 1 le calcul dynamique diverge au 1er incrément. Pourquoi ?
Pistes :
1. J'ai suspecté l'algo combiner mai sj'ai le même comportement avec les 2 algo séparé.
2. J'ai incriminé la libération brutale du ballon par suppression de la CL, bien que le calcul tourne sans RESTART. Pour adoucir le lâcher j'ai appliqué au noeud bas la force de réaction puis je réduits cette force à zéro via un fct nD. Le calcul inextenso se déroule très bien et montre l'intérêt d'un lâcher "progressif". Lorsque je fais une reprise le calcul dynamique diverge quelques incréments après que la force commence à décroitre.
3. J'ai pensé que dans le RESTART entre RD et dynamique_implicte il y pouvait y avoir un pb. J'ai le même comprtement si je fais un RESTART depuis un incrément dynamique_implicite.
4. J'ai mis en cause l'algo implicite. Même comportement avec l'algo de dynamique explicite.
5. Lorsque je maintien la CL de blocage pendant la phase de calcul dynamique pas de divergence, ce qui est normal puisu'il ne se passe rien !
Question : y a t'il des grandeurs de calcul qui ne sont pas correctement relues ou initialisées lors du RESTART ? Ce problème de RESTART est bloquant pour la suite car l'objectif est de faire un couplage avec un calcul CFD pour mettre à jour les forces de pression aérodynamique. Je peux essayer de passer par une fonction nD externe (donc sans RESTART) mais je voudrais au moins pouvoir démarrer après la fin du calcul statique d'équilibre au sol du ballon avant le lâcher (qui dure 1h30 !).
Voilà une semaine que je bataille avec ce problème :-(
Merci
Frank
Ci-joint mon calcul élémentaire (10 sec).
Fichiers