5 create Cylinder height 0.8 major radius 1 minor radius 1
6 move Volume 2 location surface 4 include_merged
7 webcut volume 2 with plane yplane imprint merge
9 chop volume 1 with volume 2
10 chop volume 5 with volume 3
16 block 1 name 'MAT_ELASTIC'
17 block 1 attribute count 2
18 block 1 attribute index 1 10.0
19 block 1 attribute index 2 0.3
21 set duplicate block elements on
24 block 2 name "BODY_FORCES")
25 block 2 attribute count 4
26 block 2 attribute index 1 1.
27 block 2 attribute index 2 0.
28 block 2 attribute index 3 0.
29 block 2 attribute index 4 0.
32 sideset 100 curve 32 29 33 48 30 28 49 47 27 1 5 31 remove
33 nodeset 101 vertex all
34 nodeset 101 vertex 17 20 18 19 29 30 remove
35 Sideset 102 surface all
36 sideset 102 surface 29 18 19 6 remove
37 Sideset 10200 surface 6
39 Sideset 200 surface 19
42 create Displacement name '' on surface 3 dof 1 dof 2 dof 3 fix 0
43 #create Displacement name '' on surface 6 dof 3 fix 0
44 create pressure name '' on surface 5 magnitude 1
46 volume all scheme Tetmesh
47 volume all size auto factor 9
52 ##refine curve 27 29 numsplit 5 bias 1.0 depth 1 smoot
53 #save as "/Users/likask/MyBuild/mofem-bitbucket/mofem_v0.2/meshes/crack1.cub" overwrite