框架beam189特征值屈曲分析

xtddpla

昨天做一个矩形框架的特征值屈曲分析，采用beam189单元，扩展5阶屈曲模态，查看结果发现5个是一样的，和理论解相去甚远，原因何在？
                把命令流也附上：
fini
/clear
/filname,wending
/title,dancengkuangjia
/prep7
et,1,beam189
mp,ex,1,206e3
mp,prxy,1,0.3
sectype,1,beam,i
secdata,600,600,600,12,12,14
secplot,1
sectype,2,beam,i
secdata,600,600,600,14,14,16
secplot,2
k,1
k,2,6000
k,3,6000,6000
k,4,,6000
k,101,6000,,6000
k,102,,6000,6000
!dk,1,ux,,,,uy,uz
!dk,2,ux,,,,uy,uz
!dk,1,rotx
!dk,2,rotx
l,1,2
l,2,3
l,3,4
l,4,1
lesize,1,,,20
latt,1,,1,,101,,1
lmesh,1
lesize,2,,,20
latt,1,,1,,101,,2
lmesh,2
lesize,3,,,20
latt,1,,1,,102,,1
lmesh,3
lesize,4,,,20
latt,1,,1,,102,,2
lmesh,4

/solu
dianh1=node(0,0,0)
dianh2=node(6000,0,0)
d,dianh1,ux,,,,,uy,uz,rotx
d,dianh2,uy,,,,,uz,rotx
fk,3,fy,-1
fk,4,fy,-1
pstres,on
allsel,all
solve
fini

/solu
antype,buckle
bucopt,lanb,5
mxpand,5
outres,all,all
solve
fini
/post1
set,list

xtddpla

得出来的结果是：
*****  INDEX OF DATA SETS ON RESULTS FILE  *****
   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1 0.18136E+07         1         1         1
     2 0.18136E+07         1         2         2
     3 0.18136E+07         1         3         3
     4 0.18136E+07         1         4         4
     5 0.18136E+07         1         5         5
5个都一样，很奇怪，为什么？

koolin

你这个算法取得有点怪，计算模态可以用这个算法吗？我也是初学者，探讨一下。

xtddpla

我是参照王新敏书上做的，改用beam4做出结果了和理论解基本相同，但189就没有做出结果。 3# koolin