From 949c4028b995e362a7f2d241920d07054f34d7df Mon Sep 17 00:00:00 2001 From: rio <> Date: Thu, 12 Jan 2017 20:02:21 +0000 Subject: [PATCH] =?UTF-8?q?=20mise=20=C3=A0=20jour=20du=20test=20de=20Zhai?= =?UTF-8?q?=20(cf.=20readme),=20initialement=20l=20algorithme=20ne=20fonct?= =?UTF-8?q?ionnait=20pas.=20La=20methode=20originale=20propose=20par=20Zha?= =?UTF-8?q?i=20(prediction=20-=20correction)=20a=20ete=20mise=20a=20jour,?= =?UTF-8?q?=20et=20les=20resultats=20sont=20surprenants=20et=20interessant?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../Test_R_zhai_dynamique_explicite/README | 27 ++++++++++++++++++- .../zhai_dynamique_explicite.info | 16 ++++++----- .../zhai_dynamique_explicite_princ.maple.ref1 | 26 +++++++++++++++--- 3 files changed, 58 insertions(+), 11 deletions(-) diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/README b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/README index 9c0d1f3..5767c56 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/README +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/README @@ -2,6 +2,10 @@ Auteur ------------------------------------------------------ Julien Troufflard (troufflard@univ-ubs.fr) +Gérard Rio (gerard.rio@univ-ubs.fr) +NB: GR: modification du test janvier 2017: +- remplacement du maple.ref1 complet: les données initiales étaient fausses +- inclusion de commentaires, ref, figures dans le Readme ------------------------------------------------------ Mots-cles @@ -20,14 +24,35 @@ test simple du mode de calcul : dynamique explicite Zhai dynamique_explicite_zhai PARA_TYPE_DE_CALCUL - phi_minus= 0.5 grand_phi= 0.5 + phi_minus= 0.5 grand_phi= 0.16667 gamma= 0.5 beta= 0.16667 "" +L'algorithme proposé par Zhai (cf. ref qui suit) est un schéma d'avancement +temporel de type prédiction-correction, qui permet d'atténuer les hautes fréquences numériques qui apparaissent +pendant le calcul. Au final, ce schéma demande environ 2 fois plus de temps cpu qu'un +calcul classique DFC de plus le rayon de stabilité est en général plus faible. +Par exemple avec les paramètres de ce test, le pas de temps critique est environ de 5% plus +faible qu'en DFC. Par contre l'amortissement est efficace et on n'observe pas d'over-shoot +au début de l'application des efforts contrairements aux schémas d' Chung Lee et Tchamwa. + +W.-M. ZHAI. Two simple fast integration methods for large-scale dynamic problems in engineering. + International Journal for Numerical Methods in Engineering, + 39(24):4199 4214, December 1996. + ------------------------------------------------------ Description du calcul ------------------------------------------------------ traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE) - loi elastique ISOELAS + +\figures: acceleration.pdf vitesse.pdf \legende: La première figure montre l'évolution de +l'accélération gamma2 du noeud 8 en début de chargement c'est-à-dire juste après +le changement de vitesse et d'accélération. On observe bien que le schéma +conduit à un amortissement très rapide des oscillations numériques, qui demeurent +présentes dans le cas du schéma classique DFC. La figure 2 montre le même type +d'évolution pour la vitesse. En comparaison avec Chung Lee, on observe un amortissement +assez semblable. Par contre pour Zhai, il n'y a pas d'over-shoot en début d'application +du déplacement imposé. \fin_legende ------------------------------------------------------------- Grandeurs de comparaison diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite.info b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite.info index c37977a..e4291b8 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite.info +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite.info @@ -1,7 +1,7 @@ dimension 3 -niveau_commentaire 3 +niveau_commentaire 0# 3 TYPE_DE_CALCUL ---------------------------------------------- @@ -9,7 +9,7 @@ TYPE_DE_CALCUL ---------------------------------------------- dynamique_explicite_zhai avec plus lectureCommandesVisu PARA_TYPE_DE_CALCUL - phi_minus= 0.5 grand_phi= 0.5 + phi_minus= 0.5 grand_phi= 0.16667 gamma= 0.5 beta= 0.16667 < cube_1elt_npti8.her @@ -55,12 +55,9 @@ N_gauche 'UX= COURBE_CHARGE: rampe ECHELLE: 0.01' controle ---------------------------------------------- TEMPSFIN 0.001 -DELTAt COEF_PASCRITIQUE 0.95 -ITERATIONS 15 -PRECISION 1e-3 +DELTAt COEF_PASCRITIQUE 0.9 SAUVEGARDE DERNIER_CALCUL MAXINCRE 99999999 -NORME Residu/Reaction_et_VarRes para_pilotage_equi_global ---------------------------------------------- @@ -69,7 +66,12 @@ para_pilotage_equi_global ---------------------------------------------- para_syteme_lineaire ---------------------------------------------- -para_affichage ---------------------------------------------- +para_affichage #------------ +#----------------------------------# +# PARAMETRE | VALEUR # +#----------------------------------# +FREQUENCE_AFFICHAGE_INCREMENT 100 +FREQUENCE_AFFICHAGE_ITERATION 100000 resultats ---------------------------------------------- COPIE 0 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite_princ.maple.ref1 index cdfea62..23af7f9 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite_princ.maple.ref1 +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_zhai_dynamique_explicite/zhai_dynamique_explicite_princ.maple.ref1 @@ -1,7 +1,7 @@ #fichier au format maple6 ############################################################################################### -# Visualisation elements finis : Herezh++ V6.713 # -# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# Visualisation elements finis : Herezh++ V6.785 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # # http://www-lg2m.univ-ubs.fr # ############################################################################################### @@ -15,6 +15,20 @@ # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # # puis pour chaque maillage # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees @@ -47,6 +61,12 @@ #---------------------------------- torseurs de reactions -------------------------------- #0 0 (nombre de torseurs et nombre total de grandeurs associees) # +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# #---------------------------------- grandeurs aux noeuds -------------------------------- #1 2 (nombre de noeuds, nombre total de grandeurs associees) # noeud_8 [2]X [3]Y [4]Z [5]V2 [6]GAMMA2 ; @@ -77,4 +97,4 @@ # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs # etc. pour tous les points d'integration - element - 1.000000000000e-03 1.010000000000e+00 9.970583790889e-01 9.970583790889e-01 -4.259581600240e-01 -8.796369095293e+07 1.000000000000e-03 2.134381140592e-01 7.863551513269e-01 7.863551513269e-01 9.851975296540e-03 -2.954651707768e-03 -2.954651707765e-03 2.069142795025e+03 + 1.000012854150e-03 1.010000128541e+00 9.970574093877e-01 9.970574093877e-01 -2.886140462174e+00 1.107763349258e+02 1.000012854150e-03 2.134381412233e-01 7.863543865476e-01 7.863543865476e-01 9.852100057629e-03 -2.955630018474e-03 -2.955630018473e-03 2.068941011815e+03