马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
模型很简单,就是模拟一个平板在空气中的振动模态,加上了空气阻尼。用modal projection方法。
命令流如下:
/PREP7
lbeam=2e-4
wbeam=2e-5
tbeam=2e-6
d_el=1e-5
pamb=1e5
visc=18.3e-6
pref=1e5
mfp=64e-9
Knud=mfp/d_el
ET, 1, FLUID136,1
ET, 2, SOLID45
MP,VISC,1,VISC
R,1,d_el,,,pamb
RMORE,pref,mfp
MP,EX, 2, 1.7e11
MP,NUXY, 2, 0.3
MP,DENS, 2, 2330
rectng,0,lbeam,0,wbeam
type,1
mat,1
real,1
smrtsize,4
amesh,all
esize,,2
type,2
mat,2
real,2
vext,all,,,,,tbeam
nsel,s,loc,x,0
nsel,a,loc,x,lbeam
nsel,a,loc,y,0
nsel,a,loc,y,wbeam
nsel,r,loc,z,-1e-9,1e-9
d,all,pres
nsel,all
esel,s,type,,1
nsle,s,1
cm,FLUN,node
allsel
nsel,s,loc,x,0
nsel,a,loc,x,lbeam
d,all,ux
d,all,uy
d,all,uz
allsel
fini
/solu
antype,modal
modopt,lanb,2
eqslv,sparse
mxpand,2
solve
fini
/post1
RMFLVEC
fini
/solu
abextract,1,2
finish
结果很奇怪,出现了负的BETA值,如下:
Mass Matrix Multiplier: ALPHADMP = 1446.86181.
Stiffness Matrix Multiplier: BETADMP= -9.581867706E-13.
请教高手为什么会出现这种结果。我的命令流哪里处理问题。敬请指教。 |