2
0
Fork 0

ajout d un test de validation des sorties des parametres materiau pour la loi de Maxwell 3D, permettant d observer leurs evolutions lorsque les parametres changent pendant le calcul

This commit is contained in:
rio 2018-06-23 13:11:10 +00:00
parent 3cf0598827
commit 5a381a4500
5 changed files with 752 additions and 0 deletions

View file

@ -0,0 +1,34 @@
------------------------------------------------------
Auteur
------------------------------------------------------
Gérard Rio (gerard.rio@univ-ubs.fr)
------------------------------------------------------
Mots-cles
------------------------------------------------------
MAXWELL3D
E_YOUNG
NU_YOUNG
MU_VISCO
MU_VISCO_SPHERIQUE
les_courbes_1D
COURBEPOLYLINEAIRE_1_D
------------------------------------------------------
But du test
------------------------------------------------------
Test de sortie des paramètres matériau pour une loi de Maxwell isotrope 3D. Les paramètres matériau dépendent de la déformation duale de Mises et de $Q_eps$ c'est-à-dire de l'intensité du déviateur des déformations et du déviateur des vitesses de déformation. Aussi pendant le calcul, les paramètres matériau changent et l'objectif du test est de montrer la mise en données et de vérifier la bonne sortie (en .mapple mais également en format gmsh) des paramètres matériaux au cours du calcul.
------------------------------------------------------
Description du calcul
------------------------------------------------------
On considère un hexaèdre linéaire à 8 points d'intégration. Les conditions limites cinématiques sont isostatiques et on applique un déplacement uniforme perpendiculairement à une face. On regarde la sortie des valeurs du module d'Young, du coefficient de viscosité $\mu$ , qui varient pendant le calcul. On sort également $\nu$ et la viscosité sphérique, qui eux sont fixes.
-------------------------------------------------------------
Grandeurs de comparaison
-------------------------------------------------------------
On regarde les paramètres matériaux. L'idée est plus ici de servir de référence.

View file

@ -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 8 x 7 ( 1 x 1 x 1 elements )
# geometrie prismatique, decoupage en hexaedre(s), interpolation lineaire.
nom_maillage hexaedre
noeuds ------------
8 NOEUDS
#---------------------------------------------------------------
#|NO DU| X | Y | Z |
#|NOEUD| | | |
#---------------------------------------------------------------
1 0 0 0
2 0 0 7
3 0 8 0
4 0 8 7
5 10 0 0
6 10 0 7
7 10 8 0
8 10 8 7
# 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

View file

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

View file

@ -0,0 +1,153 @@
#######################################################################
# #
# | | ==== === ==== ==== | | | | #
# | | | | | | / | | | | #
# |====| |=== === |=== / |====| ------- ------- #
# | | | | \ | / | | | | #
# | | ==== | \ ==== ==== | | | | #
# #
#######################################################################
# 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.847
#######################################################################
# fichier de commande #
# --- l'objectif est de tester les sorties de paramètres matériaux qui évoluent
# pendant le calcul
#----------------------------------------------------
# definition de la dimension de l'espace de travail |
#----------------------------------------------------
dimension 3
#---------------------------------------------------------------
# definition facultative du niveau d'impression (entre 0 et 10)|
#---------------------------------------------------------------
niveau_commentaire 3
#-------------------------------
# definition du type de calcul |
#-------------------------------
TYPE_DE_CALCUL
non_dynamique avec # plus visualisation
#--------------------------------------
#| definition du (ou des) maillage(s) |
#--------------------------------------
# -- def maillage
< hexaedre.her
#----------------------------------------
# definition facultative de courbes
#----------------------------------------
les_courbes_1D
fct_Mises_module COURBEPOLYLINEAIRE_1_D
Debut_des_coordonnees_des_points
Coordonnee dim= 2 0 250
Coordonnee dim= 2 0.5 350
Fin_des_coordonnees_des_points
fct_Mises COURBEPOLYLINEAIRE_1_D
Debut_des_coordonnees_des_points
Coordonnee dim= 2 0 0.06
Coordonnee dim= 2 0.5 0.09
Fin_des_coordonnees_des_points
fct_Q_Deps COURBEPOLYLINEAIRE_1_D
Debut_des_coordonnees_des_points
Coordonnee dim= 2 0. 1000000
Coordonnee dim= 2 1. 100
Fin_des_coordonnees_des_points
#------------- fin def des fonctions nD ------------------------
#======================================
#| definition des lois de comportement|
#--------------------------------------
choix_materiaux ----------
#----------------------------------------
# Elements | Nom Materiau |
#----------------------------------------
E_tout lune
materiaux #----------
#------------------------------------------
# Nom Materiau | Type loi |
#------------------------------------------
lune MAXWELL3D
E= 1. nu= 0.3 mu= 0.4 \
fac_mu_cissionD= fct_Q_Deps # mu dépend de la vitesse de déformation
fac_mu_Mises_Eps= fct_Mises # et de la déformation duale de Mises
fac_E_Mises_Eps= fct_Mises_module # E dépend de la déformation duale de Mises
fin_coeff_MAXWELL3D
#--------------------------------------- fin def des lois de comportement ------------------------
# --- divers stockages (1) -------
masse_volumique #--------#
E_tout 1
charges #------------#
blocages #------------#
#-----------------------------------------------------------
# nom du maillage | Ref noeud | Blocages
#-----------------------------------------------------------
N_droit UX
N_arriere UY
N_bas UZ
N_gauche ' UX = 10 '
controle #------------
#---------------------------
# PARAMETRE | VALEUR |
#---------------------------
SAUVEGARDE 1
#DELTAtMAXI 2.
TEMPSFIN 2.
DELTAt 0.001
NORME Residu/Reaction_et_VarRes
ITERATIONS 16
PRECISION 0.001
MULTIPLICATEUR 1
#RESTART 4
DELTAtMAXI 0.1 #0.01
#MAXINCRE 2
para_affichage #------------
#---------------------------
# PARAMETRE | VALEUR |
#---------------------------
FREQUENCE_AFFICHAGE_INCREMENT 1
FREQUENCE_AFFICHAGE_ITERATION 1#0
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
# ------------------------------------------------------------------------------------
resultats pas_de_sortie_finale_
COPIE 0
#
_fin_point_info_

View file

@ -0,0 +1,211 @@
#fichier au format maple6
###############################################################################################
# Visualisation elements finis : Herezh++ V6.848 #
# 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 --------------------------------
#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 1 41 (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] def_duale_mises
# [6]X [7]Y [8]Z
# E_YOUNG E_YOUNG SCALAIRE TABLEAU_T SCALAIRE 1 [9...9]
# NU_YOUNG NU_YOUNG SCALAIRE TABLEAU_T SCALAIRE 1 [10...10]
# MU_VISCO MU_VISCO SCALAIRE TABLEAU_T SCALAIRE 1 [11...11]
# MU_VISCO_SPHERIQUE MU_VISCO_SPHERIQUE SCALAIRE TABLEAU_T SCALAIRE 1 [12...12]
#
# [13]X [14]Y [15]Z
# Q_DEPS Q_DEPS SCALAIRE TYPE_SIMPLE SCALAIRE 1 [16...16]
#
#
#
#====================================================================
#|| 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-03 2.115361902706e+00 6.307511169766e+00 5.519072273546e+00 8.654499277207e-04 2.115361902706e+00 6.307511169766e+00 5.519072273546e+00 2.501730899855e+02 3.000000000000e-01 -1.437627876352e+03 3.594069690880e-13 2.115361902706e+00 6.307511169766e+00 5.519072273546e+00 1.059955360422e+00
2.000000000000e-03 2.117475151360e+00 6.305631371807e+00 5.517427450331e+00 1.728017429224e-03 2.117475151360e+00 6.305631371807e+00 5.517427450331e+00 2.503456034858e+02 3.000000000000e-01 -1.389594832319e+03 3.473987080797e-13 2.117475151360e+00 6.305631371807e+00 5.517427450331e+00 1.057905696006e+00
3.000000000000e-03 2.119588400014e+00 6.303758891829e+00 5.515789030350e+00 2.588003489418e-03 2.119588400014e+00 6.303758891829e+00 5.515789030350e+00 2.505176006979e+02 3.000000000000e-01 -1.350188470523e+03 3.375471176306e-13 2.119588400014e+00 6.303758891829e+00 5.515789030350e+00 1.056218255073e+00
4.000000000000e-03 2.121701648668e+00 6.301886411850e+00 5.514150610368e+00 3.446190151155e-03 2.121701648668e+00 6.301886411850e+00 5.514150610368e+00 2.506892380302e+02 3.000000000000e-01 -1.333574912338e+03 3.333937280845e-13 2.121701648668e+00 6.301886411850e+00 5.514150610368e+00 1.055480337328e+00
5.414213562373e-03 2.124690233575e+00 6.299271972209e+00 5.511862975683e+00 4.653220851858e-03 2.124690233575e+00 6.299271972209e+00 5.511862975683e+00 2.509306441704e+02 3.000000000000e-01 -1.232002498825e+03 3.080006247063e-13 2.124690233575e+00 6.299271972209e+00 5.511862975683e+00 1.051200798057e+00
6.828427124746e-03 2.127678818483e+00 6.296657532569e+00 5.509575340998e+00 5.856682629746e-03 2.127678818483e+00 6.296657532569e+00 5.509575340998e+00 2.511713365259e+02 3.000000000000e-01 -1.208385374624e+03 3.020963436561e-13 2.127678818483e+00 6.296657532569e+00 5.509575340998e+00 1.050161243297e+00
8.242640687119e-03 2.130667403390e+00 6.294075382357e+00 5.507315959562e+00 7.053161217608e-03 2.130667403390e+00 6.294075382357e+00 5.507315959562e+00 2.514106322435e+02 3.000000000000e-01 -1.113140475021e+03 2.782851187553e-13 2.130667403390e+00 6.294075382357e+00 5.507315959562e+00 1.046160628707e+00
9.656854249492e-03 2.133655988297e+00 6.291493232145e+00 5.505056578127e+00 8.246107254530e-03 2.133655988297e+00 6.291493232145e+00 5.505056578127e+00 2.516492214509e+02 3.000000000000e-01 -1.089405496164e+03 2.723513740411e-13 2.133655988297e+00 6.291493232145e+00 5.505056578127e+00 1.045125163077e+00
1.165685424949e-02 2.137882485605e+00 6.287888312806e+00 5.501902273705e+00 9.922205185389e-03 2.137882485605e+00 6.287888312806e+00 5.501902273705e+00 2.519844410371e+02 3.000000000000e-01 -9.769716043334e+02 2.442429010833e-13 2.137882485605e+00 6.287888312806e+00 5.501902273705e+00 1.040411254861e+00
1.365685424949e-02 2.142108982913e+00 6.284283393467e+00 5.498747969283e+00 1.159132734220e-02 2.142108982913e+00 6.284283393467e+00 5.498747969283e+00 2.523182654684e+02 3.000000000000e-01 -9.432663726435e+02 2.358165931609e-13 2.142108982913e+00 6.284283393467e+00 5.498747969283e+00 1.038956310104e+00
1.565685424949e-02 2.146335480221e+00 6.280758644914e+00 5.495663814300e+00 1.324494771832e-02 2.146335480221e+00 6.280758644914e+00 5.495663814300e+00 2.526489895437e+02 3.000000000000e-01 -7.827318428414e+02 1.956829607104e-13 2.146335480221e+00 6.280758644914e+00 5.495663814300e+00 1.032290734038e+00
1.765685424949e-02 2.150561977529e+00 6.277233896362e+00 5.492579659317e+00 1.489168569650e-02 2.150561977529e+00 6.277233896362e+00 5.492579659317e+00 2.529783371393e+02 3.000000000000e-01 -7.487356255167e+02 1.871839063792e-13 2.150561977529e+00 6.277233896362e+00 5.492579659317e+00 1.030842638218e+00
2.048528137424e-02 2.156539147344e+00 6.272323608596e+00 5.488283157521e+00 1.720088497856e-02 2.156539147344e+00 6.272323608596e+00 5.488283157521e+00 2.534401769957e+02 3.000000000000e-01 -6.096520554589e+02 1.524130138647e-13 2.156539147344e+00 6.272323608596e+00 5.488283157521e+00 1.025075125363e+00
2.331370849898e-02 2.162516317158e+00 6.267471992189e+00 5.484037993166e+00 1.949024618284e-02 2.162516317158e+00 6.267471992189e+00 5.484037993166e+00 2.538980492366e+02 3.000000000000e-01 -4.951869226518e+02 1.237967306630e-13 2.162516317158e+00 6.267471992189e+00 5.484037993166e+00 1.020340372286e+00
2.614213562373e-02 2.168493486973e+00 6.262691631515e+00 5.479855177575e+00 2.175853183615e-02 2.168493486973e+00 6.262691631515e+00 5.479855177575e+00 2.543517063672e+02 3.000000000000e-01 -3.661001302914e+02 9.152503257286e-14 2.168493486973e+00 6.262691631515e+00 5.479855177575e+00 1.015030834839e+00
2.897056274848e-02 2.174470656787e+00 6.257998411207e+00 5.475748609806e+00 2.400413283208e-02 2.174470656787e+00 6.257998411207e+00 5.475748609806e+00 2.548008265664e+02 3.000000000000e-01 -2.185157288982e+02 5.462893222454e-14 2.174470656787e+00 6.257998411207e+00 5.475748609806e+00 1.008992291115e+00
3.030389608181e-02 2.177288321659e+00 6.255786005936e+00 5.473812755194e+00 2.505815601478e-02 2.177288321659e+00 6.255786005936e+00 5.473812755194e+00 2.550116312030e+02 3.000000000000e-01 -2.085222850034e+02 5.213057125086e-14 2.177288321659e+00 6.255786005936e+00 5.473812755194e+00 1.008576892430e+00
3.107369644073e-02 2.178915101232e+00 6.254508673158e+00 5.472695089013e+00 2.566537461438e-02 2.178915101232e+00 6.254508673158e+00 5.472695089013e+00 2.551330749229e+02 3.000000000000e-01 -2.001185653825e+02 5.002964134562e-14 2.178915101232e+00 6.254508673158e+00 5.472695089013e+00 1.008230446764e+00
3.184349679965e-02 2.180541880804e+00 6.253310480397e+00 5.471646670349e+00 2.626304012509e-02 2.180541880804e+00 6.253310480397e+00 5.471646670349e+00 2.552526080250e+02 3.000000000000e-01 1.439361421100e+02 -3.598403552749e-14 2.180541880804e+00 6.253310480397e+00 5.471646670349e+00 9.942555638756e-01
3.246935074999e-02 2.181864465822e+00 6.252251429005e+00 5.470720000378e+00 2.675745143969e-02 2.181864465822e+00 6.252251429005e+00 5.470720000378e+00 2.553514902879e+02 3.000000000000e-01 -2.829328373940e+02 7.073320934849e-14 2.181864465822e+00 6.252251429005e+00 5.470720000378e+00 1.011582806851e+00
3.309520470033e-02 2.183187050840e+00 6.251192377613e+00 5.469793330408e+00 2.725123179806e-02 2.183187050840e+00 6.251192377613e+00 5.469793330408e+00 2.554502463596e+02 3.000000000000e-01 -2.721890084769e+02 6.804725211924e-14 2.183187050840e+00 6.251192377613e+00 5.469793330408e+00 1.011141459814e+00
3.372105865067e-02 2.184509635859e+00 6.250246819968e+00 5.468965967471e+00 2.773204682578e-02 2.184509635859e+00 6.250246819968e+00 5.468965967471e+00 2.555464093652e+02 3.000000000000e-01 3.230035946355e+02 -8.075089865888e-14 2.184509635859e+00 6.250246819968e+00 5.468965967471e+00 9.870033773385e-01
3.434691260101e-02 2.185832220877e+00 6.249094569510e+00 5.467957748323e+00 2.823469939506e-02 2.185832220877e+00 6.249094569510e+00 5.467957748323e+00 2.556469398790e+02 3.000000000000e-01 -7.308846395649e+02 1.827211598912e-13 2.185832220877e+00 6.249094569510e+00 5.467957748323e+00 1.029720263357e+00
3.523200374565e-02 2.187702638547e+00 6.247739613356e+00 5.466772161688e+00 2.891462643935e-02 2.187702638547e+00 6.247739613356e+00 5.466772161688e+00 2.557829252879e+02 3.000000000000e-01 2.874780549880e+02 -7.186951374700e-14 2.187702638547e+00 6.247739613356e+00 5.466772161688e+00 9.884572070625e-01
3.611709489029e-02 2.189573056217e+00 6.246278333692e+00 5.465493541980e+00 2.960486343769e-02 2.189573056217e+00 6.246278333692e+00 5.465493541980e+00 2.559209726875e+02 3.000000000000e-01 -8.492357565275e+01 2.123089391319e-14 2.189573056217e+00 6.246278333692e+00 5.465493541980e+00 1.003537091864e+00
3.700218603493e-02 2.191443473887e+00 6.244875081306e+00 5.464265696141e+00 3.028752922478e-02 2.191443473887e+00 6.244875081306e+00 5.464265696141e+00 2.560575058450e+02 3.000000000000e-01 1.424190009143e+02 -3.560475022858e-14 2.191443473887e+00 6.244875081306e+00 5.464265696141e+00 9.943397553715e-01
3.751319364554e-02 2.192523360032e+00 6.244064913163e+00 5.463556799016e+00 3.068109622075e-02 2.192523360032e+00 6.244064913163e+00 5.463556799016e+00 2.561362192442e+02 3.000000000000e-01 1.481260497820e+02 -3.703151244551e-14 2.192523360032e+00 6.244064913163e+00 5.463556799016e+00 9.941112167485e-01
3.802420125615e-02 2.193603246178e+00 6.243116882103e+00 5.462727271835e+00 3.108928222695e-02 2.193603246178e+00 6.243116882103e+00 5.462727271835e+00 2.562178564454e+02 3.000000000000e-01 -7.161817725628e+02 1.790454431407e-13 2.193603246178e+00 6.243116882103e+00 5.462727271835e+00 1.029044052087e+00
3.853520886677e-02 2.194683132323e+00 6.242372941452e+00 5.462076323768e+00 3.147479757243e-02 2.194683132323e+00 6.242372941452e+00 5.462076323768e+00 2.562949595145e+02 3.000000000000e-01 5.856681995566e+02 -1.464170498891e-13 2.194683132323e+00 6.242372941452e+00 5.462076323768e+00 9.764394394918e-01
3.904621647738e-02 2.195763018469e+00 6.241508162810e+00 5.461319642457e+00 3.187308184541e-02 2.195763018469e+00 6.241508162810e+00 5.461319642457e+00 2.563746163691e+02 3.000000000000e-01 -1.712989870819e+02 4.282474677047e-14 2.195763018469e+00 6.241508162810e+00 5.461319642457e+00 1.007017693714e+00
3.976889037078e-02 2.197290208101e+00 6.240442234712e+00 5.460386955373e+00 3.241849161019e-02 2.197290208101e+00 6.240442234712e+00 5.460386955373e+00 2.564836983220e+02 3.000000000000e-01 5.476294844726e+02 -1.369073711181e-13 2.197290208101e+00 6.240442234712e+00 5.460386955373e+00 9.779963985396e-01
4.049156426419e-02 2.198817397734e+00 6.239197637719e+00 5.459297933002e+00 3.298258958771e-02 2.198817397734e+00 6.239197637719e+00 5.459297933002e+00 2.565965179175e+02 3.000000000000e-01 -2.416979575432e+02 6.042448938579e-14 2.198817397734e+00 6.239197637719e+00 5.459297933002e+00 1.009850179566e+00
4.121423815760e-02 2.200344587366e+00 6.238088873608e+00 5.458327764404e+00 3.353102699439e-02 2.200344587366e+00 6.238088873608e+00 5.458327764404e+00 2.567062053989e+02 3.000000000000e-01 3.809724671409e+02 -9.524311678523e-14 2.200344587366e+00 6.238088873608e+00 5.458327764404e+00 9.847396193335e-01
4.193691205100e-02 2.201871776999e+00 6.236980109497e+00 5.457357595807e+00 3.407864033867e-02 2.201871776999e+00 6.236980109497e+00 5.457357595807e+00 2.568157280677e+02 3.000000000000e-01 3.937758796184e+02 -9.844396990461e-14 2.201871776999e+00 6.236980109497e+00 5.457357595807e+00 9.842318075508e-01
4.295892727223e-02 2.204031549289e+00 6.235342890033e+00 5.455925028829e+00 3.485925494608e-02 2.204031549289e+00 6.235342890033e+00 5.455925028829e+00 2.569718509892e+02 3.000000000000e-01 1.942138497340e+02 -4.855346243350e-14 2.204031549289e+00 6.235342890033e+00 5.455925028829e+00 9.922795716124e-01
4.354898803532e-02 2.205278494403e+00 6.234342207092e+00 5.455049431224e+00 3.531527422262e-02 2.205278494403e+00 6.234342207092e+00 5.455049431224e+00 2.570630548445e+02 3.000000000000e-01 -1.050228729759e+02 2.625571824397e-14 2.205278494403e+00 6.234342207092e+00 5.455049431224e+00 1.004327119162e+00
4.413904879841e-02 2.206525439516e+00 6.233490504346e+00 5.454304191104e+00 3.575443478099e-02 2.206525439516e+00 6.233490504346e+00 5.454304191104e+00 2.571508869562e+02 3.000000000000e-01 7.272981796666e+02 -1.818245449166e-13 2.206525439516e+00 6.233490504346e+00 5.454304191104e+00 9.708390959138e-01
4.472910956150e-02 2.207772384630e+00 6.232485709310e+00 5.453424996493e+00 3.620982945399e-02 2.207772384630e+00 6.232485709310e+00 5.453424996493e+00 2.572419658908e+02 3.000000000000e-01 -1.074967043120e+02 2.687417607799e-14 2.207772384630e+00 6.232485709310e+00 5.453424996493e+00 1.004422954222e+00
4.531917032459e-02 2.209019329743e+00 6.231678272266e+00 5.452718494401e+00 3.664304821324e-02 2.209019329743e+00 6.231678272266e+00 5.452718494401e+00 2.573286096426e+02 3.000000000000e-01 9.930752400434e+02 -2.482688100108e-13 2.209019329743e+00 6.231678272266e+00 5.452718494401e+00 9.601805089727e-01
4.615364225838e-02 2.210782776434e+00 6.230318171178e+00 5.451528401325e+00 3.727872109724e-02 2.210782776434e+00 6.230318171178e+00 5.451528401325e+00 2.574557442194e+02 3.000000000000e-01 1.571089817752e+02 -3.927724544381e-14 2.210782776434e+00 6.230318171178e+00 5.451528401325e+00 9.937884351783e-01
4.698811419217e-02 2.212546223125e+00 6.228958070089e+00 5.450338308248e+00 3.791332901791e-02 2.212546223125e+00 6.228958070089e+00 5.450338308248e+00 2.575826658036e+02 3.000000000000e-01 1.715245986321e+02 -4.288114965801e-14 2.212546223125e+00 6.228958070089e+00 5.450338308248e+00 9.932135264828e-01
4.782258612596e-02 2.214309669816e+00 6.227597969001e+00 5.449148215172e+00 3.854687576304e-02 2.214309669816e+00 6.227597969001e+00 5.449148215172e+00 2.577093751526e+02 3.000000000000e-01 1.859323645831e+02 -4.648309114576e-14 2.214309669816e+00 6.227597969001e+00 5.449148215172e+00 9.926396277745e-01
4.865705805975e-02 2.216073116507e+00 6.226237867913e+00 5.447958122095e+00 3.917936510580e-02 2.216073116507e+00 6.226237867913e+00 5.447958122095e+00 2.578358730212e+02 3.000000000000e-01 2.003322490024e+02 -5.008306225060e-14 2.216073116507e+00 6.226237867913e+00 5.447958122095e+00 9.920667367055e-01
4.983717958593e-02 2.218567006734e+00 6.224446219615e+00 5.446390433529e+00 4.005753409564e-02 2.218567006734e+00 6.224446219615e+00 5.446390433529e+00 2.580115068191e+02 3.000000000000e-01 5.895938922058e+02 -1.473984730514e-13 2.218567006734e+00 6.224446219615e+00 5.446390433529e+00 9.764774042732e-01
5.101730111211e-02 2.221060896960e+00 6.222654571317e+00 5.444822744962e+00 4.093358038312e-02 2.221060896960e+00 6.222654571317e+00 5.444822744962e+00 2.581867160766e+02 3.000000000000e-01 6.104474906004e+02 -1.526118726501e-13 2.221060896960e+00 6.222654571317e+00 5.444822744962e+00 9.756624701105e-01
5.219742263830e-02 2.223554787187e+00 6.220862923020e+00 5.443255056395e+00 4.180751444558e-02 2.223554787187e+00 6.220862923020e+00 5.443255056395e+00 2.583615028891e+02 3.000000000000e-01 6.312843562716e+02 -1.578210890679e-13 2.223554787187e+00 6.220862923020e+00 5.443255056395e+00 9.748495245760e-01
5.337754416448e-02 2.226048677414e+00 6.219071274722e+00 5.441687367828e+00 4.267934670307e-02 2.226048677414e+00 6.219071274722e+00 5.441687367828e+00 2.585358693406e+02 3.000000000000e-01 6.521044106737e+02 -1.630261026684e-13 2.226048677414e+00 6.219071274722e+00 5.441687367828e+00 9.740385611282e-01
5.504648803205e-02 2.229575570796e+00 6.216537501401e+00 5.439470321396e+00 4.390873582250e-02 2.229575570796e+00 6.216537501401e+00 5.439470321396e+00 2.587817471645e+02 3.000000000000e-01 6.857125999614e+02 -1.714281499904e-13 2.229575570796e+00 6.216537501401e+00 5.439470321396e+00 9.727276787112e-01
5.671543189963e-02 2.233102464178e+00 6.214003728079e+00 5.437253274963e+00 4.513397116774e-02 2.233102464178e+00 6.214003728079e+00 5.437253274963e+00 2.590267942335e+02 3.000000000000e-01 7.150842483962e+02 -1.787710620990e-13 2.233102464178e+00 6.214003728079e+00 5.437253274963e+00 9.715886825978e-01
5.838437576721e-02 2.236629357560e+00 6.211469954757e+00 5.435036228531e+00 4.635508163373e-02 2.236629357560e+00 6.211469954757e+00 5.435036228531e+00 2.592710163267e+02 3.000000000000e-01 7.444217111524e+02 -1.861054277881e-13 2.236629357560e+00 6.211469954757e+00 5.435036228531e+00 9.704535984627e-01
6.005331963478e-02 2.240156250941e+00 6.208936181436e+00 5.432819182098e+00 4.757209589366e-02 2.240156250941e+00 6.208936181436e+00 5.432819182098e+00 2.595144191787e+02 3.000000000000e-01 7.737247826275e+02 -1.934311956569e-13 2.240156250941e+00 6.208936181436e+00 5.432819182098e+00 9.693224082383e-01
6.241356268715e-02 2.245144031395e+00 6.205352884841e+00 5.429683804964e+00 4.928627654054e-02 2.245144031395e+00 6.205352884841e+00 5.429683804964e+00 2.598572553081e+02 3.000000000000e-01 8.209771151743e+02 -2.052442787936e-13 2.245144031395e+00 6.205352884841e+00 5.429683804964e+00 9.674961328286e-01
6.477380573951e-02 2.250131811849e+00 6.202021535765e+00 5.426768834634e+00 5.096437447593e-02 2.250131811849e+00 6.202021535765e+00 5.426768834634e+00 2.601928748952e+02 3.000000000000e-01 1.217041473646e+03 -3.042603684114e-13 2.250131811849e+00 6.202021535765e+00 5.426768834634e+00 9.518442040509e-01
6.713404879188e-02 2.255119592302e+00 6.198690186690e+00 5.423853864304e+00 5.263440163762e-02 2.255119592302e+00 6.198690186690e+00 5.423853864304e+00 2.605268803275e+02 3.000000000000e-01 1.259199566584e+03 -3.147998916460e-13 2.255119592302e+00 6.198690186690e+00 5.423853864304e+00 9.502518426642e-01
6.949429184425e-02 2.260107372756e+00 6.195208636800e+00 5.420807551260e+00 5.431319429663e-02 2.260107372756e+00 6.195208636800e+00 5.420807551260e+00 2.608626388593e+02 3.000000000000e-01 1.088937889682e+03 -2.722344724205e-13 2.260107372756e+00 6.195208636800e+00 5.420807551260e+00 9.570606645669e-01
7.283217957940e-02 2.267161159520e+00 6.190497742595e+00 5.416685515817e+00 5.665016596485e-02 2.267161159520e+00 6.190497742595e+00 5.416685515817e+00 2.613300331930e+02 3.000000000000e-01 1.369224088567e+03 -3.423060221418e-13 2.267161159520e+00 6.190497742595e+00 5.416685515817e+00 9.461022809521e-01
7.617006731455e-02 2.274214946283e+00 6.185786848389e+00 5.412563480375e+00 5.897151970277e-02 2.274214946283e+00 6.185786848389e+00 5.412563480375e+00 2.617943039406e+02 3.000000000000e-01 1.428428551580e+03 -3.571071378950e-13 2.274214946283e+00 6.185786848389e+00 5.412563480375e+00 9.438909358293e-01
7.950795504970e-02 2.281268733047e+00 6.181240670725e+00 5.408585580026e+00 6.125895611753e-02 2.281268733047e+00 6.181240670725e+00 5.408585580026e+00 2.622517912235e+02 3.000000000000e-01 1.653721566737e+03 -4.134303916842e-13 2.281268733047e+00 6.181240670725e+00 5.408585580026e+00 9.351658522451e-01
8.284584278485e-02 2.288322519811e+00 6.176694493060e+00 5.404607679678e+00 6.353111091418e-02 2.288322519811e+00 6.176694493060e+00 5.404607679678e+00 2.627062221828e+02 3.000000000000e-01 1.713201170316e+03 -4.283002925790e-13 2.288322519811e+00 6.176694493060e+00 5.404607679678e+00 9.329740741295e-01
8.756632888959e-02 2.298298080718e+00 6.170265226949e+00 5.398982079055e+00 6.671873278677e-02 2.298298080718e+00 6.170265226949e+00 5.398982079055e+00 2.633437465574e+02 3.000000000000e-01 1.808563645018e+03 -4.521409112544e-13 2.298298080718e+00 6.170265226949e+00 5.398982079055e+00 9.294493758071e-01
9.228681499432e-02 2.308273641625e+00 6.164226507448e+00 5.393698186534e+00 6.983251511297e-02 2.308273641625e+00 6.164226507448e+00 5.393698186534e+00 2.639665030226e+02 3.000000000000e-01 2.173867385945e+03 -5.434668464862e-13 2.308273641625e+00 6.164226507448e+00 5.393698186534e+00 9.154261333384e-01
9.700730109905e-02 2.318249202533e+00 6.158187787946e+00 5.388414294013e+00 7.291699923429e-02 2.318249202533e+00 6.158187787946e+00 5.388414294013e+00 2.645833998469e+02 3.000000000000e-01 2.258314208263e+03 -5.645785520657e-13 2.318249202533e+00 6.158187787946e+00 5.388414294013e+00 9.123897435078e-01
1.017277872038e-01 2.328224763440e+00 6.152149068444e+00 5.383130401492e+00 7.597272848861e-02 2.328224763440e+00 6.152149068444e+00 5.383130401492e+00 2.651945456977e+02 3.000000000000e-01 2.342432363026e+03 -5.856080907565e-13 2.328224763440e+00 6.152149068444e+00 5.383130401492e+00 9.093810661227e-01
1.084035626741e-01 2.342332336967e+00 6.143974460719e+00 5.375977647467e+00 8.020430316672e-02 2.342332336967e+00 6.143974460719e+00 5.375977647467e+00 2.660408606333e+02 3.000000000000e-01 2.665785605643e+03 -6.664464014108e-13 2.342332336967e+00 6.143974460719e+00 5.375977647467e+00 8.972625111498e-01
1.150793381444e-01 2.356439910495e+00 6.135799852993e+00 5.368824893442e+00 8.438055502180e-02 2.356439910495e+00 6.135799852993e+00 5.368824893442e+00 2.668761110044e+02 3.000000000000e-01 2.784516600516e+03 -6.961291501291e-13 2.356439910495e+00 6.135799852993e+00 5.368824893442e+00 8.930959441341e-01
1.217551136147e-01 2.370547484022e+00 6.128170709606e+00 5.362149362276e+00 8.843999989595e-02 2.370547484022e+00 6.128170709606e+00 5.362149362276e+00 2.676879999792e+02 3.000000000000e-01 3.187481015194e+03 -7.968702537986e-13 2.370547484022e+00 6.128170709606e+00 5.362149362276e+00 8.780675677182e-01
1.284308890850e-01 2.384655057549e+00 6.120541566218e+00 5.355473831109e+00 9.244643592926e-02 2.384655057549e+00 6.120541566218e+00 5.355473831109e+00 2.684892871859e+02 3.000000000000e-01 3.306450925747e+03 -8.266127314368e-13 2.384655057549e+00 6.120541566218e+00 5.355473831109e+00 8.739770603552e-01
1.378718612945e-01 2.404606179364e+00 6.110156296440e+00 5.346386755076e+00 9.797748236318e-02 2.404606179364e+00 6.110156296440e+00 5.346386755076e+00 2.695954964726e+02 3.000000000000e-01 3.646477185730e+03 -9.116192964324e-13 2.404606179364e+00 6.110156296440e+00 5.346386755076e+00 8.617076048555e-01
1.473128335039e-01 2.424557301179e+00 6.100124273018e+00 5.337608730753e+00 1.033670382968e-01 2.424557301179e+00 6.100124273018e+00 5.337608730753e+00 2.706734076594e+02 3.000000000000e-01 3.945837585092e+03 -9.864593962729e-13 2.424557301179e+00 6.100124273018e+00 5.337608730753e+00 8.510776662967e-01
1.567538057134e-01 2.444508422993e+00 6.090461891094e+00 5.329154148117e+00 1.086160612780e-01 2.444508422993e+00 6.090461891094e+00 5.329154148117e+00 2.717232122556e+02 3.000000000000e-01 4.251723516371e+03 -1.062930879093e-12 2.444508422993e+00 6.090461891094e+00 5.329154148117e+00 8.402855620689e-01
1.661947779228e-01 2.464459544808e+00 6.081148185984e+00 5.321004655923e+00 1.137297604068e-01 2.464459544808e+00 6.081148185984e+00 5.321004655923e+00 2.727459520814e+02 3.000000000000e-01 4.549931472306e+03 -1.137482868076e-12 2.464459544808e+00 6.081148185984e+00 5.321004655923e+00 8.298617307552e-01
1.795463288635e-01 2.492674691863e+00 6.068556042551e+00 5.309986530549e+00 1.207376151845e-01 2.492674691863e+00 6.068556042551e+00 5.309986530549e+00 2.741475230369e+02 3.000000000000e-01 4.969682989217e+03 -1.242420747304e-12 2.492674691863e+00 6.068556042551e+00 5.309986530549e+00 8.153191738550e-01
1.928978798041e-01 2.520889838917e+00 6.056590694567e+00 5.299516851086e+00 1.274956779188e-01 2.520889838917e+00 6.056590694567e+00 5.299516851086e+00 2.754991355838e+02 3.000000000000e-01 5.371708521606e+03 -1.342927130402e-12 2.520889838917e+00 6.056590694567e+00 5.299516851086e+00 8.015683640343e-01
2.062494307447e-01 2.549104985972e+00 6.045223089286e+00 5.289570196512e+00 1.340142902389e-01 2.549104985972e+00 6.045223089286e+00 5.289570196512e+00 2.768028580478e+02 3.000000000000e-01 5.765321137489e+03 -1.441330284372e-12 2.549104985972e+00 6.045223089286e+00 5.289570196512e+00 7.882457673307e-01
2.196009816853e-01 2.577320133027e+00 6.034426266043e+00 5.280122976213e+00 1.403032247417e-01 2.577320133027e+00 6.034426266043e+00 5.280122976213e+00 2.780606449483e+02 3.000000000000e-01 6.150775971551e+03 -1.537693992888e-12 2.577320133027e+00 6.034426266043e+00 5.280122976213e+00 7.753282498384e-01
2.384829261042e-01 2.617222376656e+00 6.020089680750e+00 5.267578464134e+00 1.488221165338e-01 2.617222376656e+00 6.020089680750e+00 5.267578464134e+00 2.797644233068e+02 3.000000000000e-01 6.692424676861e+03 -1.673106169215e-12 2.617222376656e+00 6.020089680750e+00 5.267578464134e+00 7.573479667730e-01
2.573648705231e-01 2.657124620285e+00 6.006735383635e+00 5.255893454207e+00 1.569299521748e-01 2.657124620285e+00 6.006735383635e+00 5.255893454207e+00 2.813859904350e+02 3.000000000000e-01 7.200775038875e+03 -1.800193759719e-12 2.657124620285e+00 6.006735383635e+00 5.255893454207e+00 7.407391913513e-01
2.762468149421e-01 2.697026863914e+00 5.994305042131e+00 5.245016905437e+00 1.646494300915e-01 2.697026863914e+00 5.994305042131e+00 5.245016905437e+00 2.829298860183e+02 3.000000000000e-01 7.694437219292e+03 -1.923609304823e-12 2.697026863914e+00 5.994305042131e+00 5.245016905437e+00 7.247951793104e-01
2.951287593610e-01 2.736929107543e+00 5.982742149007e+00 5.234899373996e+00 1.720020377149e-01 2.736929107543e+00 5.982742149007e+00 5.234899373996e+00 2.844004075430e+02 3.000000000000e-01 8.173299183221e+03 -2.043324795805e-12 2.736929107543e+00 5.982742149007e+00 5.234899373996e+00 7.094962563442e-01
3.218318612422e-01 2.793359401653e+00 5.967793186703e+00 5.221819032034e+00 1.818096892248e-01 2.793359401653e+00 5.967793186703e+00 5.221819032034e+00 2.863619378450e+02 3.000000000000e-01 8.847071225812e+03 -2.211767806453e-12 2.793359401653e+00 5.967793186703e+00 5.221819032034e+00 6.881506008982e-01
3.485349631234e-01 2.849789695762e+00 5.954267376847e+00 5.209983948458e+00 1.909849275262e-01 2.849789695762e+00 5.954267376847e+00 5.209983948458e+00 2.881969855052e+02 3.000000000000e-01 9.462740134073e+03 -2.365685033518e-12 2.849789695762e+00 5.954267376847e+00 5.209983948458e+00 6.690124690647e-01
3.752380650046e-01 2.906219989872e+00 5.942062205026e+00 5.199304423159e+00 1.995731367718e-01 2.906219989872e+00 5.942062205026e+00 5.199304423159e+00 2.899146273544e+02 3.000000000000e-01 1.005470743230e+04 -2.513676858075e-12 2.906219989872e+00 5.942062205026e+00 5.199304423159e+00 6.508190624931e-01
4.019411668859e-01 2.962650283981e+00 5.931068668449e+00 5.189685078695e+00 2.076177461243e-01 2.962650283981e+00 5.931068668449e+00 5.189685078695e+00 2.915235492249e+02 3.000000000000e-01 1.062155026208e+04 -2.655387565520e-12 2.962650283981e+00 5.931068668449e+00 5.189685078695e+00 6.335859731219e-01
4.397050557237e-01 3.042454771240e+00 5.917465114601e+00 5.177781969127e+00 2.181345773202e-01 3.042454771240e+00 5.917465114601e+00 5.177781969127e+00 2.936269154640e+02 3.000000000000e-01 1.141956607806e+04 -2.854891519515e-12 3.042454771240e+00 5.917465114601e+00 5.177781969127e+00 6.094511864668e-01
4.774689445616e-01 3.122259258498e+00 5.905735483707e+00 5.167518542137e+00 2.277572252273e-01 3.122259258498e+00 5.905735483707e+00 5.167518542137e+00 2.955514450455e+02 3.000000000000e-01 1.212545597493e+04 -3.031363993733e-12 3.122259258498e+00 5.905735483707e+00 5.167518542137e+00 5.885545923136e-01
5.152328333994e-01 3.202063745757e+00 5.895724462413e+00 5.158758898542e+00 2.365670034158e-01 3.202063745757e+00 5.895724462413e+00 5.158758898542e+00 2.973134006832e+02 3.000000000000e-01 1.279606236869e+04 -3.199015592173e-12 3.202063745757e+00 5.895724462413e+00 5.158758898542e+00 5.688879639509e-01
5.529967222373e-01 3.281868233015e+00 5.887241650278e+00 5.151336437956e+00 2.446429006647e-01 3.281868233015e+00 5.887241650278e+00 5.151336437956e+00 2.989285801329e+02 3.000000000000e-01 1.342768644104e+04 -3.356921610259e-12 3.281868233015e+00 5.887241650278e+00 5.151336437956e+00 5.505390242489e-01
6.064029259997e-01 3.394728821234e+00 5.877690473912e+00 5.142979158672e+00 2.549251867451e-01 3.394728821234e+00 5.877690473912e+00 5.142979158672e+00 3.009850373490e+02 3.000000000000e-01 1.431703452006e+04 -3.579258630015e-12 3.394728821234e+00 5.877690473912e+00 5.142979158672e+00 5.246909803787e-01
6.598091297621e-01 3.507589409453e+00 5.870339720884e+00 5.136547249801e+00 2.640665111881e-01 3.507589409453e+00 5.870339720884e+00 5.136547249801e+00 3.028133022376e+02 3.000000000000e-01 1.507434169973e+04 -3.768585424932e-12 3.507589409453e+00 5.870339720884e+00 5.136547249801e+00 5.031638272124e-01
7.132153335246e-01 3.620449997672e+00 5.864994318857e+00 5.131870023049e+00 2.721946357241e-01 3.620449997672e+00 5.864994318857e+00 5.131870023049e+00 3.044389271448e+02 3.000000000000e-01 1.578383603836e+04 -3.945959009590e-12 3.620449997672e+00 5.864994318857e+00 5.131870023049e+00 4.830992390316e-01
7.666215372870e-01 3.733310585891e+00 5.861355283332e+00 5.128685866982e+00 2.794364602990e-01 3.733310585891e+00 5.861355283332e+00 5.128685866982e+00 3.058872920598e+02 3.000000000000e-01 1.643895695914e+04 -4.109739239785e-12 3.733310585891e+00 5.861355283332e+00 5.128685866982e+00 4.646884724024e-01
8.421493149627e-01 3.892919560408e+00 5.858980589923e+00 5.126608010262e+00 2.883280523790e-01 3.892919560408e+00 5.858980589923e+00 5.126608010262e+00 3.076656104758e+02 3.000000000000e-01 1.736203881728e+04 -4.340509704320e-12 3.892919560408e+00 5.858980589923e+00 5.126608010262e+00 4.385267121282e-01
9.176770926384e-01 4.052528534924e+00 5.858854548552e+00 5.126497724069e+00 2.959296327650e-01 4.052528534924e+00 5.858854548552e+00 5.126497724069e+00 3.091859265530e+02 3.000000000000e-01 1.811460784529e+04 -4.528651961322e-12 4.052528534924e+00 5.858854548552e+00 5.126497724069e+00 4.176217697608e-01
9.932048703141e-01 4.212137509441e+00 5.860789832842e+00 5.128191097822e+00 3.024135295020e-01 4.212137509441e+00 5.860789832842e+00 5.128191097822e+00 3.104827059004e+02 3.000000000000e-01 1.880883353773e+04 -4.702208384433e-12 4.212137509441e+00 5.860789832842e+00 5.128191097822e+00 3.983097508944e-01
1.068732647990e+00 4.226497308104e+00 5.860155047148e+00 5.127635660341e+00 3.030663871420e-01 4.226497308104e+00 5.860155047148e+00 5.127635660341e+00 3.106132774284e+02 3.000000000000e-01 3.009037434557e+04 -7.522593586392e-12 4.226497308104e+00 5.860155047148e+00 5.127635660341e+00 3.783822564105e-02
1.168732647990e+00 4.226497308104e+00 5.860223896475e+00 5.127695903508e+00 3.030573078965e-01 4.226497308104e+00 5.860223896475e+00 5.127695903508e+00 3.106114615793e+02 3.000000000000e-01 3.127037575283e+04 -7.817593938207e-12 4.226497308104e+00 5.860223896475e+00 5.127695903508e+00 9.592621982170e-05
1.268732647990e+00 4.226497308104e+00 5.859637732286e+00 5.127183009844e+00 3.031346163224e-01 4.226497308104e+00 5.859637732286e+00 5.127183009844e+00 3.106269232645e+02 3.000000000000e-01 3.124968721564e+04 -7.812421803910e-12 4.226497308104e+00 5.859637732286e+00 5.127183009844e+00 8.168166969049e-04
1.368732647990e+00 4.226497308104e+00 5.859051568098e+00 5.126670116181e+00 3.032119479523e-01 4.226497308104e+00 5.859051568098e+00 5.126670116181e+00 3.106423895905e+02 3.000000000000e-01 3.125153910317e+04 -7.812884775791e-12 4.226497308104e+00 5.859051568098e+00 5.126670116181e+00 8.168984187748e-04
1.468732647990e+00 4.226497308104e+00 5.859029899787e+00 5.126651156415e+00 3.032148070595e-01 4.226497308104e+00 5.859029899787e+00 5.126651156415e+00 3.106429614119e+02 3.000000000000e-01 3.127621100798e+04 -7.819052751994e-12 4.226497308104e+00 5.859029899787e+00 5.126651156415e+00 3.019634892395e-05
1.568732647990e+00 4.226497308104e+00 5.859008231477e+00 5.126632196650e+00 3.032176661984e-01 4.226497308104e+00 5.859008231477e+00 5.126632196650e+00 3.106435332397e+02 3.000000000000e-01 3.127627962175e+04 -7.819069905437e-12 4.226497308104e+00 5.859008231477e+00 5.126632196650e+00 3.019646059958e-05
1.668732647990e+00 4.226497308104e+00 5.858072919254e+00 5.125813798454e+00 3.033411111410e-01 4.226497308104e+00 5.858072919254e+00 5.125813798454e+00 3.106682222282e+02 3.000000000000e-01 3.123940952565e+04 -7.809852381413e-12 4.226497308104e+00 5.858072919254e+00 5.125813798454e+00 1.303739657514e-03
1.768732647990e+00 4.226497308104e+00 5.857911200721e+00 5.125672294743e+00 3.033624611696e-01 4.226497308104e+00 5.857911200721e+00 5.125672294743e+00 3.106724922339e+02 3.000000000000e-01 3.127364872282e+04 -7.818412180706e-12 4.226497308104e+00 5.857911200721e+00 5.125672294743e+00 2.254121730554e-04
1.868732647990e+00 4.226497308104e+00 5.857749482189e+00 5.125530791032e+00 3.033838129665e-01 4.226497308104e+00 5.857749482189e+00 5.125530791032e+00 3.106767625933e+02 3.000000000000e-01 3.127416085580e+04 -7.818540213950e-12 4.226497308104e+00 5.857749482189e+00 5.125530791032e+00 2.254183962364e-04
1.968732647990e+00 4.226497308104e+00 5.856966938964e+00 5.124846065710e+00 3.034871576166e-01 4.226497308104e+00 5.856966938964e+00 5.124846065710e+00 3.106974315233e+02 3.000000000000e-01 3.124956514425e+04 -7.812391286062e-12 4.226497308104e+00 5.856966938964e+00 5.124846065710e+00 1.090985418048e-03
2.000000000000e+00 4.226497308104e+00 5.856722258419e+00 5.124631970233e+00 3.035194792537e-01 4.226497308104e+00 5.856722258419e+00 5.124631970233e+00 3.107038958507e+02 3.000000000000e-01 3.125034015850e+04 -7.812585039625e-12 4.226497308104e+00 5.856722258419e+00 5.124631970233e+00 1.090980904928e-03