CalculiX-Examples/Linear/BiaxBending/run.fbl at master · rahul-e/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
# th=10.
# he=50.
# lx=200.


# rotated rectangle section
pnt ! 0 5.0 25.0
swep all new tra 0 -10.0 0 2
swep all new tra 0 0 -50.0 8
move all rot x 45
swep all load1 tra 200.0 0 0 10

# z section
# y0z=-45
# z0z=20
# la=40
# lb=10
# th=3
del se new
seto zprof
pnt p0z 0 -45 20
swep zprof new tra 0 -3 0 2
swep zprof new1 tra 0 0 -3 2
swep new1 new2 tra 0 0 -34 6
swep new2 new3 tra 0 0 -3 2
seta fu l L00F
swep fu new tra 0 10 0 3
seta fo l L00M
swep fo new tra 0 -10 0 3
swep zprof load2 tra 200.0 0 0 10
setc

# angle section
# y0a=-90
# z0a=20
# la=40
# lb=20
# tha=3
del se new
seto angle
pnt p0a 0 -90 20
swep angle new tra 0 -3 0 2
swep new new1 tra 0 -17 0 4
swep angle new tra 0 0 -3 2
seta tmp l L01N
swep tmp new tra 0 0 -37 8
swep angle load3 tra 200.0 0 0 10
setc

# ref nodes
node 1 200.0 0 0
node 2 200.0 -45 20
node 3 200.0 -90 20

seta ref n all

# meshing
elty all he8i
mesh all

# quick and dirty support
seta nodes n all
enq nodes support rec 0 _ _ 0.01

# load application surfaces
send load1 abq surf
send load2 abq surf
send load3 abq surf

# write mesh and sets
send all abq
send all abq nam

# plot of mesh and ref node
ulin support (red), load (yellow) with ref node (blue)
frame
zoom 1.5
rot -x
rot c 90
rot l 10
rot u 10
view elem
plot n ref r
plus f all n
plus f load1 y
plus f load2 y
plus f load3 y
plus n ref b 10
plus n support r 10
hcpy png Refs/mesh

# Solving
sys ccx Biegung
read Biegung.frd new

# display of stress distribution and deformed geometry
ulin Stress distribution
view disp
view sh
scal d 1000
# sxx
ds 2 e 1
hcpy png Refs/sxx

# neutral fiber
ulin Neutral fiber
steps 2
mm 1
hcpy png Refs/neutral

#
ulin displacement
rot -x
rot c 90
plot n all r 1
plus f all n
hcpy png Refs/disp