小弟摆不平一个很简单的ANSYS求解，麻烦过路的同仁、高手指教！感激不尽！

taoboy

程序是一个2层框架结构的简化模型，按照书本上的例子自己编的代码，具体参数如下：
弹性模量E：3e6，截面惯性矩：1000（自己按照这个参数设定了梁的截面尺寸），跨度240，层高：120，集中质量：题目未说，我自己按照它的公式计算了出来！
下面是代码，烦请各位高手运行一下，为什么明明按照书本上说的自振频率是4.023rad,10.71rad，而我自己计算的却相差太远了呢！
finish
/clear
/filename,framework
/title,framework
e0=3e6$h=120$l=240$zb1=10$zh1=10.6266$lb1=5$lh1=13.39
$zb2=12$zh2=12.6$lb2=6.2$lh2=15.7015
/prep7
et,1,beam3,,,,,,1
et,2,mass21,,,4   
r,1,zb1*zh1,(zb1*zh1**3)/12,zh1
r,2,zb2*zh2,(zb2*zh2**3)/12,zh2
r,3,lb1*lh1,(lb1*lh1**3)/12,lh1
r,4,lb2*lh2,(lb2*lh2**3)/12,lh2
r,5,8.96e5$r,6,17.92e5
mp,dens,1,0
mp,prxy,1,0.2
mp,ex,1,e0
n,1,
n,2,l,
n,3,0,h,
n,4,l,h
n,5,0,2*h
n,6,l,2*h
n,7,l/2,h
n,8,l/2,2*h
type,1$mat,1,$real,1
e,4,6$e,3,5
type,1$mat,1,$real,2
e,1,3$e,2,4
type,1$mat,1,$real,3
e,5,8$e,8,6
type,1$mat,1,$real,4
e,3,7$e,7,4
type,2$real,5
e,8
type,2$real,6
e,7
finish
/solu
antype,modal,
modopt,lanb,10
mxpand,10,,,yes
d,1,ux$d,1,uy$d,1,rotz
d,2,ux$d,2,uy$d,2,rotz
solve
finish
/post1
set,list
set,first$pldisp,2

[ 本帖最后由 taoboy 于 2010-6-22 21:56 编辑 ]