隧道开挖
麻烦高手们帮我看看，出现这种错误怎么回事？下面是命令流。*** ERROR *** CP = 180.562 TIME= 10:39:48
Element 22 requires material properties and references material number
3 for which no properties are defined. Checking terminates.
fini
/cle
/UNITS,SI
/prep7
et,1,solid45
et,2,mesh200,6
!et,3,solid65
et,3,shell63
r,1,0.35
!r,2,0.2
!!!!!!!!!!定义材料参数!!!!!!!!!!!
!干沙材料参数
mp,ex,1,3.e7 !土层
mp,prxy,1,0.30
mp,dens,1,1530
TB,DP,1,,,
TBMODIF,1,1,0
TBMODIF,1,2,30
TBMODIF,1,3,10
!开挖干沙
mp,ex,2,3.e7 !土层
mp,prxy,2,0.30
mp,dens,2,1530
TB,DP,2,,,
TBMODIF,1,1,0
TBMODIF,1,2,30
TBMODIF,1,3,10
!管片材料参数，c50,管片衬砌按各向同性计算
mp,ex,3,2.07e9 !考虑刚度折减3.45e9*0.6=2.07e9
mp,prxy,3,0.17
mp,dens,3,2450
!!!!!!!!!!建立平面内模型并划分单元!!!!!!!!!!!
!在隧道中心线定义局部坐标，便于后面的实体选取
local,11,0,0,0,0 !卡式坐标
local,12,1,0,0,0 !圆柱坐标
csys,11 !声明坐标系统卡式坐标
wpcsys,-1 !观察的方向垂直工作面
cyl4,,,,,2.5,90
arsym,y,all !选择所有面沿y轴做对称面
nummrg,all
numcmp,all
rectng,0,4,0,4
aovlap,all !作面递加
nummrg,all !合并所有元素
numcmp,all !压缩没有用的编号，并重新编号
rectng,4,17.5,0,4
nummrg,all
numcmp,all
arsym,y,all !选择所有面沿y轴做对称面
nummrg,all
numcmp,all
/pnum,line,1
/pnum,area,1
aplot,all
rectng,0,4.0,4,12.5
rectng,4,17.5,4,12.5
nummrg,all
numcmp,all
rectng,0,4,-4,-8.5
rectng,4,17.5,-4,-8.5
nummrg,all
numcmp,all
meshkey,1
type,2
lesize,1,,,6 !设定线中单元尺寸
lesize,2,,,5
lesize,3,,,5
lesize,4,,,6
lesize,5,,,5
lesize,8,,,3 !设定线中单元尺寸
lesize,9,,,3
lesize,15,,,3
lesize,19,,,8
lesize,7,,,8
lesize,6,,,8 !设定线中单元尺寸
lesize,13,,,8
lesize,14,,,8
lesize,23,,,8
lesize,25,,,6 !设定线中单元尺寸
lesize,24,,,6
lesize,27,,,6
lesize,16,,,8
lesize,11,,,8
lesize,20,,,8
lesize,18,,,8 !设定线中单元尺寸
lesize,21,,,8
lesize,22,,,10
lesize,12,,,10
lesize,10,,,10
lesize,17,,,10 !设定线中单元尺寸
lesize,26,,,10
amesh,all
save
!!!!!!!!!!建立立体模型!!!!!!!!!!!
!沿隧道轴线定义一系列关键点
!首先拉伸成壳单元
k,1000,0,0,-45
l,3,1000
/view,1,1,1,1
/replot
EXTOPT,ESIZE,10,0,
LSEL,S,line,,1,4,3
ADRAG,all,,,,,,28
gplot
type,3
real,1
mat,3
ASEL,S,loc,z,-22.5
APLOT
lsel,s,loc,z,-22.5
lesize,all,,,10
MSHAPE,0,2D
MSHKEY,1
amesh,all
!延伸土的实体
ASEL,invert
aplot
EXTOPT,ESIZE,10,0,
EXTOPT,ACLEAR,1
TYPE,1
MAT,1
asel,r,area,,3,10,1
VDRAG,all,,,,,,28
allsel
!延伸开挖土的实体
MAT,2
VDRAG,1,,,,,,28
EPLOT
nummrg,all
numcmp,all
!施加边界条件，并进行求解设置
fini
/solu
!约束两侧面的X方向的约束
asel,s,loc,x,17.5
da,all,ux,0
alls
!约束地面的Y方向的约束
asel,s,loc,y,-21
da,all,uy,0
NSEL,S,LOC,X,0
DSYM,SYMM,X,0,
alls
asel,s,loc,z,-45
asel,a,loc,z,0
da,all,uz,0
allsel
acel,,10
fini
/solu
antype,static
deltim,0.1,0.05,0.2
autots,on !使用自动时间步
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
nlgeom,on !打开大位移效果
nropt,full !设定牛顿-拉普森选项
cnvtol,f,,0.02,2,0.5
esel,s,type,,3 !选择壳单元，杀死
ekill,all
esel,all
esel,s,live !选择活的单元
nsle,s !选择活单元上的节点
nsel,invert !反向选择，即选择了死单元上的节点
d,all,all,0 !将死单元上的节点约束所有位移
nsel,all
esel,all
/PBC,ALL,,1
gplot
step=0 !用于记录计算步的参数
!求解自重应力场
time,1
solve
save 麻烦大侠帮我解决一下，我真的很着急，谢谢了 我也出现过这个问题啊，好像是单元没有被赋予材料参数属性，是不是单元网格划分有歧异网格啊，我后来重新划分了网格，没事了，本人也很困惑，希望高手解答呀 同样的问题 我觉得是单元类型的问题 期待解决—
页:
[1]