tigerlily 发表于 2007-5-25 09:34

请大家帮我看看哪里出错了

模型很简单,就是模拟一个平板在空气中的振动模态,加上了空气阻尼。用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.

请教高手为什么会出现这种结果。我的命令流哪里处理问题。敬请指教。

ly030604 发表于 2008-12-11 11:22

回复 楼主 tigerlily 的帖子

我看不出来啊!不知您现在有没有解决啊?我很想知道啊!
ly030604@163.com
谢谢啦!
页: [1]
查看完整版本: 请大家帮我看看哪里出错了