diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/README b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/README new file mode 100644 index 0000000..4c51e6f --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/README @@ -0,0 +1,32 @@ +------------------------------------------------------ +Auteur +------------------------------------------------------ +Gerard Rio (gerard.rio@univ-ubs.fr) + +------------------------------------------------------ +Mots-cles +------------------------------------------------------ +_princ.maple +sortie tenseur + +------------------------------------------------------------ +But du test +------------------------------------------------------------ + Test de bon fonctionnement de la sortie des composantes du tenseur de déformation et de contrainte. +L'idée est de vérifier certaines méthodes de changement de base. +"" +TYPE_DE_CALCUL +non_dynamique +"" + +------------------------------------------------------------ +Description du calcul +------------------------------------------------------------ +Le calcul s'effectue ici en 3D sur 3 maillages en parallèle: un hexaèdre, une plaque et une barre. +On impose un déplacement à tous les noeuds ce qui permet d'imposer le tenseur de déformation. +La longueur où on test est 10 et on déplace de 5 d'où une déformation de 0.27778 + +------------------------------------------------------------ +Grandeurs de comparaison +------------------------------------------------------------ +On vérifie le déplacement la composantes 11 du tenseur de déformation et de contrainte pour les 3 maillages. diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/barre.her b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/barre.her new file mode 100644 index 0000000..4a26aa5 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/barre.her @@ -0,0 +1,61 @@ + + + nom_maillage barre +############################################################################# +# # +# 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 d'element lineaire 1D | + # | pour une utilisation par herezh++ | + # ------------------------------------------------------------------ + + # barre de dimension : 10 ( 1 elements ) + + noeuds ------------ + 2 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 10 0 0 + + # references des noeuds + + N_deb 1 + N_fi 2 + N_to 1 2 + + + # les elements + + elements ---------- + 1 ELEMENTS + #-------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #-------------------------------------------------------------------- + 1 POUT BIE1 1 2 + + + + # references des arretes + + A_to 1 1 + + + # references des elements + + E_to 1 + diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/cube.her b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/cube.her new file mode 100644 index 0000000..4e036cc --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/cube.her @@ -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 : 10 x 20 x 30 ( 1 x 1 x 1 elements ) + # geometrie prismatique, decoupage en hexaedre(s), interpolation lineaire. + + nom_maillage cube + + noeuds ------------ + 8 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 0 30 + 3 0 20 0 + 4 0 20 30 + 5 10 0 0 + 6 10 0 30 + 7 10 20 0 + 8 10 20 30 + + # 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 + diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/plaque.her b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/plaque.her new file mode 100644 index 0000000..76ff04a --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/plaque.her @@ -0,0 +1,81 @@ + +############################################################################# +# # +# 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 : 10 x 5 ( 1 x 1 elements ) + # geometrie rectangulaire, decoupage triangulaire, interpolation lineaire. + + nom_maillage plaque + + noeuds ------------ + 4 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 5 0 + 3 10 0 0 + 4 10 5 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 ---------- + 2 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 TRIANGLE LINEAIRE 1 3 2 + 2 TRIANGLE LINEAIRE 3 4 2 + + + # references des elements + +A_O 1 3 +A_S 1 1 +A_E 2 1 +A_N 2 2 +F_to 1 1 2 1 +E_O 1 +E_S 1 +E_E 2 +E_N 2 +E_SO 1 +E_SE 1 2 +E_NE 2 +E_NO 1 2 +F_O 1 1 2 1 +F_S 1 1 2 1 +F_E 1 1 2 1 +F_N 1 1 2 1 + E_to 1 2 + diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.CVisu1 b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.CVisu1 new file mode 100644 index 0000000..e358b42 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.CVisu1 @@ -0,0 +1,377 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.842 # +# Copyright (c) 1997-2018, 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: + # 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 , sur une meme ligne + 1 2 3 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 , + # 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 , puis + debut_list_torseur_reaction 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 , puis + debut_list_SM_sur_refN fin_list__SM_sur_refN + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 + 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 , le pti associe + # puis + 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 , puis + debut_list_SM_sur_refpti fin_list__SM_sur_refpti + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + 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 + 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 + 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 + deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + 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 + deb_list_ddl_element SIG11 EPS11 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 + 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 + 2 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction 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 , puis + debut_list_SM_sur_refN fin_list__SM_sur_refN + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 + 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 , le pti associe + # puis + 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 , puis + debut_list_SM_sur_refpti fin_list__SM_sur_refpti + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + 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 + 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 + 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 + deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + 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 + deb_list_ddl_element SIG11 EPS11 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 + 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 + 3 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction 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 , puis + debut_list_SM_sur_refN fin_list__SM_sur_refN + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 + 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 , le pti associe + # puis + 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 , puis + debut_list_SM_sur_refpti fin_list__SM_sur_refpti + # debut de la liste des ddl a considerer , (une liste de ddl), puis + 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 + 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 , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + 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 + 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 + 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 + deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + 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 + deb_list_ddl_element SIG11 EPS11 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 + 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 (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 || + # ================================================================================= + + + # ================================================================================= + # || ***** demande d'une visualisation Gmsh: ***** || + # ================================================================================= + # un mot cle de debut (debut_visualisation_Gmsh) + # un mot cle de fin ( fin_visualisation_Gmsh) apres tous les ordres particuliers + # la seule presence du premier mots cle suffit a activer la visualisation Gmsh + # la presence du second permet une meilleur lisibilite du fichier, mais n'est pas indispensable + debut_visualisation_Gmsh + + # ----------------------------- definition des parametres du maillage initial: ---------------- + debut_maillage_initial # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + pseudo-homothetie_sur_les_maillages_ 0 # 0 = aucune homothetie, 1 = il y en a + # --- def eventuelle de la pseudo-homothetie: une par maillage, + # pseudo-homothetie pour le maillage : 1 +# maillage_ 1 + # mot cle: maillage_ puis le numero du maillage, +# pseudo-homothetie_ 0 # 0 = non active, 1 = active + # ensuite si c'est active, on trouve : + # mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie + # puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes. +# centre_homothetie_ +# fact_mult_ + # pseudo-homothetie pour le maillage : 2 +# maillage_ 2 + # mot cle: maillage_ puis le numero du maillage, +# pseudo-homothetie_ 0 # 0 = non active, 1 = active + # ensuite si c'est active, on trouve : + # mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie + # puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes. +# centre_homothetie_ +# fact_mult_ + # pseudo-homothetie pour le maillage : 3 +# maillage_ 3 + # mot cle: maillage_ puis le numero du maillage, +# pseudo-homothetie_ 0 # 0 = non active, 1 = active + # ensuite si c'est active, on trouve : + # mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie + # puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes. +# centre_homothetie_ +# fact_mult_ + visualisation_references_sur_les_maillages_ 1 # 0 = pas de visualisation des reference, 1 = sortie des ref dans fichier unique 2= sortie des ref dans plusieurs fichiers + fin_maillage_initial # le mot cle de fin + + # ----------------------------- definition des parametres pour les isovaleurs : ---------------- + debut_isovaleur_Gmsh # mot cle de debut des parametres pour les isovaleurs + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + ancien_format_ 0 # 1 = ancien format, 0 = nouveau format +# 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_ 1 + 1 # le numero de maillage + # tableau des ddl aux noeuds a visualiser, un par maillage + debut_tableau_ddl_aux_noeuds fin_tableau_ddl_aux_noeuds + # tableau des choix_var aux noeuds a visualiser, un par maillage + # choix_var (=1 ou 0) indique si oui ou non il s'agit de la variation + debut_tableau_choix_var_ddl_aux_noeuds fin_tableau_choix_var_ddl_aux_noeuds + # tableau des ddl_etendu aux noeuds a visualiser, un par maillage + debut_tableau_ddl_etendu_aux_noeuds fin_tableau_ddl_etendu_aux_noeuds + # tableau de grandeurs evoluees aux noeuds a visualiser, un par maillage + deb_list_GrandEvoluee_noeud fin_list_GrandEvoluee_noeud + # tableau de ddl aux elements a visualiser, un par maillage + debut_tableau_ddl_aux_elements SIG11 EPS11 fin_tableau_ddl_aux_elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element CONTRAINTE_COURANTE DEFORMATION_COURANTE fin_list_GrandEvoluee_element + # tableau de grandeurs particulieres aux elements a visualiser, un par maillage + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element + 2 # le numero de maillage + # tableau des ddl aux noeuds a visualiser, un par maillage + debut_tableau_ddl_aux_noeuds fin_tableau_ddl_aux_noeuds + # tableau des choix_var aux noeuds a visualiser, un par maillage + # choix_var (=1 ou 0) indique si oui ou non il s'agit de la variation + debut_tableau_choix_var_ddl_aux_noeuds fin_tableau_choix_var_ddl_aux_noeuds + # tableau des ddl_etendu aux noeuds a visualiser, un par maillage + debut_tableau_ddl_etendu_aux_noeuds fin_tableau_ddl_etendu_aux_noeuds + # tableau de grandeurs evoluees aux noeuds a visualiser, un par maillage + deb_list_GrandEvoluee_noeud fin_list_GrandEvoluee_noeud + # tableau de ddl aux elements a visualiser, un par maillage + debut_tableau_ddl_aux_elements SIG11 EPS11 fin_tableau_ddl_aux_elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element CONTRAINTE_COURANTE DEFORMATION_COURANTE fin_list_GrandEvoluee_element + # tableau de grandeurs particulieres aux elements a visualiser, un par maillage + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element + 3 # le numero de maillage + # tableau des ddl aux noeuds a visualiser, un par maillage + debut_tableau_ddl_aux_noeuds fin_tableau_ddl_aux_noeuds + # tableau des choix_var aux noeuds a visualiser, un par maillage + # choix_var (=1 ou 0) indique si oui ou non il s'agit de la variation + debut_tableau_choix_var_ddl_aux_noeuds fin_tableau_choix_var_ddl_aux_noeuds + # tableau des ddl_etendu aux noeuds a visualiser, un par maillage + debut_tableau_ddl_etendu_aux_noeuds fin_tableau_ddl_etendu_aux_noeuds + # tableau de grandeurs evoluees aux noeuds a visualiser, un par maillage + deb_list_GrandEvoluee_noeud fin_list_GrandEvoluee_noeud + # tableau de ddl aux elements a visualiser, un par maillage + debut_tableau_ddl_aux_elements SIG11 EPS11 fin_tableau_ddl_aux_elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element CONTRAINTE_COURANTE DEFORMATION_COURANTE fin_list_GrandEvoluee_element + # tableau de grandeurs particulieres aux elements a visualiser, un par maillage + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element + fin_isovaleur_Gmsh # mot cle de fin des parametres pour les isovaleurs + + # ----------------------------- definition des parametres de deformee: ---------------- + debut_deformee # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # definition des alertes: deb_list_alerte + # un mot clef de debut + # puis deux nombres: un mini et un maxi, et un nom + # un mot clef de fin: fin_list_alerte + deb_list_alerte + fin_list_alerte + avec_vitesse= 0 avec_acceleration= 0 # 1 ou 0 pour la sortie si disponible + fin_deformee # un mot cle de fin + + # ----------------------------- definition de la liste des increments a balayer: ---------------- + debut_list_increment # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments) + # un mot cle de fin de liste ( fin_list_increment) + tous_les_increments_moins_zero dernier_increment fin_list_increment + + # ----------------------------- choix des maillages a visualiser: ---------------- + # la liste est facultative, par defaut la visualisation concerne le premier maillage + debut_choix_maillage # un mot cle de debut, + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entiers , puis , sur une meme ligne + 1 2 3 fin_choix_maillage + + fin_visualisation_Gmsh + # ================================================================================= + # || fin de la visualisation Gmsh || + # ================================================================================= + + + + + + fin_fichier_commande_visu # <<<<<< le mot cle 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 + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.info b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.info new file mode 100644 index 0000000..ddb6c3d --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur.info @@ -0,0 +1,159 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2018 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.842 + + +####################################################################### +# fichier de commande # + + +# --- test élémentaire pour sortir les composantes du tenseur des contraintes et def ---- +# vérification de certaines méthodes de changements de base interne + + +#---------------------------------------------------- +# 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) | +#-------------------------------------- + +# un maillage 3D +< cube.her + +# un maillage 2D +< plaque.her +def_mouvement_solide_initiaux_ + mouvement_solide_ # def de mouvements solides + translation_= 0. 0. 35 + fin_mouvement_solide_ + +# un maillage 1D +< barre.her +def_mouvement_solide_initiaux_ + mouvement_solide_ # def de mouvements solides + translation_= 0. 0. 40 + fin_mouvement_solide_ + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +nom_mail= cube E_tout acier +nom_mail= plaque E_tout acier2D +nom_mail= barre E_tout acier1D + + materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ + acier ISOELAS + 200000 0.3 + + acier2D ISOELAS2D_C + 200000 0.3 + + acier1D ISOELAS1D + 200000 0.3 + + +#--------------------------------------- fin def des lois de comportement ------------------------ + +# --- divers stockages (1) ------- + masse_volumique #--------# +nom_mail= cube E_tout 1 +nom_mail= plaque E_tout 1 +nom_mail= barre E_tout 1 + +epaisseurs #--------# +nom_mail= plaque E_tout 1 + +sections #--------# +nom_mail= barre E_tout 1 + + + charges #------------# + + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +nom_mail= cube N_tout UX UY UZ +nom_mail= plaque N_tout UX UY UZ +nom_mail= barre N_tout UX UY UZ + +nom_mail= cube N_gauche 'UX= 5' +nom_mail= plaque N_E 'UX= 5' +nom_mail= barre N_fi 'UX= 5' + + + + controle ------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +TEMPSFIN 1.#2.#3. +#DELTAt 0.1 +#DELTAtMAXI 0.1 +#DELTAtMINI 0.1 +SAUVEGARDE 1#INTER_TEMPS 0.1 +ITERATIONS 15 +PRECISION 1e-3 +NORME min(Res,Res/Reaction) +MAXINCRE 3#400. + + + + 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_ + + + diff --git a/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur_princ.maple.ref1 new file mode 100644 index 0000000..28e2b60 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/SORTIES/Test_R_sortie_tenseur/sortie_tenseur_princ.maple.ref1 @@ -0,0 +1,123 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.842 # +# Copyright (c) 1997-2018, 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 -------------------------------- +#3 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#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) +# +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#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) +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# maillage_1 : +# element_1 pt_integ_1: [2]X [3]Y [4]Z [5] SIG11 [6] EPS11 +# +# +#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# maillage_2 : +# element_1 pt_integ_1: [8]X [9]Y [10]Z [11] SIG11 [12] EPS11 +# +# +#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# maillage_3 : +# element_1 pt_integ_1: [14]X [15]Y [16]Z [17] SIG11 [18] EPS11 +# +# +#==================================================================== +#|| 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 3.169872981078e+00 1.577350269190e+01 2.366025403784e+01 7.478632478632e+04 2.777777777778e-01 1.000000000000e+00 5.000000000000e+00 1.666666666667e+00 3.500000000000e+01 6.105006105006e+04 2.777777777778e-01 1.000000000000e+00 7.500000000000e+00 0.000000000000e+00 4.000000000000e+01 5.555555555556e+04 2.777777777778e-01