rio/Verif_Herezh
2
0
Fork 0
Code Releases Activity

ajout test académique entre une barre constituée de 100 biellettes et

Browse source 
un massif constitué de 20 triangles, le tout en 2D.
This commit is contained in:
Gérard Rio 2024-10-11 12:23:06 +02:00
parent 3a5500a180
commit 1b13ae9dba
10 changed files with 1381 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

72
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/README Normal file
View file

@ -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

61
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel.her Normal file
View file

@ -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

277
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel100.her Normal file
View file

@ -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

274
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.CVisu1 Normal file
View file

@ -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: <debut_fichier_commande_visu>
# permet au programme de se positionner au debut du fichier, il est indispensable
# =================================================================================
# || ***** demande d'une visualisation maple: ***** ||
# =================================================================================
# un mot cle de debut (debut_visualisation_maple)
# un mot cle de fin ( fin_visualisation_maple)
# la seule presence de ces deux mots cle suffit a activer la visualisation maple
debut_visualisation_maple
# ----------------------------- definition de la liste des increments a balayer: ----------------
debut_list_increment # un mot cle de debut de liste
actif 1 # <0 ou 1> indique si l'ordre est actif ou non
# une liste d'entier separee par des blancs, ou le mot cle (tous_les_increments)
# un mot cle de fin de liste ( fin_list_increment)
fin_list_increment
# ----------------------------- choix des maillages a visualiser: ----------------
# la liste est facultative, par defaut la visualisation concerne le premier maillage
debut_choix_maillage # un mot cle de debut,
actif 0 # <0 ou 1> indique si l'ordre est actif ou non
# une liste d'entiers , puis <fin_choix_maillage>, sur une meme ligne
1 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 <un entier>,
# 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 <chaine de caracteres >, puis <fin_list_torseur_reaction>
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 <chaine de caracteres >, puis <fin_list__SM_sur_refN>
debut_list_SM_sur_refN fin_list__SM_sur_refN
# debut de la liste des ddl a considerer <deb_list_ddl_SM_sur_refN>, (une liste de ddl), puis <fin_list_ddl_SM_sur_refN>
deb_list_ddl_SM_sur_refN fin_list_ddl_SM_sur_refN
# debut de la liste des ddl etendus a considerer, (une liste de ddl), puis <fin_list_ddl_etendu__SM_sur_refN>
deb_list_ddl_etendu__SM_sur_refN fin_list_ddl_etendu__SM_sur_refN # fin des ddl etendu pour _SM_sur_refNs
# debut de la liste des grandeurs particulieres a considerer, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refN>
deb_list_GrandParticuliere__SM_sur_refN fin_list_GrandParticuliere__SM_sur_refN # fin des grandeurs particulieres aux _SM_sur_refNs
# les SM sur ref d'elements: un mot cle de debut: debut_list_SM_sur_refE
# une liste de nom de ref d'elements <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refE>
debut_list_SM_sur_refE fin_list__SM_sur_refE
# les SMs sur ref de pti: un mot cle de debut: debut_list_SM_sur_refpti
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refpti>
debut_list_SM_sur_refpti fin_list__SM_sur_refpti
# debut de la liste des ddl a considerer , (une liste de ddl), puis <fin_list_ddl__SM_sur_refEpti>
deb_list_ddl__SM_sur_refEpti fin_list_ddl__SM_sur_refEpti # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les _SM_sur_refEptis, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refEpti>
deb_list_GrandParticuliere__SM_sur_refEpti fin_list_GrandParticuliere__SM_sur_refEpti # fin des grandeurs particulieres aux _SM_sur_refEptis
# tableau de grandeurs evoluees aux _SM_sur_refEptis a visualiser, un par maillage
deb_list_GrandEvoluee__SM_sur_refEpti fin_list_GrandEvoluee__SM_sur_refEpti
# 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 <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refFE>
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 <chaine de caracteres >, puis <fin_list__SM_sur_refFpti>
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), <deb_list_Grandeur_SM_sur_refFEpti> , (la liste), puis <fin_list_Grandeur_SM_sur_refFEpti>
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 <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refAE>
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 <chaine de caracteres >, puis <fin_list__SM_sur_refApti>
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), <deb_list_Grandeur_SM_sur_refAEpti> , (la liste), puis <fin_list_Grandeur_SM_sur_refAEpti>
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 <entier>, puis <fin_list_noeud>
deb_list_noeud fin_list_noeud
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
deb_list_ref_noeud N_deb N_fi fin_list_ref_noeud
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
deb_list_ddl_noeud 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 <fin_list_ddl_etendu_noeud>
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
# debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis <fin_list_NbElement_NbPtInteg>
deb_list_NbElement_NbPtInteg 1 1 fin_list_NbElement_NbPtInteg
# debut de la liste des ref d'elements, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_element>
deb_list_ref_element fin_list_ref_element
# debut de la liste des ref de ptinteg d'elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
deb_list_ddl_element 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 <fin_list_GrandParticuliere_element>
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
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 <fin_list_NbFaceElement_NbPtInteg>
deb_list_NbFaceElement_NbPtInteg fin_list_NbFaceElement_NbPtInteg
# debut de la liste des ref de faces d'element, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_F_element>
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 <chaine de caracteres > puis <fin_list_ref_ptinteg_F_element>
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 <fin_list_Grandeur_F_element>
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 <fin_list_NbAreteElement_NbPtInteg>
deb_list_NbAreteElement_NbPtInteg fin_list_NbAreteElement_NbPtInteg
# debut de la liste des ref d'arete d'element, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_A_element>
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 <chaine de caracteres > puis <fin_list_ref_ptinteg_A_element>
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 <fin_list_Grandeur_A_element>
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 <chaine de caracteres >, puis <fin_list_torseur_reaction>
debut_list_torseur_reaction fin_list_torseur_reaction
# les sommes, moyennes etc. sur ref de noeuds: un mot cle de debut: debut_list_SM_sur_refN
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refN>
debut_list_SM_sur_refN fin_list__SM_sur_refN
# debut de la liste des ddl a considerer <deb_list_ddl_SM_sur_refN>, (une liste de ddl), puis <fin_list_ddl_SM_sur_refN>
deb_list_ddl_SM_sur_refN fin_list_ddl_SM_sur_refN
# debut de la liste des ddl etendus a considerer, (une liste de ddl), puis <fin_list_ddl_etendu__SM_sur_refN>
deb_list_ddl_etendu__SM_sur_refN fin_list_ddl_etendu__SM_sur_refN # fin des ddl etendu pour _SM_sur_refNs
# debut de la liste des grandeurs particulieres a considerer, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refN>
deb_list_GrandParticuliere__SM_sur_refN fin_list_GrandParticuliere__SM_sur_refN # fin des grandeurs particulieres aux _SM_sur_refNs
# les SM sur ref d'elements: un mot cle de debut: debut_list_SM_sur_refE
# une liste de nom de ref d'elements <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refE>
debut_list_SM_sur_refE fin_list__SM_sur_refE
# les SMs sur ref de pti: un mot cle de debut: debut_list_SM_sur_refpti
# une liste de nom de ref <chaine de caracteres >, puis <fin_list__SM_sur_refpti>
debut_list_SM_sur_refpti fin_list__SM_sur_refpti
# debut de la liste des ddl a considerer , (une liste de ddl), puis <fin_list_ddl__SM_sur_refEpti>
deb_list_ddl__SM_sur_refEpti fin_list_ddl__SM_sur_refEpti # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les _SM_sur_refEptis, (une liste de string), puis <fin_list_GrandParticuliere__SM_sur_refEpti>
deb_list_GrandParticuliere__SM_sur_refEpti fin_list_GrandParticuliere__SM_sur_refEpti # fin des grandeurs particulieres aux _SM_sur_refEptis
# tableau de grandeurs evoluees aux _SM_sur_refEptis a visualiser, un par maillage
deb_list_GrandEvoluee__SM_sur_refEpti fin_list_GrandEvoluee__SM_sur_refEpti
# 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 <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refFE>
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 <chaine de caracteres >, puis <fin_list__SM_sur_refFpti>
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), <deb_list_Grandeur_SM_sur_refFEpti> , (la liste), puis <fin_list_Grandeur_SM_sur_refFEpti>
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 <chaine de caracteres >, le pti associe
# puis <fin_list__SM_sur_refAE>
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 <chaine de caracteres >, puis <fin_list__SM_sur_refApti>
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), <deb_list_Grandeur_SM_sur_refAEpti> , (la liste), puis <fin_list_Grandeur_SM_sur_refAEpti>
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 <entier>, puis <fin_list_noeud>
deb_list_noeud fin_list_noeud
# debut de la liste des ref de noeuds, puis une liste de nom <chaine de caracteres >, puis <fin_list_ref_noeud>
deb_list_ref_noeud fin_list_ref_noeud
# debut de la liste des ddl a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_noeud>
deb_list_ddl_noeud fin_list_ddl_noeud
type_sortie_ddl_retenue= 1
# debut de la liste des ddl etendus a considerer aux noeuds, (une liste de ddl), puis <fin_list_ddl_etendu_noeud>
deb_list_ddl_etendu_noeud fin_list_ddl_etendu_noeud # fin des ddl etendu aux noeuds
# debut de la liste des grandeurs particulieres a considerer aux noeuds, (une liste de string), puis <fin_list_GrandParticuliere_noeud>
deb_list_GrandParticuliere_noeud fin_list_GrandParticuliere_noeud # fin des grandeurs particulieres aux noeuds
fin_liste_ddl_et_noeud # fin des grandeurs aux noeuds
debut_liste_ddl_ptinteg # ** debut des grandeurs aux elements
# debut de la liste des elements et points d'integration, une liste de (un element, un numero de pt d'integ), puis <fin_list_NbElement_NbPtInteg>
deb_list_NbElement_NbPtInteg fin_list_NbElement_NbPtInteg
# debut de la liste des ref d'elements, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_element>
deb_list_ref_element fin_list_ref_element
# debut de la liste des ref de ptinteg d'elements, puis une liste de: nom <chaine de caracteres > puis <fin_list_ref_ptinteg_element>
deb_list_ref_ptinteg_element fin_list_ref_ptinteg_element
# debut de la liste des ddl a considerer pour les elements, (une liste de ddl), puis <fin_list_ddl_element>
deb_list_ddl_element fin_list_ddl_element # fin de la liste de ddl a considerer pour les elements
# debut de la liste des grandeurs particulieres a considerer pour les elements, (une liste de string), puis <fin_list_GrandParticuliere_element>
deb_list_GrandParticuliere_element fin_list_GrandParticuliere_element # fin des grandeurs particulieres aux elements
# tableau de grandeurs evoluees aux elements a visualiser, un par maillage
deb_list_GrandEvoluee_element fin_list_GrandEvoluee_element
fin_liste_ddl_ptinteg # fin des grandeurs aux elements
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 <fin_list_NbFaceElement_NbPtInteg>
deb_list_NbFaceElement_NbPtInteg fin_list_NbFaceElement_NbPtInteg
# debut de la liste des ref de faces d'element, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_F_element>
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 <chaine de caracteres > puis <fin_list_ref_ptinteg_F_element>
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 <fin_list_Grandeur_F_element>
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 <fin_list_NbAreteElement_NbPtInteg>
deb_list_NbAreteElement_NbPtInteg fin_list_NbAreteElement_NbPtInteg
# debut de la liste des ref d'arete d'element, puis une liste de: nom <chaine de caracteres > + numero d'integ, puis <fin_list_ref_A_element>
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 <chaine de caracteres > puis <fin_list_ref_ptinteg_A_element>
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 <fin_list_Grandeur_A_element>
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 <un reel> (indique le temps en seconde du cycle de l'animation)
fin_animation # un mot cle de fin
fin_visualisation_maple
# =================================================================================
# || fin de la visualisation maple ||
# =================================================================================
fin_fichier_commande_visu # <<<<<< le mot cle <fin_fichier_commande_visu> permet
# l'arret de la lecture des commandes, apres ce mot cle, aucune commande n'est lu, de plus
# sans le mot cle de fin de fichier, le fichier n'est pas valide
###############################################################################################

177
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact.info Normal file
View file

@ -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_

156
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/biel_contact_princ.maple.ref1 Normal file
View file

@ -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

77
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec100x100.her Normal file
View file

@ -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

84
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_100x100.her Normal file
View file

@ -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

122
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100.her Normal file
View file

@ -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

81
Batterie/FONCTION_UNIQUE/CONTACT/Test_R_contact_100biel_20triangles/rec_tri_50x100_1.her Normal file
View file

@ -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