求助！请各位看看我的命令流，为何计算得到的残余应力完全不正确？

我想通过ANSYS模拟带圆孔平板试件的反复加卸载过程，并获得某一级荷载卸载之后的残余应力分布图。根据荷载步文件法编写了如下命令流，但计算结果貌似完全不正确。以某一级荷载卸载之后的应力分布图为例（如图所示），缺陷处的应力最大值与没卸载之前相差不大，而且缺陷形状从圆形变成了矩形，观察卸载过程中的应力变化，发现应力最大值总是先减小后增大。一直找不到问题所在很是困扰，请各位能够不吝赐教，指出问题所在，谢谢！
/PREP7
ET, 1, PLANE183
KEYOPT,1,3,3
R,1,4
MP,EX,1,206000
MP,PRXY,1,0.28
TB,KINH,1,1,11
TBTEMP,0
TBPT, ,0.00194,399.64
TBPT, ,0.01111,400
TBPT, ,0.01268,420
TBPT, ,0.01668,440
TBPT, ,0.02146,460
TBPT, ,0.02726,480
TBPT, ,0.03467,500
TBPT, ,0.04487,520
TBPT, ,0.06138,540
TBPT, ,0.11673,560
TBPT, ,0.15646,560
RECTNG,0,130,0,30
CYL4,65,15,2
ASBA,1,2
SMRT,1
MSHAPE,0,2D
MSHKEY,0
AMESH,ALL
FINISH
/SOLU
ANTYPE,0
DL,4, ,ALL
OUTRES,ALL,ALL
PRED,ON, ,ON
LNSRCH,1
AUTOTS,1
NSUBST,10
TIME,1
SFL,2,PRES,-41.6667
LSWRITE,1
TIME,2
SFL,2,PRES,0
LSWRITE,2
TIME,3
SFL,2,PRES,-83.3333
LSWRITE,3
TIME,4
SFL,2,PRES,0
LSWRITE,4
TIME,5
SFL,2,PRES,-125
LSWRITE,5
TIME,6
SFL,2,PRES,0
LSWRITE,6
TIME,7
SFL,2,PRES,-166.6667
LSWRITE,7
TIME,8
SFL,2,PRES,0
LSWRITE,8
TIME,9
SFL,2,PRES,-208.3333
LSWRITE,9
TIME,10
SFL,2,PRES,0
LSWRITE,10
TIME,11
SFL,2,PRES,-250
LSWRITE,11
TIME,12
SFL,2,PRES,0
LSWRITE,12
TIME,13
SFL,2,PRES,-291.6667
LSWRITE,13
TIME,14
SFL,2,PRES,0
LSWRITE,14
TIME,15
SFL,2,PRES,-333.3333
LSWRITE,15
TIME,16
SFL,2,PRES,0
LSWRITE,16
TIME,17
SFL,2,PRES,-375
LSWRITE,17
TIME,18
SFL,2,PRES,0
LSWRITE,18
TIME,19
SFL,2,PRES,-416.6667
LSWRITE,19
TIME,20
SFL,2,PRES,0
LSWRITE,20
TIME,21
SFL,2,PRES,-458.3333
LSWRITE,21
TIME,22
SFL,2,PRES,0
LSWRITE,22
CNVTOL,F, ,0.02,2
LSSOLVE,1,22
FINISH
/POST1
SET, , ,1, ,16
PLNSOL,S,EQV,0,1.0
FINISH