v0.14.0
dome_thin_symm.jou
Go to the documentation of this file.
1 reset
2 set duplicate block elements on
3 
4 create sphere radius 1.02
5 create sphere radius 0.98
6 
7 subtract volume 2 from volume 1
8 
9 brick x 4 y 4 z 4
10 move Volume 4 x 0 y -2 z 0 include_merged
11 subtract volume 4 from volume 3
12 
13 brick x 4 y 4 z 4
14 move Volume 5 x -2 y 0 z 0 include_merged
15 subtract volume 5 from volume 3
16 
17 brick x 4 y 4 z 4
18 move Volume 6 x 0 y 0 z -2 include_merged
19 subtract volume 6 from volume 3
20 
21 create Cylinder height 4 radius 0.2
22 #brick x 0.4 y 0.4 z 4
23 
24 #rotate Volume 5 angle 90 about X include_merged
25 ##move Volume 5 x 0.2 y 0 z 0 include_merged
26 
27 chop volume 3 with volume 7 keep
28 
29 delete volume 3
30 delete volume 7
31 
32 merge volume all with volume all
33 
34 #modify curve 74 69 chamfer radius 0.02
35 
36 #create Displacement on curve 79 dof 1 dof 2 dof 3 fix 0
37 create Displacement on surface 45 dof 3 fix 0
38 create Displacement on surface 41 47 dof 1 fix 0
39 create Displacement on surface 39 44 dof 2 fix 0
40 create pressure on surface 40 magnitude 1
41 #create pressure on surface 46 magnitude 1
42 
43 
44 
45 nodeset 4 vertex 40
46 nodeset 4 name 'LoadPath'
47 
48 volume all scheme Tetmesh
49 volume all size auto factor 8
50 mesh volume all
51 
52 block 3 tet all
53 block 3 element type tetra10
54 
55 block 1 tet all
56 block 1 name "MAT_ELASTIC_1"
57 block 1 attribute count 10
58 block 1 attribute index 1 1
59 block 1 attribute index 2 0.25
60 block 1 attribute index 3 0
61 block 1 attribute index 4 0
62 block 1 attribute index 5 0
63 block 1 attribute index 6 0
64 block 1 attribute index 7 0
65 block 1 attribute index 8 0
66 block 1 attribute index 9 0
67 block 1 attribute index 10 0
68 
69 block 2 tet all
70 block 2 name "BODY_FORCES")
71 block 2 attribute count 4
72 block 2 attribute index 1 1.
73 block 2 attribute index 2 0.
74 block 2 attribute index 3 0.
75 block 2 attribute index 4 0.
76 
77 save as "/Users/likask/MyBuild/mofem-cephas/mofem/meshes/dome_thin_symm.cub" overwrite
78 
79 
80 
81