[Gmsh] 2D unstructured mesh problem
Umberto Paoletti
paoletti at circuit.kuee.kyoto-u.ac.jp
Mon Nov 10 11:49:21 CET 2008
Hello,
I found a problem in this very simple example.
Some corners of the internal hole seem to be skipped by the mesh algorithm.
Why?
Thanks
Umberto
unit = 0.001000;
lc = 0.1*unit;
edge = 2*unit;
holedge = 0.5*unit;
z = 0*unit;
//external edge
x1 = 0*unit;
y1 = 0*unit;
x2 = x1+edge;
y2 = y1;
x3 = x2;
y3 = y2+edge;
x4 = x3-edge;
y4 = y3;
Point(1) = {x1,y1,z, 1*lc};
Point(2) = {x2,y2,z, 1*lc};
Point(3) = {x3,y3,z, 1*lc};
Point(4) = {x4,y4,z, 1*lc};
Line(1) = {1:4,1};
Line Loop(2) = {1};
x5 = 1*unit;
y5 = 1*unit;
x6 = x5+holedge;
y6 = y5;
x7 = x6;
y7 = y6+holedge;
x8 = x7-holedge;
y8 = y7;
//hole edge
Point(5) = {x5,y5,z, 1*lc};
Point(6) = {x6,y6,z, 1*lc};
Point(7) = {x7,y7,z, 1*lc};
Point(8) = {x8,y8,z, 1*lc};
Line(3) = {5:8,5};
Line Loop(4) = {3};
Plane Surface(1) = {2,4};
--
---------------------------------------------------------------------------------------
Umberto Paoletti, Dr. Eng.,
Building A1, Room 402,
Graduate School of Engineering,
Department of Electrical Engineering,
Kyoto University, Nishikyo-ku,
Katsura, 615-8510, Kyoto, Japan
Tel.: +81-(0)75-383-2246
Fax.: +81-(0)75-383-2245
E-mail: paoletti at circuit.kuee.kyoto-u.ac.jp
---------------------------------------------------------------------------------------