From 8bd0e7c60455de7b8b73a2c83799c39ed1affdec Mon Sep 17 00:00:00 2001 From: troufflard <> Date: Tue, 11 Oct 2016 21:28:51 +0000 Subject: [PATCH] ajout 2 tests en attente debug (2 cas tests en lien avec le probleme de numerotation des increments lors d un _suite_point_info_) : Test_R_non_dynamique_suite_point_info/ et Test_R_sortie_suite_point_info/; voir les README (attention : Test_R_sortie_suite_point_info/ est un test deja existant dans Batterie/; il s agit de le modifier. C est explique dans son README) --- .../README | 61 ++++++ .../cube_1elt_npti8.her | 18 ++ .../cube_1elt_npti8.lis | 158 ++++++++++++++ .../non_dynamique_suite_point_info.CVisu1 | 116 +++++++++++ .../non_dynamique_suite_point_info.info | 82 ++++++++ .../non_dynamique_suite_point_info.verif1 | 61 ++++++ ...ynamique_suite_point_info_princ.maple.ref1 | 196 ++++++++++++++++++ .../Test_R_sortie_suite_point_info/README | 72 +++++++ .../sortie_suite_point_info.CVisu1 | 136 ++++++++++++ .../sortie_suite_point_info.info | 108 ++++++++++ .../sortie_suite_point_info.verif1 | 61 ++++++ .../sortie_suite_point_info_princ.maple.ref1 | 104 ++++++++++ .../tranche_sphere.her | 140 +++++++++++++ 13 files changed, 1313 insertions(+) create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/README create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.her create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.lis create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.CVisu1 create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.info create mode 100755 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.verif1 create mode 100644 Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info_princ.maple.ref1 create mode 100755 Tests_en_attente_debug/Test_R_sortie_suite_point_info/README create mode 100755 Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.CVisu1 create mode 100755 Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.info create mode 100755 Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.verif1 create mode 100644 Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info_princ.maple.ref1 create mode 100755 Tests_en_attente_debug/Test_R_sortie_suite_point_info/tranche_sphere.her diff --git a/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/README b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/README new file mode 100644 index 0000000..d6e8381 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/README @@ -0,0 +1,61 @@ +------------------------------------------------------ +Auteur +------------------------------------------------------ +Julien Troufflard (troufflard@univ-ubs.fr) + + +------------------------------------------------------ +Mots-cles +------------------------------------------------------ +_suite_point_info_ +non_dynamique +traction uniaxiale +ISOELAS + + + +------------------------------------------------------ +But du test +------------------------------------------------------ +Vérification du nombre et de la liste des incréments calculés (INCREMENT DE CHARGE) dans le cas d une continuation de calcul (_suite_point_info_) + +@@@@@@@@@@@@@@@@@ +@@@ A EFFACER @@@ +@@@@@@@@@@@@@@@@@ +@ +@ probleme similaire que Test_R_sortie_suite_point_info/ (mais pour 20 increments) +@ +@ contrairement a Test_R_sortie_suite_point_info/, le cas ici présent de non_dynamique +@ conduit a avoir une liste d increments finale : +@ 1, 2, ..., 10 et ensuite 10, 11, ..., 19 (le 10 est répété 2 fois et le dernier +@ increment est egal a 19) +@ +@ +@ ATTENTION : +@ contrairement a Test_R_sortie_suite_point_info/, ce cas test est nouveau et a +@ ajouter completement a la batterie (dans FONCTION_UNIQUE/SORTIES par exemple, car +@ on peut partir du principe que l affichage Herezh est une sorte de "sortie".) +@@@@@@@@@@@@@@@@@ + +------------------------------------------------------ +Description du calcul +------------------------------------------------------ +traction uniaxiale suivant X sur un cube 1x1x1 (1 element HEXAEDRE LINEAIRE). Loi ISOELAS + +Le calcul se fait en 2 étapes : +1) premier calcul entre le temps 0 et 1 (traction UX de 0.01) +2) poursuite du calcul (_suite_point_info_) jusqu a un temps de 2 (traction UX de 0.02). Le déplacement imposé est piloté par une courbe de charge. + + +------------------------------------------------------------ +Informations sur les fichiers facultatifs +------------------------------------------------------------ +Le script non_dynamique_suite_point_info.verif1 lit le fichier de redirection de l affichage produit par Herezh. Il renvoie OK si il trouve un affichage de 20 incréments de calcul (INCREMENT DE CHARGE) et que ces incréments sont numerotés dans l ordre de 1 a 20. + + +------------------------------------------------------------- +Grandeurs de comparaison +------------------------------------------------------------- +pour le point d integration 1 de l element 1 : + - contrainte : SIG11 + - déformations : EPS11, EPS22, EPS33 diff --git a/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.her b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.her new file mode 100644 index 0000000..3bc9991 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.her @@ -0,0 +1,18 @@ + nom_maillage cube_1elt_npti8 + + noeuds + 8 NOEUDS + + 1 0 0 0 + 2 0 0 1 + 3 0 1 0 + 4 0 1 1 + 5 1 0 0 + 6 1 0 1 + 7 1 1 0 + 8 1 1 1 + + elements + 1 ELEMENTS + + 1 HEXAEDRE LINEAIRE 5 7 3 1 6 8 4 2 diff --git a/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.lis b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.lis new file mode 100644 index 0000000..8224c23 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/cube_1elt_npti8.lis @@ -0,0 +1,158 @@ +# CORRESPONDANCE label => position +# +# gauche => X=1 +# +# droit => X=0 +# +# avant => Y=1 +# +# arriere => Y=0 +# +# haut => Z=1 +# +# bas => Z=0 +# + + 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 + + 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 + + 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/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.CVisu1 b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.CVisu1 new file mode 100644 index 0000000..63617b9 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.CVisu1 @@ -0,0 +1,116 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.687 # +# Copyright (c) 1997-2015, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + debut_fichier_commande_visu # >>>>>> le mot cle: + # 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) + 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 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 + 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 EPS11 EPS22 EPS33 SIG11 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 || + # ================================================================================= + + + + + + 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/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.info b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.info new file mode 100644 index 0000000..6cb81b6 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.info @@ -0,0 +1,82 @@ +dimension 3 + + +niveau_commentaire 3 + + +TYPE_DE_CALCUL + +non_dynamique avec plus lectureCommandesVisu + + +< cube_1elt_npti8.her +< cube_1elt_npti8.lis + + +les_courbes_1D + evol_UX COURBEPOLYLINEAIRE_1_D + Debut_des_coordonnees_des_points + Coordonnee dim= 2 0. 0. + Coordonnee dim= 2 1. 0.01 + Coordonnee dim= 2 2. 0.02 + Fin_des_coordonnees_des_points + + +choix_materiaux + +E_tout MAT + +materiaux + +MAT ISOELAS +1000. 0.4 + + +masse_volumique +E_tout 1. + + +charges + + +blocages + +#traction uniaxiale suivant X par deplacement impose +N_droit UX +N_arriere UY +N_bas UZ +N_gauche 'UX= COURBE_CHARGE: evol_UX ECHELLE: 1' + + +controle + +DELTAtMINI 0.1 +DELTAtMAXI 0.1 +TEMPSFIN 1 +DELTAt 0.1 +ITERATIONS 15 +PRECISION 1e-3 +SAUVEGARDE DERNIER_CALCUL +MAXINCRE 99999999 +NORME Residu/Reaction_et_VarRes + + +para_pilotage_equi_global + + +para_syteme_lineaire + + +para_affichage + + +resultats pas_de_sortie_finale_ +COPIE 0 + + +_suite_point_info_ + +controle +TEMPSFIN 2. + +_fin_point_info_ diff --git a/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.verif1 b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.verif1 new file mode 100755 index 0000000..f4fc853 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info.verif1 @@ -0,0 +1,61 @@ +#!/usr/bin/env perl +use strict; +use warnings; +use English; +use File::Basename; +my $NOM_PROG = basename $PROGRAM_NAME; + + +#executable Herezh +my $exeHZ = shift(@ARGV); + +#fichier redirection affichage +my $flog = shift(@ARGV); + +#fichier .info +my $finfo = shift(@ARGV); + + +# +# +# verification des numeros d increment dans le fichier $flog (INCREMENT_DE_CHARGE) +# +# pour ce calcul, on s attend a trouver 20 increments de numeros : 1, 2, 3, ..., 20 +# +# ce script renvoie OK si il y a 20 increments numerotes de 1 a 20 +# sinon ECHEC +# +# + +#liste des increments +my @liste_increments; + +open(FIC, "<$flog"); +while() { + next if(not /^\s*INCREMENT\s+DE\s+CHARGE\s+:\s+(\d+)\b/o); + push(@liste_increments, $1);#ajout du numero increment dans la liste +} +close(FIC); + +#nombre d increments +my $nb_incr = $#liste_increments + 1; + + +#ECHEC si le nombre d increments n est pas egal a 20 +if($nb_incr != 20) { + print "\nErreur (prog:$NOM_PROG) : $nb_incr increments ont ete calcules (on s attendait a 20 increments)\n\n"; + print "resultat verification : ECHEC $nb_incr increments ont ete calcules (on s attendait a 20 increments)\n"; + exit; +} + +#ECHEC si la liste n est pas dans l ordre exact 1 2 3 ... 20 +for(1 .. 20) { + if($liste_increments[$_-1] != $_) { + print "\nErreur (prog:$NOM_PROG) : probleme dans les numeros des increments\n\n"; + print "resultat verification : ECHEC probleme dans les numeros des increments\n"; + exit; + } +} + +#OK : aucun probleme n a ete rencontre +print "resultat verification : OK\n"; diff --git a/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info_princ.maple.ref1 b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info_princ.maple.ref1 new file mode 100644 index 0000000..3bfd28a --- /dev/null +++ b/Tests_en_attente_debug/Test_R_non_dynamique_suite_point_info/non_dynamique_suite_point_info_princ.maple.ref1 @@ -0,0 +1,196 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.775 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + # entete des donnees : informations gererales: on trouve successivement: + # >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs + # precedes du numero de colonne entre crochet + # >> le nombre de maillages m, et dimension de l'espace de travail + # puis pour chaque maillage, + # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit + # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions + # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # puis pour chaque maillage + # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , + # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees + # et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere) + # precedes du numero de colonne correspondant entre crochet + # puis pour chaque maillage + # >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant + # a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques, + # on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl + # chacun precedes du numero de colonne entre crochet + # puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres), + # suivi entre crochet, de la plage des numeros de colonnes, correspondant + # chacun sur une ligne differentes + # ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D: + # tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2) + # en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3) + # non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3) + # ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements, + # ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique + # ** ( cf: parametre_style_de_sortie = 0) + +#==================================================================== +#|| recapitulatif des differentes grandeurs par colonne || +#==================================================================== +#---------------------------------- grandeur globales ------------------------------------ +#0 (nombre de grandeurs globales) +#---------------------------------- maillage et dimension -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#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_1 pt_integ_1: [2]X [3]Y [4]Z [5] EPS11 [6] EPS22 [7] EPS33 [8] SIG11 +# +# +#==================================================================== +#|| 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.134381140592e-01 7.855853835505e-01 7.855853835505e-01 9.851975296540e-03 -3.940790130020e-03 -3.940790130020e-03 9.851975263956e+00 +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.775 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + # entete des donnees : informations gererales: on trouve successivement: + # >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs + # precedes du numero de colonne entre crochet + # >> le nombre de maillages m, et dimension de l'espace de travail + # puis pour chaque maillage, + # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit + # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions + # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # puis pour chaque maillage + # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , + # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees + # et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere) + # precedes du numero de colonne correspondant entre crochet + # puis pour chaque maillage + # >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant + # a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques, + # on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl + # chacun precedes du numero de colonne entre crochet + # puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres), + # suivi entre crochet, de la plage des numeros de colonnes, correspondant + # chacun sur une ligne differentes + # ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D: + # tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2) + # en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3) + # non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3) + # ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements, + # ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique + # ** ( cf: parametre_style_de_sortie = 0) + +#==================================================================== +#|| recapitulatif des differentes grandeurs par colonne || +#==================================================================== +#---------------------------------- grandeur globales ------------------------------------ +#0 (nombre de grandeurs globales) +#---------------------------------- maillage et dimension -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#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_1 pt_integ_1: [2]X [3]Y [4]Z [5] EPS11 [6] EPS22 [7] EPS33 [8] SIG11 +# +# +#==================================================================== +#|| 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 + + 2.000000000000e+00 2.155513627133e-01 7.826231499345e-01 7.826231499345e-01 1.941560938101e-02 -7.762847935520e-03 -7.762847935522e-03 1.942531171496e+01 diff --git a/Tests_en_attente_debug/Test_R_sortie_suite_point_info/README b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/README new file mode 100755 index 0000000..645feeb --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/README @@ -0,0 +1,72 @@ +------------------------------------------------------ +Auteur +------------------------------------------------------ +Gérard Rio (gerard.rio@univ-ubs.fr) +Frank Petitjean (frank.petitjean@rtime.fr) +Julien Troufflard (troufflard@univ-ubs.fr) + +------------------------------------------------------ +Mots-cles +------------------------------------------------------ +_suite_point_info_ +FREQUENCE_SORTIE_FIL_DU_CALCUL +CAL_VOL_TOTAL_ENTRE_SURFACE_ET_PLANS_REF +relaxation dynamique +dynamique_relaxation_dynam + + +------------------------------------------------------ +But du test +------------------------------------------------------ +Vérification du fonctionnement des grandeurs en sortie et vérification du nombre et de la liste des incréments calculés (INCREMENT DE CHARGE) dans le cas d une continuation de calcul (_suite_point_info_) + +@@@@@@@@@@@@@@@@@ +@@@ A EFFACER @@@ +@@@@@@@@@@@@@@@@@ +@ +@ pas de probleme sur la comparaison .maple +@ +@ probleme sur la verif .verif (ce calcul fait 6 increments. Le .verif lit le .log +@ et constate que les numeros d increment 1 2 3 4 5 6 +@ ne sont pas tous presents => ECHEC. +@ dans le cas de la relaxation dynamique, on obtient les +@ increments : 1 2, 1 2, 1 2 !!!????!!!) +@ +@ +@ ATTENTION : +@ il ne s agit pas d un nouveau test. Il est deja present dans +@ FONCTION_UNIQUE/SORTIES/Test_R_sortie_suite_point_info/. Il n est pas question de remplacer entierement +@ ce repertoire de test. +@ le but est de regler le bug decrit ci-dessus et ensuite faire les 2 modifs +@ suivantes : +@ 1) remplacer le README du repertoire FONCTION_UNIQUE/SORTIES/Test_R_sortie_suite_point_info/ par ce README +@ 2) ajouter le .verif1 dans le repertoire FONCTION_UNIQUE/SORTIES/Test_R_sortie_suite_point_info/ +@ +@@@@@@@@@@@@@@@@@ + + + +------------------------------------------------------ +Description du calcul +------------------------------------------------------ +Le maillage est une tranche de sphère de rayon 10m. Tous les noeuds sont bloqués selon UX, UY et UZ (aucune calcul mécanique). On calcule le volume inclus dans cette tranche de sphère (mot-clé para_calculs_geometriques). Le premier calcul (temps=1) est suivi de deux _suite_point_info_ (temps=2 et temps=3) avec à chaque fois 2 incréments. + +\latex_brut +Test implanté suite à la modif Herezh en réponse à la demande \#115 sur le site d Herezh : + + \url{http://kleger.univ-ubs.fr/Herezh/issues/115} (voir ce lien pour plus d infos) + +\fin_latex_brut + + +------------------------------------------------------------ +Informations sur les fichiers facultatifs +------------------------------------------------------------ +Le script non_dynamique_suite_point_info.verif1 lit le fichier de redirection de l affichage produit par Herezh. Il renvoie OK si il trouve un affichage de 20 incréments de calcul (INCREMENT DE CHARGE) et que ces incréments sont numerotés dans l ordre de 1 a 20. + + +------------------------------------------------------------- +Grandeurs de comparaison +------------------------------------------------------------- +Diverses grandeurs globales (notamment les volumes par rapport aux plan xy, xz et yz) +chargement extérieur sur les noeuds 1 23 44 (pression_ext) diff --git a/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.CVisu1 b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.CVisu1 new file mode 100755 index 0000000..58486ac --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.CVisu1 @@ -0,0 +1,136 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.765 # +# Copyright (c) 1997-2016, 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) + 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 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 energie_cinetique energie_elastique energie_interne vol_total2D_avec_plan_xy vol_total2D_avec_plan_xz vol_total2D_avec_plan_yz 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 1 23 44 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 pression_ext 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 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 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 || + # ================================================================================= + + + + 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/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.info b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.info new file mode 100755 index 0000000..d160e45 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.info @@ -0,0 +1,108 @@ +dimension 3 + +niveau_commentaire 3 + + +TYPE_DE_CALCUL +dynamique_relaxation_dynam + +PARA_TYPE_DE_CALCUL +typeCalRelaxation= 1 lambda= 0.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.005 +nb_max_dX_OK_pourRelaxDyn_ 10 +nb_deb_testfin_pourRelaxDyn_ 1 #250 +nb_deb_test_amort_cinetique_ 1 +fi_parametre_amortissement_cinetique_ + +ARRET_A_EQUILIBRE_STATIQUE_ 2 + + + + +< tranche_sphere.her + + + +choix_materiaux +E_biel ELAS1D +E_2D ELAS2D + + +materiaux +ELAS1D ISOELAS1D + 344e6 0.3 + + +ELAS2D ISOELAS2D_C + 250e6 0.49 + + +masse_volumique +E_biel 1. +E_2D 1. + + +sections +E_biel 5e-6 + + +epaisseurs +E_2D 25e-6 + + +charges +F_to PRESSION 10 + + +blocages +N_to UX UY UZ + + +controle +TEMPSFIN 1. +DELTAt 0.5 +DELTAtMAXI 0.5 +ITERATIONS 100000 +NORME E_cin/E_stat_ET_min(Res,Res/Reac_et_Fext) +PRECISION 5.e-3 +SAUVEGARDE DERNIER_CALCUL + + +para_affichage +FREQUENCE_AFFICHAGE_INCREMENT 1 +FREQUENCE_AFFICHAGE_ITERATION 1 +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + + +para_calculs_geometriques +CAL_VOL_TOTAL_ENTRE_SURFACE_ET_PLANS_REF 1 + + +resultats pas_de_sortie_finale_ +COPIE 0 + + + +_suite_point_info_ + +controle +TEMPSFIN 2. +DELTAt 0.5 +DELTAtMAXI 0.5 + +_suite_point_info_ + +controle +TEMPSFIN 3. +DELTAt 0.5 +DELTAtMAXI 0.5 + + +_fin_point_info_ diff --git a/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.verif1 b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.verif1 new file mode 100755 index 0000000..03c3730 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info.verif1 @@ -0,0 +1,61 @@ +#!/usr/bin/env perl +use strict; +use warnings; +use English; +use File::Basename; +my $NOM_PROG = basename $PROGRAM_NAME; + + +#executable Herezh +my $exeHZ = shift(@ARGV); + +#fichier redirection affichage +my $flog = shift(@ARGV); + +#fichier .info +my $finfo = shift(@ARGV); + + +# +# +# verification des numeros d increment dans le fichier $flog (INCREMENT_DE_CHARGE) +# +# pour ce calcul, on s attend a trouver 6 increments de numeros : 1, 2, 3, 4, 5, 6 +# +# ce script renvoie OK si il y a 6 increments numerotes de 1 a 6 +# sinon ECHEC +# +# + +#liste des increments +my @liste_increments; + +open(FIC, "<$flog"); +while() { + next if(not /^\s*INCREMENT\s+DE\s+CHARGE\s+:\s+(\d+)\b/o); + push(@liste_increments, $1);#ajout du numero increment dans la liste +} +close(FIC); + +#nombre d increments +my $nb_incr = $#liste_increments + 1; + + +#ECHEC si le nombre d increments n est pas egal a 6 +if($nb_incr != 6) { + print "\nErreur (prog:$NOM_PROG) : $nb_incr increments ont ete calcules (on s attendait a 6 increments)\n\n"; + print "resultat verification : ECHEC $nb_incr increments ont ete calcules (on s attendait a 6 increments)\n"; + exit; +} + +#ECHEC si la liste n est pas dans l ordre exact 1 2 3 4 5 6 +for(1 .. 6) { + if($liste_increments[$_-1] != $_) { + print "\nErreur (prog:$NOM_PROG) : probleme dans les numeros des increments\n\n"; + print "resultat verification : ECHEC probleme dans les numeros des increments\n"; + exit; + } +} + +#OK : aucun probleme n a ete rencontre +print "resultat verification : OK\n"; diff --git a/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info_princ.maple.ref1 b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info_princ.maple.ref1 new file mode 100644 index 0000000..c83b739 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/sortie_suite_point_info_princ.maple.ref1 @@ -0,0 +1,104 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.769 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + # entete des donnees : informations gererales: on trouve successivement: + # >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs + # precedes du numero de colonne entre crochet + # >> le nombre de maillages m, et dimension de l'espace de travail + # puis pour chaque maillage, + # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit + # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions + # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # puis pour chaque maillage + # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , + # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees + # et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere) + # precedes du numero de colonne correspondant entre crochet + # puis pour chaque maillage + # >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant + # a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques, + # on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl + # chacun precedes du numero de colonne entre crochet + # puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres), + # suivi entre crochet, de la plage des numeros de colonnes, correspondant + # chacun sur une ligne differentes + # ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D: + # tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2) + # en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3) + # non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3) + # ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements, + # ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique + # ** ( cf: parametre_style_de_sortie = 0) + +#==================================================================== +#|| recapitulatif des differentes grandeurs par colonne || +#==================================================================== +#---------------------------------- grandeur globales ------------------------------------ +#6 (nombre de grandeurs globales) [2]energie_cinetique [3]energie_elastique [4]energie_interne [5]vol_total2D_avec_plan_xy [6]vol_total2D_avec_plan_xz [7]vol_total2D_avec_plan_yz +#---------------------------------- maillage et dimension -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#3 1 (nombre de noeuds, nombre total de grandeurs associees) +# noeud_1 [9]X [10]Y [11]Z [12]pression_ext ; +# noeud_23 [13]X [14]Y [15]Z [16]pression_ext ; +# noeud_44 [17]X [18]Y [19]Z [20]pression_ext ; +# +#---------------------------------- grandeurs aux elements ------------------------------ +#0 0 (nombre total d'elements, nombre totale de grandeurs associees) +#==================================================================== +#|| fin du recapitulatif des differentes grandeurs || +#==================================================================== + + # ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant + # a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun + # correspondant a un maillage. On trouve pour chaque enregistrement successivement : + # s'il y a des grandeurs globales: le temps puis les grandeurs globales, + # puis s'il y a des torseurs de reaction : + # de nouveau le temps, les composantes de la resultante puis les composantes du moments + # donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels + # puis s'il y a des grandeurs aux noeuds: de nouveau le temps + # les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud + # puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds + # puis s'il y a des grandeur aux elements: + # le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques) + # suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration + # correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes + # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs + # etc. pour tous les points d'integration - element + + 5.000000000000e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 5.000000000000e-01 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 +1.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 1.000000000000e+00 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 +1.500000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 1.500000000000e+00 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 +2.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 2.000000000000e+00 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 +2.500000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 2.500000000000e+00 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 +3.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.778094181861e+01 -2.198739281631e-01 -1.153420097091e+02 3.000000000000e+00 1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.099933241728e-16 9.961946980917e+00 0.000000000000e+00 0.000000000000e+00 -6.099933241728e-16 -9.961946980917e+00 0.000000000000e+00 diff --git a/Tests_en_attente_debug/Test_R_sortie_suite_point_info/tranche_sphere.her b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/tranche_sphere.her new file mode 100755 index 0000000..8adbf83 --- /dev/null +++ b/Tests_en_attente_debug/Test_R_sortie_suite_point_info/tranche_sphere.her @@ -0,0 +1,140 @@ +nom_maillage tranche_sphere + + noeuds ------------ +44 NOEUDS + +1 1.000000000000000e+1 0.000000000000000e+0 0.000000000000000e+0 +2 9.904184449275734e+0 8.45611829727469e-17 1.380988919117279e+0 +3 9.601347043876229e+0 1.71167500740341e-16 2.795377424078756e+0 +4 9.076542633423378e+0 2.57003545583734e-16 4.197186417547809e+0 +5 8.328684852428456e+0 3.38903210237086e-16 5.534709443948146e+0 +6 7.372773368101242e+0 4.13679692645907e-16 6.755902076156604e+0 +7 6.239540917964028e+0 4.78506854904277e-16 7.814609979586482e+0 +8 4.972287670860614e+0 5.31263706391684e-16 8.676194748748296e+0 +9 3.621552702019683e+0 5.70757529321800e-16 9.321177823992734e+0 +10 2.238940391334517e+0 5.96778646833825e-16 9.746134922319248e+0 +11 8.715574274765830e-1 6.09993324172810e-16 9.961946980917457e+0 +12 9.961946980917457e+0 8.715574274765816e-1 -5.33675006916149e-17 +13 9.866496037291203e+0 8.632065519864325e-1 1.380988919117279e+0 +14 9.564811019648355e+0 8.368125329870651e-1 2.795377424078756e+0 +15 9.042003648420060e+0 7.910728147967998e-1 4.197186417547810e+0 +16 8.296991692066261e+0 7.258927144245720e-1 5.534709443948144e+0 +17 7.344717739534476e+0 6.425795390070172e-1 6.755902076156602e+0 +18 6.215797581002266e+0 5.438118231095600e-1 7.814609979586482e+0 +19 4.953366615098298e+0 4.333634251088806e-1 8.676194748748298e+0 +20 3.607771600611843e+0 3.156391156443143e-1 9.321177823992736e+0 +21 2.230420547190903e+0 1.951365127744928e-1 9.746134922319250e+0 +22 8.682408883346531e-1 7.596123493896040e-2 9.961946980917457e+0 +23 0.000000000000000e+0 6.09993324172782e-16 9.961946980917000e+0 +24 9.904184449275734e+0 -8.45611829727469e-17 -1.380988919117279e+0 +25 9.601347043876229e+0 -1.71167500740341e-16 -2.795377424078756e+0 +26 9.076542633423378e+0 -2.57003545583734e-16 -4.197186417547809e+0 +27 8.328684852428456e+0 -3.38903210237086e-16 -5.534709443948146e+0 +28 7.372773368101242e+0 -4.13679692645907e-16 -6.755902076156604e+0 +29 6.239540917964028e+0 -4.78506854904277e-16 -7.814609979586482e+0 +30 4.972287670860614e+0 -5.31263706391684e-16 -8.676194748748296e+0 +31 3.621552702019683e+0 -5.70757529321800e-16 -9.321177823992734e+0 +32 2.238940391334517e+0 -5.96778646833825e-16 -9.746134922319248e+0 +33 8.715574274765830e-1 -6.09993324172810e-16 -9.961946980917457e+0 +34 9.866496037291203e+0 8.632065519864321e-1 -1.380988919117279e+0 +35 9.564811019648355e+0 8.368125329870644e-1 -2.795377424078756e+0 +36 9.042003648420060e+0 7.910728147967994e-1 -4.197186417547810e+0 +37 8.296991692066261e+0 7.258927144245712e-1 -5.534709443948144e+0 +38 7.344717739534476e+0 6.425795390070164e-1 -6.755902076156602e+0 +39 6.215797581002266e+0 5.438118231095592e-1 -7.814609979586482e+0 +40 4.953366615098298e+0 4.333634251088795e-1 -8.676194748748298e+0 +41 3.607771600611843e+0 3.156391156443132e-1 -9.321177823992736e+0 +42 2.230420547190903e+0 1.951365127744916e-1 -9.746134922319250e+0 +43 8.682408883346531e-1 7.596123493895918e-2 -9.961946980917457e+0 +44 0.000000000000000e+0 -6.09993324172782e-16 -9.961946980917000e+0 + + elements ------------ +44 ELEMENTS + +1 QUADRANGLE LINEAIRE 1 2 13 12 +2 QUADRANGLE LINEAIRE 2 3 14 13 +3 QUADRANGLE LINEAIRE 3 4 15 14 +4 QUADRANGLE LINEAIRE 4 5 16 15 +5 QUADRANGLE LINEAIRE 5 6 17 16 +6 QUADRANGLE LINEAIRE 6 7 18 17 +7 QUADRANGLE LINEAIRE 7 8 19 18 +8 QUADRANGLE LINEAIRE 8 9 20 19 +9 QUADRANGLE LINEAIRE 9 10 21 20 +10 QUADRANGLE LINEAIRE 10 11 22 21 +11 TRIANGLE LINEAIRE 11 23 22 +12 QUADRANGLE LINEAIRE 12 34 24 1 +13 QUADRANGLE LINEAIRE 34 35 25 24 +14 QUADRANGLE LINEAIRE 35 36 26 25 +15 QUADRANGLE LINEAIRE 36 37 27 26 +16 QUADRANGLE LINEAIRE 37 38 28 27 +17 QUADRANGLE LINEAIRE 38 39 29 28 +18 QUADRANGLE LINEAIRE 39 40 30 29 +19 QUADRANGLE LINEAIRE 40 41 31 30 +20 QUADRANGLE LINEAIRE 41 42 32 31 +21 QUADRANGLE LINEAIRE 42 43 33 32 +22 TRIANGLE LINEAIRE 33 43 44 +23 POUT BIE1 12 13 +24 POUT BIE1 13 14 +25 POUT BIE1 14 15 +26 POUT BIE1 15 16 +27 POUT BIE1 16 17 +28 POUT BIE1 17 18 +29 POUT BIE1 18 19 +30 POUT BIE1 19 20 +31 POUT BIE1 20 21 +32 POUT BIE1 21 22 +33 POUT BIE1 43 42 +34 POUT BIE1 42 41 +35 POUT BIE1 41 40 +36 POUT BIE1 40 39 +37 POUT BIE1 39 38 +38 POUT BIE1 38 37 +39 POUT BIE1 37 36 +40 POUT BIE1 36 35 +41 POUT BIE1 35 34 +42 POUT BIE1 34 12 +43 POUT BIE1 22 23 +44 POUT BIE1 44 43 + +N_y=0 2 3 4 5 6 7 8 9 10 11 + 24 25 26 27 28 29 30 31 32 33 + + +N_y=5 13 14 15 16 17 18 19 20 21 22 + 34 35 36 37 38 39 40 41 42 43 + + +N_pole 23 + +N_crochet 44 + +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 + +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 + +E_biel 23 24 25 26 27 28 29 30 31 32 + 33 34 35 36 37 38 39 40 41 42 + 43 44 + +E_quad 1 2 3 4 5 6 7 8 9 10 + 12 13 14 15 16 17 18 19 20 21 + + +E_tri 11 22 + +E_2D 1 2 3 4 5 6 7 8 9 10 + 11 12 13 14 15 16 17 18 19 20 + 21 22 + +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 +