[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