[Gmsh] Bug in 3D mesher of Gmsh - cylinder inside a block...

[David Colignon]

Hi,

your definition of the two box faces with circle holes is wrong. You should
carrefully read tutorial t4.geo 

also, these circle surfaces 

>Ruled Surface(1000430) = {1000412};
>Ruled Surface(1000431) = {1000417};

are wrong, they are not Ruled Surface but Plane Surface ans should read :

Plane Surface(1000430) = {1000412};
Plane Surface(1000431) = {1000417};


This is probably what you want : 

lc = 0.475;
lc2 = 0.25;
lent = 1.0;
Point(1) = {0, 0, -lent, lc};
Point(2) = {lent, 0, -lent, lc};
Point(3) = {lent, lent, -lent, lc};
Point(4) = {0, lent, -lent, lc};
Point(5) = {0, 0, lent, lc};
Point(6) = {lent, 0, lent, lc};
Point(7) = {lent, lent, lent, lc};
Point(8) = {0, lent, lent, lc};
Point(9) = {0.5, lent, 0, lc2};
Point(10) = {0.5, lent, 0.3, lc2};
Point(11) = {0.8, lent, 0, lc2};
Point(12) = {0.2, lent, 0, lc2};
Point(13) = {0.5, lent, -0.3, lc2};
Point(14) = {0.5, 0, 0.3, lc2};
Point(15) = {0.8, 0, 0, lc2};
Point(16) = {0.2, 0, 0, lc2};
Point(17) = {0.5, 0, -0.3, lc2};
Point(18) = {0.5, 0, 0, lc2};
Line (101) = {1, 2};
Line (102) = {2, 3};
Line (103) = {3, 4};
Line (104) = {4, 1};
Line (105) = {1, 5};
Line (106) = {4, 8};
Line (107) = {8, 5};
Line (108) = {8, 7};
Line (109) = {5, 6};
Line (110) = {7, 6};
Line (111) = {6, 2};
Line (112) = {3, 7};

Line Loop (1000300) = {101, 102, 103, 104};
Plane Surface (300) = {1000300};
Line Loop (1000301) = {108, 110, -109, -107};
Plane Surface (301) = {1000301};
Line Loop (1000302) = {104, 105, -107, -106};
Plane Surface (302) = {1000302};

Line Loop (1000304) = {102, 112, 110, 111};
Plane Surface (304) = {1000304};

Circle(1000408) = {10,9,11};
Circle(1000409) = {11,9,13};
Circle(1000410) = {13,9,12};
Circle(1000411) = {12,9,10};

Line Loop(1000412) = {1000408,1000409,1000410,1000411};
Plane Surface(1000430) = {1000412}; // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Circle(1000413) = {16,18,14};
Circle(1000414) = {14,18,15};
Circle(1000415) = {15,18,17};
Circle(1000416) = {17,18,16};

Line Loop(1000417) = {1000413,1000414,1000415,1000416};
Plane Surface(1000431) = {1000417}; // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Line(1000418) = {17,13};
Line(1000419) = {16,12};
Line(1000420) = {14,10};
Line(1000421) = {15,11};

Line Loop(1000422) = {1000420,-1000411,-1000419,1000413};
Ruled Surface(1000423) = {1000422};
Line Loop(1000424) = {-1000408,-1000420,1000414,1000421};
Ruled Surface(1000425) = {1000424};
Line Loop(1000426) = {1000418,-1000409,-1000421,1000415};
Ruled Surface(1000427) = {1000426};
Line Loop(1000428) = {1000410,-1000419,-1000416,1000418};
Ruled Surface(1000429) = {1000428};


Line Loop (1000303) = {105, 109, 111, -101};
Plane Surface (303) = {1000303,1000417}; //
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Line Loop (1000305) = {103, 106, 108, -112};
Plane Surface (305) = {1000305,1000412}; //
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


Surface Loop(1000432) = {1000431,1000423,1000425,1000430,1000427,1000429};
Volume(1000433) = {1000432};

Surface Loop(1000434) =
{1000427,1000429,305,300,303,302,301,304,1000423,1000425};
Volume(1000435) = {1000434};



Dave