/com  File: stub_demval.mac
/com  Stub Shaft used on M1A2 and Dem Val AAAV DAHA - P/N: A340A355
/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

/prep7

/title, M1A2/Dem Val AAAV DAHA Stub Shaft - 100 lb Lateral Loading

et,1,shell93
et,2,solid95

mp,EX,1,28.5e6      ! [psi] steel modulus of elasticity
mp,NUXY,1,.28       ! nu of steel

*get,ikp,kp,0,num,max
*get,ia,area,0,num,max

F = 100         ! [lbf] Applied load at CG of Az gimbal assy
nApress=ia+24       ! Area # to apply pressure to
Kaxis1=ikp+27       ! Kepoints that define axis of revolution
Kaxis2=ikp+26

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

K,ikp+1,0.3850000,0.7330242,0.
K,ikp+2,0.8561720,0.5330250,0.
K,ikp+3,1.530900,0.9502500E-01,0.
K,ikp+4,1.530900,0.6880234,0.
K,ikp+5,1.680000,0.,0.
K,ikp+6,1.680000,0.9502500E-01,0.
K,ikp+7,1.159998,0.4710245,0.
K,ikp+8,1.160000,0.,0.
K,ikp+9,0.3850000,1.205024,0.
K,ikp+10,0.3150020,1.245439,0.
K,ikp+11,0.3150009,2.111000,0.
K,ikp+12,0.4998000,1.796024,0.
K,ikp+13,0.4998000,2.342025,0.
K,ikp+14,0.6715000,0.7810242,0.
K,ikp+15,0.6715000,0.8760238,0.
K,ikp+16,1.350000,0.,0.
K,ikp+17,1.097998,0.5330242,0.
K,ikp+18,0.4689980,2.877025,0.
K,ikp+19,0.3150000,2.877024,0.
K,ikp+20,0.5882000,1.786024,0.
K,ikp+21,0.5098000,1.786024,0.
K,ikp+22,0.5982000,0.8760242,0.
K,ikp+23,1.332500,0.7190232,0.
K,ikp+24,0.7335000,0.7190242,0.
K,ikp+25,0.7335000,0.7810242,0.
K,ikp+26,0.,2.877022,0.
K,ikp+27,0.,0.,0.
K,ikp+28,0.5882000,0.8860242,0.
K,ikp+29,0.3150003,2.611000,0.
K,ikp+30,0.5982000,0.8860242,0.
K,ikp+31,0.5098000,1.796024,0.
K,ikp+32,1.097998,0.4710242,0.

L,ikp+1,ikp+2
L,3,4
L,5,6
L,7,8
L,9,1
L,10,9
L,11,10
L,12,13
L,14,15
L,6,3
L,16,5
L,2,17
L,18,19
L,20,21
L,15,22
L,23,24
!L,26,27            ! Centerline
L,20,28
L,19,29
L,29,11
L,8,16
L,23,4
L,18,13

LARC,14,24,25,.062
LARC,7,17,32,.062
LARC,28,22,30,.01
LARC,12,21,31,0.01

*go,:skip

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



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