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
Mis à jour par Frank Petitjean il y a plus d'un an
Attention, dans le fichier .info jiont le noeud bas du ballon est toujours bloqué (N_deb UY #TEMPS_MAXI= 1)
Mis à jour par Frank Petitjean il y a plus d'un an
- Fichier champ_force.txt champ_force.txt ajouté
- Fichier maillage.her maillage.her ajouté
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
Gérard,
J'ai simplifié mon cas de calcul. Il n'y a que l'algo dynamique_explicite_tchamwa en lice. Lorsque je fais un RESTART les valeurs relues sont nulles et le calcul diverge.
Il me semble qu'il y a vraiment un problème de relecture !
Ci-joint les fichiers
Frank
Mis à jour par Gérard Rio il y a plus d'un an
- Statut changé de Nouveau à En cours
Bonjour Frank,
je reprends ton exemple simplifié.
1) Je calcule en dynamique explicite de 0 à 1s --> ok (avec HZpp)
2) ensuite je fais un restart à l'incrément 1070 (cf. la ligne commentée de ton fichier) et le calcul reprend et s'arrête à 1s avec le même résultat qu'en 1)
3) je refais le même restart avec HZppfast et idem, j'obtiens le même résultat qu'en 1)
du coup, quelle est la donnée à modifier dans le .info pour obtenir une erreur ?
Mis à jour par Frank Petitjean il y a plus d'un an
- Fichier modele.log modele.log ajouté
- Fichier modele_princ.maple modele_princ.maple ajouté
- Fichier modele_princ1.maple modele_princ1.maple ajouté
- Fichier modele1.log modele1.log ajouté
Gérard,
Je procède comme toi avec la version HZppfast 7.008 et j'ai bien un problème de RESTART ! Ci-joint les 2 fichiers log et Maple.
Le comportement est le même avec la version HZpp.
Comment un RESTART peut-il dépendre de la machine ? La seule situation où le restart fonctionne, et que j'utilise depuis toujours, c'est en RD seul.
As-tu un moyen de faire tourner Herezh sur Linux ?
Frank
Mis à jour par Gérard Rio il y a plus d'un an
j'ai créer les exécutables linux de la version actuelle (7.009) et j'ai testé le restart sous linux: dans mon cas cela fonctionne sans pb.
Peux-tu regarder si cette version répond à ton pb par contre a priori je n'ai pas modifié la partie restart de l'algo ...
Mis à jour par Frank Petitjean il y a plus d'un an
Bonne nouvelle, le RESTART fonctionne avec la v7.010. Par contre que s'est-il passé avec la version V7.008 !!!
Je l'ai réinstallée et de nouveau le RESTART ne fonctionne pas. Mystère que je trouve un peu inquiétant quand même...
Donc c'est ok pour moi, merci.
Tu peux fermer le ticket.
Frank
Mis à jour par Gérard Rio il y a plus d'un an
- Statut changé de En cours à Résolu
- % réalisé changé de 0 à 100
ok je ferme le ticket
NB: également mystère pour moi...