/COM
/OUT,SCRATCH
/PREP7
*DIM,SPIN,,10 ! SPIN VELOCITY (RPM)
SPIN(1) = 1000.
SPIN(2) = 2000.
SPIN(3) = 5000.
SPIN(4) = 10000.
SPIN(5) = 15000.
SPIN(6) = 20000.
SPIN(7) = 25000.
SPIN(8) = 30000.
SPIN(9) = 40000.
SPIN(10) =50000.

RO = 7850 ! MATERIAL : STEEL
PEX = 2.059E+11
PGXY = 7.93E+10 ! NO SHEAR
MP,EX,1,PEX
MP,PRXY,1,0.3
MP,DENS,1,RO
MP,GXY,1,PGXY


ET,1,BEAM188 ! ELEMENT TYPE # 1: SHAFT
Konstante1 = 6 ! SHAFT SECTION PROPERTIES


*DIM,DIAM,ARRAY,Konstante1
DIAM(1) = 2.5E-2
DIAM(2) = 2.5E-2
DIAM(3) = 2.5E-2
DIAM(4) = 2.5E-2
DIAM(5) = 2.5E-2
DIAM(6) = 2.5E-2


*DO,I,1,Konstante1
SECTYPE,   I, BEAM, CSOLID, , 0 
SECOFFSET, CENT 
SECDATA,DIAM(I),1,1,0,0,0,0,0,0,0,0,0
*ENDDO


ET,2,MASS21 ! ELEMENT TYPE # 2: MASS
R,2,11087

ET,3,COMBI214 ! ELEMENT TYPE # 3: BEARING ELEMENT
KEYOPT,3,2,1 ! YZ PLANE
KEYOPT,3,3,0 ! ELEMENT IS SYMMETRIC


R,30,1E+8,1E+8,0,0,,    !-8.756E+6,-8.756E+6,, ! REAL CONSTANTS ( K11,K22,K12,K21) NO DAMPING
ET,4,COMBI214
KEYOPT,4,2,1 ! YZ PLANE
KEYOPT,4,3,0 ! ELEMENT IS SYMMETRIC
R,40,1E+8,1E+8,0,0,,    !,-8.756E+6,-8.756E+6,, ! REAL CONSTANTS ( K11,K22,K12,K21) NO DAMPING
/COM, NODES


N,1 ,0.
N,2 ,0.1
N,3 ,0.3
N,4 ,0.5
N,5 ,0.7
N,6 ,0.9
N,7 ,1.0

!BRG = 0.1 ! BEARING "LENGTH" FOR VISUALISATION
!N,6,0.1,BRG
!N,7,0.3,BRG

/COM, ELEMENTS
/COM, SHAFT
TYPE,1
MAT,1
*DO,I,1,Konstante1
REAL,I
E,I,I+1
*ENDDO



/COM, DISK
TYPE,2
REAL,2
E,4

!/COM, BEARINGS
!TYPE,3
!REAL,30
!E,2,6

!TYPE,4
!REAL,40
!E,4,7
!ALLSEL,ALL

! RB werden festgelegt
FINI
/SOLU
D,ALL,UX ! NO TRACTION AND TORSION
D,ALL,ROTX
!D,ALL,ROTZ
!D,ALL,ROTY
!D,ALL,UZ
!D,2,ALL ! SECOND NODES OF BEARING
!D,6,ALL
D,2,UY
D,2,UZ
D,6,UY
D,6,UZ

/COM, MODAL
RATIO = 4*ATAN(1)/30 ! in rad/sec 1/60!ist in u/min    ! 
ANTYPE,MODAL
CORIOLIS,ON,,,ON
NBF = 20
MODOPT,DAMP,NBF,,,ON
!MODOPT,DAMP,NBF,50000,1000000,ON ! 0 ist die untere Grenze und 10000 die Obere Grenze der Frequenz
*DO,I,1,10
OMEGA,SPIN(I)*RATIO ! UNIT FOR OMEGA IS RAD/SEC
MXPAND,NBF
SOLVE
*ENDDO
FINI