[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