ajout d'un test de flexion utilisant la loi hypo-élastique orthotrope entraînée
This commit is contained in:
rio
parent
b867faa2ad
commit
1284ce2c09
7 changed files with 795 additions and 0 deletions
|
|
@ -0,0 +1,35 @@
|
|||
------------------------------------------------------
|
||||
Auteur
|
||||
------------------------------------------------------
|
||||
Gérard Rio (gerard.rio@univ-ubs.fr)
|
||||
|
||||
------------------------------------------------------
|
||||
Mots-cles
|
||||
------------------------------------------------------
|
||||
HYPO_ORTHO3D
|
||||
ORTHOELA3D
|
||||
non_dynamique
|
||||
|
||||
------------------------------------------------------
|
||||
But du test
|
||||
------------------------------------------------------
|
||||
test de fonctionnement de la loi HYPO_ORTHO3D avec une comparaison avec une loi ORTHOELA3D de mêmes coefficient.
|
||||
""
|
||||
TYPE_DE_CALCUL
|
||||
non_dynamique
|
||||
""
|
||||
|
||||
------------------------------------------------------
|
||||
Description du calcul
|
||||
------------------------------------------------------
|
||||
On considère une poutre de 100 x 10 x 15 mm ( 4 x 1 x 2 éléments ), géométrie prismatique, découpage en hexaèdres, interpolation quadratique complète.
|
||||
|
||||
La poutre est encastrée à une extrémité et soumise à un déplacement imposé à l'autre de 1 mm.
|
||||
|
||||
\figures:comparaison_hypo_ortho.pdf[0.8] \legende: Comparaison entre hypo-élastique orthotrope entraînée et ortho-élastique entraînée en implicite et en explicite (relaxation visqueuse avec la même précision que l'implicite) \fin_legende
|
||||
|
||||
|
||||
-------------------------------------------------------------
|
||||
Grandeurs de comparaison
|
||||
-------------------------------------------------------------
|
||||
On regarde en particulier le torseur de réaction à l'encastrement.
|
||||
|
|
@ -0,0 +1,9 @@
|
|||
plot 'flexion_poutre_hypo_implicit.maple' u ($1*10):(-$3) w lp t 'hypo ortho en implicit '
|
||||
replot 'flexion_poutre_ortho_elas_imp.maple' u ($1*10):(-$3) w lp t ' ortho elas en implicit '
|
||||
replot 'flexion_poutre_hypo_relaxation.maple' u ($1*10):(-$3) w lp t 'hypo ortho en relaxation (precision 1.e-3 comme implicite) '
|
||||
replot 'flexion_poutre_ortho_elas_relaxation.maple' u ($1*10):(-$3) w lp t ' ortho elas en relaxation (precision 1.e-3 comme implicite) '
|
||||
set xlabel 'deplacement en mm'
|
||||
replot
|
||||
set ylabel ' resultante opposee au deplacement impose en N '
|
||||
replot
|
||||
|
||||
Binary file not shown.
|
|
@ -0,0 +1,140 @@
|
|||
###############################################################################################
|
||||
# Fichier de commande pour la visualisation elements finis #
|
||||
# Herezh++ V6.912 #
|
||||
# Copyright (c) 1997-2019, 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 N_droit N_gauche fin_list_torseur_reaction
|
||||
# les sommes, moyennes etc. sur ref de noeuds: un mot cle de debut: debut_list_SM_sur_refN
|
||||
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refN>
|
||||
debut_list_SM_sur_refN fin_list__SM_sur_refN
|
||||
# debut de la liste des ddl a considerer <deb_list_ddl_SM_sur_refN>, (une liste de ddl), puis <fin_list_ddl_SM_sur_refN>
|
||||
deb_list_ddl_SM_sur_refN fin_list_ddl_SM_sur_refN
|
||||
# debut de la liste des ddl etendus a considerer, (une liste de ddl), puis <fin_list_ddl_etendu__SM_sur_refN>
|
||||
deb_list_ddl_etendu__SM_sur_refN fin_list_ddl_etendu__SM_sur_refN # fin des ddl etendu pour _SM_sur_refNs
|
||||
# debut de la liste des grandeurs particulieres a considerer, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refN>
|
||||
deb_list_GrandParticuliere__SM_sur_refN fin_list_GrandParticuliere__SM_sur_refN # fin des grandeurs particulieres aux _SM_sur_refNs
|
||||
# les SM sur ref d'elements: un mot cle de debut: debut_list_SM_sur_refE
|
||||
# une liste de nom de ref d'elements <chaine de caracteres >, le pti associe
|
||||
# puis <fin_list__SM_sur_refE>
|
||||
debut_list_SM_sur_refE fin_list__SM_sur_refE
|
||||
# les SMs sur ref de pti: un mot cle de debut: debut_list_SM_sur_refpti
|
||||
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refpti>
|
||||
debut_list_SM_sur_refpti fin_list__SM_sur_refpti
|
||||
# debut de la liste des ddl a considerer , (une liste de ddl), puis <fin_list_ddl__SM_sur_refEpti>
|
||||
deb_list_ddl__SM_sur_refEpti fin_list_ddl__SM_sur_refEpti # fin de la liste de ddl a considerer pour les elements
|
||||
# debut de la liste des grandeurs particulieres a considerer pour les _SM_sur_refEptis, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refEpti>
|
||||
deb_list_GrandParticuliere__SM_sur_refEpti fin_list_GrandParticuliere__SM_sur_refEpti # fin des grandeurs particulieres aux _SM_sur_refEptis
|
||||
# tableau de grandeurs evoluees aux _SM_sur_refEptis a visualiser, un par maillage
|
||||
deb_list_GrandEvoluee__SM_sur_refEpti fin_list_GrandEvoluee__SM_sur_refEpti
|
||||
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 EPS_TRANSPORTEE_ANISO SIGMA_DANS_ANISO PARA_ORTHO 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,225 @@
|
|||
|
||||
#######################################################################
|
||||
# #
|
||||
# | | ==== === ==== ==== | | | | #
|
||||
# | | | | | | / | | | | #
|
||||
# |====| |=== === |=== / |====| ------- ------- #
|
||||
# | | | | \ | / | | | | #
|
||||
# | | ==== | \ ==== ==== | | | | #
|
||||
# #
|
||||
#######################################################################
|
||||
# Copyright (c) 1997-2019 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 avec le plus de verifications pendant le calcul et les I/O )
|
||||
version 6.911
|
||||
|
||||
|
||||
#######################################################################
|
||||
# fichier de commande #
|
||||
|
||||
|
||||
# --- test de flexion de la loi hypo-elast orthotrope en 3D ----
|
||||
# le fichier contient egalement une loi ortho-elas equivalente
|
||||
# pour comparaison manuel eventuel
|
||||
|
||||
|
||||
#----------------------------------------------------
|
||||
# definition de la dimension de l'espace de travail |
|
||||
#----------------------------------------------------
|
||||
dimension 3
|
||||
|
||||
#---------------------------------------------------------------
|
||||
# definition facultative du niveau d'impression (entre 0 et 10)|
|
||||
#---------------------------------------------------------------
|
||||
niveau_commentaire 3#10#3
|
||||
|
||||
#-------------------------------
|
||||
# definition du type de calcul |
|
||||
#-------------------------------
|
||||
TYPE_DE_CALCUL
|
||||
|
||||
non_dynamique #avec plus visualisation
|
||||
|
||||
|
||||
#--------------------------------------
|
||||
#| definition du (ou des) maillage(s) |
|
||||
#--------------------------------------
|
||||
|
||||
< poutre.her
|
||||
|
||||
#----------------------------------------
|
||||
# definition facultative de fonction nD|
|
||||
#----------------------------------------
|
||||
|
||||
les_fonctions_nD #------------
|
||||
|
||||
fct_repere1 FONCTION_EXPRESSION_LITTERALE_nD
|
||||
# un_argument= i un_argument= j
|
||||
fct= 1. , 1. , 0. , -1. , 1. , 0.
|
||||
fin_parametres_fonction_expression_litterale_
|
||||
|
||||
# f(norme_de_convergence)
|
||||
fonct_lambda FONCTION_EXPRESSION_LITTERALE_nD
|
||||
un_argument= amor_cinet_visqueux
|
||||
fct= (amor_cinet_visqueux == 1) \ on test si cinétique ou pas
|
||||
? 2.5 \ si cinétique on prend la valeur 2.5
|
||||
: 5 # sinon 5 en visqueux
|
||||
fin_parametres_fonction_expression_litterale_
|
||||
|
||||
fct_essai FONCTION_EXPRESSION_LITTERALE_nD
|
||||
un_argument= compteur_iteration_algo_global
|
||||
fct= (compteur_iteration_algo_global < 2000000 ) ? 1 : 0 # juste pour voir si ça marche
|
||||
fin_parametres_fonction_expression_litterale_
|
||||
|
||||
|
||||
|
||||
#------------- fin def des fonctions nD ------------------------
|
||||
|
||||
#======================================
|
||||
#| definition des lois de comportement|
|
||||
#--------------------------------------
|
||||
|
||||
choix_materiaux ----------
|
||||
#----------------------------------------
|
||||
# Elements | Nom Materiau |
|
||||
#----------------------------------------
|
||||
E_tout acier_hypo_ortho
|
||||
#E_tout acier_elas_ortho
|
||||
#E_tout acier_ortho_essai
|
||||
#E_tout acier_iso
|
||||
|
||||
materiaux #----------
|
||||
|
||||
#------------------------------------------
|
||||
# Nom Materiau | Type loi |
|
||||
#------------------------------------------
|
||||
acier_hypo_ortho HYPO_ORTHO3D
|
||||
# ....... loi de comportement hypo elastique orthotrope 3D ........
|
||||
# 9 parametres materiau: 3 modules tangent de traction, 3 coef type Poisson
|
||||
# , 3 modules type cisaillement tangent
|
||||
# et un nom de repere associe
|
||||
E1= 100000 E2= 50000 E3= 20000 \
|
||||
nu12= 0.2 nu13= 0.3 nu23= 0.1\
|
||||
G12= 20000 G13= 10000 G23= 2000
|
||||
nom_repere_associe_ repere1 avec_parametres_de_reglage_
|
||||
sortie_post_ 1
|
||||
type_transport_ 0
|
||||
permet_affichage_ 3
|
||||
fin_parametres_reglage_
|
||||
|
||||
acier_elas_ortho ORTHOELA3D
|
||||
# ....... loi de comportement hypo elastique orthotrope 3D ........
|
||||
# 9 parametres materiau: 3 modules tangent de traction, 3 coef type Poisson
|
||||
# , 3 modules type cisaillement tangent
|
||||
# et un nom de repere associe
|
||||
E1= 100000 E2= 50000 E3= 20000 \
|
||||
nu12= 0.2 nu13= 0.3 nu23= 0.1\
|
||||
G12= 20000 G13= 10000 G23= 2000
|
||||
nom_repere_associe_ repere1 avec_parametres_de_reglage_
|
||||
sortie_post_ 1
|
||||
type_transport_ 0
|
||||
permet_affichage_ 3
|
||||
fin_parametres_reglage_
|
||||
|
||||
|
||||
acier_ortho_essai ORTHOELA3D
|
||||
# en fait on utilise une loi iso
|
||||
# ....... loi de comportement orthotrope elastique 3D ........
|
||||
# 9 parametres materiau:
|
||||
# 3 modules, 3 coef de Poisson, 3 modules de cisaillement
|
||||
E1= 100000 E2= 100000 E3= 100000 \
|
||||
nu12= 0.3 nu13= 0.3 nu23= 0.3\
|
||||
G12= 50000 G13= 50000 G23= 50000
|
||||
nom_repere_associe_ repere1 avec_parametres_de_reglage_
|
||||
sortie_post_ 1
|
||||
type_transport_ 0
|
||||
permet_affichage_ 3
|
||||
fin_parametres_reglage_
|
||||
|
||||
acier_iso ISOELAS
|
||||
100000 0.3
|
||||
#--------------------------------------- fin def des lois de comportement ------------------------
|
||||
|
||||
# --- divers stockages (1) -------
|
||||
masse_volumique #--------#
|
||||
E_tout 1
|
||||
|
||||
repere_anisotropie_ #--------#
|
||||
#----------------------------------------------------
|
||||
E_tout repere1 orthonorme_ par_fonction_nD_ une_fonction_nD_ fct_repere1
|
||||
|
||||
charges #------------#
|
||||
#N_gauche PONCTUELLE 50 0 0
|
||||
#N_gauche PONCTUELLE 0 50 0
|
||||
blocages #------------#
|
||||
#-----------------------------------------------------------
|
||||
# nom du maillage | Ref noeud | Bloquages
|
||||
#-----------------------------------------------------------
|
||||
N_droit UX UY UZ
|
||||
|
||||
N_gauche 'UY= -1.'
|
||||
|
||||
|
||||
|
||||
controle ------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
TEMPSFIN 1.
|
||||
SAUVEGARDE 0
|
||||
ITERATIONS 25
|
||||
PRECISION 1e-3
|
||||
|
||||
|
||||
|
||||
para_pilotage_equi_global ------------
|
||||
#-----------------------------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#-----------------------------------------------
|
||||
FACTEUR_DIMINUTION 1.732
|
||||
FACTEUR_AUGMENTATION 1.414
|
||||
NB_BONNE_CONVERGENCE 3
|
||||
#INIT_COMP_TANGENT_SIMPLE -1
|
||||
#SUR_SOUS_RELAXATION 1.2
|
||||
NORME_MAXI_INCREMENT 1#0.2
|
||||
|
||||
#para_syteme_lineaire ------------
|
||||
#-----------------------------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#-----------------------------------------------
|
||||
#TYPE_MATRICE CARREE #_LAPACK
|
||||
#TYPE_MATRICE BANDE_NON_SYMETRIQUE_LAPACK
|
||||
#SYMETRIE_MATRICE 0
|
||||
#TYPE_MATRICE BANDE_SYMETRIQUE_LAPACK
|
||||
#SYMETRIE_MATRICE 1
|
||||
|
||||
para_dedies_dynamique ------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
#BULK_VISCOSITY 1
|
||||
#AMORTISSEMENT_VISCOSITE_ARTIFICIELLE 1
|
||||
|
||||
|
||||
para_affichage #------------
|
||||
#---------------------------
|
||||
# PARAMETRE | VALEUR |
|
||||
#---------------------------
|
||||
FREQUENCE_AFFICHAGE_ITERATION 1#20
|
||||
FREQUENCE_AFFICHAGE_INCREMENT 1#000
|
||||
FREQUENCE_SORTIE_FIL_DU_CALCUL 1#000
|
||||
|
||||
|
||||
# ------------------------------------------------------------------------------------
|
||||
resultats pas_de_sortie_finale_
|
||||
COPIE 0
|
||||
#
|
||||
|
||||
|
||||
_fin_point_info_
|
||||
|
||||
|
||||
|
||||
|
|
@ -0,0 +1,105 @@
|
|||
#fichier au format maple6
|
||||
###############################################################################################
|
||||
# Visualisation elements finis : Herezh++ V6.915 #
|
||||
# Copyright (c) 1997-2019, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ #
|
||||
# http://www-lg2m.univ-ubs.fr #
|
||||
###############################################################################################
|
||||
|
||||
|
||||
|
||||
# entete des donnees : informations gererales: on trouve successivement:
|
||||
# >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs
|
||||
# precedes du numero de colonne entre crochet
|
||||
# >> le nombre de maillages m, et dimension de l'espace de travail
|
||||
# puis pour chaque maillage,
|
||||
# >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit
|
||||
# correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions
|
||||
# des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment
|
||||
# ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds
|
||||
# >> le nombre de ref de noeuds, nombre total de grandeurs associees
|
||||
# puis pour chaque maillage
|
||||
# le numero de maillage puis pour chaque reference de noeuds
|
||||
# le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites
|
||||
# puis entre crochet la position suivi de la signification de chaque grandeur
|
||||
#
|
||||
# ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti
|
||||
# >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees
|
||||
# puis pour chaque maillage
|
||||
# le numero de maillage puis pour chaque reference d'element et de pti
|
||||
# le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites
|
||||
# puis entre crochet la position suivi de la signification de chaque grandeur
|
||||
#
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie ,
|
||||
# puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees
|
||||
# et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere)
|
||||
# precedes du numero de colonne correspondant entre crochet
|
||||
# puis pour chaque maillage
|
||||
# >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie ,
|
||||
# les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant
|
||||
# a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques,
|
||||
# on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl
|
||||
# chacun precedes du numero de colonne entre crochet
|
||||
# puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres),
|
||||
# suivi entre crochet, de la plage des numeros de colonnes, correspondant
|
||||
# chacun sur une ligne differentes
|
||||
# ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D:
|
||||
# tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2)
|
||||
# en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3)
|
||||
# non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3)
|
||||
# ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements,
|
||||
# ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique
|
||||
# ** ( cf: parametre_style_de_sortie = 0)
|
||||
|
||||
#====================================================================
|
||||
#|| recapitulatif des differentes grandeurs par colonne ||
|
||||
#====================================================================
|
||||
#---------------------------------- grandeur globales ------------------------------------
|
||||
#0 (nombre de grandeurs globales)
|
||||
#---------------------------------- maillage et dimension --------------------------------
|
||||
#1 3 (nombre de maillages et dimension)
|
||||
#---------------------------------- torseurs de reactions --------------------------------
|
||||
#2 12 (nombre de torseurs et nombre total de grandeurs associees)
|
||||
# N_droit [2]Rx [3]Ry [4]Rz [5]Mx [6]My [7]Mz ;
|
||||
# N_gauche [8]Rx [9]Ry [10]Rz [11]Mx [12]My [13]Mz ;
|
||||
#
|
||||
#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref ---------------
|
||||
#0 0 (nombre de ref de noeud et nombre total de grandeurs associees)
|
||||
#
|
||||
#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref ---------------
|
||||
#0 0 (nombre de ref d'element et nombre total de grandeurs associees)
|
||||
#
|
||||
#---------------------------------- grandeurs aux noeuds --------------------------------
|
||||
#0 0 (nombre de noeuds, nombre total de grandeurs associees)
|
||||
#---------------------------------- grandeurs aux elements ------------------------------
|
||||
#1 12 30 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles)
|
||||
# element_1 pt_integ_1: [15]X [16]Y [17]Z [18] SIG11 [19] SIG22 [20] SIG33 [21] SIG12 [22] SIG23 [23] SIG13 [24] EPS11 [25] EPS22 [26] EPS33 [27] EPS12 [28] EPS23 [29] EPS13
|
||||
# [30]X [31]Y [32]Z
|
||||
# EPS_TRANSPORTEE_ANISO EPS_TRANSPORTEE_ANISO TENSEURHH TABLEAU_T TENSEURHH 6 [33...38]
|
||||
# SIGMA_DANS_ANISO SIGMA_DANS_ANISO TENSEURHH TABLEAU_T TENSEURHH 6 [39...44]
|
||||
# PARA_ORTHO PARA_ORTHO VECTEUR TABLEAU_T VECTEUR 9 [45...53]
|
||||
#
|
||||
#
|
||||
#
|
||||
#====================================================================
|
||||
#|| fin du recapitulatif des differentes grandeurs ||
|
||||
#====================================================================
|
||||
|
||||
# ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant
|
||||
# a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun
|
||||
# correspondant a un maillage. On trouve pour chaque enregistrement successivement :
|
||||
# s'il y a des grandeurs globales: le temps puis les grandeurs globales,
|
||||
# puis s'il y a des torseurs de reaction :
|
||||
# de nouveau le temps, les composantes de la resultante puis les composantes du moments
|
||||
# donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels
|
||||
# puis s'il y a des grandeurs aux noeuds: de nouveau le temps
|
||||
# les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud
|
||||
# puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds
|
||||
# puis s'il y a des grandeur aux elements:
|
||||
# le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques)
|
||||
# suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration
|
||||
# correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes
|
||||
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
|
||||
# etc. pour tous les points d'integration - element
|
||||
|
||||
1.000000000000e+00 2.215234918395e-02 -2.017905789551e+02 -4.511946372077e-12 1.513429342164e+03 1.661426183139e-01 -2.014452650275e+04 0.000000000000e+00 2.017952550984e+02 0.000000000000e+00 -1.513464413230e+03 0.000000000000e+00 2.014300631399e+04 1.000000000000e+00 2.215989056159e+01 1.061173242519e+00 1.415559434970e+01 -5.133100620078e+01 -2.516931632485e-02 -4.172265714612e-02 -1.670489143671e+00 -1.204232320746e-01 -2.266976987164e-01 -9.634384841508e-04 3.165771111012e-04 1.284421080200e-04 1.026189098455e-04 4.734003642581e-06 -2.211483746432e-05 2.215989056159e+01 1.061173242519e+00 1.415559434970e+01 -2.292839827067e-04 6.416025895934e-04 -4.192235945520e-04 -1.240645113752e-05 1.897033230560e-05 1.284417799186e-04 -2.768776411507e+01 2.566410358373e+01 -2.373413570679e+01 -2.481290227505e-01 7.588132922240e-02 -4.179567680323e-02 Vecteur taille= 9 1.00000000000000000e+05 5.00000000000000000e+04 2.00000000000000000e+04 2.00000000000000011e-01 2.99999999999999989e-01 1.00000000000000006e-01 2.00000000000000000e+04 1.00000000000000000e+04 2.00000000000000000e+03
|
||||
|
|
@ -0,0 +1,281 @@
|
|||
|
||||
#############################################################################
|
||||
# #
|
||||
# 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 10 x 15 ( 4 x 1 x 2 elements )
|
||||
# geometrie prismatique, decoupage en hexaedre(s), interpolation quadratique complete.
|
||||
|
||||
nom_maillage poutre
|
||||
|
||||
noeuds ------------
|
||||
135 NOEUDS
|
||||
|
||||
|
||||
#---------------------------------------------------------------
|
||||
#|NO DU| X | Y | Z |
|
||||
#|NOEUD| | | |
|
||||
#---------------------------------------------------------------
|
||||
1 0 0 15
|
||||
2 0 5 15
|
||||
3 0 10 15
|
||||
4 0 0 11.25
|
||||
5 0 5 11.25
|
||||
6 0 10 11.25
|
||||
7 0 0 7.5
|
||||
8 0 5 7.5
|
||||
9 0 10 7.5
|
||||
10 0 0 3.75
|
||||
11 0 5 3.75
|
||||
12 0 10 3.75
|
||||
13 0 0 0
|
||||
14 0 5 0
|
||||
15 0 10 0
|
||||
16 12.5 0 15
|
||||
17 12.5 5 15
|
||||
18 12.5 10 15
|
||||
19 12.5 0 11.25
|
||||
20 12.5 5 11.25
|
||||
21 12.5 10 11.25
|
||||
22 12.5 0 7.5
|
||||
23 12.5 5 7.5
|
||||
24 12.5 10 7.5
|
||||
25 12.5 0 3.75
|
||||
26 12.5 5 3.75
|
||||
27 12.5 10 3.75
|
||||
28 12.5 0 0
|
||||
29 12.5 5 0
|
||||
30 12.5 10 0
|
||||
31 25 0 15
|
||||
32 25 5 15
|
||||
33 25 10 15
|
||||
34 25 0 11.25
|
||||
35 25 5 11.25
|
||||
36 25 10 11.25
|
||||
37 25 0 7.5
|
||||
38 25 5 7.5
|
||||
39 25 10 7.5
|
||||
40 25 0 3.75
|
||||
41 25 5 3.75
|
||||
42 25 10 3.75
|
||||
43 25 0 0
|
||||
44 25 5 0
|
||||
45 25 10 0
|
||||
46 37.5 0 15
|
||||
47 37.5 5 15
|
||||
48 37.5 10 15
|
||||
49 37.5 0 11.25
|
||||
50 37.5 5 11.25
|
||||
51 37.5 10 11.25
|
||||
52 37.5 0 7.5
|
||||
53 37.5 5 7.5
|
||||
54 37.5 10 7.5
|
||||
55 37.5 0 3.75
|
||||
56 37.5 5 3.75
|
||||
57 37.5 10 3.75
|
||||
58 37.5 0 0
|
||||
59 37.5 5 0
|
||||
60 37.5 10 0
|
||||
61 50 0 15
|
||||
62 50 5 15
|
||||
63 50 10 15
|
||||
64 50 0 11.25
|
||||
65 50 5 11.25
|
||||
66 50 10 11.25
|
||||
67 50 0 7.5
|
||||
68 50 5 7.5
|
||||
69 50 10 7.5
|
||||
70 50 0 3.75
|
||||
71 50 5 3.75
|
||||
72 50 10 3.75
|
||||
73 50 0 0
|
||||
74 50 5 0
|
||||
75 50 10 0
|
||||
76 62.5 0 15
|
||||
77 62.5 5 15
|
||||
78 62.5 10 15
|
||||
79 62.5 0 11.25
|
||||
80 62.5 5 11.25
|
||||
81 62.5 10 11.25
|
||||
82 62.5 0 7.5
|
||||
83 62.5 5 7.5
|
||||
84 62.5 10 7.5
|
||||
85 62.5 0 3.75
|
||||
86 62.5 5 3.75
|
||||
87 62.5 10 3.75
|
||||
88 62.5 0 0
|
||||
89 62.5 5 0
|
||||
90 62.5 10 0
|
||||
91 75 0 15
|
||||
92 75 5 15
|
||||
93 75 10 15
|
||||
94 75 0 11.25
|
||||
95 75 5 11.25
|
||||
96 75 10 11.25
|
||||
97 75 0 7.5
|
||||
98 75 5 7.5
|
||||
99 75 10 7.5
|
||||
100 75 0 3.75
|
||||
101 75 5 3.75
|
||||
102 75 10 3.75
|
||||
103 75 0 0
|
||||
104 75 5 0
|
||||
105 75 10 0
|
||||
106 87.5 0 15
|
||||
107 87.5 5 15
|
||||
108 87.5 10 15
|
||||
109 87.5 0 11.25
|
||||
110 87.5 5 11.25
|
||||
111 87.5 10 11.25
|
||||
112 87.5 0 7.5
|
||||
113 87.5 5 7.5
|
||||
114 87.5 10 7.5
|
||||
115 87.5 0 3.75
|
||||
116 87.5 5 3.75
|
||||
117 87.5 10 3.75
|
||||
118 87.5 0 0
|
||||
119 87.5 5 0
|
||||
120 87.5 10 0
|
||||
121 100 0 15
|
||||
122 100 5 15
|
||||
123 100 10 15
|
||||
124 100 0 11.25
|
||||
125 100 5 11.25
|
||||
126 100 10 11.25
|
||||
127 100 0 7.5
|
||||
128 100 5 7.5
|
||||
129 100 10 7.5
|
||||
130 100 0 3.75
|
||||
131 100 5 3.75
|
||||
132 100 10 3.75
|
||||
133 100 0 0
|
||||
134 100 5 0
|
||||
135 100 10 0
|
||||
|
||||
# references des noeuds
|
||||
|
||||
N_bas 13 14 15 28 29 30 43 44 45 58 59 60 73 74 75 88
|
||||
89 90 103 104 105 118 119 120 133 134 135
|
||||
N_haut 1 2 3 16 17 18 31 32 33 46 47 48 61 62 63 76
|
||||
77 78 91 92 93 106 107 108 121 122 123
|
||||
N_avant 3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48
|
||||
51 54 57 60 63 66 69 72 75 78 81 84 87 90 93 96
|
||||
99 102 105 108 111 114 117 120 123 126 129 132 135
|
||||
N_arriere 1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46
|
||||
49 52 55 58 61 64 67 70 73 76 79 82 85 88 91 94
|
||||
97 100 103 106 109 112 115 118 121 124 127 130 133
|
||||
N_droit 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
|
||||
N_gauche 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135
|
||||
N_avant_droit 3 6 9 12 15
|
||||
N_avant_gauche 123 126 129 132 135
|
||||
N_haut_avant 3 18 33 48 63 78 93 108 123
|
||||
N_bas_avant 15 30 45 60 75 90 105 120 135
|
||||
N_arriere_droit 1 4 7 10 13
|
||||
N_arriere_gauche 121 124 127 130 133
|
||||
N_haut_arriere 1 16 31 46 61 76 91 106 121
|
||||
N_bas_arriere 13 28 43 58 73 88 103 118 133
|
||||
N_haut_droit 1 2 3
|
||||
N_haut_gauche 121 122 123
|
||||
N_bas_droit 13 14 15
|
||||
N_bas_gauche 133 134 135
|
||||
N_haut_avant_droit 3
|
||||
N_haut_avant_gauche 123
|
||||
N_bas_avant_gauche 135
|
||||
N_bas_avant_droit 15
|
||||
N_haut_arriere_droit 1
|
||||
N_haut_arriere_gauche 121
|
||||
N_bas_arriere_gauche 133
|
||||
N_bas_arriere_droit 13
|
||||
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
|
||||
|
||||
|
||||
# les elements
|
||||
|
||||
elements ----------
|
||||
8 ELEMENTS
|
||||
#----------------------------------------------------------------------
|
||||
#| NO | | |
|
||||
#|ELTS | type element | Noeuds |
|
||||
#----------------------------------------------------------------------
|
||||
1 HEXAEDRE QUADRACOMPL _cm27pti 31 37 7 1 33 39 9 3 34 22 4 16 32 38 8 2 36 24 6 18 19 35 23 5 17 21 20
|
||||
2 HEXAEDRE QUADRACOMPL _cm27pti 37 43 13 7 39 45 15 9 40 28 10 22 38 44 14 8 42 30 12 24 25 41 29 11 23 27 26
|
||||
3 HEXAEDRE QUADRACOMPL _cm27pti 61 67 37 31 63 69 39 33 64 52 34 46 62 68 38 32 66 54 36 48 49 65 53 35 47 51 50
|
||||
4 HEXAEDRE QUADRACOMPL _cm27pti 67 73 43 37 69 75 45 39 70 58 40 52 68 74 44 38 72 60 42 54 55 71 59 41 53 57 56
|
||||
5 HEXAEDRE QUADRACOMPL _cm27pti 91 97 67 61 93 99 69 63 94 82 64 76 92 98 68 62 96 84 66 78 79 95 83 65 77 81 80
|
||||
6 HEXAEDRE QUADRACOMPL _cm27pti 97 103 73 67 99 105 75 69 100 88 70 82 98 104 74 68 102 90 72 84 85 101 89 71 83 87 86
|
||||
7 HEXAEDRE QUADRACOMPL _cm27pti 121 127 97 91 123 129 99 93 124 112 94 106 122 128 98 92 126 114 96 108 109 125 113 95 107 111 110
|
||||
8 HEXAEDRE QUADRACOMPL _cm27pti 127 133 103 97 129 135 105 99 130 118 100 112 128 134 104 98 132 120 102 114 115 131 119 101 113 117 116
|
||||
|
||||
|
||||
|
||||
# references des elements
|
||||
|
||||
E_bas 2 4 6 8
|
||||
E_haut 1 3 5 7
|
||||
E_avant 1 2 3 4 5 6 7 8
|
||||
E_arriere 1 2 3 4 5 6 7 8
|
||||
E_droit 1 2
|
||||
E_gauche 7 8
|
||||
E_avant_droit 1 2
|
||||
E_avant_gauche 7 8
|
||||
E_haut_avant 1 3 5 7
|
||||
E_bas_avant 2 4 6 8
|
||||
E_arriere_droit 1 2
|
||||
E_arriere_gauche 7 8
|
||||
E_haut_arriere 1 3 5 7
|
||||
E_bas_arriere 2 4 6 8
|
||||
E_haut_droit 1
|
||||
E_haut_gauche 7
|
||||
E_bas_droit 2
|
||||
E_bas_gauche 8
|
||||
E_haut_avant_droit 1
|
||||
E_haut_avant_gauche 7
|
||||
E_bas_avant_gauche 8
|
||||
E_bas_avant_droit 2
|
||||
E_haut_arriere_droit 1
|
||||
E_haut_arriere_gauche 7
|
||||
E_bas_arriere_gauche 8
|
||||
E_bas_arriere_droit 2
|
||||
E_to 1 2 3 4 5 6 7 8
|
||||
|
||||
# references des faces et arretes d'elements
|
||||
|
||||
F_bas 2 5 4 5 6 5 8 5
|
||||
F_haut 1 2 3 2 5 2 7 2
|
||||
F_avant 1 4 2 4 3 4 4 4 5 4 6 4 7 4 8 4
|
||||
F_arriere 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1
|
||||
F_droit 1 6 2 6
|
||||
F_gauche 7 3 8 3
|
||||
A_avant_droit 1 11 2 11
|
||||
A_avant_gauche 7 9 8 9
|
||||
A_haut_avant 1 12 3 12 5 12 7 12
|
||||
A_bas_avant 2 10 4 10 6 10 8 10
|
||||
A_arriere_droit 1 3 2 3
|
||||
A_arriere_gauche 7 1 8 1
|
||||
A_haut_arriere 1 4 3 4 5 4 7 4
|
||||
A_bas_arriere 2 2 4 2 6 2 8 2
|
||||
A_haut_droit 1 8
|
||||
A_haut_gauche 7 5
|
||||
A_bas_droit 2 7
|
||||
A_bas_gauche 8 6
|
||||
|
||||
Loading…
Reference in a new issue