CalculiX-Examples/Linear/StressConc1/par.run1.fbl 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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
# <L=100>
# <W=60>
# <R=15>
# <M=W*W/6.>
# <M>

pnt p1 <L> 0 0
pnt p2 <R> 0 0
pnt p3 0 0 0
pnt p4 0 <R> 0
pnt p5 0 <W/2.> 0
pnt p6 <L> <W/2.> 0

# line name p1 p2 division growth-rate
line l1 p1 p2 20 0.9
line l2 p2 p4 p3 30 0.9
line l3 p4 p5 20 1
line l4 p5 p6 30 1.1
line l5 p6 p1 8

# plot of the line numbers
rot -z
plot la all
plus pa all
frame
hcpy png Refs/lines

# generate the surface from the edges
surf ! l1 l2 l3 l4 l5
flip all

# meshing
asgn n 2
elty all tr6u
mesh all
elty all tr6s

# double and merge the mesh
copy all new mir y
frame
merg n all
#merg p all # this leads to mesh error with CGX 2.10 under windows
#merg l all

# symmetry boundary
seta nodes n all
enq nodes sym rec 0 _ 0 0.0001 a

# load application surface
enq nodes load rec <L> _ 0 0.0001 a
comp load do
send load abq sur

# surface ref node and load
node 1 <L> 0 0
seta nref n 1
sys echo 1,6,<M> > load.inp

# write mesh and set definitions
send all abq
send all abq nam

# images
ulin Sets: sym (red), nref (big blue), load (blue)
plot f all n
plus n nref b 10
plus n sym r 6
plus f load b 6
frame
hcpy png Refs/mesh1

# solution
sys ccx Stress1

# post-processing
read Stress1.frd new
ulin Load applied via kinematic coupling
# vector plot of worst principal stress
view elem
view disp
scal d 10000
ds 2 e 23
scal v 3
plot fv all
frame
hcpy png Refs/vecplot1

# expanded deformation plot
# view disp
# scal d 10000
# seta base all
# copy base new mir x
# frame
# zoom 1.2
# plot fv base
# plus e new n
# plus fv base
 hcpy png Refs/D1
#
# debug plot
view elem
ds 3 e 4
plot fv all
hcpy png Refs/RF1