frp加固抗剪剥离命令流运行出错！

我在输入命令流计算前出现错误提示不能计算下去
Input displacement must be in ascending order for element 4160!
我是U型条加固的 弹簧单元有切向和法向两个方向的
命令流如下
大家帮我看看 为何不能计算啊！
急！


Finish
/clear,all
/com,structural
/prep7
!定义单元类型
et,1,solid65 ! 混凝土
et,2,link8 ! 钢筋
et,3,shell41 ! CFRP
et,4,solid45 ! 垫块
et,5,combin39 ! 弹簧单元 （沿切向Z向）
et,6,combin39 ! 弹簧单元 （沿法向X向）
et,7,combin39 ! 弹簧单元 （沿法向Y向）

keyopt,1,1,1 ! 不考虑大变形
keyopt,1,5,1 !混凝土的线性解输出控制：给出每个积分点的解
keyopt,1,6,3 !混凝土的线非性解输出控制：给出每个积分点的解
keyopt,1,7,1 ! 开裂后考虑应力松弛
keyopt,3,1,2 ! 作用在拉力状态下的刚度，压力状态下崩溃
keyopt,5,3,3 ！沿Z 向平动

keyopt,5,1,0 ! 卸载路径与加载路径相同
keyopt,5,2,1 ！不能承受压力
keyopt,5,4,0 ! 控制三唯轴向单元
keyopt,6,3,1 ！沿X 向平动
keyopt,6,4,0 ! 控制三唯轴向单元
keyopt,6,1,0 ! 卸载路径与加载路径相同
keyopt,6,2,1 ！不能承受压力
keyopt,7,3,2 ！沿Y 向平动
keyopt,7,4,0 ! 控制三唯轴向单元
keyopt,7,1,0 ! 卸载路径与加载路径相同
keyopt,7,2,1 ！不能承受压力

!定义实常数
r,1 ! 混凝土实常数
r,2,28.3e-6 ！箍筋d6的截面积
r,3,380.1e-6 ！纵筋d22的截面积
r,4,1.11e-5 ！cfrp的厚度

r,5, 0,0,2.59e-6,710.54,5.18e-6,142.018 !建切（z）向弹簧单元位移-荷载曲线
rmore,7.770e-6,2131.61,1.036e-5,2842.15,1.439e-5,2842.15
rmore,1.842e-5,2842.15,2.245e-5,2842.15,2.648e-5,2842.15
rmore,3.050e-5,2842.15,3.453e-5,2842.15,3.885e-5,2486.88
rmore,4.317e-5,2131.61,4.748e-5,1776.34,5.18e-5,1421.08
rmore,5.612e-5,1065.81,6.043e-5,710.54,6.475e-5,355.27
rmore.6.9066e-5,0.00,

r,6,0,0,8.74e-7,671.19,1.75e-6,1342.37 !建法（x）向弹簧单元位移-荷载曲线
rmore,2.62e-6,2013.56,3.50-6,2684.74,4.85-6,2684.74
rmore,6.21e-6,2684.74,7.57e-6,2684.74,8.93e-6,2684.74
rmore,1.03e-5,2684.74,1.17e-5,2684.74,1.31e-5,2349.15
rmore,1.46e-5,2013.56,1.60e-5,1677.97,1.75e-5,1342.37
rmore,1.89e-5,1006.78,2.04e-5,671.19,2.18e-5,335.59
rmore,2.33e-5,0.00,

！定义材料属性
!混凝土
mp,ex,1,2.99e10 !混凝土弹模
mp,prxy,1,0.2 !混凝土泊松比
tb,conc,1 ！定义非线性材料特性
tbdata,,0.5,0.9,2.14e6,-1 ！张开系数，闭合，单轴抗拉强度，单轴抗压
tb,miso,1,1,8 ！多线性等向强化
tbpt,,0.0001,2.99e6 ！应变应力关系
tbpt,,0.0003,3.9683e6
tbpt,,0.0006,7.2930e6
tbpt,,0.0009,9.9743e6
tbpt,,0.0012,12.0120e6
tbpt,,0.0016,13.7280e6
tbpt,,0.0020,14.3e6
tbpt,,0.0033,14.3e6
！箍筋
mp,ex,2,2.1e11 ！箍筋弹模
mp,prxy,2,0.25 ！箍筋泊松比
tb,bkin,2,1,2,1 ！双线性随动强化
tbdata,,200e6 ！d6的屈服应力
!纵向钢筋
mp,ex,3,2e11 ！纵筋弹模
mp,prxy,3,0.3 ！纵筋泊松比
tb,bkin,3,1,2,1 ！双线性随动强化
tbdata,,350e6 ！d22的屈服应力
!纤维布
mp,ex,4,2.35e11
mp,prxy,4,0 ！（0.3？）
tb,bkin,4,1,2,1
tbdata,,3550e6 ！抗拉强度
! 垫块
mp,ex,5,2e5
mp,prxy,5,0.3
tb,bkin,5,1,2
tbdata,,335,2000

!建立所有节点(150*260*1500)
n,1
n,7,0.15 ！x方向单元尺寸0.025
fill,1,7
ngen,11,7,1,7,1,,0.026 ！y方向单元尺寸0.026
ngen,31,1000,1,77,1,,,-0.05 ！z方向单元尺寸0.05
/view,1,1,1,1
nplot
!箍筋
type,2
real,2
mat,2
*do,ii,1065,1068,1
e,ii,ii+1
*enddo
*do,ii,1009,1058,7
e,ii,ii+7
*enddo
*do,ii,1009,1012,1
e,ii,ii+1
*enddo
*do,ii,1013,1062,7
e,ii,ii+7
*enddo
egen,12,4000,1,24,1

!纵向钢筋
type,2
real,3
mat,3

!受拉钢筋（3 d22）
*do,ii,1009,29009,1000
e,ii,ii+1000
*enddo
*do,ii,1011,29011,1000
e,ii,ii+1000
*enddo
*do,ii,1013,29013,1000
e,ii,ii+1000
*enddo

!受压钢筋（3 d22）
*do,ii,1065,29065,1000
e,ii,ii+1000
*enddo
*do,ii,1067,29067,1000
e,ii,ii+1000
*enddo
*do,ii,1069,29069,1000
e,ii,ii+1000
*enddo
!混凝土单元
type,1
real,1
mat,1
block,0,0.15,0,0.26,0,-1.5 ! x\y\z方向的长度
/view,1,1,2,3
lplot
/replot
lsel,s,line,,9,12,1
lesize,all,0.05
lsel,all
lsel,s,line,,5,7,2
lsel,a,line,,2,4,2
lesize,all,0.025
lsel,all
lsel,s,line,,1,3,2
lsel,a,line,,6,8,2
lesize,all,0.026
lsel,all
vmesh,all
/view,1,1,1,1
/eshape,1
eplot
nummrg,node
numcmp,node
eplot

! 垫块
wpoffs,,,-0.125 !工作平面沿z轴反向移动0.125
blc4,,,0.15,-0.020,0.05 ！以新工作平面坐标建长方体
vgen,2,2,,,,,-1.3 ！沿z轴1.3复制长方体
vgen,2,2,,,,0.28,-0.5 ！向y向0.28 z反向0.5的位置复制长方体、
vgen,2,4,,,,,-0.3

VSEL,s,volu,,2,5,1 ! 选4个垫块体
VATT,5,,4,0 ! 给垫块体赋值
LSEL,s,line,,19,55,12 ！划分尺寸
LSEL,a,line,,23,59,12
LESIZE,all,0.025,,,,,,,0
LSEL,s,line,,18,54,12
LESIZE,all,0.01,,,,,,,0
MSHAPE,0,3D ！指定单元形状为6面体
MSHKEY,1
VMESH,all

! cfrp单元
wpave,0,0,0 !工作平面移回到原点
wpoffs,,,-0.175 !工作平面移到（0，0，负0.175）
wprota,,90 ！绕x轴旋转90
blc4,,,0.15,0.05 ！底面建立长0.15宽0.05 0的体（实际是个面）

wpave,0,0,0 !工作平面移回到原点
wprota,,-90 ！绕x轴旋反转90
wpoffs,,,-0.175 !工作平面移到（0，0，负0.175）
wprota,,,-90 ！绕y轴旋反转90
blc4,,,0.05,0.26 ！侧面建立宽0.05长0.260的体（实际是个面）
agen,2,32,,,0.15 ！将此面复制到另一面（以坐标原点沿z 轴）
!复制U 型frp( 12次13条，全长加固)
*do,ii,1,12
agen,2,31,33,1,,,ii*(-0.1)
ii=ii+1
*end do

/replot
asel,s,area,,31,69,1 ！选择31至69号面
aplot
cm,frp,area
aatt,4,4,3,0 ！赋予属性
lsel,s,line,,70,214,12 ！两侧面划分尺寸
lsel,a,line,,66,210,12
lesize,all,0.026,,,,,,,0
lsel,s,line,,67,211,12
lsel,a,line,,71,215,12
lesize,all,0.025,,,,,,,0
lsel,s,line,,63,207,12 ！底面划分尺寸
lsel,a,line,,62,206,12
lsel,a,line,,64,208,12
lesize,all,0.025,,,,,,,0
mshape,0,2D ！指定单元形状为2D四边形
mshkey,1
amesh,all

allsel,all
eplot

！侧面切向弹簧单元
nsel,r,loc,z,-0.125,-0.175
nsel,a,loc,z,-0.225,-0.275
nsel,a,loc,z,-0.325,-0.375
nsel,a,loc,z,-0.425,-0.475
nsel,a,loc,z,-0.525,-0.575
nsel,a,loc,z,-0.625,-0.675
nsel,a,loc,z,-0.725,-0.775
nsel,a,loc,z,-0.825,-0.875
nsel,a,loc,z,-0.925,-0.975
nsel,a,loc,z,-1.025,-1.075
nsel,a,loc,z,-1.125,-1.175
nsel,a,loc,z,-1.225,-1.275
nsel,a,loc,z,-1.325,-1.375
nsel,u,loc,y,0.265,0.28
nsel,u,loc,y,-0.005,-0.02
nsel,u,loc,x,0.025,0.145
nplot
cm,cemian,node
type,5
real,5

secnum,none, !设界面编号
tshap,line ！定义接触目标为2D或3D的简单图形
eintf,0.0001,,low,,,,, ！用二维线单元连接重合的节点从低到高排列
allsel,all
eplot
！侧面法向弹簧单元
cmsel,s,cemian,node
nplot
type,6
real,6

secnum,none, !设界面编号
tshap,line ！定义接触目标为2D或3D的简单图形
eintf,0.0001,,low,,,,, ！用二维线单元连接重合的节点从低到高排列
allsel,all
eplot

!耦合侧面重合节点y方向位移
cmsel,s,cemian,node
nplot
cpintf,uy,0.0001

！底面切向弹簧单元
nsel,r,loc,z,-0.125,-0.175
nsel,a,loc,z,-0.225,-0.275
nsel,a,loc,z,-0.325,-0.375
nsel,a,loc,z,-0.425,-0.475
nsel,a,loc,z,-0.525,-0.575
nsel,a,loc,z,-0.625,-0.675
nsel,a,loc,z,-0.725,-0.775
nsel,a,loc,z,-0.825,-0.875
nsel,a,loc,z,-0.925,-0.975
nsel,a,loc,z,-1.025,-1.075
nsel,a,loc,z,-1.125,-1.175
nsel,a,loc,z,-1.225,-1.275
nsel,a,loc,z,-1.325,-1.375
nsel,u,loc,y,0.005,0.28
nsel,u,loc,y,-0.005,-0.02
nplot
cm,dimian,node
type,5
real,5

secnum,none, !设界面编号
tshap,line ！定义接触目标为2D或3D的简单图形
eintf,0.0001,,low,,,,, ！用二维线单元连接重合的节点从低到高排列
allsel,all
eplot

！底面法向弹簧单元
cmsel,s,dimian,node
nplot
type,7
real,6

secnum,none, !设界面编号
tshap,line ！定义接触目标为2D或3D的简单图形
eintf,0.0001,,low,,,,, ！用二维线单元连接重合的节点从低到高排列
allsel,all
eplot

!耦合底面重合节点x方向位移
cmsel,s,dimian,node
nplot
cpintf,ux,0.0001

！垫块与混凝土的约束设置
VSEL,s,volu,,1
ESLV,S
VSEL,s,volu,,2,5,1
NSLV,S,1
nsel,r,loc,y,0,0.260
nsel,u,loc,z,-0.100
nsel,u,loc,z,-0.60
nsel,u,loc,z,-0.90
nsel,u,loc,z,-1.4
CEINTF,0.25,ux,uy,uz, , , ,0,
eplot

allsel,all
eplot

!求解设置
/solu
antype,0
nsel,s,loc,z,-0.1
nsel,a,loc,z,-1.4
nsel,r,loc,y,-0.02
d,all,uy
d,all,ux
nsel,r,loc,z,-1.4
d,all,uz

*do,ii,0,60
allsel,all
nsel,s,loc,z,-0.6
nsel,a,loc,z,-0.9
nsel,r,loc,y,-0.28
f,all,fy,-(0.500+ii*0.200)
allsel,all
autots,on
time,ii+1
nsubst,200
lnsrch,on
pred,on
neqit,3000
outres,basic,all
ncnv,2
lswrite,ii+1
ii=ii+1
*enddo

allsel,all
EPLOT

/solu
allsel,all
LSSOLVE,1,61,1