材料属性？？

在执行 nummrg,all 和 numcmp,all 两个语句后，材料 1 与材料 2 合并，而材料 3～5 的编号压缩为 2～4，但是后面又使用了语句 mat,5 - 指定使用 (已不存在的) 材料编号 5，由此造成了相应的错误。
可能的解决办法有多种，可以试试看：
1 因为材料 5 的编号已经改为 4 ，故应将 mat,5 改为 mat,4；同时还应该查一下其余的 mat,n 语句，将其中的材料编号 n 改为正确的值。
2 求解前再次定义材料 5，保证没有不存在的材料编号，不过对于 mat,n 语句仍需要进行检查。
3 执行完 nummrg 语句后，不执行 numcmp,all，即保持未合并的材料编号不变。
建议按第一种办法进行修改。

你好，我也遇到这样的问题，我按照论坛高手所说的把nummrg,all命令中的all改为node，这样的结果就更糟糕了。出现的错误更多了，不知道能不能帮忙解决，我的命令流是：
谢谢，不胜感激
!材料、实常数和单元类型定义
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,type,,1 !选择管片的单元
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 !求解

以后多宣传一下咱们的论坛家园

o(∩_∩)o...，非常感谢，你太厉害了

/prep7
et,1,solid65
mp,ex,1,3.3e10
mp,prxy,1,0.17
mp,ex,2,2.1e11
mp,prxy,2,0.333
mp,ex,3,4.99e10
mp,prxy,3,0.186
mp,ex,4,5.22e10
mp,prxy,4,0.188
mp,ex,5,4.3e10
mp,prxy,5,0.179
mp,ex,6,4.23e10
mp,prxy,6,0.179
mp,ex,7,7e8
mp,nuxy,7,0.27

mp,ex,8,6e8
mp,nuxy,8,0.29
mp,ex,9,5.5e8
mp,nuxy,9,0.3
mp,ex,10,5.2e8
mp,nuxy,10,0.30
mp,ex,11,0
mp,nuxy,11,0.31
mp,ex,12,6e8
mp,nuxy,12,0.29
mp,ex,13,10e8
mp,nuxy,13,0.26
mp,ex,14,5.5e8
mp,nuxy,14,0.3
r,1
r,2,0,0,0.0001,2e6,0.0002,2e6
r,3,0,0,0.0001,2e5,0.0002,2e5
r,4,0,0,0.0001,2e3,0.0002,2e3
r,5,0,0,0.0001,1e2,0.0002,1e2
r,6,0,0,0.0001,2e4,0,0002,2e4
r,7,0,0,0.0001,2e2,0,0002,2e2
r,8,2,0.10,0,90
r,9,2,0.11,0,90
r,10,2,0.05,0,90
r,11,2,0.05,0,90
r,12
r,13
r,14
blc4,0,0,3.36,2.70,34
cyl4,0,0,0,0,1.425,90,34
VSBV,1,2
vgen,,3,,,,,63,,,1
numcmp,all
blc4,0,0,3.36,2.70,5.6
cyl4,0,0,0,0,1.425,90,5.6
VSBV,2,3
vgen,,4,,,,,111.4,,,1
numcmp,all
blc4,0,0,3.36,2.70,116
cyl4,0,0,0,0,1.25,90,116
VSBV,3,4
vgen,,5,,,,,-53,,,1
numcmp,all
cyl4,0,0,1.32,0,1.4,90,54
cyl4,0,0,1.25,0,1.32,90,54
cyl4,0,0,1.20,0,1.25,90,54
cyl4,0,0,1.17,0,1.20,90,54
cyl4,0,0,1.11,0,1.17,90,54
cyl4,0,0,1.05,0,1.11,90,54
cyl4,0,0,0,0,1.05,90,54
vgen,,4,10,,,,63,,,1
cyl4,0,0,0,0,1.05,90,10
cyl4,0,0,1.05,0,1.11,90,10
cyl4,0,0,1.11,0,1.17,90,10
cyl4,0,0,1.17,0,1.25,90,10
vgen,,11,14,,,,53,,,1
cyl4,0,0,0,0,1.05,90,20
cyl4,0,0,1.05,0,1.11,90,20
cyl4,0,0,1.11,0,1.17,90,20
cyl4,0,0,1.17,0,1.25,90,20
vgen,,15,18,,,,33,,,1
cyl4,0,0,0,0,1.05,90,33
cyl4,0,0,1.05,0,1.11,90,33
cyl4,0,0,1.11,0,1.17,90,33
cyl4,0,0,1.17,0,1.25,90,33
nummrg,node
numcmp,all

问题出在nummrg,all ，这个命令把所有的项目都压缩合并，包括材料号和实常数号。如过想压缩节点可以：NUMMRG, node
NUMMRG, Label, TOLER, GTOLER, Action, Switch
Merges coincident or equivalently defined items.
PREP7: Database
MP ME ST PR PRN <> <> FL EM <> <> PP <>
Label
Items to be merged:
NODE — Nodes

ELEM — Elements

KP — Keypoints (will also merge lines, areas, and volumes)

MAT — Materials

TYPE — Element types

REAL — Real constants

CP — Coupled sets

CE — Constraint equations

ALL — All items

/prep7
et,1,solid65
mp,ex,1,3.3e10
mp,prxy,1,0.17
mp,ex,2,2.1e11
mp,prxy,2,0.333
mp,ex,3,4.99e10
mp,prxy,3,0.186
mp,ex,4,5.22e10
mp,prxy,4,0.188
mp,ex,5,4.3e10
mp,prxy,5,0.179
mp,ex,6,4.23e10
mp,prxy,6,0.179
mp,ex,7,7e8
mp,nuxy,7,0.27

mp,ex,8,6e8
mp,nuxy,8,0.29
mp,ex,9,5.5e8
mp,nuxy,9,0.3
mp,ex,10,5.2e8
mp,nuxy,10,0.30
mp,nuxy,11,0.31
mp,ex,12,6e8
mp,nuxy,12,0.29
mp,ex,13,10e8
mp,nuxy,13,0.26
mp,ex,14,5.5e8
mp,nuxy,14,0.3
r,1
r,2,0,0,0.0001,2e6,0.0002,2e6
r,3,0,0,0.0001,2e5,0.0002,2e5
r,4,0,0,0.0001,2e3,0.0002,2e3
r,5,0,0,0.0001,1e2,0.0002,1e2
r,6,0,0,0.0001,2e4,0,0002,2e4
r,7,0,0,0.0001,2e2,0,0002,2e2
r,8,2,0.10,0,90
r,9,2,0.11,0,90
r,10,2,0.05,0,90
r,11,2,0.05,0,90
r,12
r,13
r,14
blc4,0,0,3.36,2.70,34
cyl4,0,0,0,0,1.425,90,34
VSBV,1,2
vgen,,3,,,,,63,,,1
numcmp,all
blc4,0,0,3.36,2.70,5.6
cyl4,0,0,0,0,1.425,90,5.6
VSBV,2,3
vgen,,4,,,,,111.4,,,1
numcmp,all
blc4,0,0,3.36,2.70,116
cyl4,0,0,0,0,1.25,90,116
VSBV,3,4
vgen,,5,,,,,-53,,,1
numcmp,all
cyl4,0,0,1.32,0,1.4,90,54
cyl4,0,0,1.25,0,1.32,90,54
cyl4,0,0,1.20,0,1.25,90,54
cyl4,0,0,1.17,0,1.20,90,54
cyl4,0,0,1.11,0,1.17,90,54
cyl4,0,0,1.05,0,1.11,90,54
cyl4,0,0,0,0,1.05,90,54
vgen,,4,10,,,,63,,,1
cyl4,0,0,0,0,1.05,90,10
cyl4,0,0,1.05,0,1.11,90,10
cyl4,0,0,1.11,0,1.17,90,10
cyl4,0,0,1.17,0,1.25,90,10
vgen,,11,14,,,,53,,,1
cyl4,0,0,0,0,1.05,90,20
cyl4,0,0,1.05,0,1.11,90,20
cyl4,0,0,1.11,0,1.17,90,20
cyl4,0,0,1.17,0,1.25,90,20
vgen,,15,18,,,,33,,,1
cyl4,0,0,0,0,1.05,90,33
cyl4,0,0,1.05,0,1.11,90,33
cyl4,0,0,1.11,0,1.17,90,33
cyl4,0,0,1.17,0,1.25,90,33
nummrg,all
numcmp,all
再查看的时候13，14号材料已经没有了，实常数11，12，13，14也已经没了，请问是怎么回事？