introduction des cas test de Bertille.
This commit is contained in:
rio
parent
0d06d7eb8f
commit
24a005eeb7
30 changed files with 3097 additions and 0 deletions
|
|
@ -0,0 +1,19 @@
|
|||
# ----- KDR
|
||||
typeCalRelaxation= 1 lambda= 1.6 type_calcul_mass= 2 option_recalcul_mass= 1
|
||||
parametre_calcul_de_la_masse_ casMass_relax= 3
|
||||
|
||||
# ----- Amortissement cinetique
|
||||
avec_amortissement_cinetique_
|
||||
max_nb_decroit_pourRelaxDyn_ 1
|
||||
coef_arret_pourRelaxDyn_ 0.
|
||||
coef_redemarrage_pourRelaxDyn_ 0.
|
||||
max_deltaX_pourRelaxDyn_ 0 #0.005
|
||||
nb_max_dX_OK_pourRelaxDyn_ 10 #1
|
||||
nb_deb_testfin_pourRelaxDyn_ 250 #1
|
||||
nb_deb_test_amort_cinetique_ 250 #1
|
||||
fi_parametre_amortissement_cinetique_
|
||||
|
||||
#mode_debug_= 5
|
||||
|
||||
ARRET_A_EQUILIBRE_STATIQUE_ 2
|
||||
|
||||
|
|
@ -0,0 +1,44 @@
|
|||
------------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------------
|
||||
Le but du calcul est de trouver la forme d’équilibre d’un tunnel gonflé
|
||||
dont le matériau suit une loi de comportement isotrope élastique avec prise
|
||||
en compte de la condition de contraintes planes (ISOELAS + LOI_CONTRAINTES_PLANES).
|
||||
|
||||
------------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------------
|
||||
Le calcul s’effectue avec la méthode de relaxation dynamique.
|
||||
Les paramètres de contrôle de cette méthode sont définis dans KDR.algo
|
||||
|
||||
La loi de comportement testée est la loi ISOELAS.
|
||||
Les paramètres utilisés pour la définir sont des paramètres cohérents
|
||||
avec des tests réalisés sur des membranes:
|
||||
E= 125 MPa
|
||||
nu = 0,4
|
||||
pression = 0,39e-3 MPa
|
||||
épaisseur = 25e-3 mm
|
||||
masse volumique = 0,9e-9 tonnes.mm^-3
|
||||
La condition de contraintes planes est ajoutée avec la loi LOI_CONTRAINTES_PLANES,
|
||||
utilisant la méthode de NEWTON_LOCAL.
|
||||
|
||||
La structure est une membrane de dimension 1000x2000mm,
|
||||
discrétisée en 5x10 élément membranes (interpolation linéaire - 8noeuds).
|
||||
|
||||
------------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------------
|
||||
tunnel
|
||||
isoelastique
|
||||
contraintes planes
|
||||
relaxation dynamique
|
||||
ISOELAS
|
||||
LOI_CONTRAINTES_PLANES
|
||||
|
||||
------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
------------------------------------------------------------
|
||||
Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII
|
||||
(point d’intégration 1 de l’élément 23, qui est situé environ au centre de la membrane)
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,198 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# STAMM #
|
||||
# (maillage automatique de pieces simples) #
|
||||
#############################################################################
|
||||
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
|
||||
# #
|
||||
#############################################################################
|
||||
# version 03.29 , compatible avec la version Herezh++6.576
|
||||
|
||||
|
||||
# ----------------------------------------------------------------
|
||||
# | construction automatique d'un maillage 2D |
|
||||
# | pour une utilisation par herezh++ |
|
||||
# ----------------------------------------------------------------
|
||||
# rectangle de dimension : 2000 x 1000 ( 10 x 5 elements )
|
||||
# geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire.
|
||||
|
||||
nom_maillage tunnel_intermediaire
|
||||
|
||||
noeuds ------------
|
||||
66 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 0
|
||||
2 0 200 0
|
||||
3 0 400 0
|
||||
4 0 600 0
|
||||
5 0 800 0
|
||||
6 0 1000 0
|
||||
7 200 0 0
|
||||
8 200 200 0
|
||||
9 200 400 0
|
||||
10 200 600 0
|
||||
11 200 800 0
|
||||
12 200 1000 0
|
||||
13 400 0 0
|
||||
14 400 200 0
|
||||
15 400 400 0
|
||||
16 400 600 0
|
||||
17 400 800 0
|
||||
18 400 1000 0
|
||||
19 600 0 0
|
||||
20 600 200 0
|
||||
21 600 400 0
|
||||
22 600 600 0
|
||||
23 600 800 0
|
||||
24 600 1000 0
|
||||
25 800 0 0
|
||||
26 800 200 0
|
||||
27 800 400 0
|
||||
28 800 600 0
|
||||
29 800 800 0
|
||||
30 800 1000 0
|
||||
31 1000 0 0
|
||||
32 1000 200 0
|
||||
33 1000 400 0
|
||||
34 1000 600 0
|
||||
35 1000 800 0
|
||||
36 1000 1000 0
|
||||
37 1200 0 0
|
||||
38 1200 200 0
|
||||
39 1200 400 0
|
||||
40 1200 600 0
|
||||
41 1200 800 0
|
||||
42 1200 1000 0
|
||||
43 1400 0 0
|
||||
44 1400 200 0
|
||||
45 1400 400 0
|
||||
46 1400 600 0
|
||||
47 1400 800 0
|
||||
48 1400 1000 0
|
||||
49 1600 0 0
|
||||
50 1600 200 0
|
||||
51 1600 400 0
|
||||
52 1600 600 0
|
||||
53 1600 800 0
|
||||
54 1600 1000 0
|
||||
55 1800 0 0
|
||||
56 1800 200 0
|
||||
57 1800 400 0
|
||||
58 1800 600 0
|
||||
59 1800 800 0
|
||||
60 1800 1000 0
|
||||
61 2000 0 0
|
||||
62 2000 200 0
|
||||
63 2000 400 0
|
||||
64 2000 600 0
|
||||
65 2000 800 0
|
||||
66 2000 1000 0
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_O 1 2 3 4 5 6
|
||||
N_S 1 7 13 19 25 31 37 43 49 55 61
|
||||
N_E 61 62 63 64 65 66
|
||||
N_N 6 12 18 24 30 36 42 48 54 60 66
|
||||
N_SO 1
|
||||
N_SE 61
|
||||
N_NE 66
|
||||
N_NO 6
|
||||
N_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 26 27 28 29 30 31 32
|
||||
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
|
||||
49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
|
||||
65 66
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
50 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 QUADRANGLE LINEAIRE 1 7 8 2
|
||||
2 QUADRANGLE LINEAIRE 2 8 9 3
|
||||
3 QUADRANGLE LINEAIRE 3 9 10 4
|
||||
4 QUADRANGLE LINEAIRE 4 10 11 5
|
||||
5 QUADRANGLE LINEAIRE 5 11 12 6
|
||||
6 QUADRANGLE LINEAIRE 7 13 14 8
|
||||
7 QUADRANGLE LINEAIRE 8 14 15 9
|
||||
8 QUADRANGLE LINEAIRE 9 15 16 10
|
||||
9 QUADRANGLE LINEAIRE 10 16 17 11
|
||||
10 QUADRANGLE LINEAIRE 11 17 18 12
|
||||
11 QUADRANGLE LINEAIRE 13 19 20 14
|
||||
12 QUADRANGLE LINEAIRE 14 20 21 15
|
||||
13 QUADRANGLE LINEAIRE 15 21 22 16
|
||||
14 QUADRANGLE LINEAIRE 16 22 23 17
|
||||
15 QUADRANGLE LINEAIRE 17 23 24 18
|
||||
16 QUADRANGLE LINEAIRE 19 25 26 20
|
||||
17 QUADRANGLE LINEAIRE 20 26 27 21
|
||||
18 QUADRANGLE LINEAIRE 21 27 28 22
|
||||
19 QUADRANGLE LINEAIRE 22 28 29 23
|
||||
20 QUADRANGLE LINEAIRE 23 29 30 24
|
||||
21 QUADRANGLE LINEAIRE 25 31 32 26
|
||||
22 QUADRANGLE LINEAIRE 26 32 33 27
|
||||
23 QUADRANGLE LINEAIRE 27 33 34 28
|
||||
24 QUADRANGLE LINEAIRE 28 34 35 29
|
||||
25 QUADRANGLE LINEAIRE 29 35 36 30
|
||||
26 QUADRANGLE LINEAIRE 31 37 38 32
|
||||
27 QUADRANGLE LINEAIRE 32 38 39 33
|
||||
28 QUADRANGLE LINEAIRE 33 39 40 34
|
||||
29 QUADRANGLE LINEAIRE 34 40 41 35
|
||||
30 QUADRANGLE LINEAIRE 35 41 42 36
|
||||
31 QUADRANGLE LINEAIRE 37 43 44 38
|
||||
32 QUADRANGLE LINEAIRE 38 44 45 39
|
||||
33 QUADRANGLE LINEAIRE 39 45 46 40
|
||||
34 QUADRANGLE LINEAIRE 40 46 47 41
|
||||
35 QUADRANGLE LINEAIRE 41 47 48 42
|
||||
36 QUADRANGLE LINEAIRE 43 49 50 44
|
||||
37 QUADRANGLE LINEAIRE 44 50 51 45
|
||||
38 QUADRANGLE LINEAIRE 45 51 52 46
|
||||
39 QUADRANGLE LINEAIRE 46 52 53 47
|
||||
40 QUADRANGLE LINEAIRE 47 53 54 48
|
||||
41 QUADRANGLE LINEAIRE 49 55 56 50
|
||||
42 QUADRANGLE LINEAIRE 50 56 57 51
|
||||
43 QUADRANGLE LINEAIRE 51 57 58 52
|
||||
44 QUADRANGLE LINEAIRE 52 58 59 53
|
||||
45 QUADRANGLE LINEAIRE 53 59 60 54
|
||||
46 QUADRANGLE LINEAIRE 55 61 62 56
|
||||
47 QUADRANGLE LINEAIRE 56 62 63 57
|
||||
48 QUADRANGLE LINEAIRE 57 63 64 58
|
||||
49 QUADRANGLE LINEAIRE 58 64 65 59
|
||||
50 QUADRANGLE LINEAIRE 59 65 66 60
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
A_O 1 4 2 4 3 4 4 4 5 4
|
||||
A_S 1 1 6 1 11 1 16 1 21 1 26 1 31 1 36 1 41 1 46 1
|
||||
A_E 46 2 47 2 48 2 49 2 50 2
|
||||
A_N 5 3 10 3 15 3 20 3 25 3 30 3 35 3 40 3 45 3 50 3
|
||||
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 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1
|
||||
34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 1
|
||||
45 1 46 1 47 1 48 1 49 1 50 1
|
||||
E_O 1 2 3 4 5
|
||||
E_S 1 6 11 16 21 26 31 36 41 46
|
||||
E_E 46 47 48 49 50
|
||||
E_N 5 10 15 20 25 30 35 40 45 50
|
||||
E_SO 1
|
||||
E_SE 46
|
||||
E_NE 50
|
||||
E_NO 5
|
||||
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 26 27 28 29 30 31 32
|
||||
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
|
||||
49 50
|
||||
|
||||
|
|
@ -0,0 +1,114 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.703 #
|
||||
# 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 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 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 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 23 1 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 Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII 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
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,139 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2015 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
||||
# #
|
||||
#######################################################################
|
||||
# (version la plus rapide )
|
||||
version 6.695
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- commentaire sur l'objet du travail ???
|
||||
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 4
|
||||
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
#non_dynamique
|
||||
|
||||
dynamique_relaxation_dynam #avec plus visualisation
|
||||
|
||||
|
||||
|
||||
PARA_TYPE_DE_CALCUL
|
||||
|
||||
|
||||
< KDR.algo
|
||||
|
||||
|
||||
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
# -- def maillage
|
||||
< tunnel_intermediaire.her
|
||||
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_to elas_iso
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
elas_iso LOI_CONTRAINTES_PLANES
|
||||
# ....... loi de comportement LoiContraintesPlanes ........
|
||||
NEWTON_LOCAL # paramètres par défaut
|
||||
ISOELAS
|
||||
# ....... loi de comportement isoelastique 3D ........
|
||||
# module d'young : coefficient de poisson
|
||||
125 0.4
|
||||
fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane
|
||||
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
epaisseurs #--------#
|
||||
E_to 25e-3
|
||||
masse_volumique #--------#
|
||||
E_to 0.9e-9
|
||||
|
||||
charges #------------#
|
||||
F_to PRESSION -0.390e-3
|
||||
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
N_O UZ UX UY
|
||||
N_E UZ UY
|
||||
|
||||
|
||||
controle #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
SAUVEGARDE 0
|
||||
ITERATIONS 300000
|
||||
PRECISION 0.005
|
||||
MULTIPLICATEUR 1
|
||||
DELTAt 0.1
|
||||
TEMPSFIN 1
|
||||
RESTART 0
|
||||
DELTAtMAXI 0.1
|
||||
NORME E_cinetique/E_statique_ET_ResSurReact
|
||||
|
||||
|
||||
para_affichage #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
FREQUENCE_AFFICHAGE_ITERATION 100
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
|
@ -0,0 +1,87 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.703 #
|
||||
# Copyright (c) 1997-2015, 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
|
||||
# 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)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_23 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| 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-01 4.882353771716e+02 4.424495084397e+02 6.181246032206e+02 7.970860513281e-03 -3.211642487148e-03 9.949699737870e-01 -3.467794108425e-03
|
||||
2.000000000000e-01 4.959653502342e+02 4.426378607101e+02 6.223508900465e+02 1.624443498337e-02 -6.521727255003e-03 2.029127668444e+00 -3.565388553768e-03
|
||||
3.000000000000e-01 4.972811916717e+02 4.428312561206e+02 6.291126103184e+02 2.456031467167e-02 -9.954539800997e-03 3.062275759675e+00 -1.940767538481e-02
|
||||
4.000000000000e-01 5.019141162088e+02 4.430175911776e+02 6.340647017512e+02 3.297754273946e-02 -1.329514671476e-02 4.115993687453e+00 -1.549644238806e-02
|
||||
5.000000000000e-01 5.035063352809e+02 4.432172415557e+02 6.398823579697e+02 4.153055631016e-02 -1.691100713672e-02 5.173533896005e+00 -4.446284032331e-02
|
||||
6.000000000000e-01 5.056522149909e+02 4.434135298018e+02 6.479514578033e+02 5.035724872586e-02 -2.050339030892e-02 6.273197256534e+00 -5.364551442781e-02
|
||||
7.000000000000e-01 5.154781672705e+02 4.436434523049e+02 6.515363839418e+02 6.045772732791e-02 -2.475911482453e-02 7.522927180923e+00 -8.571943983130e-02
|
||||
8.000000000000e-01 5.145020074068e+02 4.438367404669e+02 6.616539343554e+02 6.950550761708e-02 -2.837728697643e-02 8.653956103341e+00 -8.557912822277e-02
|
||||
9.000000000000e-01 5.271557119695e+02 4.440942687122e+02 6.678559043492e+02 8.054214347634e-02 -3.325641313039e-02 1.000588772921e+01 -1.546976821082e-01
|
||||
1.000000000000e+00 5.356044582301e+02 4.442871686731e+02 6.730004909504e+02 9.044645294194e-02 -3.695551097365e-02 1.125955889461e+01 -1.156164550013e-01
|
||||
|
|
@ -0,0 +1,20 @@
|
|||
# ----- KDR
|
||||
#typeCalRelaxation= 1 lambda= 1.6 type_calcul_mass= 2 option_recalcul_mass= 1
|
||||
typeCalRelaxation= 1 lambda= 2. type_calcul_mass= 2 option_recalcul_mass= 1
|
||||
parametre_calcul_de_la_masse_ casMass_relax= 3
|
||||
|
||||
# ----- Amortissement cinetique
|
||||
avec_amortissement_cinetique_
|
||||
max_nb_decroit_pourRelaxDyn_ 1
|
||||
coef_arret_pourRelaxDyn_ 0.
|
||||
coef_redemarrage_pourRelaxDyn_ 0.
|
||||
max_deltaX_pourRelaxDyn_ 0 #0.005
|
||||
nb_max_dX_OK_pourRelaxDyn_ 10 #1
|
||||
nb_deb_testfin_pourRelaxDyn_ 250 #1
|
||||
nb_deb_test_amort_cinetique_ 250 #1
|
||||
fi_parametre_amortissement_cinetique_
|
||||
|
||||
|
||||
|
||||
ARRET_A_EQUILIBRE_STATIQUE_ 2
|
||||
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
------------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------------
|
||||
Le but du calcul est de trouver la forme d’équilibre d’un tunnel gonflé
|
||||
dont le matériau suit une loi de comportement isotrope hyperélastique
|
||||
ISOHYOER3DORGEAS1 + LOI_CONTRAINTES_PLANES.
|
||||
|
||||
------------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------------
|
||||
Le calcul s’effectue avec la méthode de relaxation dynamique.
|
||||
Les paramètres de contrôle de la méthode de relaxation dynamique sont définis dans KDR.algo
|
||||
|
||||
La loi de comportement testée est:
|
||||
- ISOHYPER3DORGEAS1
|
||||
- LOI_CONTRAINTES_PLANES avec prise en compte de la condition de contraintes planes par la méthode de NEWTON_LOCAL
|
||||
Les paramètres utilisés pour définir ces lois sont des paramètres cohérents avec des tests réalisés sur des membranes:
|
||||
K = 624; Q0s = 2; mu01 = 30; mu02 = 14.6; mu03 = 90; alpha1 = 1e-3; alpha2 = 1e-3; Q0e =0.12
|
||||
|
||||
La structure est une membrane de dimension 1000x2000mm,
|
||||
discrétisée en 1x12 éléments membranes (interpolation linéaire - 8noeuds).
|
||||
|
||||
------------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------------
|
||||
tunnel
|
||||
hyperelastique
|
||||
contraintes planes
|
||||
relaxation dynamique
|
||||
isohyper3dorgeas1
|
||||
|
||||
------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
------------------------------------------------------------
|
||||
Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII
|
||||
(point d’intégration 1 de l’élément 6, qui est situé environ au centre de la membrane)
|
||||
|
|
@ -0,0 +1,115 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# STAMM #
|
||||
# (maillage automatique de pieces simples) #
|
||||
#############################################################################
|
||||
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
|
||||
# #
|
||||
#############################################################################
|
||||
# version 03.29 , compatible avec la version Herezh++6.576
|
||||
|
||||
|
||||
# ----------------------------------------------------------------
|
||||
# | construction automatique d'un maillage 2D |
|
||||
# | pour une utilisation par herezh++ |
|
||||
# ----------------------------------------------------------------
|
||||
# rectangle de dimension : 2000 x 1000 ( 12 x 1 elements )
|
||||
# geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire.
|
||||
|
||||
nom_maillage tunnel_12x1
|
||||
|
||||
noeuds ------------
|
||||
26 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 0
|
||||
2 0 1000 0
|
||||
3 166.6666666666667 0 0
|
||||
4 166.6666666666667 1000 0
|
||||
5 333.3333333333333 0 0
|
||||
6 333.3333333333333 1000 0
|
||||
7 500 0 0
|
||||
8 500 1000 0
|
||||
9 666.6666666666666 0 0
|
||||
10 666.6666666666666 1000 0
|
||||
11 833.3333333333333 0 0
|
||||
12 833.3333333333333 1000 0
|
||||
13 999.9999999999999 0 0
|
||||
14 999.9999999999999 1000 0
|
||||
15 1166.666666666667 0 0
|
||||
16 1166.666666666667 1000 0
|
||||
17 1333.333333333333 0 0
|
||||
18 1333.333333333333 1000 0
|
||||
19 1500 0 0
|
||||
20 1500 1000 0
|
||||
21 1666.666666666667 0 0
|
||||
22 1666.666666666667 1000 0
|
||||
23 1833.333333333333 0 0
|
||||
24 1833.333333333333 1000 0
|
||||
25 2000 0 0
|
||||
26 2000 1000 0
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_O 1 2
|
||||
N_S 1 3 5 7 9 11 13 15 17 19 21 23 25
|
||||
N_E 25 26
|
||||
N_N 2 4 6 8 10 12 14 16 18 20 22 24 26
|
||||
N_SO 1
|
||||
N_SE 25
|
||||
N_NE 26
|
||||
N_NO 2
|
||||
N_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 26
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
12 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 QUADRANGLE LINEAIRE 1 3 4 2
|
||||
2 QUADRANGLE LINEAIRE 3 5 6 4
|
||||
3 QUADRANGLE LINEAIRE 5 7 8 6
|
||||
4 QUADRANGLE LINEAIRE 7 9 10 8
|
||||
5 QUADRANGLE LINEAIRE 9 11 12 10
|
||||
6 QUADRANGLE LINEAIRE 11 13 14 12
|
||||
7 QUADRANGLE LINEAIRE 13 15 16 14
|
||||
8 QUADRANGLE LINEAIRE 15 17 18 16
|
||||
9 QUADRANGLE LINEAIRE 17 19 20 18
|
||||
10 QUADRANGLE LINEAIRE 19 21 22 20
|
||||
11 QUADRANGLE LINEAIRE 21 23 24 22
|
||||
12 QUADRANGLE LINEAIRE 23 25 26 24
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
A_O 1 4
|
||||
A_S 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1
|
||||
12 1
|
||||
A_E 12 2
|
||||
A_N 1 3 2 3 3 3 4 3 5 3 6 3 7 3 8 3 9 3 10 3 11 3
|
||||
12 3
|
||||
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
|
||||
E_O 1
|
||||
E_S 1 2 3 4 5 6 7 8 9 10 11
|
||||
12
|
||||
E_E 12
|
||||
E_N 1 2 3 4 5 6 7 8 9 10 11
|
||||
12
|
||||
E_SO 1
|
||||
E_SE 12
|
||||
E_NE 12
|
||||
E_NO 1
|
||||
E_to 1 2 3 4 5 6 7 8 9 10 11 12
|
||||
|
||||
|
|
@ -0,0 +1,113 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.703 #
|
||||
# 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 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 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 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 6 1 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 Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII 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
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,145 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2015 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
||||
# #
|
||||
#######################################################################
|
||||
# (version la plus rapide )
|
||||
version 6.695
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- commentaire sur l'objet du travail ???
|
||||
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 4
|
||||
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
#non_dynamique
|
||||
|
||||
dynamique_relaxation_dynam #avec plus visualisation
|
||||
|
||||
PARA_TYPE_DE_CALCUL
|
||||
|
||||
< KDR.algo
|
||||
#mode_debug_= 5
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
# -- def maillage
|
||||
< tunnel_12x1.her
|
||||
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_to orgeas2D
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
orgeas2D LOI_CONTRAINTES_PLANES
|
||||
|
||||
# ....... loi de comportement LoiContraintesPlanes ........
|
||||
NEWTON_LOCAL
|
||||
ISOHYPER3DORGEAS1
|
||||
|
||||
# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e |
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# parametres coherents avec des membranes, ajustes sur resultats experimentaux
|
||||
624 2 30 14.6 90 1e-3 1e-3 0.12
|
||||
|
||||
|
||||
fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane
|
||||
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
epaisseurs #--------#
|
||||
E_to 25e-3
|
||||
masse_volumique #--------#
|
||||
E_to 0.9e-9
|
||||
|
||||
charges #------------#
|
||||
F_to PRESSION -0.390e-3
|
||||
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
N_O UZ UX UY
|
||||
N_E UZ UY
|
||||
|
||||
controle #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
SAUVEGARDE 1#0
|
||||
ITERATIONS 300000
|
||||
PRECISION 0.005
|
||||
MULTIPLICATEUR 1
|
||||
DELTAt 0.1
|
||||
TEMPSFIN 1
|
||||
RESTART 0
|
||||
DELTAtMAXI 0.1
|
||||
NORME E_cinetique/E_statique_ET_ResSurReact
|
||||
MAXINCRE 1
|
||||
|
||||
para_pilotage_equi_global #------ parametre pour le jacobien
|
||||
|
||||
CAS_JACOBIEN_NEGATIF 2
|
||||
|
||||
para_affichage #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
FREQUENCE_AFFICHAGE_INCREMENT 1
|
||||
FREQUENCE_AFFICHAGE_ITERATION 100
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
|
@ -0,0 +1,78 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.710 #
|
||||
# Copyright (c) 1997-2015, 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
|
||||
# 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)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_6 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| 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-01 5.155709198288e+02 2.122580326756e+02 6.250076281113e+02 8.047703236407e-03 -3.248351140868e-03 1.001048740192e+00 2.258452053693e-03
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
------------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------------
|
||||
Ce programme teste la loi de comportement isotrope hyperélastique
|
||||
ISOHYPER3DORGEAS1 sur une menbrane carrée en traction simple
|
||||
en utilisant une condition de contrainte plane appliquée à la loi 3D.
|
||||
|
||||
------------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------------
|
||||
Le calcul est de type statique.
|
||||
|
||||
La loi de comportement testée est la loi ISOHYPER3DORGEAS1 (loi 3D).
|
||||
Les paramètres définissant cette loi sont des paramètres ajustés sur ses essais réalisés sur des membranes:
|
||||
K= 624; Q0s= 2; mu01=30; mu02= 14.6; mu03= 90; alpha01= 1e-3; alpha02= 1e-3; Q0e= 0.12
|
||||
Nous lui associons la loi LOI_CONTRAINTES_PLANES afin de faire le calcul en 2D prenant en compte
|
||||
l’hypothèse des contraintes planes (par la méthode de Newton).
|
||||
|
||||
La structure étudiée est une membrane carrée de dimensions 100*100 mm^2, discrétisée
|
||||
avec un seul élément de type membrane (interpolation linéaire).
|
||||
|
||||
L'essai simulé est un essai de traction simple dans la direction 11.
|
||||
|
||||
------------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------------
|
||||
quadrangle
|
||||
isohyper3dorgeas1
|
||||
contraintes planes
|
||||
|
||||
|
||||
------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
------------------------------------------------------------
|
||||
SIG11 SIG22 SIG12
|
||||
EPS11 EPS22 EPS12
|
||||
|
|
@ -0,0 +1,76 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# STAMM #
|
||||
# (maillage automatique de pieces simples) #
|
||||
#############################################################################
|
||||
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
|
||||
# #
|
||||
#############################################################################
|
||||
# version 03.29 , compatible avec la version Herezh++6.576
|
||||
|
||||
|
||||
# ----------------------------------------------------------------
|
||||
# | construction automatique d'un maillage 2D |
|
||||
# | pour une utilisation par herezh++ |
|
||||
# ----------------------------------------------------------------
|
||||
# rectangle de dimension : 100 x 100 ( 1 x 1 elements )
|
||||
# geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire.
|
||||
|
||||
nom_maillage quadrangle1
|
||||
|
||||
noeuds ------------
|
||||
4 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 0
|
||||
2 0 100 0
|
||||
3 100 0 0
|
||||
4 100 100 0
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_O 1 2
|
||||
N_S 1 3
|
||||
N_E 3 4
|
||||
N_N 2 4
|
||||
N_SO 1
|
||||
N_SE 3
|
||||
N_NE 4
|
||||
N_NO 2
|
||||
N_to 1 2 3 4
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
1 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 QUADRANGLE LINEAIRE 1 3 4 2
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
A_O 1 4
|
||||
A_S 1 1
|
||||
A_E 1 2
|
||||
A_N 1 3
|
||||
F_to 1 1
|
||||
E_O 1
|
||||
E_S 1
|
||||
E_E 1
|
||||
E_N 1
|
||||
E_SO 1
|
||||
E_SE 1
|
||||
E_NE 1
|
||||
E_NO 1
|
||||
E_to 1
|
||||
|
||||
|
|
@ -0,0 +1,112 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.699 #
|
||||
# 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 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 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 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 1 1 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 SIG11 SIG22 SIG12 EPS11 EPS22 EPS12 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
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,149 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2015 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
||||
# #
|
||||
#######################################################################
|
||||
# (version la plus rapide )
|
||||
version 6.697
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- commentaire sur l'objet du travail ???
|
||||
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 3
|
||||
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
non_dynamique #avec plus visualisation
|
||||
|
||||
|
||||
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
# -- def maillage
|
||||
< quadrangle1.her
|
||||
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_to orgeas2D
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
orgeas2D LOI_CONTRAINTES_PLANES
|
||||
|
||||
|
||||
NEWTON_LOCAL
|
||||
|
||||
ISOHYPER3DORGEAS1
|
||||
|
||||
# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e |
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# parametres coherents avec des membranes, ajustes sur resultats experimentaux
|
||||
624 2 30 14.6 90 1e-3 1e-3 0.12
|
||||
|
||||
fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane
|
||||
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
epaisseurs #--------#
|
||||
E_to 25e-3
|
||||
masse_volumique #--------#
|
||||
E_to 0.9e-9
|
||||
|
||||
charges #------------#
|
||||
# ici pas de chargement a definir
|
||||
|
||||
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
|
||||
# conditions limites pour un essai de traction selon l axe des X
|
||||
N_N UZ
|
||||
N_S UZ
|
||||
N_E UZ
|
||||
N_O UZ UX
|
||||
N_E 'UX=20'
|
||||
N_SO UY
|
||||
N_SE UY
|
||||
|
||||
|
||||
controle #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
SAUVEGARDE 0
|
||||
ITERATIONS 300000
|
||||
PRECISION 0.005
|
||||
MULTIPLICATEUR 1
|
||||
DELTAt 0.005
|
||||
TEMPSFIN 1
|
||||
RESTART 0
|
||||
DELTAtMAXI 0.1
|
||||
|
||||
|
||||
|
||||
para_affichage #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
|
||||
|
||||
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,113 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.703 #
|
||||
# Copyright (c) 1997-2015, 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
|
||||
# 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)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 6 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_1 pt_integ_1: [2]X [3]Y [4]Z [5] SIG11 [6] SIG22 [7] SIG12 [8] EPS11 [9] EPS22 [10] EPS12
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| 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
|
||||
|
||||
5.000000000000e-03 2.115361902706e+01 2.112406169668e+01 0.000000000000e+00 1.245967678291e-01 9.603156326874e-05 0.000000000000e+00 9.985019975028e-04 -3.989064216068e-04 0.000000000000e+00
|
||||
1.000000000000e-02 2.117475151360e+01 2.111563685285e+01 0.000000000000e+00 2.487309504414e-01 2.189700783719e-04 0.000000000000e+00 1.994015960096e-03 -7.982904127122e-04 0.000000000000e+00
|
||||
1.500000000000e-02 2.119588400014e+01 2.110721200902e+01 0.000000000000e+00 3.723318674740e-01 2.947116133150e-04 0.000000000000e+00 2.986553798226e-03 -1.198152735947e-03 0.000000000000e+00
|
||||
2.000000000000e-02 2.121701648668e+01 2.109878716518e+01 0.000000000000e+00 4.954025576742e-01 3.229168731565e-04 0.000000000000e+00 3.976127363058e-03 -1.598494155465e-03 0.000000000000e+00
|
||||
2.707106781187e-02 2.124690233575e+01 2.108687263677e+01 0.000000000000e+00 6.684817550666e-01 2.114892521729e-04 0.000000000000e+00 5.370558286617e-03 -2.165481776918e-03 0.000000000000e+00
|
||||
3.414213562373e-02 2.127678818483e+01 2.107495810836e+01 0.000000000000e+00 8.405845945444e-01 7.308843204531e-05 0.000000000000e+00 6.759117391982e-03 -2.733431293449e-03 0.000000000000e+00
|
||||
4.121320343560e-02 2.130667403390e+01 2.106304357995e+01 0.000000000000e+00 1.011647903613e+00 -1.642749678023e-04 0.000000000000e+00 8.141837600510e-03 -3.302344882098e-03 0.000000000000e+00
|
||||
4.828427124746e-02 2.133655988297e+01 2.105112905154e+01 0.000000000000e+00 1.181679194540e+00 -5.017109596818e-04 0.000000000000e+00 9.518751603152e-03 -3.872224726065e-03 0.000000000000e+00
|
||||
5.828427124746e-02 2.137882485605e+01 2.103427936387e+01 0.000000000000e+00 1.420256734858e+00 -1.286104655104e-03 0.000000000000e+00 1.145615326926e-02 -4.679810286139e-03 0.000000000000e+00
|
||||
6.828427124746e-02 2.142108982913e+01 2.101742967620e+01 0.000000000000e+00 1.656938955111e+00 -2.140067561486e-03 0.000000000000e+00 1.338209845641e-02 -5.489338950370e-03 0.000000000000e+00
|
||||
7.828427124746e-02 2.146335480221e+01 2.100057998853e+01 0.000000000000e+00 1.891602939365e+00 -3.203540448401e-03 0.000000000000e+00 1.529667731467e-02 -6.300816957409e-03 0.000000000000e+00
|
||||
8.828427124746e-02 2.150561977529e+01 2.098373030087e+01 0.000000000000e+00 2.124258568724e+00 -4.478167388988e-03 0.000000000000e+00 1.719997910912e-02 -7.114250570967e-03 0.000000000000e+00
|
||||
1.024264068712e-01 2.156539147344e+01 2.095990124404e+01 0.000000000000e+00 2.449586900641e+00 -6.901881440541e-03 0.000000000000e+00 1.987257471975e-02 -8.267970071699e-03 0.000000000000e+00
|
||||
1.165685424949e-01 2.162516317158e+01 2.093607218722e+01 0.000000000000e+00 2.771063439458e+00 -9.469614229263e-03 0.000000000000e+00 2.252303990196e-02 -9.425631241595e-03 0.000000000000e+00
|
||||
1.307106781187e-01 2.168493486973e+01 2.091224313040e+01 0.000000000000e+00 3.087874889780e+00 -1.236182894675e-02 0.000000000000e+00 2.515161831774e-02 -1.058725205668e-02 0.000000000000e+00
|
||||
1.448528137424e-01 2.174470656787e+01 2.088841407358e+01 0.000000000000e+00 3.394533227239e+00 -1.454248524462e-02 0.000000000000e+00 2.775855028480e-02 -1.175285059557e-02 0.000000000000e+00
|
||||
1.648528137424e-01 2.182923651403e+01 2.085471469824e+01 0.000000000000e+00 3.635970900778e+00 8.900885125220e-03 0.000000000000e+00 3.140881587810e-02 -1.340808160172e-02 0.000000000000e+00
|
||||
1.848528137424e-01 2.191376646020e+01 2.081559778866e+01 0.000000000000e+00 3.775275349391e+00 -8.523720758650e-05 0.000000000000e+00 3.501692149978e-02 -1.533949935420e-02 0.000000000000e+00
|
||||
2.048528137424e-01 2.199829640636e+01 2.077648087908e+01 0.000000000000e+00 3.910439062755e+00 -9.243746091742e-03 0.000000000000e+00 3.858351391342e-02 -1.728183650063e-02 0.000000000000e+00
|
||||
2.248528137424e-01 2.208282635252e+01 2.073736396950e+01 0.000000000000e+00 4.042992206342e+00 -1.891739832852e-02 0.000000000000e+00 4.210922752774e-02 -1.923517550773e-02 0.000000000000e+00
|
||||
2.531370849898e-01 2.220236974881e+01 2.068765464959e+01 0.000000000000e+00 4.241640540628e+00 -1.196950925605e-04 0.000000000000e+00 4.702676176270e-02 -2.173346133597e-02 0.000000000000e+00
|
||||
2.814213562373e-01 2.232191314510e+01 2.063794532968e+01 0.000000000000e+00 4.435204617174e+00 1.547688925023e-02 0.000000000000e+00 5.186550124481e-02 -2.424982128783e-02 0.000000000000e+00
|
||||
3.097056274848e-01 2.244145654139e+01 2.058332849849e+01 0.000000000000e+00 4.611002260797e+00 1.945234469155e-04 0.000000000000e+00 5.662712042139e-02 -2.703565370811e-02 0.000000000000e+00
|
||||
3.379898987322e-01 2.256099993768e+01 2.052871166730e+01 0.000000000000e+00 4.782487246584e+00 -1.663624311373e-02 0.000000000000e+00 6.131324949586e-02 -2.984375089417e-02 0.000000000000e+00
|
||||
3.779898987322e-01 2.273005983000e+01 2.045898895472e+01 0.000000000000e+00 5.039077829464e+00 -6.024431578780e-06 0.000000000000e+00 6.781464266963e-02 -3.346124066146e-02 0.000000000000e+00
|
||||
4.179898987322e-01 2.289911972233e+01 2.038926624214e+01 0.000000000000e+00 5.285582564498e+00 1.055045129288e-02 0.000000000000e+00 7.417257236372e-02 -3.711590472384e-02 0.000000000000e+00
|
||||
4.579898987322e-01 2.306817961465e+01 2.031954352955e+01 0.000000000000e+00 5.524282424941e+00 1.692089993712e-02 0.000000000000e+00 8.039122873540e-02 -4.080825418225e-02 0.000000000000e+00
|
||||
4.979898987322e-01 2.323723950697e+01 2.024982081697e+01 0.000000000000e+00 5.756098513475e+00 1.931845591194e-02 0.000000000000e+00 8.647465007341e-02 -4.453880895168e-02 0.000000000000e+00
|
||||
5.545584412272e-01 2.347632629955e+01 2.015121801123e+01 0.000000000000e+00 6.090949936850e+00 1.429757509892e-02 0.000000000000e+00 9.485458158242e-02 -4.988086006574e-02 0.000000000000e+00
|
||||
6.111269837221e-01 2.371541309213e+01 2.008095327997e+01 0.000000000000e+00 7.892359737027e+00 -1.364063141236e-03 0.000000000000e+00 1.029823461463e-01 -5.373573961986e-02 0.000000000000e+00
|
||||
6.676955262170e-01 2.395449988471e+01 2.001957237612e+01 0.000000000000e+00 9.920994044886e+00 -1.486750209981e-04 0.000000000000e+00 1.108679607376e-01 -5.713650214774e-02 0.000000000000e+00
|
||||
7.242640687119e-01 2.419358667729e+01 1.995819147227e+01 0.000000000000e+00 1.187929967617e+01 6.835735725519e-04 0.000000000000e+00 1.185209498365e-01 -6.056868977309e-02 0.000000000000e+00
|
||||
8.042640687119e-01 2.453170646193e+01 1.987138576558e+01 0.000000000000e+00 1.453402508573e+01 1.752436737719e-04 0.000000000000e+00 1.289643090895e-01 -6.547693781182e-02 0.000000000000e+00
|
||||
8.842640687119e-01 2.486982624658e+01 1.978458005889e+01 0.000000000000e+00 1.706330747541e+01 -4.133812178877e-03 0.000000000000e+00 1.389846198943e-01 -7.044993293332e-02 0.000000000000e+00
|
||||
9.642640687119e-01 2.520794603123e+01 1.969777435220e+01 0.000000000000e+00 1.947653199461e+01 -1.444581530112e-02 0.000000000000e+00 1.486044270979e-01 -7.548881897902e-02 0.000000000000e+00
|
||||
1.000000000000e+00 2.535898384862e+01 1.965899831758e+01 0.000000000000e+00 2.051931548721e+01 -2.156859637703e-02 0.000000000000e+00 1.527777777778e-01 -7.776128287454e-02 0.000000000000e+00
|
||||
BIN
Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/.DS_Store
vendored
Normal file
BIN
Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/.DS_Store
vendored
Normal file
Binary file not shown.
|
|
@ -0,0 +1,19 @@
|
|||
# ----- KDR
|
||||
typeCalRelaxation= 1 lambda= 1.6 type_calcul_mass= 2 option_recalcul_mass= 1
|
||||
parametre_calcul_de_la_masse_ casMass_relax= 3
|
||||
|
||||
# ----- Amortissement cinetique
|
||||
avec_amortissement_cinetique_
|
||||
max_nb_decroit_pourRelaxDyn_ 1
|
||||
coef_arret_pourRelaxDyn_ 0.
|
||||
coef_redemarrage_pourRelaxDyn_ 0.
|
||||
max_deltaX_pourRelaxDyn_ 0 #0.005
|
||||
nb_max_dX_OK_pourRelaxDyn_ 10 #1
|
||||
nb_deb_testfin_pourRelaxDyn_ 250 #1
|
||||
nb_deb_test_amort_cinetique_ 250 #1
|
||||
fi_parametre_amortissement_cinetique_
|
||||
|
||||
#mode_debug_= 5
|
||||
|
||||
ARRET_A_EQUILIBRE_STATIQUE_ 2
|
||||
|
||||
|
|
@ -0,0 +1,38 @@
|
|||
------------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------------
|
||||
Le but du calcul est de trouver la forme d’équilibre d’un tunnel gonflé
|
||||
dont le matériau suit une loi de comportement isotrope élastique (ISOELAST2D_C).
|
||||
|
||||
------------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------------
|
||||
Le calcul s’effectue avec la méthode de relaxation dynamique.
|
||||
Les paramètres de contrôle de cette méthode sont définis dans KDR.algo
|
||||
|
||||
La loi de comportement testée est la loi ISOELAST2D_C.
|
||||
Les paramètres utilisés pour la définir sont des paramètres cohérents avec des tests réalisés sur des membranes:
|
||||
E= 125 MPa
|
||||
nu = 0,4
|
||||
pression = 0,39e-3 MPa
|
||||
épaisseur = 25e-3 mm
|
||||
masse volumique = 0,9e-9 tonnes.mm^-3
|
||||
|
||||
La structure est une membrane de dimension 1000x2000mm,
|
||||
discrétisée en 10x20 élément membranes (interpolation linéaire - 8noeuds).
|
||||
|
||||
------------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------------
|
||||
tunnel
|
||||
isoelastique
|
||||
relaxation dynamique
|
||||
ISOELAST2D_C
|
||||
|
||||
------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
------------------------------------------------------------
|
||||
Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII
|
||||
(point d’intégration 1 de l’élément 116, qui est situé environ au centre de la membrane)
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,552 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# STAMM #
|
||||
# (maillage automatique de pieces simples) #
|
||||
#############################################################################
|
||||
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
|
||||
# #
|
||||
#############################################################################
|
||||
# version 03.29 , compatible avec la version Herezh++6.576
|
||||
|
||||
|
||||
# ----------------------------------------------------------------
|
||||
# | construction automatique d'un maillage 2D |
|
||||
# | pour une utilisation par herezh++ |
|
||||
# ----------------------------------------------------------------
|
||||
# rectangle de dimension : 2000 x 1000 ( 20 x 10 elements )
|
||||
# geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire.
|
||||
|
||||
nom_maillage tunnel_bis
|
||||
|
||||
noeuds ------------
|
||||
231 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 0
|
||||
2 0 100 0
|
||||
3 0 200 0
|
||||
4 0 300 0
|
||||
5 0 400 0
|
||||
6 0 500 0
|
||||
7 0 600 0
|
||||
8 0 700 0
|
||||
9 0 800 0
|
||||
10 0 900 0
|
||||
11 0 1000 0
|
||||
12 100 0 0
|
||||
13 100 100 0
|
||||
14 100 200 0
|
||||
15 100 300 0
|
||||
16 100 400 0
|
||||
17 100 500 0
|
||||
18 100 600 0
|
||||
19 100 700 0
|
||||
20 100 800 0
|
||||
21 100 900 0
|
||||
22 100 1000 0
|
||||
23 200 0 0
|
||||
24 200 100 0
|
||||
25 200 200 0
|
||||
26 200 300 0
|
||||
27 200 400 0
|
||||
28 200 500 0
|
||||
29 200 600 0
|
||||
30 200 700 0
|
||||
31 200 800 0
|
||||
32 200 900 0
|
||||
33 200 1000 0
|
||||
34 300 0 0
|
||||
35 300 100 0
|
||||
36 300 200 0
|
||||
37 300 300 0
|
||||
38 300 400 0
|
||||
39 300 500 0
|
||||
40 300 600 0
|
||||
41 300 700 0
|
||||
42 300 800 0
|
||||
43 300 900 0
|
||||
44 300 1000 0
|
||||
45 400 0 0
|
||||
46 400 100 0
|
||||
47 400 200 0
|
||||
48 400 300 0
|
||||
49 400 400 0
|
||||
50 400 500 0
|
||||
51 400 600 0
|
||||
52 400 700 0
|
||||
53 400 800 0
|
||||
54 400 900 0
|
||||
55 400 1000 0
|
||||
56 500 0 0
|
||||
57 500 100 0
|
||||
58 500 200 0
|
||||
59 500 300 0
|
||||
60 500 400 0
|
||||
61 500 500 0
|
||||
62 500 600 0
|
||||
63 500 700 0
|
||||
64 500 800 0
|
||||
65 500 900 0
|
||||
66 500 1000 0
|
||||
67 600 0 0
|
||||
68 600 100 0
|
||||
69 600 200 0
|
||||
70 600 300 0
|
||||
71 600 400 0
|
||||
72 600 500 0
|
||||
73 600 600 0
|
||||
74 600 700 0
|
||||
75 600 800 0
|
||||
76 600 900 0
|
||||
77 600 1000 0
|
||||
78 700 0 0
|
||||
79 700 100 0
|
||||
80 700 200 0
|
||||
81 700 300 0
|
||||
82 700 400 0
|
||||
83 700 500 0
|
||||
84 700 600 0
|
||||
85 700 700 0
|
||||
86 700 800 0
|
||||
87 700 900 0
|
||||
88 700 1000 0
|
||||
89 800 0 0
|
||||
90 800 100 0
|
||||
91 800 200 0
|
||||
92 800 300 0
|
||||
93 800 400 0
|
||||
94 800 500 0
|
||||
95 800 600 0
|
||||
96 800 700 0
|
||||
97 800 800 0
|
||||
98 800 900 0
|
||||
99 800 1000 0
|
||||
100 900 0 0
|
||||
101 900 100 0
|
||||
102 900 200 0
|
||||
103 900 300 0
|
||||
104 900 400 0
|
||||
105 900 500 0
|
||||
106 900 600 0
|
||||
107 900 700 0
|
||||
108 900 800 0
|
||||
109 900 900 0
|
||||
110 900 1000 0
|
||||
111 1000 0 0
|
||||
112 1000 100 0
|
||||
113 1000 200 0
|
||||
114 1000 300 0
|
||||
115 1000 400 0
|
||||
116 1000 500 0
|
||||
117 1000 600 0
|
||||
118 1000 700 0
|
||||
119 1000 800 0
|
||||
120 1000 900 0
|
||||
121 1000 1000 0
|
||||
122 1100 0 0
|
||||
123 1100 100 0
|
||||
124 1100 200 0
|
||||
125 1100 300 0
|
||||
126 1100 400 0
|
||||
127 1100 500 0
|
||||
128 1100 600 0
|
||||
129 1100 700 0
|
||||
130 1100 800 0
|
||||
131 1100 900 0
|
||||
132 1100 1000 0
|
||||
133 1200 0 0
|
||||
134 1200 100 0
|
||||
135 1200 200 0
|
||||
136 1200 300 0
|
||||
137 1200 400 0
|
||||
138 1200 500 0
|
||||
139 1200 600 0
|
||||
140 1200 700 0
|
||||
141 1200 800 0
|
||||
142 1200 900 0
|
||||
143 1200 1000 0
|
||||
144 1300 0 0
|
||||
145 1300 100 0
|
||||
146 1300 200 0
|
||||
147 1300 300 0
|
||||
148 1300 400 0
|
||||
149 1300 500 0
|
||||
150 1300 600 0
|
||||
151 1300 700 0
|
||||
152 1300 800 0
|
||||
153 1300 900 0
|
||||
154 1300 1000 0
|
||||
155 1400 0 0
|
||||
156 1400 100 0
|
||||
157 1400 200 0
|
||||
158 1400 300 0
|
||||
159 1400 400 0
|
||||
160 1400 500 0
|
||||
161 1400 600 0
|
||||
162 1400 700 0
|
||||
163 1400 800 0
|
||||
164 1400 900 0
|
||||
165 1400 1000 0
|
||||
166 1500 0 0
|
||||
167 1500 100 0
|
||||
168 1500 200 0
|
||||
169 1500 300 0
|
||||
170 1500 400 0
|
||||
171 1500 500 0
|
||||
172 1500 600 0
|
||||
173 1500 700 0
|
||||
174 1500 800 0
|
||||
175 1500 900 0
|
||||
176 1500 1000 0
|
||||
177 1600 0 0
|
||||
178 1600 100 0
|
||||
179 1600 200 0
|
||||
180 1600 300 0
|
||||
181 1600 400 0
|
||||
182 1600 500 0
|
||||
183 1600 600 0
|
||||
184 1600 700 0
|
||||
185 1600 800 0
|
||||
186 1600 900 0
|
||||
187 1600 1000 0
|
||||
188 1700 0 0
|
||||
189 1700 100 0
|
||||
190 1700 200 0
|
||||
191 1700 300 0
|
||||
192 1700 400 0
|
||||
193 1700 500 0
|
||||
194 1700 600 0
|
||||
195 1700 700 0
|
||||
196 1700 800 0
|
||||
197 1700 900 0
|
||||
198 1700 1000 0
|
||||
199 1800 0 0
|
||||
200 1800 100 0
|
||||
201 1800 200 0
|
||||
202 1800 300 0
|
||||
203 1800 400 0
|
||||
204 1800 500 0
|
||||
205 1800 600 0
|
||||
206 1800 700 0
|
||||
207 1800 800 0
|
||||
208 1800 900 0
|
||||
209 1800 1000 0
|
||||
210 1900 0 0
|
||||
211 1900 100 0
|
||||
212 1900 200 0
|
||||
213 1900 300 0
|
||||
214 1900 400 0
|
||||
215 1900 500 0
|
||||
216 1900 600 0
|
||||
217 1900 700 0
|
||||
218 1900 800 0
|
||||
219 1900 900 0
|
||||
220 1900 1000 0
|
||||
221 2000 0 0
|
||||
222 2000 100 0
|
||||
223 2000 200 0
|
||||
224 2000 300 0
|
||||
225 2000 400 0
|
||||
226 2000 500 0
|
||||
227 2000 600 0
|
||||
228 2000 700 0
|
||||
229 2000 800 0
|
||||
230 2000 900 0
|
||||
231 2000 1000 0
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_O 1 2 3 4 5 6 7 8 9 10 11
|
||||
N_S 1 12 23 34 45 56 67 78 89 100 111 122 133 144 155 166
|
||||
177 188 199 210 221
|
||||
N_E 221 222 223 224 225 226 227 228 229 230 231
|
||||
N_N 11 22 33 44 55 66 77 88 99 110 121 132 143 154 165 176
|
||||
187 198 209 220 231
|
||||
N_SO 1
|
||||
N_SE 221
|
||||
N_NE 231
|
||||
N_NO 11
|
||||
N_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 26 27 28 29 30 31 32
|
||||
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
|
||||
49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
|
||||
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
|
||||
81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
|
||||
97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112
|
||||
113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
|
||||
129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
|
||||
145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
|
||||
161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176
|
||||
177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
|
||||
193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208
|
||||
209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224
|
||||
225 226 227 228 229 230 231
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
200 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 QUADRANGLE LINEAIRE 1 12 13 2
|
||||
2 QUADRANGLE LINEAIRE 2 13 14 3
|
||||
3 QUADRANGLE LINEAIRE 3 14 15 4
|
||||
4 QUADRANGLE LINEAIRE 4 15 16 5
|
||||
5 QUADRANGLE LINEAIRE 5 16 17 6
|
||||
6 QUADRANGLE LINEAIRE 6 17 18 7
|
||||
7 QUADRANGLE LINEAIRE 7 18 19 8
|
||||
8 QUADRANGLE LINEAIRE 8 19 20 9
|
||||
9 QUADRANGLE LINEAIRE 9 20 21 10
|
||||
10 QUADRANGLE LINEAIRE 10 21 22 11
|
||||
11 QUADRANGLE LINEAIRE 12 23 24 13
|
||||
12 QUADRANGLE LINEAIRE 13 24 25 14
|
||||
13 QUADRANGLE LINEAIRE 14 25 26 15
|
||||
14 QUADRANGLE LINEAIRE 15 26 27 16
|
||||
15 QUADRANGLE LINEAIRE 16 27 28 17
|
||||
16 QUADRANGLE LINEAIRE 17 28 29 18
|
||||
17 QUADRANGLE LINEAIRE 18 29 30 19
|
||||
18 QUADRANGLE LINEAIRE 19 30 31 20
|
||||
19 QUADRANGLE LINEAIRE 20 31 32 21
|
||||
20 QUADRANGLE LINEAIRE 21 32 33 22
|
||||
21 QUADRANGLE LINEAIRE 23 34 35 24
|
||||
22 QUADRANGLE LINEAIRE 24 35 36 25
|
||||
23 QUADRANGLE LINEAIRE 25 36 37 26
|
||||
24 QUADRANGLE LINEAIRE 26 37 38 27
|
||||
25 QUADRANGLE LINEAIRE 27 38 39 28
|
||||
26 QUADRANGLE LINEAIRE 28 39 40 29
|
||||
27 QUADRANGLE LINEAIRE 29 40 41 30
|
||||
28 QUADRANGLE LINEAIRE 30 41 42 31
|
||||
29 QUADRANGLE LINEAIRE 31 42 43 32
|
||||
30 QUADRANGLE LINEAIRE 32 43 44 33
|
||||
31 QUADRANGLE LINEAIRE 34 45 46 35
|
||||
32 QUADRANGLE LINEAIRE 35 46 47 36
|
||||
33 QUADRANGLE LINEAIRE 36 47 48 37
|
||||
34 QUADRANGLE LINEAIRE 37 48 49 38
|
||||
35 QUADRANGLE LINEAIRE 38 49 50 39
|
||||
36 QUADRANGLE LINEAIRE 39 50 51 40
|
||||
37 QUADRANGLE LINEAIRE 40 51 52 41
|
||||
38 QUADRANGLE LINEAIRE 41 52 53 42
|
||||
39 QUADRANGLE LINEAIRE 42 53 54 43
|
||||
40 QUADRANGLE LINEAIRE 43 54 55 44
|
||||
41 QUADRANGLE LINEAIRE 45 56 57 46
|
||||
42 QUADRANGLE LINEAIRE 46 57 58 47
|
||||
43 QUADRANGLE LINEAIRE 47 58 59 48
|
||||
44 QUADRANGLE LINEAIRE 48 59 60 49
|
||||
45 QUADRANGLE LINEAIRE 49 60 61 50
|
||||
46 QUADRANGLE LINEAIRE 50 61 62 51
|
||||
47 QUADRANGLE LINEAIRE 51 62 63 52
|
||||
48 QUADRANGLE LINEAIRE 52 63 64 53
|
||||
49 QUADRANGLE LINEAIRE 53 64 65 54
|
||||
50 QUADRANGLE LINEAIRE 54 65 66 55
|
||||
51 QUADRANGLE LINEAIRE 56 67 68 57
|
||||
52 QUADRANGLE LINEAIRE 57 68 69 58
|
||||
53 QUADRANGLE LINEAIRE 58 69 70 59
|
||||
54 QUADRANGLE LINEAIRE 59 70 71 60
|
||||
55 QUADRANGLE LINEAIRE 60 71 72 61
|
||||
56 QUADRANGLE LINEAIRE 61 72 73 62
|
||||
57 QUADRANGLE LINEAIRE 62 73 74 63
|
||||
58 QUADRANGLE LINEAIRE 63 74 75 64
|
||||
59 QUADRANGLE LINEAIRE 64 75 76 65
|
||||
60 QUADRANGLE LINEAIRE 65 76 77 66
|
||||
61 QUADRANGLE LINEAIRE 67 78 79 68
|
||||
62 QUADRANGLE LINEAIRE 68 79 80 69
|
||||
63 QUADRANGLE LINEAIRE 69 80 81 70
|
||||
64 QUADRANGLE LINEAIRE 70 81 82 71
|
||||
65 QUADRANGLE LINEAIRE 71 82 83 72
|
||||
66 QUADRANGLE LINEAIRE 72 83 84 73
|
||||
67 QUADRANGLE LINEAIRE 73 84 85 74
|
||||
68 QUADRANGLE LINEAIRE 74 85 86 75
|
||||
69 QUADRANGLE LINEAIRE 75 86 87 76
|
||||
70 QUADRANGLE LINEAIRE 76 87 88 77
|
||||
71 QUADRANGLE LINEAIRE 78 89 90 79
|
||||
72 QUADRANGLE LINEAIRE 79 90 91 80
|
||||
73 QUADRANGLE LINEAIRE 80 91 92 81
|
||||
74 QUADRANGLE LINEAIRE 81 92 93 82
|
||||
75 QUADRANGLE LINEAIRE 82 93 94 83
|
||||
76 QUADRANGLE LINEAIRE 83 94 95 84
|
||||
77 QUADRANGLE LINEAIRE 84 95 96 85
|
||||
78 QUADRANGLE LINEAIRE 85 96 97 86
|
||||
79 QUADRANGLE LINEAIRE 86 97 98 87
|
||||
80 QUADRANGLE LINEAIRE 87 98 99 88
|
||||
81 QUADRANGLE LINEAIRE 89 100 101 90
|
||||
82 QUADRANGLE LINEAIRE 90 101 102 91
|
||||
83 QUADRANGLE LINEAIRE 91 102 103 92
|
||||
84 QUADRANGLE LINEAIRE 92 103 104 93
|
||||
85 QUADRANGLE LINEAIRE 93 104 105 94
|
||||
86 QUADRANGLE LINEAIRE 94 105 106 95
|
||||
87 QUADRANGLE LINEAIRE 95 106 107 96
|
||||
88 QUADRANGLE LINEAIRE 96 107 108 97
|
||||
89 QUADRANGLE LINEAIRE 97 108 109 98
|
||||
90 QUADRANGLE LINEAIRE 98 109 110 99
|
||||
91 QUADRANGLE LINEAIRE 100 111 112 101
|
||||
92 QUADRANGLE LINEAIRE 101 112 113 102
|
||||
93 QUADRANGLE LINEAIRE 102 113 114 103
|
||||
94 QUADRANGLE LINEAIRE 103 114 115 104
|
||||
95 QUADRANGLE LINEAIRE 104 115 116 105
|
||||
96 QUADRANGLE LINEAIRE 105 116 117 106
|
||||
97 QUADRANGLE LINEAIRE 106 117 118 107
|
||||
98 QUADRANGLE LINEAIRE 107 118 119 108
|
||||
99 QUADRANGLE LINEAIRE 108 119 120 109
|
||||
100 QUADRANGLE LINEAIRE 109 120 121 110
|
||||
101 QUADRANGLE LINEAIRE 111 122 123 112
|
||||
102 QUADRANGLE LINEAIRE 112 123 124 113
|
||||
103 QUADRANGLE LINEAIRE 113 124 125 114
|
||||
104 QUADRANGLE LINEAIRE 114 125 126 115
|
||||
105 QUADRANGLE LINEAIRE 115 126 127 116
|
||||
106 QUADRANGLE LINEAIRE 116 127 128 117
|
||||
107 QUADRANGLE LINEAIRE 117 128 129 118
|
||||
108 QUADRANGLE LINEAIRE 118 129 130 119
|
||||
109 QUADRANGLE LINEAIRE 119 130 131 120
|
||||
110 QUADRANGLE LINEAIRE 120 131 132 121
|
||||
111 QUADRANGLE LINEAIRE 122 133 134 123
|
||||
112 QUADRANGLE LINEAIRE 123 134 135 124
|
||||
113 QUADRANGLE LINEAIRE 124 135 136 125
|
||||
114 QUADRANGLE LINEAIRE 125 136 137 126
|
||||
115 QUADRANGLE LINEAIRE 126 137 138 127
|
||||
116 QUADRANGLE LINEAIRE 127 138 139 128
|
||||
117 QUADRANGLE LINEAIRE 128 139 140 129
|
||||
118 QUADRANGLE LINEAIRE 129 140 141 130
|
||||
119 QUADRANGLE LINEAIRE 130 141 142 131
|
||||
120 QUADRANGLE LINEAIRE 131 142 143 132
|
||||
121 QUADRANGLE LINEAIRE 133 144 145 134
|
||||
122 QUADRANGLE LINEAIRE 134 145 146 135
|
||||
123 QUADRANGLE LINEAIRE 135 146 147 136
|
||||
124 QUADRANGLE LINEAIRE 136 147 148 137
|
||||
125 QUADRANGLE LINEAIRE 137 148 149 138
|
||||
126 QUADRANGLE LINEAIRE 138 149 150 139
|
||||
127 QUADRANGLE LINEAIRE 139 150 151 140
|
||||
128 QUADRANGLE LINEAIRE 140 151 152 141
|
||||
129 QUADRANGLE LINEAIRE 141 152 153 142
|
||||
130 QUADRANGLE LINEAIRE 142 153 154 143
|
||||
131 QUADRANGLE LINEAIRE 144 155 156 145
|
||||
132 QUADRANGLE LINEAIRE 145 156 157 146
|
||||
133 QUADRANGLE LINEAIRE 146 157 158 147
|
||||
134 QUADRANGLE LINEAIRE 147 158 159 148
|
||||
135 QUADRANGLE LINEAIRE 148 159 160 149
|
||||
136 QUADRANGLE LINEAIRE 149 160 161 150
|
||||
137 QUADRANGLE LINEAIRE 150 161 162 151
|
||||
138 QUADRANGLE LINEAIRE 151 162 163 152
|
||||
139 QUADRANGLE LINEAIRE 152 163 164 153
|
||||
140 QUADRANGLE LINEAIRE 153 164 165 154
|
||||
141 QUADRANGLE LINEAIRE 155 166 167 156
|
||||
142 QUADRANGLE LINEAIRE 156 167 168 157
|
||||
143 QUADRANGLE LINEAIRE 157 168 169 158
|
||||
144 QUADRANGLE LINEAIRE 158 169 170 159
|
||||
145 QUADRANGLE LINEAIRE 159 170 171 160
|
||||
146 QUADRANGLE LINEAIRE 160 171 172 161
|
||||
147 QUADRANGLE LINEAIRE 161 172 173 162
|
||||
148 QUADRANGLE LINEAIRE 162 173 174 163
|
||||
149 QUADRANGLE LINEAIRE 163 174 175 164
|
||||
150 QUADRANGLE LINEAIRE 164 175 176 165
|
||||
151 QUADRANGLE LINEAIRE 166 177 178 167
|
||||
152 QUADRANGLE LINEAIRE 167 178 179 168
|
||||
153 QUADRANGLE LINEAIRE 168 179 180 169
|
||||
154 QUADRANGLE LINEAIRE 169 180 181 170
|
||||
155 QUADRANGLE LINEAIRE 170 181 182 171
|
||||
156 QUADRANGLE LINEAIRE 171 182 183 172
|
||||
157 QUADRANGLE LINEAIRE 172 183 184 173
|
||||
158 QUADRANGLE LINEAIRE 173 184 185 174
|
||||
159 QUADRANGLE LINEAIRE 174 185 186 175
|
||||
160 QUADRANGLE LINEAIRE 175 186 187 176
|
||||
161 QUADRANGLE LINEAIRE 177 188 189 178
|
||||
162 QUADRANGLE LINEAIRE 178 189 190 179
|
||||
163 QUADRANGLE LINEAIRE 179 190 191 180
|
||||
164 QUADRANGLE LINEAIRE 180 191 192 181
|
||||
165 QUADRANGLE LINEAIRE 181 192 193 182
|
||||
166 QUADRANGLE LINEAIRE 182 193 194 183
|
||||
167 QUADRANGLE LINEAIRE 183 194 195 184
|
||||
168 QUADRANGLE LINEAIRE 184 195 196 185
|
||||
169 QUADRANGLE LINEAIRE 185 196 197 186
|
||||
170 QUADRANGLE LINEAIRE 186 197 198 187
|
||||
171 QUADRANGLE LINEAIRE 188 199 200 189
|
||||
172 QUADRANGLE LINEAIRE 189 200 201 190
|
||||
173 QUADRANGLE LINEAIRE 190 201 202 191
|
||||
174 QUADRANGLE LINEAIRE 191 202 203 192
|
||||
175 QUADRANGLE LINEAIRE 192 203 204 193
|
||||
176 QUADRANGLE LINEAIRE 193 204 205 194
|
||||
177 QUADRANGLE LINEAIRE 194 205 206 195
|
||||
178 QUADRANGLE LINEAIRE 195 206 207 196
|
||||
179 QUADRANGLE LINEAIRE 196 207 208 197
|
||||
180 QUADRANGLE LINEAIRE 197 208 209 198
|
||||
181 QUADRANGLE LINEAIRE 199 210 211 200
|
||||
182 QUADRANGLE LINEAIRE 200 211 212 201
|
||||
183 QUADRANGLE LINEAIRE 201 212 213 202
|
||||
184 QUADRANGLE LINEAIRE 202 213 214 203
|
||||
185 QUADRANGLE LINEAIRE 203 214 215 204
|
||||
186 QUADRANGLE LINEAIRE 204 215 216 205
|
||||
187 QUADRANGLE LINEAIRE 205 216 217 206
|
||||
188 QUADRANGLE LINEAIRE 206 217 218 207
|
||||
189 QUADRANGLE LINEAIRE 207 218 219 208
|
||||
190 QUADRANGLE LINEAIRE 208 219 220 209
|
||||
191 QUADRANGLE LINEAIRE 210 221 222 211
|
||||
192 QUADRANGLE LINEAIRE 211 222 223 212
|
||||
193 QUADRANGLE LINEAIRE 212 223 224 213
|
||||
194 QUADRANGLE LINEAIRE 213 224 225 214
|
||||
195 QUADRANGLE LINEAIRE 214 225 226 215
|
||||
196 QUADRANGLE LINEAIRE 215 226 227 216
|
||||
197 QUADRANGLE LINEAIRE 216 227 228 217
|
||||
198 QUADRANGLE LINEAIRE 217 228 229 218
|
||||
199 QUADRANGLE LINEAIRE 218 229 230 219
|
||||
200 QUADRANGLE LINEAIRE 219 230 231 220
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
A_O 1 4 2 4 3 4 4 4 5 4 6 4 7 4 8 4 9 4 10 4
|
||||
A_S 1 1 11 1 21 1 31 1 41 1 51 1 61 1 71 1 81 1 91 1 101 1
|
||||
111 1 121 1 131 1 141 1 151 1 161 1 171 1 181 1 191 1
|
||||
A_E 191 2 192 2 193 2 194 2 195 2 196 2 197 2 198 2 199 2 200 2
|
||||
A_N 10 3 20 3 30 3 40 3 50 3 60 3 70 3 80 3 90 3 100 3 110 3
|
||||
120 3 130 3 140 3 150 3 160 3 170 3 180 3 190 3 200 3
|
||||
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 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1
|
||||
34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 1
|
||||
45 1 46 1 47 1 48 1 49 1 50 1 51 1 52 1 53 1 54 1 55 1
|
||||
56 1 57 1 58 1 59 1 60 1 61 1 62 1 63 1 64 1 65 1 66 1
|
||||
67 1 68 1 69 1 70 1 71 1 72 1 73 1 74 1 75 1 76 1 77 1
|
||||
78 1 79 1 80 1 81 1 82 1 83 1 84 1 85 1 86 1 87 1 88 1
|
||||
89 1 90 1 91 1 92 1 93 1 94 1 95 1 96 1 97 1 98 1 99 1
|
||||
100 1 101 1 102 1 103 1 104 1 105 1 106 1 107 1 108 1 109 1 110 1
|
||||
111 1 112 1 113 1 114 1 115 1 116 1 117 1 118 1 119 1 120 1 121 1
|
||||
122 1 123 1 124 1 125 1 126 1 127 1 128 1 129 1 130 1 131 1 132 1
|
||||
133 1 134 1 135 1 136 1 137 1 138 1 139 1 140 1 141 1 142 1 143 1
|
||||
144 1 145 1 146 1 147 1 148 1 149 1 150 1 151 1 152 1 153 1 154 1
|
||||
155 1 156 1 157 1 158 1 159 1 160 1 161 1 162 1 163 1 164 1 165 1
|
||||
166 1 167 1 168 1 169 1 170 1 171 1 172 1 173 1 174 1 175 1 176 1
|
||||
177 1 178 1 179 1 180 1 181 1 182 1 183 1 184 1 185 1 186 1 187 1
|
||||
188 1 189 1 190 1 191 1 192 1 193 1 194 1 195 1 196 1 197 1 198 1
|
||||
199 1 200 1
|
||||
E_O 1 2 3 4 5 6 7 8 9 10
|
||||
E_S 1 11 21 31 41 51 61 71 81 91 101
|
||||
111 121 131 141 151 161 171 181 191
|
||||
E_E 191 192 193 194 195 196 197 198 199 200
|
||||
E_N 10 20 30 40 50 60 70 80 90 100 110
|
||||
120 130 140 150 160 170 180 190 200
|
||||
E_SO 1
|
||||
E_SE 191
|
||||
E_NE 200
|
||||
E_NO 10
|
||||
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 26 27 28 29 30 31 32
|
||||
33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
|
||||
49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
|
||||
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
|
||||
81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
|
||||
97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112
|
||||
113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
|
||||
129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144
|
||||
145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
|
||||
161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176
|
||||
177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
|
||||
193 194 195 196 197 198 199 200
|
||||
|
||||
|
|
@ -0,0 +1,115 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.703 #
|
||||
# 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 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 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 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 116 1 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 Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII 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
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,136 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2015 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
||||
# #
|
||||
#######################################################################
|
||||
# (version la plus rapide )
|
||||
version 6.695
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- commentaire sur l'objet du travail ???
|
||||
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 4
|
||||
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
#non_dynamique
|
||||
|
||||
dynamique_relaxation_dynam
|
||||
|
||||
|
||||
PARA_TYPE_DE_CALCUL
|
||||
#TYPE 1
|
||||
|
||||
< KDR.algo
|
||||
|
||||
|
||||
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
# -- def maillage
|
||||
< tunnel_bis.her
|
||||
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_tout membrane
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
membrane ISOELAS2D_C
|
||||
# ....... loi de comportement isoelastique 3D ........
|
||||
# module d'young : coefficient de poisson
|
||||
125 0.4
|
||||
|
||||
# -- definition du type de deformation (par defaut: DEFORMATION_STANDART) --
|
||||
type_de_deformation DEFORMATION_STANDART
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
epaisseurs #--------#
|
||||
E_tout 25e-3
|
||||
masse_volumique #--------#
|
||||
E_tout 0.9e-9
|
||||
|
||||
charges #------------#
|
||||
F_to PRESSION -0.390e-3
|
||||
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
N_O UZ UX UY
|
||||
N_E UZ UY
|
||||
|
||||
|
||||
controle #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
SAUVEGARDE 0
|
||||
ITERATIONS 300000
|
||||
PRECISION 0.005
|
||||
MULTIPLICATEUR 1
|
||||
DELTAt 0.1
|
||||
TEMPSFIN 1
|
||||
DELTAtMAXI 0.1
|
||||
NORME E_cinetique/E_statique_ET_ResSurReact
|
||||
|
||||
|
||||
|
||||
|
||||
para_affichage #------------
|
||||
#FREQUENCE_AFFICHAGE_INCREMENT 1
|
||||
#FREQUENCE_AFFICHAGE_ITERATION 100
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,87 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.703 #
|
||||
# Copyright (c) 1997-2015, 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
|
||||
# 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)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_116 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| 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-01 7.730408920278e+02 5.210652559614e+02 6.245103853949e+02 8.033427076661e-03 -3.196644897709e-03 1.005173975830e+00 2.488978118371e-03
|
||||
2.000000000000e-01 7.681447945156e+02 5.209963438606e+02 6.345274366115e+02 1.612793207041e-02 -6.505133685726e-03 2.012779552994e+00 -8.029889518028e-03
|
||||
3.000000000000e-01 7.809829180392e+02 5.209274263835e+02 6.366382254352e+02 2.455888703561e-02 -9.846662840037e-03 3.068485401726e+00 -3.438694314249e-03
|
||||
4.000000000000e-01 7.845634757231e+02 5.208549836996e+02 6.418662485518e+02 3.316246027256e-02 -1.339468098880e-02 4.137587481702e+00 -1.930013091984e-02
|
||||
5.000000000000e-01 7.962474522385e+02 5.207802318866e+02 6.475992468084e+02 4.199020586806e-02 -1.709493636559e-02 5.230986803842e+00 -4.447232416248e-02
|
||||
6.000000000000e-01 8.016561095824e+02 5.207021944225e+02 6.546612524038e+02 5.127867716650e-02 -2.100063340888e-02 6.380717827819e+00 -7.279204498310e-02
|
||||
7.000000000000e-01 8.138162142176e+02 5.206263072773e+02 6.607739699857e+02 6.088097138966e-02 -2.484116374907e-02 7.581027662207e+00 -7.273440375093e-02
|
||||
8.000000000000e-01 8.212930601219e+02 5.205427028046e+02 6.707023601256e+02 7.136238854896e-02 -2.912196824722e-02 8.885952566976e+00 -8.586500411190e-02
|
||||
9.000000000000e-01 8.349537584645e+02 5.204682193168e+02 6.789366509884e+02 8.112411765326e-02 -3.298031591607e-02 1.010892727483e+01 -7.896857957788e-02
|
||||
1.000000000000e+00 8.429195604637e+02 5.203823488565e+02 6.882414598862e+02 9.263139024120e-02 -3.748051654705e-02 1.155344994381e+01 -6.368459085936e-02
|
||||
Binary file not shown.
|
|
@ -0,0 +1,31 @@
|
|||
------------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------------
|
||||
Ce programme teste la loi de comportement isotrope hyperélastique
|
||||
ISOHYPER3DORGEAS1 sur un cube en traction simple.
|
||||
|
||||
------------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------------
|
||||
Le calcul est de type statique.
|
||||
|
||||
La loi de comportement testée est la loi ISOHYPER3DORGEAS1 (loi 3D).
|
||||
Les paramètres définissant cette loi ont été ajustés sur des essais réalisés sur des membranes:
|
||||
K= 624; Q0s= 0.4; mu01= 44; mu02= 0.6; mu03= 44; alpha01= 1e-3; alpha02= 1e-3; Q0e= 0.01
|
||||
|
||||
La structure étudiée est un cube de dimensions 100*100*100 mm^3, discrétisé
|
||||
avec un seul élément de type volumique (interpolation linéaire, 8noeuds)
|
||||
|
||||
L'essai simulé est un essai de traction simple dans la direction 11.
|
||||
|
||||
------------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------------
|
||||
cube
|
||||
isohyper3dorgeas1
|
||||
|
||||
------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
------------------------------------------------------------
|
||||
SIG11 SIG22 SIG33 SIG12 SIG23 SIG13
|
||||
EPS11 EPS22 EPS33 EPS12 EPS23 EPS13
|
||||
|
|
@ -0,0 +1,133 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# STAMM #
|
||||
# (maillage automatique de pieces simples) #
|
||||
#############################################################################
|
||||
# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio #
|
||||
# #
|
||||
#############################################################################
|
||||
# version 03.29 , compatible avec la version Herezh++6.576
|
||||
|
||||
|
||||
# ----------------------------------------------------------------
|
||||
# | construction automatique d'un maillage 3D |
|
||||
# | pour une utilisation par herezh++ |
|
||||
# ----------------------------------------------------------------
|
||||
# prisme de dimension : 100 x 100 x 100 ( 1 x 1 x 1 elements )
|
||||
# geometrie prismatique, decoupage en hexaedre(s), interpolation lineaire.
|
||||
|
||||
nom_maillage cube1ele
|
||||
|
||||
noeuds ------------
|
||||
8 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 0
|
||||
2 0 0 100
|
||||
3 0 100 0
|
||||
4 0 100 100
|
||||
5 100 0 0
|
||||
6 100 0 100
|
||||
7 100 100 0
|
||||
8 100 100 100
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_avant 3 4 7 8
|
||||
N_arriere 1 2 5 6
|
||||
N_haut 2 4 6 8
|
||||
N_bas 1 3 5 7
|
||||
N_droit 1 2 3 4
|
||||
N_gauche 5 6 7 8
|
||||
N_haut_droit 2 4
|
||||
N_haut_gauche 6 8
|
||||
N_haut_arriere 2 6
|
||||
N_haut_avant 4 8
|
||||
N_bas_droit 1 3
|
||||
N_bas_gauche 5 7
|
||||
N_bas_arriere 1 5
|
||||
N_bas_avant 3 7
|
||||
N_arriere_droit 1 2
|
||||
N_arriere_gauche 5 6
|
||||
N_avant_droit 3 4
|
||||
N_avant_gauche 7 8
|
||||
N_haut_arriere_droit 2
|
||||
N_haut_arriere_gauche 6
|
||||
N_haut_avant_gauche 8
|
||||
N_haut_avant_droit 4
|
||||
N_bas_arriere_droit 1
|
||||
N_bas_arriere_gauche 5
|
||||
N_bas_avant_gauche 7
|
||||
N_bas_avant_droit 3
|
||||
N_to 1 2 3 4 5 6 7 8
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
1 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 HEXAEDRE LINEAIRE 5 7 3 1 6 8 4 2
|
||||
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
E_avant 1
|
||||
E_derriere 1
|
||||
E_haut 1
|
||||
E_bas 1
|
||||
E_droit 1
|
||||
E_gauche 1
|
||||
E_haut_droit 1
|
||||
E_haut_gauche 1
|
||||
E_haut_arriere 1
|
||||
E_haut_avant 1
|
||||
E_bas_droit 1
|
||||
E_bas_gauche 1
|
||||
E_bas_arriere 1
|
||||
E_bas_avant 1
|
||||
E_arriere_droit 1
|
||||
E_arriere_gauche 1
|
||||
E_avant_droit 1
|
||||
E_avant_gauche 1
|
||||
E_haut_arriere_droit 1
|
||||
E_haut_arriere_gauche 1
|
||||
E_haut_avant_gauche 1
|
||||
E_haut_avant_droit 1
|
||||
E_bas_arriere_droit 1
|
||||
E_bas_arriere_gauche 1
|
||||
E_bas_avant_gauche 1
|
||||
E_bas_avant_droit 1
|
||||
E_to 1
|
||||
|
||||
# references des faces et arretes d'elements
|
||||
|
||||
F_avant 1 5
|
||||
F_arriere 1 2
|
||||
F_haut 1 4
|
||||
F_bas 1 1
|
||||
F_droit 1 6
|
||||
F_gauche 1 3
|
||||
A_haut_droit 1 11
|
||||
A_haut_gauche 1 9
|
||||
A_haut_arriere 1 12
|
||||
A_haut_avant 1 10
|
||||
A_bas_droit 1 3
|
||||
A_bas_gauche 1 1
|
||||
A_bas_arriere 1 4
|
||||
A_bas_avant 1 2
|
||||
A_arriere_droit 1 8
|
||||
A_arriere_gauche 1 5
|
||||
A_avant_droit 1 7
|
||||
A_avant_gauche 1 6
|
||||
|
||||
|
|
@ -0,0 +1,114 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.699 #
|
||||
# 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)
|
||||
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 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 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 1 1 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 SIG11 SIG22 SIG33 SIG12 SIG23 SIG13 EPS11 EPS22 EPS33 EPS12 EPS23 EPS13 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
|
||||
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,165 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2015 Gerard Rio, gerard.rio@univ-ubs.fr #
|
||||
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
|
||||
# #
|
||||
#######################################################################
|
||||
# (version la plus rapide )
|
||||
version 6.697
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- commentaire sur l'objet du travail ???
|
||||
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 3
|
||||
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
non_dynamique #avec plus visualisation
|
||||
|
||||
|
||||
#-----------------------------------------------------------------------
|
||||
#| parametres (falcultatifs ) associes au calcul implicite statique |
|
||||
#-----------------------------------------------------------------------
|
||||
|
||||
## PARA_TYPE_DE_CALCUL
|
||||
# ................................................................
|
||||
# / acceleration de convergence par extrapolation methode MMPE /
|
||||
#..................................................................
|
||||
## acceleration_convergence_= 1 cas_acceleration_convergence_= 1 nb_vec_cycle_= 8
|
||||
#
|
||||
# acceleration_convergence_= indique si oui (=1) ou non (=0) on veut une acceleration
|
||||
# cas_acceleration_convergence_= 1 : on projette sur les vecteurs positions
|
||||
# 2 : on projette sur les vecteurs residus
|
||||
# nb_vec_cycle_= indique le nombre de vecteur maxi que l'on prendre en compte
|
||||
# pour la projection
|
||||
#..................................................................
|
||||
# / mode debug /
|
||||
#..................................................................
|
||||
#mode_debug_= 5
|
||||
#
|
||||
# mode_debug_= parametre de debug (=0) calcul classique
|
||||
# 1 : a chaque iteration on sort une visualisation
|
||||
# **important**: ce parametre doit etre mis sur une ligne differente de
|
||||
# l'acceleration de convergence
|
||||
#
|
||||
|
||||
|
||||
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
# -- def maillage
|
||||
< cube1ele.her
|
||||
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_to orgeas
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
orgeas ISOHYPER3DORGEAS1
|
||||
# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e |
|
||||
#-----------------------------------------------------------------------------------------
|
||||
# parametres coherents avec des membranes, ajustes sur resultats experimentaux
|
||||
624 2 30 14.6 90 1e-3 1e-3 0.12
|
||||
|
||||
|
||||
|
||||
# -- definition du type de deformation (par defaut: DEFORMATION_STANDART) --
|
||||
type_de_deformation DEFORMATION_STANDART
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
masse_volumique #--------#
|
||||
E_to 1
|
||||
|
||||
|
||||
charges
|
||||
# ici pas de chargement a definir
|
||||
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
|
||||
# conditions limites pour l essai de traction
|
||||
N_droit UX
|
||||
N_arriere UY
|
||||
N_bas UZ
|
||||
N_gauche 'UX=20'
|
||||
|
||||
|
||||
|
||||
controle #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
SAUVEGARDE 0
|
||||
ITERATIONS 300000
|
||||
PRECISION 0.005
|
||||
MULTIPLICATEUR 1
|
||||
DELTAt 0.005
|
||||
TEMPSFIN 1
|
||||
RESTART 0
|
||||
DELTAtMAXI 0.1
|
||||
|
||||
|
||||
|
||||
para_affichage #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,113 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.703 #
|
||||
# Copyright (c) 1997-2015, 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
|
||||
# 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)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 12 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_1 pt_integ_1: [2]X [3]Y [4]Z [5] SIG11 [6] SIG22 [7] SIG33 [8] SIG12 [9] SIG23 [10] SIG13 [11] EPS11 [12] EPS22 [13] EPS33 [14] EPS12 [15] EPS23 [16] EPS13
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| 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
|
||||
|
||||
5.000000000000e-03 2.115361902706e+01 7.883599917067e+01 7.883599917067e+01 1.246754248083e-01 9.260057771884e-05 9.260057771069e-05 2.698641420015e-34 3.179201731063e-15 -1.264538710493e-18 9.985019975030e-04 -3.998248006416e-04 -3.998248006417e-04 2.017792062434e-36 3.558397234536e-17 -1.419327565492e-20
|
||||
1.000000000000e-02 2.117475151360e+01 7.880448488185e+01 7.880448488185e+01 2.487203362453e-01 4.378397638141e-05 4.378397638141e-05 2.258352856448e-31 9.564229201673e-15 -2.525876536287e-18 1.994015960096e-03 -8.001293727149e-04 -8.001293727149e-04 2.531420612556e-33 1.068658098527e-16 -2.838090712908e-20
|
||||
1.500000000000e-02 2.119588400014e+01 7.877297059304e+01 7.877297059304e+01 3.721373352617e-01 -1.471408770181e-04 -1.471408770345e-04 -3.151671773921e-30 2.877267804918e-14 3.152852087291e-15 2.986553798226e-03 -1.200914484132e-03 -1.200914484132e-03 -1.261424708586e-32 3.209396208290e-16 3.546338609505e-17
|
||||
2.000000000000e-02 2.121701648668e+01 7.874145630423e+01 7.874145630423e+01 4.949289904455e-01 -4.808510016829e-04 -4.808510016829e-04 -4.052626957671e-30 3.847550531748e-14 3.147627423472e-15 3.976127363058e-03 -1.602180904342e-03 -1.602180904342e-03 -1.513412655257e-32 4.284334921657e-16 3.544224318952e-17
|
||||
2.707106781187e-02 2.124690233575e+01 7.869688836958e+01 7.869688836958e+01 6.675199999052e-01 -1.197819588030e-03 -1.197819588014e-03 -1.441168458111e-29 5.795042623729e-14 6.287355878473e-15 5.370558286617e-03 -2.170480461629e-03 -2.170480461629e-03 -7.071036449980e-32 6.437343313157e-16 7.090164628709e-17
|
||||
3.414213562373e-02 2.127678818483e+01 7.865232043493e+01 7.865232043493e+01 8.388724111112e-01 -2.203487543149e-03 -2.203487543132e-03 -2.922081750873e-30 8.081701077360e-14 1.562194203704e-15 6.759117391982e-03 -2.739746364703e-03 -2.739746364703e-03 -1.170820325417e-33 8.955862966472e-16 1.764296174247e-17
|
||||
4.121320343560e-02 2.130667403390e+01 7.860775250028e+01 7.860775250028e+01 1.008993068946e+00 -3.499578063698e-03 -3.499578063698e-03 -1.080627929555e-29 1.038700094798e-13 3.125685958814e-15 8.141837600510e-03 -3.309980805726e-03 -3.309980805726e-03 -4.032467466340e-32 1.148294254774e-15 3.535313377434e-17
|
||||
4.828427124746e-02 2.133655988297e+01 7.856318456563e+01 7.856318456563e+01 1.177888593891e+00 -5.087707111604e-03 -5.087707111604e-03 -2.928091900597e-29 1.205921932226e-13 6.248979357515e-15 9.518751603152e-03 -3.881185983081e-03 -3.881185983081e-03 -1.414695120625e-31 1.329977917946e-15 7.078411739725e-17
|
||||
5.828427124746e-02 2.137882485605e+01 7.850153129136e+01 7.850153129136e+01 1.420924356942e+00 1.778836690065e-07 1.778836690065e-07 3.122677735298e-15 1.520889211799e-17 -1.449068593345e-17 1.145615326926e-02 -4.672969880109e-03 -4.672969880109e-03 3.544598261520e-17 1.671791241742e-19 -1.644859461076e-19
|
||||
6.828427124746e-02 2.142108982913e+01 7.843987801710e+01 7.843987801710e+01 1.661553887406e+00 4.518916700280e-03 4.518916700280e-03 9.344644646674e-15 -1.252314128054e-13 -1.574330285611e-15 1.338209845641e-02 -5.466621522972e-03 -5.466621522972e-03 1.062950350269e-16 -1.372036130272e-15 -1.790795682236e-17
|
||||
7.828427124746e-02 2.146335480221e+01 7.837822474283e+01 7.837822474283e+01 1.899791408606e+00 8.466112807827e-03 8.466112807818e-03 1.242853128120e-14 -2.452641324545e-13 -1.927644176175e-17 1.529667731467e-02 -6.262146790748e-03 -6.262146790748e-03 1.416702459251e-16 -2.678327586566e-15 -2.197281547730e-19
|
||||
8.828427124746e-02 2.150561977529e+01 7.831450148208e+01 7.831450148208e+01 2.126190274600e+00 4.056138742312e-07 4.056138656686e-07 3.099636839317e-15 -6.644225028668e-15 1.528084042314e-15 1.719997910912e-02 -7.086356845695e-03 -7.086356845695e-03 3.540550292127e-17 -7.231573320331e-17 1.745449122873e-17
|
||||
1.024264068712e-01 2.156539147344e+01 7.822665044411e+01 7.822665044411e+01 2.452639937988e+00 4.893950551095e-07 4.893950551095e-07 4.538919521320e-31 -6.693345758016e-15 3.063706421145e-15 1.987257471975e-02 -8.225944875921e-03 -8.225944875921e-03 -6.704925374937e-49 -7.251909107094e-17 3.509917242802e-17
|
||||
1.165685424949e-01 2.162516317158e+01 7.813879940613e+01 7.813879940613e+01 2.774254243306e+00 -1.099583339563e-03 -1.099583339519e-03 -3.077833676010e-15 -3.355380061689e-15 3.049395342465e-15 2.252303990196e-02 -9.369378764813e-03 -9.369378764812e-03 -3.536830484083e-17 -3.619287119000e-17 3.504151146736e-17
|
||||
1.307106781187e-01 2.168493486973e+01 7.805094836815e+01 7.805094836815e+01 3.090634186812e+00 -3.070516967621e-03 -3.070516967568e-03 -9.198359692841e-15 -1.361332398737e-14 7.633554271603e-15 2.515161831774e-02 -1.051667583709e-02 -1.051667583709e-02 -1.060507807566e-16 -1.462360687732e-16 8.800964709843e-17
|
||||
1.448528137424e-01 2.174470656787e+01 7.796309733017e+01 7.796309733017e+01 3.397608088517e+00 -3.640918407576e-03 -3.640918407496e-03 -1.219232805421e-14 -1.368260537136e-14 9.109292047365e-15 2.775855028480e-02 -1.166785351515e-02 -1.166785351515e-02 -1.413303313424e-16 -1.466900077400e-16 1.055925708056e-16
|
||||
1.648528137424e-01 2.182923651403e+01 7.782695993902e+01 7.782695993902e+01 3.630891115443e+00 1.691978495235e-05 1.691978496919e-05 2.858117939270e-15 4.313592403843e-17 1.391349847677e-15 3.140881587810e-02 -1.345947045607e-02 -1.345947045607e-02 3.531900230236e-17 4.877875694900e-19 1.719351318512e-17
|
||||
1.848528137424e-01 2.191376646020e+01 7.768538164037e+01 7.768538164037e+01 3.774726117038e+00 3.036236455006e-06 3.036236455006e-06 1.952285839935e-31 -6.110243048725e-15 1.283961988856e-15 3.501692149978e-02 -1.533269185833e-02 -1.533269185833e-02 -1.908966389196e-35 -7.379037166445e-17 1.712628670278e-17
|
||||
2.048528137424e-01 2.199829640636e+01 7.754380334172e+01 7.754380334172e+01 3.911298424153e+00 -3.353151710913e-03 -3.353151710898e-03 -2.476766896937e-15 -1.165918454092e-14 2.435196911312e-15 3.858351391342e-02 -1.721618292187e-02 -1.721618292187e-02 -3.530399559257e-17 -1.488837228210e-16 3.471145432795e-17
|
||||
2.248528137424e-01 2.208282635252e+01 7.740222504307e+01 7.740222504307e+01 4.041777089805e+00 -1.066073173812e-02 -1.066073173811e-02 -7.009631604381e-15 -1.117294446732e-14 1.124861022343e-15 4.210922752774e-02 -1.911001885331e-02 -1.911001885331e-02 -1.058928929021e-16 -1.496347314162e-16 1.699301682192e-17
|
||||
2.531370849898e-01 2.220236974881e+01 7.720702231727e+01 7.720702231727e+01 4.242156394540e+00 2.577368087696e-06 2.577368087696e-06 -2.176771897216e-15 -2.617435738778e-15 -4.583015611096e-17 4.702676176271e-02 -2.173827173030e-02 -2.173827173030e-02 -3.528532883763e-17 -3.707173389460e-17 -7.429038068358e-19
|
||||
2.814213562373e-01 2.232191314510e+01 7.701181959148e+01 7.701181959148e+01 4.431240939974e+00 2.669684945765e-03 2.669684945752e-03 4.097350408477e-15 1.292539703571e-14 -3.121216608497e-15 5.186550124481e-02 -2.438653551191e-02 -2.438653551191e-02 7.054900914095e-17 1.916864592340e-16 -5.374173968333e-17
|
||||
3.097056274848e-01 2.244145654139e+01 7.681661686568e+01 7.681661686568e+01 4.609222836896e+00 -2.680564849878e-03 -2.680564849898e-03 7.774154339240e-15 1.753198745518e-14 -4.909386874531e-15 5.662712042139e-02 -2.705501385940e-02 -2.705501385940e-02 1.410605927022e-16 2.701453759769e-16 -8.907991687667e-17
|
||||
3.379898987322e-01 2.256099993768e+01 7.662141413988e+01 7.662141413988e+01 4.776266882282e+00 -1.605509988623e-02 -1.605509988626e-02 1.113349612981e-14 3.659356340441e-14 -3.764011254849e-15 6.131324949586e-02 -2.974391303156e-02 -2.974391303156e-02 2.115435002830e-16 5.822630148238e-16 -7.151860535730e-17
|
||||
3.779898987322e-01 2.273005983000e+01 7.635399097803e+01 7.635399097803e+01 5.039088489648e+00 7.786784323796e-06 7.786784310695e-06 1.753430675707e-15 -4.667414834065e-15 -2.686028178850e-15 6.781464266963e-02 -3.346117460201e-02 -3.346117460201e-02 3.524057961813e-17 -7.705718981576e-17 -5.398399332503e-17
|
||||
4.179898987322e-01 2.289911972233e+01 7.608656781617e+01 7.608656781617e+01 5.281101461860e+00 4.188695379502e-04 4.188695379874e-04 -8.349974158114e-15 -4.180499680814e-14 -2.563877898162e-15 7.417257236372e-02 -3.721770044001e-02 -3.721770044001e-02 -1.761336499809e-16 -7.108068172303e-16 -5.408222393956e-17
|
||||
4.579898987322e-01 2.306817961465e+01 7.581914465432e+01 7.581914465432e+01 5.502768903361e+00 -1.486716933725e-02 -1.486716933719e-02 -1.920581604346e-14 -7.341157130583e-14 -1.662346231887e-15 8.039122873540e-02 -4.101404548037e-02 -4.101404548037e-02 -4.225881034147e-16 -1.278036045639e-15 -3.657682338320e-17
|
||||
4.979898987322e-01 2.323723950697e+01 7.555978607164e+01 7.555978607164e+01 5.747205570340e+00 6.846022119505e-06 6.846022125279e-06 -4.625456062906e-15 9.551384603328e-17 -4.304665733468e-15 8.647465007341e-02 -4.473447565687e-02 -4.473447565687e-02 -1.055996277891e-16 1.694392636938e-18 -9.827595225819e-17
|
||||
5.545584412272e-01 2.347632629955e+01 7.519299764650e+01 7.519299764650e+01 6.072495284825e+00 -1.349997809418e-02 -1.349997809426e-02 1.477529737856e-14 1.046096116767e-13 -4.870811226558e-15 9.485458158242e-02 -5.006182316306e-02 -5.006182316306e-02 3.518213344976e-16 1.887643342859e-15 -1.159811042653e-16
|
||||
6.111269837221e-01 2.371541309213e+01 7.494407477561e+01 7.494407477561e+01 7.893443913008e+00 1.174474057125e-03 1.174474057110e-03 7.071438838960e-31 -2.672724509365e-15 5.209333833797e-15 1.029823461463e-01 -5.372189528702e-02 -5.372189528702e-02 4.982209748012e-33 -3.934430703275e-17 1.034334563857e-16
|
||||
6.676955262170e-01 2.395449988471e+01 7.471413002634e+01 7.471413002634e+01 9.921062208520e+00 3.190389927665e-05 3.190389928598e-05 2.062306067108e-15 1.790013034956e-16 4.069697219451e-16 1.108679607376e-01 -5.713547648234e-02 -5.713547648234e-02 3.492323853827e-17 2.200646573397e-18 6.891654397969e-18
|
||||
7.242640687119e-01 2.419358667729e+01 7.448418527706e+01 7.448418527706e+01 1.187711853859e+01 -3.449522271665e-03 -3.449522271601e-03 9.231584851297e-15 -7.648278494211e-15 -7.628931811219e-15 1.185209498365e-01 -6.058072128757e-02 -6.058072128757e-02 1.391677793116e-16 -8.200736945582e-17 -1.150075004227e-16
|
||||
8.042640687119e-01 2.453170646193e+01 7.415899429404e+01 7.415899429404e+01 1.452614570673e+01 -1.360483111091e-02 -1.360483111078e-02 1.811545479203e-14 -4.344610922162e-15 -2.279882166040e-14 1.289643090895e-01 -6.550785078148e-02 -6.550785078148e-02 2.422979411704e-16 -4.018174953881e-17 -3.049389382075e-16
|
||||
8.842640687119e-01 2.486982624658e+01 7.383380331102e+01 7.383380331102e+01 1.704568019323e+01 -3.271895550780e-02 -3.271895550757e-02 3.369765411616e-14 -1.540453391274e-14 -3.668502104018e-14 1.389846198943e-01 -7.050022616153e-02 -7.050022616153e-02 4.133377972658e-16 -1.271366858575e-16 -4.499810502276e-16
|
||||
9.642640687119e-01 2.520794603123e+01 7.350861232800e+01 7.350861232800e+01 1.944401396488e+01 -6.385033599079e-02 -6.385033599044e-02 4.475374837916e-14 -1.125837325714e-14 -4.985158812666e-14 1.486044270979e-01 -7.555900452951e-02 -7.555900452950e-02 5.142620258237e-16 -8.477323560279e-17 -5.728409268280e-16
|
||||
1.000000000000e+00 2.535898384862e+01 7.337486990659e+01 7.337486990659e+01 2.052550812407e+01 1.213215295226e-05 1.213215295226e-05 3.046647892073e-15 4.224030430256e-16 1.197441103707e-14 1.527777777778e-01 -7.765909162409e-02 -7.765909162409e-02 3.420428927494e-17 3.072528279066e-18 1.344350359864e-16
|
||||
Loading…
Reference in a new issue