# vim:fileencoding=utf-8:ft=cgx
# All values in SI units: m,s,kg,N,Pa,J
# (length,time,mass,force,stress,energy).

rot y
rot u 30
rot r 20
zoom 0.8

valu divp 8

pnt RP1 0 0.2071 0.0101
pnt RP2 0 0.2249 0.0101
pnt RP3 0 0.2249 0.0279
pnt RP4 0 0.2071 0.0279
pnt CP  0 0.2160 0.0190

seta pincontour p RP1 RP2 RP3 RP4
seto pincontour
pnt PP1 0 0.2071 0.0051
pnt PP2 0 0.2249 0.0051
pnt PP3 0 0.2299 0.0101
pnt PP4 0 0.2299 0.0279
pnt PP5 0 0.2249 0.0328
pnt PP6 0 0.2071 0.0328
pnt PP7 0 0.2021 0.0279
pnt PP8 0 0.2021 0.0101
line LP01 PP1 PP2 divp
line LP02 PP2 PP3 RP2 divp
line LP03 PP3 PP4 divp
line LP04 PP4 PP5 RP3 divp
line LP05 PP5 PP6 divp
line LP06 PP6 PP7 RP4 divp
line LP07 PP7 PP8 divp
line LP08 PP8 PP1 RP1 divp
setc pincontour
comp pincontour do

line LP09 RP1 RP2 divp
line LP10 RP2 RP3 divp
line LP11 RP3 RP4 divp
line LP12 RP4 RP1 divp
line LP13 RP1 PP8 divp
line LP14 RP1 PP1 divp
line LP15 RP2 PP2 divp
line LP16 RP2 PP3 divp
line LP17 RP3 PP4 divp
line LP18 RP3 PP5 divp
line LP19 RP4 PP6 divp
line LP20 RP4 PP7 divp

seto pin
surf SP1 blend LP09 LP10 LP11 LP12
surf SP2 blend LP01 LP15 LP09 LP14
surf SP3 blend LP15 LP02 LP16
surf SP4 blend LP16 LP03 LP17 LP10
surf SP5 blend LP17 LP04 LP18
surf SP6 blend LP18 LP05 LP19 LP11
surf SP7 blend LP19 LP06 LP20
surf SP8 blend LP20 LP07 LP13 LP12
surf SP9 blend LP13 LP08 LP14
setc pin
comp pin do

valu dl 64
valu lp 0.16

seto pin
swep pin IP1 tra lp 0 0 dl
setc pin
comp pin do

seta pintop s A00D A00E A00F
comp pintop do
seta pinbot s A00A A00B A00H
comp pinbot do

# solid shell
seto shell
pnt PP1a 0 0.2071 0.0050
swep shell new tra 0 0 -0.002 2
swep shell J01 tra 0 0.0178 0 divp
swep J01 J02 rot RP2 x 90 divp
swep J02 J03 tra 0 0 0.0178 divp
swep J03 J04 rot RP3 x 90 divp
swep J04 J05 tra 0 -0.0178 0 divp
swep J05 J06 rot RP4 x 90 divp
swep J06 J07 tra 0 0 -0.0178 divp
swep J07 new rot RP1 x 90 divp
setc shell
merg p shell
merg l shell
seto shell
swep shell new tra 0.08 0 0 32
setc shell
seta shelltop s A01L A01O A01R
comp shelltop do
seta shellbot s A01C A01F A01W
comp shellbot do

# Reserve nodes for e.g. reference nodes for rigid bodies.
asgn n 2
# 0.08 + 0.448 = 0.528
node 1 0.528 0.216 0.019
elty all he20r
mesh pin
mesh shell

# Node sets for fixation and rigid body load.
seta nodes n pin
enq nodes load rec lp _ _ 1e-6

# Shell top surface
seta tsnodes n shelltop
enq tsnodes topindep rec 0.010000 _ _ 0.02 a
comp topindep f

# Shell bottom surface
seta bsnodes n shellbot
enq bsnodes botindep rec 0.070000 _ _ 0.02 a
comp botindep f

seta tpnodes n pintop
enq tpnodes topdep rec 0.010000 _ _ 0.02 a
comp topdep f

seta bpnodes n pinbot
enq bpnodes botdep rec 0.070000 _ _ 0.02 a
comp botdep f

valu output yes
# Write data
if output eq yes
send all abq
send pin abq nam
send shell abq nam
send load abq nam
send topindep abq sur
send topdep abq sur
send botindep abq sur
send botdep abq sur
endif