diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README index 8e54e6f..2e5b85f 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/README +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_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, figures dans le Readme ------------------------------------------------------ Mots-cles @@ -19,12 +23,33 @@ test simple du mode de calcul : dynamique explicite (DFC) TYPE_DE_CALCUL dynamique_explicite "" +L'objectif est ici d'observer le fonctionnement de l'algorithme classique DFC. +On observe le déplacement, vitesse, accélération transversalement au sens imposé. +Ces mouvements dépendent donc ici du coefficient de poisson, c-a-d -0.3 ce qui +se passe dans le sens imposé. + +Remarques: On observe l'apparition de fréquences dites numériques. +Pour supprimer ou atténuer ces fréquences numériques on pourra utiliser: + 1) soit des schémas numériques adaptés, type : + - en explicite Tchamwa et Chung lee, + - en implicite HHT, + 2) ou une loi de comportement modifiée par l'ajout d'un bulk viscosity. + +(cf. documentions théoriques et utilisateur d'herezh) + ------------------------------------------------------ Description du calcul ------------------------------------------------------ traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE) - loi elastique ISOELAS + - chargement suivant une vitesse constante (une rampe de déplacement) + +\figures: acceleration.pdf vitesse.pdf \legende: La première figure montre l'évolution de +l'accélération gamma2 du noeud 8. On observe que globalement on a bien une vitesse +transversale globalement = -0.3 la vitesse imposée. On observe également des oscillations +numériques, ce qui est classique en DFC.. La figure 2 montre le même type +d'évolution pour la vitesse \fin_legende ------------------------------------------------------------- Grandeurs de comparaison diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf new file mode 100644 index 0000000..e274cd2 Binary files /dev/null and b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/acceleration.pdf differ diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info index ca7a623..a90c94b 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite.info @@ -1,7 +1,7 @@ dimension 3 -niveau_commentaire 3 +niveau_commentaire 0#3 TYPE_DE_CALCUL ---------------------------------------------- @@ -66,7 +66,12 @@ para_pilotage_equi_global ---------------------------------------------- para_syteme_lineaire ---------------------------------------------- -para_affichage ---------------------------------------------- +para_affichage #------------ +#----------------------------------# +# PARAMETRE | VALEUR # +#----------------------------------# +FREQUENCE_AFFICHAGE_INCREMENT 100000 +FREQUENCE_AFFICHAGE_ITERATION 100000 resultats ---------------------------------------------- COPIE 0 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 index 969b3c3..63fadb1 100644 --- a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/cube_dynamique_explicite_princ.maple.ref1 @@ -1,7 +1,7 @@ #fichier au format maple6 ############################################################################################### -# Visualisation elements finis : Herezh++ V6.687 # -# 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.970574516706e-01 9.970574516706e-01 -3.777977376754e+00 -4.888446928246e+05 1.000000000000e-03 2.134381140592e-01 7.863544198951e-01 7.863544198951e-01 9.851975296540e-03 -2.955587360158e-03 -2.955587360161e-03 2.068916079252e+03 + 1.000073916208e-03 1.010000739162e+00 9.970574516706e-01 9.970574516706e-01 -5.723375540897e+00 -2.061055598212e+07 1.000073916208e-03 2.134382702626e-01 7.863544198951e-01 7.863544198951e-01 9.852692719186e-03 -2.955587360158e-03 -2.955587360161e-03 2.069118889116e+03 diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu new file mode 100644 index 0000000..a45fbcd --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/plotPdf.gnu @@ -0,0 +1,41 @@ +set term X11 font "Bold,20" +show output +# set key 0.02,1000. + set xtics 1. + set ytics 1. +set size 0.98, 1. +set terminal pdf lw 2 font "Times-Roman,10" #"Helvetica" 14 + set key center top +set ylabel "Acceleration / 10000000 " +set xlabel "temps x 1000000" +set xrange [0.:10] +#set yrange [-6.:6.] +##set y2range [0:9.] +##set y2label "Axial strain [%]" +set grid + set output "acceleration.pdf" +# set title " Strain = 2 x (strain_max in log)-0.0035 " font "Times-Roman,10" + show title + set style data line +plot 'cube_dynamique_explicite_princ.maple' u ($1*1000000.):($6/10000000.) \ +t 'acceleration ' lt 1 lw 2, \ +'cube_dynamique_explicite_princ.maple.ref1' u ($1*1000000.):($6/10000000.) \ + t 'acceleration ref' lt 6 lw 1\ + +#--------------- + set output "vitesse.pdf" +# set title " Strain = 2 x (strain_max in log)-0.0035 " font "Times-Roman,10" + show title +set autoscale +set xrange [0.:10] +set yrange [-14.:8.] + set style data line +plot 'cube_dynamique_explicite_princ.maple' u ($1*1000000.):($5/1.) \ +t 'vitesse ' lt 1 lw 2, \ +'cube_dynamique_explicite_princ.maple.ref1' u ($1*1000000.):($5/1.) \ + t 'vitesse ref' lt 6 lw 1\ + + + +# pause -1 "Hit return to continue" + diff --git a/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf new file mode 100644 index 0000000..21ae96a Binary files /dev/null and b/Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_explicite/vitesse.pdf differ