恳请高手指点

小弟正在研究盾构施工地表沉降数值模拟，在我建立ANSYS模型处理的时候，遇到了这样的问题： Element 65 requires material properties and references material number 5 for which no properties are defined. Checking terminates.
我按照论坛上的一些高手改了nummrg,all 命令中all，但是还是不行，这样改的话问题更多，更不能处理了，请问哪位高手能给予帮助，小弟不胜感激。希望多多帮忙啊。下面是小弟的命令流，希望高手帮忙改改。后处理部分也是不对的。
!材料、实常数和单元类型定义
Fini
/clear !更新数据库
/filename,exm
/title,the analysis of the arch tunnel
*set,h,12 !隧道埋深
*set,D,6 !隧道直径
*set,th,0.35 !管片厚度
*set,m,-1.5*D !围岩顶面坐标
/prep7 !进入前处理器
et,1,plane42 !设置单元类型
keyopt,1,3,2 !定义为 平面应变问题
r,1,th,th*th*th/12,th,, !定义实常数
save !保持数据
!定义模型中的材料参数
!被开挖土体材料参数
mp,ex,1,5e6 !土体弹性模量
mp,prxy,1,0.35 !土体泊松比
mp,dens,1,1830
tb,dp,1 !定义DP材料
tbdata,1, 2.0e3,35 !定义DP材料参数
!隧道上方土体材料参数
mp,ex,2,5e6 !土体弹性模量
mp,prxy,2,0.35 !土体泊松比
mp,dens,2,1830 !土体密度
tb,dp,2 !定义DP材料
tbdata,1, 2.0e3,35 !定义DP材料参数
!基岩材料参数
Mp,ex,3,5.0e8
Mp,prxy,3,0.33
Mp,dens,3,2160
tb,dp,3 !定义DP材料
tbdata,1, 2.7e6,33 !定义DP材料参数
!管片材料参数，管片衬砌按各向同性计算
mp,ex,4,2.07e9 !管片衬砌弹性模量
mp,prxy,4,0.2 !管片衬砌泊松比
mp,dens,4,2500 !管片衬砌密度
!等代层，参数按水泥土取值
mp,ex,5,1.2e9 !等代层弹性模量
mp,prxy,5,0.2 !等代层泊松比
mp,dens,5,2100 !等代层密度
tb,dp,4 !定义DP材料
tbdata,1, 1.8e3,30 !定义DP材料参数
save !保持数据
!********************************************
!建立平面内模型并划分单元
!在隧道中心线定义局部坐标，便于后来的实体选取。
local,11,0,0,0,0 !局部笛卡儿坐标
local,12,1,0,0,0 !局部极坐标
csys,11 !将当前坐标转换为局部坐标
wpcsys,-1 !同时将工作平面转换到局坐标
cyl4,,,,,0.5*D,90 !画隧道内经为D/2
cyl4,0,0,0.5*D,0,0.5*D+0.35,90 !画隧道内经为D/2
cyl4,0,0,0.5*D+0.35,0,0.5*D+0.35+0.2,90 !画管片
rectng,0,1.5*D,0,1.5*D !画外边界矩形
aovlap,all !做面递加
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
rectng,1.5*D,5.5*D,0,1.5*D !画矩形面
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
save !保持数据
!划分被开挖土体单元
meshkey,1 !选择划分方式为映射划分
type,1 !采用plane42 进行平面内的单元划分
mat,1 !指定为1号材料号
lesize,1,,,8 !对线设置单元数
lesize,2,,,8
lesize,3,,,8
amesh,1 !对面1划分单元
!划分管片的有限元网格
type,1 !采用plane42 进行平面内的单元划分
mat,4 !指定为4号材料号
lesize,4,,,8 !对线设置单元数
lesize,8,,,2
lesize,9,,,2
amesh,2 !对面2划分单元
save
!划分等代层的有限元网格
type,1 !采用plane42 进行平面内的单元划分
mat,5 !指定为3号材料号
lesize,5,,,8 !对线设置单元数
lesize,10,,,1
lesize,11,,,1
amesh,3 !对面2划分单元
save
!划分隧道右边土体有限元网格
type,1 !采用plane42 进行平面内的单元划分
mat,2 !指定为2号材料号
lesize,12,,,4 !对线设置单元数
lesize,13,,,4
lesize,6,,,4
lesize,7,,,4
lesize,14,,,8,2
lesize,16,,,8,0.5
amap,4,7,6,8,10
amap,5,9,8,11,12
save !保持数据
!利用对称性划分单元得到下半部分模型
arsym,y,all !以y轴为对称轴进行镜像操作
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
allsel,all !选择所有元素
save !保持数据
!建立隧道下方地层模型
rectng,0,1.5*D,-1.5*D,-5.5*D !绘制下边界矩形面
rectng,1.5*D,5.5*D,-1.5*D,-5.5*D
save
!**************************************
!划分隧道下方围岩的有限元网格
type,1 !采用plane42 进行平面内的单元划分
mat,3
lesize,28,,,4 !对线设置单元数
lesize,29,,,8,0.5
lesize,30,,,4
lesize,31,,,8,2
lesize,32,,,8,2
lesize,33,,,8,0.5
lesize,34,,,8,0.5
lesize,35,,,8,2
asel,s,area,,11,12 !对面11和12划分单元
amesh,all
save !保持数据
!建立隧道上方土层模型
rectng,0,1.5*D,1.5*D,h+0.5*D !绘制上边界矩形面
rectng,1.5*D,5.5*D,1.5*D,h+0.5*D
type,1 !采用plane42 进行平面内的单元划分
mat,2 !指定为2号材料号
lesize,36,,,4 !对线设置单元数
lesize,37,,,6
lesize,38,,,4
lesize,39,,,6
lesize,40,,,8,2
lesize,41,,,6
lesize,42,,,8,0.5
lesize,43,,,6
amesh,13 !对面13和14划分单元
amesh,14
nummrg,all !合并所有元素
numcmp,all !压缩所有元素编号
save !保持数据
!***********************************
!施加边界条件
!施加约束条件
nsel,s,loc,x,5.5*D !选择节点
d,all,Ux,0 !施加X方向约束
nsel,s,loc,x,0
d,all,Ux,0 !施加X方向约束
nsel,s,loc,y,-5.5*D !选择节点
d,all,uy !施加Y方向约束
allsel,all !选中所有元素
save
!**********************************
!施加荷载
acel,0,10,0 !施加重力
lsel,s,loc,y,h+0.5*D !选择地面所有节点
sfl,all,pres,20e3 !施加地面荷载
save
fini
/solu !进入求解器
antype,static !设置分析类型为静态分析
autots,on !使用自动时间步
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
nlgeom,on !打开大变形求解
nropt,full !设定牛顿-拉普森选项
outres,all,all !输出所有项，每一步都输出
!初始地应力的计算
esel,s,mat,,4 !选择管片的单元
ekill,all !杀死单元即在隧道未修建前的自重应力场中不存在壳单元
esel,all !选择所有单元
esel,s,live !选择所有活单元
nsle,s !选择当前单元中所有节点
nsel,inve !反向选择
d,all,all !约束选择集所有节点、所有自由度
allsel !选择所有元素
solve !求解
esel,s,mat,,1,5,4 !选择挖去部分土体和等代层单元
ekill,all !杀死单元
esel,s,mat,,4 !选择管片的单元
ealive,all !激活单元即考虑隧道开挖地层损失对地表沉降的影响
nsle,s !选择当前单元中所有节点
ddele,all,all !将位移约束删除
esel,all !选择所有单元
esel,s,live !选择所有活单元
nsle,s !选择当前单元中所有节点
nsel,inve !反向选择
d,all,all,0 !将非活动节点约束
nsel,all !选择所有节点
esel,all !选择所有单元
solve !求解