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

new def assi_parameterrad=3.0 ;隧道洞径c1=10.0 ;埋深s=1.0 ;最小净距beta=15 kd=4.0b=8.0*kd ;宽度long=10long1=long long2=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*kd x10=-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*rad z0_4=2.0*z3_2/3 x2_4=-1.0*(rad+long*tan(beta*degrad));;;;;;;;;;;z1_5=-6.0*kd ;决定模型底部的高度rad1=2.0*rad+s/2.0int_rad1=rad1+0.1 int_rad2=rad1-0.1;;;;;;;;;;;;;x3_6=-1.0*rad;;;;;;;;;;;;; xz=-1.0*(2.0*(rad+long*tan(beta*degrad)+s/2.0))-b xz1=xz-0.1xz2=xz+0.1xy1=b-0.1xy2=b+0.1zx1=z1_5-0.1 zx2=z1_5+0.1;;;;;;;;;;;;n1=b/2.0 n2=long n3=rad n3_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 &。