关于特征值稳定的问题，急

小妹是个初学者，做关于钢桁架整体稳定性的研究，遇到了一些问题，每次运行特征值稳定性计算时都会遇到这个问题，请论坛里的诸位大神帮我看看我的命令流到底是哪里出错了呗，谢谢啦。
如果和小妹研究方向差不多的我们也可以一起研究

/FILNAME GHJ 100
/PREP7
ET,1,LINK8
R,1,0.00722 !SXG 465*10
R,2,0.001306 !HG 108*8
R,3,0.00162 !ZC 133*8
R,4,0.001946 !XFG 159*8
R,5,0.011595 !XXG 500*15
MP,EX,1,2.06E11
MP,DENS,1,7.8E3
MP,PRXY,1,0.3
acel,0,9.81,0
SECTYPE, 1, BEAM, CTUBE, SXG, 0
SECOFFSET, CENT
SECDATA,0.2275,0.2325,0,0,0,0,0,0,0,0
SECTYPE, 2, BEAM, CTUBE, HG, 0
SECOFFSET, CENT
SECDATA,0.05,0.054,0,0,0,0,0,0,0,0
SECTYPE, 3, BEAM, CTUBE, ZC, 0
SECOFFSET, CENT
SECDATA,0.0625,0.0665,0,0,0,0,0,0,0,0
SECTYPE, 4, BEAM, CTUBE, XFG, 0
SECOFFSET, CENT
SECDATA,0.0755,0.0795,0,0,0,0,0,0,0,0
SECTYPE, 5, BEAM, CTUBE, XXG, 0
SECOFFSET, CENT
SECDATA,0.2425,0.25,0,0,0,0,0,0,0,0
k,1,-45,-7,0
*do,i,1,9,1
k,i+1,-45+i*10,-7,0
*enddo
*do,i,1,9,1
l,i,i+1,
*enddo
k,11,-50.0,0,2.5
*do,i,1,20,1
k,i+11,-50+i*5,0,2.5
*enddo
*do,i,11,30,1
l,i,i+1,
*enddo
k,32,-50.0,0,-2.5
*do,i,1,20,1
k,i+32,-50+i*5,0,-2.5
*enddo
*do,i,32,51,1
l,i,i+1,
*enddo
*do,i,1,10,1
l,i,2*i+9,
*enddo
*do,i,1,10,1
l,i,2*i+11,
*enddo
*do,i,1,10,1
l,i,2*i+30,
*enddo
*do,i,1,10,1
l,i,2*i+32,
*enddo
*do,i,11,31,1
l,i,i+21,
*enddo
*do,i,11,30,1
l,i,i+22,
*enddo
*do,i,12,31,1
l,i,i+20,
*enddo
*do,i,1,10,1
l,i,2*i+10,
*enddo
*do,i,1,10,1
l,i,2*i+31,
*enddo
LPLOT
*do,i,10,49,1
ESIZE, , 1
MAT,1
REAL,1
LMESH,i
*enddo
LPLOT
*do,i,90,110,1
ESIZE, , 1
MAT,1
REAL,2
LMESH,i
*enddo
LPLOT
*do,i,111,150,1
ESIZE, , 1
MAT,1
REAL,3
LMESH,i
*enddo
LPLOT
*do,i,50,89,1
ESIZE, , 1
MAT,1
REAL,4
LMESH,i
*enddo
LPLOT
*do,i,1,9,1
ESIZE, , 1
MAT,1
REAL,5
LMESH,i
*enddo
/SOLU
ANTYPE,STATIC
EPLOT
D,43,UX,,,,,UY,UZ
D,52,UY,,,,,UZ
D,1,UZ
D,21,UZ
*do,i,1,42,1
F,i,FY,-1
*enddo
*do,i,1,21,1
F,i,FZ,-1
*enddo
*do,i,43,52,1
F,i,FZ,-1
*enddo
PSTRES,ON
SOLVE
FINISH
/solu
antype,buckle
bucopt,lanb,1
mxpand,5
outres,all,all
solve
finish
/post1
set,list
set,1,1
pldisp
finish