!   TITLE:  tc_stiff.mac
!
!---------------------------------------------------------------------------
!       Rev -   03/25/2001      D. Edwards
!---------------------------------------------------------------------------
!
!   Resume the db file saved by tc.mac.  Apply constraints and loads
!   so that a load verses deflection solve can be used to compute
!   structure stiffness.
!
!   CS10:           Cartesian Coordinate System
!   CS11:           Cylindrical Coordinate System
!   CS12:           Spherical Coordinate System
!
!   USAGE:  tc_stiff
!
!******************************* TABLE OF CONTENTS **************************
!
!       Model Assembly Preparation
!
!   1.  Parameter List
!   2.  Model Coordinate System
!   3.  Material and Element Type
!   4.  Solid Model
!       4.1  Outer perimeter
!       4.2  Surface Cut Patterns
!       4.3  Inner Cylinder
!       4.4  Create Separate Areas with Coincident Entities
!       4.5  Create Section Symmetry
!   5.  FE Model
!       5.1  Mesh Areas
!       5.2  Apply Constraints & Coupling
!           5.2.1  Apply Constraints
!           5.2.2  Apply Coupling
!   6.  Solve
!   7.  Post Processing
!
!***************************** END Table of Contents ************************
!
resume,tc,db

fini

/prep7

! Parameter list
uy0=.005            ! Arbitrary deflection for measuring stiffness

tstart

/title,Stiffness of Lower Azimuth Bearing w/ Epoxied Sleeve
!
!


! Delete lower BC's and reapply with displacement uy0
csys,CS10
cmsel,s,nbottom
ddele,all,uy
d,all,uy,uy0
alls
ddele,ncorner,ux
d,ncorner,ux,-uy0       ! ux is vertical due to node cs being rotated

! Select the inner set of surface nodes of the outer race and radially
! constrain them.
csys,CS11
nsel,s,loc,x,rracei
nsel,r,loc,y,-90,0

nrotat,all
d,all,ux

alls

fini

/solu
antype,0
tunif,Troom
tref,Troom
solv
fini

/post1

csys,CS11
nsel,s,loc,x,rracei

nplot

/PSF,DEFA, ,1,0
/PBF,DEFA, ,1
/PSYMB,DOT,1
/PSYMB,PCONV,
/PBC,ALL, ,1
/REP

fsum

tstop
:skip