From c78c1831515992cc63937cc9ca3459f6ddb56489 Mon Sep 17 00:00:00 2001 From: rio <> Date: Tue, 10 Jan 2017 10:08:00 +0000 Subject: [PATCH] mise a jour du test de Chung Lee cf. le readme --- .../README | 22 ++++++++++++++++ .../chung_lee_dynamique_explicite.info | 9 +++++-- ...g_lee_dynamique_explicite_princ.maple.ref1 | 26 ++++++++++++++++--- 3 files changed, 52 insertions(+), 5 deletions(-) diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/README b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/README index afddab8..bee716a 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/README +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_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 @@ -22,12 +26,29 @@ test simple du mode de calcul : dynamique explicite Chung-Lee PARA_TYPE_DE_CALCUL beta= 1.037 "" +L'algorithme proposé par Chung et Lee (cf. ref qui suit) est un schéma d'avancement +temporel qui permet d'atténuer les hautes fréquences numériques qui apparaissent +pendant le calcul. Deux grandes familles de méthodes sont disponibles: soit via une intervention +au niveau du schéma d'avancement temporel (algo de Chung Lee , Tchamwa, HHT etc.) +soit au niveau de la loi de comportement via par exemple le Bulk Viscosity. + + + J. Chung and J. Lee. A new family of explicit time integration methods for linear and non-linear structural dynamics. International Journal for Numerical Methods in Engineering, + 37:3961 3976, 1994. ------------------------------------------------------ Description du calcul ------------------------------------------------------ traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE) - loi elastique ISOELAS + - vitesse de traction uniforme + +\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 de Chung Lee +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 \fin_legende ------------------------------------------------------------- Grandeurs de comparaison @@ -39,3 +60,4 @@ pour le noeud 8 : pour le point d integration 1 de l element 1 : - deformations : EPS11 EPS22 EPS33 - contrainte : SIG11 + diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite.info b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite.info index d3391cc..894f9e8 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite.info +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite.info @@ -1,7 +1,7 @@ dimension 3 -niveau_commentaire 3 +niveau_commentaire 0#3 TYPE_DE_CALCUL ---------------------------------------------- @@ -69,7 +69,12 @@ para_pilotage_equi_global ---------------------------------------------- para_syteme_lineaire ---------------------------------------------- -para_affichage ---------------------------------------------- +para_affichage #------------ +#----------------------------------# +# PARAMETRE | VALEUR # +#----------------------------------# +FREQUENCE_AFFICHAGE_INCREMENT 100000 +FREQUENCE_AFFICHAGE_ITERATION 100 resultats ---------------------------------------------- COPIE 0 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite_princ.maple.ref1 index 7346203..8e54fec 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_dynamique_explicite_princ.maple.ref1 +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_chung_lee_dynamique_explicite/chung_lee_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 1.000000000000e+00 1.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.000000000000e-03 2.134381140592e-01 7.886751345948e-01 7.886751345948e-01 9.851975296540e-03 0.000000000000e+00 0.000000000000e+00 2.785077631906e+03 + 1.000159262703e-03 1.010001592627e+00 9.970564659791e-01 9.970564659791e-01 -2.886104243346e+00 1.108081251954e+02 1.000159262703e-03 2.134384506209e-01 7.863536425047e-01 7.863536425047e-01 9.853521080987e-03 -2.956059413922e-03 -2.956059413922e-03 2.069239426720e+03