恳请高手指点

小弟正在研究盾构施工地表沉降数值模拟,在我建立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 !求解
邀请回答 举报

共1个回答

默认排序
  • 邹正刚
    2009年5月24日
    在执行 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,即保持未合并的材料编号不变。
    建议按第一种办法进行修改。检查一下所有 mat,n 语句,看材料号 n 是否正确;至少有一句 mat,5 是不对的,应该改为 mat,4;其余你自己检查一下。
    0评论 采纳答案 举报

  • 恳请高手指点的相关案例教程

    点击蓝字 关注我们 好久没有更新算例了,这期更新一个小算例,用强度折减法计算边坡的稳定性,用ANSYS软件计算,ABAQUS的算例我前面写了两篇,感兴趣的可以看看。 边坡的尺寸及参数与之前的ABAQUS的边坡算例相同。模型如下。 接下来讲解一下建模的过程。先定义一下单元的类型。这里采用平面应变单元,PLANE42单元。 !定义单元类型ET,1,PLANE82     !定义边坡围岩单元KEYOPT
    本文以二维静态磁场为例,介绍一下使用远场单元注意事项,并给出一个简单的APDL算例,软件版本ANSYS19.0。 一、问题介绍及注意事项 对于ANSYS二维静态磁场分析,磁力线总是平行或垂直于边界的,有时与实际情况是不符的,这时候就要引入infin110等远场单元。下面就以infin110为例,列出几项远场单元的注意事项: 1)infin110只需在最外层划分一层单元。 2)最外层的infin11
    !ANSYS命令流学习笔记6 !热应力分析 !学习重点: !1、 理解热力耦合的直接法和间接法 !间接法:先进行热分析,然后将求得的节点温度作为载荷施加到结构应力分析中。 !直接法:直接采用具有温度和位移自由度的耦合单元,同时得到热分析和结构应力分析的结果。直接法又分弱耦合和强耦合选择强耦合时,形成不对称矩阵,线性系统可以直接求解。选择弱耦合时,对称矩阵,还是把热和结构分别进行求解,并将热结果施加
    !多工况下的拓扑优化 !学习重点: !1、 何为拓扑优化 !区分尺寸优化、形状优化、拓扑优化。拓扑优化是形状优化的一种特殊形式。网上资料和ansys help文件都有详细说明。 !2、单一工况载荷下的拓扑优化。单一工况拓扑操作流程很简单,APDL命令也很简单。Workbench也可以轻松实现。 !3、多工况下的拓扑优化。需要apdl用到lswrite 命令。目前不清楚如何完全由workbench完
    背景说明:两个半径均为100mm,的正交圆柱体发生正接触,作用在两圆柱接触体法线方向上的压力总和1000N,两圆柱体均为钢制分析两圆柱的接触情况。 分析思路:对模型进行对称约束,施加载荷后要进行节点自由度耦合,使加载面的所有节点能在承载后具有相同的Y方向位移,并且在变形后仍然保持为水平平面。 总结: 1. 如果不进行节点耦合,结构变形会产生刚性位移,造成不收敛。但是节点耦合造成了变形不符合实际情况