/com  File: stub2.mac
/com  Revised Stub Shaft with single thread retainers for AAAV - P/N: 
/com  Apply a pressure at the axial CG to give 100 lbf lateral load.
/com  Mounting surface nodes fixed in all directions.

fini
/clear,all

/title, Revised Stub Shaft - 100 lb Lateral Loading

/prep7

/pnum,KP,1
/pnum,LINE,1

F = 100         ! [lbf] Applied load to resolver cover
nApress=4       ! Area # to apply pressure to
Kaxis1=21       ! KP's that define axis of revolution
Kaxis2=22

K,1,0.8561700,0.5330250,0.
K,2,0.2500000,0.7903281,0.
K,3,1.530900,0.6114967,0.
K,4,1.530900,0.9502500E-01,0.
K,5,1.160000,0.,0.
K,6,1.160000,0.4710242,0.
K,7,0.4998000,2.787025,0.
K,8,0.4998000,1.956024,0.
K,9,0.6715020,1.036024,0.
K,10,0.6715013,1.004027,0.
K,11,0.7965011,0.8790241,0.
K,12,1.148223,0.8790235,0.
K,13,0.5979250,1.036024,0.
K,14,0.5098000,1.946024,0.
K,15,0.5879250,1.946024,0.
K,16,0.2500000,3.037025,0.
K,17,0.4690000,3.037025,0.
K,18,1.098000,0.5330242,0.
K,19,1.350000,0.,0.
K,20,1.630000,0.9502500E-01,0.
K,21,0.,3.037025,0.
K,22,0.,0.,0.
K,23,0.2500000,2.111000,0.
K,24,0.7965013,1.004024,0.
K,25,1.236612,0.9156350,0.
K,26,1.400000,1.079023,0.
K,27,1.148224,1.004024,0.
K,28,1.459900,1.079023,0.
K,29,1.459900,0.6638573,0.
K,30,1.479900,0.6114967,0.
K,31,1.479900,0.6414967,0.
K,32,1.630000,0.,0.
K,33,0.5879250,1.046024,0.
K,34,0.2500000,2.611000,0.
K,35,0.5979250,1.046024,0.
K,36,0.5098000,1.956024,0.
K,37,1.098000,0.4710242,0.
L,1,2
L,3,4
L,5,6
L,7,8
L,9,10
L,11,12
L,13,9
L,14,15
L,16,17
L,18,1
L,19,5
L,4,20
!L,21,22            ! Centerline
L,23,2
L,25,26
L,26,28
L,28,29
L,3,30
L,20,32
L,15,33
L,32,19
L,16,34
L,34,23
L,7,17

larc,10,11,24,.125
larc,12,25,27,.125
larc,29,30,31,0.30000E-01
larc,33,13,35,0.10000E-01
larc,8,14,36,0.10000E-01
larc,6,18,37,0.62000E-01

*go,:skip

et,1,shell93
et,2,solid95

mp,EX,1,28.5e6
mp,NUXY,1,.28

al,all

desize,,,,,,,.04

!lesize,9,,,7,1,

amesh,all

type,2
esize,,9
vrotat,all,,,,,,Kaxis1,Kaxis2,180,2

aclear,1

numcmp,all

asel,s,loc,z,0
da,all,symm

asel,s,loc,y,0
da,all,all,0


/com Apply 100 equivalent pounds to the part in the lateral direction
! Calculate the pressure required to generate F/2 pounds on the model half.
! Apress = (pi/2)*(rL)
! Aproj = rL

asel,s,area,,nApress
asum,fine
*get,Apress,area,0,area
Aproj = Apress*(2/3.1415)
press = (F/2)/Aproj
/PSF,PRES,NORM,2,0
/PBF,DEFA, ,1
sfa,all,1,PRES,press

alls

!*go,:skip

static

fini
/post1
! The CG is mid way between the applied pressure length. "uxcg" is thus the
! deflection at the CG of the Azimuth gimbal assy (24 lb) located at y = 2.361 inches
uxcg = (UX(ikp+11) + UX(ikp+29))/2
klat=F/uxcg
*msg,info,F,uxcg,klat
Stub Shaft Lateral Stiffness:  k = F/uxcg = %g / %g = %g

:skip