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 1 gen 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 2 gen 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_5p5=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 4gen 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_tunnel\1.savplot 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_tunnel\1.savdef 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 2500step 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 2500step 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 2500step 300end_commandcut1_1=cut1_1+s1cut2_1=cut2_1+s1chu1_1=cut1_1chu2_1=cut2_1+0.4end;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;def main_function ;;;开挖主函数s=1.0 ;;;;;开挖步距s1=1.0 ;;;;大跨度开挖步距;xx=s*tan(beta*degrad)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_1-0.5chu2_1=cut2_1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;左边与右边同时开挖10米，然后开挖大跨度10米。