碳纤维加固，求助I！！！！！！！

*SET,H,80
*SET,B,250
*SET,L,1800
*SET,LF,100
*SET,LS,50
*SET,A,20
*SET,PI,ACOS(-1)
*SET,SR,PI*4**2
*SET,RO,PI*3**2/200/H
*SET,CB,300
*SET,CH,0.167
*SET,CL,L-2*LS
*SET,F1,0.002
*SET,F2,30/LF
!***************************定义完参数,开始建立模型
/PREP7
ET,1,SOLID65
ET,2,LINK8
ET,3,SHELL41
KEYOPT,3,1,2
ET,4,SOLID45
MP,EX,1,18000
MP,PRXY,1,0.2
!**************************定义MISES本构关系
TB,KINH,1,1,7
TBPT,,0.0001,1.8
TBPT,,0.0004,6.66
TBPT,,0.0008,11.84
TBPT,,0.0012,15.54
TBPT,,0.0016,17.76
TBPT,,0.002,18.5
TBPT,,0.0033,18.5
TB,CONC,1,1,9
TBDATA,,0.3,0.5,1.75,-1
!**************************钢筋双线性随动强化
MP,EX,2,2.1E5
MP,PRXY,2,0.3
TB,BKIN,2,1,2,1
TBDATA,,235,0
!**************************碳纤维线弹性
MP,EX,3,2.35E5
MP,PRXY,3,0
!**************************支座及整体钢筋
MP,EX,4,2.1E5
MP,PRXY,4,0.3
R,1,2,RO,0,0,,,
RMORE,,,,,,,
R,2,SR,,
R,3,SR/2,,
R,4,CH,,,,,,
RMORE
!**************************建立板模型
BLOCK,0,-B,0,H,-LF/2,-(L/2-LS)
LPLOT
/PNUM,LINE,1
/REPLOT
LGEN,2,8,,,50,,,,0
LGEN,2,13,,,100,,,,0
ADRAG,13,,,,,,9
ADRAG,14,,,,,,9
/PNUM,LINE,0
/REPLOT
VSBA,1,7
VSBA,3,8
!**************************建立受力钢筋模型
LSEL,S,LOC,Z,-LF/2
LSEL,R,LOC,Y,0
LGEN,2,ALL,,,,A,,,0
LSEL,S,LOC,Z,-50
LSEL,R,LOC,Y,A
ADRAG,ALL,,,,,,9
ALLSEL,ALL
VSBA,2,15
VSBA,4,14
VSBA,1,13
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,A
LSEL,A,LOC,Y,H
LESIZE,ALL,50
LSEL,S,LOC,Z,-LF/2
LSEL,A,LOC,Z,-(L/2-LS)
LSEL,U,LOC,Y,0
LSEL,U,LOC,Y,A
LSEL,U,LOC,Y,H
LESIZE,ALL,20
ALLSEL,ALL
TYPE,1
MAT,1
REAL,1
ESYS,O
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL
ASEL,S,LOC,Z,-LF/2
!**************************完成另一半板模型
TYPE,1
EXTOPT,ESIZE,LF/50,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,O
VEXT,ALL,,,0,0,LF
ASEL,S,LOC,Z,LF/2
TYPE,1
EXTOPT,ESIZE,(L/2-LF/2-LS)/50,0
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,0
VEXT,ALL,,,0,0,L/2-LF/2-LS,
ALLSEL,ALL
!**************************生成支座模型
ASEL,S,LOC,Z,L/2-LS
TYPE,1
EXTOPT,ESIZE,1,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,0
VEXT,ALL,,,0,0,25
ASEL,S,LOC,Z,L/2-LS+25
TYPE,1
EXTOPT,ESIZE,2,0
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,0
VEXT,ALL,,,0,0,LS
ASEL,S,LOC,Z,-(L/2-LS)
TYPE,1
EXTOPT,ESIZE,1,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,0
VEXT,ALL,,,0,0,-25
ASEL,S,LOC,Z,-(L/2-LS+25)
TYPE,1
EXTOPT,ESIZE,2,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,1
REAL,1
ESYS,0
VEXT,ALL,,,0,0,-LS
ALLSEL,ALL
EPLOT
!**************************划分钢筋单元
LSEL,S,LOC,X,-200
LSEL,A,LOC,X,-100
LSEL,R,LOC,Y,A
TYPE,2
MAT,2
REAL,2
ESYS,0
LMESH,ALL
LSEL,S,LOC,X,0
LSEL,R,LOC,Y,A
TYPE,2
MAT,2
REAL,3
ESYS,0
LMESH,ALL
ALLSEL,ALL
!**************************划分碳纤维单元
KSEL,S,LOC,Z,850
KSEL,R,LOC,Y,0
KSEL,R,LOC,X,0
KGEN,2,ALL,,,-CB/2,,,,0
LSTR,41,97
ALLSEL,ALL
LSTR,41,4,
ADRAG,221,,,,,,222
ASEL,S,,,168
LSLA,S
LESIZE,ALL,50,,,,,,,1
TYPE,3
MAT,3
REAL,4
ESYS,0
MSHAPE,1,2D
AMESH,ALL
ALLSEL,ALL
NUMMRG,NODE,,,,LOW
NUMCMP,NODE
!**************************划分支座单元
ASEL,S,LOC,Y,0
ASEL,R,LOC,Z,L/2-25,L/2+25
TYPE,4
EXTOPT,ESIZE,1,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,4
REAL,4
VEXT,ALL,,,0,-20
ASEL,S,LOC,Y,0
ASEL,R,LOC,Z,-(L/2-25),-(L/2+25)
TYPE,4
EXTOPT,ESIZE,1,0,
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
MAT,4
REAL,4
ESYS,0
VEXT,ALL,,,0,-20,
ALLSEL,ALL
FINISH
!**************************设置边界
/SOLU
NSEL,S,LOC,X,0
/GO
DSYM,SYMM,X
NSEL,S,LOC,Z,L/2
NSEL,R,LOC,Y,-20
/GO
D,ALL,,,,,,UY
NSEL,S,LOC,Z,-L/2
NSEL,R,LOC,Z,-20
/GO
D,ALL,,,,,,UY,UZ
!**************************加载1
ANTYPE,0
NLGEOM,1
NROPT,FULL
EQSLV,SPAR,,0
ASEL,S,LOC,Y,H
/GO
SFA,ALL,1,PRES,F1
ALLSEL,ALL
OUTRES,ALL,ALL
TIME,1
AUTOTS,1
NSUBST,10,20,5,1
CNVTOL,F,,0,05,2,,
NEQIT,30,
PRED,ON,,ON
ESEL,MAT,3
EKILL,ALL
ALLSEL,ALL
SOLVE
!**************************加载2
ASEL,S,LOC,Y,H
ASEL,R,LOC,Z,-LF/2,LF/2
/GO
SFA,ALL,1,PRES,F2
OUTRES,ALL,ALL
TIME,2
AUTOTS,0
NSUBST,30,,,1
PRED,-1
ARCLEN,1,10,0,
ESEL,MAT,3
EALIVE,ALL
ALLSEL,ALL
SOLVE
FINISH
!**************************查看结果
/POST1
ESEL,MAT,1
SET,1,LAST,1,
PLNSOL,S,Z,0,1
SET,2,LAST,1,
PLNSOL,S,Z,0,1
FINISH
!**************************查看碳纤维结果
/POST26
ALLSEL,ALL
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,0
*GET,N1,NODE,0,NUM,MAX
NSOL,2,N1,U,Y
PROD,3,2,,,,,,-1,1,1,
XVAR,3
PLVAR,1
FINISH



就是过不去！！各位高手看看问题出在哪里？急需帮助！

是不是施加约束的问题，我粗略的看了一下

恩，我运行了一下，警告说约束有问题

我现在也在做这方面的东西，30米T梁加固，也出了很多问题，不知道是什么原因，楼主可以给我发一份你修正过的命令流吗，我的邮箱是jiandazhxch@163.com，小弟不胜感激