2
0
Fork 0

ajout d un test de sortie des differentes deformations en 1D: Almansi, logarithmique, Green-Lagrange

This commit is contained in:
rio 2017-10-22 10:06:28 +00:00
parent 92ad8c8ba8
commit e0229c254e
5 changed files with 520 additions and 0 deletions

View file

@ -0,0 +1,51 @@
------------------------------------------------------
Auteur
------------------------------------------------------
Gérard Rio (gerard.rio@univ-ubs.fr)
------------------------------------------------------
Mots-cles
------------------------------------------------------
Green-Lagrange11
Almansi11
logarithmique11
------------------------------------------------------
But du test
------------------------------------------------------
Test de sortie de résultats concernant le calcul des différentes déformations suivantes, accessibles en post-traitement:
- déformation d'Almansi
- déformation de Green-Lagrange
- déformation logarithmique
------------------------------------------------------
Description du calcul
------------------------------------------------------
Étude 1D de la traction d'une barre maillée en 2 éléments.
-------------------------------------------------------------
Grandeurs de comparaison
-------------------------------------------------------------
\latex_brut
Les 3 déformations que l'on peut calculer. La loi de comportement est Hooke 1D de module 10000 MPa
Longueur: L = 400mm, force: $F_centre = 1000N$, $F_bout = 2000N$, $section: 3mm^2$
Dans ce calcul, on impose une section constante d'où une contrainte qui est exactement:
$\sigma = F / section$
La déformation utilisée pour le calcul est celle d'Almansi. D'où sur la première barre:
$\sigma = 3000/3=1000MPa$ D'où $\epsilon_A = 0.1= 10$\%
Sur la seconde barre:
$\sigma = 2000/3 MPa$ D'où $\epsilon_A = 0.2/3 = 0.0666 = 6.66 $\%
À partir de ces grandeurs, sachant que l'élongation vaut :
$\lambda = 1./\sqrt{(1.-2.\epsilon_A)} $, on en déduit :
- la déformation logarithmique : $\epsilon_L = log(\lambda) $
- la déformation de Green-Lagrange: $\epsilon_G = 0.5 (\lambda \lambda - 1.)$
D'où pour le premier tronçon:
$\epsilon_L = 0.11157$ et $\epsilon_G = 0.125$
et pour le deuxième tronçon:
$\epsilon_L = 0.071550$ et $\epsilon_G = 0.076923$
On retrouve ces différentes valeurs dans le fichier .maple
\fin_latex_brut

View file

@ -0,0 +1,138 @@
###############################################################################################
# Fichier de commande pour la visualisation elements finis #
# Herezh++ V6.811 #
# Copyright (c) 1997-2017, 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 2 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 EPS11 Green-Lagrange11 Almansi11 logarithmique11 Def_principaleI fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
# informations particuliere dans le cas ou il y a une animation
# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements)
# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent
# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme.
debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier
# un parametre de pilotage du style de sortie
parametre_style_de_sortie 1
# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc
# (tangent pour les coques, suivant la fibre moyenne pour les element 1D )
tenseur_en_absolue_ 0
fin_grandeurs_maple # fin des grandeurs a visualiser au format maple
# ----------------------------- definition des parametres d'animation: ----------------
debut_animation # un mot cle de debut de liste (debut_animation)
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
# des parametres avec des valeurs: (sur une meme ligne)
cycleInterval 8 # cycleInterval <un reel> (indique le temps en seconde du cycle de l'animation)
fin_animation # un mot cle de fin
fin_visualisation_maple
# =================================================================================
# || fin de la visualisation maple ||
# =================================================================================
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
# l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus
# sans le mot cle de fin de fichier, le fichier n'est pas valide
###############################################################################################

View file

@ -0,0 +1,104 @@
#############################################################################
# #
# msh2her.pl #
# transformation d'un maillage genere par le programme gmsh #
# en un maillage utilisable par Herezh++ #
#############################################################################
# licence GPL, gerard.rio@univ-ubs.fr #
# http://www-lg2m.univ-ubs.fr/~rio #
# version 2.2 #
#############################################################################
nom_maillage barre1d
noeuds ------------
3 NOEUDS
#---------------------------------------------------------------
#|NO DU| X | Y | Z |
#|NOEUD| | | |
#---------------------------------------------------------------
1 0 0 0
2 200 0 0
3 400 0 0
#---------------------------------------------------------------
# references de noeuds |
#---------------------------------------------------------------
N_1_2 2
N_encastrement 1
N_LIGNE_poutre 1 2 3
N_1_3 3
N_milieu 2
N_extremite 3
N_1_1 1
#---------------------------------------------------------------
# les elements |
#---------------------------------------------------------------
elements ------------
2 ELEMENTS
#----------------------------------------------------------------------
#| NO | | |
#|ELTS | type element | Noeuds |
#----------------------------------------------------------------------
1 POUT BIE1 1 2
2 POUT BIE1 2 3
#---------------------------------------------------------------
# references elements |
#---------------------------------------------------------------
E_2_1 1
E_2_2 2
E_poutre 1 2
#---------------------------------------------------------------
# references aretes elements |
#---------------------------------------------------------------
A_poutre 1 1 2 1
#---------------------------------------------------------------
# references de faces elements |
#---------------------------------------------------------------

View file

@ -0,0 +1,119 @@
#######################################################################
# #
# | | ==== === ==== ==== | | | | #
# | | | | | | / | | | | #
# |====| |=== === |=== / |====| ------- ------- #
# | | | | \ | / | | | | #
# | | ==== | \ ==== ==== | | | | #
# #
#######################################################################
# Copyright (c) 1997-2017 Gerard Rio, gerard.rio@univ-ubs.fr #
# All rights reserved. http://http://kleger.univ-ubs.fr/Herezh/ #
# Certification IDDN.FR.010.0106078.000.R.P.2006.035.20600 #
# #
#######################################################################
# (version la plus rapide )
version 6.811
#######################################################################
# fichier de commande #
# --- traction d'une barre soumise a une force en bout et une force au milieu
# objectifs: sortie en post-traitement des differentes deformations
#----------------------------------------------------
# definition de la dimension de l'espace de travail |
#----------------------------------------------------
dimension 1
#---------------------------------------------------------------
# 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
< barre1d.her
#======================================
#| definition des lois de comportement|
#--------------------------------------
choix_materiaux ----------
#----------------------------------------
# Elements | Nom Materiau |
#----------------------------------------
E_tout acier
materiaux #----------
#------------------------------------------
# Nom Materiau | Type loi |
#---------------------------
acier ISOELAS1D
# ....... loi de comportement isoelastique 2D, contraintes planes ........
# module d'young : coefficient de poisson
10000 0.3
#--------------------------------------- fin def des lois de comportement ------------------------
# --- divers stockages (1) -------
sections #--------#
E_tout 3
variation_section
E_tout 0
masse_volumique #--------#
E_tout 1
charges #------------#
N_milieu PONCTUELLE 1e+03
N_extremite PONCTUELLE 2e+03
blocages #------------#
#-----------------------------------------------------------
# nom du maillage | Ref noeud | Blocages
#-----------------------------------------------------------
N_encastrement UX
controle ------------
#---------------------------
# PARAMETRE | VALEUR |
#---------------------------
SAUVEGARDE 0
DELTAtMAXI 0.1
TEMPSFIN 1.
DELTAt 0.1
MULTIPLICATEUR 1.
ITERATIONS 15
PRECISION 1e-3
para_affichage #------------
#---------------------------
# PARAMETRE | VALEUR |
#---------------------------
FREQUENCE_SORTIE_FIL_DU_CALCUL 1
# ------------------------------------------------------------------------------------
resultats pas_de_sortie_finale_
COPIE 0
#
_fin_point_info_

View file

@ -0,0 +1,108 @@
#fichier au format maple6
###############################################################################################
# Visualisation elements finis : Herezh++ V6.811 #
# Copyright (c) 1997-2017, 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 1 (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 ------------------------------
#2 5 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] EPS11 [4] Green-Lagrange11 [5] Almansi11 [6] logarithmique11 [7] Def_principaleI
# element_2 pt_integ_1: [8]X [9] EPS11 [10] Green-Lagrange11 [11] Almansi11 [12] logarithmique11 [13] Def_principaleI
#
#
#====================================================================
#|| fin du recapitulatif des differentes grandeurs ||
#====================================================================
# ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant
# a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun
# correspondant a un maillage. On trouve pour chaque enregistrement successivement :
# s'il y a des grandeurs globales: le temps puis les grandeurs globales,
# puis s'il y a des torseurs de reaction :
# de nouveau le temps, les composantes de la resultante puis les composantes du moments
# donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels
# puis s'il y a des grandeurs aux noeuds: de nouveau le temps
# les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud
# puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds
# puis s'il y a des grandeur aux elements:
# le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques)
# suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration
# correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes
# puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs
# etc. pour tous les points d'integration - element
1.000000000000e-01 1.010152510000e+02 9.999966479154e-03 1.020404672965e-02 9.999966479154e-03 1.010131945382e-02 9.999966479154e-03 3.027039096049e+02 6.666660029993e-03 6.756749939502e-03 6.666660029993e-03 6.711503439712e-03 6.666660029993e-03
2.000000000000e-01 1.020620579690e+02 1.999986222984e-02 2.083318384319e-02 1.999986222984e-02 2.041085374957e-02 1.999986222984e-02 3.054847198832e+02 1.333330630490e-02 1.369860160725e-02 1.333330630490e-02 1.351430842502e-02 1.333330630490e-02
3.000000000000e-01 1.031421085400e+02 2.999985339907e-02 3.191472770384e-02 2.999985339907e-02 3.093754590063e-02 2.999985339907e-02 3.083462867023e+02 1.999997184190e-02 2.083330277985e-02 1.999997184190e-02 2.041096792878e-02 1.999997184190e-02
4.000000000000e-01 1.042571893575e+02 3.999984406529e-02 4.347807663675e-02 3.999984406529e-02 4.169063497531e-02 3.999984406529e-02 3.112926820050e+02 2.666663733784e-02 2.816898135794e-02 2.666663733784e-02 2.740408726634e-02 2.666663733784e-02
5.000000000000e-01 1.054092358756e+02 4.999983381926e-02 5.555535039422e-02 4.999983381926e-02 5.268007318368e-02 4.999983381926e-02 3.143283022619e+02 3.333330275930e-02 3.571425061654e-02 3.333330275930e-02 3.449640298558e-02 3.333330275930e-02
6.000000000000e-01 1.066003366948e+02 5.999982265495e-02 6.818158917229e-02 5.999982265495e-02 6.391648422652e-02 5.999982265495e-02 3.174578768042e+02 3.999996810923e-02 4.347822319144e-02 3.999996810923e-02 4.169076980565e-02 3.999996810923e-02
7.000000000000e-01 1.078327494386e+02 6.999981046772e-02 8.139509257410e-02 6.999981046772e-02 7.541122448097e-02 6.999981046772e-02 3.206865013239e+02 4.666663338246e-02 5.147054774575e-02 4.666663338246e-02 4.899016746958e-02 4.666663338246e-02
8.000000000000e-01 1.091089187679e+02 7.999979713756e-02 9.523780773477e-02 7.999979713756e-02 8.717645206954e-02 7.999979713756e-02 3.240196757977e+02 5.333329857331e-02 5.970144898081e-02 5.333329857331e-02 5.639770816220e-02 5.333329857331e-02
9.000000000000e-01 1.104314967872e+02 8.999978252725e-02 1.097557741335e-01 8.999978252725e-02 9.922520415131e-02 8.999978252725e-02 3.274633473521e+02 5.999996367559e-02 6.818177127530e-02 5.999996367559e-02 6.391664447721e-02 5.999996367559e-02
1.000000000000e+00 1.118033662395e+02 9.999976647933e-02 1.249996351242e-01 9.999976647933e-02 1.115714837563e-01 9.999976647933e-02 3.310239588770e+02 6.666662868257e-02 7.692302635253e-02 6.666662868257e-02 7.155037799253e-02 6.666662868257e-02