v0.14.0
elasticty_atom_test_thermal_expansion.jou
Go to the documentation of this file.
1 #!python
2 cubit.cmd('reset')
3 
4 cubit.cmd('brick x 1 y 1 z 1')
5 cubit.cmd('set duplicate block elements on')
6 cubit.cmd('block 1 volume 1 ')
7 cubit.cmd('block 1 attribute count 3')
8 
9 cubit.cmd('block 1 name "MAT_ELASTIC1"')
10 #!Young's Modulus
11 cubit.cmd('block 1 attribute index 1 1.')
12 #!Poisson ratio
13 cubit.cmd('block 1 attribute index 2 0.25')
14 #!Ternak expansion:
15 cubit.cmd('block 2 attribute index 3 1.0')
16 
17 cubit.cmd('block 2 volume 1 ')
18 cubit.cmd('block 2 attribute count 2')
19 
20 cubit.cmd('block 2 name "MAT_THERMAL"')
21 cubit.cmd('block 2 attribute count 2')
22 #!Thermal expansion
23 cubit.cmd('block 2 attribute index 1 1.')
24 #!Heat capacity
25 cubit.cmd('block 2 attribute index 2 1.')
26 
27 #!Displacement bcs
28 cubit.cmd('create Displacement on vertex 7 dof 1 dof 2 dof 3 fix 0')
29 cubit.cmd('create Displacement on vertex 4 dof 1 dof 2 fix 0')
30 cubit.cmd('create Displacement on vertex 8 dof 1 fix 0')
31 
32 #!Mesh
33 cubit.cmd('volume 1 size 1')
34 cubit.cmd('volume 1 scheme Tetmesh')
35 cubit.cmd('mesh volume 1')
36 
37 cubit.cmd('save as "/Users/likask/MyBuild/mofem-bitbucket/mofem/meshes/elasticity_atom_test_thermal_expansion.cub" overwrite')
38 
39