|
| 1 | +SetFactory("OpenCASCADE"); |
| 2 | + |
| 3 | +dx=0.01; |
| 4 | +w=1.0; |
| 5 | + |
| 6 | +n=81; |
| 7 | + |
| 8 | +Point(1)={0.0,0.0,0.0,dx}; |
| 9 | +Point(2)={ w,0.0,0.0,dx}; |
| 10 | +Point(3)={ w, w,0.0,dx}; |
| 11 | +Point(4)={0.0, w,0.0,dx};//+ |
| 12 | + |
| 13 | +//+ |
| 14 | +Line(1) = {1, 2}; |
| 15 | +//+ |
| 16 | +Line(2) = {2, 3}; |
| 17 | +//+ |
| 18 | +Line(3) = {3, 4}; |
| 19 | +//+ |
| 20 | +Line(4) = {4, 1}; |
| 21 | +//+ |
| 22 | +Line Loop(1) = {4, 1, 2, 3}; |
| 23 | +//+ |
| 24 | +Plane Surface(1) = {1}; |
| 25 | + |
| 26 | +Circle(5) = {0.25, 0.15, 0, 0.08, 0, 2*Pi}; |
| 27 | +//+ |
| 28 | +Circle(6) = {0.5, 0.15, 0, 0.08, 0, 2*Pi}; |
| 29 | +//+ |
| 30 | +Circle(7) = {0.75, 0.15, 0, 0.08, 0, 2*Pi}; |
| 31 | + |
| 32 | +//+ |
| 33 | +Circle(8) = {0.25, 0.4, 0, 0.08, 0, 2*Pi}; |
| 34 | +//+ |
| 35 | +Circle(9) = {0.5, 0.4, 0, 0.08, 0, 2*Pi}; |
| 36 | +//+ |
| 37 | +Circle(10) = {0.75, 0.4, 0, 0.08, 0, 2*Pi}; |
| 38 | + |
| 39 | +//+ |
| 40 | +Circle(11) = {0.75, 0.625, 0, 0.08, 0, 2*Pi}; |
| 41 | +//+ |
| 42 | +Circle(12) = {0.5, 0.625, 0, 0.08, 0, 2*Pi}; |
| 43 | +//+ |
| 44 | +Circle(13) = {0.25, 0.625, 0, 0.08, 0, 2*Pi}; |
| 45 | +//+ |
| 46 | + |
| 47 | +Circle(14) = {0.25, 0.85, 0, 0.08, 0, 2*Pi}; |
| 48 | +//+ |
| 49 | +Circle(15) = {0.5, 0.85, 0, 0.08, 0, 2*Pi}; |
| 50 | +//+ |
| 51 | +Circle(16) = {0.75, 0.85, 0, 0.08, 0, 2*Pi}; |
| 52 | +//+ |
| 53 | +Line Loop(2) = {14}; |
| 54 | +//+ |
| 55 | +Plane Surface(2) = {2}; |
| 56 | +//+ |
| 57 | +Line Loop(3) = {15}; |
| 58 | +//+ |
| 59 | +Plane Surface(3) = {3}; |
| 60 | +//+ |
| 61 | +Line Loop(4) = {16}; |
| 62 | +//+ |
| 63 | +Plane Surface(4) = {4}; |
| 64 | +//+ |
| 65 | +Line Loop(5) = {13}; |
| 66 | +//+ |
| 67 | +Plane Surface(5) = {5}; |
| 68 | +//+ |
| 69 | +Line Loop(6) = {12}; |
| 70 | +//+ |
| 71 | +Plane Surface(6) = {6}; |
| 72 | +//+ |
| 73 | +Line Loop(7) = {11}; |
| 74 | +//+ |
| 75 | +Plane Surface(7) = {7}; |
| 76 | +//+ |
| 77 | +Line Loop(8) = {8}; |
| 78 | +//+ |
| 79 | +Plane Surface(8) = {8}; |
| 80 | +//+ |
| 81 | +Line Loop(9) = {9}; |
| 82 | +//+ |
| 83 | +Plane Surface(9) = {9}; |
| 84 | +//+ |
| 85 | +Line Loop(10) = {10}; |
| 86 | +//+ |
| 87 | +Plane Surface(10) = {10}; |
| 88 | +//+ |
| 89 | +Line Loop(11) = {5}; |
| 90 | +//+ |
| 91 | +Plane Surface(11) = {11}; |
| 92 | +//+ |
| 93 | +Line Loop(12) = {6}; |
| 94 | +//+ |
| 95 | +Plane Surface(12) = {12}; |
| 96 | +//+ |
| 97 | +Line Loop(13) = {7}; |
| 98 | +//+ |
| 99 | +Plane Surface(13) = {13}; |
| 100 | +//+ |
| 101 | +BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Surface{3}; Surface{4}; Surface{5}; Surface{6}; Surface{7}; Surface{8}; Surface{9}; Surface{10}; Surface{11}; Surface{12}; Surface{13}; Delete; } |
| 102 | +//+ |
| 103 | +Physical Line("left") = {17}; |
| 104 | +//+ |
| 105 | +Physical Line("bottom") = {19}; |
| 106 | +//+ |
| 107 | +Physical Line("right") = {20}; |
| 108 | +//+ |
| 109 | +Physical Line("top") = {18}; |
| 110 | +//+ |
| 111 | +Physical Surface("block") = {1}; |
| 112 | +//+ |
| 113 | +Transfinite Line {14, 15, 16, 13, 12, 11, 8, 9, 10, 5, 6, 7} = n Using Progression 1; |
0 commit comments