2
0
Fork 0

mise a jour du test DFC cf. Readme

This commit is contained in:
rio 2017-01-10 10:59:05 +00:00
parent c78c183151
commit 8265f7f7af
6 changed files with 96 additions and 5 deletions

View file

@ -2,6 +2,10 @@
Auteur Auteur
------------------------------------------------------ ------------------------------------------------------
Julien Troufflard (troufflard@univ-ubs.fr) 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 Mots-cles
@ -19,12 +23,33 @@ test simple du mode de calcul : dynamique explicite (DFC)
TYPE_DE_CALCUL TYPE_DE_CALCUL
dynamique_explicite 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 Description du calcul
------------------------------------------------------ ------------------------------------------------------
traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE) traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE)
- loi elastique ISOELAS - 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 Grandeurs de comparaison

View file

@ -1,7 +1,7 @@
dimension 3 dimension 3
niveau_commentaire 3 niveau_commentaire 0#3
TYPE_DE_CALCUL ---------------------------------------------- TYPE_DE_CALCUL ----------------------------------------------
@ -66,7 +66,12 @@ para_pilotage_equi_global ----------------------------------------------
para_syteme_lineaire ---------------------------------------------- para_syteme_lineaire ----------------------------------------------
para_affichage ---------------------------------------------- para_affichage #------------
#----------------------------------#
# PARAMETRE | VALEUR #
#----------------------------------#
FREQUENCE_AFFICHAGE_INCREMENT 100000
FREQUENCE_AFFICHAGE_ITERATION 100000
resultats ---------------------------------------------- resultats ----------------------------------------------
COPIE 0 COPIE 0

View file

@ -1,7 +1,7 @@
#fichier au format maple6 #fichier au format maple6
############################################################################################### ###############################################################################################
# Visualisation elements finis : Herezh++ V6.687 # # Visualisation elements finis : Herezh++ V6.785 #
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # # Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
# http://www-lg2m.univ-ubs.fr # # 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 # >> 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 # 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 # 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 # puis pour chaque maillage
# >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , # >> 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 # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees
@ -47,6 +61,12 @@
#---------------------------------- torseurs de reactions -------------------------------- #---------------------------------- torseurs de reactions --------------------------------
#0 0 (nombre de torseurs et nombre total de grandeurs associees) #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 -------------------------------- #---------------------------------- grandeurs aux noeuds --------------------------------
#1 2 (nombre de noeuds, nombre total de grandeurs associees) #1 2 (nombre de noeuds, nombre total de grandeurs associees)
# noeud_8 [2]X [3]Y [4]Z [5]V2 [6]GAMMA2 ; # 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 # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
# etc. pour tous les points d'integration - element # 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

View file

@ -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"