求助：计算过程错误，那个高手帮我看看，谢谢

!***************************************************************************
!**********自然边坡稳定性分析***********************************
!***************************************************************************
!清理内存
fini
/clear
/prep7

!定义关键点
k,1,-10,,-10
k,2,-10,0,0
k,3,0,0,0
k,4,22 ,, 18.7
k,5,45.0  ,, 38.0000
k,6,67.0  ,, 49.2000
k,7,112.0  , , 66.5000
k,8,135.0  ,, 81.7000
k,9,145.0  ,, 84.2000
k,10,167.12 , , 92.0000
k,11 ,177.12 ,, 92.0000
k,12 ,177.12 ,,-10.0000
k,13,-10,6.8,-10
l,1,13
a,1,2,3,4,5,6,7,8,9,10,11,12
!生成斜坡模型（体1）
vdrag,1,,,,,,1

MSHAPE,1,3D
MSHKEY,0
et,1,solid95   !岩土单元
!定义参数
!岩土参数
mp,ex,1,6e9
mp,nuxy,1,0.3
mp,dens,1,2500  
tb,dp,1
tbdata, 1,700000,39
!定义单元类型
!划分网格
type,1
mat,1
vmesh,all

fini
!进入求解器
/solu
!声明分析类型
antype,STATIC
NLGEOM,ON
NROPT,FULL
SSTIF,ON
da,1,uy,,
da,14,uy,,
da,2,ux,,
da,12,ux,,
da,13,all

!工况1：自然状态下变形内力场计算
time,1
DELTIM,1,,1
autots,on
acel,,,10
LNSRCH,ON
ncnv,2
PRED,ON
!输出所有结果
outres,all
outpr,all
*do,j,5,6,1
i=5.2+0.02*j
tbdata, 1,700000/i,180*atan(tan(39*3.1415926/180)/i)/3.1415926
solve
*enddo