转子动力学

不懂

自己想做一个转子模态分析，以下是我给出的程序，我的模态分析用qrdamp算法，出现问题，始终解决不了，主要是我要在模态分析的基础上做模态叠加分析，而damp算法无法实现，所以请求高手帮忙看下！谢谢

/filename,zhuanzi
/title,motai
/prep7
k,1,0,
k,2,125e-3
k,3,250e-3,
k,4,750e-3
k,5,1000e-3
k,6,1250e-3
k,7,1500e-3
k,8,1750e-3
k,9,2000e-3
k,10,2250e-3,   ！轴承位置
k,11,2500e-3
k,12,3000e-3
k,13,3500e-3
k,14,4500e-3
k,15,5500e-3
k,16,6500e-3
k,17,7500e-3
k,18,7720e-3,   ！轴承位置
k,19,8500e-3
k,20,9500e-3
k,21,10500e-3
k,22,11500e-3
k,23,12050e-3,  ！轴承位置
k,24,13500e-3
k,25,14150e-3,
!建立直线
*do,i,1,24,1
l,i,i+1
*enddo

!* 设置单元属性
ET,1,188
ET,2,14
keypoint,2,2,2
r,2,5e11,0.1e4
et,3,14
keypoint,2,2,3
r,3,5e11,0.1e4
ET,4,14
keypoint,2,2,1
r,4,5e11,0.1e4
mp,ex,1,2.06e11  
mp,PRXY,1,0.3
mp,dens,1,7.850e3

mp,ex,2,2.06e11  
mp,PRXY,2,0.3
mp,dens,2,7.01250e3
LESIZE,all,20e-3
SECTYPE,   1, BEAM, CSOLID, sec01, 0
SECOFFSET, CENT
SECDATA,8800e-3/2,20,20  
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,   1
!*  
LMESH, 1  ,2

SECTYPE,   2, BEAM, CSOLID, sec02, 0
SECOFFSET, CENT
SECDATA,470e-3/2,20,10
TYPE,   1   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,   2
lmesh,3,7,1

SECTYPE,   3, BEAM,CSOLID, sec03, 0
SECOFFSET, CENT
SECDATA,500e-3/2,20,10
TYPE,   1   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,   3
LMESH,8,11,1


SECTYPE, 4, BEAM, CSOLID, sec04, 0
SECOFFSET, CENT
SECDATA,420e-3/2,20,10
TYPE,   1   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,  4
LMESH,12,16,1

SECTYPE,5,BEAM, CSOLID, sec05, 0
SECOFFSET, CENT
SECDATA,470e-3/2,20,10
TYPE,   1   
MAT,      2
REAL,   
ESYS,       0   
SECNUM,   5  
LMESH,17,19,1

SECTYPE,   6, BEAM, CTUBE, sec06, 0
SECOFFSET, CENT
SECDATA,570e-3/2,20,10
TYPE,   1   
MAT,      2
REAL,   
ESYS,       0   
SECNUM,  6
LMESH, 20,24,1

N,1000,2250e-3,-0.1
N,1001,7720e-3,-0.1
N,1002,12050e-3,-0.1

N,1010,2250e-3,0.00,0.1
N,1011,7720e-3,0.00,0.1
N,1012,12050e-3,0.00,0.1

type,2
real,2
e,106,1000
e,382,1001
e,582,1002
e,106,1010
e,382,1011
e,582,1012

type,3
real,3
e,106,1000
e,382,1001
e,582,1002
e,106,1010
e,382,1011
e,582,1012
type,4
real,4
e,106,1000
e,382,1001
e,582,1002
e,106,1010
e,382,1011
e,582,1012
!显示梁模型
/eshape,0
eplot

!加载
/solu
d,all,ux
d,1000,uy,-0.01
d,1000,ux,0,,,,uz,rotx,roty,rotz
d,1001,all
d,1002,all
d,1010,all
d,1011,all
d,1012,all

f,625,fy,-9e4
f,700,mz,-3e4
acel,,9.8
/solu
ANTYPE,static
OMEGA,3000
coriolis,on,,,on
pstres,on
solve
finish
/solu
ANTYPE,MODAL
modopt,QRDAMP,40,400,10000,on
MXPAND,40
betad,5e-5
coriolis,on,,,on
OMEGA,3000
pstres,on
solve
fini

不懂

我的问题是这样的“Number of matrix rows for BEAM188 elements in EMAT file does not match with that on data base.  Please rerun the QRdamp eigenproblem with EMATWRITE,YES to force writing of element matrix file.”高手帮忙看下怎么回事