CalculiX-Examples/Drahtbiegen/Zug/par.pre.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
# run param.py pre.fbd.par
#
# <Dia=1.>         Diameter
# <DivR=10>      Radial division
# <Len=2>          Length
# <DivL=50>      Longitudinal division
# <BiasL=4>      Longitudinal bias
# <Etyp='qu8cr'>  Element type

# geometry
pnt p1 0 0 0
swep all rad tra <Dia/2.> 0 0 <DivR>
seta bottom se all
swep all top tra 0 <Len> 0 <DivL>

# meshing
enq all vertlines rect _ <Len/2.> 0 0.1
comp vertlines up
bia vertlines <BiasL>
flip all
elty all <Etyp>
mesh all
send all abq

# symmetry axis
seta nodes n all
enq nodes sym rec  0 _ 0 0.001
comp sym down

send all abq nam

# display of sets
rot -z
frame
tra l <Dia/2.>
plot la all
plus sa all
plus pa all
hcpy png zug-geo

ulin Sets: top (blue), sym (green), bottom (turq)
view elem
plot f all n
plus n top b 8
plus n bottom t 8
plus n sym g 6


hcpy png zug-sets
# END OF SCRIPT