[Gmsh] 3D mesh

Philippe Ackerer ackerer at unistra.fr
Mon Oct 7 22:53:39 CEST 2013


Hi
I missunderstood something. I plan to discretize a colmun.
Could you explain me why only plan z=0 (nodes and elements) is stored in 
the .msh file ?
I want also to identify the top surfaces, as I have done for the bottom 
surfaces (physical surface 101 and 'Dirichlet'). Any tips are welcome.
I attached the file.
Thanks and thanks a lot for this very helpfull mesh generator.

Best regards
Philippe

-- 
-------------------------------------------------------------
Philippe ACKERER
LHyGeS UMR 7517 CNRS-UdS Engees
1, rue Blessig
67000 STRASBOURG
Tel: +33 (0)388 560 561

-------------- next part --------------
/********************************************************************* 
 *
 *  Column discretization
 *
 *********************************************************************/

// As usual, we start by defining some variables:

r1=0.01;
r2=0.10;
Lcw=0.001;
Lc1 =0.005;
teta=3.14159/3.0;

// Then we define some points and some lines using these variables:

Point(1) = { 0.0, 0.0  , 0.0 , Lcw};
Point(2) = { r1, 0.0   , 0.0 , Lcw};
Point(3) = { r2 , 0.0   , 0.0 , Lc1};
x1=r1*Cos(teta);
y1=r1*Sin(teta);
Point(4) = { x1 , y1 , 0.0 , Lcw};
x1=r2*Cos(teta);
y1=r2*Sin(teta);
Point(5) = { x1 , y1 , 0.0 , Lc1};


// circle (n1,n2,n3), circle starting at point n1 and ending at point n3, with centre n2

Circle(1) = {2,1,4};
Circle(2) = {3,1,5};

Line(10) = {1,2};
Line(11) = {2,3};
Line(12) = {1,4};
Line(13) = {4,5};

Line Loop(21) = {10,1,-12};     Plane Surface(31) = {21};
Line Loop(22) = {-1,11,2,-13};  Plane Surface(32) = {22};

Physical Surface(101)={32};
Physical Surface("Dirichlet") = {31} ;

//Extrusion de la surface sur une distance h
  h = 0.200;
Extrude {0,0,h} {
Surface{31,32}; Layers{ {20}, {1} };Recombine;}