ansys竖井开挖命令流.doc

本人在做一个地铁竖井施工的模型。计算过程中发现开挖后,应力释放的结果不理想,请各位帮帮忙!!
附:命令流:
finish
/clear,start
/filname,process of shujing
/title,process of shujing
/prep7
r1=       !竖井半径(井壁内径)
r2=r1+   !二支半径(井壁外径),
r3=r2+   !初支半径,
r4=r2+    !围岩中锚杆区的半径,(包括初支和锚杆区)
r5=12*r1     !边界范围半径
h=50         !基岩上部竖井深度
et,1,solid45
et,2,154      !定义三维表面效应单元  (井壁摩擦)
KEYOPT,2,2,0    !指定表面效应单元的K2=0,所加荷载与单元坐标系方向相同
KEYOPT,2,4,1    !指定表面效应单元的K4=0,去掉边中点,成为四结点表面单元
!围岩参数
ex1=1500e6   !弹性模量(Pa)
nuxy1=    !泊松比
dens1=1900   !密度(kg*m3)
c1=100e3     !粘聚力(Pa)
jiaodu1=25   !内摩擦家(°度)
mp,ex,1,ex1
mp,nuxy,1,nuxy1
mp,dens,1,dens1
tb,dp,1,,,0     !德鲁克——普拉格准则
tbdata,1,c1,jiaodu1,jiaodu1
!锚杆区
ex2=1600e6
nuxy2=
dens2=1900
c2=120e3
jiaodu2=26
mp,ex,2,ex2
mp,nuxy,2,nuxy2
mp,dens,2,dens2
tb,dp,2,,,0
tbdata,1,c2,jiaodu2,jiaodu2
!初支
ex3=2500e6
nuxy3=
dens3=2245
mp,ex,3,ex3
mp,nuxy,3,nuxy3
mp,dens,3,dens3
!井壁
ex4=2800e6
nuxy4=
dens4=2500
mp,ex,4,ex4
mp,nuxy,4,nuxy4
mp,dens,4,dens4
/view,1,1,1,1
!建立模型
cyl4,0,0,0,0,r1,90,10
cyl4,0,0,0,0,r2,90,10
cyl4,0,0,0,0,r3,90,10
cyl4,0,0,0,0,r4,90,10
vptn,all
vgen,5,all,,,,,10,,,
cyl4,0,0,0,0,r5,90,50
vptn,all
nummrg,all
numcmp,all
allsel
!划分网格
local,11,1,,,          !摩擦力
csys,11
asel,s,loc,x,r2
aatt,4,,2,0
lesize,70,,,10  !圆弧
lesize,96,,,4
lesize,22,,,4
lesize,99,,,4
lesize,111,,,4
lesize,