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