声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1196|回复: 1

[结构分析] 转子动力学

[复制链接]
发表于 2010-8-24 10:37 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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


本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2010-8-25 15:55 | 显示全部楼层
我的问题是这样的“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.”高手帮忙看下怎么回事     
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 15:34 , Processed in 0.115082 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表