请问关于ansys算出的滞回曲线怎么会这样?

本人前天发过请教位移加载错误的帖子,后把钢材本构关系中的TBDATA,,325,0,,,, 改为TBDATA,,325,2e3,,,, 就没有错误了,请问是什么原因?但是现在画出的滞回曲线非常奇怪,后面有两段位移不变,力反而急剧变化,请问是什么原因.命令流如下:!注：所用的单位为：mm,n,e3kg,e3kg/mm2,mpa,9810mm/s2 /filname,experiment
/title,frame model
!截面尺寸为150mm×150mm,其中钢为5mm的围板
/PREP7

!生成截面面积
blc5,,,120,120
blc5,,,112,112
aptn,all
numcmp,all
!定义单元与材料号，仅材料号有用；
!且该材料号的性质，将套用将来使用截面时所定义的性质；
!当然也可以通过GUI方式对材料号另赋值，但比较麻烦(不知命令方式?)
et,1,plane82
mp,ex,1,3e10 !混凝土
mp,ex,2,2e11 !钢
!divide面积，以划分截面
wpoff,-56,-56
wprota,,90
asbw,all
wpoff,,,-112
asbw,all
wprota,,,90
asbw,all
wpoff,,,112
asbw,all
numcmp,all
wpcsys,-1
asel,s,loc,x,0
asel,r,loc,y,0
*get,a1,area,0,num,min
aatt,1,,1 !关联材料号
asel,s,,,all
asel,u,,,a1
aatt,2,,1
allsel,all
esize,11.2
mshape,0,2d
mshkey,1
amesh,all !剖分面积
secwrite,xxjm,sect !写入文件xxjm.sect
!以上定义了一个复合截面，将来在188/189中使用。

finish
/clear,start
/config,nres,3000
/prep7

mp,ex,1,3.45e4!混凝土
mp,prxy,1,0.3
mp,dens,1,2.5e-9

TB,BISO,1,1,2,
TBTEMP,0
TBDATA,,17.5,0,,,,
tbplot,biso,1,,,yes

mp,ex,2,2.1e5　　!钢
mp,prxy,2,0.3
mp,dens,2,7.8e-9

TB,BISO,2,1,2,
TBTEMP,0
TBDATA,,325,2e3,,,,
tbplot,biso,2,,,yes

et,1,beam188
sectype,1,beam,mesh
secread,'xxjm','sect',' ',mesh
secplot,1,1

sectype,2,beam,I,I-shape
secoffset,cent
secdata,120,120,180,6,6,5
secplot,2,1

k,1
k,2,,1000
k,3,,2000
k,4,,3000
*do,i,1,3
lstr,i,i+1
*enddo
lgen,3,1,3,1,1500
*do,i,2,4
lstr,i,i+4
*enddo
*do,i,6,8
lstr,i,i+4
*enddo

lsel,s,line,,1,9
latt,1,1,1,,,,1
lesize,all,50
lmesh,all

lsel,s,line,,10,15
latt,2,1,1,,1,,2
lesize,all,50
lmesh,all

/VIEW, 1 ,1,1,1 !三维显示
/ANG, 1
/eshape,1
/DEVICE,VECTOR,0
eplot

!合并节点，整理编号
nummrg,node
numcmp,node

/solu
nsel,s,loc,y,0
d,all,,,,,,all
ANTYPE,STATIC
ALLSEL
nlgeom,off　　　　　　!在静态或完全瞬态中是否包含大变形
pred,off 　　　　　　　!在非线性分析中激活一个预测器
nropt,full,,on　　　!在静态或完全瞬态分析中，指定Newton-Raphson选项
lnsrch,on　　　　　!激活一个线性搜索与Newton-Raphson一起使用
sstif,on　　　　　　!在非线性分析中，是否激活应力刚度效应
NSUBST,50,1000,5　　!指定荷载步中所需只步数
cnvtol,u,,0.03,0　!为非线性分析指定收敛值
autots,off　　　　　!指定是否使用自动时间步长跟踪或荷载步跟踪
rescontrl,,all,1,1!输出子步
OUTRES,nsol,all　　　!控制结果的输出
OUTRES,rsol,all

*dim,disp1,,39
*dim,disp2,,39
*dim,disp3,,39

*create,ff,TXT,
*VREAD,disp1(1),disp1,txt,,IJK,39
(E3.0)！读荷载数组
*end
/input,ff,TXT,

*create,ff,TXT,
*VREAD,disp2(1),disp2,txt,,IJK,39
(E3.0)！读荷载数组
*end
/input,ff,TXT,

*create,ff,TXT,
*VREAD,disp3(1),disp3,txt,,IJK,39
(E3.0)！读荷载数组
*end
/input,ff,TXT,

*DO,I,1,39
time,i
DK,12,ux,-disp1(i)！位移加载
DK,11,ux,-disp2(i)
DK,10,ux,-disp3(i)
ACEL,0,-9810,0,
SOLVE
*ENDDO
save

!画出荷载位移曲线
FINISH
/POST26
NSOL,2,164,U,X,
RFORCE,3,164,F,X,
XVAR,2
PLVAR,3, , , , , , , , , ,