2
0
Fork 0

mise a jour du test de Chung Lee cf. le readme

This commit is contained in:
rio 2017-01-10 10:08:00 +00:00
parent e969bcb4fb
commit c78c183151
3 changed files with 52 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, ref, figures dans le Readme
------------------------------------------------------ ------------------------------------------------------
Mots-cles Mots-cles
@ -22,12 +26,29 @@ test simple du mode de calcul : dynamique explicite Chung-Lee
PARA_TYPE_DE_CALCUL PARA_TYPE_DE_CALCUL
beta= 1.037 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 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
- 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 Grandeurs de comparaison
@ -39,3 +60,4 @@ pour le noeud 8 :
pour le point d integration 1 de l element 1 : pour le point d integration 1 de l element 1 :
- deformations : EPS11 EPS22 EPS33 - deformations : EPS11 EPS22 EPS33
- contrainte : SIG11 - contrainte : SIG11

View file

@ -1,7 +1,7 @@
dimension 3 dimension 3
niveau_commentaire 3 niveau_commentaire 0#3
TYPE_DE_CALCUL ---------------------------------------------- TYPE_DE_CALCUL ----------------------------------------------
@ -69,7 +69,12 @@ para_pilotage_equi_global ----------------------------------------------
para_syteme_lineaire ---------------------------------------------- para_syteme_lineaire ----------------------------------------------
para_affichage ---------------------------------------------- para_affichage #------------
#----------------------------------#
# PARAMETRE | VALEUR #
#----------------------------------#
FREQUENCE_AFFICHAGE_INCREMENT 100000
FREQUENCE_AFFICHAGE_ITERATION 100
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.713 # # 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 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