怎么输出屈曲分析的 临界载荷？？？谢谢

bluecmy

命令流如下，，，
t=1
L=1000
R=100
s=30
/PREP7
ET,1,SHELL93
ET,2,mass21
R,2,1,1,1
R,1,t
MP,EX,1,206e3
MP,PRXY,1,0.3
/view,1,1,1,1
/ang,1,60,zs,1
/ang,1,60,xs,1
K,1,0,0,0,
K,2,0,0,L
LSTR, 1, 2
circle,1,R
lsel,s,loc,z,0
ADRAG,all,0,0,0,,,1
ldele,1
aesize,all,s
MSHAPE,0,2D
MSHKEY,0
amesh,all


ksel,s,,,2
type,2
real,2
kmesh,all
allsel
nsel,s,loc,z,L
CERIG,node(0,0,L),ALL,ALL, , , ,
allsel
/SOL 静力及特征值屈曲分析
!加约束
lsel,s,loc,z,0
dl,all,,ux
dl,all,,uy
dl,all,,uz
!加荷载
f,node(0,0,L),mx,1
allsel

PSTRES,ON
SOLVE
FINISH
/SOLUTION
ANTYPE,1
MXPAND,3,0,0,0,0.001,
BUCOPT,LANB,3
SOLVE
FINISH
/POST1
SET,FIRST
PLDISP,2

jxxansys

/POST1
SET,FIRST
SET,LIST
PLDISP,2
FINISH