FLAC3D学习（SIMWE论坛君之枫整理） zz=gp_zpos(p_gp)-z0 if xx>zz then
ang2=angmax*zz/xx ; remap direction from uniform lengths to uniform angles rr=xx*sqrt2
gp_xpos(p_gp)=x0+rr*cos(ang2) gp_zpos(p_gp)=z0+rr*sin(ang2) else
if zz>xx then
ang2=angmax*xx/zz ; remap direction from uniform lengths to uniform angles rr=zz*sqrt2
gp_xpos(p_gp)=x0+rr*sin(ang2) gp_zpos(p_gp)=z0+rr*cos(ang2) end_if end_if end_if
p_gp=gp_next(p_gp) end_loop end def makeall tunnel_data make1 ; outside top make2 ; outside lower-right make3 ; outside bottom make4 ; outside bottom-right
make5 ; inside bottom (call make5 before make6 so common points mapped with map2) make6 ; inside top make7 ; far top and right
; remap inner tunnel walls into outer arcs
map1a ; inside bottom (adjust right side to vertical) map1 ; inside bottom (adjust right side to sector) map2 ; inside top end makeall gen merge 1e-4
; gen zone reflect orig x0 y1 z0 norm 1 0 0 plo crea qqq plo add block group plo sho
26
FLAC3D学习（SIMWE论坛君之枫整理） pause
group 'top heading' range group g5 group invert range group g6
group rock range group 'top heading' not group invert not pause
14.我的例子，出了问题，但不知道原因 ;梯形溶洞 y 方向 范围 2 10 ;create gong ding
gen zone radcyl size 1 12 10 10 p0 0 0 0 p1 7 0 0 p2 0 12 0 p3 0 0 7＆ p4 7 12 0 p5 0 12 7 p6 7 0 7 p7 7 12 7 p8 5 0 0 p9 0 0 5 ＆ p10 5 12 0 p11 0 12 5 fill
gen zone reflect norm 1 0 0 ori 0 0 0 gen zone reflect norm 0 0 1 ori 0 0 0 group 2 range cyl end1 0 0 0 end2 0 12 0 ra 5 group 1 range x -7 7 y 0 12 z -7 7 group 2 not ;create rongdong
gen zone brick size 13 12 10 p0 7 0 -7 p1 20 0 -7 p2 7 12 -7 p3 7 0 7 group 3 gen zone brick size 16 12 3 p0 -7 0 7 p1 9 0 7 p2 -7 12 7 p3 -7 0 10 ＆ p4 9 12 7 p5 -7 12 10 p6 12 0 10 p7 12 12 10 group 4 ;rongdong position
gen zone brick size 10 2 3 p0 9 0 7 p1 19 0 7 p2 9 2 7 p3 12 0 10 ＆ p4 19 2 7 p5 12 2 10 p6 16 0 10 p7 16 2 10 group 7
gen zone brick size 10 10 3 p0 9 2 7 p1 19 2 7 p2 9 10 7 p3 12 2 10 ＆ p4 19 10 7 p5 12 10 10 p6 16 2 10 p7 16 10 10 group 8
gen zone brick size 10 2 3 p0 9 10 7 p1 19 10 7 p2 9 12 7 p3 12 10 10 ＆ p4 19 12 7 p5 12 12 10 p6 16 10 10 p7 16 12 10 group 9 ;rongdong dingban
gen zone brick size 20 12 2 p0 -7 0 10 p1 20 0 10 p2 -7 12 10 p3 -7 0 12 ＆group 5
;rongdong youbianqiang
gen zone brick size 1 12 3 p0 19 0 7 p1 20 0 7 p2 19 12 7 p3 16 0 10 ＆ p4 20 12 7 p5 16 12 10 p6 20 0 10 p7 20 12 10 group 6 plot add surface yellow plot show ;create macro
macro rock 'prop bulk 2.5e8 shear 3.3e8'
mo ela range group 1 group 2 group 3 group 4 group 5 group 6 group 7 group 9 mo nu range group 8
27
FLAC3D学习（SIMWE论坛君之枫整理）
macro rock range group 1 group 2 group 3 group 4 group 5 group 6 group 7 group 9 plot add surface yellow plot show set grav 0 0 -10 fix x range x -7.1 -6.9 fix x range x 20.1 19.9 fix y range y -0.1 0.1 fix y range y 12.1 11.9 fix z range z -6.9 -7.1 fix z range z 12.1 11.9
ini szz -5e4 grad 0 0 -1e4 range z 12 -7 ini sxx -5e4 grad 0 0 -1e2 range z 12 -7 ini sxx -5e4 grad 0 0 -1e2 range z 12 -7 app szz -5e4 range z 12.1 11.9 app syy -5e4 range z 12.1 11.9 app sxx -5e4 range z 12.1 11.9 set large 我的诊断：
（1）没有定义强度参数
（2）几个group取并集的时候语法错误，具体参考command reference的range这条命令。 没有定义材料参数 修改如下：
gen zone radcyl size 1 12 10 10 p0 0 0 0 p1 7 0 0 p2 0 12 0 p3 0 0 7 ＆
p4 7 12 0 p5 0 12 7 &p6 7 0 7 p7 7 12 7 p8 5 0 0 p9 0 0 5 p10 5 12 0 p11 0 12 5 fill gen zone reflect norm 1 0 0 ori 0 0 0 gen zone reflect norm 0 0 1 ori 0 0 0 group 2 range cyl end1 0 0 0 end2 0 12 0 ra 5 group 1 range x -7 7 y 0 12 z -7 7 group 2 not ;create rongdong
gen zone brick size 13 12 10 p0 7 0 -7 p1 20 0 -7 p2 7 12 -7 p3 7 0 7 group 3 gen zone brick size 16 12 3 p0 -7 0 7 p1 9 0 7 p2 -7 12 7 p3 -7 0 10 p4 9 12 7 p5 -7 12 10 p6 12 0 10 p7 12 12 10 group 4 ;rongdong position
gen zone brick size 10 2 3 p0 9 0 7 p1 19 0 7 p2 9 2 7 p3 12 0 10 p4 19 2 7 p5 12 2 10 p6 16 0 10 p7 16 2 10 group 7
gen zone brick size 10 10 3 p0 9 2 7 p1 19 2 7 p2 9 10 7 p3 12 2 10 p4 19 10 7 p5 12 10 10 p6 16 2 10 p7 16 10 10 group 8
gen zone brick size 10 2 3 p0 9 10 7 p1 19 10 7 p2 9 12 7 p3 12 10 10 p4 19 12 7
28
FLAC3D学习（SIMWE论坛君之枫整理） p5 12 12 10 p6 16 10 10 p7 16 12 10 group 9 ;rongdong dingban
gen zone brick size 20 12 2 p0 -7 0 10 p1 20 0 10 p2 -7 12 10 p3 -7 0 12 group 5 ;rongdong youbianqiang
gen zone brick size 1 12 3 p0 19 0 7 p1 20 0 7 p2 19 12 7 p3 16 0 10 p4 20 12 7 p5 16 12 10 p6 20 0 10 p7 20 12 10 group 6 plot add surface yellow plot show ;create macro
;prop bulk 2.5e8 shear 3.3e8 'prop bulk 2.5e8 shear 3.3e8' mo ela
prop bulk 2.5e8 shear 3.3e8 range group 1 prop bulk 2.5e8 shear 3.3e8 range group 2 prop bulk 2.5e8 shear 3.3e8 range group 3 prop bulk 2.5e8 shear 3.3e8 range group 4 prop bulk 2.5e8 shear 3.3e8 range group 5 prop bulk 2.5e8 shear 3.3e8 range group 6 prop bulk 2.5e8 shear 3.3e8 range group 7 prop bulk 2.5e8 shear 3.3e8 range group 9 mo nu range group 8 plot show set grav 0 0 -10 fix x range x -7.1 -6.9 fix x range x 20.1 19.9 fix y range y -0.1 0.1 fix y range y 12.1 11.9 fix z range z -6.9 -7.1 fix z range z 12.1 11.9
ini szz -5e4 grad 0 0 -1e4 range z 12 -7 ini sxx -5e4 grad 0 0 -1e2 range z 12 -7 ini sxx -5e4 grad 0 0 -1e2 range z 12 -7 app szz -5e4 range z 12.1 11.9 app syy -5e4 range z 12.1 11.9 app sxx -5e4 range z 12.1 11.9 set large
15.建立interface的几种方法 首先给出一个简单的例子 new
29
FLAC3D学习（SIMWE论坛君之枫整理）
gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1 gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &p5 0 10 90 p6 47 0 90 p7 47 10 90 size 11 1 9 group 3
gen zone wedge p0 66.4 0 90.1 p1 47 0 90.1 p2 66.4 10 90.1 p3 111.4 0 45.1 &p4 47 10 90.1 p5 111.4 10 45.1 size 5 1 9 group 4
;interface 1 face rang plane norm 45,0,64.4 origin 47 5 90
interface 1 face range plane norm 45,0,64.4 origin 47 5 90 dist 0.1 ini z add -0.1 range group 4
gen zone brick p0 111.4 0 0 p1 128 0 0 p2 111.4 10 0 p3 111.4 0 45 & size 2 1 5 group 2 model mohr
pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10 range group 1 pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10 range group 2 int 1 prop ks 2e12 kn 2e12 fric 36 pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10 range group 3 pro bulk=1e8 she=3e7 dens 2e3 coh 12e3 &
friction 20 dilation 0. tension 1.0e10 range group 4 fix z range z -0.1 0.1 fix x range x -0.1 0.1 fix x range x 127.9 128.1 fix y
set gravity 0 0 -10 step 1000 ini state 0
ini xdis 0 ydis 0 zdis 0 ;ini xvel 0 yvel 0 zvel 0 ;set large step 1000
;solve fos file cable.sav associated 修改上面蓝色字体部分
第一种建立interface的方法（书上的移来移去式）：
gen zone brick p0 0 0 0 p1 111.4 0 0 p2 0 10 0 p3 0 0 45 size 11 1 5 group 1 gen zone brick p0 0 0 45 p1 111.4 0 45 p2 0 10 45 p3 0 0 90 p4 111.4 10 45 &p5 0 10 90 p6 47 0 90 p7 47 10 90 size 11 1 9 group 3
gen zone wedge p0 66.4 0 90.1 p1 47 0 90.1 p2 66.4 10 90.1 p3 111.4 0 45.1 &p4 47
30