我想做一个循环改变臂架约束，然后绘出其最大应力曲线

yujin0728

老是出现错误，不知道哪错了，还请高手指教！以下是源程序！
FINISH
/CLEAR,START
*DIM,vs,,126
*DIM,num,,126
*DIM,x1,,126
*DIM,y1,,126
dl=630*2  !每个斜拉杆的节距630毫米
A1=0
A2=630
B1=900      !前臂宽度
H1=750      !前臂高度
/PREP7
!(1)定义节点
K,(10:310:20),(A1::dl),H1,0
K,(20:300:20),(A2::dl),0,0
L,(10:290:10),(20:300:10)
L,(10:290:20),(30:310:20)
L,(20:280:20),(40:300:20)
!前补充弦杆
K,305,15*dl-0.5*A2,0.5*H1,0
K,320,15*dl+0.5*A2,0.5*H1,0
K,330,15*dl+1.5*A2,0.5*H1,0
L,300,305
L,305,310
L,310,320
L,305,320
L,320,330
NUMSTR,LINE,501
LGEN,2,ALL,,,,,B1,500,0
/VIEW,1,1,1,1
!底面横杆
L,(20:300:20),(520:800:20)
L,10,510    ！尾部弯杆
L,330,830    !头部横杆
/VIEW,1,1,2,3
!建立方向关键点
K,831,,H1/2,H1/2
K,832,,H1/2,B1+H1/2
LPLOT
ET,1,BEAM188
MP,EX,1,210E3
MP,PRXY,1,0.3
MP,DENS,1,7800E-9
R,10
!定义前臂型材
FRL1=80 !方钢管长度80mm
FRW1=80 !宽度80mm
FRT1=6    !厚度6mm
FCI1=20 !圆钢管内径40mm
FCO1=24   !外径48mm
FRL2=70   !矩形钢管长度70mm
FRW2=50  !矩形钢管宽度50mm
FRT2=4     !矩形钢管厚度4mm
SECTYPE,11,BEAM,HREC
SECOFFSET,CENT
SECDATA,FRL1,FRW1,FRT1,FRT1,FRT1,FRT1
SECTYPE,12,BEAM,CTUBE
SECOFFSET,CENT
SECDATA,FCI1,FCO1
SECTYPE,13,BEAM,HREC
SECOFFSET,CENT
SECDATA,FRL2,FRW2,FRT2,FRT2,FRT2,FRT2
！划分前臂网格
LSEL,S,,,30,58     !划分方钢管
LESIZE,ALL,50
LATT,1,10,1,,831,,11
LMESH,ALL
ALLSEL,ALL
LSEL,S,,,530,558     !划分方钢管
LESIZE,ALL,50
LATT,1,10,1,,832,,11
LMESH,ALL
ALLSEL,ALL
LSEL,S,,,1,29          !划分弦杆圆钢管48×4
LSEL,A,,,501,529
LSEL,A,,,59,61
LSEL,A,,,559,561
LESIZE,ALL,50
LATT,1,10,1,,,,12
LMESH,ALL
ALLSEL,ALL
LSEL,S,,,564,579          !划分底面横杆圆钢管48×4
LESIZE,ALL,50
LATT,1,10,1,,,,12
LMESH,ALL
ALLSEL,ALL
LSEL,S,,,62,63       !划分头部矩形钢管70×50×4
LSEL,A,,,562,563
LSEL,A,,,580
LESIZE,ALL,50
LATT,1,10,1,,,,13
LMESH,ALL
ALLSEL,ALL
/SOLU
ACEL,0,9.8,0
FK,320,FY,-2940
FK,820,FY,-2940
FK,(30:290:20),FY,-595.35
FK,(530:790:20),FY,-595.35
SAVE,front,DB,,ALL
*DO,i,1,126
num(i)=i
RESUME,front,DB,,0
NSEL,S,LOC,Y,H1
NSEL,R,LOC,X,A2-100+50*i,A2+300+50*i
NSEL,U,LOC,Z,10,B1-10
D,ALL,UY
ALLSEL,ALL
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,3*dl+250+50*i,3*dl-450+50*i
NSEL,U,LOC,Z,10,B1-10
D,ALL,UY
D,ALL,UZ
ALLSEL,ALL
SOLVE
FINISH
/POST1
/ESHAPE,1
PLNSOL,S,EQV
*GET,vs(i),plnsol,0,max
FINISH
*ENDO
/SOLU
SMOOTH,num,vs,,1,x1,y1,1
FINISH
/IMAGE,SAVE,Von_stre,BMP

经过测试，如果去掉循环的部分，程序没有问题！

yujin0728

求高手指导啊！！:'(