fini
/clear

resume,rotor_housing_model_189,db
allsel
/eshape,1
eplo
finish
!
!##################
! Static analysis
!##################
!
/solu
antype,static
irlf,1	                        ! get accurate moment of inertia (0.147e06)
solve
fini
/post1
eplot
irlist

/eof

/solu
antype,static
irlf,0
coriolis,on,,,off		! static must be done in rotating reference frame
cmomega,rotshaft,200,0,0        ! 200 rads/sec on the rotor-shaft component 			
!!
cmsel,,rotshaft
nsle
ic,all,rotz,0.1	                ! initial velocities for wz = 0.1 rads/sec
!!
*get,n_count,node,,count
*do,iii,1,n_count,1
/gopr
ic,iii,uy,,0.1*nx(iii)
*enddo
!!
allsel
solve
finish
!*
/post1
set,last
prrsol                       	! Expected MY = 2 * 0.1 * 200 * 0.147e06 = 0.588e07 
				! ANSYS MY = 0.563e07 
fini
