flac3d 隧道 蠕变 求指点

fangaochenabc

generate zone radcyl p0 0 0 0 p1 12 0 0 p2 0 1 0 p3 0 0 20 &
size 4 2 10 8 dim 6 6 6 6 rat 1 1 1 1.2  fill group inrockys
generate zone radcyl p0 0 0 0 p1 0 0 -18 p2 0 1 0 p3 12 0 0 &
size 4 2 10 8 dim 3 6 3 6 rat 1 1 1 1.2  fill group inrockyx
generate zone brick p0 12 0 0 p1 36 0 0 p2 12 1 0 p3 12 0 20 &
size 10 2 5
generate zone brick p0 12 0 0 p1 12 0 -18 p2 12 1 0 p3 36 0 0 &
size 5 2 10
gen zon reflect dip 90 dd 270 orig 0 0 0
inter 1 face range x 11.9 12.1 y 0 1 z -18 20
inter 2 face range x -12.1 y 0 1 z -18 20
set gravity 0 0 -9.8
fix z range z -18.01 -17.99
fix x range x 35.99 36.01
fix x range x -36.01 -35.99
fix y range y -0.01 0.01
fix y range y 0.99 1.01
model moh
ini density 2000
prop bulk=1.47e9 shear=5.6e8 fric=30 coh=4e5 tension 5e5

hist unbal
hist gp zdisp 0 0 6
hist gp xdisp 6 0 0
hist gp xdisp -6 0 0
plot add surface yellow
plot show
plot cont zdisp
solve
model null range group inrockys any
solve
sel shell id=1 range x -6 6 z 0 6 cyl end1 0 0 0 end2 0 1 0 rad 6.0
sel node fix x yr zr range x -0.1 0.1
sel node fix y xr zr range y -0.1 0.1
sel node fix z xr yr range z -0.1 0.1
sel shell id=1 prop iso=(10.5e9,0.25) thick=0.2
plot add sel geom shell on
solve
pause

config creep
def burcal
bkk=bme/(3*(1-2*b1))
bk1=bke/(2*(1+b1))
bk2=bme/(2*(1+b1))
bvis1=vis1
bvis2=vis2
end
set bme=6.35714e9     ;Em
set bke=3.25714e9     ;Ek
set b1=0.35        ;泊松比
set vis1=2.4e6
set vis2=4.5e6
burcal
model burger
prop  bulk bkk ksh bk1 msh bk2 kvis bvis1 mvis bvis2
set gravity 0 0 -9.8
ini xdis 0 ydis 0 zdis 0
set creep dt 5
set creep on
plot cont zdisp
hist id=98 unbal
hist id=99 crtime
hist dt
hist id=100 gp zdisp 0 0 6
hist id=101 gp xdisp 6 0 0
hist id=102 gp xdisp -6 0 0
solve age 86400


帮忙看下 有什么问题 我运行后 开挖后的隧道怎么又被恢复了，蠕变计算时不收敛，并且位移过大。

lg14507

model burger 这步应该设置范围，这个范围应该不包括已经开挖的部分吧。

fangaochenabc

开挖了以后不是不存在了吗，这样定义难道把已经开挖的又给加上了？

幻想成为学术姐

model burger 和prop 之后都得加range