Как я могу получить более регулярную поверхностную четверку, используя gmsh? - PullRequest
0 голосов
/ 20 июня 2019

Я хочу иметь возможность создавать четырехугольную поверхностную сетку, которая является очень правильной (каждая грань имеет, насколько это возможно, одну и ту же область) и выровнена по границе поверхности.

Следующий тестовый файл .geo упрощает тип предполагаемого варианта использования:

lc = 0.1; 

// vertices.
Point(1) = {0, 0, 0, lc};
Point(2) = {0.5, 0, 0, lc};
Point(3) = {1.0, 0, 0, lc};
Point(4) = {1.0, 0.5, 0.5, lc};
Point(5) = {1.0, 1.0, 1.0, lc};
Point(6) = {1.0, 1.5, 0.5, lc};
Point(7) = {1.0, 2.0, 0.0, lc};
Point(8) = {0.5, 2.0, 0.0, lc};
Point(9) = {0.0, 2.0, 0.0, lc};
Point(10) = {0.0, 1.5, 0.5, lc};
Point(11) = {0.0, 1.0, 1.0, lc};
Point(12) = {0.0, 0.5, 0.5, lc};

// curves.
Spline(1) = {1,2,3};
Spline(2) = {3,4,5,6,7};
Spline(3) = {7,8,9};
Spline(4) = {9,10,11,12,1};
Physical Line("bottom") = {1};
Physical Line("top") = {3};
Curve Loop(1) = {2, 3, 4, 1};

//surface.
Transfinite Curve{1} = 20
Transfinite Surface(1) = {2,3,4,1};
Physical Surface("mysurface") = {1};

Когда я загружаю этот файл .geo в графический интерфейс gmsh (v 4.3.0) и запускаю сетку 1D, затем 2D (опция Frontal-Delaunay) и, наконец, команда 2D рекомбинации (опция Blossom), получающаяся поверхностная сетка не такая регулярная:

enter image description here

Журнал консоли показывает:

Info    : Meshing 1D...
Info    : Meshing curve 1 (Nurb)
Info    : Meshing curve 2 (Nurb)
Info    : Meshing curve 3 (Nurb)
Info    : Meshing curve 4 (Nurb)
Info    : Done meshing 1D (0.008326 s)
Info    : 70 vertices 74 elements
Info    : Meshing 2D...
Info    : Meshing surface 1 (Surface, Frontal)
Info    : Done meshing 2D (0.013711 s)
Info    : 272 vertices 538 elements
Info    : Recombining 2D mesh...
Info    : Blossom: 665 internal 62 closed
Info    : Blossom recombination completed (0.012128 s): 230 quads, 0 triangles, 0 invalid quads, 0 quads with Q < 0.1, avg Q = 0.799983, min Q = 0.502415
Info    : Done recombining 2D mesh (0.012205 s)

Я подозреваю, что это может быть из-за моей относительной неопытности с geo / gmsh. Совет ценится.

...