diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/README b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/README new file mode 100644 index 0000000..e942a8d --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/README @@ -0,0 +1,72 @@ +------------------------------------------------------ +Auteur +------------------------------------------------------ +Gérard Rio (gerardrio56@free.fr) + +------------------------------------------------------ +Mots-cles +------------------------------------------------------ +contact solide-déformable +ISOELAS1D +espace 2D +comparaison analytique + +------------------------------------------------------ +But du test +------------------------------------------------------ +test d'un contact basique en 2D entre une barre représentée par 100 biellettes déformables +et un quadrangle indéformable de 50x100mm représenté par 20 triangles dans la hauteur. +Par rapport au test Test_R_contact_biel_quadrangle l'idée est ici d'étudier l'impact du nombre d'éléments sur la convergence et la précision du résultat. +: +"" + TYPE_DE_CALCUL + non_dynamique + + para_contact + CONTACT_TYPE 2 +"" + +------------------------------------------------------ +Description du calcul +------------------------------------------------------ +Les biellettes sont horizontales, constituées de 10 éléments. Seule les biellettes sont déformables. Sa déformation transversale est nulle. +Le maître de forme quadrangulaire est entièrement bloqué, et maillé en triangles linéaires: 2 triangles suivant x et 10 triangles suivant y. + +Cas académique très simple: +- on déplace le noeud sommet de la biel de 0.3mm maxi, sa longueur est de 100 -> 0.3% de def nominale + On est limité sur le déplacement du noeud à chaque incrément, sous peine de retourner la première biellette si le déplacement est supérieur à sa longueur (ici 1mm). +- le module d'Young: E=200000 MPa -> contrainte de 600 MPa +- la section de la biel est 50 mm^2 -> une réaction de 30000 N +-> et on trouve avec Herezh réaction = -2.890709228294e+04N +NB: PENALISATION_PENETRATION 0.9, TYPE_PENALISATION_PENETRATION 2, DISTANCE_MAXI_AU_PT_PROJETE 0.1 + +plus précisémment: def d'Almansi = 0.5*(1.-10000./99.7/99.7) = -3.0136e-03 +sig = -602.71 MPa -> Réaction = -3.0136e+04 N +d'où une différence avec le calcul d'Herezh : 4.0779% ce qui est moins précis qu'avec un seul élément et un quadrangle. +Le résultat est néanmoins cohérent. + +Avec un facteur de pénalisation de 10 on obtient une différence de 0.3792 pour une réaction de -3.002172277419e+04 + +On peut observer que l'épaisseur du quadrangle n'influe pas sur le calcul du facteur de pénalisation et donc sur le résultat. + +Avec 2 triangles au lieu de 20 on obtient une différence de 0.1187 pour une réaction de -3.010022585416e+04 +Cela provient de l'utilisation du longueur caractéristique (= sqrt (surface)) qui diminue à mesure que les triangles s'applatissent. + +Les calcul s'effectuent en 5 itérations. + +Le test utilise les paramètres: +CONTACT_TYPE 2 +PENALISATION_PENETRATION 10. + +------------------------------------------------------------- +Grandeurs de comparaison +------------------------------------------------------------- +biellette : + . [2]compteur_iteration_algo_global + pour la ref de noeuds N_deb (== le premier noeud auquel on impose un déplacement de 0.3mm) : + - N_deb [4]Rx [5]Ry [6]Mz + pour la ref de noeuds N_fi (== le dernier noeud) ce qui permet de connaître la pénétration + reference N_fi : noeud_101 [11]X [12]Y [13]delta_0_t_X1 + pour le point d integration 1 de la biellette : + - element_1 pt_integ_1: [15]X [16]Y [17] SIG11 [18] EPS11 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel.her new file mode 100644 index 0000000..409ec4c --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel.her @@ -0,0 +1,61 @@ + + + nom_maillage biel +############################################################################# +# # +# 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 d'element lineaire 1D | + # | pour une utilisation par herezh++ | + # ------------------------------------------------------------------ + + # barre de dimension : 100 ( 1 elements ) + + noeuds ------------ + 2 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 100 0 0 + + # references des noeuds + + N_deb 1 + N_fi 2 + N_to 1 2 + + + # les elements + + elements ---------- + 1 ELEMENTS + #-------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #-------------------------------------------------------------------- + 1 POUT BIE1 1 2 + + + + # references des arretes + + A_to 1 1 + + + # references des elements + + E_to 1 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel100.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel100.her new file mode 100644 index 0000000..68fc527 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel100.her @@ -0,0 +1,277 @@ + + + nom_maillage biel +############################################################################# +# # +# 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 d'element lineaire 1D | + # | pour une utilisation par herezh++ | + # ------------------------------------------------------------------ + + # barre de dimension : 100 ( 100 elements ) + + noeuds ------------ + 101 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 1 0 0 + 3 2 0 0 + 4 3 0 0 + 5 4 0 0 + 6 5 0 0 + 7 6 0 0 + 8 7 0 0 + 9 8 0 0 + 10 9 0 0 + 11 10 0 0 + 12 11 0 0 + 13 12 0 0 + 14 13 0 0 + 15 14 0 0 + 16 15 0 0 + 17 16 0 0 + 18 17 0 0 + 19 18 0 0 + 20 19 0 0 + 21 20 0 0 + 22 21 0 0 + 23 22 0 0 + 24 23 0 0 + 25 24 0 0 + 26 25 0 0 + 27 26 0 0 + 28 27 0 0 + 29 28 0 0 + 30 29 0 0 + 31 30 0 0 + 32 31 0 0 + 33 32 0 0 + 34 33 0 0 + 35 34 0 0 + 36 35 0 0 + 37 36 0 0 + 38 37 0 0 + 39 38 0 0 + 40 39 0 0 + 41 40 0 0 + 42 41 0 0 + 43 42 0 0 + 44 43 0 0 + 45 44 0 0 + 46 45 0 0 + 47 46 0 0 + 48 47 0 0 + 49 48 0 0 + 50 49 0 0 + 51 50 0 0 + 52 51 0 0 + 53 52 0 0 + 54 53 0 0 + 55 54 0 0 + 56 55 0 0 + 57 56 0 0 + 58 57 0 0 + 59 58 0 0 + 60 59 0 0 + 61 60 0 0 + 62 61 0 0 + 63 62 0 0 + 64 63 0 0 + 65 64 0 0 + 66 65 0 0 + 67 66 0 0 + 68 67 0 0 + 69 68 0 0 + 70 69 0 0 + 71 70 0 0 + 72 71 0 0 + 73 72 0 0 + 74 73 0 0 + 75 74 0 0 + 76 75 0 0 + 77 76 0 0 + 78 77 0 0 + 79 78 0 0 + 80 79 0 0 + 81 80 0 0 + 82 81 0 0 + 83 82 0 0 + 84 83 0 0 + 85 84 0 0 + 86 85 0 0 + 87 86 0 0 + 88 87 0 0 + 89 88 0 0 + 90 89 0 0 + 91 90 0 0 + 92 91 0 0 + 93 92 0 0 + 94 93 0 0 + 95 94 0 0 + 96 95 0 0 + 97 96 0 0 + 98 97 0 0 + 99 98 0 0 + 100 99 0 0 + 101 100 0 0 + + # references des noeuds + + N_deb 1 + N_fi 101 + 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 67 68 69 70 71 72 73 74 75 76 77 78 79 80 + 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 + 97 98 99 100 101 + + + # les elements + + elements ---------- + 100 ELEMENTS + #-------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #-------------------------------------------------------------------- + 1 POUT BIE1 1 2 + 2 POUT BIE1 2 3 + 3 POUT BIE1 3 4 + 4 POUT BIE1 4 5 + 5 POUT BIE1 5 6 + 6 POUT BIE1 6 7 + 7 POUT BIE1 7 8 + 8 POUT BIE1 8 9 + 9 POUT BIE1 9 10 + 10 POUT BIE1 10 11 + 11 POUT BIE1 11 12 + 12 POUT BIE1 12 13 + 13 POUT BIE1 13 14 + 14 POUT BIE1 14 15 + 15 POUT BIE1 15 16 + 16 POUT BIE1 16 17 + 17 POUT BIE1 17 18 + 18 POUT BIE1 18 19 + 19 POUT BIE1 19 20 + 20 POUT BIE1 20 21 + 21 POUT BIE1 21 22 + 22 POUT BIE1 22 23 + 23 POUT BIE1 23 24 + 24 POUT BIE1 24 25 + 25 POUT BIE1 25 26 + 26 POUT BIE1 26 27 + 27 POUT BIE1 27 28 + 28 POUT BIE1 28 29 + 29 POUT BIE1 29 30 + 30 POUT BIE1 30 31 + 31 POUT BIE1 31 32 + 32 POUT BIE1 32 33 + 33 POUT BIE1 33 34 + 34 POUT BIE1 34 35 + 35 POUT BIE1 35 36 + 36 POUT BIE1 36 37 + 37 POUT BIE1 37 38 + 38 POUT BIE1 38 39 + 39 POUT BIE1 39 40 + 40 POUT BIE1 40 41 + 41 POUT BIE1 41 42 + 42 POUT BIE1 42 43 + 43 POUT BIE1 43 44 + 44 POUT BIE1 44 45 + 45 POUT BIE1 45 46 + 46 POUT BIE1 46 47 + 47 POUT BIE1 47 48 + 48 POUT BIE1 48 49 + 49 POUT BIE1 49 50 + 50 POUT BIE1 50 51 + 51 POUT BIE1 51 52 + 52 POUT BIE1 52 53 + 53 POUT BIE1 53 54 + 54 POUT BIE1 54 55 + 55 POUT BIE1 55 56 + 56 POUT BIE1 56 57 + 57 POUT BIE1 57 58 + 58 POUT BIE1 58 59 + 59 POUT BIE1 59 60 + 60 POUT BIE1 60 61 + 61 POUT BIE1 61 62 + 62 POUT BIE1 62 63 + 63 POUT BIE1 63 64 + 64 POUT BIE1 64 65 + 65 POUT BIE1 65 66 + 66 POUT BIE1 66 67 + 67 POUT BIE1 67 68 + 68 POUT BIE1 68 69 + 69 POUT BIE1 69 70 + 70 POUT BIE1 70 71 + 71 POUT BIE1 71 72 + 72 POUT BIE1 72 73 + 73 POUT BIE1 73 74 + 74 POUT BIE1 74 75 + 75 POUT BIE1 75 76 + 76 POUT BIE1 76 77 + 77 POUT BIE1 77 78 + 78 POUT BIE1 78 79 + 79 POUT BIE1 79 80 + 80 POUT BIE1 80 81 + 81 POUT BIE1 81 82 + 82 POUT BIE1 82 83 + 83 POUT BIE1 83 84 + 84 POUT BIE1 84 85 + 85 POUT BIE1 85 86 + 86 POUT BIE1 86 87 + 87 POUT BIE1 87 88 + 88 POUT BIE1 88 89 + 89 POUT BIE1 89 90 + 90 POUT BIE1 90 91 + 91 POUT BIE1 91 92 + 92 POUT BIE1 92 93 + 93 POUT BIE1 93 94 + 94 POUT BIE1 94 95 + 95 POUT BIE1 95 96 + 96 POUT BIE1 96 97 + 97 POUT BIE1 97 98 + 98 POUT BIE1 98 99 + 99 POUT BIE1 99 100 + 100 POUT BIE1 100 101 + + + + # references des arretes + + A_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 51 1 52 1 53 1 54 1 55 1 56 1 57 1 58 1 59 1 60 1 61 1 62 1 63 1 64 1 + 65 1 66 1 67 1 68 1 69 1 70 1 71 1 72 1 73 1 74 1 75 1 76 1 77 1 78 1 79 1 80 1 + 81 1 82 1 83 1 84 1 85 1 86 1 87 1 88 1 89 1 90 1 91 1 92 1 93 1 94 1 95 1 96 1 + 97 1 98 1 99 1 100 1 + + + # references des elements + + 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 51 52 53 54 55 56 57 58 59 60 61 62 63 64 + 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 + 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 + 97 98 99 100 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.CVisu1 b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.CVisu1 new file mode 100644 index 0000000..433f9ee --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.CVisu1 @@ -0,0 +1,274 @@ +############################################################################################### +# Fichier de commande pour la visualisation elements finis # +# Herezh++ V7.036 # +# Copyright (C) 1997-2022 Université Bretagne Sud (France), AUTHOR : Gérard Rio (gerardrio56@free.fr)/ # +# 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 2 fin_choix_maillage + + # ----------------------------- definition des grandeurs a visualiser (maple): ---------------- + debut_grandeurs_maple # un mot cle de debut (debut_grandeurs_maple), + actif 1 # <0 ou 1> indique si l'ordre est actif ou non + # les grandeurs globales (energie, puissance ....) un mot cle de debut, le nom des grandeurs puis un mot de fin + debut_grandeurs_globales compteur_iteration_algo_global 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 de sommes, moyennes etc. sur des références + # - de noeuds + # - d'elements + # - de face d'element + # - d'arete d'element + # 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 + # les infos pour la visualisation eventuelle aux faces d'element de grandeurs quelconques, + # les infos pour la visualisation eventuelle aux aretes d'element de grandeurs quelconques, + # 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_deb 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 + # les SM sur ref de faces d'element: un mot cle de debut: debut_list_SM_sur_refFE + # une liste de nom de ref de faces d'element , le pti associe + # puis + debut_list_SM_sur_refFE fin_list__SM_sur_refFE + # les SMs sur ref de pti de face: un mot cle de debut: debut_list_SM_sur_refFpti + # une liste de nom de ref , puis + debut_list_SM_sur_refFpti fin_list__SM_sur_refFpti + # debut de la liste des grandeurs a considerer pour les SM sur face, (une liste de string), , (la liste), puis + deb_list_Grandeur_SM_sur_refFEpti fin_list_Grandeur_SM_sur_refFEpti # fin des grandeurs aux _SM_sur_refFEptis + # les SM sur ref d'aretes d'element: un mot cle de debut: debut_list_SM_sur_refAE + # une liste de nom de ref d'aretes d'element , le pti associe + # puis + debut_list_SM_sur_refAE fin_list__SM_sur_refAE + # les SMs sur ref de pti d'aretes: un mot cle de debut: debut_list_SM_sur_refApti + # une liste de nom de ref , puis + debut_list_SM_sur_refApti fin_list__SM_sur_refApti + # debut de la liste des grandeurs a considerer pour les SM sur aretes, (une liste de string), , (la liste), puis + deb_list_Grandeur_SM_sur_refAEpti fin_list_Grandeur_SM_sur_refAEpti # fin des grandeurs aux _SM_sur_refAEptis + debut_liste_ddl_et_noeud # ** debut des ddl principaux aux noeuds + # debut de la liste de noeuds, puis une liste de numero de noeud , puis + deb_list_noeud fin_list_noeud + # debut de la liste des ref de noeuds, puis une liste de nom , puis + deb_list_ref_noeud N_deb N_fi fin_list_ref_noeud + # debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis + deb_list_ddl_noeud X1 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 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 EPS11 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 + debut_liste_grandeurs_ptinteg_F_E # ** debut des grandeurs aux faces d'elements + # debut de la liste des faces d'element et points d'integration, une liste de (un element, une face, un numero de pt d'integ), puis + deb_list_NbFaceElement_NbPtInteg fin_list_NbFaceElement_NbPtInteg + # debut de la liste des ref de faces d'element, puis une liste de: nom + numero d'integ, puis + deb_list_ref_F_element fin_list_ref_F_element + # debut de la liste des ref de ptinteg de faces d'element, puis une liste de: nom puis + deb_list_ref_ptinteg_F_element fin_list_ref_ptinteg_F_element + # debut de la liste des grandeurs a considerer pour les faces d'elements, (une liste de string), puis + deb_list_Grandeur_F_element fin_list_Grandeur_F_element + fin_liste_grandeurs_ptinteg_F_E # fin des grandeurs aux faces d'elements + debut_liste_grandeurs_ptinteg_A_E # ** debut des grandeurs aux aretes d'elements + # debut de la liste des aretes d'element et points d'integration, une liste de (un element, une arete, un numero de pt d'integ), puis + deb_list_NbAreteElement_NbPtInteg fin_list_NbAreteElement_NbPtInteg + # debut de la liste des ref d'arete d'element, puis une liste de: nom + numero d'integ, puis + deb_list_ref_A_element fin_list_ref_A_element + # debut de la liste des ref de ptinteg d'aretes d'element, puis une liste de: nom puis + deb_list_ref_ptinteg_A_element fin_list_ref_ptinteg_A_element + # debut de la liste des grandeurs a considerer pour les aretes d'elements, (une liste de string), puis + deb_list_Grandeur_A_element fin_list_Grandeur_A_element + fin_liste_grandeurs_ptinteg_A_E # fin des grandeurs aux aretes d'elements + 2 # le numero de maillage + # les torseurs de reaction: un mot cle de debut: debut_list_torseur_reaction + # une liste de nom , puis + debut_list_torseur_reaction fin_list_torseur_reaction + # 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 + # les SM sur ref de faces d'element: un mot cle de debut: debut_list_SM_sur_refFE + # une liste de nom de ref de faces d'element , le pti associe + # puis + debut_list_SM_sur_refFE fin_list__SM_sur_refFE + # les SMs sur ref de pti de face: un mot cle de debut: debut_list_SM_sur_refFpti + # une liste de nom de ref , puis + debut_list_SM_sur_refFpti fin_list__SM_sur_refFpti + # debut de la liste des grandeurs a considerer pour les SM sur face, (une liste de string), , (la liste), puis + deb_list_Grandeur_SM_sur_refFEpti fin_list_Grandeur_SM_sur_refFEpti # fin des grandeurs aux _SM_sur_refFEptis + # les SM sur ref d'aretes d'element: un mot cle de debut: debut_list_SM_sur_refAE + # une liste de nom de ref d'aretes d'element , le pti associe + # puis + debut_list_SM_sur_refAE fin_list__SM_sur_refAE + # les SMs sur ref de pti d'aretes: un mot cle de debut: debut_list_SM_sur_refApti + # une liste de nom de ref , puis + debut_list_SM_sur_refApti fin_list__SM_sur_refApti + # debut de la liste des grandeurs a considerer pour les SM sur aretes, (une liste de string), , (la liste), puis + deb_list_Grandeur_SM_sur_refAEpti fin_list_Grandeur_SM_sur_refAEpti # fin des grandeurs aux _SM_sur_refAEptis + 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 fin_list_ref_element + # debut de la liste des ref de ptinteg d'elements, puis une liste de: nom puis + deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element + # debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis + deb_list_ddl_element fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements + # debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis + deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements + # tableau de grandeurs evoluees aux elements a visualiser, un par maillage + deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element + fin_liste_ddl_ptinteg # fin des grandeurs aux elements + debut_liste_grandeurs_ptinteg_F_E # ** debut des grandeurs aux faces d'elements + # debut de la liste des faces d'element et points d'integration, une liste de (un element, une face, un numero de pt d'integ), puis + deb_list_NbFaceElement_NbPtInteg fin_list_NbFaceElement_NbPtInteg + # debut de la liste des ref de faces d'element, puis une liste de: nom + numero d'integ, puis + deb_list_ref_F_element fin_list_ref_F_element + # debut de la liste des ref de ptinteg de faces d'element, puis une liste de: nom puis + deb_list_ref_ptinteg_F_element fin_list_ref_ptinteg_F_element + # debut de la liste des grandeurs a considerer pour les faces d'elements, (une liste de string), puis + deb_list_Grandeur_F_element fin_list_Grandeur_F_element + fin_liste_grandeurs_ptinteg_F_E # fin des grandeurs aux faces d'elements + debut_liste_grandeurs_ptinteg_A_E # ** debut des grandeurs aux aretes d'elements + # debut de la liste des aretes d'element et points d'integration, une liste de (un element, une arete, un numero de pt d'integ), puis + deb_list_NbAreteElement_NbPtInteg fin_list_NbAreteElement_NbPtInteg + # debut de la liste des ref d'arete d'element, puis une liste de: nom + numero d'integ, puis + deb_list_ref_A_element fin_list_ref_A_element + # debut de la liste des ref de ptinteg d'aretes d'element, puis une liste de: nom puis + deb_list_ref_ptinteg_A_element fin_list_ref_ptinteg_A_element + # debut de la liste des grandeurs a considerer pour les aretes d'elements, (une liste de string), puis + deb_list_Grandeur_A_element fin_list_Grandeur_A_element + fin_liste_grandeurs_ptinteg_A_E # fin des grandeurs aux aretes d'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/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.info b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.info new file mode 100644 index 0000000..aadae08 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.info @@ -0,0 +1,177 @@ + +####################################################################### +# # +# | | ==== === ==== ==== | | | | # +# | | | | | | / | | | | # +# |====| |=== === |=== / |====| ------- ------- # +# | | | | \ | / | | | | # +# | | ==== | \ ==== ==== | | | | # +# # +####################################################################### +# Herezh++ is distributed under GPL 3 license ou ultérieure. # +# Copyright (C) 1997-2022 Université Bretagne Sud (France) # +# AUTHOR : Gérard Rio # +# E-MAIL : gerardrio56@free.fr # +# 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 7.036 + +####################################################################### +# fichier de commande # + +# --- test académique : contact 2D entre une barre composée de 100 biellettes et un quadrangle +# on compare avec un calcul théorique + +#---------------------------------------------------- +# definition de la dimension de l'espace de travail | +#---------------------------------------------------- + dimension 2 + +#--------------------------------------------------------------- +# definition facultative du niveau d'impression (entre 0 et 10)| +#--------------------------------------------------------------- + niveau_commentaire 3#0 + +#------------------------------- +# definition du type de calcul | +#------------------------------- + TYPE_DE_CALCUL + +non_dynamique #avec plus visualisation +#-------------------------------------- +#| definition du (ou des) maillage(s) | +#-------------------------------------- +#< biel.her +< biel100.her + +#< rec100x100.her +#< rec_tri_100x100.her +< rec_tri_50x100_1.her +#< rec_tri_50x100.her + + def_mouvement_solide_initiaux_ + mouvement_solide_ + translation_= 100.000 -25. 0. + centre_= 0. 0. 0. + rotation_= 0. 0. 0. + fin_mouvement_solide_ + + domaine_esclave +#-------------------------------------------------------------------- +# definition du nombre de domaine esclave | +#-------------------------------------------------------------------- +1 + +#----------------------------- +# Definition des courbes +#----------------------------- + les_courbes_1D + +charge_triangle COURBEPOLYLINEAIRE_1_D + Debut_des_coordonnees_des_points + Coordonnee dim= 2 0. 0. + Coordonnee dim= 2 1 1. + Coordonnee dim= 2 2. 0. + Fin_des_coordonnees_des_points + +#====================================== +#| definition des lois de comportement| +#-------------------------------------- + + choix_materiaux ---------- +#---------------------------------------- +# Elements | Nom Materiau | +#---------------------------------------- +nom_mail= biel E_tout acier_1D +nom_mail= rec100x100 E_tout acier + + materiaux #---------- + +#------------------------------------------ +# Nom Materiau | Type loi | +#------------------------------------------ + acier_1D ISOELAS1D +# module d'young : coefficient de poisson + 200000 0. + + acier ISOELAS2D_C +# module d'young : coefficient de poisson + 200000 0.3 + +#--------------------------------------- fin def des lois de comportement ------------------------ + +# --- divers stockages (1) ------- + masse_volumique #--------# +nom_mail= biel E_tout 8.9e-12 # en tonne/mm3 +nom_mail= rec100x100 E_tout 8.9e-9 # en tonne/mm3 + sections ------------ +nom_mail= biel E_tout 50. + variation_section # pas de variation de section +nom_mail= biel E_tout 0 + epaisseurs ------------ +nom_mail= rec100x100 E_tout 1. + + charges #------------# + + blocages #------------# +#----------------------------------------------------------- +# nom du maillage | Ref noeud | Bloquages +#----------------------------------------------------------- +# la biel: dep vertical bloqué +nom_mail= biel N_tout UY + +# déplacement maxi suivant x de 0.3 mm -> environ 3% +nom_mail= biel N_deb 'UX = COURBE_CHARGE: charge_triangle ECHELLE: 0.3' + +# bloc inférieur bloqué +nom_mail= rec100x100 N_tout UY UX + + controle ------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +TEMPSFIN 1.#0.1#1.#2.#3. +#DELTAt 0.1 +SAUVEGARDE 1#INTER_TEMPS 0.1 +PRECISION 1.e-4 + + para_contact ------------ +#------------------------------------------ +# PARAMETRE | VALEUR | +#------------------------------------------ +CONTACT_TYPE 2 +PENALISATION_PENETRATION 10.#0.9 +TYPE_PENALISATION_PENETRATION 2 +NB_DECOLLEMENT_MAXI 1 +DISTANCE_MAXI_AU_PT_PROJETE 0.1 +TYPE_DE_DECOLLEMENT 0 + +para_pilotage_equi_global ------------ +#----------------------------------------------- +# PARAMETRE | VALEUR | +#----------------------------------------------- +FACTEUR_DIMINUTION 1.732 +FACTEUR_AUGMENTATION 1.414 +NB_BONNE_CONVERGENCE 3 +NORME_MAXI_INCREMENT 0.1 + + para_affichage #------------ +#--------------------------- +# PARAMETRE | VALEUR | +#--------------------------- +FREQUENCE_AFFICHAGE_ITERATION 1#00 +FREQUENCE_AFFICHAGE_INCREMENT 1#00 +FREQUENCE_SORTIE_FIL_DU_CALCUL 1 + +# ------------------------------------------------------------------------------------ + resultats pas_de_sortie_finale_ + COPIE 0 +# + + + _fin_point_info_ + + + \ No newline at end of file diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact_princ.maple.ref1 b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact_princ.maple.ref1 new file mode 100644 index 0000000..61c1636 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact_princ.maple.ref1 @@ -0,0 +1,156 @@ +#fichier au format maple6 +############################################################################################### +# Visualisation elements finis : Herezh++ V7.036 # +# Copyright (C) 1997-2022 Université Bretagne Sud (France), AUTHOR : Gérard Rio (gerardrio56@free.fr)/ # +# 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 + # + # ensuite pour les moyennes, sommes, maxi etc. calculees sur des references de faces d'element ou de pti + # >> le nombre de ref de face + ref de pti de face, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference de face et de pti de face + # le nom de la reference, le nombre de face 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'arete d'element ou de pti + # >> le nombre de ref d'arete + ref de pti d'arete, nombre total de grandeurs associees + # puis pour chaque maillage + # le numero de maillage puis pour chaque reference d'arete et de pti d'arete + # le nom de la reference, le nombre d'arete 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 + # puis pour chaque maillage + # >> le nombre de triplets element-nb_face-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # puis le nombre grandeurs, suivi des identificateurs+categorie+type (chaines de caracteres), + # suivi entre crochet, de la plage des numeros de colonnes, correspondant + # chacun sur une ligne differentes + # puis pour chaque maillage + # >> le nombre de triplets element-nb_arete-pt_integ (peut etre nulle) ou il y a des grandeurs en sortie , + # puis le nombre grandeurs, suivi des 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 ------------------------------------ +#1 (nombre de grandeurs globales) [2]compteur_iteration_algo_global +#---------------------------------- maillage et dimension -------------------------------- +#2 2 (nombre de maillages et dimension) +#---------------------------------- torseurs de reactions -------------------------------- +# maillage_1 : 1 3 (nombre de torseurs et nombre total de grandeurs associees) +# N_deb [4]Rx [5]Ry [6]Mz ; +# +# maillage_2 : 0 0 (nombre de torseurs et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux noeuds pour des ref --------------- +# maillage_1 : 0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +# maillage_2 : 0 0 (nombre de ref de noeud et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux elements pour des ref --------------- +# maillage_1 : 0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +# maillage_2 : 0 0 (nombre de ref d'element et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux faces d'elements pour des ref --------------- +# maillage_1 : 0 0 (nombre de ref de faces d'element et nombre total de grandeurs associees) +# +# maillage_2 : 0 0 (nombre de ref de faces d'element et nombre total de grandeurs associees) +# +#-------- moyenne, somme, maxi etc. de grandeurs aux aretes d'element pour des ref --------------- +# maillage_1 : 0 0 (nombre de ref d'arete d'element et nombre total de grandeurs associees) +# +# maillage_2 : 0 0 (nombre de ref d'arete d'element et nombre total de grandeurs associees) +# +#---------------------------------- grandeurs aux noeuds -------------------------------- +#2 1 (nombre de noeuds, nombre total de grandeurs associees) +# maillage_1 : +# reference N_deb : noeud_1 [8]X [9]Y [10]delta_0_t_X1 ; +# reference N_fi : noeud_101 [11]X [12]Y [13]delta_0_t_X1 ; +# +# maillage_2 : 0 0 (nombre de noeuds, nombre total de grandeurs associees) +#---------------------------------- grandeurs aux elements ------------------------------ +#1 2 00 (nombre total d'elements, nombre totale de grandeurs associees, nombre de grandeurs particulieres, nombre de grandeurs tensorielles) +# maillage_1 : +# element_1 pt_integ_1: [15]X [16]Y [17] SIG11 [18] EPS11 +# +# +# maillage_2 : 0 0 (nombre total d'elements, nombre totale de grandeurs associees) +#---------------------------------- grandeurs aux faces d'elements ------------------------------ +# maillage_1 : +#0 0 (nombre total de faces d'elements, nombre totale de grandeurs associees) +# maillage_2 : +#0 0 (nombre total de faces d'elements, nombre totale de grandeurs associees) +#---------------------------------- grandeurs aux aretes d'elements ------------------------------ +# maillage_1 : +#0 0 (nombre total d'aretes d'elements, nombre totale de grandeurs associees) +# maillage_2 : +#0 0 (nombre total d'aretes d'elements, nombre totale de grandeurs associees) +#==================================================================== +#|| fin du recapitulatif des differentes grandeurs || +#==================================================================== + + # ensuite les donnees sont organisees sur differentes lignes, chaques lignes correspondant + # a un calcul (par exemple un pas de temps), sur chaque ligne il y a m enregistrement, chacun + # correspondant a un maillage. On trouve pour chaque enregistrement successivement : + # s'il y a des grandeurs globales: le temps puis les grandeurs globales, + # puis s'il y a des torseurs de reaction : + # de nouveau le temps, les composantes de la resultante puis les composantes du moments + # donc en 1D -> 1 reels (resultante), en 2D -> 3 reels (resultante 2, moment 1) et en 3D 6 reels + # puis s'il y a des grandeurs aux noeuds: de nouveau le temps + # les coordonnees a t du premier noeud suivi des p1 grandeurs correspondant au premier noeud + # puis les coordonnees du second noeud, les p1 grandeurs etc. pour tous les noeuds + # puis s'il y a des grandeur aux elements: + # le temps, puis les coordonnees a t du point d'integration d'un element (pour les grandeurs generiques) + # suivi des p2 grandeurs correspondantes puis les coordonnees a t du point d'integration + # correspondant aux grandeurs specifiques suivi des p3 grandeurs correspondantes + # puis les coordonnees d'un second point d'integration d'un element, les p2 grandeurs + # etc. pour tous les points d'integration - element + + 1.000000000000e+00 5 1.000000000000e+00 -3.010022585416e+04 0.000000000000e+00 0.000000000000e+00 1.000000000000e+00 3.000000000000e-01 0.000000000000e+00 3.000000000000e-01 1.000003611890e+02 0.000000000000e+00 3.611890172692e-04 1.000000000000e+00 7.985017499735e-01 0.000000000000e+00 -6.020045170833e+02 -3.010022585416e-03 diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec100x100.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec100x100.her new file mode 100644 index 0000000..3c66a36 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec100x100.her @@ -0,0 +1,77 @@ + +############################################################################# +# # +# 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 rec100x100 + + 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 + N_contact 1 2 + + + # 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/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_100x100.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_100x100.her new file mode 100644 index 0000000..0a1340d --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_100x100.her @@ -0,0 +1,84 @@ + +############################################################################# +# # +# 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 triangulaire, interpolation lineaire. + + nom_maillage rec100x100 + + 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 + + N_contact 1 2 + + + + # les elements + + elements ---------- + 2 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 TRIANGLE LINEAIRE 1 3 2 + 2 TRIANGLE LINEAIRE 3 4 2 + + + # references des elements + +A_O 1 3 +A_S 1 1 +A_E 2 1 +A_N 2 2 +F_to 1 1 2 1 +E_O 1 +E_S 1 +E_E 2 +E_N 2 +E_SO 1 +E_SE 1 2 +E_NE 2 +E_NO 1 2 +F_O 1 1 2 1 +F_S 1 1 2 1 +F_E 1 1 2 1 +F_N 1 1 2 1 + E_to 1 2 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100.her new file mode 100644 index 0000000..400c169 --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100.her @@ -0,0 +1,122 @@ + +############################################################################# +# # +# 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 50 ( 10 x 1 elements ) + # geometrie rectangulaire, decoupage triangulaire, interpolation lineaire. + + nom_maillage rec100x100 + + noeuds ------------ + 22 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 50 0 0 + 2 0 0 0 + 3 50 10 0 + 4 0 10 0 + 5 50 20 0 + 6 0 20 0 + 7 50 30 0 + 8 0 30 0 + 9 50 40 0 + 10 0 40 0 + 11 50 50 0 + 12 0 50 0 + 13 50 60 0 + 14 0 60 0 + 15 50 70 0 + 16 0 70 0 + 17 50 80 0 + 18 0 80 0 + 19 50 90 0 + 20 0 90 0 + 21 50 100 0 + 22 0 100 0 + + # references des noeuds + +N_O 2 4 6 8 10 12 14 16 18 20 22 +N_S 1 2 +N_E 1 3 5 7 9 11 13 15 17 19 21 +N_N 21 22 +N_SO 2 +N_SE 1 +N_NE 21 +N_NO 22 + N_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 + 17 18 19 20 21 22 + + + # les elements + + elements ---------- + 20 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 TRIANGLE LINEAIRE 1 3 2 + 2 TRIANGLE LINEAIRE 3 4 2 + 3 TRIANGLE LINEAIRE 3 5 4 + 4 TRIANGLE LINEAIRE 5 6 4 + 5 TRIANGLE LINEAIRE 5 7 6 + 6 TRIANGLE LINEAIRE 7 8 6 + 7 TRIANGLE LINEAIRE 7 9 8 + 8 TRIANGLE LINEAIRE 9 10 8 + 9 TRIANGLE LINEAIRE 9 11 10 + 10 TRIANGLE LINEAIRE 11 12 10 + 11 TRIANGLE LINEAIRE 11 13 12 + 12 TRIANGLE LINEAIRE 13 14 12 + 13 TRIANGLE LINEAIRE 13 15 14 + 14 TRIANGLE LINEAIRE 15 16 14 + 15 TRIANGLE LINEAIRE 15 17 16 + 16 TRIANGLE LINEAIRE 17 18 16 + 17 TRIANGLE LINEAIRE 17 19 18 + 18 TRIANGLE LINEAIRE 19 20 18 + 19 TRIANGLE LINEAIRE 19 21 20 + 20 TRIANGLE LINEAIRE 21 22 20 + + + # references des elements + +A_O 2 2 4 2 6 2 8 2 10 2 12 2 14 2 16 2 18 2 20 2 +A_S 1 3 +A_E 1 1 3 1 5 1 7 1 9 1 11 1 13 1 15 1 17 1 19 1 +A_N 20 1 +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 +E_O 2 4 6 8 10 12 14 16 18 20 +E_S 1 +E_E 1 3 5 7 9 11 13 15 17 19 +E_N 20 +E_SO 1 2 +E_SE 1 +E_NE 19 20 +E_NO 20 +F_O 1 1 2 1 +F_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 13 1 14 1 15 1 16 1 17 1 18 1 19 1 20 1 +F_E 10 1 11 1 +F_N 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 + E_to 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 + 17 18 19 20 + diff --git a/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100_1.her b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100_1.her new file mode 100644 index 0000000..eee3eff --- /dev/null +++ b/Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100_1.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 : 50 x 100 ( 1 x 1 elements ) + # geometrie rectangulaire, decoupage triangulaire, interpolation lineaire. + + nom_maillage rec100x100 + + noeuds ------------ + 4 NOEUDS + + +#--------------------------------------------------------------- +#|NO DU| X | Y | Z | +#|NOEUD| | | | +#--------------------------------------------------------------- + 1 0 0 0 + 2 0 100 0 + 3 50 0 0 + 4 50 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 ---------- + 2 ELEMENTS + #---------------------------------------------------------------------- + #| NO | | | + #|ELTS | type element | Noeuds | + #---------------------------------------------------------------------- + 1 TRIANGLE LINEAIRE 1 3 2 + 2 TRIANGLE LINEAIRE 3 4 2 + + + # references des elements + +A_O 1 3 +A_S 1 1 +A_E 2 1 +A_N 2 2 +F_to 1 1 2 1 +E_O 1 +E_S 1 +E_E 2 +E_N 2 +E_SO 1 +E_SE 1 2 +E_NE 2 +E_NO 1 2 +F_O 1 1 2 1 +F_S 1 1 2 1 +F_E 1 1 2 1 +F_N 1 1 2 1 + E_to 1 2 +