烦请帮我看下该命令流哪里有错误，谢谢

!!!!!!!!!!!!!!!!!!!!!!受弯构件的临界弯矩和极限承载力!!!!!!!
finish
/clear
/unis,SI
/filname,ansys4
!!!!!参数输入!!!!!
/prep7
H=1.232 !!! 外径 !!!!!
B1=0.4 !!!上翼缘宽度!!!!!
B2=0.2 !!!下翼缘宽度!!!!!
Tw=0.012 !!! 腹板厚度 !!!!!
Tf1=0.02 !!!上翼缘厚度!!!!!
Tf2=0.012 !!!下翼缘厚度!!!!!
L=16 !!! 计算跨度 !!!!!
TEC=L/1000 !!! 初始弯曲 !!!!!
Uymax=0.2 !!!弧长法控制参数!!!
FX=235.0e6 !!!fy (MPa)
E0=2.06e11 !!!弹性模量(MPa)
ET=E0*0.0 !!!切线弹性模量(MPa)
NU=0.3 !!!柏松比
FF=623214.1 !上边缘屈服弯矩对应的集中力
NN=40 !!!分段总数
DL=L/NN
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!定义单元类型
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/prep7
et,1,beam188
keyopt,1,1,1 !!考虑翘曲自由度
keyopt,1,3,2 !!定义形函数精度
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!定义截面
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
sectype,1,beam,I,,2
secoffset,cent
secdata,B2,B1,H,Tf2,Tf1,Tw,0,0,0,0,0,0,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!定义本构关系
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
mp,ex,1,E0
mp,prxy,1,NU
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!定义单元和节点
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,i,1,NN+1
n,i,(i-1)*DL,0,0
*enddo
n,NN+2,0,1e10,0
*do,ii,1,NN
e,ii,ii+1,NN+2
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!定义边界条件,施加荷载
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
d,1,ux,0
d,1,uy,0
d,1,uz,0
d,1,rotx,0 !简支梁不能绕x轴扭转
d,NN+1,uy,0
d,NN+1,uz,0
d,NN+1,rotx,0
f,NN/2+1,fy,-FF
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!首先进行静力解
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu
antype,0
pstres,on
solve
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!然后进行屈曲分析
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu
antype,1
bucopt,lanb,4
mxpand,4,,,1
solve
!!!显示设置
/PSF,PRES,NORM,2,0,1 !荷载显示箭头
/PBC,pres,,1
/eshape,1
/VIEW,1,1,1,1
/ANG,1
/REP,FAST !!构件尺寸

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!观察屈曲模态
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
set,first
plnsol,u,sum,0,1.0
nsort,u,sum
*get,usmax,sort,0,max !a1=最大合位移
*get,nmax,sort,0,imax
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!施加初始缺陷
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/prep7
upgeom,TEC/usmax,1,1,'ansys4','rst'
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!再进行非线性分析
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!随动强化多线模型
/prep7
tb,bkin,1,1,2,1
tbdata,,FX,ET
finish
/solu
antype,0
time,1
pstres,on
nlgeom,on
nsubst,20
arclen,on
ARCTRM,U,Uymax,NN/2+1,UY
outres,all,none
outres,nsol,all
!outres,all,all
solve
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!后处理，找寻最大稳定系数
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
*get,nset,active,0,set,nset
*do,iii,1,nset-1
set, , ,1, , , ,iii,
*get,fi0,active,0,set,time
*if,fi0,ge,fimax,then
fimax=fi0
*endif
*enddo
!!!后处理
/POST1
SET,LAST
SET,PREVIOUS
PLNSOL, U,x, 0,1.0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!后处理，画荷载位移曲线
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post26
NSOL,2,NN/2+1,U,Y,UX_2
XVAR,1
PLVAR,2,
PRVAR,2,