From 24a005eeb781d23190e9563ab433ad6ed4982c3f Mon Sep 17 00:00:00 2001 From: rio <> Date: Mon, 22 Jun 2015 13:40:00 +0000 Subject: [PATCH] introduction des cas test de Bertille. --- .../Test_L_tunnel_ISOELAS_CP/KDR.algo | 19 + .../Test_L_tunnel_ISOELAS_CP/README.txt | 44 ++ .../tunnel_intermediaire.her | 198 +++++++ .../tunnel_isoelast_CP.CVisu1 | 114 ++++ .../tunnel_isoelast_CP.info | 139 +++++ .../tunnel_isoelast_CP_princ.maple.ref1 | 87 +++ .../KDR.algo | 20 + .../README.txt | 36 ++ .../tunnel_12x1.her | 115 ++++ .../tunnel_hyperelast.CVisu1 | 113 ++++ .../tunnel_hyperelast.info | 145 +++++ .../tunnel_hyperelast_princ.maple.ref1 | 78 +++ .../README.txt | 36 ++ .../quadrangle1.her | 76 +++ .../quadrangle_traction.CVisu1 | 112 ++++ .../quadrangle_traction.info | 149 +++++ .../quadrangle_traction_princ.maple.ref1 | 113 ++++ .../Test_L_tunnel_ISOELAS2D_C/.DS_Store | Bin 0 -> 6148 bytes .../Test_L_tunnel_ISOELAS2D_C/KDR.algo | 19 + .../Test_L_tunnel_ISOELAS2D_C/README.txt | 38 ++ .../Test_L_tunnel_ISOELAS2D_C/tunnel_bis.her | 552 ++++++++++++++++++ .../tunnel_isoelast.CVisu1 | 115 ++++ .../tunnel_isoelast.info | 136 +++++ .../tunnel_isoelast_princ.maple.ref1 | 87 +++ .../.DS_Store | Bin 0 -> 6148 bytes .../README.txt | 31 + .../cube1ele.her | 133 +++++ .../cube_traction.CVisu1 | 114 ++++ .../cube_traction.info | 165 ++++++ .../cube_traction_princ.maple.ref1 | 113 ++++ 30 files changed, 3097 insertions(+) create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/KDR.algo create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/README.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_intermediaire.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP_princ.maple.ref1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/KDR.algo create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/README.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_12x1.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast_princ.maple.ref1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/README.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle1.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction_princ.maple.ref1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/.DS_Store create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/KDR.algo create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/README.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_bis.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast_princ.maple.ref1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/.DS_Store create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/README.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube1ele.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction_princ.maple.ref1 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/KDR.algo b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/KDR.algo new file mode 100644 index 0000000..6c9a4be --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/KDR.algo @@ -0,0 +1,19 @@ +# ----- KDR +typeCalRelaxation= 1 lambda= 1.6 type_calcul_mass= 2 option_recalcul_mass= 1 +parametre_calcul_de_la_masse_ casMass_relax= 3 + +# ----- Amortissement cinetique +avec_amortissement_cinetique_ + max_nb_decroit_pourRelaxDyn_ 1 + coef_arret_pourRelaxDyn_ 0. + coef_redemarrage_pourRelaxDyn_ 0. + max_deltaX_pourRelaxDyn_ 0 #0.005 +nb_max_dX_OK_pourRelaxDyn_ 10 #1 +nb_deb_testfin_pourRelaxDyn_ 250 #1 +nb_deb_test_amort_cinetique_ 250 #1 +fi_parametre_amortissement_cinetique_ + +#mode_debug_= 5 + +ARRET_A_EQUILIBRE_STATIQUE_ 2 + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/README.txt b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/README.txt new file mode 100644 index 0000000..0716396 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/README.txt @@ -0,0 +1,44 @@ +------------------------------------------------------------ +But du test +------------------------------------------------------------ +Le but du calcul est de trouver la forme d’équilibre d’un tunnel gonflé +dont le matériau suit une loi de comportement isotrope élastique avec prise +en compte de la condition de contraintes planes (ISOELAS + LOI_CONTRAINTES_PLANES). + +------------------------------------------------------------ +Description du calcul +------------------------------------------------------------ +Le calcul s’effectue avec la méthode de relaxation dynamique. +Les paramètres de contrôle de cette méthode sont définis dans KDR.algo + +La loi de comportement testée est la loi ISOELAS. +Les paramètres utilisés pour la définir sont des paramètres cohérents +avec des tests réalisés sur des membranes: +E= 125 MPa +nu = 0,4 +pression = 0,39e-3 MPa +épaisseur = 25e-3 mm +masse volumique = 0,9e-9 tonnes.mm^-3 +La condition de contraintes planes est ajoutée avec la loi LOI_CONTRAINTES_PLANES, +utilisant la méthode de NEWTON_LOCAL. + +La structure est une membrane de dimension 1000x2000mm, +discrétisée en 5x10 élément membranes (interpolation linéaire - 8noeuds). + +------------------------------------------------------------ +Mots-cles +------------------------------------------------------------ +tunnel +isoelastique +contraintes planes +relaxation dynamique +ISOELAS +LOI_CONTRAINTES_PLANES + +------------------------------------------------------------ +Grandeurs de comparaison +------------------------------------------------------------ +Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII +(point d’intégration 1 de l’élément 23, qui est situé environ au centre de la membrane) + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_intermediaire.her b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_intermediaire.her new file mode 100644 index 0000000..f8f703f --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_intermediaire.her @@ -0,0 +1,198 @@ + +############################################################################# +# # +# 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 : 2000 x 1000 ( 10 x 5 elements ) + # geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire. + + nom_maillage tunnel_intermediaire + + noeuds ------------ + 66 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 200 0 + 3 0 400 0 + 4 0 600 0 + 5 0 800 0 + 6 0 1000 0 + 7 200 0 0 + 8 200 200 0 + 9 200 400 0 + 10 200 600 0 + 11 200 800 0 + 12 200 1000 0 + 13 400 0 0 + 14 400 200 0 + 15 400 400 0 + 16 400 600 0 + 17 400 800 0 + 18 400 1000 0 + 19 600 0 0 + 20 600 200 0 + 21 600 400 0 + 22 600 600 0 + 23 600 800 0 + 24 600 1000 0 + 25 800 0 0 + 26 800 200 0 + 27 800 400 0 + 28 800 600 0 + 29 800 800 0 + 30 800 1000 0 + 31 1000 0 0 + 32 1000 200 0 + 33 1000 400 0 + 34 1000 600 0 + 35 1000 800 0 + 36 1000 1000 0 + 37 1200 0 0 + 38 1200 200 0 + 39 1200 400 0 + 40 1200 600 0 + 41 1200 800 0 + 42 1200 1000 0 + 43 1400 0 0 + 44 1400 200 0 + 45 1400 400 0 + 46 1400 600 0 + 47 1400 800 0 + 48 1400 1000 0 + 49 1600 0 0 + 50 1600 200 0 + 51 1600 400 0 + 52 1600 600 0 + 53 1600 800 0 + 54 1600 1000 0 + 55 1800 0 0 + 56 1800 200 0 + 57 1800 400 0 + 58 1800 600 0 + 59 1800 800 0 + 60 1800 1000 0 + 61 2000 0 0 + 62 2000 200 0 + 63 2000 400 0 + 64 2000 600 0 + 65 2000 800 0 + 66 2000 1000 0 + + # references des noeuds + +N_O 1 2 3 4 5 6 +N_S 1 7 13 19 25 31 37 43 49 55 61 +N_E 61 62 63 64 65 66 +N_N 6 12 18 24 30 36 42 48 54 60 66 +N_SO 1 +N_SE 61 +N_NE 66 +N_NO 6 + N_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 + 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 + 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 + 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 + 65 66 + + + # les elements + + elements ---------- + 50 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 QUADRANGLE LINEAIRE 1 7 8 2 + 2 QUADRANGLE LINEAIRE 2 8 9 3 + 3 QUADRANGLE LINEAIRE 3 9 10 4 + 4 QUADRANGLE LINEAIRE 4 10 11 5 + 5 QUADRANGLE LINEAIRE 5 11 12 6 + 6 QUADRANGLE LINEAIRE 7 13 14 8 + 7 QUADRANGLE LINEAIRE 8 14 15 9 + 8 QUADRANGLE LINEAIRE 9 15 16 10 + 9 QUADRANGLE LINEAIRE 10 16 17 11 + 10 QUADRANGLE LINEAIRE 11 17 18 12 + 11 QUADRANGLE LINEAIRE 13 19 20 14 + 12 QUADRANGLE LINEAIRE 14 20 21 15 + 13 QUADRANGLE LINEAIRE 15 21 22 16 + 14 QUADRANGLE LINEAIRE 16 22 23 17 + 15 QUADRANGLE LINEAIRE 17 23 24 18 + 16 QUADRANGLE LINEAIRE 19 25 26 20 + 17 QUADRANGLE LINEAIRE 20 26 27 21 + 18 QUADRANGLE LINEAIRE 21 27 28 22 + 19 QUADRANGLE LINEAIRE 22 28 29 23 + 20 QUADRANGLE LINEAIRE 23 29 30 24 + 21 QUADRANGLE LINEAIRE 25 31 32 26 + 22 QUADRANGLE LINEAIRE 26 32 33 27 + 23 QUADRANGLE LINEAIRE 27 33 34 28 + 24 QUADRANGLE LINEAIRE 28 34 35 29 + 25 QUADRANGLE LINEAIRE 29 35 36 30 + 26 QUADRANGLE LINEAIRE 31 37 38 32 + 27 QUADRANGLE LINEAIRE 32 38 39 33 + 28 QUADRANGLE LINEAIRE 33 39 40 34 + 29 QUADRANGLE LINEAIRE 34 40 41 35 + 30 QUADRANGLE LINEAIRE 35 41 42 36 + 31 QUADRANGLE LINEAIRE 37 43 44 38 + 32 QUADRANGLE LINEAIRE 38 44 45 39 + 33 QUADRANGLE LINEAIRE 39 45 46 40 + 34 QUADRANGLE LINEAIRE 40 46 47 41 + 35 QUADRANGLE LINEAIRE 41 47 48 42 + 36 QUADRANGLE LINEAIRE 43 49 50 44 + 37 QUADRANGLE LINEAIRE 44 50 51 45 + 38 QUADRANGLE LINEAIRE 45 51 52 46 + 39 QUADRANGLE LINEAIRE 46 52 53 47 + 40 QUADRANGLE LINEAIRE 47 53 54 48 + 41 QUADRANGLE LINEAIRE 49 55 56 50 + 42 QUADRANGLE LINEAIRE 50 56 57 51 + 43 QUADRANGLE LINEAIRE 51 57 58 52 + 44 QUADRANGLE LINEAIRE 52 58 59 53 + 45 QUADRANGLE LINEAIRE 53 59 60 54 + 46 QUADRANGLE LINEAIRE 55 61 62 56 + 47 QUADRANGLE LINEAIRE 56 62 63 57 + 48 QUADRANGLE LINEAIRE 57 63 64 58 + 49 QUADRANGLE LINEAIRE 58 64 65 59 + 50 QUADRANGLE LINEAIRE 59 65 66 60 + + + # references des elements + +A_O 1 4 2 4 3 4 4 4 5 4 +A_S 1 1 6 1 11 1 16 1 21 1 26 1 31 1 36 1 41 1 46 1 +A_E 46 2 47 2 48 2 49 2 50 2 +A_N 5 3 10 3 15 3 20 3 25 3 30 3 35 3 40 3 45 3 50 3 +F_to 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 + 12 1 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 21 1 22 1 + 23 1 24 1 25 1 26 1 27 1 28 1 29 1 30 1 31 1 32 1 33 1 + 34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 42 1 43 1 44 1 + 45 1 46 1 47 1 48 1 49 1 50 1 +E_O 1 2 3 4 5 +E_S 1 6 11 16 21 26 31 36 41 46 +E_E 46 47 48 49 50 +E_N 5 10 15 20 25 30 35 40 45 50 +E_SO 1 +E_SE 46 +E_NE 50 +E_NO 5 + E_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 + 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 + 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 + 49 50 + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.CVisu1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.CVisu1 new file mode 100644 index 0000000..22d895e --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.CVisu1 @@ -0,0 +1,114 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.703 # +# 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 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales fin_grandeurs_globales + # ensuite pour chaque maillage:, + # le numero du maillage , + # les infos pour la visualisation eventuelle des torseurs de reaction, + # les infos pour la visualisation eventuelle aux noeud, + # - tout d'abord les ddl principaux (position, deplacements, temperature etc.) + # - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses + # les infos pour la visualisation eventuelle aux elements, + # - tout d'abord les grandeurs generiques (deformations, contraintes etc.) + # - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def + # enfin un mot cle de fin ( fin_grandeurs_maple) + 1 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud fin_list_ddl_noeud + type_sortie_ddl_retenue= 0 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg 23 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements +# informations particuliere dans le cas ou il y a une animation +# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements) +# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent +# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme. + debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier +# un parametre de pilotage du style de sortie + parametre_style_de_sortie 1 +# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc +# (tangent pour les coques, suivant la fibre moyenne pour les element 1D ) + tenseur_en_absolue_ 0 + fin_grandeurs_maple # fin des grandeurs a visualiser au format maple + + # ----------------------------- definition des parametres d'animation: ---------------- + debut_animation # un mot cle de debut de liste (debut_animation) + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # des parametres avec des valeurs: (sur une meme ligne) + cycleInterval 8 # cycleInterval (indique le temps en seconde du cycle de l'animation) + fin_animation # un mot cle de fin + + fin_visualisation_maple + # ================================================================================= + # || fin de la visualisation maple || + # ================================================================================= + + + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.info new file mode 100644 index 0000000..7a88349 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP.info @@ -0,0 +1,139 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2015 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.695 + + +####################################################################### +# fichier de commande # + + +# --- commentaire sur l'objet du travail ??? + + + +#---------------------------------------------------- +# definition de la dimension de l'espace de travail | +#---------------------------------------------------- +dimension 3 + + +#--------------------------------------------------------------- +# definition facultative du niveau d'impression (entre 0 et 10)| +#--------------------------------------------------------------- +niveau_commentaire 4 + + +#------------------------------- +# definition du type de calcul | +#------------------------------- +TYPE_DE_CALCUL + +#non_dynamique + +dynamique_relaxation_dynam #avec plus visualisation + + + +PARA_TYPE_DE_CALCUL + + +< KDR.algo + + + + +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- + +# -- def maillage +< tunnel_intermediaire.her + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + +choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +E_to elas_iso + +materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ +elas_iso LOI_CONTRAINTES_PLANES + # ....... loi de comportement LoiContraintesPlanes ........ + NEWTON_LOCAL # paramètres par défaut + ISOELAS + # ....... loi de comportement isoelastique 3D ........ + # module d'young : coefficient de poisson + 125 0.4 +fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane + +#--------------------------------------- fin def des lois de comportement ------------------------ + + +# --- divers stockages (1) ------- +epaisseurs #--------# +E_to 25e-3 +masse_volumique #--------# +E_to 0.9e-9 + +charges #------------# +F_to PRESSION -0.390e-3 + +blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +N_O UZ UX UY +N_E UZ UY + + +controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +SAUVEGARDE 0 +ITERATIONS 300000 +PRECISION 0.005 +MULTIPLICATEUR 1 +DELTAt 0.1 +TEMPSFIN 1 +RESTART 0 +DELTAtMAXI 0.1 +NORME E_cinetique/E_statique_ET_ResSurReact + + +para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +FREQUENCE_AFFICHAGE_ITERATION 100 +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + +# ------------------------------------------------------------------------------------ +resultats pas_de_sortie_finale_ +COPIE 0 +# + + +_fin_point_info_ diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP_princ.maple.ref1 new file mode 100644 index 0000000..743dfe2 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_ISOELAS_CP/tunnel_isoelast_CP_princ.maple.ref1 @@ -0,0 +1,87 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.703 # +# 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 -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# element_23 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII +# +# +#==================================================================== +#|| 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 4.882353771716e+02 4.424495084397e+02 6.181246032206e+02 7.970860513281e-03 -3.211642487148e-03 9.949699737870e-01 -3.467794108425e-03 +2.000000000000e-01 4.959653502342e+02 4.426378607101e+02 6.223508900465e+02 1.624443498337e-02 -6.521727255003e-03 2.029127668444e+00 -3.565388553768e-03 +3.000000000000e-01 4.972811916717e+02 4.428312561206e+02 6.291126103184e+02 2.456031467167e-02 -9.954539800997e-03 3.062275759675e+00 -1.940767538481e-02 +4.000000000000e-01 5.019141162088e+02 4.430175911776e+02 6.340647017512e+02 3.297754273946e-02 -1.329514671476e-02 4.115993687453e+00 -1.549644238806e-02 +5.000000000000e-01 5.035063352809e+02 4.432172415557e+02 6.398823579697e+02 4.153055631016e-02 -1.691100713672e-02 5.173533896005e+00 -4.446284032331e-02 +6.000000000000e-01 5.056522149909e+02 4.434135298018e+02 6.479514578033e+02 5.035724872586e-02 -2.050339030892e-02 6.273197256534e+00 -5.364551442781e-02 +7.000000000000e-01 5.154781672705e+02 4.436434523049e+02 6.515363839418e+02 6.045772732791e-02 -2.475911482453e-02 7.522927180923e+00 -8.571943983130e-02 +8.000000000000e-01 5.145020074068e+02 4.438367404669e+02 6.616539343554e+02 6.950550761708e-02 -2.837728697643e-02 8.653956103341e+00 -8.557912822277e-02 +9.000000000000e-01 5.271557119695e+02 4.440942687122e+02 6.678559043492e+02 8.054214347634e-02 -3.325641313039e-02 1.000588772921e+01 -1.546976821082e-01 +1.000000000000e+00 5.356044582301e+02 4.442871686731e+02 6.730004909504e+02 9.044645294194e-02 -3.695551097365e-02 1.125955889461e+01 -1.156164550013e-01 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/KDR.algo b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/KDR.algo new file mode 100644 index 0000000..021312e --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/KDR.algo @@ -0,0 +1,20 @@ +# ----- KDR +#typeCalRelaxation= 1 lambda= 1.6 type_calcul_mass= 2 option_recalcul_mass= 1 +typeCalRelaxation= 1 lambda= 2. type_calcul_mass= 2 option_recalcul_mass= 1 +parametre_calcul_de_la_masse_ casMass_relax= 3 + +# ----- Amortissement cinetique +avec_amortissement_cinetique_ + max_nb_decroit_pourRelaxDyn_ 1 + coef_arret_pourRelaxDyn_ 0. + coef_redemarrage_pourRelaxDyn_ 0. + max_deltaX_pourRelaxDyn_ 0 #0.005 +nb_max_dX_OK_pourRelaxDyn_ 10 #1 +nb_deb_testfin_pourRelaxDyn_ 250 #1 +nb_deb_test_amort_cinetique_ 250 #1 +fi_parametre_amortissement_cinetique_ + + + +ARRET_A_EQUILIBRE_STATIQUE_ 2 + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/README.txt b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/README.txt new file mode 100644 index 0000000..334584c --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/README.txt @@ -0,0 +1,36 @@ +------------------------------------------------------------ +But du test +------------------------------------------------------------ +Le but du calcul est de trouver la forme d’équilibre d’un tunnel gonflé +dont le matériau suit une loi de comportement isotrope hyperélastique +ISOHYOER3DORGEAS1 + LOI_CONTRAINTES_PLANES. + +------------------------------------------------------------ +Description du calcul +------------------------------------------------------------ +Le calcul s’effectue avec la méthode de relaxation dynamique. +Les paramètres de contrôle de la méthode de relaxation dynamique sont définis dans KDR.algo + +La loi de comportement testée est: +- ISOHYPER3DORGEAS1 +- LOI_CONTRAINTES_PLANES avec prise en compte de la condition de contraintes planes par la méthode de NEWTON_LOCAL +Les paramètres utilisés pour définir ces lois sont des paramètres cohérents avec des tests réalisés sur des membranes: +K = 624; Q0s = 2; mu01 = 30; mu02 = 14.6; mu03 = 90; alpha1 = 1e-3; alpha2 = 1e-3; Q0e =0.12 + +La structure est une membrane de dimension 1000x2000mm, +discrétisée en 1x12 éléments membranes (interpolation linéaire - 8noeuds). + +------------------------------------------------------------ +Mots-cles +------------------------------------------------------------ +tunnel +hyperelastique +contraintes planes +relaxation dynamique +isohyper3dorgeas1 + +------------------------------------------------------------ +Grandeurs de comparaison +------------------------------------------------------------ +Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII +(point d’intégration 1 de l’élément 6, qui est situé environ au centre de la membrane) \ No newline at end of file diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_12x1.her b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_12x1.her new file mode 100644 index 0000000..aa26a49 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_12x1.her @@ -0,0 +1,115 @@ + +############################################################################# +# # +# 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 : 2000 x 1000 ( 12 x 1 elements ) + # geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire. + + nom_maillage tunnel_12x1 + + noeuds ------------ + 26 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 1000 0 + 3 166.6666666666667 0 0 + 4 166.6666666666667 1000 0 + 5 333.3333333333333 0 0 + 6 333.3333333333333 1000 0 + 7 500 0 0 + 8 500 1000 0 + 9 666.6666666666666 0 0 + 10 666.6666666666666 1000 0 + 11 833.3333333333333 0 0 + 12 833.3333333333333 1000 0 + 13 999.9999999999999 0 0 + 14 999.9999999999999 1000 0 + 15 1166.666666666667 0 0 + 16 1166.666666666667 1000 0 + 17 1333.333333333333 0 0 + 18 1333.333333333333 1000 0 + 19 1500 0 0 + 20 1500 1000 0 + 21 1666.666666666667 0 0 + 22 1666.666666666667 1000 0 + 23 1833.333333333333 0 0 + 24 1833.333333333333 1000 0 + 25 2000 0 0 + 26 2000 1000 0 + + # references des noeuds + +N_O 1 2 +N_S 1 3 5 7 9 11 13 15 17 19 21 23 25 +N_E 25 26 +N_N 2 4 6 8 10 12 14 16 18 20 22 24 26 +N_SO 1 +N_SE 25 +N_NE 26 +N_NO 2 + N_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 + 17 18 19 20 21 22 23 24 25 26 + + + # les elements + + elements ---------- + 12 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 QUADRANGLE LINEAIRE 1 3 4 2 + 2 QUADRANGLE LINEAIRE 3 5 6 4 + 3 QUADRANGLE LINEAIRE 5 7 8 6 + 4 QUADRANGLE LINEAIRE 7 9 10 8 + 5 QUADRANGLE LINEAIRE 9 11 12 10 + 6 QUADRANGLE LINEAIRE 11 13 14 12 + 7 QUADRANGLE LINEAIRE 13 15 16 14 + 8 QUADRANGLE LINEAIRE 15 17 18 16 + 9 QUADRANGLE LINEAIRE 17 19 20 18 + 10 QUADRANGLE LINEAIRE 19 21 22 20 + 11 QUADRANGLE LINEAIRE 21 23 24 22 + 12 QUADRANGLE LINEAIRE 23 25 26 24 + + + # references des elements + +A_O 1 4 +A_S 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 + 12 1 +A_E 12 2 +A_N 1 3 2 3 3 3 4 3 5 3 6 3 7 3 8 3 9 3 10 3 11 3 + 12 3 +F_to 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 11 1 + 12 1 +E_O 1 +E_S 1 2 3 4 5 6 7 8 9 10 11 + 12 +E_E 12 +E_N 1 2 3 4 5 6 7 8 9 10 11 + 12 +E_SO 1 +E_SE 12 +E_NE 12 +E_NO 1 + E_to 1 2 3 4 5 6 7 8 9 10 11 12 + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.CVisu1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.CVisu1 new file mode 100644 index 0000000..b96ef40 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.CVisu1 @@ -0,0 +1,113 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.703 # +# 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 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales fin_grandeurs_globales + # ensuite pour chaque maillage:, + # le numero du maillage , + # les infos pour la visualisation eventuelle des torseurs de reaction, + # les infos pour la visualisation eventuelle aux noeud, + # - tout d'abord les ddl principaux (position, deplacements, temperature etc.) + # - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses + # les infos pour la visualisation eventuelle aux elements, + # - tout d'abord les grandeurs generiques (deformations, contraintes etc.) + # - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def + # enfin un mot cle de fin ( fin_grandeurs_maple) + 1 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud fin_list_ddl_noeud + type_sortie_ddl_retenue= 0 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg 6 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements +# informations particuliere dans le cas ou il y a une animation +# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements) +# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent +# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme. + debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier +# un parametre de pilotage du style de sortie + parametre_style_de_sortie 1 +# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc +# (tangent pour les coques, suivant la fibre moyenne pour les element 1D ) + tenseur_en_absolue_ 0 + fin_grandeurs_maple # fin des grandeurs a visualiser au format maple + + # ----------------------------- definition des parametres d'animation: ---------------- + debut_animation # un mot cle de debut de liste (debut_animation) + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # des parametres avec des valeurs: (sur une meme ligne) + cycleInterval 8 # cycleInterval (indique le temps en seconde du cycle de l'animation) + fin_animation # un mot cle de fin + + fin_visualisation_maple + # ================================================================================= + # || fin de la visualisation maple || + # ================================================================================= + + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.info new file mode 100644 index 0000000..e7fdfea --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast.info @@ -0,0 +1,145 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2015 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.695 + + +####################################################################### +# fichier de commande # + + +# --- commentaire sur l'objet du travail ??? + + + +#---------------------------------------------------- +# definition de la dimension de l'espace de travail | +#---------------------------------------------------- +dimension 3 + + +#--------------------------------------------------------------- +# definition facultative du niveau d'impression (entre 0 et 10)| +#--------------------------------------------------------------- +niveau_commentaire 4 + + +#------------------------------- +# definition du type de calcul | +#------------------------------- +TYPE_DE_CALCUL + +#non_dynamique + +dynamique_relaxation_dynam #avec plus visualisation + +PARA_TYPE_DE_CALCUL + +< KDR.algo +#mode_debug_= 5 + +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- + +# -- def maillage +< tunnel_12x1.her + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + +choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +E_to orgeas2D + +materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ +orgeas2D LOI_CONTRAINTES_PLANES + + # ....... loi de comportement LoiContraintesPlanes ........ + NEWTON_LOCAL + ISOHYPER3DORGEAS1 + +# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........ +#----------------------------------------------------------------------------------------- +# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e | +#----------------------------------------------------------------------------------------- +# parametres coherents avec des membranes, ajustes sur resultats experimentaux +624 2 30 14.6 90 1e-3 1e-3 0.12 + + +fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane + +#--------------------------------------- fin def des lois de comportement ------------------------ + + +# --- divers stockages (1) ------- +epaisseurs #--------# +E_to 25e-3 +masse_volumique #--------# +E_to 0.9e-9 + +charges #------------# +F_to PRESSION -0.390e-3 + +blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +N_O UZ UX UY +N_E UZ UY + +controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +SAUVEGARDE 1#0 +ITERATIONS 300000 +PRECISION 0.005 +MULTIPLICATEUR 1 +DELTAt 0.1 +TEMPSFIN 1 +RESTART 0 +DELTAtMAXI 0.1 +NORME E_cinetique/E_statique_ET_ResSurReact +MAXINCRE 1 + +para_pilotage_equi_global #------ parametre pour le jacobien + +CAS_JACOBIEN_NEGATIF 2 + +para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +FREQUENCE_AFFICHAGE_INCREMENT 1 +FREQUENCE_AFFICHAGE_ITERATION 100 +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + +# ------------------------------------------------------------------------------------ +resultats pas_de_sortie_finale_ +COPIE 0 +# + + +_fin_point_info_ diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast_princ.maple.ref1 new file mode 100644 index 0000000..59d0b41 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_L_tunnel_isohyper3dorgeas1_CP/tunnel_hyperelast_princ.maple.ref1 @@ -0,0 +1,78 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.710 # +# 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 -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# element_6 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII +# +# +#==================================================================== +#|| 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 5.155709198288e+02 2.122580326756e+02 6.250076281113e+02 8.047703236407e-03 -3.248351140868e-03 1.001048740192e+00 2.258452053693e-03 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/README.txt b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/README.txt new file mode 100644 index 0000000..342fec7 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/README.txt @@ -0,0 +1,36 @@ +------------------------------------------------------------ +But du test +------------------------------------------------------------ +Ce programme teste la loi de comportement isotrope hyperélastique +ISOHYPER3DORGEAS1 sur une menbrane carrée en traction simple +en utilisant une condition de contrainte plane appliquée à la loi 3D. + +------------------------------------------------------------ +Description du calcul +------------------------------------------------------------ +Le calcul est de type statique. + +La loi de comportement testée est la loi ISOHYPER3DORGEAS1 (loi 3D). +Les paramètres définissant cette loi sont des paramètres ajustés sur ses essais réalisés sur des membranes: +K= 624; Q0s= 2; mu01=30; mu02= 14.6; mu03= 90; alpha01= 1e-3; alpha02= 1e-3; Q0e= 0.12 +Nous lui associons la loi LOI_CONTRAINTES_PLANES afin de faire le calcul en 2D prenant en compte +l’hypothèse des contraintes planes (par la méthode de Newton). + +La structure étudiée est une membrane carrée de dimensions 100*100 mm^2, discrétisée +avec un seul élément de type membrane (interpolation linéaire). + +L'essai simulé est un essai de traction simple dans la direction 11. + +------------------------------------------------------------ +Mots-cles +------------------------------------------------------------ +quadrangle +isohyper3dorgeas1 +contraintes planes + + +------------------------------------------------------------ +Grandeurs de comparaison +------------------------------------------------------------ +SIG11 SIG22 SIG12 +EPS11 EPS22 EPS12 \ No newline at end of file diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle1.her b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle1.her new file mode 100644 index 0000000..3ff0bea --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle1.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 : 100 x 100 ( 1 x 1 elements ) + # geometrie rectangulaire, decoupage rectangulaire, interpolation lineaire. + + nom_maillage quadrangle1 + + noeuds ------------ + 4 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 100 0 + 3 100 0 0 + 4 100 100 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 QUADRANGLE 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/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.CVisu1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.CVisu1 new file mode 100644 index 0000000..6adf36a --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.CVisu1 @@ -0,0 +1,112 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.699 # +# 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 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales fin_grandeurs_globales + # ensuite pour chaque maillage:, + # le numero du maillage , + # les infos pour la visualisation eventuelle des torseurs de reaction, + # les infos pour la visualisation eventuelle aux noeud, + # - tout d'abord les ddl principaux (position, deplacements, temperature etc.) + # - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses + # les infos pour la visualisation eventuelle aux elements, + # - tout d'abord les grandeurs generiques (deformations, contraintes etc.) + # - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def + # enfin un mot cle de fin ( fin_grandeurs_maple) + 1 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud fin_list_ddl_noeud + type_sortie_ddl_retenue= 0 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element SIG11 SIG22 SIG12 EPS11 EPS22 EPS12 fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements +# informations particuliere dans le cas ou il y a une animation +# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements) +# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent +# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme. + debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier +# un parametre de pilotage du style de sortie + parametre_style_de_sortie 1 +# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc +# (tangent pour les coques, suivant la fibre moyenne pour les element 1D ) + tenseur_en_absolue_ 0 + fin_grandeurs_maple # fin des grandeurs a visualiser au format maple + + # ----------------------------- definition des parametres d'animation: ---------------- + debut_animation # un mot cle de debut de liste (debut_animation) + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # des parametres avec des valeurs: (sur une meme ligne) + cycleInterval 8 # cycleInterval (indique le temps en seconde du cycle de l'animation) + fin_animation # un mot cle de fin + + fin_visualisation_maple + # ================================================================================= + # || fin de la visualisation maple || + # ================================================================================= + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.info new file mode 100644 index 0000000..5988570 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction.info @@ -0,0 +1,149 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2015 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.697 + + +####################################################################### +# fichier de commande # + + +# --- commentaire sur l'objet du travail ??? + + + +#---------------------------------------------------- +# 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 + < quadrangle1.her + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + +choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +E_to orgeas2D + +materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ +orgeas2D LOI_CONTRAINTES_PLANES + + +NEWTON_LOCAL + +ISOHYPER3DORGEAS1 + +# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........ +#----------------------------------------------------------------------------------------- +# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e | +#----------------------------------------------------------------------------------------- +# parametres coherents avec des membranes, ajustes sur resultats experimentaux +624 2 30 14.6 90 1e-3 1e-3 0.12 + +fin_loi_contrainte_plane # ----- fin de Loi de contrainte plane + +#--------------------------------------- fin def des lois de comportement ------------------------ + + +# --- divers stockages (1) ------- + epaisseurs #--------# +E_to 25e-3 + masse_volumique #--------# +E_to 0.9e-9 + + charges #------------# +# ici pas de chargement a definir + + + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- + +# conditions limites pour un essai de traction selon l axe des X +N_N UZ +N_S UZ +N_E UZ +N_O UZ UX +N_E 'UX=20' +N_SO UY +N_SE UY + + +controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +SAUVEGARDE 0 +ITERATIONS 300000 +PRECISION 0.005 +MULTIPLICATEUR 1 +DELTAt 0.005 +TEMPSFIN 1 +RESTART 0 +DELTAtMAXI 0.1 + + + +para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + + + +# ------------------------------------------------------------------------------------ + resultats pas_de_sortie_finale_ + COPIE 0 +# + + + _fin_point_info_ + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction_princ.maple.ref1 new file mode 100644 index 0000000..1748acb --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/CP_venant_du_3D/Test_R_quadrangle_isohyper3dorgeas1_CP_traction/quadrangle_traction_princ.maple.ref1 @@ -0,0 +1,113 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.703 # +# 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 -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 6 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# element_1 pt_integ_1: [2]X [3]Y [4]Z [5] SIG11 [6] SIG22 [7] SIG12 [8] EPS11 [9] EPS22 [10] EPS12 +# +# +#==================================================================== +#|| fin du recapitulatif des differentes grandeurs || +#==================================================================== + + # ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant + # a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun + # correspondant a un maillage. On trouve pour chaque enregistrement successivement : + # s'il y a des grandeurs globales: le temps puis les grandeurs globales, + # puis s'il y a des torseurs de reaction : + # de nouveau le temps, les composantes de la resultante puis les composantes du moments + # donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels + # puis s'il y a des grandeurs aux noeuds: de nouveau le temps + # les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud + # puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds + # puis s'il y a des grandeur aux elements: + # le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques) + # suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration + # correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes + # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs + # etc. pour tous les points d'integration - element + + 5.000000000000e-03 2.115361902706e+01 2.112406169668e+01 0.000000000000e+00 1.245967678291e-01 9.603156326874e-05 0.000000000000e+00 9.985019975028e-04 -3.989064216068e-04 0.000000000000e+00 +1.000000000000e-02 2.117475151360e+01 2.111563685285e+01 0.000000000000e+00 2.487309504414e-01 2.189700783719e-04 0.000000000000e+00 1.994015960096e-03 -7.982904127122e-04 0.000000000000e+00 +1.500000000000e-02 2.119588400014e+01 2.110721200902e+01 0.000000000000e+00 3.723318674740e-01 2.947116133150e-04 0.000000000000e+00 2.986553798226e-03 -1.198152735947e-03 0.000000000000e+00 +2.000000000000e-02 2.121701648668e+01 2.109878716518e+01 0.000000000000e+00 4.954025576742e-01 3.229168731565e-04 0.000000000000e+00 3.976127363058e-03 -1.598494155465e-03 0.000000000000e+00 +2.707106781187e-02 2.124690233575e+01 2.108687263677e+01 0.000000000000e+00 6.684817550666e-01 2.114892521729e-04 0.000000000000e+00 5.370558286617e-03 -2.165481776918e-03 0.000000000000e+00 +3.414213562373e-02 2.127678818483e+01 2.107495810836e+01 0.000000000000e+00 8.405845945444e-01 7.308843204531e-05 0.000000000000e+00 6.759117391982e-03 -2.733431293449e-03 0.000000000000e+00 +4.121320343560e-02 2.130667403390e+01 2.106304357995e+01 0.000000000000e+00 1.011647903613e+00 -1.642749678023e-04 0.000000000000e+00 8.141837600510e-03 -3.302344882098e-03 0.000000000000e+00 +4.828427124746e-02 2.133655988297e+01 2.105112905154e+01 0.000000000000e+00 1.181679194540e+00 -5.017109596818e-04 0.000000000000e+00 9.518751603152e-03 -3.872224726065e-03 0.000000000000e+00 +5.828427124746e-02 2.137882485605e+01 2.103427936387e+01 0.000000000000e+00 1.420256734858e+00 -1.286104655104e-03 0.000000000000e+00 1.145615326926e-02 -4.679810286139e-03 0.000000000000e+00 +6.828427124746e-02 2.142108982913e+01 2.101742967620e+01 0.000000000000e+00 1.656938955111e+00 -2.140067561486e-03 0.000000000000e+00 1.338209845641e-02 -5.489338950370e-03 0.000000000000e+00 +7.828427124746e-02 2.146335480221e+01 2.100057998853e+01 0.000000000000e+00 1.891602939365e+00 -3.203540448401e-03 0.000000000000e+00 1.529667731467e-02 -6.300816957409e-03 0.000000000000e+00 +8.828427124746e-02 2.150561977529e+01 2.098373030087e+01 0.000000000000e+00 2.124258568724e+00 -4.478167388988e-03 0.000000000000e+00 1.719997910912e-02 -7.114250570967e-03 0.000000000000e+00 +1.024264068712e-01 2.156539147344e+01 2.095990124404e+01 0.000000000000e+00 2.449586900641e+00 -6.901881440541e-03 0.000000000000e+00 1.987257471975e-02 -8.267970071699e-03 0.000000000000e+00 +1.165685424949e-01 2.162516317158e+01 2.093607218722e+01 0.000000000000e+00 2.771063439458e+00 -9.469614229263e-03 0.000000000000e+00 2.252303990196e-02 -9.425631241595e-03 0.000000000000e+00 +1.307106781187e-01 2.168493486973e+01 2.091224313040e+01 0.000000000000e+00 3.087874889780e+00 -1.236182894675e-02 0.000000000000e+00 2.515161831774e-02 -1.058725205668e-02 0.000000000000e+00 +1.448528137424e-01 2.174470656787e+01 2.088841407358e+01 0.000000000000e+00 3.394533227239e+00 -1.454248524462e-02 0.000000000000e+00 2.775855028480e-02 -1.175285059557e-02 0.000000000000e+00 +1.648528137424e-01 2.182923651403e+01 2.085471469824e+01 0.000000000000e+00 3.635970900778e+00 8.900885125220e-03 0.000000000000e+00 3.140881587810e-02 -1.340808160172e-02 0.000000000000e+00 +1.848528137424e-01 2.191376646020e+01 2.081559778866e+01 0.000000000000e+00 3.775275349391e+00 -8.523720758650e-05 0.000000000000e+00 3.501692149978e-02 -1.533949935420e-02 0.000000000000e+00 +2.048528137424e-01 2.199829640636e+01 2.077648087908e+01 0.000000000000e+00 3.910439062755e+00 -9.243746091742e-03 0.000000000000e+00 3.858351391342e-02 -1.728183650063e-02 0.000000000000e+00 +2.248528137424e-01 2.208282635252e+01 2.073736396950e+01 0.000000000000e+00 4.042992206342e+00 -1.891739832852e-02 0.000000000000e+00 4.210922752774e-02 -1.923517550773e-02 0.000000000000e+00 +2.531370849898e-01 2.220236974881e+01 2.068765464959e+01 0.000000000000e+00 4.241640540628e+00 -1.196950925605e-04 0.000000000000e+00 4.702676176270e-02 -2.173346133597e-02 0.000000000000e+00 +2.814213562373e-01 2.232191314510e+01 2.063794532968e+01 0.000000000000e+00 4.435204617174e+00 1.547688925023e-02 0.000000000000e+00 5.186550124481e-02 -2.424982128783e-02 0.000000000000e+00 +3.097056274848e-01 2.244145654139e+01 2.058332849849e+01 0.000000000000e+00 4.611002260797e+00 1.945234469155e-04 0.000000000000e+00 5.662712042139e-02 -2.703565370811e-02 0.000000000000e+00 +3.379898987322e-01 2.256099993768e+01 2.052871166730e+01 0.000000000000e+00 4.782487246584e+00 -1.663624311373e-02 0.000000000000e+00 6.131324949586e-02 -2.984375089417e-02 0.000000000000e+00 +3.779898987322e-01 2.273005983000e+01 2.045898895472e+01 0.000000000000e+00 5.039077829464e+00 -6.024431578780e-06 0.000000000000e+00 6.781464266963e-02 -3.346124066146e-02 0.000000000000e+00 +4.179898987322e-01 2.289911972233e+01 2.038926624214e+01 0.000000000000e+00 5.285582564498e+00 1.055045129288e-02 0.000000000000e+00 7.417257236372e-02 -3.711590472384e-02 0.000000000000e+00 +4.579898987322e-01 2.306817961465e+01 2.031954352955e+01 0.000000000000e+00 5.524282424941e+00 1.692089993712e-02 0.000000000000e+00 8.039122873540e-02 -4.080825418225e-02 0.000000000000e+00 +4.979898987322e-01 2.323723950697e+01 2.024982081697e+01 0.000000000000e+00 5.756098513475e+00 1.931845591194e-02 0.000000000000e+00 8.647465007341e-02 -4.453880895168e-02 0.000000000000e+00 +5.545584412272e-01 2.347632629955e+01 2.015121801123e+01 0.000000000000e+00 6.090949936850e+00 1.429757509892e-02 0.000000000000e+00 9.485458158242e-02 -4.988086006574e-02 0.000000000000e+00 +6.111269837221e-01 2.371541309213e+01 2.008095327997e+01 0.000000000000e+00 7.892359737027e+00 -1.364063141236e-03 0.000000000000e+00 1.029823461463e-01 -5.373573961986e-02 0.000000000000e+00 +6.676955262170e-01 2.395449988471e+01 2.001957237612e+01 0.000000000000e+00 9.920994044886e+00 -1.486750209981e-04 0.000000000000e+00 1.108679607376e-01 -5.713650214774e-02 0.000000000000e+00 +7.242640687119e-01 2.419358667729e+01 1.995819147227e+01 0.000000000000e+00 1.187929967617e+01 6.835735725519e-04 0.000000000000e+00 1.185209498365e-01 -6.056868977309e-02 0.000000000000e+00 +8.042640687119e-01 2.453170646193e+01 1.987138576558e+01 0.000000000000e+00 1.453402508573e+01 1.752436737719e-04 0.000000000000e+00 1.289643090895e-01 -6.547693781182e-02 0.000000000000e+00 +8.842640687119e-01 2.486982624658e+01 1.978458005889e+01 0.000000000000e+00 1.706330747541e+01 -4.133812178877e-03 0.000000000000e+00 1.389846198943e-01 -7.044993293332e-02 0.000000000000e+00 +9.642640687119e-01 2.520794603123e+01 1.969777435220e+01 0.000000000000e+00 1.947653199461e+01 -1.444581530112e-02 0.000000000000e+00 1.486044270979e-01 -7.548881897902e-02 0.000000000000e+00 +1.000000000000e+00 2.535898384862e+01 1.965899831758e+01 0.000000000000e+00 2.051931548721e+01 -2.156859637703e-02 0.000000000000e+00 1.527777777778e-01 -7.776128287454e-02 0.000000000000e+00 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/.DS_Store b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/.DS_Store new file mode 100644 index 0000000000000000000000000000000000000000..5008ddfcf53c02e82d7eee2e57c38e5672ef89f6 GIT binary patch literal 6148 zcmeH~Jr2S!425mzP>H1@V-^m;4Wg<&0T*E43hX&L&p$$qDprKhvt+--jT7}7np#A3 zem<@ulZcFPQ@L2!n>{z**++&mCkOWA81W14cNZlEfg7;MkzE(HCqgga^y>{tEnwC%0;vJ&^%eQ zLs35+`xjp>T0>>>>> 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 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales fin_grandeurs_globales + # ensuite pour chaque maillage:, + # le numero du maillage , + # les infos pour la visualisation eventuelle des torseurs de reaction, + # les infos pour la visualisation eventuelle aux noeud, + # - tout d'abord les ddl principaux (position, deplacements, temperature etc.) + # - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses + # les infos pour la visualisation eventuelle aux elements, + # - tout d'abord les grandeurs generiques (deformations, contraintes etc.) + # - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def + # enfin un mot cle de fin ( fin_grandeurs_maple) + 1 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud fin_list_ddl_noeud + type_sortie_ddl_retenue= 0 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg 116 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element Def_principaleI Def_principaleII Sigma_principaleI Sigma_principaleII fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements +# informations particuliere dans le cas ou il y a une animation +# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements) +# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent +# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme. + debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier +# un parametre de pilotage du style de sortie + parametre_style_de_sortie 1 +# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc +# (tangent pour les coques, suivant la fibre moyenne pour les element 1D ) + tenseur_en_absolue_ 0 + fin_grandeurs_maple # fin des grandeurs a visualiser au format maple + + # ----------------------------- definition des parametres d'animation: ---------------- + debut_animation # un mot cle de debut de liste (debut_animation) + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # des parametres avec des valeurs: (sur une meme ligne) + cycleInterval 8 # cycleInterval (indique le temps en seconde du cycle de l'animation) + fin_animation # un mot cle de fin + + fin_visualisation_maple + # ================================================================================= + # || fin de la visualisation maple || + # ================================================================================= + + + + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast.info new file mode 100644 index 0000000..2cd3386 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast.info @@ -0,0 +1,136 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2015 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.695 + + +####################################################################### +# fichier de commande # + + +# --- commentaire sur l'objet du travail ??? + + + +#---------------------------------------------------- +# definition de la dimension de l'espace de travail | +#---------------------------------------------------- + dimension 3 + + +#--------------------------------------------------------------- +# definition facultative du niveau d'impression (entre 0 et 10)| +#--------------------------------------------------------------- + niveau_commentaire 4 + + +#------------------------------- +# definition du type de calcul | +#------------------------------- + TYPE_DE_CALCUL + +#non_dynamique + +dynamique_relaxation_dynam + + +PARA_TYPE_DE_CALCUL +#TYPE 1 + +< KDR.algo + + + + +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- + + # -- def maillage + < tunnel_bis.her + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +E_tout membrane + + materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ + membrane ISOELAS2D_C +# ....... loi de comportement isoelastique 3D ........ +# module d'young : coefficient de poisson + 125 0.4 + +# -- definition du type de deformation (par defaut: DEFORMATION_STANDART) -- + type_de_deformation DEFORMATION_STANDART +#--------------------------------------- fin def des lois de comportement ------------------------ + +# --- divers stockages (1) ------- + epaisseurs #--------# +E_tout 25e-3 + masse_volumique #--------# +E_tout 0.9e-9 + + charges #------------# +F_to PRESSION -0.390e-3 + + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +N_O UZ UX UY +N_E UZ UY + + + controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- + SAUVEGARDE 0 + ITERATIONS 300000 + PRECISION 0.005 + MULTIPLICATEUR 1 + DELTAt 0.1 + TEMPSFIN 1 + DELTAtMAXI 0.1 + NORME E_cinetique/E_statique_ET_ResSurReact + + + + + para_affichage #------------ +#FREQUENCE_AFFICHAGE_INCREMENT 1 +#FREQUENCE_AFFICHAGE_ITERATION 100 +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + +# ------------------------------------------------------------------------------------ + resultats pas_de_sortie_finale_ + COPIE 0 +# + + + _fin_point_info_ + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast_princ.maple.ref1 new file mode 100644 index 0000000..89b1dd8 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/ISOELAS2D_C/Test_L_tunnel_ISOELAS2D_C/tunnel_isoelast_princ.maple.ref1 @@ -0,0 +1,87 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.703 # +# 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 -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 4 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# element_116 pt_integ_1: [2]X [3]Y [4]Z [5] Def_principaleI [6] Def_principaleII [7] Sigma_principaleI [8] Sigma_principaleII +# +# +#==================================================================== +#|| 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 7.730408920278e+02 5.210652559614e+02 6.245103853949e+02 8.033427076661e-03 -3.196644897709e-03 1.005173975830e+00 2.488978118371e-03 +2.000000000000e-01 7.681447945156e+02 5.209963438606e+02 6.345274366115e+02 1.612793207041e-02 -6.505133685726e-03 2.012779552994e+00 -8.029889518028e-03 +3.000000000000e-01 7.809829180392e+02 5.209274263835e+02 6.366382254352e+02 2.455888703561e-02 -9.846662840037e-03 3.068485401726e+00 -3.438694314249e-03 +4.000000000000e-01 7.845634757231e+02 5.208549836996e+02 6.418662485518e+02 3.316246027256e-02 -1.339468098880e-02 4.137587481702e+00 -1.930013091984e-02 +5.000000000000e-01 7.962474522385e+02 5.207802318866e+02 6.475992468084e+02 4.199020586806e-02 -1.709493636559e-02 5.230986803842e+00 -4.447232416248e-02 +6.000000000000e-01 8.016561095824e+02 5.207021944225e+02 6.546612524038e+02 5.127867716650e-02 -2.100063340888e-02 6.380717827819e+00 -7.279204498310e-02 +7.000000000000e-01 8.138162142176e+02 5.206263072773e+02 6.607739699857e+02 6.088097138966e-02 -2.484116374907e-02 7.581027662207e+00 -7.273440375093e-02 +8.000000000000e-01 8.212930601219e+02 5.205427028046e+02 6.707023601256e+02 7.136238854896e-02 -2.912196824722e-02 8.885952566976e+00 -8.586500411190e-02 +9.000000000000e-01 8.349537584645e+02 5.204682193168e+02 6.789366509884e+02 8.112411765326e-02 -3.298031591607e-02 1.010892727483e+01 -7.896857957788e-02 +1.000000000000e+00 8.429195604637e+02 5.203823488565e+02 6.882414598862e+02 9.263139024120e-02 -3.748051654705e-02 1.155344994381e+01 -6.368459085936e-02 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/.DS_Store b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/.DS_Store new file mode 100644 index 0000000000000000000000000000000000000000..5008ddfcf53c02e82d7eee2e57c38e5672ef89f6 GIT binary patch literal 6148 zcmeH~Jr2S!425mzP>H1@V-^m;4Wg<&0T*E43hX&L&p$$qDprKhvt+--jT7}7np#A3 zem<@ulZcFPQ@L2!n>{z**++&mCkOWA81W14cNZlEfg7;MkzE(HCqgga^y>{tEnwC%0;vJ&^%eQ zLs35+`xjp>T0>>>>> le mot cle: + # permet au programme de se positionner au debut du fichier, il est indispensable + + + # ================================================================================= + # || ***** demande d'une visualisation maple: ***** || + # ================================================================================= + # un mot cle de debut (debut_visualisation_maple) + # un mot cle de fin ( fin_visualisation_maple) + # la seule presence de ces deux mots cle suffit a activer la visualisation maple + debut_visualisation_maple + + # ----------------------------- definition de la liste des increments a balayer: ---------------- + debut_list_increment # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments) + # un mot cle de fin de liste ( fin_list_increment) + fin_list_increment + + # ----------------------------- choix des maillages a visualiser: ---------------- + # la liste est facultative, par defaut la visualisation concerne le premier maillage + debut_choix_maillage # un mot cle de debut, + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entiers , puis , sur une meme ligne + 1 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales fin_grandeurs_globales + # ensuite pour chaque maillage:, + # le numero du maillage , + # les infos pour la visualisation eventuelle des torseurs de reaction, + # les infos pour la visualisation eventuelle aux noeud, + # - tout d'abord les ddl principaux (position, deplacements, temperature etc.) + # - puis les ddl etendus et particulieres qui peuvent representer des grandeurs diverses + # les infos pour la visualisation eventuelle aux elements, + # - tout d'abord les grandeurs generiques (deformations, contraintes etc.) + # - puis les grandeurs particulieres, par exemple specifique a une loi de comp ou une mesure de def + # enfin un mot cle de fin ( fin_grandeurs_maple) + 1 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud fin_list_ddl_noeud + type_sortie_ddl_retenue= 0 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element SIG11 SIG22 SIG33 SIG12 SIG23 SIG13 EPS11 EPS22 EPS33 EPS12 EPS23 EPS13 fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements +# informations particuliere dans le cas ou il y a une animation +# type_xi indique si oui ou non les grandeurs a tracer sont aux noeuds (sinon c'est au elements) +# x1 et x2 indiquent les noms des ddls des grandeurs en x et y. accroi_x1 et accroi_x2 indiquent +# si oui ou non x1 et x2 represente l'accroissement entre 0 et t de la grandeur ou bien la grandeur elle meme. + debut_info_particulier grandeur_au_noeud? 1 x1= NU_DDL x2= NU_DDL accroi_x1= 0 accroi_x2= 1 fin_info_particulier +# un parametre de pilotage du style de sortie + parametre_style_de_sortie 1 +# un parametre indiquant si les tenseurs sont en absolue (rep 1) ou suivant un repere ad hoc +# (tangent pour les coques, suivant la fibre moyenne pour les element 1D ) + tenseur_en_absolue_ 0 + fin_grandeurs_maple # fin des grandeurs a visualiser au format maple + + # ----------------------------- definition des parametres d'animation: ---------------- + debut_animation # un mot cle de debut de liste (debut_animation) + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # des parametres avec des valeurs: (sur une meme ligne) + cycleInterval 8 # cycleInterval (indique le temps en seconde du cycle de l'animation) + fin_animation # un mot cle de fin + + fin_visualisation_maple + # ================================================================================= + # || fin de la visualisation maple || + # ================================================================================= + + + + fin_fichier_commande_visu # <<<<<< le mot cle permet + # l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus + # sans le mot cle de fin de fichier, le fichier n'est pas valide + + ############################################################################################### + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction.info new file mode 100644 index 0000000..778e2c3 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction.info @@ -0,0 +1,165 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Copyright (c) 1997-2015 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.697 + + +####################################################################### +# fichier de commande # + + +# --- commentaire sur l'objet du travail ??? + + + +#---------------------------------------------------- +# 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 + + +#----------------------------------------------------------------------- +#| parametres (falcultatifs ) associes au calcul implicite statique | +#----------------------------------------------------------------------- + + ## PARA_TYPE_DE_CALCUL + # ................................................................ + # / acceleration de convergence par extrapolation methode MMPE / + #.................................................................. + ## acceleration_convergence_= 1 cas_acceleration_convergence_= 1 nb_vec_cycle_= 8 + # + # acceleration_convergence_= indique si oui (=1) ou non (=0) on veut une acceleration + # cas_acceleration_convergence_= 1 : on projette sur les vecteurs positions + # 2 : on projette sur les vecteurs residus + # nb_vec_cycle_= indique le nombre de vecteur maxi que l'on prendre en compte + # pour la projection + #.................................................................. + # / mode debug / + #.................................................................. +#mode_debug_= 5 + # + # mode_debug_= parametre de debug (=0) calcul classique + # 1 : a chaque iteration on sort une visualisation + # **important**: ce parametre doit etre mis sur une ligne differente de + # l'acceleration de convergence + # + + + + +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- + + # -- def maillage + < cube1ele.her + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +E_to orgeas + + materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ + orgeas ISOHYPER3DORGEAS1 +# ....... loi de comportement 3D hyperelastique isotrope Orgeas 1 ........ +#----------------------------------------------------------------------------------------- +# K | Q0s | mu01 | mu02 | mu03 | alpha1 | alpha2 | Q0e | +#----------------------------------------------------------------------------------------- +# parametres coherents avec des membranes, ajustes sur resultats experimentaux +624 2 30 14.6 90 1e-3 1e-3 0.12 + + + +# -- definition du type de deformation (par defaut: DEFORMATION_STANDART) -- + type_de_deformation DEFORMATION_STANDART +#--------------------------------------- fin def des lois de comportement ------------------------ + +# --- divers stockages (1) ------- + masse_volumique #--------# +E_to 1 + + + charges +# ici pas de chargement a definir + + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- + +# conditions limites pour l essai de traction +N_droit UX +N_arriere UY +N_bas UZ +N_gauche 'UX=20' + + + +controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +SAUVEGARDE 0 +ITERATIONS 300000 +PRECISION 0.005 +MULTIPLICATEUR 1 +DELTAt 0.005 +TEMPSFIN 1 +RESTART 0 +DELTAtMAXI 0.1 + + + + para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- + + FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + +# ------------------------------------------------------------------------------------ + resultats pas_de_sortie_finale_ + COPIE 0 +# + + + _fin_point_info_ + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction_princ.maple.ref1 new file mode 100644 index 0000000..1e77eff --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Orgeas/Test_R_cube_isohyperd3dorgeas1_traction/cube_traction_princ.maple.ref1 @@ -0,0 +1,113 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.703 # +# 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 -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 12 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# element_1 pt_integ_1: [2]X [3]Y [4]Z [5] SIG11 [6] SIG22 [7] SIG33 [8] SIG12 [9] SIG23 [10] SIG13 [11] EPS11 [12] EPS22 [13] EPS33 [14] EPS12 [15] EPS23 [16] EPS13 +# +# +#==================================================================== +#|| fin du recapitulatif des differentes grandeurs || +#==================================================================== + + # ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant + # a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun + # correspondant a un maillage. On trouve pour chaque enregistrement successivement : + # s'il y a des grandeurs globales: le temps puis les grandeurs globales, + # puis s'il y a des torseurs de reaction : + # de nouveau le temps, les composantes de la resultante puis les composantes du moments + # donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels + # puis s'il y a des grandeurs aux noeuds: de nouveau le temps + # les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud + # puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds + # puis s'il y a des grandeur aux elements: + # le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques) + # suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration + # correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes + # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs + # etc. pour tous les points d'integration - element + + 5.000000000000e-03 2.115361902706e+01 7.883599917067e+01 7.883599917067e+01 1.246754248083e-01 9.260057771884e-05 9.260057771069e-05 2.698641420015e-34 3.179201731063e-15 -1.264538710493e-18 9.985019975030e-04 -3.998248006416e-04 -3.998248006417e-04 2.017792062434e-36 3.558397234536e-17 -1.419327565492e-20 +1.000000000000e-02 2.117475151360e+01 7.880448488185e+01 7.880448488185e+01 2.487203362453e-01 4.378397638141e-05 4.378397638141e-05 2.258352856448e-31 9.564229201673e-15 -2.525876536287e-18 1.994015960096e-03 -8.001293727149e-04 -8.001293727149e-04 2.531420612556e-33 1.068658098527e-16 -2.838090712908e-20 +1.500000000000e-02 2.119588400014e+01 7.877297059304e+01 7.877297059304e+01 3.721373352617e-01 -1.471408770181e-04 -1.471408770345e-04 -3.151671773921e-30 2.877267804918e-14 3.152852087291e-15 2.986553798226e-03 -1.200914484132e-03 -1.200914484132e-03 -1.261424708586e-32 3.209396208290e-16 3.546338609505e-17 +2.000000000000e-02 2.121701648668e+01 7.874145630423e+01 7.874145630423e+01 4.949289904455e-01 -4.808510016829e-04 -4.808510016829e-04 -4.052626957671e-30 3.847550531748e-14 3.147627423472e-15 3.976127363058e-03 -1.602180904342e-03 -1.602180904342e-03 -1.513412655257e-32 4.284334921657e-16 3.544224318952e-17 +2.707106781187e-02 2.124690233575e+01 7.869688836958e+01 7.869688836958e+01 6.675199999052e-01 -1.197819588030e-03 -1.197819588014e-03 -1.441168458111e-29 5.795042623729e-14 6.287355878473e-15 5.370558286617e-03 -2.170480461629e-03 -2.170480461629e-03 -7.071036449980e-32 6.437343313157e-16 7.090164628709e-17 +3.414213562373e-02 2.127678818483e+01 7.865232043493e+01 7.865232043493e+01 8.388724111112e-01 -2.203487543149e-03 -2.203487543132e-03 -2.922081750873e-30 8.081701077360e-14 1.562194203704e-15 6.759117391982e-03 -2.739746364703e-03 -2.739746364703e-03 -1.170820325417e-33 8.955862966472e-16 1.764296174247e-17 +4.121320343560e-02 2.130667403390e+01 7.860775250028e+01 7.860775250028e+01 1.008993068946e+00 -3.499578063698e-03 -3.499578063698e-03 -1.080627929555e-29 1.038700094798e-13 3.125685958814e-15 8.141837600510e-03 -3.309980805726e-03 -3.309980805726e-03 -4.032467466340e-32 1.148294254774e-15 3.535313377434e-17 +4.828427124746e-02 2.133655988297e+01 7.856318456563e+01 7.856318456563e+01 1.177888593891e+00 -5.087707111604e-03 -5.087707111604e-03 -2.928091900597e-29 1.205921932226e-13 6.248979357515e-15 9.518751603152e-03 -3.881185983081e-03 -3.881185983081e-03 -1.414695120625e-31 1.329977917946e-15 7.078411739725e-17 +5.828427124746e-02 2.137882485605e+01 7.850153129136e+01 7.850153129136e+01 1.420924356942e+00 1.778836690065e-07 1.778836690065e-07 3.122677735298e-15 1.520889211799e-17 -1.449068593345e-17 1.145615326926e-02 -4.672969880109e-03 -4.672969880109e-03 3.544598261520e-17 1.671791241742e-19 -1.644859461076e-19 +6.828427124746e-02 2.142108982913e+01 7.843987801710e+01 7.843987801710e+01 1.661553887406e+00 4.518916700280e-03 4.518916700280e-03 9.344644646674e-15 -1.252314128054e-13 -1.574330285611e-15 1.338209845641e-02 -5.466621522972e-03 -5.466621522972e-03 1.062950350269e-16 -1.372036130272e-15 -1.790795682236e-17 +7.828427124746e-02 2.146335480221e+01 7.837822474283e+01 7.837822474283e+01 1.899791408606e+00 8.466112807827e-03 8.466112807818e-03 1.242853128120e-14 -2.452641324545e-13 -1.927644176175e-17 1.529667731467e-02 -6.262146790748e-03 -6.262146790748e-03 1.416702459251e-16 -2.678327586566e-15 -2.197281547730e-19 +8.828427124746e-02 2.150561977529e+01 7.831450148208e+01 7.831450148208e+01 2.126190274600e+00 4.056138742312e-07 4.056138656686e-07 3.099636839317e-15 -6.644225028668e-15 1.528084042314e-15 1.719997910912e-02 -7.086356845695e-03 -7.086356845695e-03 3.540550292127e-17 -7.231573320331e-17 1.745449122873e-17 +1.024264068712e-01 2.156539147344e+01 7.822665044411e+01 7.822665044411e+01 2.452639937988e+00 4.893950551095e-07 4.893950551095e-07 4.538919521320e-31 -6.693345758016e-15 3.063706421145e-15 1.987257471975e-02 -8.225944875921e-03 -8.225944875921e-03 -6.704925374937e-49 -7.251909107094e-17 3.509917242802e-17 +1.165685424949e-01 2.162516317158e+01 7.813879940613e+01 7.813879940613e+01 2.774254243306e+00 -1.099583339563e-03 -1.099583339519e-03 -3.077833676010e-15 -3.355380061689e-15 3.049395342465e-15 2.252303990196e-02 -9.369378764813e-03 -9.369378764812e-03 -3.536830484083e-17 -3.619287119000e-17 3.504151146736e-17 +1.307106781187e-01 2.168493486973e+01 7.805094836815e+01 7.805094836815e+01 3.090634186812e+00 -3.070516967621e-03 -3.070516967568e-03 -9.198359692841e-15 -1.361332398737e-14 7.633554271603e-15 2.515161831774e-02 -1.051667583709e-02 -1.051667583709e-02 -1.060507807566e-16 -1.462360687732e-16 8.800964709843e-17 +1.448528137424e-01 2.174470656787e+01 7.796309733017e+01 7.796309733017e+01 3.397608088517e+00 -3.640918407576e-03 -3.640918407496e-03 -1.219232805421e-14 -1.368260537136e-14 9.109292047365e-15 2.775855028480e-02 -1.166785351515e-02 -1.166785351515e-02 -1.413303313424e-16 -1.466900077400e-16 1.055925708056e-16 +1.648528137424e-01 2.182923651403e+01 7.782695993902e+01 7.782695993902e+01 3.630891115443e+00 1.691978495235e-05 1.691978496919e-05 2.858117939270e-15 4.313592403843e-17 1.391349847677e-15 3.140881587810e-02 -1.345947045607e-02 -1.345947045607e-02 3.531900230236e-17 4.877875694900e-19 1.719351318512e-17 +1.848528137424e-01 2.191376646020e+01 7.768538164037e+01 7.768538164037e+01 3.774726117038e+00 3.036236455006e-06 3.036236455006e-06 1.952285839935e-31 -6.110243048725e-15 1.283961988856e-15 3.501692149978e-02 -1.533269185833e-02 -1.533269185833e-02 -1.908966389196e-35 -7.379037166445e-17 1.712628670278e-17 +2.048528137424e-01 2.199829640636e+01 7.754380334172e+01 7.754380334172e+01 3.911298424153e+00 -3.353151710913e-03 -3.353151710898e-03 -2.476766896937e-15 -1.165918454092e-14 2.435196911312e-15 3.858351391342e-02 -1.721618292187e-02 -1.721618292187e-02 -3.530399559257e-17 -1.488837228210e-16 3.471145432795e-17 +2.248528137424e-01 2.208282635252e+01 7.740222504307e+01 7.740222504307e+01 4.041777089805e+00 -1.066073173812e-02 -1.066073173811e-02 -7.009631604381e-15 -1.117294446732e-14 1.124861022343e-15 4.210922752774e-02 -1.911001885331e-02 -1.911001885331e-02 -1.058928929021e-16 -1.496347314162e-16 1.699301682192e-17 +2.531370849898e-01 2.220236974881e+01 7.720702231727e+01 7.720702231727e+01 4.242156394540e+00 2.577368087696e-06 2.577368087696e-06 -2.176771897216e-15 -2.617435738778e-15 -4.583015611096e-17 4.702676176271e-02 -2.173827173030e-02 -2.173827173030e-02 -3.528532883763e-17 -3.707173389460e-17 -7.429038068358e-19 +2.814213562373e-01 2.232191314510e+01 7.701181959148e+01 7.701181959148e+01 4.431240939974e+00 2.669684945765e-03 2.669684945752e-03 4.097350408477e-15 1.292539703571e-14 -3.121216608497e-15 5.186550124481e-02 -2.438653551191e-02 -2.438653551191e-02 7.054900914095e-17 1.916864592340e-16 -5.374173968333e-17 +3.097056274848e-01 2.244145654139e+01 7.681661686568e+01 7.681661686568e+01 4.609222836896e+00 -2.680564849878e-03 -2.680564849898e-03 7.774154339240e-15 1.753198745518e-14 -4.909386874531e-15 5.662712042139e-02 -2.705501385940e-02 -2.705501385940e-02 1.410605927022e-16 2.701453759769e-16 -8.907991687667e-17 +3.379898987322e-01 2.256099993768e+01 7.662141413988e+01 7.662141413988e+01 4.776266882282e+00 -1.605509988623e-02 -1.605509988626e-02 1.113349612981e-14 3.659356340441e-14 -3.764011254849e-15 6.131324949586e-02 -2.974391303156e-02 -2.974391303156e-02 2.115435002830e-16 5.822630148238e-16 -7.151860535730e-17 +3.779898987322e-01 2.273005983000e+01 7.635399097803e+01 7.635399097803e+01 5.039088489648e+00 7.786784323796e-06 7.786784310695e-06 1.753430675707e-15 -4.667414834065e-15 -2.686028178850e-15 6.781464266963e-02 -3.346117460201e-02 -3.346117460201e-02 3.524057961813e-17 -7.705718981576e-17 -5.398399332503e-17 +4.179898987322e-01 2.289911972233e+01 7.608656781617e+01 7.608656781617e+01 5.281101461860e+00 4.188695379502e-04 4.188695379874e-04 -8.349974158114e-15 -4.180499680814e-14 -2.563877898162e-15 7.417257236372e-02 -3.721770044001e-02 -3.721770044001e-02 -1.761336499809e-16 -7.108068172303e-16 -5.408222393956e-17 +4.579898987322e-01 2.306817961465e+01 7.581914465432e+01 7.581914465432e+01 5.502768903361e+00 -1.486716933725e-02 -1.486716933719e-02 -1.920581604346e-14 -7.341157130583e-14 -1.662346231887e-15 8.039122873540e-02 -4.101404548037e-02 -4.101404548037e-02 -4.225881034147e-16 -1.278036045639e-15 -3.657682338320e-17 +4.979898987322e-01 2.323723950697e+01 7.555978607164e+01 7.555978607164e+01 5.747205570340e+00 6.846022119505e-06 6.846022125279e-06 -4.625456062906e-15 9.551384603328e-17 -4.304665733468e-15 8.647465007341e-02 -4.473447565687e-02 -4.473447565687e-02 -1.055996277891e-16 1.694392636938e-18 -9.827595225819e-17 +5.545584412272e-01 2.347632629955e+01 7.519299764650e+01 7.519299764650e+01 6.072495284825e+00 -1.349997809418e-02 -1.349997809426e-02 1.477529737856e-14 1.046096116767e-13 -4.870811226558e-15 9.485458158242e-02 -5.006182316306e-02 -5.006182316306e-02 3.518213344976e-16 1.887643342859e-15 -1.159811042653e-16 +6.111269837221e-01 2.371541309213e+01 7.494407477561e+01 7.494407477561e+01 7.893443913008e+00 1.174474057125e-03 1.174474057110e-03 7.071438838960e-31 -2.672724509365e-15 5.209333833797e-15 1.029823461463e-01 -5.372189528702e-02 -5.372189528702e-02 4.982209748012e-33 -3.934430703275e-17 1.034334563857e-16 +6.676955262170e-01 2.395449988471e+01 7.471413002634e+01 7.471413002634e+01 9.921062208520e+00 3.190389927665e-05 3.190389928598e-05 2.062306067108e-15 1.790013034956e-16 4.069697219451e-16 1.108679607376e-01 -5.713547648234e-02 -5.713547648234e-02 3.492323853827e-17 2.200646573397e-18 6.891654397969e-18 +7.242640687119e-01 2.419358667729e+01 7.448418527706e+01 7.448418527706e+01 1.187711853859e+01 -3.449522271665e-03 -3.449522271601e-03 9.231584851297e-15 -7.648278494211e-15 -7.628931811219e-15 1.185209498365e-01 -6.058072128757e-02 -6.058072128757e-02 1.391677793116e-16 -8.200736945582e-17 -1.150075004227e-16 +8.042640687119e-01 2.453170646193e+01 7.415899429404e+01 7.415899429404e+01 1.452614570673e+01 -1.360483111091e-02 -1.360483111078e-02 1.811545479203e-14 -4.344610922162e-15 -2.279882166040e-14 1.289643090895e-01 -6.550785078148e-02 -6.550785078148e-02 2.422979411704e-16 -4.018174953881e-17 -3.049389382075e-16 +8.842640687119e-01 2.486982624658e+01 7.383380331102e+01 7.383380331102e+01 1.704568019323e+01 -3.271895550780e-02 -3.271895550757e-02 3.369765411616e-14 -1.540453391274e-14 -3.668502104018e-14 1.389846198943e-01 -7.050022616153e-02 -7.050022616153e-02 4.133377972658e-16 -1.271366858575e-16 -4.499810502276e-16 +9.642640687119e-01 2.520794603123e+01 7.350861232800e+01 7.350861232800e+01 1.944401396488e+01 -6.385033599079e-02 -6.385033599044e-02 4.475374837916e-14 -1.125837325714e-14 -4.985158812666e-14 1.486044270979e-01 -7.555900452951e-02 -7.555900452950e-02 5.142620258237e-16 -8.477323560279e-17 -5.728409268280e-16 +1.000000000000e+00 2.535898384862e+01 7.337486990659e+01 7.337486990659e+01 2.052550812407e+01 1.213215295226e-05 1.213215295226e-05 3.046647892073e-15 4.224030430256e-16 1.197441103707e-14 1.527777777778e-01 -7.765909162409e-02 -7.765909162409e-02 3.420428927494e-17 3.072528279066e-18 1.344350359864e-16