求助！使用ansys进行路面结构分析出错

目的：将实际分层路面用ansys进行建模，模拟受力分析。将ansys计算的结果与实际测得数据进行比对。
过程：用命令流进行建模，但得出的数据不准确，且云图不连续，因自己使用ANSYS对路面分析不懂，实在找不出问题所在，切盼高手指点！命令流如下：
/filename,Pavement
/Title,Pavement.txt
/PREP7 !进入前处理模块PREP7
ET,1,PLANE42,,,1
MP,EX,7,1.5E9 !AC-13的材料特性
MP,NUXY,7,0.35
MP,DENS,7,2700
MP,EX,6,1.3E9 !AC-20C的材料特性
MP,NUXY,6,0.35
MP,DENS,6,2600
MP,EX,5,1.0E9 !AC-25C的材料特性
MP,NUXY,5,0.35
MP,DENS,5,2500
MP,EX,4,1.2E9 !ATB-30的材料特性
MP,NUXY,4,0.35
MP,DENS,4,2300
MP,EX,3,1.5E9 !水泥碎石的材料特性
MP,NUXY,3,0.3
MP,DENS,3,2300
MP,EX,2,1.2E9 !水泥砂砾的材料特性
MP,NUXY,2,0.35
MP,DENS,2,2200
MP,EX,1,2.675E9 !土的材料特性
MP,NUXY,1,0.3375
MP,DENS,1,1000
!定义节点
n,100,0,0
n,1000,0,0.5
n,2000,0,1
n,3000,0,1.5
n,4000,0,2
n,5000,0,2.5
n,200,0,3
n,300,0,3.15
n,400,0,3.46
n,500,0,3,56
n,600,0,3.64
n,700,0,3.70
n,800,0,3.74

*Do,i,100,119,1 !水平方向结点的循环
n,i,(5/20)*(i-100),0
*enddo
*Do,i,1000,1019,1 !水平方向结点的循环
n,i,(5/20)*(i-1000),0.5
*enddo
*Do,i,2000,2019,1 !水平方向结点的循环
n,i,(5/20)*(i-2000),1
*enddo
*Do,i,3000,3019,1 !水平方向结点的循环
n,i,(5/20)*(i-3000),1.5
*enddo
*Do,i,4000,4019,1 !水平方向结点的循环
n,i,(5/20)*(i-4000),2
*enddo
*Do,i,5000,5019,1 !水平方向结点的循环
n,i,(5/20)*(i-5000),2.5
*enddo
*Do,i,200,219,1 !水平方向结点的循环
n,i,(5/20)*(i-200),3
*enddo
*Do,i,300,319,1 !水平方向结点的循环
n,i,(5/20)*(i-300),3.15
*enddo
*Do,i,400,419,1 !水平方向结点的循环
n,i,(5/20)*(i-400),3.46
*enddo
*Do,i,500,519,1 !水平方向结点的循环
n,i,(5/20)*(i-500),3.56
*enddo
*Do,i,600,619,1 !水平方向结点的循环
n,i,(5/20)*(i-600),3.64
*enddo
*Do,i,700,719,1 !水平方向结点的循环
n,i,(5/20)*(i-700),3.70
*enddo
*Do,i,800,819,1 !水平方向结点的循环
n,i,(5/20)*(i-800),3.74
*enddo
!定义单元
TYPE,1
MAT,1
*do,i,100,118,1
e,i,i+900,i+901,i+1
*enddo
TYPE,1
MAT,1
*do,i,1000,1018,1
e,i,i+1000,i+1001,i+1
*enddo
TYPE,1
MAT,1
*do,i,2000,2018,1
e,i,i+1000,i+1001,i+1
*enddo
TYPE,1
MAT,1
*do,i,3000,3018,1
e,i,i+1000,i+1001,i+1
*enddo
TYPE,1
MAT,1
*do,i,4000,4018,1
e,i,i+1000,i+1001,i+1
*enddo
TYPE,1
MAT,1
*do,i,5000,5018,1
e,i,i-4800,i-4799,i+1
*enddo
type,1
mat,2
*do,i,200,218,1
e,i,i+100,i+101,i+1
*enddo
type,1
mat,3
*do,i,300,318,1
e,i,i+100,i+101,i+1
*enddo
type,1
mat,4
*do,i,400,418,1
e,i,i+100,i+101,i+1
*enddo
type,1
mat,5
*do,i,500,518,1
e,i,i+100,i+101,i+1
*enddo
type,1
mat,6
*do,i,600,618,1
e,i,i+100,i+101,i+1
*enddo
type,1
mat,7
*do,i,700,718,1
e,i,i+100,i+101,i+1
*enddo
!定义约束
*do,i,100,119,1
d,i,uy
*enddo
F,806,FY,-1,,815,1 !定义节点间分布力
/SOLU !进入求解模块SOLUTION
SOLVE !开始求解
/POST1 !进入后处理模块POST1
SET,1 !读入第一个荷载步的计算结果
pldisp,2 !显示变形图(保留原结构轮廓)
PLNSOL,S,y,0,1 !在图形窗口显示结点上的X方向应力
PLNSOL,S,1,0,1 !在图形窗口显示结点上的主应力
PLNSOL,S,EQV,0,1 !在图形窗口显示结点Mises应力
FINISH !退出后模块