[Gmsh] Identification numbers
Philippe Ackerer
ackerer at unistra.fr
Wed Feb 12 10:19:03 CET 2014
Hi,
I want to mesh a cylinder (see attached file) and need the identity
numbers of the bottom, top and sides (will be used to define boundary
conditions).
Can you explain the numbering that appears in out[] (see attached file) ??
Thanks
Best regards and congratulations for Gmsh,
Philippe
--
--------------------------------------------------------------
Philippe ACKERER
LHyGeS UMR 7517 UdS-CNRS-Engees
1, rue Blessig
67000 STRASBOURG
tel +33 (0)368 850 561
-------------------------------------------------------------
-------------- next part --------------
// Meshing a cylinder
r1=19.0;
Lcw=0.20;
Lc1 = 0.40;
teta=2.0*3.14159/8.0;
// Then we define some points and some lines using these variables:
Point(1) = { 0.0, 0.0 , 0.0 , Lcw};
Point(2) = { r1, 0.0 , 0.0 , Lc1};
x1=r1*Cos(teta);
y1=r1*Sin(teta);
Point(3) = { x1 , y1 , 0.0 , Lc1};
Circle(1) = {2,1,3};
Line(10) = {1,2};
Line(11) = {3,1};
Line Loop(21) = {10,1,11}; Plane Surface(22) = {21};
//Extrusion over h
h = 15.0;
out[]=Extrude {0,0,h} {Surface{22}; Layers{ {3,5,3}, {0.20,0.80,1.00} };Recombine; };
Printf(" Top=%g", out[0]);
Printf(" Lateral=%g", out[1]);
Printf(" Lateral=%g", out[2]);
Printf(" Lateral=%g", out[3]);