!   TITLE:  retainer_od_race.mac (AAAV)
!
!-----------------------------------------------------------------------
!       Rev -   10/31/2002      D. Edwards
!-----------------------------------------------------------------------
!
!   Analysis of upper azimuth bearing retainer to determine proper clearance
!   and thickness to achieve desired preload on outer DT bearing races.
!
!   Assembly: CU01715097 prototype machined and sleeved to use DT bearings.
!
!   CS_0:       Cartesian ring Coordinate System
!   CS_1:       Cylindrical ring Coordinate System
!   CS_2:       Spherical ring Coordinate System
!
!
!************************** TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Solid Model
!
!********************* Model Assembly Preparation **********************
!
mod_assemb_prep,arg1,arg2,arg3      ! prepare database for assembly
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>
!
retainer_od_race.parm
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate System >>>>>>>>>>>>>>>>>>>>>
!
mod_coord_system,acs,xc,yc,zc,rotxyc,rotyzc,rotzxc
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 3.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!
csys,CS_0

k,,r1
k,,r2
k,,r2,t1
k,,ro,t1
k,,ri,t1
k,,r3,t1
k,,r3,t1-t2
k,,r1,t1-t2

create_area

view,'f'


! Apply pressure load along annulus about bolt pattern
sfl,4,pres,Pt

!MSH_AREA, mstyp, mshdens, etpa, elmnl, elmxl, lszmsh, mat, real, ecs
msh_area,0,0,1,.01,.01,,c416

! Apply contact elements at bearing race contact
lsel,s,line,,7
nsll,s,0                    ! 0=interior nodes only
ngen,2,nm+1,all,,,,-.01     ! Copy nodes down y=-.01

type,2
mat,2
real,2

eintf,,,low,0,0,-.01,0

nsll,u,1
d,all,uy

init_entity_num

! Apply contact elements at sleeve surface

lsel,s,line,,1
nsll,s,1                    ! All nodes
ngen,2,nm+1,all,,,,-gapS    ! Copy nodes down the correct gap distance

type,3
mat,2
real,2

eintf,,,low,0,0,-gapS,0

nsll,u,1
d,all,uy

init_entity_num,1

fini

/sol

antype,0

solcontrol
!autots,on
!nsub,20
solv

fini

/post1

/dscale,,1

avprin,0,0,
plnsol,s,eqv,2,1

lsel,s,line,,1
nsll

!fcn1='gapS+'
rplot_nodal,0,'x','U','Y'