!   TITLE:  ring.mac (AAAV)
!
!---------------------------------------------------------------------------
!       Rev -   12/11/2000  D. Edwards
!---------------------------------------------------------------------------
!
!   Axisymmetric model of the outer race of the lower azimuth axis
!   bearing 490A045-102.  NHBB Outer Ring, P/N SSMER-1907DU20A7.
!
!   CS10:   Cartesian 1 Coordinate System
!   CS11:   Cylindrical 1 Coordinate System
!   CS12:   Spherical 1 Coordinate System
!
!
!******************************* TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Material and Element Type
!   4.  Solid Model
!
!
!************************** Model Assembly Preparation **********************
!
!
mod_assemb_prep
tstart
!
!
/title,NHBB Outer Race, SSMER-1907DU20A7 (490A045-102)
!
!
/units,bin
!
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>
!
!
ring.parm                   ! name of parameter file
!
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate Systems >>>>>>>>>>>>>>>>>>>>
!
!
mod_coord_system,0,0,0,0,0,0,0
!
!
!<<<<<<<<<<<<<<<<<<<< 3.  Material and Element Type >>>>>>>>>>>>>>>>>>>
!
!
!---------------------- Material Properties --------------------
!

alum=1
cres=2

mp,ex,alum,10.4E6
mp,dens,alum,.097/386.4
mp,nuxy,alum,0.33
mp,alpx,alum,aalum

mp,ex,cres,28.0E6
mp,dens,cres,.283/386.4
mp,nuxy,cres,0.27
mp,alpx,cres,acres
!
!------------------------- Element Type ------------------------
!
!et,1,plane82,,,1       ! Axisymmetric elements
et,1,plane42,,,1        ! Axisymmetric elements

!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 4.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!

!
! *** Create and mesh the outer race
!

!*get,km,kp,0,num,max
init_entity_num

k,km+1,1.0141,0
k,km+2,1.01985,.16841364

k,km+3,1.0141,.236855
kcontact=km+3

k,km+4,.97,.2965488

l,km+1,km+2
larc,km+2,km+4,km+3

ldiv,lm+2,0.0,km+3,2,0      ! Divide the track at the ball contact point

k,km+5,orL,orW-orr1
k,km+6,orL+orr1,orW
k,km+7,orD-orrc,orW
k,km+8,orD,orW-orrc
k,km+9,orD,orCH
k,km+10,orD-orCH,0

*do,i,km+4,km+9
    l,i,i+1
*enddo

l,km+10,km+1

al,all

dk,km+1,uy      ! Constrain a point on the outer race

/view,1,,,1
/ang,1
lplot


!                  (1)      (2)     (3)     (4)      (5)       (6)     (7)   (8)
! MSH_ALGORITHM, perflsz, perfmsh, mstyp, mshdens, elmszovr, atargov, etpa, etpv,
!
!                elmnl, elmxl, mtrl, rl, ecs, volsel
!                 (9)    (10)  (11) (12) (13)  (14)

msh_algorithm,1,1,0,0,,,1,,.01,.03,cres,CS10

!
!*** Create and mesh the aluminum ring
!

init_entity_num

k,,orD,0
k,,orD,orCH
k,,orD,orW-orrc
k,,orD,orW
k,,orD+alT,orW
k,,orD+alT,0

create_area

init_entity_num

asel,s,area,,am
msh_algorithm,1,1,0,0,,,1,,.01,.03,alum,CS10

eplot

!
!----------------------- Apply Boundary conditions
!

! Couple the interface nodes in x direction only
nsel,s,loc,x,orD
cpintf,ux
alls

dk,km,uy            ! Tie down the aluminum sleeve

!----------------------- 6  Solve -----------------------
!

fini
/solu

antype,0
tunif,Ttest
tref,Troom

solv

fini

/post1

ksel,s,kp,,kcontact
nslk,s
*get,ncontact,node,0,num,max
*get,uxct,node,ncontact,u,x         ! Total movement of node

! Compensate for thermal shrinkage of race by itself (no alum ring)
uxc=uxct-nx(ncontact)*(Ttest-Troom)*acres

/out,ring,sum
*msg
Net Radial deflection of ball contact point with outer race.

*msg

*msg
For use in press.tkw to determine equivalent ring thickness:

*msg,,uxc
dri =  %g

/out
*uilist,ring,sum

tstop

:skip