-
Notifications
You must be signed in to change notification settings - Fork 0
/
circ.geo
104 lines (102 loc) · 5 KB
/
circ.geo
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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
SetFactory("OpenCASCADE");
Point(1) = {1.0},0.0,0};
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{1};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{2};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{3};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{4};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{5};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{6};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{7};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{8};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{9};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{10};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{11};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{12};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{13};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{14};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{15};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{16};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{17};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{18};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{19};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{20};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{21};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{22};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{23};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{24};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{25};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{26};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{27};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{28};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{29};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{30};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{31};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{32};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{33};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{34};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{35};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{36};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{37};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{38};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{39};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{40};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{41};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{42};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{43};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{44};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{45};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{46};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{47};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{48};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{49};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{50};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{51};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{52};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{53};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{54};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{55};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{56};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{57};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{58};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{59};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{60};}
Extrude {0, 0, 1}, {0, 0, 0}, -Pi/{Nel} {Point{61};}
Extrude {0, 0, 1}, {0, 0, 0}, -{gaplengths[ii // 2]} {Point{62};}
Extrude {0, 0, 1}, {0, 0, 0}, -{0.95 * gaplengths[-1]} {Point{63};}
Line(65) = {65,1};
Curve Loop(1) = {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, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65};
Plane Surface(1) = {1};
Physical Curve(1) = {1};
Physical Curve(2) = {3};
Physical Curve(3) = {5};
Physical Curve(4) = {7};
Physical Curve(5) = {9};
Physical Curve(6) = {11};
Physical Curve(7) = {13};
Physical Curve(8) = {15};
Physical Curve(9) = {17};
Physical Curve(10) = {19};
Physical Curve(11) = {21};
Physical Curve(12) = {23};
Physical Curve(13) = {25};
Physical Curve(14) = {27};
Physical Curve(15) = {29};
Physical Curve(16) = {31};
Physical Curve(17) = {33};
Physical Curve(18) = {35};
Physical Curve(19) = {37};
Physical Curve(20) = {39};
Physical Curve(21) = {41};
Physical Curve(22) = {43};
Physical Curve(23) = {45};
Physical Curve(24) = {47};
Physical Curve(25) = {49};
Physical Curve(26) = {51};
Physical Curve(27) = {53};
Physical Curve(28) = {55};
Physical Curve(29) = {57};
Physical Curve(30) = {59};
Physical Curve(31) = {61};
Physical Curve(32) = {63};
Physical Surface(50) = {1};
Mesh.SecondOrderLinear = 1;