[Gmsh] Now, I am lost :(

[Philippe Ackerer]

Hi,
Following my previous mail, I try to generate a prismatic mesh of a 
cylinder without the 'extrude' command to prescribe the identity value 
of the sides.
I used Transfinite instructions without success (the 3D elements are not 
prisms - see attached file).

Any help is very welcomed. Thanks,
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;
h = 15.0 ;
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 , Lcw};
Point(4) = { 0.0, 0.0  , h , Lcw};
Point(5) = { r1, 0.0   , h , Lc1};
x1=r1*Cos(teta);
y1=r1*Sin(teta);
Point(6) = { x1 , y1 , h , Lcw};

//Horizontal lines
Circle(1) = {2,1,3};
Circle(2) = {5,4,6};
Line(11) = {1,2}; 
Line(12) = {3,1}; 
Line(13) = {4,5}; 
Line(14) = {6,4}; 

//Vertical lines 
Line(21) = {1,4}; 
Line(22) = {2,5}; 
Line(23) = {3,6}; 

Transfinite Line{21,22,23} = 6;
Transfinite Line{11,-12,13,-14} = 15 Using Progression 1.1;
Transfinite Line{1,2} = 6;

// *** Horizontal Surfaces
Line Loop(101) = {11,1,12}; Plane Surface(31) = {101}; Transfinite Surface {31}={1,2,3};
Line Loop(102) = {13,2,14}; Plane Surface(32) = {102}; Transfinite Surface {32}={4,5,6};

// ****Vertical Surfaces
Line Loop(103) = {11,22,-13,-21};  Plane Surface(33) = {103};  Transfinite Surface {33}={1,2,5,4};
Line Loop(104) = {12,21,-14,-23} ; Plane Surface(34) = {104};  Transfinite Surface {34}={3,1,4,6};
Line Loop(105) = {1,23,-2,-22};    Ruled Surface(35) = {105};  Transfinite Surface {35}={2,3,6,5};                

Recombine Surface {33,34,35} ;

Surface Loop(201)= {31,32,33,34,35};
Volume(301)={201};
Recombine Volume {301};
Physical Volume(501)={301};