v0.14.0
Strip_footing.jou
Go to the documentation of this file.
1 #!python
2 cubit.cmd('new')
3 
4 SigY=848.7e-3;
5 sqrt3=1.7321;
6 c=SigY/sqrt3
7 Pres=5.14*c;
8 
9 
10 cubit.cmd('set duplicate block elements on ')
11 cubit.cmd('brick x 5 y 5 z 0.1')
12 
13 cubit.cmd('brick x 0.5 y 5 z 0.1')
14 cubit.cmd('move surface 10 location surface 4 include_merged')
15 cubit.cmd('chop volume 1 with volume 2')
16 
17 cubit.cmd('brick x 1.5 y 1.5 z 0.1')
18 cubit.cmd('move volume 5 x -1.75 y 1.75 include_merged')
19 cubit.cmd('chop volume 3 with volume 5 keep')
20 cubit.cmd('chop volume 4 with volume 5 keep')
21 cubit.cmd('delete volume 3 4 5')
22 
23 cubit.cmd('imprint volume all')
24 cubit.cmd('merge volume all')
25 
26 str1='create pressure on surface 31 magnitude '+str(Pres); cubit.cmd(str1)
27 
28 cubit.cmd('create Displacement on surface 33 39 dof 1 fix 0')
29 cubit.cmd('create Displacement on surface 35 48 dof 2 fix 0')
30 cubit.cmd('create Displacement on surface 53 45 29 36 49 38 43 32 dof 3 fix 0')
31 
32 cubit.cmd('volume 6 8 size auto factor 7')
33 cubit.cmd('volume 6 8 scheme tetmesh')
34 cubit.cmd('mesh volume 6 8')
35 
36 cubit.cmd('volume 7 9 size auto factor 8')
37 cubit.cmd('volume 7 9 scheme tetmesh')
38 cubit.cmd('mesh volume 7 9')
39 
40 
41 cubit.cmd('block 1 vertex 38')
42 cubit.cmd('block 1 name "LoadPath"')
43 
44 cubit.cmd('save as "/Users/zahur/Desktop/Plasticity/Strip_footing.cub" overwrite')
45 
46