[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
//