CalculiX-Examples/Contact/Tubes/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
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
# parameters
valu R1 50
valu T1 5
valu L1 500
valu R2 45
valu T2 5
valu L2 500
valu OVL 150
valu divt 1
valu divrot 18
valu divl 50

# dependent values
valu T1neg * T1 -1
valu T2neg * T2 -1
valu x21 - L1 OVL
valu x22 + x21 L2

# sketch of tube 1
seto tube1
pnt p1 0 0 R1
swep all new tra 0 0 T1neg divt
swep all new tra L1 0 0 divl
swep all new rot x 180 divrot
setc
# sketch of tube 2
seto tube2
pnt p2 x21 0 R2
swep tube2 new tra 0 0 T2neg divt
swep tube2 new tra L2 0 0 divl
swep tube2 new rot x 180 divrot
setc

# meshing
node 1 x22 0 0

elty all he8i
mesh all
send all abq

# node sets for picking
seta nodes n all
seta n1 n tube1
seta n2 n tube2

# fix
enq n1 fix rec 0 _ _ 0.1
send fix abq nam

# symmetry
enq nodes ysym rec _ 0 _ 0.1
send ysym abq nam

# load
enq n2 load rec x22 _ _ 0.1
send abq load
comp load do
send load abq surf

# indep
valu R1i - R1 T1
enq n1 ind cx R1i _ _ 0.1
comp ind do
send ind abq surf

# dep
enq n2 dep cx R2 _ _ 0.1
comp dep do
send dep abq surf


rot -y
plot f all n
frame
rot l 80
rot u 10
plus f ind b
plus f dep r
view elem
zoom 4
hcpy png contact

plot n fix m 6
plus n load t 6
plus n ysym r 4
ulin fix (magenta), ysym (red), load (turq)
hcpy png sets