Skip to content
Snippets Groups Projects
Commit e280b539 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

simpler

parent 29acc340
Branches
Tags
No related merge requests found
...@@ -5,54 +5,52 @@ Mesh.CharacteristicLengthMin = 0.1; ...@@ -5,54 +5,52 @@ Mesh.CharacteristicLengthMin = 0.1;
Mesh.CharacteristicLengthMax = 0.1; Mesh.CharacteristicLengthMax = 0.1;
// build volume through (closed) line loops // build volume through (closed) line loops
Disk(1) = {0,0,0, 0.2}; Circle(1) = {0,0,0, 0.2};
Disk(2) = {0.1,0.04,1, 0.3}; Circle(2) = {0.1,0.04,1, 0.3};
Disk(3) = {0.03,-0.08,2, 0.25}; Circle(3) = {0.03,-0.08,2, 0.25};
For i In {1:3} For i In {1:3}
Line Loop(i) = Boundary{ Surface{i}; }; Line Loop(i) = i;
Delete { Surface{i}; }
EndFor EndFor
ThruSections(1) = {1:3}; ThruSections(1) = {1:3};
Disk(4) = {1,0,0, 0.2}; Circle(4) = {1,0,0, 0.2};
Disk(5) = {1.1,0.04,1, 0.3}; Circle(5) = {1.1,0.04,1, 0.3};
Disk(6) = {1.03,-0.08,2, 0.25}; Circle(6) = {1.03,-0.08,2, 0.25};
For i In {4:6} For i In {4:6}
Line Loop(i) = Boundary{ Surface{i}; }; Line Loop(i) = i;
Delete { Surface{i}; }
EndFor EndFor
Ruled ThruSections(2) = {4:6}; Ruled ThruSections(2) = {4:6};
// build surfaces through (closed) line loops // build surfaces through (closed) line loops
Disk(7) = {2,0,0, 0.2}; Circle(7) = {2,0,0, 0.2};
Disk(8) = {2.1,0.04,1, 0.3}; Circle(8) = {2.1,0.04,1, 0.3};
Disk(9) = {2.03,-0.08,2, 0.25}; Circle(9) = {2.03,-0.08,2, 0.25};
For i In {7:9} For i In {7:9}
Line Loop(i) = Boundary{ Surface{i}; }; Line Loop(i) = i;
Delete { Surface{i}; }
EndFor EndFor
a() = ThruSections{7:9}; ThruSections{7:9}
Disk(10) = {3,0,0, 0.2}; Circle(10) = {3,0,0, 0.2};
Disk(11) = {3.1,0.04,1, 0.3}; Circle(11) = {3.1,0.04,1, 0.3};
Disk(12) = {3.03,-0.08,2, 0.25}; Circle(12) = {3.03,-0.08,2, 0.25};
For i In {10:12} For i In {10:12}
Line Loop(i) = Boundary{ Surface{i}; }; Line Loop(i) = i;
Delete { Surface{i}; }
EndFor EndFor
b() = Ruled ThruSections{10:12}; // builds ruled surfaces Ruled ThruSections{10:12}
// build surfaces through wires // build surfaces through (open) wires
n = 100; Circle(13) = {4,0,0, 0.2, Pi/3};
Point(n+1) = {0,-2,0}; Circle(14) = {4.1,0.04,1, 0.3, Pi/2};
Point(n+2) = {0.5,-1.5,0}; Circle(15) = {4.03,-0.08,2, 0.25, Pi/3};
Point(n+3) = {0,-1.5,0}; For i In {13:15}
Circle(n+1) = {n+1,n+3,n+2}; Wire(i) = i;
Point(n+4) = {0,-2,1}; EndFor
Point(n+5) = {0.7,-1.3,1}; ThruSections{13:15}
Point(n+6) = {0,-1.3,1};
Circle(n+2) = {n+4,n+6,n+5};
Wire(n+1) = {n+1}; Circle(16) = {5,0,0, 0.2, Pi/3};
Wire(n+2) = {n+2}; Circle(17) = {5.1,0.04,1, 0.3, Pi/2};
c() = ThruSections{n+1,n+2}; Circle(18) = {5.03,-0.08,2, 0.25, Pi/3};
For i In {16:18}
Wire(i) = i;
EndFor
Ruled ThruSections{16:18}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment