BEAM189压杆屈曲分析特征值比非线性值小

aslkjddfkljakl

我用BEAM189做压杆屈曲分析。
先做特征值分析，结果与理论值一样。然后接着做非线性屈曲分析，结果竟然比特征值分析所得结果大，前者为42346，后者为time=43833。
按照我的理解，42346N就是临界屈曲力的大小，而43833就是非线性分析时导致压杆屈曲的力的大小。
若是如此，则后者应该比前者小。但为什么大呢。下面是我的命令。求助呀。


/FILNAME,file,0
/TITLE,file
/PREP7  
!*  
ET,1,BEAM189
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,206e9  
MPDATA,PRXY,1,,0.3  
SECTYPE,   1, BEAM, RECT, , 0   
SECOFFSET, CENT
SECDATA,0.05,0.05,0,0,0,0,0,0,0,0,0,0   
K,1,0,0,0,  
K,2,5,0,0,  
LSTR,       1,       2  
FLST,5,1,4,ORDE,1   
FITEM,5,1   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,100, , , , ,1
!*  
LMESH,       1  
FLST,2,201,1,ORDE,2
FITEM,2,1   
FITEM,2,-201
!*  
/GO
D,P51X, , , , , ,UZ,ROTX,ROTY, , ,  
FLST,2,1,1,ORDE,1   
FITEM,2,2   
!*  
/GO
D,P51X, , , , , ,UY,UZ,ROTX,ROTY, ,
FLST,2,1,1,ORDE,1   
FITEM,2,1   
!*  
/GO
D,P51X, , , , , ,UX,UY,UZ,ROTX,ROTY,
FLST,2,1,1,ORDE,1   
FITEM,2,2   
!*  
/GO
F,P51X,FX,-1
SAVE
FINISH  
/SOL
!*  
ANTYPE,0
PSTRES,1
/STATUS,SOLU
SOLVE   
FINISH  
/SOLU   
!*  
ANTYPE,1
!*  
BUCOPT,LANB,5,0,0,  
MXPAND,5,0,1000000,1,0.001,
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
SET,LIST
SAVE
!*  
*GET,load,ACTIVE, ,SET,FREQ
FINISH  
/PREP7  
FINISH  
/SOL
FINISH  
/PREP7  
SAVE
SAVE
/BATCH  
/COM,ANSYS RELEASE 13.0    UP20101012       01:00:07    02/27/2012
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/COM,ANSYS RELEASE 13.0    UP20101012       01:00:07    02/27/2012
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
/SOLU   
FINISH  
/POST1  
SET,LIST
SET,LIST,999
SET,,, ,,, ,1   
!*  
*GET,solution,ACTIVE, ,SET,FREQ
FINISH  
/PREP7  
UPGEOM,0.001,1,1,'file','rst',' '   
!*  
FINISH  
/SOL
!*  
ANTYPE,0
ANTYPE,0
NLGEOM,1
NSUBST,50,10000,50  
OUTRES,ERASE
OUTRES,ALL,ALL  
LNSRCH,1
TIME,1.15*solution  
DOFSEL,S,FX
FSCALE,1.15*solution,   
DOFSEL,ALL  
!*  
/STATUS,SOLU
SOLVE

wulei

不是很懂啊