2
0
Fork 0

mise à jour du test de Zhai (cf. readme), initialement l algorithme ne fonctionnait pas. La methode originale propose par Zhai (prediction - correction) a ete mise a jour, et les resultats sont surprenants et interessant

This commit is contained in:
rio 2017-01-12 20:02:21 +00:00
parent 3e6c561328
commit 949c4028b9
3 changed files with 58 additions and 11 deletions

View file

@ -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,8 +24,20 @@ 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
@ -29,6 +45,15 @@ 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
-------------------------------------------------------------

View file

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

View file

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