[Gmsh] RE : 3D Mesh problem

Yann COLLETTE ycollette at free.fr
Fri Dec 8 08:11:42 CET 2006


Hello,

Thanks a lot. It works fine now.

YC

Patrick Lefevre a écrit :
> Hello Collette,
>  
> I think there is a mistake in the physical volume decleration. Indeed, with the extrude command with surface, the volumes are automatically created. You can see that in using the very useful visibility option (ctrl+shift+v). The number of the 2 volumes created are 1 and 2. So I have modified your geo file and it seems to work now. 
>  
> Regards,
> Patrick
>
> ________________________________
>
> De: gmsh-bounces at geuz.org de la part de Collette Yann
> Date: jeu. 30/11/2006 18:45
> À: gmsh at geuz.org
> Objet : [Gmsh] 3D Mesh problem
>
>
>
> Hello,
>
> I've got a problem with a 3d volume: when I produce my volume and mesh it
> under gmsh, no problem.
> But when I use gmsh -3 myfile.geo I've got only meshing of surfaces, no 3d
> meshes ...
> I've tested various combination of physical volume with no success.
> What kind of mistake did I made ?
>
> Yann COLLETTE
>
>
>   
> ------------------------------------------------------------------------
>
> lo = 1;
> la = 1;
> ep = 0.5;
> rho = 0.2;
> int = 0.2;
> ext = 0.3;
> mesh = 0.1;
> Point(1) = {0,0,0,mesh};
> Point(2) = {int,0,0,mesh};
> Point(3) = {-la/2+rho,lo/2,0,mesh};
> Point(4) = {la/2-rho,lo/2,0,mesh};
> Point(5) = {la/2-rho,-lo/2,0,mesh};
> Point(6) = {-la/2+rho,-lo/2,0,mesh};
> Point(7) = {-la/2,-lo/2+rho,0,mesh};
> Point(8) = {-la/2,lo/2-rho,0,mesh};
> Point(9) = {la/2,lo/2-rho,0,mesh};
> Point(10) = {la/2,-lo/2+rho,0,mesh};
> Line(2) = {4,3};
> Line(3) = {8,7};
> Line(4) = {6,5};
> Line(5) = {10,9};
> Point(11) = {la/2-rho,lo/2-rho,0,mesh};
> Point(12) = {-la/2+rho,lo/2-rho,0,mesh};
> Point(13) = {-la/2+rho,-lo/2+rho,0,mesh};
> Point(14) = {la/2-rho,-lo/2+rho,0,mesh};
> Circle(6) = {9,11,4};
> Circle(7) = {3,12,8};
> Circle(8) = {7,13,6};
> Circle(9) = {5,14,10};
> Point(15) = {-int,0,0,mesh};
> Point(16) = {0,int,0,mesh};
> Point(17) = {0,-int,0,mesh};
> Circle(10) = {2,1,16};
> Circle(11) = {16,1,15};
> Circle(12) = {15,1,17};
> Circle(13) = {17,1,2};
> Line Loop(14) = {6,2,7,3,8,4,9,5};
> Line Loop(15) = {10,11,12,13};
> Point(18) = {ext,0,0,mesh};
> Point(19) = {-ext,0,0,mesh};
> Point(21) = {0,ext,0,mesh};
> Point(22) = {0,-ext,0,mesh};
> Circle(16) = {18,1,21};
> Circle(17) = {21,1,19};
> Circle(18) = {19,1,22};
> Circle(19) = {22,1,18};
> Line Loop(20) = {16,17,18,19};
> Plane Surface(21) = {14,20};
> Plane Surface(22) = {20,15};
> Extrude {0,0,ep} {
>   Surface{21,22};
> }
>
> Physical Surface(127) = {84};
> Physical Surface(128) = {21};
> Physical Surface(129) = {126};
> Physical Surface(130) = {22};
> Physical Surface(131) = {47,51,55,59,63,67,39,43};
> Physical Surface(132) = {83,71,75,79};
> Physical Surface(133) = {121,117,113,125};
> Physical Line(134) = {25,24,31,30,29,28,27,26};
> Physical Line(135) = {7,3,8,4,9,5,6,2};
> Physical Line(136) = {35,34,33,32};
> Physical Line(137) = {19,16,17,18};
> Physical Line(138) = {92,91,90,93};
> Physical Line(139) = {10,11,12,13};
> /*Surface Loop(140) = {84,39,-21,43,47,51,55,59,63,67,71,75,79,83};
> Volume(141) = {140};
> Surface Loop(142) = {126,-83,-22,-79,-75,-71,113,117,121,125};
> Volume(143) = {142};*/
> Physical Volume(144) = {1, 2};
>
> Coherence;
>