[Gmsh] error message

Erik Arlemark Erik.Arlemark.6091 at student.uu.se
Fri Jun 16 21:37:19 CEST 2006



Hello I got an error message when trying to do a 3D mesh of the following code.
The error message told me to send information of my error to this adress.

//radie för bränslestavar
r=7.25;

//tjocklek av finare mesh
f=1;

//distans mellan bränslestavar
d=3.825;

//antal cellager för finmesh
fl=3;

//fjärdedels båge (cellager)
fb=5;

//3/4 Pi båge (cellager)
tfb=15;

//antal cellager för block
bl=3;

//för extrude
//längd av figuren
L=10;

//antal cellager på längden
cl=4;

//volym 1

//mitt för volym 1
Point(1) = {0, 0, 0, 0.1};
//figur 1
Point(2) = {r, 0, 0, 0.1};
Point(3) = {r+f, 0, 0, 0.1};
Point(4) = {0, r+f, 0, 0.1};
Point(5) = {0, r, 0, 0.1};
Line(1) = {2,3};
Line(2) = {4,5};
Circle(3) = {5,1,2};
Circle(4) = {3,1,4};
Line Loop(5) = {4,2,3,1};
Plane Surface(1) = {5};
Transfinite Line {3,4} = fb Using Progression 1;
Transfinite Line {1,2} = fl Using Progression 1;
Transfinite Surface {1} = {2,3,4,5};
Recombine Surface {1};

Extrude {0,0,L} {
Surface{1}; Layers{ {cl}, {1}, {1} }; Recombine;
}

//volym 2

//centrum för högerkanten
Point(20) = {2*r+d, 0 , 0, 0.1};
//figur 2
Point(21) = {r+d-f, 0, 0, 0.1};
Point(22) = {2*r+d-(r+f)*Cos(Pi/4), (r+f)*Sin(Pi/4), 0, 0.1};
Point(23) = {(r+f)*Cos(Pi/4), (r+f)*Sin(Pi/4) , 0, 0.1};
Line(28) = {3,21};
Line(29) = {22,23};
Circle(30) = {21,20,22};
Circle(31) = {23,1,3};
Line Loop(32) = {30,29,31,28};
Plane Surface(2) = {32};
Transfinite Line {28,29,30,31} = bl Using Progression 1;
Transfinite Surface {2} = {3,21,22,23};
Recombine Surface {2};

Extrude {0,0,L} {
Surface{2}; Layers{ {cl}, {1}, {1} }; Recombine;
}

//volym 3

Point(31) = {r+d, 0, 0, 0.1};
Point(32) = {2*r+d+r*Cos(Pi/4), r*Sin(Pi/4), 0, 0.1};
Point(33) = {2*r+d+(r+f)*Cos(Pi/4), (r+f)*Sin(Pi/4), 0, 0.1};
Line(55) = {21,31};
Line(56) = {32,33};
Circle(57) = {31,20,32};
Circle(58) = {33,20,21};
Line Loop(58) = {57,56,58,55};
Plane Surface(3) = {58};
Transfinite Line {55,56} = fl Using Progression 1;
Transfinite Line {57,58} = tfb Using Progression 1;
Transfinite Surface {3} = {21,31,32,33};
Recombine Surface {3};

Extrude {0,0,L} {
Surface{3}; Layers{ {cl}, {1}, {1} }; Recombine;
}

//volym 4

Point(41) = {3*r+d, 0, 0, 0.1};
Point(42) = {3*r+d+f, 0, 0, 0.1};
Line(41) = {41,42}; 
Circle(81) = {32,20,41};
Circle(82) = {33,20,42};
Line Loop(83) = {81,41,-82,-56};
Plane Surface(4) = {83};
Transfinite Line {41} = fl Using Progression 1;
Transfinite Line {81,82} = fb Using Progression 1;
Transfinite Surface {4} = {41,42,33,32};
Recombine Surface {4};

Extrude {0,0,L} {
Surface{4}; Layers{ {cl}, {1}, {1} }; Recombine;
}

Thanks
/Erik