- ;1/2建模隧道部分
- new
- gen zone radcy p0 9 0 0 p1 18 0 0 p2 9 36 0 p3 9 0 9 size 4 28 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outr
- gen zone cshell p0 9 0 0 p1 12 0 0 p2 9 36 0 p3 9 0 3 size 1 24 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group liny fill group inr
- gen zone reflect dip 90 dd 90 orig 9 0 0
- gen zone reflect dip 0 dd 0 ori 9 0 0
- ;1/2地层
- gen zone brick p0 18 0 -9 p1 36 0 -9 p2 18 36 -9 p3 18 0 9 size 10 12 8 group outr1
- gen zone brick p0 18 0 9 p1 36 0 9 p2 18 36 9 p3 18 0 13 size 10 12 2 group outr2
- gen zone brick p0 0 0 9 p1 18 0 9 p2 0 36 9 p3 0 0 13 size 8 12 2 group outr2
- gen zone brick p0 18 0 -23 p1 36 0 -23 p2 18 36 -23 p3 18 0 -9 size 10 12 8 group outr3
- gen zone brick p0 0 0 -23 p1 18 0 -23 p2 0 36 -23 p3 0 0 -9 size 8 12 8 group outr3
- ;隧道全形成
- gen zone reflect dip 90 dd 270 orig 0 0 0
- group linz range group liny x -12.1 -5.99
- plot add axes red
- plot set rotation 20 0 45
- plot block group
根据隧道及地下工程FLAC解析方法一书当中的例子，建立上述模型
赋予重力
然后对左边的隧道求解
第一步开挖当中的第5行里面的inrz是哪个组？前面建模的时候没有inrz这个组，在3.0运行的时候也会提示不正确，但5.0能顺利运行，求解
- ;左开挖
- ;step1
- rest Gset.sav
- initial xdisp=0 ydisp=0 zdisp=0
- model null range group inrz y 0 6
- model elas range group linz y 0 6
- ini den 2400 range group linz y 0 6
- prop bulk=16.2e9 shear=11.8e9 range group linz y 0 6
- step 1000
- save l_st1.sav
