[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;}