[Gmsh] how to load a pre-defined 2D mesh

Christophe Geuzaine cag32 at case.edu
Fri Feb 17 15:37:24 CET 2006


Romain Quey wrote:
> Hello,
> 
> I want to mesh volumes with pre-defined 2D-meshes on some of their
> faces, but it doesn't work ...
> 
> Please find attached a simple example based on a cube:
> 'cube.geo' first defines 5 of the faces, then load the mesh
> of the 6th one (from file 'z1.msh', obtained from 'z1.geo')
> and finally defines the volume ...


Hello Romain - Meshing with partially pre-meshed parts is not supported
at the moment. (It might be added in the future, but it's not high on my
own priority list--maybe on someone else's...)



> 
> When I run the 3D meshing (Netgen algorithm), Gmsh returns:
> Info    : 'gmsh -3 cube.geo ' started on Thu Feb 16 16:34:07 2006
> Info    : Parsing file 'cube.geo'
> Info    : Reading mesh file 'z1.msh'
> Info    : Detected mesh file format 2
> Info    : 134 nodes
> Info    : 266 elements
> Info    : Read mesh file 'z1.msh'
> Info    : Parsed file 'cube.geo'
> Info    : Mesh 3D...
> Info    : Meshing volume 1
> Info    : CalcLocalH: 134 Points 0 Elements 230 Surface Elements
> Info    : Check subdomain 1 / 1
> Info    : 230 open elements
> Fatal   : ERROR: Edge 2 - 5 multiple times in surface mesh
> Aborted
> 
> Where am I wrong?
> 
> Thank you for your help ... and for Gmsh too!
> 
> -- 
> Romain QUEY
> 
> 
> ------------------------------------------------------------------------
> 
> Point(1) = {0,0,0, 0.2};
> Point(2) = {1,0,0, 0.2};
> Point(3) = {1,1,0, 0.2};
> Point(4) = {0,1,0, 0.2};
> Point(5) = {0,0,1, 0.2};
> Point(6) = {0,1,1, 0.2};
> Point(7) = {1,1,1, 0.2};
> Point(8) = {1,0,1, 0.2};
> Line(1) = {4,6} ;
> Line(2) = {6,5} ;
> Line(3) = {5,1} ;
> Line(4) = {1,4} ;
> Line(5) = {2,8} ;
> Line(6) = {8,7} ;
> Line(7) = {7,3} ;
> Line(8) = {3,2} ;
> Line(9) = {5,8} ;
> Line(10) = {2,1} ;
> Line(11) = {4,3} ;
> Line(12) = {7,6} ;
> 
> Line Loop(1) = {1,2,3,4};       Plane Surface(1) = {1};
> Line Loop(2) = {5,6,7,8};       Plane Surface(2) = {2};
> Line Loop(3) = {-3,9,-5,10};    Plane Surface(3) = {3};
> Line Loop(4) = {11,-7,12,-1};   Plane Surface(4) = {4};
> Line Loop(5) = {-10,-8,-11,-4}; Plane Surface(5) = {5};
> Merge "z1.msh"; // mesh file of face 6
> 
> Surface Loop(1) = {-1,-2,-3,-4,-5,-6};
> Volume(1) = {1};
> 
> 
> ------------------------------------------------------------------------
> 
> Point(5) = {0, 0, 1, 0.1};
> Point(6) = {0, 1, 1, 0.1};
> Point(7) = {1, 1, 1, 0.1};
> Point(8) = {1, 0, 1, 0.1};
> Line(2) = {6,5} ;
> Line(6) = {8,7} ;
> Line(9) = {5,8} ;
> Line(12) = {7,6} ;
> 
> Line Loop(6) = {-2,-12,-6,-9};
> Plane Surface(6) = {6};
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> 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