-
Notifications
You must be signed in to change notification settings - Fork 101
/
Copy pathpar.pre.fbl
134 lines (120 loc) · 2.34 KB
/
par.pre.fbl
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
132
133
134
# test of geometry generation for ccx
# <xh=100>
# <rh=5>
# <lx=xh+2*rh>
# <wi=4*rh>
# <th=rh>
# <angle=-60>
# <dz=10>
# outer rectangle
pnt p1 0 <-wi/2> 0
pnt p2 <lx> <-wi/2> 0
pnt p3 <lx> <wi/2> 0
pnt p4 0 <wi/2> 0
line r1 p1 p2
line r2 p2 p3
line r3 p3 p4
line r4 p4 p1
# circle
pnt pc <xh> 0 0
pnt pc1 <xh-rh> 0 0
pnt pc2 <xh> <-rh> 0
pnt pc3 <xh+rh> 0 0
pnt pc4 <xh> <rh> 0
line c1 pc1 pc2 pc
line c2 pc2 pc3 pc
line c3 pc3 pc4 pc
line c4 pc4 pc1 pc
# create part
surf s1 all
swep all new tra 0 0 <th>
# duplicate part
seto part2
copy all new1 rot pc y 180
setc
move part2 rot pc z <angle>
move part2 tra 0 0 <-dz>
# meshing
div all auto <rh/2>
elty all te10
mesh all
seta parts e all
# display
view elem
plot f all
# pin
seto ref1
pnt ref1 <xh> 0 <th/2.>
setc ref1
seto ref2
pnt ref2 <xh> 0 <-th/2.-dz>
setc ref2
seto pin1
line ! ref1 ref2 2
setc pin1
seto pin2
swep ref1 pin2 tra 0 0 <th/2.> 2
swep ref2 pin2 tra 0 0 <-th/2.> 2
setc
seta pin se pin1 pin2
elty pin be3
mesh pin
stack on
prnt se ref1
valu nref1 pop 18
prnt se ref2
valu nref2 pop 18
# coupling surfaces
seta s1 s A006 A007 A008 A009
seta s2 s A00G A00H A00I A00J
comp s* do
comp s* do
send s1 abq sur
send s2 abq sur
# write coupling cards, kinematic
sys echo *coupling, ref node= nref1 ,surface= Ss1 ,constraint name=p11 > kpin.inc
sys echo *kinematic >> kpin.inc
sys echo 1,3 >> kpin.inc
sys echo *coupling, ref node= nref2 ,surface= Ss2 ,constraint name=p12 >> kpin.inc
sys echo *kinematic >> kpin.inc
sys echo 1,3 >> kpin.inc
# support
seta nodes n all
enq nodes fix rec 0 _ _ 0.1
# load
seta load s A00F
comp load do
comp load do
send load abq sur
# write mesh and sets
send all abq
send all abq nam
# plot of the model
frame
rot -z
rot u 90
rot l 45
rot u 45
tra l <0.3*xh>
plot f all w
hcpy png Refs/mesh
# plot of support and load application sets
ulin nodeset 'fix' and surface 'load'
plot n fix r
plus f load b
plus l all
hcpy png Refs/fix_load
# plot of coupling details
ulin Coupling surfaces and reference nodes
frame pin
zoom 0.5
#rot z
view ill
plot l all r 2
plus l part2 b 2
plus f s1 r
plus pa ref1 r 10
plus f s2 b
plus pa ref2 b 10
plus e pin k 6
hcpy png Refs/coupling