ANSYS边坡稳定性分析求助？

/CLEAR FINISH !定义工作名称和分析标题 /FILNAME,Slope,1 /TITLE,Slope stability Analysis !定义单元类型 /prep7 ET,1,PLANE183 KEYOPT,1,3,2 !平面应变选项 !定义模型中的材料参数 MP,EX,1,31E9 !围岩1弹性模量 MP,PRXY,1,0.24 !围岩1泊松比 MP,DENS,1,2700 !围岩1密度 *DO,I,2,13 XS=1.0+(I-2)*0.2 MP,EX,I,30E9 !围岩2弹性模量 MP,PRXY,I,0.25 !围岩2泊松比 MP,DENS,I,2500 !围岩2密度 TB,DP,I TBMODIF,1,1,0.9E6/XS !围岩2折减XS后的内聚力 TBMODIF,1,2,42/XS !围岩2折减XS后的摩擦角 *ENDDO SAVE !存储数据库 !建立几何模型并划分单元 K,1, !根据坐标生成关键点 K,2,-800 K,3,-800,-800 K,4,-800,-1200 K,5,1200,-1200 K,6,1200,-800 k,7,1200 K,8,1200,378 K,9,430,378 SAVE *DO,I,1,8 !由关键点创建线 L,I,I+1 *ENDDO L,9,1 L,1,7 L,3,6 SAVE AL,3,4,5,11 AL,1,2,11,6,10 AL,7,8,9,10 SAVE !划分网格 LSEL,S,,,3,5,2 LESIZE,ALL,,,5 LSEL,S,,,4,11,7 LESIZE,ALL,,,26 MAT,1 TYPE,1 MSHAPE,0 !网格形状为四边形，1为三角形网格 MSHKEY,0 !网格划分采用扫略方式；0为自由划分；2为以扫略为主，自由为辅 AMESH,1 LSEL,S,,,2,6,4 LSEL,A,,,1 LESIZE,ALL,,,10 LSEL,S,,,8,10,2 LESIZE,ALL,,,16 LSEL,S,,,7,9,2 LESIZE,ALL,,,12 MAT,2 TYPE,1 MSHAPE,0 MSHKEY,0 AMESH,2 AMESH,3 ALLSEL SAVE,Slope_Grid.db !指定边界条件 NSEL,S,LOC,X,-800 NSEL,A,LOC,X,1200 D,ALL,UX ALLSEL NSEL,S,LOC,Y,-1200 D,ALL,UX D,ALL,UY ACEL,9.8 SAVE !求解设置 /SOLU ANTYPE,STATIC !静力求解 NSUBST,5,100,1 !荷载步5，最大荷载子步数100，最小荷载子步数1 PRED,ON !打开时间步长预测器 NROPT,FULL !设定牛顿-拉普森选项 NLGEOM,ON !打开大变形 LNSRCH,ON !打开线性搜索 OUTRES,ALL,ALL !输出所有项 CNVTOL,F,,0.005,2,0.5 !力收敛准则 CNVTOL,U,,0.005,2,1 !位移收敛准则 !边坡在强度折减系数XS=1.0时求解 ALLSEL SOLVE SAVE,F1.0,DB