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

newdef assi_parameterrad=3.0 ;隧道洞径c1=10.0 ;埋深s=1.0 ;最小净距beta=15kd=4.0b=8.0*kd ;宽度long=10long1=longlong2=long1+10.0x1=kdy2=longx2=-1.0*long*tan(beta*degrad)z3=8.0c=c1+z3int_x1=-1.0*(2.0*(rad+long*tan(beta*degrad)+s/2.0)+kd)int_x2=-1.0*(s+3.0*rad+long*tan(beta*degrad)int_x3=-1.0*(long*tan(beta*degrad)-rad);x1_1=-1.0*(rad+long*tan(beta*degrad)+s/2.0)z3_1=-1.0*z3;z3_2=-1.0*kdx10=-1.0*long*tan(beta*degrad)-rad;x1_2=-1.0*(rad+long*tan(beta*degrad);x2_3=rad-long*tan(beta*degrad);x0_4=-1.0*radz0_4=2.0*z3_2/3x2_4=-1.0*(rad+long*tan(beta*degrad);z1_5=-6.0*kd ;决定模型底部的高度 rad1=2.0*rad+s/2.0int_rad1=rad1+0.1int_rad2=rad1-0.1;x3_6=-1.0*rad;xz=-1.0*(2.0*(rad+long*tan(beta*degrad)+s/2.0)-bxz1=xz-0.1xz2=xz+0.1xy1=b-0.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=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 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 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 ratio=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 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 group 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 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 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=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 ratio=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 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 zone 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 long 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_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=x2_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 range 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 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_x1 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 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.0*p_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 ten=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 szz 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=r-hdrr=2.0*r+dd/2.0 ;大跨度隧道的洞径rr1=rr+0.1rr2=rr-2*hdc1=10.0 ;埋深b=4.0*kd ;宽度length=10long1=lengthlong2=length+10.0o2=-1.0*length*tan(beta*degrad)o1_z=-1.0*(2.0*(r+length*tan(beta*degrad)+dd)o2_z=-1.0*(2.0*r+dd+length*tan(beta*degrad);x0_1=r-hdz0_1=-2.0*kd/3x2_1=r-long*tan(beta*degrad)-hd;x0_2=-1.0*rx1_2=x0_2+hdx2_2=-1.0*(r+length*tan(beta*degrad);x1_3=-1.0*(rad+length*tan(beta*degrad)+dd/2.0);x0_4=r-long*tan(beta*degrad)-hd2x1_4=r-long*tan(beta*degrad)x1_4_1=x1_4-0.1x1_4_2=x1_4+0.1;x0_5=-1.0*(3.0*r+dd+length*tan(beta*degrad)x1_5=-1.0*(3.0*r+dd+length*tan(beta*degrad)+hd2x0_5_1=x0_5-0.1x0_5_2=x0_5+0.1;z3_6=z0_1+0.1x4_6=r-long*tan(beta*degrad);x0_7=x0_5-0.1x1_7=x1_4+0.1;x0_8=o1_z-rx1_8=x0_8+hdx2_8=o2_z-r;x0_9=o1_z+r-hdx1_9=o1_z+rx2_9=o2_z+r-hd;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 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 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=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=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_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 5gen 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 density 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 chu1_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 zone 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 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+s1ch