分岔隧道FLAC3D数值模拟命令流.doc

《分岔隧道FLAC3D数值模拟命令流.doc》由会员分享，可在线阅读，更多相关《分岔隧道FLAC3D数值模拟命令流.doc（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、分岔隧道FLAC3D数值模拟命令流newdef assi_parameterrad=3。0 ；隧道洞径c1=10。0 ；埋深s=1.0 ；最小净距beta=15kd=4。0b=8.0kd ；宽度long=10long1=longlong2=long1+10.0x1=kdy2=longx2=1。0long*tan（beta*degrad）z3=8.0c=c1+z3int_x1=1。0（2。0*(rad+longtan(betadegrad）+s/2。0）+kd)int_x2=-1.0（s+3。0rad+longtan(beta*degrad)）int_x3=-1.0(long*tan（betad

2、egrad)-rad);；;；;；;；;x1_1=1。0（rad+longtan(beta*degrad)+s/2。0）z3_1=-1.0*z3;；;；;；;z3_2=-1。0*kdx10=-1.0*long*tan(betadegrad）-rad；;；;；;；x1_2=-1.0(rad+longtan(beta*degrad）;；;；x2_3=rad-longtan(beta*degrad)；;；;；x0_4=1。0radz0_4=2.0z3_2/3x2_4=1。0（rad+longtan(betadegrad)）；;；;z1_5=6.0kd ；决定模型底部的高度 rad1=2.0*rad+

3、s/2。0int_rad1=rad1+0。1int_rad2=rad10。1；;；;；;；;；x3_6=1。0rad;；;；;；xz=-1.0*（2.0(rad+long*tan(betadegrad）+s/2.0）-bxz1=xz0。1xz2=xz+0.1xy1=b0。1xy2=b+0。1zx1=z1_5-0.1zx2=z1_5+0。1；;；;；;n1=b/2.0n2=longn3=radn3_1=c1/2n4=3n5=3.0*kdendassi_parameter ;；;运行函数assi_parametergen zone radcylinder p0=0 0 0 p1=x1 0 0 p2

4、=x2 y2 0 p3=0 0 z3 p4=x1 long 0 p5=x2 y2 z3 p6=x1 0 z3 p7=x1 long z3 dimension=rad rad rad rad size=rad long 6 n4 ratio=1 1 1 1 group 1gen zone radcylinder p0=0 0 0 p1=x1_1 0 0 p2=x2 y2 0 p3=0 0 z3_1 p4= x1_1 long 0 p5=x2 long z3_1 p6=x1_1 0 z3_1 p7=x1_1 long z3_1 dimension=rad rad rad rad size=rad

5、long 6 n4 ratio=1 1 1 1 group 2gen zone reflect dip=0 dd=0 origin=0 0 0 range group 2 ;以水平面为对称面镜像del zone range z=0,10gen zone brick p0=rad 0 0 p1=rad 0 z3_2 p2=x2_3 y2 0 p3=x1 0 0 p4=x2_3 y2 z3_2 p5=x1 long 0 p6=x1 0 z3_2 p7=x1 long z3_2 size=rad long n4 ratio=1 1 1 group 2gen zone brick p0=x0_4 0

6、z0_4 p1=x0_4 0 z3_2 p2=x2_4 y2 z0_4 p3=rad 0 z0_4 p4=x2_4 y2 z3_2 p5=x2_3 long z0_4 p6=rad 0 z3_2 p7=x2_3 long z3_2 size=1 long 6 ratio=1 1 1 group 1gen zone brick p0=x1_1 0 0 p1=x1_1 0 z3_2 p2=x1_1 y2 0 p3=x0_4 0 0 p4=x1_1 y2 z3_2 p5=x2_4 long 0 p6=x0_4 0 z3_2 p7=x2_4 long z3_2 size=rad long n4 rat

7、io=1 1 1 group 2group 1 range group 2；;；;；生成洞身右侧围岩gen zone brick p0=x1 0 0 p1=b 0 0 p2=x1 long2 0 p3=x1 0 z3 size=n1 long2 rad ratio=1 1 1 group 1gen zone brick p0=x1 0 z3_2 p1=b 0 z3_2 p2=x1 long2 z3_2 p3=x1 0 0 size=n1 long2 rad ratio=1 1 1 group 1;；;；;生成洞底下方围岩gen zone brick p0=rad 0 z3_2 p1=rad 0

8、 z1_5 p2=x2_3 y2 z3_2 p3=x1 0 z3_2 p4=x2_3 y2 z1_5 p5=x1 long z3_2 p6=x1 0 z1_5 p7=x1 long z1_5 size=n5 long n4 ratio=1 1 1 group 1gen zone brick p0=x0_4 0 z3_2 p1=x0_4 0 z1_5 p2=x2_4 y2 z3_2 p3=rad 0 z3_2 p4=x2_4 y2 z1_5 p5=x2_3 long z3_2 p6=rad 0 z1_5 p7=x2_3 long z1_5 size=n5 long 6 ratio=1 1 1 g

9、roup 1gen zone brick p0=x1_1 0 z3_2 p1=x1_1 0 z1_5 p2=x1_1 y2 z3_2 p3=x0_4 0 z3_2 p4=x1_1 y2 z1_5 p5=x2_4 long z3_2 p6=x0_4 0 z1_5 p7=x2_4 long z1_5 size=n5 long n4 ratio=1 1 1 group 1；;；;生成洞顶上部围岩gen zone brick p0=0 0 z3 p1=x1 0 z3 p2=x2 y2 z3 p3=0 0 c p4=x1 long z3 p5=x2 y2 c p6=x1 0 c p7=x1 long c

10、 size=rad long n3_1 ratio=1 1 1 group 1gen zone brick p0=x1_1 0 z3 p1=0 0 z3 p2=x1_1 long z3 p3=x1_1 0 c p4=x2 long z3 p5=x1_1 long c p6=0 0 c p7=x2 long c size=rad long n3_1 ratio=1 1 1 group 1gen zone brick p0=x1 0 z1_5 p1=b 0 z1_5 p2=x1 long2 z1_5 p3=x1 0 z3_2 size=n1 long2 n5 ratio=1 1 1 1 group

11、 1gen zone brick p0=x1 0 z3 p1=b 0 z3 p2=x1 long2 z3 p3=x1 0 c size=n1 long2 n3_1 ratio=1 1 1 1 group 1;；;；;；生成大跨度段洞顶上方围岩gen zone brick p0=x2 long1 z3 p1=x1 long1 z3 p2=x2 long2 z3 p3=x2 long1 c size=rad long n3_1 ratio=1 1 1 group 1gen zone brick p0=x1_1 long1 z3 p1=x2 long1 z3 p2=x1_1 long2 z3 p3=

12、x1_1 long1 c size=rad long n3_1 ratio=1 1 1 group 1；;；;生成大跨度段洞底下方围岩gen zone brick p0=x2_3 long1 z3_2 p1=x2_3 long1 z1_5 p2=x2_3 long2 z3_2 p3=x1 long1 z3_2 size=n5 long n4 ratio=1 1 1 group 1gen zone brick p0=x2_4 long1 z3_2 p1=x2_4 long1 z1_5 p2=x2_4 long2 z3_2 p3=x2_3 long1 z3_2 size=n5 long 6 rat

13、io=1 1 1 group 1gen zone brick p0=x1_1 long1 z3_2 p1=x1_1 long1 z1_5 p2=x1_1 long2 z3_2 p3=x2_4 long1 z3_2 size=n5 long n4 ratio=1 1 1 group 1;；;；生成大跨度段洞身直边墙gen zone brick p0=x2_3 long1 0 p1=x2_3 long1 z3_2 p2=x2_3 long2 0 p3=x1 long1 0 size=n4 long n4 ratio=1 1 1 group 1gen zone brick p0=x2_4 long1

14、 z0_4 p1=x2_4 long1 z3_2 p2=x2_4 long2 z0_4 p3=x2_3 long1 z0_4 size=1 long 6 ratio=1 1 1 group 1gen zone brick p0=x1_1 long1 z0_4 p1=x1_1 long1 z3_2 p2=x1_1 long2 z0_4 p3=x2_4 long1 z0_4 size=1 long n4 ratio=1 1 1 group 1gen zone reflect dip=90 dd=90 origin=x1_1 0 0 ；以竖直面为对称面镜像；;；;；;生成大跨度段洞身拱顶gen zo

15、ne radcylinder p0=x1_1 long1 0 p1=x1 long1 0 p2=x1_1 long2 0 p3=x1_1 long1 z3 dimension=rad1 rad1 rad1 rad1 size=1 10 9 n4 ratio=1 1 0.87 1 group 2gen zone reflect dip=90 dd=90 origin=x1_1 0 0 range group 2 ；以竖直面为对称面镜像;；; 生成小净距隧道gen zone cylinder p0=0 0 0 p1= rad 0 0 p2=x2 y2 0 p3=0 0 rad size rad l

16、ong 6 group 3gen zone cylinder p0=0 0 0 p1= 0 0 rad p2=x2 y2 0 p3=x3_6 0 0 size rad long 6 group 3gen zone brick p0=x0_4 0 z0_4 p1=rad 0 z0_4 p2=x2_4 y2 z0_4 p3=x0_4 0 0 p4=x2_3 long z0_4 p5=x2_4 y2 0 p6=rad 0 0 p7=x2_3 long 0 size=6 long 2 ratio=1 1 1 group 3gen zone reflect dip=90 dd=90 origin=x1_

17、1 0 0 range group 3；以竖直面为对称面镜像;；;；; 生成大跨度隧道gen zone cylinder p0=x1_1 long1 0 p1=add(rad1 0 0） p2=x1_1 long2 0 p3=x1_1 long1 rad1 size rad1 long 9 ratio=1 1 0。87 1 group 4 gen zone brick p0=x2_4 long1 z0_4 p1=x2_3 long1 z0_4 p2=x2_4 long2 z0_4 p3=x1_1 long1 0 p4=x2_3 long2 z0_4 & p5=x1_1 long2 0 p6=x

18、2_3 long1 0 p7=x2_3 long2 0 size=6 long 2 ratio=1 1 1 group 5gen zone wedge p0=x1_1 long1 z0_4 p1=x2_4 long1 z0_4 p2=x1_1 long2 z0_4 p3=x1_1 long1 0 size=3 long 2 ratio=1 1 1 group 5gen zone reflect dip=90 dd=90 origin=x1_1 0 0 range group 4；以竖直面为对称面镜像gen zone reflect dip=90 dd=90 origin=x1_1 0 0 ra

19、nge group 5;以竖直面为对称面镜像generate merge 0.3 range x=x1 y=long1 long2 z=0 z3generate merge 0.3 range x=int_x1 y=long1 long2 z=0 z3generate merge 0。83 range x=int_x1 x1 y=long1 long2 z=z3generate merge 0.202 range x=int_x1 x1 y=long1 z=0 z3attach face range x=int_x1 x1 y=long1 z=0 z3;gen separate 4;int 1

20、 face range x=int_x2 int_x3 y=long1 long2 z=0 ；int 1 face range x=int_x2 int_x3 z=0 rad1 ；cylinder end1=x1_1 long1 0 end2=x1_1 long2 0 radius=int_rad1 cylinder end1=x1_1 long1 0 end2=x1_1 long2 0 radius=int_rad2 not；gen separate 2；int 1 face range x=int_x1 x1 y=long1 z=0 z3 ；int 1 face range x=int_x

21、1 x1 y=long1 long2 z=z3;int 1 face range x=int_x1 y=long1 long2 z=0 z3;int 1 face range x=x1 y=long1 long2 z=0 z3group 1 range group 2group 4 range group 5group 2 range x=x1_1 x1 group 3set grav 0 0 10fix x range x xz1 xz2fix x range x xy1 xy2fix y range y -0.1 0。1fix y range y 19。9 20.1fix z range

22、z zx1 zx2model mohrdef derives_mod=y_mod/(2.0*（1。0+p_ratio）) ；s_mod为剪切模量,y_mod为杨氏模量，p_ratio为泊松比b_mod=y_mod/(3。0*（1。0-2.0p_ratio)） ;b_mod为体积模量；k=10（b_mod+4.0*s_mod/3.0)/0。2endset y_mod=2e9 p_ratio=0.33derive ;执行函数deriveprop bulk=b_mod shear=s_mod fric=20 coh=2e4;int 1 prop kn=k ks=k coh=2e4 fric=89 t

23、en=10e20 dila=0。0ini dens 1500solve eini xd=0ini yd=0ini zd=0ini xv=0ini yv=0ini zv=0ini state=0save F:bifurcation_tunnel1。sav plot block groupplot add axes red;plot sketch;plot inter red ske;plot set plane ori x1_1 0 0 nrom 0 1 0；plot set plane ori x1_1 0 0 nrom 1 0 0;plot con szz plane;plot con sz

24、z ou on magf 20 ;；围岩竖向应力显示;plot con zdis ou on magf 20 ；围岩竖向位移显示；plot con xdis ou on magf 20 ;；;围岩水平位移显示；generate merge 0。2 range x=x1 y=long1 long2 z=0 z3newrest F：bifurcation_tunnel1。sav def assi_parameterdd=1。0 ;最小净距r=3.0 ；小净距隧道的洞径beta=15hd=0.4hd2=0。6kd=4。0r1=r+0。1r2=rhdrr=2.0*r+dd/2。0 ；大跨度隧道的洞径r

25、r1=rr+0.1rr2=rr2hdc1=10。0 ;埋深b=4。0*kd ;宽度length=10long1=lengthlong2=length+10。0o2=1。0lengthtan(betadegrad）o1_z=-1.0*(2。0(r+lengthtan(beta*degrad)）+dd)o2_z=-1.0*(2。0r+dd+lengthtan(beta*degrad）;；;；;；;；x0_1=rhdz0_1=-2。0kd/3x2_1=r-long*tan（betadegrad）hd;；;；;；;；;x0_2=1.0rx1_2=x0_2+hdx2_2=-1。0(r+lengthtan

26、（betadegrad）);；;；x1_3=1.0*(rad+length*tan（betadegrad)+dd/2。0)；;；;；;；;x0_4=r-longtan(beta*degrad)hd2x1_4=r-long*tan(betadegrad）x1_4_1=x1_40.1x1_4_2=x1_4+0。1;；;；;；x0_5=1。0(3.0*r+dd+length*tan（beta*degrad)x1_5=1。0*（3。0r+dd+lengthtan(betadegrad）+hd2x0_5_1=x0_5-0.1x0_5_2=x0_5+0。1;；;；;；z3_6=z0_1+0.1x4_6=r

27、-longtan（beta*degrad)；;；;；x0_7=x0_50。1x1_7=x1_4+0.1；;；;；;；;；x0_8=o1_z-rx1_8=x0_8+hdx2_8=o2_zr;；;；;；x0_9=o1_z+rhdx1_9=o1_z+rx2_9=o2_z+rhd；;；;；;；x4_10=o2_z+rendassi_parameter ；;运行函数assi_parameterdef exca_zb ；开挖左边洞室commandm n range group 3 y cut1_z cut2_zstep 100end_commandcommand ;；初期支护sel shell id=1

28、range y=chu1_z chu2_z z=0 r cylinder end1 o1_z 0 0 end2 o2_z long1 0 rad r1 cylinder end1 o1_z 0 0 end2 o2_z long1 0 rad r2 notgen zone brick p0=x0_8 0 z0_1 p1=x1_8 0 z0_1 p2=x2_8 length z0_1 p3=x0_8 0 0.1 size 1 10 2 group 6sel shell id=2 range y=chu1_z chu2_z z=z0_1 0.1 group 6del zone range group

29、 6gen zone brick p0=x0_9 0 z0_1 p1=x1_9 0 z0_1 p2=x2_9 length z0_1 p3=x0_9 0 0。1 size 1 10 2 group 6sel shell id=2 range y=chu1_z chu2_z z=z0_1 0.1 group 6del zone range group 6gen zone brick p0=x0_8 0 z0_1 p1=x1_9 0 z0_1 p2=x2_8 length z0_1 p3=x0_8 0 z3_6 p4=x4_10 length z0_1 p5=x2_8 length z3_6 p6

30、=x1_9 0 z3_6 p7=x4_10 length z3_6 & size 6 10 1 group 6sel shell id=2 range y=chu1_z chu2_z z=z0_1 group 6del zone range group 6sel shell id=1 prop isotropic 21e9，0.20 thickness 0.6 density 2500sel shell id=2 prop isotropic 21e9,0.20 thickness hd density 2500 step 300end_commandcut1_z=cut1_z+scut2_z

31、=cut2_z+schu1_z=cut1_zchu2_z=cut2_z+0。1enddef exca_yb ；;；开挖右边洞室commandm n range group 2 y cut1 cut2step 100end_commandcommand ;；;初期支护sel shell id=1 range y=chu1 chu2 z=0 r cylinder end1 0 0 0 end2 o2 long1 0 rad r1 cylinder end1 0 0 0 end2 o2 long1 0 rad r2 notgen zone brick p0=x0_1 0 z0_1 p1=r 0 z0

32、_1 p2=x2_1 length z0_1 p3=x0_1 0 0。1 size 1 10 2 group 5sel shell id=2 range y=chu1 chu2 z=z0_1 0.1 group 5del zone range group 5gen zone brick p0=x0_2 0 z0_1 p1=x1_2 0 z0_1 p2=x2_2 length z0_1 p3=x0_2 0 0。1 size 1 10 2 group 5sel shell id=2 range y=chu1 chu2 z=z0_1 0.1 group 5del zone range group 5

33、gen zone brick p0=x0_2 0 z0_1 p1=rad 0 z0_1 p2=x2_2 length z0_1 p3=x0_2 0 z3_6 p4=x4_6 length z0_1 p5=x2_2 length z3_6 p6=rad 0 z3_6 p7=x4_6 length z3_6 & size 6 10 1 group 5sel shell id=2 range y=chu1 chu2 z=z0_1 group 5del zone range group 5sel shell id=1 prop isotropic 21e9，0。20 thickness 0.6 den

34、sity 2500sel shell id=2 prop isotropic 21e9，0.20 thickness hd density 2500 step 300end_commandcut1=cut1+scut2=cut2+schu1=cut1chu2=cut2+0.1;x1=x2;x2=x2-xxenddef exca_dkd ；;开挖大跨度洞室commandm n range group 4 y cut1_1 cut2_1step 100end_commandcommand ;；;初期支护sel shell id=1 range z=0 rr cylinder end1 x1_3 c

35、hu1_1 0 end2 x1_3 chu2_1 0 rad rr1 cylinder end1 x1_3 chu1_1 0 end2 x1_3 chu2_1 0 rad rr2 notgen zone brick p0=x0_4 long1 z0_1 p1=x1_4 long1 z0_1 p2=x0_4 long2 z0_1 p3=x0_4 long1 0。1 size 1 10 2 group 5sel shell id=2 range x=x1_4_1 x1_4_2 y=chu1_1 chu2_1 z=z0_1 0。1 group 5del zone range group 5gen z

36、one brick p0=x0_5 long1 z0_1 p1=x1_5 long1 z0_1 p2=x0_5 long2 z0_1 p3=x0_5 long1 0.1 size 1 10 2 group 5sel shell id=2 range x=x0_5_1 x0_5_2 y=chu1_1 chu2_1 z=z0_1 0.1 group 5del zone range group 5gen zone brick p0=x0_7 long1 z0_1 p1=x1_7 long1 z0_1 p2=x0_7 long2 z0_1 p3=x0_7 long1 z3_6 size 13 10 1

37、 group 5sel shell id=2 range y=chu1_1 chu2_1 z=z0_1 group 5del zone range group 5sel shell id=1 prop isotropic 21e9，0。20 thickness 0.8 density 2500sel shell id=2 prop isotropic 21e9，0.20 thickness hd2 density 2500 step 300end_commandcut1_1=cut1_1+s1cut2_1=cut2_1+s1chu1_1=cut1_1chu2_1=cut2_1+0.4end；;

38、；;；;；;；;；;；;；;def main_function ；;；开挖主函数s=1.0 ；;；开挖步距s1=1。0 ；;；大跨度开挖步距；xx=s*tan（betadegrad)cut1=0cut2=schu1=cut1-0.1chu2=cut2；x1=0;x2=x1-xxcut1_z=0cut2_z=schu1_z=cut1_z-0.1chu2_z=cut2_zcut1_1=lengthcut2_1=length+s1chu1_1=cut1_10.5chu2_1=cut2_1;；;；;；;；;；;；;；;；;；;；;；;；;左边与右边同时开挖10米,然后开挖大跨度10米。exca_zb

39、；左边开挖1米exca_yb ；;；右边开挖1米exca_zb ；;；左边开挖2米exca_yb ；;；右边开挖2米exca_zb ;；左边开挖3米exca_yb ；;右边开挖3米exca_zb ；;左边开挖4米exca_yb ；;；;右边开挖4米exca_zb ；;左边开挖5米exca_yb ;；;右边开挖5米exca_zb ;；左边开挖6米exca_yb ;；;右边开挖6米exca_zb ;；左边开挖7米exca_yb ;；右边开挖7米exca_zb ;；左边开挖8米exca_yb ；;；右边开挖8米exca_zb ;；;左边开挖9米exca_yb ;；;右边开挖9米exca_zb ；;左

40、边开挖10米exca_yb ；;右边开挖10米exca_dkd ;；;；大跨度开挖1米exca_dkdexca_dkdexca_dkdexca_dkdexca_dkdexca_dkdexca_dkdexca_dkdexca_dkdendset log onmain_functionsave F：bifurcation_tunnelexcavate_1.sav plot block groupplot add axes red;plot sketch；plot inter red ske;plot set plane ori x1_1 0 0 nrom 0 1 0；plot set plane

41、ori x1_1 0 0 nrom 1 0 0；plot con szz plane；plot con szz ou on magf 20 ;；围岩竖向应力显示;plot con zdis ou on magf 20 ；;；围岩竖向位移显示;plot con xdis ou on magf 20 ;；;围岩水平位移显示;镜头转换；% ；plot set rotation 0 0 40 ;plot set center 0 0 0 ;plot set magnification 3 ；plot ske ;plot add sel geom cable ；plot add sel geom shell yellow ；set movie avi step 1 file F:excavate1.avi ;movie start;%