电机转子临界转速的计算程序（模态分析）

以下是我计算的电机转子临界转速的程序，不知道有没有问题
FINISH
/CLEAR,NOSTART
/FILNAME,BEAM4
/TITLE,BEAM4,62MOTOR
/PREP7
ET,1,BEAM4
MP,EX,1,210E9
MP,DENS,1,7800
MP,PRXY,1,0.3
R,1,0.5027E-4,0.0201E-8,0.0201E-8,0.008,0.008 ！六种截面参数
R,2,0.7854E-4,0.0491E-8,0.0491E-8,0.01,0.01
R,3,1.7671E-4,0.2485E-8,0.2485E-8,0.015,0.015
R,4,0.3848E-4,0.0118E-8,0.0118E-8,0.007,0.007
R,5,0.9195E-4,0.0673E-8,0.0673E-8,0.01082,0.01082
R,6,1.7671E-4,0.2485E-8,0.2485E-8,0.015,0.015
N,1
N,13,0.012
FILL
N,14,0.017
N,16,0.027
FILL
N,17,0.032
N,29,0.044
FILL
REAL,1
E,1,2
EGEN,3,1,1
E,26,27
EGEN,3,1,4
REAL,2
E,4,5
E,5,6
E,24,25
E,25,26
REAL,3
E,6,7
E,7,8
E,8,9
E,21,22
E,22,23
E,23,24
REAL,4
E,9,10
E,10,11
E,19,20
E,20,21
REAL,5
E,11,12
E,12,13
E,17,18
E,18,19
REAL,6
E,13,14
EGEN,4,1,25
/ESHAPE,1
EPLOT
FINISH
/SOLU
ANTYPE,MODAL
D,1,UX ！1，29为两个端点，为轴承处
D,1,UY
D,1,UZ
D,29,UX
D,29,UY
D,29,UZ
以下采用gui操作，模态扩展为四阶
ansys计算结果和理论计算误差为0.33%
没有考虑陀螺效应，不知道对不对，请高手指点。