CalculiX-Examples/Linear/PSE/par.run.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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
# <b0=1> width of the plane models
# <b1=10> width of the beam
# <b2=150> width of the plate
# <he=10> heigth
# <le=50> length
# <dist=30> distance of the models

seto PE
pnt ! 0 0 0
swep all new tra 0 <b0> 0 1 a
setc

seto PS
pnt ! 0 <b0+dist> 0
swep PS new tra 0 <b0> 0 1 a
setc

seto Beam
pnt ! 0 <2*b0+2*dist> 0
swep Beam new tra 0 <b1> 0 2 a
setc

seto Plate
pnt ! 0 <2*b0+b1+3*dist> 0
swep Plate new tra 0 <b2> 0 20 a
setc

swep all new tra 0 0 <he> 2 a
swep all new tra <le> 0 0 20 a

elty all he8i
mesh all
send all abq

# mesh plot
frame
rot y
rot l 20
rot u 25
zoom 2.3
view elem
ulin From bottom right: plane strain, plane stress, beam, plate
seta ! all
hcpy png mesh

# z support
seta nodes n all
enq nodes fixz rec <le> _ _ 0.1
send fixz abq spc 3
ulin Constraint uz=0 (support)
plot n fixz
hcpy png uz0

# y support
enq nodes fixy rec _ <b0/2.> _ <1.1*b0/2.>
seta fixy p D003 D005 D007
comp fixy do
send fixy abq spc 2
ulin Constraint uy=0 (rigid body constraint and plane strain)
plot n fixy
hcpy png uy0

# symmetry
enq nodes symx rec 0 _ _ 0.1
send symx abq spc 1
ulin Constraint ux=0 (symmetry)
plot n symx
hcpy png ux0

# load
seta load s A00D A00E A00F A00G
comp load do
comp load do
send load abq pres 1
text Pressure application area
plot f load
hcpy png load

sys ccx solve
read solve.frd

# sys wmctrl -r :ACTIVE: -e 0,200,0,1200,600

# sys wmctrl -r Calculix -e 0,200,0,1200,600
# rot y
# rot l 60
# rot u 30
# zoom 2


view disp
view elem off
view sh
scal d 200
rot y
rot l 10
zoom
zoom 3
ds 1 e 3
text Vertical displacement
plot fv all
hcpy png U3-flat

seta base se all
copy all new mir x
rot y
rot l 20
rot u 25
frame
zoom 2.3
tra l <le/2>

ds 1 e 3
plot fv base
hcpy png U3

text Transverse strain
ds 3 e 2
hcpy png EYY

text Transverse stress
ds 2 e 2
hcpy png SYY

text Longitudinal (bending) stress
ds 2 e 1
hcpy png SXX

text Equivalent (von Mises) stress
ds 2 e 7
hcpy png SE