Herez++_test

bonjour Frank,

éléments de réponse:
---------------------

- Au niveau de la loi de comportement, je suppose qu'il s'agit de la même loi
que dans le précédent test.

- a priori le résultat correspond à ce que j'imaginerais:

courbes_1.png : une sortie tous les 50 itérations, avec une ligne entre les points,
et avec maxi 1000 incréments donc 1.01 s -> le résultat paraît logique ...

courbes_2.png : une sortie tous les 200 itérations (car c'est le nombre de sortie qui
pénalise ici le temps global de calcul !), ici jusqu'à 1.18 s -> là on commence à voir
premiers rebonds. On voit aussi que les premières hautes fréquences ont disparu

courbes_3.png : une sortie tous les 1000 itérations, donc cela donne un filtrage
numérique, mais néanmoins on n'a plus d'oscillations importantes, seulement une
fréquence de base qui va en s'atténuant. Ici le calcul va jusqu'à environ 1.9s et le
temps global calcul + sortie : 1:58.67 , bon... le résultat me paraît assez logique.
On a l'impression que si l'on attend encore un peu on aura une stabilisation du
déplacement autour de -0.018mm ? (à la louche en regardant le graphique !)

@+,

NB: syntaxe du dernier ordre pour le temps:
time calcul_script

---------- réponse 2 --------
maintenant j'utilise une loi beaucoup plus simple :
#----------------------------
sangle1 LOI_ADDITIVE_EN_SIGMA

ISOELAS1D
  1.5e9 0.3

MAXWELL1D
    E= 1.5e11     mu=      10e6 type_derivee  -1
                fin_coeff_MAXWELL1D

fin_liste_lois_elementaires

et là, on a une stabilisation beaucoup plus rapide ... vers le même état a priori
c'est très étonnant ! -> courbe4.png

Nouvel essai avec DFC à la place de Tchamwa
au observe de nouveau des oscillations ==> donc le filtrage précédent était dû à Tchamwa
qui filtrait bigrement. Cependant la loi élastique que j'ai utilisée n'est pas exactement
identique à la loi non linéaire hypo-élastique d'où une fréquence d'oscillation plus
grande. -> courbe5.png

Autre essai avec Tchamwa + un amortissement numérique plus faible

dynamique_explicite_tchamwa
PARA_TYPE_DE_CALCUL
#-----------------------------------
phi= 1.01

--> on obtient toujours le même résultat: en gros ici l'algo de tchamwa filtre totalement
les fréquences, ce qui fait que l'on a une stabilisation très rapide, mais ... sans
doute pas physique ??

autre essai avec DFC et un temps plus long 4s -> 1:12.72 de calcul avec un pas de temps
0.4 le pas critique
là on obtient une stabilisation correcte à 4s -> courbe5.png
j'ai ici plus confiance dans ce dernier résultat, compte tenu du cas traité, mais il
faudrait sans doute plus approfondir

c'est intéressant, affaire à suivre !

Gérard Rio a écrit :

---------- réponse 2 --------
maintenant j'utilise une loi beaucoup plus simple :
#----------------------------
sangle1 LOI_ADDITIVE_EN_SIGMA

ISOELAS1D
1.5e9 0.3

MAXWELL1D
E= 1.5e11 mu= 10e6 type_derivee -1
fin_coeff_MAXWELL1D

fin_liste_lois_elementaires

et là, on a une stabilisation beaucoup plus rapide ... vers le même état a priori
c'est très étonnant ! -> courbe4.png

Nouvel essai avec DFC à la place de Tchamwa
au observe de nouveau des oscillations ==> donc le filtrage précédent était dû à Tchamwa
qui filtrait bigrement. Cependant la loi élastique que j'ai utilisée n'est pas exactement
identique à la loi non linéaire hypo-élastique d'où une fréquence d'oscillation plus
grande. -> courbe5.png

Autre essai avec Tchamwa + un amortissement numérique plus faible

dynamique_explicite_tchamwa
PARA_TYPE_DE_CALCUL
#-----------------------------------
phi= 1.01

--> on obtient toujours le même résultat: en gros ici l'algo de tchamwa filtre totalement
les fréquences, ce qui fait que l'on a une stabilisation très rapide, mais ... sans
doute pas physique ??

autre essai avec DFC et un temps plus long 4s -> 1:12.72 de calcul avec un pas de temps
0.4 le pas critique
là on obtient une stabilisation correcte à 4s -> courbe5.png
j'ai ici plus confiance dans ce dernier résultat, compte tenu du cas traité, mais il
faudrait sans doute plus approfondir

c'est intéressant, affaire à suivre !

Voici le meilleur choix de paramètres pour un comportement viscoélastique nonlinéaire. Pour 4 sec de temps physique le temps de calcul est de quelques secondes seulement !
J'ai maintenant les 2 méthodes de modélisation pour les sangles de chaine de vol :
a) avec 1 seul élément et une loi complexe qui est désactivée en compression
b) avec un ensemble d'éléments biellettes et une loi standard.
Il me reste à comparer plus finement les résultats.