各位大侠,我正在做一个煤矿马头门的模型,模型已经建好了,网格采用自由划分,只考虑开挖对围岩的应力和应变的影响,不考虑初期支护,在模拟开挖时却出问题了,请各位看看是什么原因,谢谢了!
! EKILL is not a recognized POST1 command,abbreviation,or macro.this command will be ignored.
命令流:
/prep7
!定义单元类型、实常数、材料属性
et,1,mesh200,2 !3-D线单元2节点
et,2,mesh200,6 !3-D面单元4节点
et,3,solid45
mp,ex,1,11.3e9 !剩余土体材料属性
mp,prxy,1,0.2
mp,dens,1,2500
mp,ex,2,11.3e9 !挖去5-5土体材料属性
mp,prxy,2,0.22
mp,dens,2,2550
mp,ex,3,11.3e9 !挖去3-3土体材料属性
mp,prxy,3,0.22
mp,dens,3,2550
mp,ex,4,11.3e9 !挖去立井土体材料属性
mp,prxy,4,0.22
mp,dens,4,2550
!5-5
type,3
mat,2
CYL4,0,0,3.35,0, ,180,21
BLOCK,-3.35,3.35,0,-5.3,0,21,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VADD,P51X
!右边3-3
type,3
mat,3
CYL4,0,-3.3,3.3,0, ,180,-26
BLOCK,-3.3,3.3,-3.3,-5.81,0,-26,
GPLOT
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VADD,P51X
!左边3-3
type,3
mat,3
FLST,2,1,8
FITEM,2,0,-3.3,21
WPAVE,P51X
CYL4,0,0,3.3,0, ,180,10.1
BLOCK,-3.3,3.3,0,-2.51,0,10.1,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VADD,P51X
!体粘结
FLST,2,2,6,ORDE,2
FITEM,2,3
FITEM,2,5
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,4
VGLUE,P51X
!立井
type,3
mat,4
FLST,2,1,8
FITEM,2,0,17.3,10.5
WPAVE,P51X
wpro,,-30.000000,
wpro,,-30.000000,
wpro,,-30.000000,
CYL4,0,0,4, , , ,-34.6
!创建围岩
type,3
mat,1
FLST,2,1,8
FITEM,2,0,0,0
WPAVE,P51X
wpro,,30.000000,
wpro,,30.000000,
wpro,,30.000000,
BLOCK,-16.75,16.75,-17.3,17.3,-26,31.1,
!体分割
FLST,2,5,6,ORDE,2
FITEM,2,1
FITEM,2,-5
VPTN,P51X
!自由划分网格
FLST,5,1,4,ORDE,1
FITEM,5,79
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,11, , , , ,1
FLST,5,3,4,ORDE,3
FITEM,5,77
FITEM,5,-78
FITEM,5,80
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,11, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,70
FITEM,5,72
FITEM,5,-73
FITEM,5,75
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,11, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,69
FITEM,5,71
FITEM,5,74
FITEM,5,76
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,12, , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,67
FITEM,5,-68
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,12, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,3
FITEM,5,-4
FITEM,5,83
FITEM,5,-84
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,93
FITEM,5,-94
FITEM,5,97
FITEM,5,-98
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,4, , , , ,1
FLST,5,2,4,ORDE,2
FITEM,5,81
FITEM,5,-82
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,2, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,17
FITEM,5,20
FITEM,5,34
FITEM,5,-35
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,13, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,33
FITEM,5,36
FITEM,5,50
FITEM,5,-51
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
FLST,5,3,4,ORDE,3
FITEM,5,11
FITEM,5,14
FITEM,5,27
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
FLST,5,3,4,ORDE,3
FITEM,5,28
FITEM,5,37
FITEM,5,42
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
FLST,5,8,4,ORDE,8
FITEM,5,38
FITEM,5,40
FITEM,5,48
FITEM,5,-49
FITEM,5,55
FITEM,5,57
FITEM,5,59
FITEM,5,-60
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1
FLST,5,8,4,ORDE,8
FITEM,5,9
FITEM,5,16
FITEM,5,21
FITEM,5,29
FITEM,5,-30
FITEM,5,32
FITEM,5,45
FITEM,5,-46
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,3, , , , ,1
FLST,5,8,4,ORDE,6
FITEM,5,12
FITEM,5,-13
FITEM,5,39
FITEM,5,41
FITEM,5,63
FITEM,5,-66
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,13, , , , ,1
FLST,5,9,4,ORDE,8
FITEM,5,22
FITEM,5,-24
FITEM,5,37
FITEM,5,48
FITEM,5,85
FITEM,5,-86
FITEM,5,88
FITEM,5,-89
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,3, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,95
FITEM,5,-96
FITEM,5,99
FITEM,5,-100
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1
FLST,5,4,4,ORDE,3
FITEM,5,87
FITEM,5,90
FITEM,5,-92
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,13, , , , ,1
MSHAPE,1,3D
MSHKEY,0
FLST,5,7,6,ORDE,4
FITEM,5,1
FITEM,5,3
FITEM,5,6
FITEM,5,-10
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!施加边界条件,并进行求解设置
/solu
antype,static
!约束两侧面X方向的约束
asel,s,loc,x,-16.75
asel,a,loc,x,16.75
da,all,ux,0
alls
!约束地面Y方向的约束
asel,s,loc,Y,-17.3
da,all,uy,0
alls
!约束岩石前后Z方向的约束
asel,s,loc,z,-26
asel,a,loc,z,31.1
da,all,uz,0
allsel
acel,,10,0, !施加重力加速度
!设定分析选项
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,all
esel,s,live ! 选择活的单元,即所有土体单元
nsel,all
esel,all
/PBC,ALL,,1
gplot
solve
!立井上部开挖模拟(24米,每天2米)
*do,ii,1,12,1
!以下步骤表示,先选择每天挖去的岩石单元为有效单元,然后将其杀死
esel,s,mat,,4
nsle,s
nsel,r,loc,y,17.3-(ii-1)*2,15.3-(ii-1)*2
esln,r,1
ekill,all
!选择活单元,此时指未挖去的岩石单元
esel,all
esel,s,live
nsle,s
nsel,invert !反向选择,将死单元上的节点约束所有自由度
d,all,all,
nsel,all
esel,all
solve
*enddo
!左右两边马头门同时开挖模拟(10.5米,每天2米)
*do,ii,1,6,1
!以下步骤表示,先选择每天挖去的岩石单元为有效单元,然后将其杀死
esel,s,mat,,2
nsle,s
nsel,r,loc,z,10.5+(ii-1)*2,12.5+(ii-1)*2
nsel,r,loc,z,10.5-(ii-1)*2,8.5-(ii-1)*2
esln,r,1
ekill,all
!选择活单元,此时指未挖去的岩石单元
esel,all
esel,s,live
nsle,s
nsel,invert !反向选择,将死单元上的节点约束所有自由度
d,all,all,
nsel,all
esel,all
solve
*enddo
!左右两边3-3同时开挖模拟(26米,每天2米)
*do,ii,1,13,1
!以下步骤表示,先选择每天挖去的岩石单元为有效单元,然后将其杀死
esel,s,mat,,3
nsle,s
nsel,r,loc,z,-(ii-1)*2,-(2+(ii-1)*2)
nsel,r,loc,z,21+(ii-1)*2,23+(ii-1)*2
esln,r,1
ekill,all
!选择活单元,此时指未挖去的岩石单元
esel,all
esel,s,live
nsle,s
nsel,invert !反向选择,将死单元上的节点约束所有自由度
d,all,all,
nsel,all
esel,all
solve
*enddo
!立井下部开挖模拟(10.6米,每天2米)
*do,ii,1,6,1
!以下步骤表示,先选择每天挖去的岩石单元为有效单元,然后将其杀死
esel,s,mat,,4
nsle,s
nsel,r,loc,y,-6.7-(ii-1)*2,-8.7-(ii-1)*2
esln,r,1
ekill,all
!选择活单元,此时指未挖去的岩石单元
esel,all
esel,s,live
nsle,s
nsel,invert !反向选择,将死单元上的节点约束所有自由度
d,all,all,
nsel,all
esel,all
solve
*enddo
finish |