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