马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
自己想做一个转子模态分析,以下是我给出的程序,我的模态分析用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
|