//parametre de maillage:
cl1 = 10;
//hauteur de la canopee sous gmsh:
H=0;
x=3.0;
//ecriture des points
//cir n°0:
//face n°1:
Point(1) = {45.0, 35.0, H, cl1};
Point(2) = {25.0, 35.0, H, cl1};
Point(3) = {25.0, 25.0, H, cl1};
Point(4) = {45.0, 25.0, H, cl1};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};
//face n°2:
Point(5) = {45.0, 5.0, H, cl1};
Point(6) = {25.0, 5.0, H, cl1};
Line(5) = {5, 4};
Line(6) = {3, 6};
Line(7) = {6, 5};
//face n°3:
Point(7) = {65.0, 55.0, H, cl1};
Point(8) = {45.0, 55.0, H, cl1};
Point(9) = {65.0, 5.0, H, cl1};
Line(8) = {7, 8};
Line(9) = {8, 1};
Line(10) = {5, 9};
Line(11) = {9, 7};
//face n°4:
Point(10) = {25.0, 55.0, H, cl1};
Line(12) = {10, 2};
Line(13) = {8, 10};
//face n°5:
Point(11) = {5.0, 55.0, H, cl1};
Point(12) = {5.0, 5.0, H, cl1};
Line(14) = {11, 12};
Line(15) = {12, 6};
Line(16) = {10, 11};
//on ecrit chaque surface du sol:
//face n1:
// les lignes de cette face sont transfinit￩es:
//creation de la lineloop regroupant toutes les lignes de cet face:
Line Loop(1) = {1,2,3,4
};
//creation du trace au sol du toit du batiment:
Plane Surface(1) = {1};
Recombine Surface {1}= 90;
//extrusion de cette surface plusieurs fois de x metre jusqu'a la hauteur du toit
//le toit est a la hauteur:h1 il faut extruder le trace au sol de 0 vers h1 ,
H1=6.0;
//calcul du nombre d'etage: nbr_etage1=H1/x
local0[]=Extrude {0,0,3.0} {
Surface{1}; Layers{1};Recombine;};
Physical Surface ('plancher_local0')={1};Physical Volume ('vol_local0')={local0[1]};
Physical Surface ('plafond_local0')={local0[0]};
Physical Surface ('surf_local0')={local0[2],local0[3],local0[4],local0[5]};
local1[]=Extrude {0,0,3.0} {
Surface{local0[0]}; Layers{1};Recombine;};
Physical Volume ('vol_local1')={local1[1]};
Physical Surface ('surf_local1')={local1[2],local1[3],local1[4],local1[5]};
Physical Surface ('plafond_local1')={local1[0]};canope0[]=Extrude {0,0,9.0} {
Surface{local1[0]}; Layers{3.0};Recombine;};
//face n2:
// les lignes de cette face sont transfinit￩es:
//creation de la lineloop regroupant toutes les lignes de cet face:
Line Loop(2) = {5,-3,6,7
};
//creation du trace au sol du toit du batiment:
Plane Surface(2) = {2};
Recombine Surface {2}= 90;
//extrusion de cette surface plusieurs fois de x metre jusqu'a la hauteur de la canopee
canope1[]=Extrude {0,0,15} {
Surface{2}; Layers{5.0};Recombine;};
//face n3:
// les lignes de cette face sont transfinit￩es:
//creation de la lineloop regroupant toutes les lignes de cet face:
Line Loop(3) = {8,9,-4,-5,10,11
};
//creation du trace au sol du toit du batiment:
Plane Surface(3) = {3};
Recombine Surface {3}= 90;
//extrusion de cette surface plusieurs fois de x metre jusqu'a la hauteur de la canopee
canope2[]=Extrude {0,0,15} {
Surface{3}; Layers{5.0};Recombine;};
//face n4:
// les lignes de cette face sont transfinit￩es:
//creation de la lineloop regroupant toutes les lignes de cet face:
Line Loop(4) = {12,-1,-9,13
};
//creation du trace au sol du toit du batiment:
Plane Surface(4) = {4};
Recombine Surface {4}= 90;
//extrusion de cette surface plusieurs fois de x metre jusqu'a la hauteur de la canopee
canope3[]=Extrude {0,0,15} {
Surface{4}; Layers{5.0};Recombine;};
//face n5:
// les lignes de cette face sont transfinit￩es:
//creation de la lineloop regroupant toutes les lignes de cet face:
Line Loop(5) = {14,15,-6,-2,-12,16
};
//creation du trace au sol du toit du batiment:
Plane Surface(5) = {5};
Recombine Surface {5}= 90;
//extrusion de cette surface plusieurs fois de x metre jusqu'a la hauteur de la canopee
canope4[]=Extrude {0,0,15} {
Surface{5}; Layers{5.0};Recombine;};
Physical Surface ('sol')={2,3,4,5};
Physical Volume ('canope')={canope0[],canope1[],canope2[],canope3[],canope4[]};
Physical Surface('sol_mineral') = {5, 4, 3, 2};
Physical Surface('toit') = {60};
Physical Surface('mur') = {55, 33, 59, 37, 47, 25, 51, 29};
