/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.06E+11
PGXY = 7.9E+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 = 4 ! SHAFT SECTION PROPERTIES


*DIM,DIAM,ARRAY,Konstante1
DIAM(1) = 1E-1
DIAM(2) = 1.5E-1
DIAM(3) = 1.5E-1
DIAM(4) = 1E-1
!DIAM(1) = 2E-1
!DIAM(2) = 3E-1
!DIAM(3) = 3E-1
!DIAM(4) = 2E-1

*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,20,1.401,1.401,1.401,0.002,0.00136,0.00136
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.2
N,4 ,0.3
N,5 ,0.4

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,20
!E,2,6
/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,6,ALL ! SECOND NODES OF BEARING
D,7,ALL


/COM, MODAL
RATIO = 4*ATAN(1)/30 ! in rad/sec 1/60!ist in u/min    ! 
ANTYPE,MODAL
CORIOLIS,ON,,,ON
NBF = 30
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
/POST1
/OUT,
PRCAMP,,1.,RPM
PLCAMP,,1.,RPM
!*GET,F1,CAMP,1,FREQ,6
!*GET,F2,CAMP,2,FREQ,6
!*GET,F3,CAMP,3,FREQ,6
!*GET,F4,CAMP,4,FREQ,6
!*DIM,LABEL,CHAR,1,4
!*DIM,VALUE,ARRAY,4,3
!LABEL(1,1) = 'WHIRL BW'
!LABEL(1,2) = 'WHIRL FW'
!LABEL(1,3) = 'WHIRL BW'
!LABEL(1,4) = 'WHIRL FW'
!*VFILL,VALUE(1,1),DATA,10747
!*VFILL,VALUE(1,2),DATA,F1*60
!*VFILL,VALUE(1,3),DATA,ABS((F1*60)/10747)
!*VFILL,VALUE(2,1),DATA,19665
!*VFILL,VALUE(2,2),DATA,F2*60
!*VFILL,VALUE(2,3),DATA,ABS((F2*60)/19665)
!*VFILL,VALUE(3,1),DATA,39077
!*VFILL,VALUE(3,2),DATA,F3*60
!*VFILL,VALUE(3,3),DATA,ABS((F3*60)/39077)
!*VFILL,VALUE(4,1),DATA,47549
!*VFILL,VALUE(4,2),DATA,F4*60
!*VFILL,VALUE(4,3),DATA,ABS((F4*60)/47549)
SAVE,TABLE_1
FINI
/CLEAR,NOSTART