|
1
|
#!/usr/bin/perl
|
|
2
|
#!/usr/local/bin/perl
|
|
3
|
use strict;
|
|
4
|
use English;
|
|
5
|
use File::Basename;
|
|
6
|
my $NOM_PROG = basename $PROGRAM_NAME;
|
|
7
|
|
|
8
|
|
|
9
|
my $VERSION = '1.01';
|
|
10
|
#####################################################################################################
|
|
11
|
# script pour visualiser un ou plusieurs maillages dans Gmsh #
|
|
12
|
# version 1.00 : version initiale #
|
|
13
|
# version 1.01 : which du package File::Which n est plus utilise pour verifier l existence #
|
|
14
|
# d une commande (pour eviter une erreur dans le cas ou ce package n est pas #
|
|
15
|
# installe). A la place, on regarde la variable $PATH via la subroutine verif_cmd() #
|
|
16
|
#####################################################################################################
|
|
17
|
|
|
18
|
|
|
19
|
|
|
20
|
|
|
21
|
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
|
|
22
|
# VARIABLES QUE L UTILISATEUR PEUT ETRE AMENE A RENSEIGNER #
|
|
23
|
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
|
|
24
|
#
|
|
25
|
# les variables $exeHZ et $exeGMSH servent a imposer le choix des executables Herezh et Gmsh
|
|
26
|
# (rq : les alias shell ne fonctionneront pas)
|
|
27
|
#
|
|
28
|
#commande HZ++ definie par l utilisateur (modifiable egalement avec l option -exeHZ)
|
|
29
|
my $exeHZ; # = 'mon_executable_a_moi';
|
|
30
|
#commande gmsh definie par l utilisateur (modifiable egalement avec l option -exeGMSH)
|
|
31
|
my $exeGMSH; # = 'mon_executable_a_moi';
|
|
32
|
|
|
33
|
|
|
34
|
#commande Herezh par defaut selon la plateforme
|
|
35
|
# rq : ces variables seront utilisees uniquement si la variable $exeHZ n a pas ete definie
|
|
36
|
my $EXE_HZ_Linux64bits = 'HZppfast64';#linux 64 bits
|
|
37
|
my $EXE_HZ_MacOSX = 'HZppfast_Vn-1';#MacOSX
|
|
38
|
my $EXE_HZ_autres = 'HZppfast';#tous les autres systemes
|
|
39
|
#commande gmsh par defaut (quelque soit la plateforme)
|
|
40
|
# rq : cette variable sera utilisee uniquement si la variable $exeGMSH n a pas ete definie
|
|
41
|
my $exeGMSH_defaut = 'gmsh';
|
|
42
|
|
|
43
|
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
|
|
44
|
# FIN VARIABLES UTILISATEUR #
|
|
45
|
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
|
|
46
|
|
|
47
|
|
|
48
|
|
|
49
|
|
|
50
|
sub affichage_aide {
|
|
51
|
use Text::Wrap;
|
|
52
|
#config du package Text::Wrap
|
|
53
|
$Text::Wrap::columns = 80;#le nombre de caracteres maximum par ligne sera egal a ($Text::Wrap::columns - 1) dans le cas ou utilise wrap
|
|
54
|
|
|
55
|
#indentation pour l affichage de Usage 2
|
|
56
|
my $indent_Usage2 = ""; $indent_Usage2 .= " " for(1 .. length(" > $NOM_PROG "));
|
|
57
|
|
|
58
|
print "--------------------------------\n";
|
|
59
|
print " script $NOM_PROG\n";
|
|
60
|
print "--------------------------------\n";
|
|
61
|
print "\n";
|
|
62
|
print "But : visualiser des maillages Herezh avec Gmsh\n";
|
|
63
|
print "\n";
|
|
64
|
print "Usage 1 : mode interactif\n";
|
|
65
|
print " > $NOM_PROG\n";
|
|
66
|
print "Usage 2 : lancement avec arguments\n";
|
|
67
|
print wrap(" ",$indent_Usage2, "> $NOM_PROG [-h|help] [-v] [-exeHZ path_HZ] [-exeGMSH path_GMSH] [-saveVisu] fher_1 [fher_2 .. fher_N]\n");
|
|
68
|
print "\n";
|
|
69
|
print "Arguments :\n";
|
|
70
|
print " - fic_1 : 1er maillage .her\n";
|
|
71
|
print " - [fher_2 .. fher_N] : maillages .her supplementaires eventuels\n";
|
|
72
|
print "\n";
|
|
73
|
print "Fonctionnement :\n";
|
|
74
|
print " $NOM_PROG lance un calcul herezh pour creer un fichier .msh\n";
|
|
75
|
print " qui est ensuite visualise dans Gmsh. La visualisation des references est\n";
|
|
76
|
print " desactivee a l ouverture (exceptee la derniere vue qui est specialement\n";
|
|
77
|
print " creee par ce script pour afficher les elements en vue solide)\n";
|
|
78
|
print " Le code couleur pour les references est :\n";
|
|
79
|
print " - reference de noeuds => rouge\n";
|
|
80
|
print " - reference d aretes => jaune\n";
|
|
81
|
print " - reference de faces => vert\n";
|
|
82
|
print " - reference d elements => bleu\n";
|
|
83
|
print " - reference de points d integration => mauve\n";
|
|
84
|
print " Dans le cas de maillages avec beaucoup d elements, la preparation des fichiers\n";
|
|
85
|
print " peut prendre du temps. Il est alors conseille d utiliser l option -saveVisu pour\n";
|
|
86
|
print " conserver les fichiers de visualisation pour pouvoir les reouvrir ulterieurment\n";
|
|
87
|
print " sans avoir a reexecuter ce script.\n";
|
|
88
|
print "\n";
|
|
89
|
print "Options :\n";
|
|
90
|
print " -v : affichage du numero de version\n";
|
|
91
|
print "\n";
|
|
92
|
print " -exeHZ path_HZ : choisir l executable Herezh++. Le calcul Herezh se\n";
|
|
93
|
print " fera avec l executable path_HZ\n";
|
|
94
|
print " par defaut : linux 64 bits => path_HZ=HZppfast64\n";
|
|
95
|
print " MacOSX (darwin) => path_HZ=HZppfast_Vn-1\n";
|
|
96
|
print " autres => path_HZ=HZppfast\n";
|
|
97
|
print " (a noter que les alias shell ne fonctionneront pas)\n";
|
|
98
|
print " (a noter que cette option fonctionne aussi en mode interactif)\n";
|
|
99
|
print "\n";
|
|
100
|
print " -exeGMSH path_GMSH : choisir l executable Gmsh. La visualisation Gmsh se\n";
|
|
101
|
print " fera avec l executable path_GMSH\n";
|
|
102
|
print " par defaut : path_GMSH=gmsh (quelque soit la plateforme)\n";
|
|
103
|
print " (a noter que les alias shell ne fonctionneront pas)\n";
|
|
104
|
print " (a noter que cette option fonctionne aussi en mode interactif)\n";
|
|
105
|
print "\n";
|
|
106
|
print " -saveVisu : sauvegarde des fichiers .geo et .msh de visu\n";
|
|
107
|
print "\n";
|
|
108
|
print "Exemples :\n";
|
|
109
|
print " o visualiser les maillages mail1.her et mail2.her :\n";
|
|
110
|
print " > $NOM_PROG mail1.her mail2.her\n";
|
|
111
|
print " o visualiser le maillage mail1.her en choisissant l executable Herezh :\n";
|
|
112
|
print " > $NOM_PROG -exeHZ HZpp_perso mail1.her\n";
|
|
113
|
print " o visualiser le maillage mail1.her et sauvegarder les fichiers de visu :\n";
|
|
114
|
print " > $NOM_PROG -saveVisu mail1.her\n";
|
|
115
|
print " o visualiser le maillage mail1.her en choisissant l executable Herezh et\n";
|
|
116
|
print " et l executable Gmsh (exemple dans le cas ou les executables se trouvent\n";
|
|
117
|
print " quelque part dans le HOME)\n";
|
|
118
|
print " > $NOM_PROG -exeHZ ~/mon_rep_HZ/HZpp -exeGMSH ~/mon_rep_GMSH/gmsh mail1.her\n";
|
|
119
|
print "\n";
|
|
120
|
print "Auteur :\n";
|
|
121
|
print " TROUFFLARD Julien\n";
|
|
122
|
print " julien.troufflard\@univ-ubs.fr\n";
|
|
123
|
print " julien.troufflard\@free.fr\n";
|
|
124
|
print "--------------------------------\n";
|
|
125
|
}
|
|
126
|
|
|
127
|
|
|
128
|
|
|
129
|
#------------------------------------
|
|
130
|
#option -h ou -help => affichage de l aide et arret
|
|
131
|
#------------------------------------
|
|
132
|
# rq : insensible a la casse
|
|
133
|
foreach my $arg (@ARGV) {
|
|
134
|
if(($arg =~ /^-h$/i) or ($arg =~ /^-help$/i)) {
|
|
135
|
affichage_aide();
|
|
136
|
exit;
|
|
137
|
}
|
|
138
|
}
|
|
139
|
|
|
140
|
#------------------------------------
|
|
141
|
#option -v => affichage de la version et arret
|
|
142
|
#------------------------------------
|
|
143
|
foreach my $arg (@ARGV) {
|
|
144
|
if($arg eq '-v') {
|
|
145
|
print "\n $NOM_PROG : version $VERSION\n\n";
|
|
146
|
exit;
|
|
147
|
}
|
|
148
|
}
|
|
149
|
|
|
150
|
|
|
151
|
|
|
152
|
#------------------------------------
|
|
153
|
#recuperation des arguments et options
|
|
154
|
#------------------------------------
|
|
155
|
my @liste_fher;#liste des maillages
|
|
156
|
my $is_opt_saveVisu = 0;#indicateur de l option -saveVisu (sauvegarde des fichiers .geo et _Gmsh.msh de visualisation)
|
|
157
|
|
|
158
|
#on parcourt la liste des arguments (on traite les options connues et on stocke les autres dans @args)
|
|
159
|
my $opt;
|
|
160
|
my @args;
|
|
161
|
while($#ARGV != -1) {
|
|
162
|
$opt = shift(@ARGV);
|
|
163
|
|
|
164
|
#option -exeHZ
|
|
165
|
if($opt eq '-exeHZ') {
|
|
166
|
$exeHZ = shift(@ARGV);
|
|
167
|
}
|
|
168
|
#option -exeGMSH
|
|
169
|
elsif($opt eq '-exeGMSH') {
|
|
170
|
$exeGMSH = shift(@ARGV);
|
|
171
|
}
|
|
172
|
#option -saveVisu
|
|
173
|
elsif($opt eq '-saveVisu') {
|
|
174
|
$is_opt_saveVisu = 1;
|
|
175
|
}
|
|
176
|
|
|
177
|
#cas d une option inconnue
|
|
178
|
elsif($opt =~ /^-/) {
|
|
179
|
warn "**Attention : option $opt inconnue (on ignore cette option)...\n";
|
|
180
|
}
|
|
181
|
|
|
182
|
#autres arguments
|
|
183
|
else {
|
|
184
|
push(@args,$opt);
|
|
185
|
}
|
|
186
|
}#while($#ARGV != -1)
|
|
187
|
|
|
188
|
|
|
189
|
#---------------------
|
|
190
|
#verif de l executable Herezh
|
|
191
|
#---------------------
|
|
192
|
#si la variable $exeHZ n a pas ete renseigne au prealable => on selectionne l executable par defaut en fonction de la plateforme
|
|
193
|
if(not defined($exeHZ)) {
|
|
194
|
#- - - - - -
|
|
195
|
#type de plateforme
|
|
196
|
#- - - - - -
|
|
197
|
my $type_OS;
|
|
198
|
#-avec uname
|
|
199
|
if(verif_cmd('uname')) {
|
|
200
|
$type_OS = qx(uname -a);
|
|
201
|
chomp($type_OS);
|
|
202
|
}
|
|
203
|
#-sinon : warning (OS inconnue)
|
|
204
|
else {warn "**Attention : impossible de saisir le type de systeme d exploitation avec uname -a ...\n";}
|
|
205
|
|
|
206
|
#selection de l executable Herezh
|
|
207
|
#-linux 64 bits
|
|
208
|
if($type_OS =~ /linux/i and $type_OS =~ /x86_64/i) {$exeHZ = $EXE_HZ_Linux64bits;}
|
|
209
|
#-MacOSX (darwin)
|
|
210
|
elsif($type_OS =~ /darwin/i) {$exeHZ = $EXE_HZ_MacOSX;}
|
|
211
|
#-tous les autres (y compris linux 32 bits)
|
|
212
|
else{$exeHZ = $EXE_HZ_autres;}
|
|
213
|
}
|
|
214
|
#verif de l executable Herezh
|
|
215
|
verif_cmd($exeHZ) or die "\n**Erreur (prog:$NOM_PROG) : executable Herezh ($exeHZ) introuvable (pour eviter cette erreur : vous pouvez soit utiliser l option -exeHZ, soit renseigner directement la commande Herezh dans la variable \$exeHZ au debut du script)...\n\n";
|
|
216
|
|
|
217
|
#---------------------
|
|
218
|
#verif de l executable Gmsh
|
|
219
|
#---------------------
|
|
220
|
#si la variable $exeGMSH n a pas ete renseigne au prealable => on selectionne l executable par defaut
|
|
221
|
$exeGMSH = $exeGMSH_defaut if(not defined($exeGMSH));
|
|
222
|
#verif de l executable Gmsh
|
|
223
|
verif_cmd($exeGMSH) or die "\n**Erreur (prog:$NOM_PROG) : executable Gmsh ($exeGMSH) introuvable (pour eviter cette erreur : vous pouvez soit utiliser l option -exeGMSH, soit renseigner directement la commande Gmsh dans la variable \$exeGMSH au debut du script)...\n\n";
|
|
224
|
|
|
225
|
|
|
226
|
|
|
227
|
|
|
228
|
|
|
229
|
#---------------------
|
|
230
|
#liste des maillages (on en profite pour faire des verifs sur la validite du fichier)
|
|
231
|
#---------------------
|
|
232
|
#
|
|
233
|
# la liste @args va servir de liste temporaire pour la saisie des noms de fichier
|
|
234
|
#
|
|
235
|
#cas avec arguments
|
|
236
|
if($#args >= 0) {
|
|
237
|
#rien a faire pour l instant (la liste @args a deja ete remplie par les arguments)
|
|
238
|
}
|
|
239
|
#cas en mode interactif (on remplit la liste @args avec un menu interactif)
|
|
240
|
else {
|
|
241
|
print "\nChoix des maillages a visualiser :\n";
|
|
242
|
print " (taper liste pour afficher la liste actuelle)\n";
|
|
243
|
my $choix = -1;
|
|
244
|
while() {
|
|
245
|
print " Donner un nom de maillage (f pour finir) : ";
|
|
246
|
$choix = <STDIN>; chomp($choix);
|
|
247
|
next if($choix eq '');#cas ou l utilisateur a tape "Entree"
|
|
248
|
|
|
249
|
#choix f => on arrete la saisie
|
|
250
|
last if($choix eq 'f');
|
|
251
|
|
|
252
|
if($choix eq 'liste') {
|
|
253
|
print " Liste actuelle :\n @args\n\n" if($choix eq 'liste');
|
|
254
|
next;
|
|
255
|
}
|
|
256
|
|
|
257
|
push(@args, $choix);
|
|
258
|
}
|
|
259
|
print "\n";
|
|
260
|
}
|
|
261
|
#-verif de la validite des maillages
|
|
262
|
foreach my $arg (@args) {
|
|
263
|
$arg .= '.her' if(not $arg =~ /\.her$/);#rajout de l extension .her si manquante
|
|
264
|
#existence du fichier
|
|
265
|
(-e $arg) or do {warn "**Attention : fichier $arg introuvable (on ignore ce fichier)...\n"; next;};
|
|
266
|
#ouverture du fichier
|
|
267
|
open(FIC, "<$arg") or do {warn "**Attention : impossible d ouvrir le fichier $arg (on ignore ce fichier)...\n"; next;};
|
|
268
|
#est-ce un maillage Herezh ? (verif par la presence des mots noeuds et NOEUDS)
|
|
269
|
my ($is_noeuds, $is_NOEUDS) = (0,0);
|
|
270
|
while(<FIC>) {
|
|
271
|
$is_noeuds = 1 if(/^\s*noeuds/);
|
|
272
|
$is_NOEUDS = 1 if(/^\s*\d+\s+NOEUDS/);
|
|
273
|
last if($is_NOEUDS);
|
|
274
|
}
|
|
275
|
close(FIC);
|
|
276
|
$is_noeuds or do {warn "**Attention : le fichier $arg ne contient pas le mot \"noeuds\" (on ignore ce fichier)...\n"; next;};
|
|
277
|
$is_NOEUDS or do {warn "**Attention : le fichier $arg ne contient pas le mot \"NOEUDS\" (on ignore ce fichier)...\n"; next;};
|
|
278
|
|
|
279
|
#fichier ok => ajout a la liste
|
|
280
|
push(@liste_fher, $arg);
|
|
281
|
}
|
|
282
|
|
|
283
|
#si a ce stade, la liste des maillages est vide => on arrete le programme
|
|
284
|
($#liste_fher > -1) or die "\nArret du programme car aucun maillage valide n a ete fourni...\n\n";
|
|
285
|
|
|
286
|
|
|
287
|
|
|
288
|
#---------------------
|
|
289
|
#verif des maillages (constitution de la liste des noms de maillage)
|
|
290
|
#---------------------
|
|
291
|
print " verification des maillages...\n";
|
|
292
|
my @liste_nom_maillage;
|
|
293
|
#-verif 1 : si il y a un seul maillage et qu il n a pas de nom, on indique "premier_maillage" dans la liste des noms de maillages
|
|
294
|
if($#liste_fher == 0) {
|
|
295
|
my $nom_maillage = '';
|
|
296
|
open(FIC, "<$liste_fher[0]") or die "\n**Erreur (prog:$NOM_PROG) : impossible d ouvrir le fichier $liste_fher[0] ...\n\n";
|
|
297
|
while(<FIC>) {
|
|
298
|
last if(/^\s*noeuds/);
|
|
299
|
next if(not /^\s*nom_maillage\s+(\S+)/);
|
|
300
|
$nom_maillage = $1;
|
|
301
|
last;
|
|
302
|
}
|
|
303
|
close(FIC);
|
|
304
|
$nom_maillage = 'premier_maillage' if($nom_maillage eq '');
|
|
305
|
push(@liste_nom_maillage, $nom_maillage);
|
|
306
|
}
|
|
307
|
|
|
308
|
#-verif 2 : si il y a plusieurs maillages, il faut qu il ait chacun un nom_maillage different (sinon arret du programme)
|
|
309
|
if($#liste_fher > 0) {
|
|
310
|
my $is_maillage_ok = 1;
|
|
311
|
foreach my $fher (@liste_fher) {
|
|
312
|
my $nom_maillage = '';
|
|
313
|
open(FIC, "<$fher") or die "\n**Erreur (prog:$NOM_PROG) : impossible d ouvrir le fichier $fher ...\n\n";;
|
|
314
|
while(<FIC>) {
|
|
315
|
next if(not /^\s*nom_maillage\s+(\S+)/);
|
|
316
|
$nom_maillage = $1;
|
|
317
|
last;
|
|
318
|
}
|
|
319
|
close(FIC);
|
|
320
|
push(@liste_nom_maillage, $nom_maillage);
|
|
321
|
}
|
|
322
|
|
|
323
|
my @maillages_deja_traites; for(my $i=0; $i<=$#liste_nom_maillage; $i++) {$maillages_deja_traites[$i] = 0;}
|
|
324
|
for(my $i=0; $i<=$#liste_nom_maillage; $i++) {
|
|
325
|
next if($maillages_deja_traites[$i]);
|
|
326
|
|
|
327
|
#cas d un maillage sans nom
|
|
328
|
if($liste_nom_maillage[$i] eq '') {
|
|
329
|
warn "**Erreur (prog:$NOM_PROG) : le maillage $liste_fher[$i] n a pas de nom (nom_maillage non specifie)...\n";
|
|
330
|
$is_maillage_ok = 0;
|
|
331
|
next;
|
|
332
|
}
|
|
333
|
|
|
334
|
#cas d un maillage ayant le meme nom qu un ou plusieurs autres maillages
|
|
335
|
my @liste_maillages_meme_nom = ();
|
|
336
|
for(my $j=$i+1; $j<=$#liste_nom_maillage; $j++) {
|
|
337
|
if($liste_nom_maillage[$i] eq $liste_nom_maillage[$j]) {
|
|
338
|
$maillages_deja_traites[$j] = 1;
|
|
339
|
push(@liste_maillages_meme_nom, $liste_fher[$j]);
|
|
340
|
}
|
|
341
|
}
|
|
342
|
if($#liste_maillages_meme_nom > -1) {
|
|
343
|
warn "**Erreur (prog:$NOM_PROG) : les maillages suivants ont le meme nom => $liste_fher[$i] @liste_maillages_meme_nom\n";
|
|
344
|
$is_maillage_ok = 0;
|
|
345
|
}
|
|
346
|
}
|
|
347
|
|
|
348
|
#arret du programme si on a trouve des maillages sans nom ou des noms utilises plusieurs fois
|
|
349
|
$is_maillage_ok or die "\nArret du programme a cause d un probleme sur les noms de maillages...\n\n";
|
|
350
|
}
|
|
351
|
|
|
352
|
|
|
353
|
#---------------------
|
|
354
|
#listes des elements 1D, 2D et 3D (pour affecter des lois de type LOI_RIEN)
|
|
355
|
#---------------------
|
|
356
|
print " preparation du calcul Herezh...\n";
|
|
357
|
#-prefixes et suffixe pour la reconnaissance des elements 1D, 2D et axisymetriques
|
|
358
|
my @PREFIXE_1D = qw(POUT);
|
|
359
|
my @PREFIXE_2D = qw(TRIA QUAD);
|
|
360
|
my @SUFFIXE_AXI = qw(_AXI);
|
|
361
|
|
|
362
|
#-table de hachage pour definir pour chaque maillage, la liste des elements 1D, 2D, 3D (rq : les elements AXI sont consideres comme 3D car il necessite une loi 3D)
|
|
363
|
my %TAB_DIM_LOI;#@{$TAB_DIM_LOI{indice maillage}{dim loi}} = (liste elements) (par exemple, pour le premier maillage => @{$TAB_DIM_LOI{0}{'3D'}} = (1,5,10)
|
|
364
|
|
|
365
|
#on boucle sur les maillages pour constituer les listes d elements par dimension de loi de comportement pour chaque maillage
|
|
366
|
# rq : on en profite pour reperer si il y a au moins 1 element 1D et au moins un element 2D (pour savoir si il faudra renseigner le mot-cle sections et epaisseurs)
|
|
367
|
my $is_elt_1D = 0;#indicateur de la presence d au moins 1 element 1D
|
|
368
|
my $is_elt_2D = 0;#indicateur de la presence d au moins 1 element 2D
|
|
369
|
my $nb_elts_tot = 0;
|
|
370
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
371
|
#saisie des elements
|
|
372
|
my ($nb_elts, $ref_elements);
|
|
373
|
($_, $_, $_, $nb_elts, $ref_elements) = lecture_mail_her($liste_fher[$no_mail]);
|
|
374
|
$nb_elts_tot += $nb_elts;
|
|
375
|
my @ELEM_1D = ();
|
|
376
|
my @ELEM_2D = ();
|
|
377
|
my @ELEM_3D = ();
|
|
378
|
ELEM:for(my $i=1; $i<=$nb_elts; $i++) {
|
|
379
|
#verif si element AXI (=> loi LOI_RIEN3D )
|
|
380
|
foreach my $suffixe (@SUFFIXE_AXI) {
|
|
381
|
@_ = split(/\s+/, $ref_elements->{$i}{'TYPE'});
|
|
382
|
if($_[0] =~ /$suffixe\s*$/) {
|
|
383
|
push(@ELEM_3D, $i);
|
|
384
|
next ELEM;
|
|
385
|
}
|
|
386
|
}
|
|
387
|
#verif si element 1D (=> loi LOI_RIEN1D )
|
|
388
|
foreach my $prefixe (@PREFIXE_1D) {
|
|
389
|
if($ref_elements->{$i}{'TYPE'} =~ /^\s*$prefixe/) {
|
|
390
|
push(@ELEM_1D, $i);
|
|
391
|
$is_elt_1D = 1;
|
|
392
|
next ELEM;
|
|
393
|
}
|
|
394
|
}
|
|
395
|
#verif si element 2D (=> loi LOI_RIEN2D_C )
|
|
396
|
foreach my $prefixe (@PREFIXE_2D) {
|
|
397
|
if($ref_elements->{$i}{'TYPE'} =~ /^\s*$prefixe/) {
|
|
398
|
push(@ELEM_2D, $i);
|
|
399
|
$is_elt_2D = 1;
|
|
400
|
next ELEM;
|
|
401
|
}
|
|
402
|
}
|
|
403
|
#sinon, c est un element 3D (=> loi LOI_RIEN3D )
|
|
404
|
push(@ELEM_3D, $i);
|
|
405
|
}#FIN BOUCLE SUR LES ELEMENTS DU MAILLAGE indice $no_mail
|
|
406
|
|
|
407
|
#remplissage de la table pour ce maillage
|
|
408
|
push(@{$TAB_DIM_LOI{$no_mail}{'1D'}}, @ELEM_1D);
|
|
409
|
push(@{$TAB_DIM_LOI{$no_mail}{'2D'}}, @ELEM_2D);
|
|
410
|
push(@{$TAB_DIM_LOI{$no_mail}{'3D'}}, @ELEM_3D);
|
|
411
|
|
|
412
|
}#FIN BOUCLE SUR LES MAILLAGES
|
|
413
|
|
|
414
|
|
|
415
|
|
|
416
|
#--------------------
|
|
417
|
#nom des fichiers .info, .CVisu, _Gmsh.msh et .geo temporaires (on s assure qu ils n existent pas deja)
|
|
418
|
#--------------------
|
|
419
|
my $racine_fic_tmp = $NOM_PROG; $racine_fic_tmp =~ s/\.\S+$//;
|
|
420
|
my $no = 0;
|
|
421
|
$racine_fic_tmp .= "_$no";
|
|
422
|
my $finfo = "$racine_fic_tmp.info";#fichier de calcul temporaire
|
|
423
|
my $fCVisu = "$racine_fic_tmp.CVisu";#.CVisu associe
|
|
424
|
my $fGmsh = "$racine_fic_tmp\_Gmsh.msh";#.msh qui va etre cree apres calcul
|
|
425
|
my $fgeo = "$racine_fic_tmp.geo";#.geo qui sera utilise pour lancer la visu
|
|
426
|
|
|
427
|
while(-e $finfo or -e $fCVisu or -e $fGmsh or -e $fgeo) {
|
|
428
|
$no++;
|
|
429
|
$racine_fic_tmp = $NOM_PROG; $racine_fic_tmp =~ s/\.\S+$//;
|
|
430
|
$racine_fic_tmp .= "_$no";
|
|
431
|
$finfo = "$racine_fic_tmp.info";
|
|
432
|
$fCVisu = "$racine_fic_tmp.CVisu";
|
|
433
|
$fGmsh = "$racine_fic_tmp\_Gmsh.msh";
|
|
434
|
$fgeo = "$racine_fic_tmp.geo";
|
|
435
|
}
|
|
436
|
#-memorisation des eventuels fichiers deja existants qui commencent comme le fichier .info (pour ne pas les effacer a la fin du script)
|
|
437
|
my @liste_fic_a_ne_pas_effacer = glob("$racine_fic_tmp*");
|
|
438
|
|
|
439
|
|
|
440
|
#---------------------
|
|
441
|
#ecriture du .info
|
|
442
|
#---------------------
|
|
443
|
open(FIC, ">$finfo");
|
|
444
|
print FIC "dimension 3\n\n";
|
|
445
|
|
|
446
|
print FIC "niveau_commentaire 1\n\n";
|
|
447
|
|
|
448
|
#on choisit de lancer le calcul en dynamique_explicite en prevision des maillages a grand nombre de noeuds (plus rapide que non_dynamique pour faire un increment)
|
|
449
|
print FIC "TYPE_DE_CALCUL\ndynamique_explicite avec plus lectureCommandesVisu\n";
|
|
450
|
|
|
451
|
#ecriture des maillages et de references d elements speciales
|
|
452
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
453
|
my $fher = $liste_fher[$no_mail];
|
|
454
|
my $flis = $fher; $flis =~ s/.her$/.lis/;
|
|
455
|
print FIC "\n< $fher\n";
|
|
456
|
print FIC "< $flis\n" if(-e $flis);#inclusion du .lis si il existe
|
|
457
|
#set d elements speciaux en fonction de la dimension de la loi de comportement
|
|
458
|
# rq : on utilise une subroutine qui ecrit les references avec 15 elements max par ligne (sinon bug Herezh si il y a trop d elements par ligne)
|
|
459
|
foreach my $dim ('1D', '2D', '3D') {
|
|
460
|
next if($#{$TAB_DIM_LOI{$no_mail}{$dim}} == -1);#pas d ecriture si aucun element de dimension $dim
|
|
461
|
ecrire_liste_N_E(*FIC, "E_tmp_visu_elem_$dim", @{$TAB_DIM_LOI{$no_mail}{$dim}});
|
|
462
|
}
|
|
463
|
}
|
|
464
|
|
|
465
|
print FIC "\nchoix_materiaux\n";
|
|
466
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
467
|
#choix materiau par dimension de loi
|
|
468
|
foreach my $dim ('1D', '2D', '3D') {
|
|
469
|
next if($#{$TAB_DIM_LOI{$no_mail}{$dim}} == -1);#pas d ecriture si aucun element de dimension $dim
|
|
470
|
print FIC "nom_mail= $liste_nom_maillage[$no_mail] E_tmp_visu_elem_$dim MAT$dim\n";
|
|
471
|
}
|
|
472
|
}
|
|
473
|
|
|
474
|
print FIC "\nmateriaux\n";
|
|
475
|
print FIC "MAT1D LOI_RIEN1D\nMAT2D LOI_RIEN2D_C\nMAT3D LOI_RIEN3D\n\n";
|
|
476
|
|
|
477
|
print FIC "masse_volumique\n";
|
|
478
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
479
|
foreach my $dim ('1D', '2D', '3D') {
|
|
480
|
next if($#{$TAB_DIM_LOI{$no_mail}{$dim}} == -1);#pas d ecriture si aucun element de dimension $dim
|
|
481
|
print FIC "nom_mail= $liste_nom_maillage[$no_mail] E_tmp_visu_elem_$dim 1.\n";
|
|
482
|
}
|
|
483
|
}
|
|
484
|
|
|
485
|
#sections pour les eventuels elements 1D
|
|
486
|
if($is_elt_1D) {
|
|
487
|
print FIC "\nsections\n";
|
|
488
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
489
|
next if($#{$TAB_DIM_LOI{$no_mail}{'1D'}} == -1);#pas d ecriture si aucun element de dimension $dim
|
|
490
|
print FIC "nom_mail= $liste_nom_maillage[$no_mail] E_tmp_visu_elem_1D 1.\n";
|
|
491
|
}
|
|
492
|
}
|
|
493
|
|
|
494
|
#epaisseurs pour les eventuels elements 2D
|
|
495
|
if($is_elt_2D) {
|
|
496
|
print FIC "\nepaisseurs\n";
|
|
497
|
for(my $no_mail=0; $no_mail<=$#liste_fher; $no_mail++) {
|
|
498
|
next if($#{$TAB_DIM_LOI{$no_mail}{'2D'}} == -1);#pas d ecriture si aucun element de dimension $dim
|
|
499
|
print FIC "nom_mail= $liste_nom_maillage[$no_mail] E_tmp_visu_elem_2D 1.\n";
|
|
500
|
}
|
|
501
|
}
|
|
502
|
|
|
503
|
print FIC "\ncharges\n\n";
|
|
504
|
print FIC "blocages\n\n";
|
|
505
|
|
|
506
|
#controle => un seul increment
|
|
507
|
print FIC "controle\n";
|
|
508
|
print FIC "DELTAt 1.\n";
|
|
509
|
print FIC "TEMPSFIN 1.\n";
|
|
510
|
print FIC "SAUVEGARDE 0\n";
|
|
511
|
print FIC "MAXINCRE 1\n";
|
|
512
|
|
|
513
|
print FIC "\npara_pilotage_equi_global\n\n";
|
|
514
|
|
|
515
|
print FIC "para_syteme_lineaire\n\n";
|
|
516
|
|
|
517
|
print FIC "para_affichage\nFREQUENCE_SORTIE_FIL_DU_CALCUL 1\n\n";
|
|
518
|
|
|
519
|
print FIC "resultats pas_de_sortie_finale_\nCOPIE 0\n\n";
|
|
520
|
|
|
521
|
print FIC "_fin_point_info_\n";
|
|
522
|
close(FIC);
|
|
523
|
|
|
524
|
|
|
525
|
#---------------------
|
|
526
|
#ecriture du .CVisu
|
|
527
|
#---------------------
|
|
528
|
open(FIC, ">$fCVisu");
|
|
529
|
print FIC "
|
|
530
|
debut_fichier_commande_visu
|
|
531
|
|
|
532
|
debut_visualisation_Gmsh
|
|
533
|
debut_maillage_initial
|
|
534
|
actif 1
|
|
535
|
pseudo-homothetie_sur_les_maillages_ 0
|
|
536
|
visualisation_references_sur_les_maillages_ 1
|
|
537
|
fin_maillage_initial
|
|
538
|
|
|
539
|
debut_choix_maillage
|
|
540
|
actif 0
|
|
541
|
1";
|
|
542
|
for(my $i=1; $i<=$#liste_fher; $i++) {$_ = $i + 1; print FIC " $_";}
|
|
543
|
print FIC " fin_choix_maillage
|
|
544
|
fin_visualisation_Gmsh
|
|
545
|
|
|
546
|
fin_fichier_commande_visu
|
|
547
|
";
|
|
548
|
close(FIC);
|
|
549
|
|
|
550
|
|
|
551
|
|
|
552
|
#---------------------
|
|
553
|
#lancement du calcul
|
|
554
|
#---------------------
|
|
555
|
#-on lance le calcul avec redirection dans un fichier .log (au cas ou le calcul a plante => on propose la lecture de ce .log a l utilisateur)
|
|
556
|
my $fredir = "$racine_fic_tmp.log";
|
|
557
|
print " creation du fichier _Gmsh.msh (calcul Herezh en cours)...\n";
|
|
558
|
system("rm -f $fGmsh $fredir");
|
|
559
|
system("echo \'#\' > $fredir; echo \'#fichier genere suite a l execution du script $NOM_PROG\' >> $fredir; echo \'#\' >> $fredir; echo \'\' >> $fredir");#affichage d un en-tete dans le fichier .log pour indiquer que ce fichier a ete cree lors de l execution de ce script
|
|
560
|
system("echo n | $exeHZ -f $finfo >> $fredir");#on lance avec "echo n" pour repondre automatiquement "non" au cas ou il y a une erreur de calcul Herezh
|
|
561
|
|
|
562
|
#cas ou la calcul n a pas fonctionne (si le fichier _Gmsh.msh n a pas ete cree)
|
|
563
|
# => on propose a l utilisateur de visualiser le .log pour savoir ce qui s est passe
|
|
564
|
if(not -e $fGmsh) {
|
|
565
|
print "\n**Erreur (prog:$NOM_PROG) : le calcul Herezh++ n a pas fonctionne (le fichier pour Gmsh n a pas ete cree)...\n\n";
|
|
566
|
my $choix = -1;
|
|
567
|
while($choix ne 'o' and $choix ne 'n') {
|
|
568
|
print " Voulez-vous visualiser l affichage Herezh++ du calcul ? (o/n) ";
|
|
569
|
$choix = <STDIN>; chomp($choix); $choix = lc($choix);
|
|
570
|
}
|
|
571
|
if($choix eq 'o') {
|
|
572
|
print " => voir fichier $fredir\n";
|
|
573
|
#on ajoute le fichier .log a la liste de fichiers a ne pas effacer
|
|
574
|
push(@liste_fic_a_ne_pas_effacer, $fredir);
|
|
575
|
}
|
|
576
|
|
|
577
|
#destruction des fichiers de calcul
|
|
578
|
efface_fic_herezh();
|
|
579
|
|
|
580
|
#arret du programme
|
|
581
|
die "\nArret du programme a cause d un probleme d execution Herezh++...\n\n";
|
|
582
|
}
|
|
583
|
|
|
584
|
|
|
585
|
#---------------------
|
|
586
|
#lecture du .msh et reecriture pour modifier les couleurs selon le type de reference (noeud, arete, face, element)
|
|
587
|
# et pour saisir la liste des types de reference dans l ordre du fichier .msh pour affecter des options Gmsh suivant le type de reference
|
|
588
|
#---------------------
|
|
589
|
|
|
590
|
#on affecte une isovaleur par couleur :
|
|
591
|
# - 0 => gris (dedie a l affichage du maillage)
|
|
592
|
# - 1 => bleu (dedie a l affichage des refs d elements)
|
|
593
|
# - 2 => vert (dedie a l affichage des refs de faces)
|
|
594
|
# - 3 => jaune (dedie a l affichage des refs d aretes)
|
|
595
|
# - 4 => rouge (dedie a l affichage des refs de noeuds)
|
|
596
|
# - 5 => mauve (dedie a l affichage des refs de points d integration)
|
|
597
|
my $couleur_RGB_maillage = '{190, 190, 190, 255}';#gris
|
|
598
|
my $couleur_RGB_element = '{ 0, 150, 255, 255}';#bleu
|
|
599
|
my $couleur_RGB_face = '{ 0, 255, 0, 255}';#vert
|
|
600
|
my $couleur_RGB_arete = '{240, 200, 0, 255}';#jaune
|
|
601
|
my $couleur_RGB_noeud = '{255, 0, 0, 255}';#rouge
|
|
602
|
my $couleur_RGB_pti = '{238, 130, 238, 255}';#mauve
|
|
603
|
|
|
604
|
#jeu de couleur par type de reference (on utilise les isovaleur gmsh qui vont du bleu fonce au rouge fonce en passant par le vert et le jaune)
|
|
605
|
# rq : pour donner une idee des couleur, si on fixe les bornes d isovaleur entre [1:4], on a : 1.7=bleu clair, 2.5=vert, 3=>jaune, 3.6=rouge fonce
|
|
606
|
my %table_couleur_type_ref = ('pt_integr'=>5, 'noeud'=>4, 'arete'=>3, 'face'=>2, 'element'=>1);
|
|
607
|
|
|
608
|
#on va lister les types de ref
|
|
609
|
my @liste_type_reference;#liste des types de references
|
|
610
|
|
|
611
|
#fichier temporaire
|
|
612
|
my $ftmp = $fGmsh.rand(9999999); while(-e $ftmp) {$ftmp = $fGmsh.rand(9999999);}
|
|
613
|
open(FIC, "<$fGmsh");
|
|
614
|
open(Ftmp, ">$ftmp");
|
|
615
|
my $is_Element_data = 0;
|
|
616
|
my $couleur_type_ref;
|
|
617
|
print " modification du fichier _Gmsh.msh...\n";
|
|
618
|
while(my $ligne = <FIC>) {
|
|
619
|
|
|
620
|
#cas d un element data => on le recopie si il s agit d une reference existant reellement dans les maillages (avec la bonne isovaleur suivant le type de reference)
|
|
621
|
if($ligne =~ /^\s*\$ElementData\s*$/i) {
|
|
622
|
my $entete = $ligne;
|
|
623
|
#on lit jusqu au nom de la reference
|
|
624
|
while($ligne = <FIC>) {
|
|
625
|
$entete .= $ligne;
|
|
626
|
last if($ligne =~ /^\s*\"/);
|
|
627
|
}
|
|
628
|
|
|
629
|
#selection de l isovaleur en fonction du type de reference
|
|
630
|
#-ref de noeuds
|
|
631
|
if($ligne =~ /^\s*\"\s*N(\S+)/) {
|
|
632
|
push(@liste_type_reference, 'noeud');
|
|
633
|
$couleur_type_ref = $table_couleur_type_ref{'noeud'};
|
|
634
|
}
|
|
635
|
#-ref d aretes
|
|
636
|
elsif($ligne =~ /^\s*\"\s*A(\S+)/) {
|
|
637
|
push(@liste_type_reference, 'arete');
|
|
638
|
$couleur_type_ref = $table_couleur_type_ref{'arete'};
|
|
639
|
}
|
|
640
|
#-ref de faces
|
|
641
|
elsif($ligne =~ /^\s*\"\s*F(\S+)/) {
|
|
642
|
push(@liste_type_reference, 'face');
|
|
643
|
$couleur_type_ref = $table_couleur_type_ref{'face'};
|
|
644
|
}
|
|
645
|
#-ref d elements
|
|
646
|
elsif($ligne =~ /^\s*\"\s*E(\S+)/) {
|
|
647
|
push(@liste_type_reference, 'element');
|
|
648
|
$couleur_type_ref = $table_couleur_type_ref{'element'};
|
|
649
|
}
|
|
650
|
#-ref de points d integration
|
|
651
|
elsif($ligne =~ /^\s*\"\s*G(\S+)/) {
|
|
652
|
push(@liste_type_reference, 'pt_integr');
|
|
653
|
$couleur_type_ref = $table_couleur_type_ref{'pt_integr'};
|
|
654
|
}
|
|
655
|
|
|
656
|
#si le nom est une reference d elements creee par ce script => on ne l ecrit pas, on lit jusqu a la fin de l element data et on supprime cette ref de la liste @liste_type_reference
|
|
657
|
if($ligne =~ /E_tmp_visu_elem_/) {
|
|
658
|
while($ligne = <FIC>) {last if($ligne =~ /^\s*\$EndElementData\s*$/i);}
|
|
659
|
pop(@liste_type_reference);#suppression du dernier element de la liste
|
|
660
|
}
|
|
661
|
#si c est une reference reellement dans le maillage, on recopie l en-tete actuel et on recopie l element data avec l isovaleur
|
|
662
|
else {
|
|
663
|
$is_Element_data = 1;
|
|
664
|
print Ftmp $entete;
|
|
665
|
while($ligne = <FIC>) {
|
|
666
|
$ligne = "$1 $couleur_type_ref\n" if ($ligne =~ /^\s*(\d+)\s+\d+\s*$/);
|
|
667
|
print Ftmp $ligne;
|
|
668
|
last if($ligne =~ /^\s*\$EndElementData\s*$/i);
|
|
669
|
}
|
|
670
|
}
|
|
671
|
}#if($ligne =~ /^\s*\$ElementData\s*$/i)
|
|
672
|
|
|
673
|
#cas general : on recopie simplement la ligne courante
|
|
674
|
else {
|
|
675
|
print Ftmp $ligne;
|
|
676
|
}
|
|
677
|
}
|
|
678
|
close(FIC);
|
|
679
|
|
|
680
|
#rajout d une ref supplementaire contenant tous les elements pour affichage des faces des elements 2D et 3D (on leur affecte l isovaleur 0, c est a dire couleur gris)
|
|
681
|
# rq : les elements 1D se retrouvent egalement dans cette ref
|
|
682
|
print Ftmp "\$ElementData
|
|
683
|
1
|
|
684
|
\"Activer/Desactiver vue elements 2D 3D\"
|
|
685
|
0
|
|
686
|
3
|
|
687
|
0
|
|
688
|
1
|
|
689
|
$nb_elts_tot\n";
|
|
690
|
for(my $i=1; $i<=$nb_elts_tot; $i++) {
|
|
691
|
print Ftmp "$i 0\n";#isovaleur 0
|
|
692
|
}
|
|
693
|
print Ftmp "\$EndElementData\n";
|
|
694
|
close(Ftmp);
|
|
695
|
system("cp -nf $fGmsh a.msh");
|
|
696
|
system("mv -f $ftmp $fGmsh");
|
|
697
|
|
|
698
|
|
|
699
|
#---------------------
|
|
700
|
#lancement de la visu Gmsh
|
|
701
|
#---------------------
|
|
702
|
#-on cree un fichier .geo pour y ecrire dans l ordre :
|
|
703
|
# - des options generales a toutes les vues
|
|
704
|
# - un Merge du .msh
|
|
705
|
# - des options au cas par cas par type de reference
|
|
706
|
open(FIC, ">$fgeo");
|
|
707
|
|
|
708
|
#variable donnant le nom du fichier .msh
|
|
709
|
print FIC "//nom du fichier .msh\n";
|
|
710
|
print FIC "fichier_msh = \"$fGmsh\";\n\n";
|
|
711
|
|
|
712
|
#-options generales
|
|
713
|
print FIC '
|
|
714
|
Geometry.Light = 0; //desactivation de la lumiere (geometrie)
|
|
715
|
Mesh.Light = 0; //desactivation de la lumiere (maillage)
|
|
716
|
View.Light = 0; //desactivation de la lumiere (vue)
|
|
717
|
Mesh.ColorCarousel = 0; //type de couleur (0=by element type, 1=by elementary entity, 2=by physical entity, 3=by partition)
|
|
718
|
Geometry.Points = 0; //affichage des points (=0 desactiver, =1 activer)
|
|
719
|
Mesh.Points = 0; //affichage des noeuds (=0 desactiver, =1 activer)
|
|
720
|
Mesh.Lines = 1; //affichage des lignes (elements 1D) (=0 desactiver, =1 activer)
|
|
721
|
Mesh.SurfaceEdges = 1; //affichage des aretes des elements 2D (=0 desactiver, =1 activer)
|
|
722
|
Mesh.SurfaceFaces = 0; //affichage des faces des elements 2D (=0 desactiver, =1 activer)
|
|
723
|
Mesh.VolumeEdges = 1; //affichage des aretes des elements 3D (=0 desactiver, =1 activer)
|
|
724
|
Mesh.VolumeFaces = 0; //affichage des faces des elements 3D (=0 desactiver, =1 activer)
|
|
725
|
View.Visible = 0; //desactivation de toutes les vues au demarrage
|
|
726
|
View.ShowScale = 0; //desactivation de la vue de l echelle d isovaleur
|
|
727
|
View.RangeType = 2; //type de bornes des isovaleurs (2=Custom)
|
|
728
|
View.CustomMin = 0; //borne mini isovaleur (rappel de la convention adoptee pour $NOM_PROG : 0=>maillage, 1=>element, 2=>face, 3=>arete, 4=>noeud, 5=>point integration)
|
|
729
|
View.CustomMax = 5; //borne maxi isovaleur
|
|
730
|
View.PointType = 1; //type d affichage des points (1=3D sphere)
|
|
731
|
View.PointSize = 4.; //taille des points
|
|
732
|
View.LineType = 1; //type d affichage des lignes (1=3D cylinder)
|
|
733
|
View.LineWidth = 3.; //taille des lignes
|
|
734
|
|
|
735
|
Mesh.PointSize = 3.; //taille des noeuds
|
|
736
|
Mesh.Color.Points = {0,0,0}; //on met toutes les couleurs des elements en noir pour que les noeuds et les aretes des elements soient en noir
|
|
737
|
Mesh.Color.Lines = {0,0,0}; // la couleur des faces des elements 2D et 3D sera geree par l isovaleur 0 (gris) de la derniere View intitulee : Activer/Desactiver vue elements 2D 3D
|
|
738
|
Mesh.Color.Triangles = {0,0,0};
|
|
739
|
Mesh.Color.Quadrangles = {0,0,0};
|
|
740
|
Mesh.Color.Tetrahedra = {0,0,0};
|
|
741
|
Mesh.Color.Hexahedra = {0,0,0};
|
|
742
|
Mesh.Color.Prisms = {0,0,0};
|
|
743
|
Mesh.Color.Pyramids = {0,0,0};
|
|
744
|
';
|
|
745
|
|
|
746
|
#echelle de couleur des isovaleurs
|
|
747
|
print FIC "//couleur isovaleur : couleur 1 (gris) => isovaleur=0 (maillage)\n";
|
|
748
|
print FIC "// couleur 2 (bleu) => isovaleur=1 (ref elements)\n";
|
|
749
|
print FIC "// couleur 3 (vert) => isovaleur=2 (ref faces)\n";
|
|
750
|
print FIC "// couleur 4 (jaune) => isovaleur=3 (ref aretes)\n";
|
|
751
|
print FIC "// couleur 5 (rouge) => isovaleur=4 (ref noeuds)\n";
|
|
752
|
print FIC "// couleur 6 (mauve) => isovaleur=5 (ref points integration)\n";
|
|
753
|
# valeur isovaleur 0 1 2 3 4 5
|
|
754
|
print FIC "View.ColorTable = {$couleur_RGB_maillage, $couleur_RGB_element, $couleur_RGB_face, $couleur_RGB_arete, $couleur_RGB_noeud, $couleur_RGB_pti};\n";
|
|
755
|
|
|
756
|
#-Merge du fichier .msh (rq : on passe par Sprintf pour utiliser la variable fichier_msh definie au debut)
|
|
757
|
print FIC "\nMerge Sprintf(fichier_msh);\n\n";
|
|
758
|
|
|
759
|
#on indique d afficher la vue speciale qui sert a afficher les faces des elements 2D 3D (la derniere qui a ete cree)
|
|
760
|
# rq : pour cette vue, on remet l affichage classique pour les points et les lignes
|
|
761
|
print FIC "//options speciales pour la derniere vue qui sert a l affichage des faces des elements 2D et 3D\n";
|
|
762
|
$_ = $#liste_type_reference + 1;#cette vue n est pas enregistree dans la liste, son numero est donc egal a la derniere + 1
|
|
763
|
print FIC "View[$_].Visible = 1;\n";
|
|
764
|
print FIC "View[$_].PointType = 0; //type d affichage des points (0=Color dot)\n";
|
|
765
|
print FIC "View[$_].PointSize = 3.; //taille des points\n";
|
|
766
|
print FIC "View[$_].LineType = 0; //type d affichage des lignes (0=Color segment)\n";
|
|
767
|
print FIC "View[$_].LineWidth = 1.; //taille des lignes\n";
|
|
768
|
print FIC "View[$_].Explode = 0.999; //on reduit legerement la taille de cette vue pour eviter un conflit de couleur quand on affiche des ref de faces ou d elements par dessus\n";
|
|
769
|
|
|
770
|
close(FIC);
|
|
771
|
|
|
772
|
print " visu Gmsh en cours (via fichiers $fgeo et $fGmsh)...\n";
|
|
773
|
system("$exeGMSH $fgeo");
|
|
774
|
|
|
775
|
#cas de l option -saveVisu (sauvegarde des fichiers .geo et _Gmsh.msh)
|
|
776
|
if($is_opt_saveVisu) {
|
|
777
|
#on ajoute les fichiers .geo et _Gmsh.msh a la liste des fichiers a ne pas effacer
|
|
778
|
push(@liste_fic_a_ne_pas_effacer, $fgeo, $fGmsh);
|
|
779
|
print "\nopt -saveVisu => Les fichiers $fgeo et $fGmsh ont ete sauvegardes...\n\n";
|
|
780
|
}
|
|
781
|
|
|
782
|
#destruction des fichiers de calcul
|
|
783
|
efface_fic_herezh();
|
|
784
|
|
|
785
|
|
|
786
|
|
|
787
|
|
|
788
|
|
|
789
|
|
|
790
|
#rq : cette subroutine utilise les variables globales $racine_fic_tmp et @liste_fic_a_ne_pas_effacer
|
|
791
|
sub efface_fic_herezh {
|
|
792
|
#on transforme la liste des fichiers a ne pas effacer en table d indicateur (liste @liste_fic_a_ne_pas_effacer)
|
|
793
|
my %NE_PAS_EFFACER;
|
|
794
|
foreach my $fic (@liste_fic_a_ne_pas_effacer) {$NE_PAS_EFFACER{$fic} = 1;}
|
|
795
|
#on saisit la liste actuelle des fichiers qui commencent comme le fichier .info
|
|
796
|
my @liste_fic_actuelle = glob("$racine_fic_tmp*");
|
|
797
|
#on efface seulement ceux qui n ont pas d indicateur %NE_PAS_EFFACER
|
|
798
|
foreach my $fic (@liste_fic_actuelle) {
|
|
799
|
system("rm -rf $fic") if(not defined($NE_PAS_EFFACER{$fic}));
|
|
800
|
}
|
|
801
|
}
|
|
802
|
|
|
803
|
|
|
804
|
#surcouche de which pour traiter le cas d une commande en chemin absolu (cas qui apparemment ne marche pas avec la sub which d origine du package File::Which ??!!?!??!!!)
|
|
805
|
# N EST PLUS UTILISE DEPUIS LA VERSION 1.01 => A LA PLACE, ON UTILISE LA SUB verif_cmd()
|
|
806
|
##sub which_absolute {
|
|
807
|
## use File::Which;
|
|
808
|
## my $cmd = shift;
|
|
809
|
## my @path;
|
|
810
|
## #cas d une commande avec chemin absolu
|
|
811
|
## if($cmd =~ /^\// and -x $cmd) {
|
|
812
|
## @path = ($cmd);
|
|
813
|
## }
|
|
814
|
## #commande which classique
|
|
815
|
## push(@path, which($cmd));
|
|
816
|
## return(@path);
|
|
817
|
##}
|
|
818
|
|
|
819
|
#cette subroutine verifie l existence d une commande dans $PATH (remplace l utilisation de which depuis la version 1.01)
|
|
820
|
sub verif_cmd {
|
|
821
|
my $cmd = shift;
|
|
822
|
|
|
823
|
#verif directe : est-ce que le fichier existe et est executable
|
|
824
|
return 1 if(-x $cmd);
|
|
825
|
|
|
826
|
#sinon, on regarde dans les path
|
|
827
|
foreach my $path (split(/\s*:\s*/, $ENV{PATH})) {
|
|
828
|
return 1 if(-x "$path/$cmd");
|
|
829
|
}
|
|
830
|
|
|
831
|
#cas ou la commande n existe pas
|
|
832
|
return 0;
|
|
833
|
}
|
|
834
|
|
|
835
|
|
|
836
|
#ecrire une liste de noeuds ou d elements
|
|
837
|
# rq : a priori, si on met trop de noeuds/elements sur une meme ligne, Herezh plante (donc on se limite 15 nombres par ligne)
|
|
838
|
sub ecrire_liste_N_E {
|
|
839
|
my $handle = shift;#le handle de fichier est passe par une variable
|
|
840
|
my $nom_liste = shift;
|
|
841
|
my @liste_no = @_;#liste des noeuds ou d elements
|
|
842
|
|
|
843
|
my $cpt; my $cpt_max = 15; my $nb_blancs;
|
|
844
|
|
|
845
|
$nb_blancs = ""; $nb_blancs .= " " for(1 .. length($nom_liste));
|
|
846
|
$_ = shift(@liste_no);
|
|
847
|
print $handle " $nom_liste $_";
|
|
848
|
$cpt = 1;
|
|
849
|
foreach my $no (@liste_no) {
|
|
850
|
$cpt++;
|
|
851
|
if($cpt == 1) {print $handle " $nb_blancs $no";}
|
|
852
|
elsif($cpt == $cpt_max) {print $handle " $no\n"; $cpt = 0;}
|
|
853
|
else {print $handle " $no";}
|
|
854
|
}
|
|
855
|
print $handle "\n" if($cpt != $cpt_max);
|
|
856
|
}#sub ecrire_liste_noeuds
|
|
857
|
|
|
858
|
|
|
859
|
#----------------
|
|
860
|
#sub qui lit un maillage herezh++ pour recuperer les noeuds, les elements et les listes de references
|
|
861
|
#et les renvoier sous forme de reference (lecture du .her et d un .lis si il existe)
|
|
862
|
#
|
|
863
|
# exemple d appel :
|
|
864
|
# my ($nom_maillage, $nb_noeuds, $ref_tab_noeuds, $nb_elts, $ref_tab_elements, @ref_listes) = lecture_mail_her("fic_her.her");
|
|
865
|
#
|
|
866
|
# avec - $nom_maillage : nom du maillage (si il y en a un. sinon $nom_maillage sera egal a undef
|
|
867
|
# - $nb_noeuds : nombre de noeuds (entier)
|
|
868
|
# - $ref_tab_noeuds : reference vers un tableau de noeuds => $ref_tab_noeuds->[no noeud][0] : coordonnee x
|
|
869
|
# $ref_tab_noeuds->[no noeud][1] : coordonnee y
|
|
870
|
# $ref_tab_noeuds->[no noeud][2] : coordonnee z)
|
|
871
|
# - $nb_elts : nombre d elements (entier)
|
|
872
|
# - $ref_tab_elements : reference vers une table de hashage => $ref_tab_elements->{no elt}{'TYPE'} : type d element
|
|
873
|
# @{$ref_tab_elements->{no elt}{'CONNEX'}} : (liste des noeuds)
|
|
874
|
# - @ref_listes : liste de references vers les tables de hashage contenant les listes de references de noeuds, aretes, faces et elements
|
|
875
|
# => $ref_listes[0] : reference vers la table de hashage des listes de noeuds => @{$ref_listes[0]->{'nom liste'}} : (liste des noeuds)
|
|
876
|
# $ref_listes[1] : reference vers la table de hashage des listes d aretes => @{$ref_listes[1]->{'nom liste'}} : (liste des aretes)
|
|
877
|
# $ref_listes[2] : reference vers la table de hashage des listes de faces => @{$ref_listes[2]->{'nom liste'}} : (liste des faces)
|
|
878
|
# $ref_listes[3] : reference vers la table de hashage des listes d elements => @{$ref_listes[3]->{'nom liste'}} : (liste des elements)
|
|
879
|
#
|
|
880
|
sub lecture_mail_her {
|
|
881
|
my $fher = shift;
|
|
882
|
|
|
883
|
$fher .= '.her' if(not $fher =~ /\.her$/);
|
|
884
|
|
|
885
|
my $nom_maillage;
|
|
886
|
|
|
887
|
#------------------------
|
|
888
|
# lecture du maillage .her
|
|
889
|
#------------------------
|
|
890
|
#-lecture de noeuds
|
|
891
|
my @tab_noeuds; my $nb_noeuds;
|
|
892
|
my $no_noeud;
|
|
893
|
open(Fher, "<$fher") or die "\nErreur (sub:lecture_mail_her) : impossible d ouvrir le fichier $fher ...\n\n";
|
|
894
|
while(<Fher>) {
|
|
895
|
if(/^\s*nom_maillage\s+(\S+)/) {$nom_maillage = $1; next;}
|
|
896
|
next if(not /(\d+)\s+NOEUDS/);
|
|
897
|
$nb_noeuds = $1;
|
|
898
|
last;
|
|
899
|
}
|
|
900
|
while(<Fher>) {
|
|
901
|
last if($no_noeud == $nb_noeuds);
|
|
902
|
next if(not /^\s*(\d+)\s+(\S+)\s+(\S+)\s+(\S+)\s*$/);
|
|
903
|
$no_noeud = $1;
|
|
904
|
@{$tab_noeuds[$no_noeud]} = ($2,$3,$4);
|
|
905
|
}
|
|
906
|
|
|
907
|
#-lecture des elements
|
|
908
|
my %tab_elements; my $nb_elts;
|
|
909
|
my $no_elt;
|
|
910
|
while(<Fher>) {
|
|
911
|
next if(not /(\d+)\s+ELEMENTS/);
|
|
912
|
$nb_elts = $1;
|
|
913
|
last;
|
|
914
|
}
|
|
915
|
while(<Fher>) {
|
|
916
|
last if($no_elt == $nb_elts);
|
|
917
|
next if(not /^\s*\d+\s+\w+\s+\w+/);
|
|
918
|
s/^\s+//;s/\s+$//;
|
|
919
|
$_ =~ /^(\d+)\s+/;
|
|
920
|
$no_elt = $1; s/^(\d+)\s+//;
|
|
921
|
$_ =~ /\s+(\d+(?:\s+\d+)*)$/;
|
|
922
|
@{$tab_elements{$no_elt}{'CONNEX'}} = split(/\s+/, $1); s/\s+(\d+(?:\s+\d+)*)$//;
|
|
923
|
$tab_elements{$no_elt}{'TYPE'} = $_; $tab_elements{$no_elt}{'TYPE'} =~ s/\s+/ /g;
|
|
924
|
}
|
|
925
|
close(Fher);
|
|
926
|
|
|
927
|
|
|
928
|
#------------------------
|
|
929
|
# lecture des references (dans le .her et dans un eventuel .lis)
|
|
930
|
#------------------------
|
|
931
|
my $flis = $fher; $flis =~ s/.her$/.lis/;
|
|
932
|
my $nom_liste;
|
|
933
|
my $is_liste_en_cours;
|
|
934
|
my %listes_NOEUDS;
|
|
935
|
my %listes_ARETES;
|
|
936
|
my %listes_FACES;
|
|
937
|
my %listes_ELEMENTS;
|
|
938
|
|
|
939
|
#-dans le .her
|
|
940
|
open(Fher, "<$fher");
|
|
941
|
$is_liste_en_cours = 0;
|
|
942
|
while(<Fher>) {
|
|
943
|
chomp;
|
|
944
|
if(/^\s*(N\S+)/) {
|
|
945
|
$nom_liste = $1;
|
|
946
|
$is_liste_en_cours = 1;
|
|
947
|
s/^\s*N\S+\s+//; s/\s+$//;
|
|
948
|
push(@{$listes_NOEUDS{$nom_liste}},split(/\s+/,$_));
|
|
949
|
}
|
|
950
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[AFE]/) {
|
|
951
|
$is_liste_en_cours = 0;
|
|
952
|
}
|
|
953
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
954
|
s/^\s+//; s/\s+$//;
|
|
955
|
push(@{$listes_NOEUDS{$nom_liste}},split(/\s+/,$_));
|
|
956
|
}
|
|
957
|
}
|
|
958
|
close(Fher);
|
|
959
|
|
|
960
|
open(Fher, "<$fher");
|
|
961
|
$is_liste_en_cours = 0;
|
|
962
|
while(<Fher>) {
|
|
963
|
chomp;
|
|
964
|
if(/^\s*(A\S+)/) {
|
|
965
|
$nom_liste = $1;
|
|
966
|
$is_liste_en_cours = 1;
|
|
967
|
s/^\s*A\S+\s+//; s/\s+$//;
|
|
968
|
push(@{$listes_ARETES{$nom_liste}},split(/\s+/,$_));
|
|
969
|
}
|
|
970
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NFE]/) {
|
|
971
|
$is_liste_en_cours = 0;
|
|
972
|
}
|
|
973
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
974
|
s/^\s+//; s/\s+$//;
|
|
975
|
push(@{$listes_ARETES{$nom_liste}},split(/\s+/,$_));
|
|
976
|
}
|
|
977
|
}
|
|
978
|
close(Fher);
|
|
979
|
|
|
980
|
open(Fher, "<$fher");
|
|
981
|
$is_liste_en_cours = 0;
|
|
982
|
while(<Fher>) {
|
|
983
|
chomp;
|
|
984
|
if(/^\s*(F\S+)/) {
|
|
985
|
$nom_liste = $1;
|
|
986
|
$is_liste_en_cours = 1;
|
|
987
|
s/^\s*F\S+\s+//; s/\s+$//;
|
|
988
|
push(@{$listes_FACES{$nom_liste}},split(/\s+/,$_));
|
|
989
|
}
|
|
990
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NAE]/) {
|
|
991
|
$is_liste_en_cours = 0;
|
|
992
|
}
|
|
993
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
994
|
s/^\s+//; s/\s+$//;
|
|
995
|
push(@{$listes_FACES{$nom_liste}},split(/\s+/,$_));
|
|
996
|
}
|
|
997
|
}
|
|
998
|
close(Fher);
|
|
999
|
|
|
1000
|
open(Fher, "<$fher");
|
|
1001
|
$is_liste_en_cours = 0;
|
|
1002
|
while(<Fher>) {
|
|
1003
|
chomp;
|
|
1004
|
if(/^\s*(E\S+)/) {
|
|
1005
|
$nom_liste = $1;
|
|
1006
|
$is_liste_en_cours = 1;
|
|
1007
|
s/^\s*E\S+\s+//; s/\s+$//;
|
|
1008
|
push(@{$listes_ELEMENTS{$nom_liste}},split(/\s+/,$_));
|
|
1009
|
}
|
|
1010
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NAF]/) {
|
|
1011
|
$is_liste_en_cours = 0;
|
|
1012
|
}
|
|
1013
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
1014
|
s/^\s+//; s/\s+$//;
|
|
1015
|
push(@{$listes_ELEMENTS{$nom_liste}},split(/\s+/,$_));
|
|
1016
|
}
|
|
1017
|
}
|
|
1018
|
close(Fher);
|
|
1019
|
|
|
1020
|
|
|
1021
|
#dans le .lis (si il existe)
|
|
1022
|
if(-e $flis) {
|
|
1023
|
|
|
1024
|
open(Flis, "<$flis") or die "\nErreur (sub:lecture_mail_her) : impossible d ouvrir le fichier $flis ...\n\n";
|
|
1025
|
$is_liste_en_cours = 0;
|
|
1026
|
while(<Flis>) {
|
|
1027
|
chomp;
|
|
1028
|
if(/^\s*(N\S+)/) {
|
|
1029
|
$nom_liste = $1;
|
|
1030
|
$is_liste_en_cours = 1;
|
|
1031
|
s/^\s*N\S+\s+//; s/\s+$//;
|
|
1032
|
push(@{$listes_NOEUDS{$nom_liste}},split(/\s+/,$_));
|
|
1033
|
}
|
|
1034
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[AFE]/) {
|
|
1035
|
$is_liste_en_cours = 0;
|
|
1036
|
}
|
|
1037
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
1038
|
s/^\s+//; s/\s+$//;
|
|
1039
|
push(@{$listes_NOEUDS{$nom_liste}},split(/\s+/,$_));
|
|
1040
|
}
|
|
1041
|
}
|
|
1042
|
close(Flis);
|
|
1043
|
|
|
1044
|
open(Flis, "<$flis");
|
|
1045
|
$is_liste_en_cours = 0;
|
|
1046
|
while(<Flis>) {
|
|
1047
|
chomp;
|
|
1048
|
if(/^\s*(A\S+)/) {
|
|
1049
|
$nom_liste = $1;
|
|
1050
|
$is_liste_en_cours = 1;
|
|
1051
|
s/^\s*A\S+\s+//; s/\s+$//;
|
|
1052
|
push(@{$listes_ARETES{$nom_liste}},split(/\s+/,$_));
|
|
1053
|
}
|
|
1054
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NFE]/) {
|
|
1055
|
$is_liste_en_cours = 0;
|
|
1056
|
}
|
|
1057
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
1058
|
s/^\s+//; s/\s+$//;
|
|
1059
|
push(@{$listes_ARETES{$nom_liste}},split(/\s+/,$_));
|
|
1060
|
}
|
|
1061
|
}
|
|
1062
|
close(Flis);
|
|
1063
|
|
|
1064
|
open(Flis, "<$flis");
|
|
1065
|
$is_liste_en_cours = 0;
|
|
1066
|
while(<Flis>) {
|
|
1067
|
chomp;
|
|
1068
|
if(/^\s*(F\S+)/) {
|
|
1069
|
$nom_liste = $1;
|
|
1070
|
$is_liste_en_cours = 1;
|
|
1071
|
s/^\s*F\S+\s+//; s/\s+$//;
|
|
1072
|
push(@{$listes_FACES{$nom_liste}},split(/\s+/,$_));
|
|
1073
|
}
|
|
1074
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NAE]/) {
|
|
1075
|
$is_liste_en_cours = 0;
|
|
1076
|
}
|
|
1077
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
1078
|
s/^\s+//; s/\s+$//;
|
|
1079
|
push(@{$listes_FACES{$nom_liste}},split(/\s+/,$_));
|
|
1080
|
}
|
|
1081
|
}
|
|
1082
|
close(Flis);
|
|
1083
|
|
|
1084
|
open(Flis, "<$flis");
|
|
1085
|
$is_liste_en_cours = 0;
|
|
1086
|
while(<Flis>) {
|
|
1087
|
chomp;
|
|
1088
|
if(/^\s*(E\S+)/) {
|
|
1089
|
$nom_liste = $1;
|
|
1090
|
$is_liste_en_cours = 1;
|
|
1091
|
s/^\s*E\S+\s+//; s/\s+$//;
|
|
1092
|
push(@{$listes_ELEMENTS{$nom_liste}},split(/\s+/,$_));
|
|
1093
|
}
|
|
1094
|
elsif(/^\s*noeuds/i or /^\s*elements/i or /^\s*[NAF]/) {
|
|
1095
|
$is_liste_en_cours = 0;
|
|
1096
|
}
|
|
1097
|
elsif($is_liste_en_cours and /^\s*\d+(\s+\d+)*\s*$/i) {
|
|
1098
|
s/^\s+//; s/\s+$//;
|
|
1099
|
push(@{$listes_ELEMENTS{$nom_liste}},split(/\s+/,$_));
|
|
1100
|
}
|
|
1101
|
}
|
|
1102
|
close(Flis);
|
|
1103
|
|
|
1104
|
}#if(-e $flis)
|
|
1105
|
|
|
1106
|
#AFFICHAGE DES LISTES DE NOEUDS
|
|
1107
|
#foreach my $nom (keys(%listes_NOEUDS)) {
|
|
1108
|
# print "$nom : @{$listes_NOEUDS{$nom}}\n";
|
|
1109
|
#}
|
|
1110
|
#AFFICHAGE DES LISTES D ARETES
|
|
1111
|
#foreach my $nom (keys(%listes_ARETES)) {
|
|
1112
|
# print "$nom : @{$listes_ARETES{$nom}}\n";
|
|
1113
|
#}
|
|
1114
|
#AFFICHAGE DES LISTES DE FACES
|
|
1115
|
#foreach my $nom (keys(%listes_FACES)) {
|
|
1116
|
# print "$nom : @{$listes_FACES{$nom}}\n";
|
|
1117
|
#}
|
|
1118
|
#AFFICHAGE DES LISTES D ELEMENTS
|
|
1119
|
#foreach my $nom (keys(%listes_ELEMENTS)) {
|
|
1120
|
# print "$nom : @{$listes_ELEMENTS{$nom}}\n";
|
|
1121
|
#}
|
|
1122
|
|
|
1123
|
return($nom_maillage, $nb_noeuds, \@tab_noeuds, $nb_elts, \%tab_elements,
|
|
1124
|
\%listes_NOEUDS, \%listes_ARETES,
|
|
1125
|
\%listes_FACES, \%listes_ELEMENTS);
|
|
1126
|
}#sub lecture_mail_her
|