[Gmsh] GMSH, extrude to angle greater than pi?
David Colignon
David.Colignon at univ.u-3mrs.fr
Tue Apr 8 18:33:20 CEST 2003
Hi,
I think you should try this :
...
Extrude Surface {53, {0,0,1}, {0,0,0}, 2*Pi/3} {
Recombine ; Layers { 6, 54, 1 } ;
} ;
Extrude Surface {310, {0,0,1}, {0,0,0}, 2*Pi/3} {
Recombine ; Layers { 6, 311 , 1 } ;
} ;
Extrude Surface {567, {0,0,1}, {0,0,0}, 2*Pi/3} {
Recombine ; Layers { 6, 568 , 1 } ;
} ;
...
Dave
A 08:44 08/04/2003 -0700, vous avez écrit :
>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
>//
>
>
>_______________________________________________
>gmsh mailing list
>gmsh at geuz.org
>http://www.geuz.org/mailman/listinfo/gmsh
>
>
Dave