!
!   TITLE:  ms_sol.mac
!
!   Read in Ansys DB file "mirror_support.db", apply boundary conditions
!   and perform a static analysis.
!
!   The applied boundary conditions are for determining the torsional
!   stiffness about the axis of rotation.
!
!   D. Edwards 10/16/00
!

!*go,:skip

FINISH

! /psearch,e:\sl\ansys\ansys\aaav\elev_axis\

/clear,start
!resume,mesh_msupport,db
resume,ms_mesh,db
!
/out,ms_sol_modal,out

/prep7

tstart

/pbc,all, ,0        ! Turn off symbols until after nrotat command

! Constrain the torquer bearing journal in all directions
csys,CSm_trq1

parm_val,'mir_supt','EXTD5','extd5'
parm_val,'mir_supt','L18','l18'
parm_val,'mir_supt','D1','d1'
asel,s,loc,z,1.05*extd5,.95*(l18+extd5)
asel,r,loc,x,0,d1/2
da,all,all,0,1
alls


! Constrain band drive journal to rotate only
csys,CSm_res1

parm_val,'mir_supt','EXTD4','extd4'
parm_val,'mir_supt','TSIDE','tside'
parm_val,'mir_supt','DOS','bRadius',2.000
bRadius=bRadius/2       ! [in]    Radius of band drive journal

tside=.045

nsel,s,loc,z,tside,tside+extd4
nsel,r,loc,x,bRadius

nrot,all
d,all,ux
alls

*get,ncnt,node,0,count
*msg,info,ncnt,ncnt*3
There are %G nodes and %G degrees of freedom in this model

/PBC,ALL, ,1
aplot

FINISH
/config,fsplit,250

/solu

!******************************
!*** Powerdynamics - PCG Solve
!******************************
antype,2

modopt,subsp,4
eqslv,pcg,1.0e-5
mxpand,4, , ,0
lumpm,1
pstres,0

modopt,subsp,4,10,1000, ,off
rigid,
subopt,0,0, ,100,

solve

tstop

finish
/post1

rsys,CSm_res1

/title,Mirror Support - Modal Frequency

/out