CalculiX-Examples/Linear/Crack1/ct.fbd at master · declmal/CalculiX-Examples · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
sys gmsh ct.geo
read gmsh.inp

# remove all surface elements
zap +CPS6

send all abq
send Symy abq nam
send Symz abq nam
comp Load do
send Load abq pres 1
# point at origin for x constraint
enq all origin rec 0 0 0 0.1
send origin abq nam

# plot of the sets
rot -z
rot l 30
rot d 30
view ill
plot fb all n 30
plus f Load
plus n Symy r
plus n Symz b
plus n origin g 8
hcpy png Refs/sets

# solver run
sys ccx solve
read solve.frd

# symmetry expansion
seta orig e all
copy all new mir y
copy all new mir z
# plot of the equivalent stress
view disp
view sh
scal d 2000
ds -1 e 7
rot -z
rot r 25
rot u 5
plot ev orig
plus e new n
min 0 f
max 10 f
frame
hcpy png Refs/se

# display of the strain energy
sys more solve.dat