[Gmsh] GMSH, extrude to angle greater than pi?

Steffens-Jr, Alfred P alfred.p.steffens-jr at boeing.com
Tue Apr 8 17:44:59 CEST 2003


Hello,

Thanks again for your help.  I have a 2D surface which I want to extrude to
2*pi, even if there must be a discontinuity between 0 and 2pi.  I notice
that when I try to do that on the display the angle pi to 2pi becomes pi to
0.  Is there a way in GMSH to extrude to a full 2pi?  I have included my
.geo file below as illustration.

Al Steffens
alfred.p.steffens-jr at boeing.com


//
//	tpcyl2.geo
//


scale_small= 4.0;
scale_big= 10.0;

Point(1) = {657.0515,  0, 155.7671, scale_big};
Point(2) = {558.1988,  0, 161.2588, scale_big};
Point(3) = {474.9661,  0, 180.0553, scale_big};
Point(4) = {426.8788,  0, 201.2593, scale_big};
Point(5) = {399.8606,  0, 223.5084, scale_big};
Point(6) = {400.0422,  0, 227.5042, scale_small};
Point(7) = {406.1267,  0, 229.2299, scale_small};
Point(8) = {411.1216,  0, 229.003, scale_small};
Point(9) = {411.3939,  0, 234.9968, scale_small};
Point(10) = {406.399,  0, 235.2237, scale_small};
Point(11) = {397.4537,  0, 236.6311, scale_small};
Point(12) = {379.7899,  0, 244.4407, scale_small};
Point(13) = {378,  0, 247.483, scale_small};
Point(14) = {381.9693,  0, 248.4366, scale_small};
Point(15) = {397.2715,  0, 248.528, scale_small};
Point(16) = {397.2715,  0, 254.6579, scale_small};
Point(17) = {388.2808,  0, 255.0664, scale_small};
Point(18) = {365.6222,  0, 263.1029, scale_small};
Point(19) = {364.8047,  0, 267.1442, scale_small};
Point(20) = {368.9821,  0, 270.9585, scale_small};
Point(21) = {385.9646,  0, 270.1871, scale_small};
Point(22) = {386.1008,  0, 276.2716, scale_small};
Point(23) = {373.2503,  0, 276.7708, scale_small};
Point(24) = {367.3472,  0, 279.041, scale_small};
Point(25) = {362.4885,  0, 282.2649, scale_small};
Point(26) = {362.8516,  0, 290.2566, scale_small};
Point(27) = {371.8877,  0, 290.8472, scale_small};
Point(28) = {369.1177,  0, 295.9781, scale_small};
Point(29) = {363.2146,  0, 297.25, scale_small};
Point(30) = {355.1775,  0, 297.6124, scale_small};
Point(31) = {352.0898,  0, 295.7506, scale_small};


Point(32) = {658.958,  0, 175.7012, scale_big};
Point(33) = {584.8981,  0, 176.0624, scale_big};
Point(34) = {506.4787,  0, 190.6361, scale_big};
Point(35) = {497.8964,  0, 200.0353, scale_big};
Point(36) = {479.8696,  0, 199.8531, scale_big};
Point(37) = {458.2099,  0, 207.8443, scale_big};
Point(38) = {434.5069,  0, 214.9273, scale_big};
Point(39) = {434.5,  0, 221.9654, scale_big};
Point(40) = {434.1433,  0, 228.9582, scale_big};
Point(41) = {438.2753,  0, 231.7736, scale_big};
Point(42) = {427.3774,  0, 234.2707, scale_big};
Point(43) = {426.5151,  0, 215.2903, scale_small};
Point(44) = {419.5677,  0, 216.6069, scale_small};
Point(45) = {418.5228,  0, 237.6761, scale_small};
Point(46) = {406.6259,  0, 240.2186, scale_small};
Point(47) = {405.5356,  0, 260.2887, scale_small};
Point(48) = {392.5944,  0, 261.8776, scale_small};
Point(49) = {391.504,  0, 281.9478, scale_small};
Point(50) = {377.5184,  0, 282.5831, scale_small};
Point(51) = {377.2456,  0, 298.612, scale_small};
Point(52) = {375.6561,  0, 307.6935, scale_small};
Point(53) = {370.7066,  0, 308.9194, scale_small};
Point(54) = {364.7128,  0, 309.1917, scale_small};
Point(55) = {360.7623,  0, 310.3721, scale_small};
Point(56) = {356.5395,  0, 305.5588, scale_small};
Point(57) = {355.359,  0, 301.6083, scale_small};
Point(58) = {350.3641,  0, 301.8352, scale_small};

//
//		lower curve
//
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,5};
CatmullRom(5) = {5,6,7};
Line(6) = {7,8};
Line(7) = {8,9};
Line(8) = {9,10};
Line(9) = {10,11};
Line(10) = {11,12};
CatmullRom(11) = {12,13,14};
Line(12) = {14,15};
Line(13) = {15,16};
Line(14) = {16,17};
Line(15) = {17,18};
CatmullRom(16) = {18,19,20};
Line(17) = {20,21};
Line(18) = {21,22};
Line(19) = {22,23};
CatmullRom(20) = {23,24,25};
Line(21) = {25,26};
Line(22) = {26,27};
CatmullRom(23) = {27,28,29};
Line(24) = {29,30};
Line(25) = {30,31};


//
//		upper curve
//

Line(26) = {32,33};
Line(27) = {33,34};
Line(28) = {34,35};
Line(29) = {35,36};
Line(30) = {36,37};
Line(31) = {37,38};
Line(32) = {38,39};
Line(33) = {39,40};
Line(34) = {40,41};
Line(35) = {41,42};
Line(36) = {42,43};
Line(37) = {43,44};
Line(38) = {44,45};
Line(39) = {45,46};
Line(40) = {46,47};
Line(41) = {47,48};
Line(42) = {48,49};
Line(43) = {49,50};
Line(44) = {50,51};
CatmullRom(45) = {51,52,53};
Line(46) = {53,54};
CatmullRom(47) = {54,55,56};
Line(48) = {56,57};
Line(49) = {57,58};

//
//	top opening
//
Line(50) = {58,31};
//
//	bottom opening
//
Line(51) = {1,32};


//
//	Q1 cross-section surface
//
Line Loop(52) =
{2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,-50,-49,-48
,-47,-46,-45,-44,-43,-42,-41,-40,-39,-38,-37,-36,-35,-34,-33,-32,-31,-30,-29
,-28,-27,-26,-51,1};
Plane Surface(53) = {52};


//
//	Q1,Q2 shell
//
Extrude Surface {53, {0,0,1}, {0,0,0}, 4.1} {
  Recombine ; Layers { 3, 54, 1 } ; 
} ;

//
//	now GMSH will define surfaces for this extrusion
//	and we must use the GUI to view the surface numbers
//


//
//	Region 1  -- bottom opening
//
Physical Surface (1) = {305};

//
//	Region 2  -- top opening
//
Physical Surface (2) = {205};

//
//	Region 3  -- the enclosing surface (other than top and bottom)
//
Physical Surface (3) = {53, 310, 109, 113, 117, 121, 125, 129, 133, 137,
141, 145, 149, 153, 157, 161, 165, 169, 173, 177, 181, 185, 189, 193, 197,
201, 209, 213, 217, 221, 225, 229, 233, 237, 241, 245, 249, 253, 257, 261,
265, 269, 273, 277, 281, 285, 289, 293, 297, 301, 309};

//
//	Region 4 -- the inside (the volume)
//
Physical Volume (4) = {54};

// end GMSH .geo file
//