[Gmsh] structured and unstructured meshes bug
Christophe Geuzaine
cag32 at case.edu
Wed Mar 22 02:11:12 CET 2006
Michel Aubourg wrote:
> Re-bonjour,
>
> Avec le fichier joint cela ira mieux.
>
Hi Michel - Unfortunately, that's inherent to the Delaunay technique as
implemented in the default "isotropic" algorithm. You can use the
combined Delaunay+advancing front technique from the "netgen" algorithm
to fix this. (Add e.g. "Mesh.Algorithm3D = 4;" at the beginning of your
file; beware that your volume boundary has to be defined with exterior
pointing normals when using this algorithm.)
Best,
Christophe
> Cordialement
> Michel Aubourg
>
>
> ------------------------------------------------------------------------
>
> wBoite = 8;
> lBoite = 12;
> hBoite = 4;
>
> hSubstrat = 0.635;
>
> wStrip = 0.63;
> lStrip = 7.2;
>
> lPlot = 1.5;
>
> lGap = ( lBoite - lStrip ) / 2 - lPlot;
>
> mSurface = 0.999999 * wStrip / 2;
> mVolume = 0.999999 * 2;
>
> Point(1) = {0,0,0,mSurface};
> Extrude {wStrip/2,0,0} {
> Point{1};
> Layers {{wStrip/2/mSurface},{1.0}};
> }
> Extrude {(wBoite-wStrip)/2,0,0} {
> Point{2};
> Layers {{(wBoite-wStrip)/2/mSurface},{1.0}};
> }
> Extrude {0,lPlot,0} {
> Line{1,2};
> Layers {{lPlot/mSurface},{1.0}};
> }
> Extrude {0,lGap,0} {
> Line{3,7};
> Layers {{lGap/mSurface},{1.0}};
> }
> Extrude {0,lStrip/2,0} {
> Line{11,15};
> Layers {{lStrip/2/mSurface},{1.0}};
> }
>
> volSubstrat[] =
> Extrude {0,0,-hSubstrat} {
> Surface{6,10,14,18,22,26};
> Layers {{hSubstrat/mSurface},{1.0}};
> }
> ;
>
> Point(64) = {0,0,hBoite,mVolume};
>
> Extrude {wBoite/2,0,0} {
> Point{64};
> Layers {{wBoite/2/mVolume},{1.0}};
> }
> Extrude {0,lBoite/2,0} {
> Line{159};
> Layers {{lBoite/2/mVolume},{1.0}};
> }
>
> Line(164) = {1,64};
> Line(165) = {65,3};
> Line(166) = {15,67};
> Line(167) = {66,12};
>
> Line Loop(168) = {9,17,25,166,-162,165};
> Plane Surface(169) = {168};
> Line Loop(170) = {23,166,-160,167,19};
> Plane Surface(171) = {170};
> Line Loop(172) = {167,-20,-12,-4,164,161};
> Plane Surface(173) = {172};
> Line Loop(174) = {1,2,-165,-159,-164};
> Plane Surface(175) = {174};
> Surface Loop(176) = {163,175,-6,-10,169,-18,-26,-171,173,-22,-14};
>
> Volume(177) = {176};
>
> // ligne
> Physical Line(178) = {33};
>
> // zMin
> Physical Surface(179) = {70,114,158,136,48,92};
> // plot
> Physical Surface(180) = {6};
> // gap
> Physical Surface(181) = {14};
> // strip
> Physical Surface(182) = {22};
> // substrat
> Physical Surface(183) = {10,18,26};
> // zMax
> Physical Surface(184) = {163};
>
> // yMin
> Physical Surface(185) = {57,175,35};
> // yMax
> Physical Surface(186) = {153,171,131};
>
> // xMin
> Physical Surface(187) = {47,91,135,173};
> // xMax
> Physical Surface(188) = {61,105,149,169};
>
> // substrat
> Physical Volume(189) = {volSubstrat[]};
> // air
> Physical Volume(190) = {177};
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
--
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine