/usr/share/doc/rheolef-doc/examples/torus.mshcad is in rheolef-doc 6.5-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | R = 1; // large radius
r = 0.6; // small radius
n = 10; // density of discretisation
h = 1.0/n;
n1 = Ceil(2*Pi*r/h);
n2 = Ceil(2*Pi*R/(3*h));
Printf("n = %g", n);
Printf("n1 = %g", n1);
Printf("n2 = %g", n2);
Point(1) = { R, 0, 0, h};
Point(2) = {R+r, 0, 0, h};
Point(3) = { R, 0, r, h};
Point(4) = { R, 0,-r, h};
Point(5) = {R-r, 0, 0, h};
Circle(1) = {2,1,3};
Circle(2) = {3,1,5};
Circle(3) = {5,1,4};
Circle(4) = {4,1,2};
Line Loop(5) = {1,2,3,4};
Plane Surface(6) = {5};
third1[]=Extrude {{0,0,1}, {0,0,0}, 2*Pi/3}{Surface{6}; Layers{n2};};
third2[]=Extrude {{0,0,1}, {0,0,0}, 2*Pi/3}{Surface{third1[0]}; Layers{n2};};
third3[]=Extrude {{0,0,1}, {0,0,0},-2*Pi/3}{Surface{6}; Layers{n2};};
Physical Surface("boundary") = {
-third1[2],-third1[3],-third1[4],-third1[5],
-third2[2],-third2[3],-third2[4],-third2[5],
third3[3], third3[4], third3[5], third3[6]
};
|