2
0
Fork 0

ajout test : Batterie/FONCTION_UNIQUE/TYPE_DE_CALCUL/Test_R_dynamique_relaxation_dynam_visco (idem Test_R_dynamique_relaxation_dynam mais avec amortissement visqueux)

This commit is contained in:
troufflard 2017-04-06 10:32:06 +00:00
parent c89ecd8c74
commit 1e66763d77
10 changed files with 433 additions and 1335 deletions

View file

@ -1,7 +1,7 @@
------------------------------------------------------
Auteur
------------------------------------------------------
GŽrard Rio (gerard.rio@univ-ubs.fr)
Gérard Rio (gerard.rio@univ-ubs.fr)
------------------------------------------------------
@ -27,29 +27,29 @@ FONCTION_EXPRESSION_LITTERALE_nD
------------------------------------------------------
But du test
------------------------------------------------------
Test d'une loi additive pondŽrŽe par des grandeurs calculŽes localement par une loi membre
ici il s'agit de la pression d'hystŽrŽsis ˆ t, de la derni<6E>re pression d'hystŽrŽsis de rŽfŽrence,
et enfin de la pression prŽcŽdente de rŽfŽrence cf. le fichier elastomere_joint.loi
Test d'une loi additive pondérée par des grandeurs calculées localement par une loi membre
ici il s'agit de la pression d'hystérésis à t, de la dernière pression d'hystérésis de référence,
et enfin de la pression précédente de référence cf. le fichier elastomere_joint.loi
Dans ce test on se sert d'une fonction littŽrale nD ˆ 3 variables
Dans ce test on se sert d'une fonction littérale nD à 3 variables
------------------------------------------------------
Description du calcul
------------------------------------------------------
Le test est issu d'un exemple de calcul fourni par ƒmilie Vieville dans le cadre de sa th<74>se de doctorat.
Le test est issu d'un exemple de calcul fourni par Émilie Vieville dans le cadre de sa thèse de doctorat.
On consid<EFBFBD>re un tube creux de section rectangulaire soumis ˆ sur la face supŽrieure ˆ un dŽplacement de compression.
Le tube est modŽlisŽ par un seul ŽlŽment quadrangulaire axisymŽtrique. Les 3 autres faces sont bloquŽes perpendiculairement
ˆ la mati<74>re.
La loi de comportement est une combinaison additive d'hyperŽlasticitŽ et d'hystŽrŽsis sphŽrique et dŽviatorique.
On considère un tube creux de section rectangulaire soumis à sur la face supérieure à un déplacement de compression.
Le tube est modélisé par un seul élément quadrangulaire axisymétrique. Les 3 autres faces sont bloquées perpendiculairement
à la matière.
La loi de comportement est une combinaison additive d'hyperélasticité et d'hystérésis sphérique et déviatorique.
La description de la loi est complexe.
-------------------------------------------------------------
Grandeurs de comparaison
-------------------------------------------------------------
Principalement:
- rŽaction due au dŽplacement
- au premier point d'intŽgration: contrainte, dŽformation, globales et locales,
pondŽration et pressions
- réaction due au déplacement
- au premier point d'intégration: contrainte, déformation, globales et locales,
pondération et pressions

View file

@ -11,6 +11,7 @@ membrane
traction uniaxiale
ISOELAS2D_C
relaxation dynamique
amortissement cinetique
------------------------------------------------------
But du test

View file

@ -0,0 +1,37 @@
------------------------------------------------------
Auteur
------------------------------------------------------
Julien Troufflard (troufflard@univ-ubs.fr)
------------------------------------------------------
Mots-cles
------------------------------------------------------
coussin circulaire
membrane
traction uniaxiale
ISOELAS2D_C
relaxation dynamique
amortissement visqueux
------------------------------------------------------
But du test
------------------------------------------------------
test simple du mode de calcul : statique explicite (relaxation dynamique avec amortissement visqueux)
""
TYPE_DE_CALCUL
dynamique_relaxation_dynam
""
------------------------------------------------------
Description du calcul
------------------------------------------------------
simulation du gonflage d un coussin circulaire de rayon 100mm
- gonflage dans la direction Z
- modelisation d un quart de cercle par symetrie
- comportement de membrane elastique (TRIANGLE LINEAIRE + ISOELAS2D_C)
-------------------------------------------------------------
Grandeurs de comparaison
-------------------------------------------------------------
pour le noeud du set N_centre :
- coordonnee : X3

View file

@ -0,0 +1,116 @@
###############################################################################################
# Fichier de commande pour la visualisation elements finis #
# Herezh++ V6.687 #
# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
# http://www-lg2m.univ-ubs.fr #
###############################################################################################
debut_fichier_commande_visu # >>>>>> le mot cle: <debut_fichier_commande_visu>
# permet au programme de se positionner au debut du fichier, il est indispensable
# =================================================================================
# || ***** demande d'une visualisation maple: ***** ||
# =================================================================================
# un mot cle de debut (debut_visualisation_maple)
# un mot cle de fin ( fin_visualisation_maple)
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
debut_visualisation_maple
# ----------------------------- definition de la liste des increments a balayer: ----------------
debut_list_increment # un mot cle de debut de liste
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
# une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments)
# un mot cle de fin de liste ( fin_list_increment)
tous_les_increments_moins_zero dernier_increment fin_list_increment
# ----------------------------- choix des maillages a visualiser: ----------------
# la liste est facultative, par defaut la visualisation concerne le premier maillage
debut_choix_maillage # un mot cle de debut,
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
# une liste d'entiers , puis <fin_choix_maillage>, sur une meme ligne
1 fin_choix_maillage
# ----------------------------- definition des grandeurs a visualiser (maple): ----------------
debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple),
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
# les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin
debut_grandeurs_globales fin_grandeurs_globales
# ensuite pour chaque maillage:,
# le numero du maillage <un entier>,
# les infos pour la visualisation eventuelle des torseurs de reaction,
# les infos pour la visualisation eventuelle aux noeud,
# - tout d'abord les ddl principaux (position, deplacements, temperature etc.)
# - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses
# les infos pour la visualisation eventuelle aux elements,
# - tout d'abord les grandeurs generiques (deformations, contraintes etc.)
# - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def
# enfin un mot cle de fin ( fin_grandeurs_maple)
1 # le numero de maillage
# les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction
# une liste de nom <chaine de caracteres >, puis <fin_list_torseur_reaction>
debut_list_torseur_reaction fin_list_torseur_reaction
debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds
# debut de la liste de noeuds, puis une liste de numero de noeud <entier>, puis <fin_list_noeud>
deb_list_noeud fin_list_noeud
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
deb_list_ref_noeud N_centre fin_list_ref_noeud
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
deb_list_ddl_noeud X3 fin_list_ddl_noeud
type_sortie_ddl_retenue= 0
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
# debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis <fin_list_NbElement_NbPtInteg>
deb_list_NbElement_NbPtInteg fin_list_NbElement_NbPtInteg
# debut de la liste des ref d'elements, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_element>
deb_list_ref_element fin_list_ref_element
# debut de la liste des ref de ptinteg d'elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
deb_list_ddl_element fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
# informations particuliere dans le cas ou il y a une animation
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements)
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
# un parametre de pilotage du style de sortie
parametre_style_de_sortie 1
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
tenseur_en_absolue_ 0
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
# ----------------------------- definition des parametres d'animation: ----------------
debut_animation # un mot cle de debut de liste (debut_animation)
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
# des parametres avec des valeurs: (sur une meme ligne)
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l'animation)
fin_animation # un mot cle de fin
fin_visualisation_maple
# =================================================================================
# || fin de la visualisation maple ||
# =================================================================================
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
# l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus
# sans le mot cle de fin de fichier, le fichier n'est pas valide
###############################################################################################

View file

@ -0,0 +1,92 @@
dimension 3
niveau_commentaire 3
TYPE_DE_CALCUL ----------------------------------------------
dynamique_relaxation_dynam avec plus lectureCommandesVisu
PARA_TYPE_DE_CALCUL
typeCalRelaxation= 2 lambda= 1.1 type_calcul_mass= 2 option_recalcul_mass= 4
parametre_calcul_de_la_masse_ casMass_relax= 3
parametre_recalcul_de_la_masse_ fac_epsilon= 0.8
parametre_calcul_de_la_viscosite_ type_calcul_visqu_critique= 3 \
opt_cal_C_critique= 1 f_= 0.9
ARRET_A_EQUILIBRE_STATIQUE_ 1
< quart_cercle_membrane_tri3.her
< quart_cercle_membrane_tri3.lis
# renumerotation_des_noeuds_
les_courbes_1D ----------------------------------------------
rampe COURBEPOLYLINEAIRE_1_D
Debut_des_coordonnees_des_points
Coordonnee dim= 2 0. 0.
Coordonnee dim= 2 1. 1.
Fin_des_coordonnees_des_points
choix_materiaux ----------------------------------------------
E_tout MAT_ELAS
materiaux ----------------------------------------------
MAT_ELAS ISOELAS2D_C
400. 0.4
masse_volumique ----------------------------------------------
E_tout 1.e-9
epaisseurs ----------------------------------------------
E_tout 0.2
charges ----------------------------------------------
F_to PRESSION -0.01
blocages ----------------------------------------------
N_X0 UX
N_Y0 UY
N_cercle UZ
controle ----------------------------------------------
DELTAtMINI 1.e-7
DELTAtMAXI 0.1
TEMPSFIN 1.
DELTAt 0.1
ITERATIONS 999999999
PRECISION 1e-3
SAUVEGARDE DERNIER_CALCUL
MAXINCRE 99999999
NORME Residu/Reaction_et_VarRes
para_pilotage_equi_global ----------------------------------------------
para_syteme_lineaire ----------------------------------------------
para_affichage ----------------------------------------------
resultats ----------------------------------------------
COPIE 0
POINTS_INTEGRATION E_tout
Green-Lagrange Almansi Cauchy_global Def_mixte_local Sigma_mixte_local
_fin_point_info_

View file

@ -0,0 +1,97 @@
#fichier au format maple6
###############################################################################################
# Visualisation elements finis : Herezh++ V6.788 #
# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
# http://www-lg2m.univ-ubs.fr #
###############################################################################################
# entete des donnees : informations gererales: on trouve successivement:
# >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs
# precedes du numero de colonne entre crochet
# >> le nombre de maillages m, et dimension de l'espace de travail
# puis pour chaque maillage,
# >> 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
# et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere)
# precedes du numero de colonne correspondant entre crochet
# puis pour chaque maillage
# >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie ,
# les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant
# a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques,
# on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl
# chacun precedes du numero de colonne entre crochet
# puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres),
# suivi entre crochet, de la plage des numeros de colonnes, correspondant
# chacun sur une ligne differentes
# ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D:
# tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2)
# en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3)
# non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3)
# ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements,
# ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique
# ** ( cf: parametre_style_de_sortie = 0)
#====================================================================
#|| recapitulatif des differentes grandeurs par colonne ||
#====================================================================
#---------------------------------- grandeur globales ------------------------------------
#0 (nombre de grandeurs globales)
#---------------------------------- maillage et dimension --------------------------------
#1 3 (nombre de maillages et dimension)
#---------------------------------- 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 1 (nombre de noeuds, nombre total de grandeurs associees)
# reference N_centre : noeud_1 [2]X [3]Y [4]Z [5]X3 ;
#
#---------------------------------- grandeurs aux elements ------------------------------
#0 0 (nombre total d'elements, nombre totale de grandeurs associees)
#====================================================================
#|| fin du recapitulatif des differentes grandeurs ||
#====================================================================
# ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant
# a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun
# correspondant a un maillage. On trouve pour chaque enregistrement successivement :
# s'il y a des grandeurs globales: le temps puis les grandeurs globales,
# puis s'il y a des torseurs de reaction :
# de nouveau le temps, les composantes de la resultante puis les composantes du moments
# donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels
# puis s'il y a des grandeurs aux noeuds: de nouveau le temps
# les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud
# puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds
# puis s'il y a des grandeur aux elements:
# le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques)
# suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration
# correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes
# 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+00 0.000000000000e+00 0.000000000000e+00 2.648470271881e+01 2.648470271881e+01

View file

@ -0,0 +1,54 @@
nom_maillage quart_cercle_membrane_tri3
noeuds
20 NOEUDS
1 0 0 0
2 100 0 0
3 0 100 0
4 24.9999999998881 0 0
5 49.99999999982357 0 0
6 74.99999999986318 0 0
7 0 75.00000000017774 0
8 0 50.00000000025827 0
9 0 25.00000000017774 0
10 95.10565162954303 30.90169943740956 0
11 80.90169943759216 58.77852522911324 0
12 58.77852522936804 80.90169943740703 0
13 30.90169943755233 95.10565162949665 0
14 37.50000000010934 37.50000000021801 0
15 26.33854535607836 62.39459629332053 0
16 62.39459629305664 26.33854535595977 0
17 21.87500000008213 21.87500000011834 0
18 53.86242505428095 53.86242505420438 0
19 17.96875000003368 41.40625000017588 0
20 41.40625000013176 17.96875000018651 0
elements
25 ELEMENTS
1 TRIANGLE LINEAIRE 6 10 16
2 TRIANGLE LINEAIRE 7 15 13
3 TRIANGLE LINEAIRE 7 13 3
4 TRIANGLE LINEAIRE 6 2 10
5 TRIANGLE LINEAIRE 10 11 16
6 TRIANGLE LINEAIRE 12 13 15
7 TRIANGLE LINEAIRE 11 18 16
8 TRIANGLE LINEAIRE 12 15 18
9 TRIANGLE LINEAIRE 1 4 17
10 TRIANGLE LINEAIRE 1 17 9
11 TRIANGLE LINEAIRE 7 8 15
12 TRIANGLE LINEAIRE 5 6 16
13 TRIANGLE LINEAIRE 14 18 15
14 TRIANGLE LINEAIRE 14 16 18
15 TRIANGLE LINEAIRE 8 19 15
16 TRIANGLE LINEAIRE 5 16 20
17 TRIANGLE LINEAIRE 4 5 20
18 TRIANGLE LINEAIRE 8 9 19
19 TRIANGLE LINEAIRE 14 15 19
20 TRIANGLE LINEAIRE 14 20 16
21 TRIANGLE LINEAIRE 11 12 18
22 TRIANGLE LINEAIRE 4 20 17
23 TRIANGLE LINEAIRE 9 17 19
24 TRIANGLE LINEAIRE 14 19 17
25 TRIANGLE LINEAIRE 14 17 20

View file

@ -0,0 +1,16 @@
N_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
17 18 19 20
N_cercle 2 3 10 11 12 13
N_centre 1
N_X0 1 3 7 8 9
N_Y0 1 2 4 5 6
F_to 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 12 1 13 1 14 1 15 1 16 1
17 1 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1
E_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
17 18 19 20 21 22 23 24 25

File diff suppressed because it is too large Load diff

Binary file not shown.