3 cubit.cmd('brick x 10 y 10 z 10')
4 cubit.cmd('volume all size 10')
5 cubit.cmd('volume all scheme tetmesh')
6 cubit.cmd('mesh volume all')
7 cubit.cmd('set duplicate block elements off')
8 cubit.cmd('block 1 volume 1')
9 cubit.cmd('block 1 attribute count 2')
10 cubit.cmd('block 1 name "MAT_ELASTIC1"')
11 cubit.cmd('block 1 attribute index 1 1')
12 cubit.cmd('block 1 attribute index 2 0')
13 cubit.cmd('create displacement on surface 4 dof 1 fix 1')
14 cubit.cmd('create displacement on surface 3 dof 2 fix 1')
15 cubit.cmd('create displacement on surface 2 dof 3 fix 1')
16 cubit.cmd('create displacement on surface 1 dof 3 fix -1')
17 cubit.cmd('create displacement on surface 5 dof 2 fix -1')
18 cubit.cmd('create displacement on surface 6 dof 1 fix -1')
19 cubit.cmd('save as "/Users/michaelcortis/Documents/PhD/moFEM/mofem/meshes/tri_axial_cube_atom.cub" ')
