Hi all,
I try to create a body of revolution with gmsh.
The shape will be quite complex. But for tests I have tried to create a
simple sphere by extrusion of 2*Pi. I have attached my example.
The problem is that I am not able to create the volume mesh.
Can somebody give me a hint, what's going wrong?
Thanks a lot,
Patrick
System: Ubuntu 11.04 64bit
Version : 2.5.0
GUI toolkit : FLTK 1.1.10
License : GNU General Public License
Build OS : Linux
Build options : Ann Bamg Blas Chaco DIntegration Dlopen FlTree Fltk Gmm
Have64BitSizeT Jpeg Kbipack Lapack MPI MathEx Mesh Metis Mpeg
NativeFileChooser Netgen OpenCascade OpenGL Parser Plugins Png Post
Solver Zlib
Build date : 20101030
Build host : crested
Packager : buildd
Web site : http://www.geuz.org/gmsh/
Mailing list : [email protected]
// Center of sphere 01
centerSphere01X = 0.0;
centerSphere01Y = 0.0;
centerSphere01Z = 0.0;
// Dimension of sphere 01
radiusSphere01AlongAxis = 0.1;
radiusSphere01PerpAxis = 0.1;
// Characteristic length for all.
cl = 0.02;
// Points of body
ipt20 = newp; Point(ipt20) = {centerSphere01X, centerSphere01Y,
centerSphere01Z, cl};
ipt21 = newp; Point(ipt21) = {centerSphere01X, centerSphere01Y, centerSphere01Z
- radiusSphere01AlongAxis, cl};
ipt22 = newp; Point(ipt22) = {centerSphere01X, centerSphere01Y +
radiusSphere01PerpAxis, centerSphere01Z, cl};
ipt23 = newp; Point(ipt23) = {centerSphere01X, centerSphere01Y, centerSphere01Z
+ radiusSphere01AlongAxis, cl};
ipt24 = newp; Point(ipt24) = {centerSphere01X, centerSphere01Y +
0.2*radiusSphere01PerpAxis, centerSphere01Z, cl};
// Create lines for ellipsoid hull cell 01.
icl20 = newreg; Ellipse(icl20) = {ipt21,ipt20,ipt24,ipt22};
icl21 = newreg; Ellipse(icl21) = {ipt22,ipt20,ipt24,ipt23};
out1[] = Extrude{{0,0,1},{0,0,0}, Pi}{ Line{icl20}; };
out2[] = Extrude{{0,0,1},{0,0,0},-Pi+0.0001}{ Line{icl20}; };
out3[] = Extrude{{0,0,1},{0,0,0}, Pi}{ Line{icl21}; };
out4[] = Extrude{{0,0,1},{0,0,0},-Pi+0.0001}{ Line{icl21}; };
// Make surface loop for volume of cubes.
ilsl1 = newreg; Surface Loop(ilsl1) = {out1[1],out2[1],out3[1],out4[1]};
// Define volume of cubes.
ivl1 = newv; Volume(ivl1) = {ilsl1};
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
