From db3f715d42ede7406cfc700199e4286294f68480 Mon Sep 17 00:00:00 2001 From: rio <> Date: Mon, 26 Sep 2016 16:47:02 +0000 Subject: [PATCH] ajout test R sur l utilisation de grandeurs locales de ponderation. Exemple egalement d utilisation de l hysteresis bulk --- ...NTRAINTES_PLANES_traction_princ.maple.ref1 | 4 +- .../README | 9 +- .../Joint17_1x1.her | 81 +++++++ .../Test_R_fct3D_variables_locales/README | 55 +++++ .../courbe_ponderation.txt | 34 +++ .../elastomere_joint.loi | 81 +++++++ .../test_fct3D.CVisu1 | 138 ++++++++++++ .../test_fct3D.info | 206 ++++++++++++++++++ .../test_fct3D_princ.maple.ref1 | 108 +++++++++ 9 files changed, 710 insertions(+), 6 deletions(-) create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/Joint17_1x1.her create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/README create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/courbe_ponderation.txt create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/elastomere_joint.loi create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.CVisu1 create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.info create mode 100644 Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D_princ.maple.ref1 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction_princ.maple.ref1 index 87566f7..be466ea 100644 --- a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction_princ.maple.ref1 +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction_princ.maple.ref1 @@ -1,6 +1,6 @@ #fichier au format maple6 ############################################################################################### -# Visualisation elements finis : Herezh++ V6.752 # +# Visualisation elements finis : Herezh++ V6.771 # # Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # # http://www-lg2m.univ-ubs.fr # ############################################################################################### @@ -98,4 +98,4 @@ # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs # etc. pour tous les points d'integration - element - 1.000000000000e+00 2.176646113673e-01 2.094640616973e-01 0.000000000000e+00 2.870204543312e-02 -8.923101412272e-03 2.947035243163e+02 2.955880224154e-02 -8.844414469719e-03 2.176646113673e-01 2.094640616973e-01 0.000000000000e+00 9.911156715684e-03 + 1.000000000000e+00 2.176646113673e-01 2.094652158455e-01 0.000000000000e+00 2.870204543312e-02 -8.917493119952e-03 2.947217216813e+02 2.955880224154e-02 -8.838904478911e-03 2.176646113673e-01 2.094652158455e-01 0.000000000000e+00 9.911133403096e-03 diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/README b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/README index 06b63b1..8fc601a 100644 --- a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/README +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/HYPO_ELAS/Test_R_HYPO_ELAS3D_LOI_CONTRAINTES_PLANES_traction/README @@ -17,10 +17,11 @@ d ------------------------------------------------------ But du test ------------------------------------------------------ -****** -****** ATTENTION : pour l instant, le fichier .maple.ref1 est celui généré par le -****** cas test Test_R_HYPO_ELAS2D_C_traction en attendant la correction du problème -****** +%****** +%****** ATTENTION : pour l instant, le fichier .maple.ref1 est celui généré par le +%****** cas test Test_R_HYPO_ELAS2D_C_traction en attendant la correction du problème +%****** +% gérard rio -> mise à jour et correction le 12 septembre 2016 test simple de la loi HYPO_ELAS3D en contraintes planes (LOI_CONTRAINTES_PLANES) en traction uniaxiale : "" LOI_CONTRAINTES_PLANES ... diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/Joint17_1x1.her b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/Joint17_1x1.her new file mode 100644 index 0000000..19e16be --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/Joint17_1x1.her @@ -0,0 +1,81 @@ + +############################################################################# +# # +# 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 : 7.96 x 7.81 ( 1 x 1 elements ) + # geometrie rectangulaire, decoupage rectangulaire axi-symetrique, interpolation quadratique complete. + + nom_maillage Joint + + noeuds ------------ + 9 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 3.905 0 + 3 0 7.81 0 + 4 3.98 0 0 + 5 3.98 3.905 0 + 6 3.98 7.81 0 + 7 7.96 0 0 + 8 7.96 3.905 0 + 9 7.96 7.81 0 + + # references des noeuds + +N_O 1 2 3 +N_S 1 4 7 +N_E 7 8 9 +N_N 3 6 9 +N_SO 1 +N_SE 7 +N_NE 9 +N_NO 3 + N_to 1 2 3 4 5 6 7 8 9 + + + # les elements + + elements ---------- + 1 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 QUAD_AXI QUADRACOMPL 1 7 9 3 4 8 6 2 5 + + + # 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/Test_R_fct3D_variables_locales/README b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/README new file mode 100644 index 0000000..5afea61 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/README @@ -0,0 +1,55 @@ +------------------------------------------------------ +Auteur +------------------------------------------------------ +GŽrard Rio (gerard.rio@univ-ubs.fr) + + +------------------------------------------------------ +Mots-cles +------------------------------------------------------ +contrainte_individuelle_a_chaque_loi_a_t +contrainte_individuelle_a_t_sans_proportion +LOI_ADDITIVE_EN_SIGMA +avec_fonction_de_ponderation_ +les_grandeur_ponderation= +def_duale_mises_maxi +def_duale_mises +ISOHYPERBULK_GENE +avec_ponder_grandeur_locale_ +PRESSION_HYST_T +PRESSION_HYST_REF +PRESSION_HYST_REF_M1 +HYSTERESIS_BULK +les_Fonctions_nD +FONCTION_EXPRESSION_LITTERALE_nD + + +------------------------------------------------------ +But du test +------------------------------------------------------ +Test d'une loi additive pondŽrŽe par des grandeurs calculŽes localement par une loi membre +ici il s'agit de la pression d'hystŽrŽsis ˆ t, de la dernire pression d'hystŽrŽsis de rŽfŽrence, +et enfin de la pression prŽcŽdente de rŽfŽrence cf. le fichier elastomere_joint.loi + +Dans ce test on se sert d'une fonction littŽrale nD ˆ 3 variables + + +------------------------------------------------------ +Description du calcul +------------------------------------------------------ +Le test est issu d'un exemple de calcul fourni par ƒmilie Vieville dans le cadre de sa thse de doctorat. + +On considre un tube creux de section rectangulaire soumis ˆ sur la face supŽrieure ˆ un dŽplacement de compression. +Le tube est modŽlisŽ par un seul ŽlŽment quadrangulaire axisymŽtrique. Les 3 autres faces sont bloquŽes perpendiculairement +ˆ la matire. +La loi de comportement est une combinaison additive d'hyperŽlasticitŽ et d'hystŽrŽsis sphŽrique et dŽviatorique. +La description de la loi est complexe. + +------------------------------------------------------------- +Grandeurs de comparaison +------------------------------------------------------------- +Principalement: +- rŽaction due au dŽplacement +- au premier point d'intŽgration: contrainte, dŽformation, globales et locales, + pondŽration et pressions + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/courbe_ponderation.txt b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/courbe_ponderation.txt new file mode 100644 index 0000000..8c86a3f --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/courbe_ponderation.txt @@ -0,0 +1,34 @@ + courbe_bulk COURBEPOLYLINEAIRE_1_D + Debut_des_coordonnees_des_points + #Coordonnee dim= 2 0.0 0.15 #parametres base + #Coordonnee dim= 2 0.08 0.15 #parametres base + #Coordonnee dim= 2 2.0 0.15 #parametres base + Coordonnee dim= 2 0.0 0.17 + Coordonnee dim= 2 2.0 0.17 + Fin_des_coordonnees_des_points + + + courbe_orgeas COURBEPOLYLINEAIRE_1_D + Debut_des_coordonnees_des_points + #Coordonnee dim= 2 0.0 0.50 #parametres base + #Coordonnee dim= 2 0.08 0.50 #parametres base + #Coordonnee dim= 2 2.0 0.50 #parametres base + Coordonnee dim= 2 0.0 3.50 + Coordonnee dim= 2 2.0 0.50 + Fin_des_coordonnees_des_points + + courbe_hysteresis_3D COURBEPOLYHERMITE_1_D + Debut_des_coordonnees_des_points + #Coordonnee dim= 3 0.0 0.05 0.6 #parametres base + #Coordonnee dim= 3 2.0 7.0 3. #parametres base + Coordonnee dim= 3 0.0 0.05 0.6 + Coordonnee dim= 3 2.0 5.5 3.5 + Fin_des_coordonnees_des_points + + courbe_hysteresis_bulk COURBEPOLYHERMITE_1_D + Debut_des_coordonnees_des_points + #Coordonnee dim= 3 0.0 0.05 0.6 #parametres base + #Coordonnee dim= 3 2.0 7.0 3. #parametres base + Coordonnee dim= 3 0.0 0.05 0.6 + Coordonnee dim= 3 2.0 5.5 3.5 + Fin_des_coordonnees_des_points \ No newline at end of file diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/elastomere_joint.loi b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/elastomere_joint.loi new file mode 100644 index 0000000..6e4bb00 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/elastomere_joint.loi @@ -0,0 +1,81 @@ +elastomere LOI_ADDITIVE_EN_SIGMA + +avec_fonction_de_ponderation_ + + #--------------------------------------- + # loi 1: hyperbulk gene + #--------------------------------------- + les_grandeur_ponderation= def_duale_mises_maxi AuPti_ fin_grandeur_ponderation_ + deb_fonct_ponder= courbe_bulk fin_fonct_ponder_ + # ....... loi de comportement 3D hyperelastique isotrope ISOHYPERBULK_GENE + ISOHYPERBULK_GENE + #-------------------------" + # K | fonction |" + #-------------------------" + omega_V= courbe_polynomiale sortie_post_ 1 + + #--------------------------------------- + # loi 2: Orgeas sans partie sphérique + #--------------------------------------- + les_grandeur_ponderation= def_duale_mises_maxi AuPti_ fin_grandeur_ponderation_ + deb_fonct_ponder= courbe_orgeas fin_fonct_ponder_ + ISOHYPER3DORGEAS1 + 0 10.0 0.1 1. 11.5 0.1 0.1 0.04 sortie_post_ 1 + + #--------------------------------------- + # loi 3: Hystérésis déviatorique + #--------------------------------------- + les_grandeur_ponderation= def_duale_mises AuPti_ fin_grandeur_ponderation_ + deb_fonct_ponder= courbe_hysteresis_3D fin_fonct_ponder_ + avec_ponder_grandeur_locale_ + les_grandeurs_= PRESSION_HYST_T PRESSION_HYST_REF PRESSION_HYST_REF_M1 fin_grandeurs_ + deb_fonct_= fct_dep_pression fin_fonct_ + HYSTERESIS_3D + np= 1.75 mu= 420 Qzero= 26 avec_parametres_de_reglage_ + type_de_resolution_ 5 type_calcul_comportement_tangent_ 2 + cas_kutta_ 5 erreurAbsolue_ 1.e-4 erreurRelative_ 1.e-3 + nbMaxiAppel_ 3000 + + nb_iteration_maxi_ 30 + nb_dichotomie_maxi_ 20 + tolerance_residu_ 1.e-4 + tolerance_residu_rel_ 1.e-3 + + tolerance_coincidence_ 1.e-3 + nb_maxInvCoinSurUnPas_ -20 + mini_Qsig_pour_phase_sigma_Oi_tdt_ 1 + mini_rayon_ 1.e-12 min_angle_trajet_neutre_ 1.e-2 + possibilite_cosAlphaNegatif_ 0 + permet_affichage_ 5 + type_de_transport_memorisation_ -1 + fin_parametres_reglage_Hysteresis_ + #-------- fin hysteresis_3D déviatorique -------- + + #--------------------------------------- + # loi 4: Hystérésis bulk + #--------------------------------------- + les_grandeur_ponderation= def_duale_mises AuPti_ fin_grandeur_ponderation_ + deb_fonct_ponder= courbe_hysteresis_bulk fin_fonct_ponder_ + avec_ponder_grandeur_locale_ + les_grandeurs_= PRESSION_HYST_T PRESSION_HYST_REF PRESSION_HYST_REF_M1 fin_grandeurs_ + deb_fonct_= fct_dep_pression fin_fonct_ + HYSTERESIS_BULK #hystérésis sphérique + np= 1.75 mu= 1000 Qzero= 26 avec_parametres_de_reglage_ + type_de_resolution_ 2 + cas_kutta_ 5 + erreurAbsolue_ 1.e-5 erreurRelative_ 1.e-4 + #erreurAbsolue_ 1.e-8 erreurRelative_ 1.e-8 + nbMaxiAppel_ 3600 + nb_iteration_maxi_ 30 nb_dichotomie_maxi_ 20 + tolerance_residu_ 1.e-3 tolerance_residu_rel_ 1.e-5 + depassement_Q0_ 1.1 + sortie_post_ 1 + permet_affichage_ 0 + depassement_Q0_ 1.1 + tolerance_coincidence_ 1.e-10 + nb_maxInvCoinSurUnPas_ -20 + sortie_post_ 1 + fin_parametres_reglage_Hysteresis_ + #-------- fin hysteresis_bulk -------- + +fin_liste_lois_elementaires diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.CVisu1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.CVisu1 new file mode 100644 index 0000000..33948d7 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.CVisu1 @@ -0,0 +1,138 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V6.756 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + debut_fichier_commande_visu # >>>>>> le mot cle: + # permet au programme de se positionner au debut du fichier, il est indispensable + + + # ================================================================================= + # || ***** demande d'une visualisation maple: ***** || + # ================================================================================= + # un mot cle de debut (debut_visualisation_maple) + # un mot cle de fin ( fin_visualisation_maple) + # la seule presence de ces deux mots cle suffit a activer la visualisation maple + debut_visualisation_maple + + # ----------------------------- definition de la liste des increments a balayer: ---------------- + debut_list_increment # un mot cle de debut de liste + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments) + # un mot cle de fin de liste ( fin_list_increment) + fin_list_increment + + # ----------------------------- choix des maillages a visualiser: ---------------- + # la liste est facultative, par defaut la visualisation concerne le premier maillage + debut_choix_maillage # un mot cle de debut, + actif 0 # <0 ou 1> indique si l'ordre est actif ou non + # une liste d'entiers , puis , sur une meme ligne + 1 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales 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 N_N fin_list_torseur_reaction + # les sommes, moyennes etc. sur ref de noeuds: un mot cle de debut: debut_list_SM_sur_refN + # une liste de nom de ref , puis + debut_list_SM_sur_refN fin_list__SM_sur_refN + # debut de la liste des ddl a considerer , (une liste de ddl), puis + deb_list_ddl_SM_sur_refN fin_list_ddl_SM_sur_refN + # debut de la liste des ddl etendus a considerer, (une liste de ddl), puis + deb_list_ddl_etendu__SM_sur_refN fin_list_ddl_etendu__SM_sur_refN # fin des ddl etendu pour _SM_sur_refNs + # debut de la liste des grandeurs particulieres a considerer, (une liste de string), puis + deb_list_GrandParticuliere__SM_sur_refN fin_list_GrandParticuliere__SM_sur_refN # fin des grandeurs particulieres aux _SM_sur_refNs + # les SM sur ref d'elements: un mot cle de debut: debut_list_SM_sur_refE + # une liste de nom de ref d'elements , le pti associe + # puis + debut_list_SM_sur_refE fin_list__SM_sur_refE + # les SMs sur ref de pti: un mot cle de debut: debut_list_SM_sur_refpti + # une liste de nom de ref , puis + debut_list_SM_sur_refpti fin_list__SM_sur_refpti + # debut de la liste des ddl a considerer , (une liste de ddl), puis + deb_list_ddl__SM_sur_refEpti fin_list_ddl__SM_sur_refEpti # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les _SM_sur_refEptis, (une liste de string), puis + deb_list_GrandParticuliere__SM_sur_refEpti fin_list_GrandParticuliere__SM_sur_refEpti # fin des grandeurs particulieres aux _SM_sur_refEptis + # tableau de grandeurs evoluees aux _SM_sur_refEptis a visualiser, un par maillage + deb_list_GrandEvoluee__SM_sur_refEpti fin_list_GrandEvoluee__SM_sur_refEpti + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud 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= 1 + # debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds + # debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis + deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds + fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds + debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements + # debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis + deb_list_NbElement_NbPtInteg fin_list_NbElement_NbPtInteg + # debut de la liste des ref d'elements, puis une liste de: nom + numero d'integ, puis + deb_list_ref_element E_to 5 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 EPS11 EPS22 EPS33 def_equivalente Q_eps 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 contrainte_individuelle_a_chaque_loi_a_t FONC_PONDERATION V_vol PRESSION_HYST_REF PRESSION_HYST POTENTIEL 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/Test_R_fct3D_variables_locales/test_fct3D.info b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.info new file mode 100644 index 0000000..a918d94 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D.info @@ -0,0 +1,206 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# 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 avec le plus de verifications pendant le calcul et les I/O ) + version 6.686 + + +####################################################################### +# fichier de commande # + + +# --- écrasement du joint, avec contact entre le joint et la chemise fine + + + +#---------------------------------------------------- +# definition de la dimension de l'espace de travail | +#---------------------------------------------------- + dimension 3 + + +#--------------------------------------------------------------- +# definition facultative du niveau d'impression (entre 0 et 10)| +#--------------------------------------------------------------- + niveau_commentaire 0 + +#------------------------------- +# definition du type de calcul | +#------------------------------- + TYPE_DE_CALCUL +non_dynamique #avec plus visualisation + +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- + # -- def maillage joint + < Joint17_1x1.her +def_mouvement_solide_initiaux_ +mouvement_solide_ +translation_= 15.015 7.6 0. +centre_= 0. 0. 0. +rotation_= 0. 0. 0. +fin_mouvement_solide_ + + +#----------------------------- +# Definition des courbes +#----------------------------- + les_courbes_1D +< courbe_ponderation.txt + +courbe_polynomiale COURBEPOLYNOMIALE # nom de la courbe puis le type de la courbe +# def des coefficients d’un polynome du quatrieme degre a0+a1x+a2x^2+a3x^3+a4x^4 +# debut_coef= a0 a1 a2 a3 a4 fin_coef +debut_coef= -171.403600607624 1489.547106 -3391.128916 2999.230915 -926.245505244461 fin_coef #w1_4(x) + +# def du chargement de traction + simple_char COURBEPOLYLINEAIRE_1_D + Debut_des_coordonnees_des_points + Coordonnee dim= 2 0.0 0.00000000 + Coordonnee dim= 2 1 0.39 + Coordonnee dim= 2 2 0.27 + Coordonnee dim= 2 3 0.39 + Coordonnee dim= 2 4 0.79 + Coordonnee dim= 2 5 0.55 + Coordonnee dim= 2 6 0.79 + Coordonnee dim= 2 7 1.17 + Coordonnee dim= 2 8 0.83 + Coordonnee dim= 2 9 1.17 + Coordonnee dim= 2 10 1.57 + Coordonnee dim= 2 11 1.10 + Coordonnee dim= 2 12 1.57 + Coordonnee dim= 2 13 1.98 + Coordonnee dim= 2 14 1.38 + Coordonnee dim= 2 15 1.98 + Fin_des_coordonnees_des_points + + +les_Fonctions_nD #------------ +# fonction Fonction_expression_litterale_nD +# f(P,PR,PRM1) +# P : la pression, +# PR : la dernière pression de référence +# PRM1: la pression de référence avant PR +fct_dep_pression FONCTION_EXPRESSION_LITTERALE_nD +un_argument= P un_argument= PR un_argument= PRM1 + +fct= (P > 0.) ? \ # d'abord le cas positif + ( ( (PR != 0.) && (P > PR) && (P <= (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ + ( 0.4/(PR- (PR+PRM1)/2.) * P + (1. - PR*0.4/(PR- (PR+PRM1)/2.))) \ # première moitiè de la remontée + : ( ( (PR != 0.) && (P > (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ + ( 0.4/(PRM1 - (PR+PRM1)/2.) * P + (1. - 0.4 * PRM1/(PRM1 - (PR+PRM1)/2.))) \ #deuxième moitié de remontée + : 1. \ # tous les autres cas + ) \ + ) \ + : ( ( (PR != 0.) && (P > PR) && (P <= (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ # sinon cas négatif + ( -0.4/(PR- (PR+PRM1)/2.) * P + (1. + PR*0.4/(PR- (PR+PRM1)/2.))) \ # première moitiè de la remontée + : ( ( (PR != 0.) && (P > (PR+PRM1)/2.) && ((PRM1 - PR) > 1.e-6) ) ? \ + ( -0.4/(PRM1 - (PR+PRM1)/2.) * P + (1. + 0.4 * PRM1/(PRM1 - (PR+PRM1)/2.))) \ #deuxième moitié de remontée + : 1. \ # tous les autres cas + ) \ + ) + + fin_parametres_fonction_expression_litterale_ + + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +nom_mail= Joint E_tout elastomere + + materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ + acier ISOELAS +# ....... loi de comportement isoelastique 3D ........ +# module d'young : coefficient de poisson + 210000 0.3 + +< elastomere_joint.loi + +#--------------------------------------- fin def des lois de comportement ------------------------ + +# --- divers stockages (1) ------- + masse_volumique #--------# +nom_mail= Joint E_tout 1 + + charges #------------# + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +# le joint +nom_mail= Joint N_S UY +nom_mail= Joint N_O UX +nom_mail= Joint N_N 'UY = COURBE_CHARGE: simple_char ECHELLE: -1', UZ +nom_mail= Joint N_tout UZ + +nom_mail= Joint N_E UX + + controle #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +SAUVEGARDE 0 +DELTAtMAXI 0.05 +DELTAtMINI 1.e-8 +TEMPSFIN 10 +DELTAt 0.01 +MAXINCRE 1000000 +ITERATIONS 16 +PRECISION 5.e-3 + + +para_pilotage_equi_global ------------ +#----------------------------------------------- +# PARAMETRE | VALEUR | +#----------------------------------------------- +FACTEUR_DIMINUTION 1.732 #1.732 +FACTEUR_AUGMENTATION 1.414 #1.414 +NB_BONNE_CONVERGENCE 3 +NORME_MAXI_INCREMENT 0.5#0.1 + +para_syteme_lineaire ------------ +#----------------------------------------------- +# PARAMETRE | VALEUR | +#----------------------------------------------- +TYPE_MATRICE BANDE_NON_SYMETRIQUE_LAPACK +SYMETRIE_MATRICE 0 + + para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +FREQUENCE_AFFICHAGE_ITERATION 20 +FREQUENCE_AFFICHAGE_INCREMENT 20 +FREQUENCE_SORTIE_FIL_DU_CALCUL DERNIER_CALCUL + +# ------------------------------------------------------------------------------------ + resultats pas_de_sortie_finale_ + COPIE 0 +# + + + _fin_point_info_ + + diff --git a/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D_princ.maple.ref1 new file mode 100644 index 0000000..fed5893 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/LOI_DE_COMPORTEMENT/Test_R_fct3D_variables_locales/test_fct3D_princ.maple.ref1 @@ -0,0 +1,108 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V6.774 # +# Copyright (c) 1997-2016, Gerard Rio (gerard.rio@univ-ubs.fr) http://kleger.univ-ubs.fr/Herezh/ # +# http://www-lg2m.univ-ubs.fr # +############################################################################################### + + + + # entete des donnees : informations gererales: on trouve successivement: + # >> le nombre de grandeurs globales (peut etre nul) suivi des identificateurs + # precedes du numero de colonne entre crochet + # >> le nombre de maillages m, et dimension de l'espace de travail + # puis pour chaque maillage, + # >> le nombre de torseurs de reaction (peut etre nul), le nombre total de reel qui va etre ecrit + # correspondant aux composantes des torseurs, puis les noms de ref associee suivi des positions + # des composantes entre crochet accolees a un identificateur: R pour reaction, M pour moment + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de noeuds + # >> le nombre de ref de noeuds, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de noeuds + # le nom de la reference, le nombre de noeud de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references d'elements ou de pti + # >> le nombre de ref d'element+ref de pti, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'element et de pti + # le nom de la reference, le nombre d'element de la ref, le nombre de grandeurs qui vont etre ecrites + # puis entre crochet la position suivi de la signification de chaque grandeur + # + # puis pour chaque maillage + # >> le nombre de noeud n (peut etre nul) ou il y a des grandeurs en sortie , + # puis le nombre des grandeurs p1 correspondantes, la position entre crochet des coordonnees + # et enfin l'idendificateur de ces grandeurs(p1 chaines de caractere) + # precedes du numero de colonne correspondant entre crochet + # puis pour chaque maillage + # >> le nombre de couples element-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # les grandeurs aux elements sont decomposees en 2 listes: la premiere de quantite P2 correspondant + # a des grandeurs generiques, la seconde de quantite P3 corresponds aux grandeurs specifiques, + # on trouve donc a la suite du nombre d'element: le nombre P2, suivi de P2 identificateurs de ddl + # chacun precedes du numero de colonne entre crochet + # puis le nombre P3, suivi de P3 identificateurs+categorie+type (chaines de caracteres), + # suivi entre crochet, de la plage des numeros de colonnes, correspondant + # chacun sur une ligne differentes + # ==== NB: pour les grandeurs specifique tensorielle: exemple d'ordre en 2D: + # tenseur symetrique, A(1,1) A(2,1) A(2,2), non symetrique A(1,1) A(1,2) A(2,1) A(2,2) + # en 3D c'est: tenseur symetrique, A(1,1) A(2,1) A(2,2) A(3,1) A(3,2) A(3,3) + # non symetrique A(1,1) A(1,2) A(2,1) A(2,2) A(2,3) A(3,1) A(3,2) A(3,3) + # ** dans le cas ou il n'y a qu'un seul increment en sortie, pour les grandeurs aux noeuds et aux elements, + # ** les informations peuvent etre decoupees selon: une ligne = un noeud, et le temps n'est pas indique + # ** ( cf: parametre_style_de_sortie = 0) + +#==================================================================== +#|| recapitulatif des differentes grandeurs par colonne || +#==================================================================== +#---------------------------------- grandeur globales ------------------------------------ +#0 (nombre de grandeurs globales) +#---------------------------------- maillage et dimension -------------------------------- +#1 3 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +#1 6 (nombre de torseurs et nombre total de grandeurs associees) +# N_N [2]Rx [3]Ry [4]Rz [5]Mx [6]My [7]Mz ; +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +#0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +#0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 8 60 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# reference E_to : element_1 pt_integ_5: [9]X [10]Y [11]Z [12] SIG11 [13] SIG22 [14] SIG33 [15] EPS11 [16] EPS22 [17] EPS33 [18] def_equivalente [19] Q_eps +# [20]X [21]Y [22]Z +# contrainte_individuelle_a_chaque_loi_a_t TABLEAU_T TENSEURHH 24 [23...46] +# FONC_PONDERATION TABLEAU_T SCALAIRE 4 [47...50] +# V_vol TABLEAU_T SCALAIRE 2 [51...52] +# PRESSION_HYST_REF TABLEAU_T SCALAIRE 1 [53...53] +# PRESSION_HYST TABLEAU_T SCALAIRE 1 [54...54] +# POTENTIEL TABLEAU_T SCALAIRE 2 [55...56] +# +# +# +#==================================================================== +#|| fin du recapitulatif des differentes grandeurs || +#==================================================================== + + # ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant + # a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun + # correspondant a un maillage. On trouve pour chaque enregistrement successivement : + # s'il y a des grandeurs globales: le temps puis les grandeurs globales, + # puis s'il y a des torseurs de reaction : + # de nouveau le temps, les composantes de la resultante puis les composantes du moments + # donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels + # puis s'il y a des grandeurs aux noeuds: de nouveau le temps + # les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud + # puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds + # puis s'il y a des grandeur aux elements: + # le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques) + # suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration + # correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes + # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs + # etc. pour tous les points d'integration - element + + 1.000000000000e-02 0.000000000000e+00 7.938767185078e+01 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 1.530036607400e+03 1.000000000000e-02 1.899500000000e+01 1.150305002486e+01 0.000000000000e+00 -4.596046510682e-02 -8.356288482106e-02 -4.596046510682e-02 0.000000000000e+00 -4.997340846383e-04 0.000000000000e+00 3.331560564256e-04 4.080311714803e-04 1.899500000000e+01 1.150305002486e+01 0.000000000000e+00 -4.952119141932e-02 0.000000000000e+00 -4.952119141932e-02 0.000000000000e+00 0.000000000000e+00 -4.952119141932e-02 1.574583935193e-03 0.000000000000e+00 -3.149955917467e-03 0.000000000000e+00 0.000000000000e+00 1.574583935193e-03 1.398994352691e-01 0.000000000000e+00 -2.797988705382e-01 0.000000000000e+00 0.000000000000e+00 1.398994352691e-01 -9.975098995414e-01 0.000000000000e+00 -9.975098995414e-01 0.000000000000e+00 0.000000000000e+00 -9.975098995414e-01 1.700000000000e-01 3.499500265915e+00 5.020008647164e-02 5.020008647164e-02 9.995006402049e-01 9.995006402049e-01 -0.000000000000e+00 9.975098995414e-01 1.147801657453e-05 7.859582947296e-07 +1.000000000000e+01 0.000000000000e+00 3.462151693250e+04 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 6.672600279171e+05 1.000000000000e+01 1.899500771143e+01 1.072033034002e+01 0.000000000000e+00 -5.363786760356e+00 -3.641099788708e+01 -5.363817712934e+00 1.894104332568e-11 -2.832544897515e-01 4.059711402842e-07 2.826528097308e-01 2.312764892791e-01 1.899500771143e+01 1.072033034002e+01 0.000000000000e+00 -5.387779170475e+01 -2.168404344971e-19 -5.387779170475e+01 0.000000000000e+00 0.000000000000e+00 -5.387779170475e+01 2.246465403689e+00 -4.229206433188e-04 -5.197979034130e+00 0.000000000000e+00 0.000000000000e+00 2.246471695490e+00 1.061403742687e+01 5.235653417365e-03 -2.122784528128e+01 0.000000000000e+00 0.000000000000e+00 1.061380785441e+01 -2.599999045748e+01 1.084202172486e-19 -2.599999045748e+01 0.000000000000e+00 0.000000000000e+00 -2.599999045748e+01 1.700000000000e-01 3.216745305882e+00 2.229871107671e-01 2.229871107671e-01 7.989759965917e-01 7.989759965917e-01 -0.000000000000e+00 2.599999045748e+01 6.204389175885e+00 4.452556487267e-01