diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/README b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/README new file mode 100755 index 0000000..2a9da2f --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/README @@ -0,0 +1,39 @@ +------------------------------------------------------------ +Auteur +------------------------------------------------------------ +Gérard Rio (gerard.rio@univ-ubs.fr) +Julien Troufflard (troufflard@univ-ubs.fr) + +------------------------------------------------------------ +Mots-cles +------------------------------------------------------------ +axisymetrique +contact rigide deformable +ISOELAS + +------------------------------------------------------------ +But du test +------------------------------------------------------------ + +L'objectif est de vérifier le bon fonctionnement d'un contact basique, en axisymétrique. +- on vérifie le déplacement +- on vérifie la force de réaction +- on vérifie que l'ensemble fonctionne et qu'il n'y a pas d'erreur dans le processus de calcul. + +------------------------------------------------------------ +Description du calcul +------------------------------------------------------------ +Le maillage est constitué de deux éléments rectangulaires axisymétriques décrivant ainsi un barreau cylindrique avec un tube creux. +L'élément externe (le plus grand) est bloqué et l'élément interne est comprimé via un déplacement nul à la base suivant y et un déplacement imposé sur la partie haute. Cela conduit à un contact d'où des réactions et une pénétration qui est gérée par la méthode de pénalisation utilisée ici pour le contact. + +\figures: maillage.png[0.2] [1] force_reaction.png[0.2] \legende: À gauche le maillage constitué de 2 éléments et à droite les forces de réaction \fin_legende + +On peut observer par exemple les forces de réaction via les grandeurs évoluées aux noeuds: VECT_REAC. + + +------------------------------------------------------------ +Grandeurs de comparaison +------------------------------------------------------------ + +- déplacements et réaction pour les noeuds en contact : noeuds 3 et 4 du maillage 1 + et pour les noeuds du solide rigide : noeuds 1 et 4 du maillage 2 diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/force_reaction.png b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/force_reaction.png new file mode 100755 index 0000000..020299c Binary files /dev/null and b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/force_reaction.png differ diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/maillage.png b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/maillage.png new file mode 100755 index 0000000..46204bb Binary files /dev/null and b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/maillage.png differ diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plateau_axi.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plateau_axi.her new file mode 100755 index 0000000..9cc444c --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plateau_axi.her @@ -0,0 +1,35 @@ + nom_maillage plateau_axi + + noeuds + 4 NOEUDS + + 1 9 -10 0 + 2 29 -10 0 + 3 29 30 0 + 4 9 30 0 + + N_O 1 4 + + N_1_4 4 + + N_1_3 3 + + N_1_2 2 + + N_SURFACE_0 1 2 3 4 + + N_LIGNE_0 1 2 3 4 + + N_0 1 2 3 4 + + + elements + 1 ELEMENTS + + 1 QUAD_AXI LINEAIRE 1 2 3 4 + + E_2_6 1 + + E_0 1 + + F_0 1 1 diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plot_axi.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plot_axi.her new file mode 100755 index 0000000..3dc39c3 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/plot_axi.her @@ -0,0 +1,76 @@ + +############################################################################# +# # +# STAMM # +# (maillage automatique de pieces simples) # +############################################################################# +# Copyright © 1997-2013 Gerard Rio, gerard.rio@univ-ubs.fr # +# All rights reserved. http://www-lg2m.univ-ubs.fr/~rio # +# # +############################################################################# +# version 03.29 , compatible avec la version Herezh++6.576 + + + # ---------------------------------------------------------------- + # | construction automatique d'un maillage 2D | + # | pour une utilisation par herezh++ | + # ---------------------------------------------------------------- + # rectangle de dimension : 9 x 30 ( 1 x 1 elements ) + # geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire. + + nom_maillage plot_axi + + noeuds ------------ + 4 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 25.26 0 + 3 8.835 0 + 4 8.835 25.26 0 + + # references des noeuds + +N_O 1 2 +N_S 1 3 +N_E 3 4 +N_N 2 4 +N_SO 1 +N_SE 3 +N_NE 4 +N_NO 2 + N_to 1 2 3 4 + + + # les elements + + elements ---------- + 1 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 QUAD_AXI LINEAIRE 1 3 4 2 + + + # references des elements + +A_O 1 4 +A_S 1 1 +A_E 1 2 +A_N 1 3 +F_to 1 1 +E_O 1 +E_S 1 +E_E 1 +E_N 1 +E_SO 1 +E_SE 1 +E_NE 1 +E_NO 1 + E_to 1 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.CVisu1 b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.CVisu1 new file mode 100755 index 0000000..4b4c799 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.CVisu1 @@ -0,0 +1,254 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.736 # +# 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) + tous_les_increments fin_list_increment + + # ----------------------------- choix des maillages a visualiser: ---------------- + # la liste est facultative, par defaut la visualisation concerne le premier maillage + debut_choix_maillage # un mot cle de debut, + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entiers , puis , sur une meme ligne + 1 2 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 N_E fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud X1 X2 X3 R_X1 R_X2 R_X3 fin_list_ddl_noeud + type_sortie_ddl_retenue= 1 + # 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 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 + 2 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + 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 N_O fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud X1 X2 X3 R_X1 R_X2 R_X3 fin_list_ddl_noeud + type_sortie_ddl_retenue= 1 + # 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 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 || + # ================================================================================= + + + # ================================================================================= + # || ***** demande d'une visualisation Gmsh: ***** || + # ================================================================================= + # un mot cle de debut (debut_visualisation_Gmsh) + # un mot cle de fin ( fin_visualisation_Gmsh) apres tous les ordres particuliers + # la seule presence du premier mots cle suffit a activer la visualisation Gmsh + # la presence du second permet une meilleur lisibilite du fichier, mais n'est pas indispensable + debut_visualisation_Gmsh + + # ----------------------------- definition des parametres du maillage initial: ---------------- + debut_maillage_initial # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + pseudo-homothetie_sur_les_maillages_ 0 # 0 = aucune homothetie, 1 = il y en a + # --- def eventuelle de la pseudo-homothetie: une par maillage, + # pseudo-homothetie pour le maillage : 1 +# maillage_ 1 + # mot cle: maillage_ puis le numero du maillage, +# pseudo-homothetie_ 0 # 0 = non active, 1 = active + # ensuite si c'est active, on trouve : + # mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie + # puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes. +# centre_homothetie_ +# fact_mult_ + # pseudo-homothetie pour le maillage : 2 +# maillage_ 2 + # mot cle: maillage_ puis le numero du maillage, +# pseudo-homothetie_ 0 # 0 = non active, 1 = active + # ensuite si c'est active, on trouve : + # mot cle: centre_homothetie_ puis les coordonnees du centre d'homothetie + # puis mot cle: fact_mult_ puis les coordonnees donnant les coefs multiplicatifs selon les axes. +# centre_homothetie_ +# fact_mult_ + 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 X1 X2 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 1 1 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 VECT_REAC fin_list_GrandEvoluee_noeud + # tableau de ddl aux elements a visualiser, un par maillage + debut_tableau_ddl_aux_elements 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 fin_list_GrandParticuliere_element + 2 # le numero de maillage + # tableau des ddl aux noeuds a visualiser, un par maillage + debut_tableau_ddl_aux_noeuds X1 X2 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 1 1 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 VECT_REAC fin_list_GrandEvoluee_noeud + # tableau de ddl aux elements a visualiser, un par maillage + debut_tableau_ddl_aux_elements 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 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 + 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) + dernier_increment fin_list_increment + + # ----------------------------- choix des maillages a visualiser: ---------------- + # la liste est facultative, par defaut la visualisation concerne le premier maillage + debut_choix_maillage # un mot cle de debut, + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entiers , puis , sur une meme ligne + 1 2 fin_choix_maillage + + fin_visualisation_Gmsh + # ================================================================================= + # || fin de la visualisation Gmsh || + # ================================================================================= + + + + + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.info b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.info new file mode 100755 index 0000000..f214586 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI.info @@ -0,0 +1,81 @@ +dimension 3 #---# + +niveau_commentaire 3 #---# + +TYPE_DE_CALCUL #---# +non_dynamique # avec plus visualisation + +< plot_axi.her + +< plateau_axi.her + +domaine_esclave #---# +1 #nb premiers maillages definis comme esclave + +choix_materiaux #---# +nom_mail= plot_axi E_tout MAT_PLOT +nom_mail= plateau_axi E_tout MAT_OUTIL + +materiaux #---# +MAT_OUTIL ISOELAS +210000. 0.3 + +MAT_PLOT ISOELAS +210. 0.3 + +masse_volumique #---# +nom_mail= plot_axi E_tout 1. +nom_mail= plateau_axi E_tout 1. + +charges #---# + +blocages #---# +nom_mail= plot_axi N_O UX +nom_mail= plot_axi N_S UY +nom_mail= plot_axi N_tout UZ +nom_mail= plot_axi N_N 'UY=-5.635' + +nom_mail= plateau_axi N_tout UX UY UZ + +controle #---# + +DELTAtMINI 1.e-6 +DELTAtMAXI 0.5 +TEMPSFIN 1. +DELTAt 0.01 +ITERATIONS 15 +PRECISION 1.e-3 +SAUVEGARDE 1 +MAXINCRE 9999999 +NORME Residu/Reaction +#RESTART 15 + +para_contact #---# + +#CONTACT_TYPE : +# > 0 : pas de contact +# > 1 : contact avec une methode sans multiplicateur de Lagrange ni penalisation +# > 2 : contact avec une methode de penalisation +CONTACT_TYPE 2 + +PENALISATION_PENETRATION 10. + +#TYPE_PENALISATION_PENETRATION (indique la methode pour calculer le facteur beta) +# > 1 : dans ce cas beta=1 et alpha = facteur de penalisation +# > 2 : beta est calcule +# > 3 : meme fonctionnement que 2 mais avec un algorithme supplementaire +# > 5 : meme fonctionnement que 2 mais avec un algorithme supplementaire +TYPE_PENALISATION_PENETRATION 4 +NB_DECOLLEMENT_MAXI 4 +DISTANCE_MAXI_AU_PT_PROJETE 2. +TYPE_DE_DECOLLEMENT 1 +PENETRATION_CONTACT_MAXI 0.01 +PENETRATION_BORNE_REGULARISATION 0.1 + +para_affichage #---# +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + + resultats pas_de_sortie_finale_ #---# + COPIE 0 + + _fin_point_info_ diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI_princ.maple.ref1 new file mode 100755 index 0000000..910281e --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_basique_axi_deformable_solide/test_AXI_princ.maple.ref1 @@ -0,0 +1,115 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.736 # +# Copyright (c) 1997-2015, 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 + # 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 -------------------------------- +#2 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#2 6 (nombre de noeuds, nombre total de grandeurs associees) +# maillage_1 : +# reference N_E : noeud_3 [2]X [3]Y [4]Z [5]delta_0_t_X1 [6]delta_0_t_X2 [7]delta_0_t_X3 [8]delta_0_t_R_X1 [9]delta_0_t_R_X2 [10]delta_0_t_R_X3 ; +# noeud_4 [11]X [12]Y [13]Z [14]delta_0_t_X1 [15]delta_0_t_X2 [16]delta_0_t_X3 [17]delta_0_t_R_X1 [18]delta_0_t_R_X2 [19]delta_0_t_R_X3 ; +# +#2 6 (nombre de noeuds, nombre total de grandeurs associees) +# maillage_2 : +# reference N_O : noeud_1 [21]X [22]Y [23]Z [24]delta_0_t_X1 [25]delta_0_t_X2 [26]delta_0_t_X3 [27]delta_0_t_R_X1 [28]delta_0_t_R_X2 [29]delta_0_t_R_X3 ; +# noeud_4 [30]X [31]Y [32]Z [33]delta_0_t_X1 [34]delta_0_t_X2 [35]delta_0_t_X3 [36]delta_0_t_R_X1 [37]delta_0_t_R_X2 [38]delta_0_t_R_X3 ; +# +#---------------------------------- grandeurs aux elements ------------------------------ +#0 0 (nombre total d'elements, nombre totale de grandeurs associees) +#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 + + 1.000000000000e-02 8.840936565084e+00 0.000000000000e+00 0.000000000000e+00 5.936565084060e-03 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -7.695567144281e+01 0.000000000000e+00 8.840936565084e+00 2.520365000000e+01 0.000000000000e+00 5.936565084060e-03 -5.635000000000e-02 0.000000000000e+00 0.000000000000e+00 7.695567144281e+01 0.000000000000e+00 1.000000000000e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +2.000000000000e-02 8.846929200292e+00 0.000000000000e+00 0.000000000000e+00 1.192920029242e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -1.546207033397e+02 0.000000000000e+00 8.846929200292e+00 2.514730000000e+01 0.000000000000e+00 1.192920029242e-02 -1.127000000000e-01 0.000000000000e+00 0.000000000000e+00 1.546207033397e+02 0.000000000000e+00 2.000000000000e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +3.000000000000e-02 8.852972565132e+00 0.000000000000e+00 0.000000000000e+00 1.797256513183e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -2.330315826996e+02 0.000000000000e+00 8.852972565132e+00 2.509095000000e+01 0.000000000000e+00 1.797256513183e-02 -1.690500000000e-01 0.000000000000e+00 0.000000000000e+00 2.330315826996e+02 0.000000000000e+00 3.000000000000e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +4.000000000000e-02 8.859069312393e+00 0.000000000000e+00 0.000000000000e+00 2.406931239320e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -3.121893400754e+02 0.000000000000e+00 8.859069312393e+00 2.503460000000e+01 0.000000000000e+00 2.406931239320e-02 -2.254000000000e-01 0.000000000000e+00 0.000000000000e+00 3.121893400754e+02 0.000000000000e+00 4.000000000000e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +5.414213562373e-02 8.867783844327e+00 0.000000000000e+00 0.000000000000e+00 3.278384432655e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -4.254300001496e+02 0.000000000000e+00 8.867783844327e+00 2.495490906576e+01 0.000000000000e+00 3.278384432655e-02 -3.050909342397e-01 0.000000000000e+00 0.000000000000e+00 4.254300001496e+02 0.000000000000e+00 5.414213562373e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +6.828427124746e-02 8.876608406252e+00 0.000000000000e+00 0.000000000000e+00 4.160840625242e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -5.402138945501e+02 0.000000000000e+00 8.876608406252e+00 2.487521813152e+01 0.000000000000e+00 4.160840625242e-02 -3.847818684794e-01 0.000000000000e+00 0.000000000000e+00 5.402138945501e+02 0.000000000000e+00 6.828427124746e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +8.242640687119e-02 8.885544970661e+00 0.000000000000e+00 0.000000000000e+00 5.054497066124e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -6.565709882637e+02 0.000000000000e+00 8.885544970661e+00 2.479552719728e+01 0.000000000000e+00 5.054497066124e-02 -4.644728027192e-01 0.000000000000e+00 0.000000000000e+00 6.565709882637e+02 0.000000000000e+00 8.242640687119e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +9.656854249492e-02 8.894595556473e+00 0.000000000000e+00 0.000000000000e+00 5.959555647312e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -7.745320136684e+02 0.000000000000e+00 8.894595556473e+00 2.471583626304e+01 0.000000000000e+00 5.959555647312e-02 -5.441637369589e-01 0.000000000000e+00 0.000000000000e+00 7.745320136684e+02 0.000000000000e+00 9.656854249492e-02 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +1.165685424949e-01 8.907593688930e+00 0.000000000000e+00 0.000000000000e+00 7.259368893045e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -9.441535902343e+02 0.000000000000e+00 8.907593688930e+00 2.460313626304e+01 0.000000000000e+00 7.259368893045e-02 -6.568637369589e-01 0.000000000000e+00 0.000000000000e+00 9.441535902343e+02 0.000000000000e+00 1.165685424949e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +1.365685424949e-01 8.920830024037e+00 0.000000000000e+00 0.000000000000e+00 8.583002403732e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -1.117138173821e+03 0.000000000000e+00 8.920830024037e+00 2.449043626304e+01 0.000000000000e+00 8.583002403732e-02 -7.695637369589e-01 0.000000000000e+00 0.000000000000e+00 1.117138173821e+03 0.000000000000e+00 1.365685424949e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +1.565685424949e-01 8.934310785322e+00 0.000000000000e+00 0.000000000000e+00 9.931078532201e-02 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -1.293581145060e+03 0.000000000000e+00 8.934310785322e+00 2.437773626304e+01 0.000000000000e+00 9.931078532201e-02 -8.822637369589e-01 0.000000000000e+00 0.000000000000e+00 1.293581145060e+03 0.000000000000e+00 1.565685424949e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +1.765685424949e-01 8.948042410615e+00 0.000000000000e+00 0.000000000000e+00 1.130424106154e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -1.473581465122e+03 0.000000000000e+00 8.948042410615e+00 2.426503626304e+01 0.000000000000e+00 1.130424106154e-01 -9.949637369589e-01 0.000000000000e+00 0.000000000000e+00 1.473581465122e+03 0.000000000000e+00 1.765685424949e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +2.048528137424e-01 8.967903098670e+00 0.000000000000e+00 0.000000000000e+00 1.329030986698e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -1.734413476677e+03 0.000000000000e+00 8.967903098670e+00 2.410565439456e+01 0.000000000000e+00 1.329030986698e-01 -1.154345605438e+00 0.000000000000e+00 0.000000000000e+00 1.734413476677e+03 0.000000000000e+00 2.048528137424e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +2.331370849898e-01 8.988298538725e+00 0.000000000000e+00 0.000000000000e+00 1.532985387249e-01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 -2.002870240885e+03 0.000000000000e+00 8.988298538725e+00 2.394627252608e+01 0.000000000000e+00 1.532985387249e-01 -1.313727473918e+00 0.000000000000e+00 0.000000000000e+00 2.002870240885e+03 0.000000000000e+00 2.331370849898e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 +2.614213562373e-01 9.000011442868e+00 0.000000000000e+00 0.000000000000e+00 1.650114428681e-01 0.000000000000e+00 0.000000000000e+00 -2.682585767084e+02 -2.315198505388e+03 0.000000000000e+00 9.000011442868e+00 2.378689065760e+01 0.000000000000e+00 1.650114428681e-01 -1.473109342397e+00 0.000000000000e+00 -2.682585767084e+02 2.315198505388e+03 0.000000000000e+00 2.614213562373e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 2.417930404297e+02 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 2.947241129872e+02 0.000000000000e+00 0.000000000000e+00 +2.897056274848e-01 9.000037670876e+00 0.000000000000e+00 0.000000000000e+00 1.650376708761e-01 0.000000000000e+00 0.000000000000e+00 -8.835928710791e+02 -2.681142759682e+03 0.000000000000e+00 9.000037670876e+00 2.362750878912e+01 0.000000000000e+00 1.650376708761e-01 -1.632491210877e+00 0.000000000000e+00 -8.835928710791e+02 2.681142759682e+03 0.000000000000e+00 2.897056274848e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 7.999411313637e+02 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.672446107944e+02 0.000000000000e+00 0.000000000000e+00 +3.297056274848e-01 9.000075010658e+00 0.000000000000e+00 0.000000000000e+00 1.650750106585e-01 0.000000000000e+00 0.000000000000e+00 -1.760733418677e+03 -3.211495375873e+03 0.000000000000e+00 9.000075010658e+00 2.340210878912e+01 0.000000000000e+00 1.650750106585e-01 -1.857891210877e+00 0.000000000000e+00 -1.760733418677e+03 3.211495375873e+03 0.000000000000e+00 3.297056274848e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.601056519888e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.920410317466e+03 0.000000000000e+00 0.000000000000e+00 +3.697056274848e-01 9.000112669301e+00 0.000000000000e+00 0.000000000000e+00 1.651126693007e-01 0.000000000000e+00 0.000000000000e+00 -2.646689657992e+03 -3.757408559212e+03 0.000000000000e+00 9.000112669301e+00 2.317670878912e+01 0.000000000000e+00 1.651126693007e-01 -2.083291210877e+00 0.000000000000e+00 -2.646689657992e+03 3.757408559212e+03 0.000000000000e+00 3.697056274848e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 2.421581504303e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 2.871797811680e+03 0.000000000000e+00 0.000000000000e+00 +4.097056274848e-01 9.000150650513e+00 0.000000000000e+00 0.000000000000e+00 1.651506505125e-01 0.000000000000e+00 0.000000000000e+00 -3.541583632722e+03 -4.319496766342e+03 0.000000000000e+00 9.000150650513e+00 2.295130878912e+01 0.000000000000e+00 1.651506505125e-01 -2.308691210877e+00 0.000000000000e+00 -3.541583632722e+03 4.319496766342e+03 0.000000000000e+00 4.097056274848e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 3.260319136385e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 3.822848129059e+03 0.000000000000e+00 0.000000000000e+00 +4.497056274848e-01 9.000188964172e+00 0.000000000000e+00 0.000000000000e+00 1.651889641715e-01 0.000000000000e+00 0.000000000000e+00 -4.445683902658e+03 -4.898405060436e+03 0.000000000000e+00 9.000188964172e+00 2.272590878912e+01 0.000000000000e+00 1.651889641715e-01 -2.534091210877e+00 0.000000000000e+00 -4.445683902658e+03 4.898405060436e+03 0.000000000000e+00 4.497056274848e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 4.117669253268e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 4.773698552047e+03 0.000000000000e+00 0.000000000000e+00 +5.062741699797e-01 9.000243736296e+00 0.000000000000e+00 0.000000000000e+00 1.652437362963e-01 0.000000000000e+00 0.000000000000e+00 -5.740557052386e+03 -5.747141922352e+03 0.000000000000e+00 9.000243736296e+00 2.240714505216e+01 0.000000000000e+00 1.652437362963e-01 -2.852854947836e+00 0.000000000000e+00 -5.740557052386e+03 5.747141922352e+03 0.000000000000e+00 5.062741699797e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 5.349351179297e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.131762925475e+03 0.000000000000e+00 0.000000000000e+00 +5.628427124746e-01 9.000299223387e+00 0.000000000000e+00 0.000000000000e+00 1.652992233871e-01 0.000000000000e+00 0.000000000000e+00 -7.055211127456e+03 -6.632915325085e+03 0.000000000000e+00 9.000299223387e+00 2.208838131521e+01 0.000000000000e+00 1.652992233871e-01 -3.171618684794e+00 0.000000000000e+00 -7.055211127456e+03 6.632915325085e+03 0.000000000000e+00 5.628427124746e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.630638213521e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 7.479784041392e+03 0.000000000000e+00 0.000000000000e+00 +6.194112549695e-01 9.000355456969e+00 0.000000000000e+00 0.000000000000e+00 1.653554569692e-01 0.000000000000e+00 0.000000000000e+00 -8.390509984291e+03 -7.557910957489e+03 0.000000000000e+00 9.000355456969e+00 2.176961757825e+01 0.000000000000e+00 1.653554569692e-01 -3.490382421753e+00 0.000000000000e+00 -8.390509984291e+03 7.557910957489e+03 0.000000000000e+00 6.194112549695e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 7.952445377385e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 8.828574591197e+03 0.000000000000e+00 0.000000000000e+00 +6.759797974645e-01 9.000412471038e+00 0.000000000000e+00 0.000000000000e+00 1.654124710379e-01 0.000000000000e+00 0.000000000000e+00 -9.747384734411e+03 -8.524478146237e+03 0.000000000000e+00 9.000412471038e+00 2.145085384129e+01 0.000000000000e+00 1.654124710379e-01 -3.809146158712e+00 0.000000000000e+00 -9.747384734411e+03 8.524478146237e+03 0.000000000000e+00 6.759797974645e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 9.316156150213e+03 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.017861331861e+04 0.000000000000e+00 0.000000000000e+00 +7.559797974645e-01 9.000494508399e+00 0.000000000000e+00 0.000000000000e+00 1.654945083987e-01 0.000000000000e+00 0.000000000000e+00 -1.170516853974e+04 -9.967327690208e+03 0.000000000000e+00 9.000494508399e+00 2.100005384129e+01 0.000000000000e+00 1.654945083987e-01 -4.259946158712e+00 0.000000000000e+00 -1.170516853974e+04 9.967327690208e+03 0.000000000000e+00 7.559797974645e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.128060632127e+04 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.212973075821e+04 0.000000000000e+00 0.000000000000e+00 +8.359797974645e-01 9.000578284021e+00 0.000000000000e+00 0.000000000000e+00 1.655782840212e-01 0.000000000000e+00 0.000000000000e+00 -1.371099508705e+04 -1.150623764702e+04 0.000000000000e+00 9.000578284021e+00 2.054925384129e+01 0.000000000000e+00 1.655782840212e-01 -4.710746158712e+00 0.000000000000e+00 -1.371099508705e+04 1.150623764702e+04 0.000000000000e+00 8.359797974645e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.336820175444e+04 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.405378841967e+04 0.000000000000e+00 0.000000000000e+00 +9.159797974645e-01 9.000663914343e+00 0.000000000000e+00 0.000000000000e+00 1.656639143429e-01 0.000000000000e+00 0.000000000000e+00 -1.576808739818e+04 -1.314992502031e+04 0.000000000000e+00 9.000663914343e+00 2.009845384129e+01 0.000000000000e+00 1.656639143429e-01 -5.161546158712e+00 0.000000000000e+00 -1.576808739818e+04 1.314992502031e+04 0.000000000000e+00 9.159797974645e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.555157033385e+04 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.598460446251e+04 0.000000000000e+00 0.000000000000e+00 +9.959797974645e-01 9.000751528731e+00 0.000000000000e+00 0.000000000000e+00 1.657515287308e-01 0.000000000000e+00 0.000000000000e+00 -1.788002890817e+04 -1.490811846214e+04 0.000000000000e+00 9.000751528731e+00 1.964765384129e+01 0.000000000000e+00 1.657515287308e-01 -5.612346158712e+00 0.000000000000e+00 -1.788002890817e+04 1.490811846214e+04 0.000000000000e+00 9.959797974645e-01 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.783601996996e+04 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.792403784638e+04 0.000000000000e+00 0.000000000000e+00 +1.000000000000e+00 9.000755931576e+00 0.000000000000e+00 0.000000000000e+00 1.657559315756e-01 0.000000000000e+00 0.000000000000e+00 -1.798635128945e+04 -1.499968481569e+04 0.000000000000e+00 9.000755931576e+00 1.962500000000e+01 0.000000000000e+00 1.657559315756e-01 -5.635000000000e+00 0.000000000000e+00 -1.798635128945e+04 1.499968481569e+04 0.000000000000e+00 1.000000000000e+00 9.000000000000e+00 -1.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.794208065507e+04 0.000000000000e+00 0.000000000000e+00 9.000000000000e+00 3.000000000000e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.803062192383e+04 0.000000000000e+00 0.000000000000e+00