title model of rear hub
para 
  r 7.0                c average radius of sinusoidal curve
  s 1                  c scaling factor for the mesh
  t 1                  c minimum number of nodes per region
r0  [%r+0.5*cos(  0)]  c radii for the sinusoidal curve
r1  [%r+0.5*cos( 10)]  c radii for the sinusoidal curve
r2  [%r+0.5*cos( 20)]  c radii for the sinusoidal curve
r3  [%r+0.5*cos( 30)]  c radii for the sinusoidal curve
r4  [%r+0.5*cos( 40)]  c radii for the sinusoidal curve
r5  [%r+0.5*cos( 50)]  c radii for the sinusoidal curve
r6  [%r+0.5*cos( 60)]  c radii for the sinusoidal curve
r7  [%r+0.5*cos( 70)]  c radii for the sinusoidal curve
r8  [%r+0.5*cos( 80)]  c radii for the sinusoidal curve
r9  [%r+0.5*cos( 90)]  c radii for the sinusoidal curve
r10 [%r+0.5*cos(100)]  c radii for the sinusoidal curve
r11 [%r+0.5*cos(110)]  c radii for the sinusoidal curve
r12 [%r+0.5*cos(120)]  c radii for the sinusoidal curve
r13 [%r+0.5*cos(130)]  c radii for the sinusoidal curve
r14 [%r+0.5*cos(140)]  c radii for the sinusoidal curve
r15 [%r+0.5*cos(150)]  c radii for the sinusoidal curve
r16 [%r+0.5*cos(160)]  c radii for the sinusoidal curve
r17 [%r+0.5*cos(170)]  c radii for the sinusoidal curve
r18 [%r+0.5*cos(180)]; c radii for the sinusoidal curve
ld 1 lp 
[%r0 *sin(  0)] [%r0 *cos(  0)]
[%r1 *sin( -1)] [%r1 *cos( -1)]
[%r2 *sin( -2)] [%r2 *cos( -2)]
[%r3 *sin( -3)] [%r3 *cos( -3)]
[%r4 *sin( -4)] [%r4 *cos( -4)]
[%r5 *sin( -5)] [%r5 *cos( -5)]
[%r6 *sin( -6)] [%r6 *cos( -6)]
[%r7 *sin( -7)] [%r7 *cos( -7)]
[%r8 *sin( -8)] [%r8 *cos( -8)]
[%r9 *sin( -9)] [%r9 *cos( -9)]
[%r10*sin(-10)] [%r10*cos(-10)]
[%r11*sin(-11)] [%r11*cos(-11)]
[%r12*sin(-12)] [%r12*cos(-12)]
[%r13*sin(-13)] [%r13*cos(-13)]
[%r14*sin(-14)] [%r14*cos(-14)]
[%r15*sin(-15)] [%r15*cos(-15)]
[%r16*sin(-16)] [%r16*cos(-16)]
[%r17*sin(-17)] [%r17*cos(-17)]
[%r18*sin(-18)] [%r18*cos(-18)] ;
sd 1 cy 0 0 0 0 0 1 6.625
sd 2 cy 0 0 0 0 0 1 3.23
sd 3 cy 0.0 5.56 0.0 0 0 1 0.507
sd 4 cy 0.0 3.5 0.0 0 0 1 0.2735
sd 5 cy 0 3.5 0 0 0 1 0.6349
sd 6 cy -1.15 3.6 0 0 0 1 0.562
sd 7 ts 0 0 5.81 0 0 1 1.2 0 3.14
sd 8 cy 0 0 0 0 0 1 2.375
sd 9 ts 0 0 5.81 0 0 1 1.59 0 2.0
sd 10 ts 0 0 5.81 0 0 1 2.15 0 1.44
sd 11 cp 1 rx -90 ;
sd 12 plan 0 0 0 1 0 0
sd 13 plan 0 0 0 [cos(18)] [sin(18)] 0
sd 14 plan 0 0 6.44 0 0 1
sd 15 plan 0 0 7.94 0 0 1
sd 16 cy 0 0 0 0 0 1 4.34
gct 1 ryz ; lev 1 grep 0 1 ; ; pslv 1
para
i2 [  1+max(%t,2*%s)]
i3 [%i2+max(%t,  %s)]
i4 [%i3+max(%t,2*%s)]
i5 [%i4+max(%t,3*%s)]
i6 [%i5+max(%t,3*%s)]
i7 [%i6+max(%t,4*%s)]
i8 [%i7+max(%t,5*%s)]
j2 [  1+max(%t,2*%s)]
j3 [%j2+max(%t,6*%s)]
k2 [  1+max(%t,2*%s)]
k3 [%k2+max(%t,2*%s)]
k4 [%k3+max(%t,3*%s)]
k5 [%k4+max(%t,3*%s)]
k6 [%k5+max(%t,3*%s)] ;
cylinder
1 %i2 %i3 %i4 %i5 %i6 %i7 %i8 ;
1 %j2 %j3 ;
1 %k2 %k3 %k4 %k5 %k6 ;
2.59 3.096 3.257 3.8873 4.34 5.053 6.067 6.625
90.0 94.5 108.0
5.81 6.44 7.25 7.94 8.9 10.0
pb 6 2 2 7 2 4 y 96.0
pa 4 0 6 x 3.5383
pa 4 0 3 x 3.4
mb 7 2 1 7 2 6 xy .6 .6
mb 4 1 3 4 3 3 yz .6 -.6
de 5 0 1 8 0 2
de 1 0 1 4 0 3
de 1 0 4 2 0 6
de 2 0 5 3 0 6
de 4 0 4 8 0 6
de 6 1 2 7 2 4
sfi -6 -7 ; 1 -2 ; 2 4 ;sd 3
sfi 1 -4 ; ; 1 -3 ; sd 10
c
c choose one of the following two commands 
c This determines the shape of the outer disc
c
c sfi -8;1 3;1 6;sd 1  c cylinder
  sfi -8;1 3;1 6;sd 11 c sinusoidal curve
sfi  ;-1; ;sd 12
sfi  ;-3; ;sd 13
sfi 4 8; ;-4;sd 15      
sfi 5 8; ;-2;sd 14
sfi -5; ;1 2;sd 16
res 8 1 2 8 3 4 j 1
res 7 1 2 8 1 4 k 1
res 5 3 2 8 3 2 i 1
res 4 3 4 8 3 4 i 1
res 1 1 3 4 3 3 i 1
res 4 1 1 4 3 3 k 1
res 5 2 2 8 3 4 k 1
res 5 1 2 6 1 4 k 1
relaxi 1 8;-3;1 6;5 0.001 1
relaxi 1 5;-1;1 4;5 0.001 1
relaxi 4 8; ;-4;5 0.001 1
relaxi 5 8; ;-2;5 0.001 1
lct 10 rz 0 ; rz 36 ; rz 72 ; rz 108 ; rz 144 ;
ryz rz 36 ; ryz rz 72 ; ryz rz 108 ; ryz rz 144 ; ryz rz 180 ;
lrep 1 2 3 4 5 6 7 8 9 10 ;
endpart
para i2 [  1+max(%t, 2*%s)]
     i3 [%i2+max(%t, 2*%s)]
     i4 [%i3+max(%t,   %s)]
     i5 [%i4+max(%t, 2*%s)]
     j2 [  1+max(%t, 4*%s)]
     j3 [%j2+max(%t, 4*%s)]
     j4 [%j3+max(%t, 2*%s)]
     k2 [  1+max(%t,10*%s)]
     k3 [%k2+max(%t, 2*%s)]
     k4 [%k3+max(%t, 3*%s)]
     k5 [%k4+max(%t, 7*%s)] ;
cylinder
1 %i2 %i3 %i4 %i5 ;
1 %j2 %j3 %j4;
1 %k2 %k3 %k4 %k5 ;
2.375 2.865 3.2265 3.7734 4.1349
90.0 96.75 103.5 112.5
0. 2.75 3.444 3.949 5.81
de 1 0 1 2 0 3
de 1 0 4 3 0 5
de 3 1 1 4 2 2
ma 5 2 4 y 2.5
ma 5 3 4 y 1.0
mb 3 2 1 3 2 2 xy .1 5
mb 4 2 1 4 2 2 xy -.1 5
sf 5 1 1 5 2 4 sd 5
sf 5 2 1 5 3 3 sd 6
sf 5 3 1 5 4 3 sd 2
sfi -3 -4 ; 1 -2 ; 1 2 ;sd 4
sf 1 0 3 1 0 4 sd 8
sfi 1 -3 ; ; -4 5 ; sd 9
sf 5 3 3 5 4 5 sd 7
sf 5 0 5 5 0 5 sd 7
lin 3 4 4 5 4 4
lin 2 4 1 5 4 1
lin 2 4 1 5 4 4
lin 2 2 1 5 4 4
lin 2 3 1 5 3 1
lin 2 2 1 5 3 1
lin 2 3 1 5 4 1
res 3 0 4 5 0 5 i 1.0
res 2 0 3 5 0 3 i 1.0
res 2 3 1 5 4 2 i 1.0
res 3 1 5 5 4 5 j 1
res 3 1 4 5 4 4 j 1
res 1 1 3 1 4 4 j 1
res 2 1 1 2 4 4 j 1
relaxi ;-1;2 5;5 0.001 1
relaxi ;-4; ;5 0.001 1
lct 8 rz 0 ; rz 45 ; rz 90 ; rz 135 ;
ryz rz 45 ; ryz rz 90 ; ryz rz 135 ; ryz rz 180 ;
lrep 1 2 3 4 5 6 7 8 ;
endpart
pplv
merge
rx 30 ry 45 zf 1.7 tvv
