mesh.geo 838 Bytes
Newer Older
dubois's avatar
maj    
dubois committed
1
d=1e-3;
dubois's avatar
sablier  
dubois committed
2
3
L = 15*d;
H = 50*d;
Jonathan Lambrechts's avatar
Jonathan Lambrechts committed
4
5
lc = 3*d;
R=2*d;
dubois's avatar
sablier  
dubois committed
6
7
8
9

Point(1) = {0., 0., 0, lc};
Point(2) = {L, 0., 0, lc};
Point(3) = {L, H, 0, lc};
Jonathan Lambrechts's avatar
Jonathan Lambrechts committed
10
11
Point(4) = {0.5*(L+R), H, 0, lc/3};
Point(5) = {0.5*(L+R), H+d, 0, lc/3};
dubois's avatar
sablier  
dubois committed
12
13
14
15
Point(6) = {L, H+d, 0, lc};
Point(7) = {L, H+d+H, 0, lc};
Point(8) = {0., H+d+H, 0, lc};
Point(9) = {0., H+d, 0, lc};
Jonathan Lambrechts's avatar
Jonathan Lambrechts committed
16
17
Point(10) = {0.5*(L-R), H+d, 0, lc/3};
Point(11) = {0.5*(L-R), H, 0, lc/3};
dubois's avatar
sablier  
dubois committed
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
Point(12) = {0., H, 0, lc};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 9};
Line(9) = {9, 10};
Line(10) = {10,11};
Line(11) = {11,12};
Line(12) = {12,1};

Line Loop(1) = {1:12};

Plane Surface(1) = {1};
Physical Line("Box") = {1,2,3,4,5,6,8,9,10,11,12};
Physical Line("Top") = {7};
Physical Surface("Domain") = {1};
Physical Point("PtFix") = {1};