为什么非线型屈曲分析得不到下降段？

youtaohit

为什么非线型屈曲分析得不到下降段？急！！！！！！

将特征值屈曲分析得到的一阶屈曲模态缩小后作为初始缺陷对结构进行非线型屈曲分析，然而得不到下降段，为什么？先谢谢了！
        问题是一个压杆的稳定，一下是命令流：

FINISH
/CLEAR,NOSTART
/FILNAME,BUCKLING
/TITLE,BUCKLING OF A BAR BY BEAM4

!特征值屈曲分析

/prep7
ET,1,BEAM4,,1,,,,1
MP,EX,1,2.0E11
MP,PRXY,1,0.3
R,1,2.16e-2,2.06e-4,2.06e-4,0.3,0.3
K,1,0,0,0,
K,2,0,10,0,
LSTR,1,2
LSEL,ALL
LESIZE,ALL,,,20
LMESH,ALL
d,1,ux,0
d,1,uy,0
d,1,uz,0
d,1,roty,0
d,2,ux,0
d,2,uz,0
d,2,roty,0
f,2,fy,-1
FINISH


/SOLU
ANTYPE,STATIC
EQSLV,SPAR
PSTRES,ON
SOLVE
FINISH


/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,4
MXPAND,4
SOLVE
FINISH
/POST1
SET,FIRST
*get,freq,mode,1,freq
PLDISP,1
finish


!非线型屈曲分析

/comfig,nres,100
/prep7
tb,biso,1,1,2,                       !修改材料模型
tbtemp,0
tbdata,,2.0e8,0,,,,
upgeom,0.01,1,1,'BUCKLING','rst',''  !引入几何缺陷
finish


/solu
antype,static
f,2,fy,-freq*1.2
nlgeom,on               !打开大变形选项
outres,all,all          !输出所有子步的结果
arclen,1,0,0,           !打开弧长法
arclen,u,0.2,11,ux      !终止条件设置
nsubst,200,,,1          !指定荷载步的子步数
solve
finish


/post26
/axlab,x,disp          !指定X轴的标注
/axlab,y,force         !指定Y轴的标注
nsol,2,11,u,x          !指定2号变量保存11号节点的X方向的位移计算结果
rforce,3,1,f,y         !指定3号变量保存1号节点的Y方向的支反力计算结果
xvar,2
plvar,3                !绘出位移——反力图
prvar,2,3              !列表显示位移和反力随时间的变化
finish

[ 本帖最后由 xinyuxf 于 2007-7-22 16:09 编辑 ]