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

new

def assi_parameter

rad=3.0 ;隧道洞径

c1=10.0 ;埋深

s=1.0 ;最小净距

beta=15

kd=4.0

b=8.0*kd ;宽度

long=10

long1=long

long2=long1+10.0

x1=kd

y2=long

x2=-1.0*long*tan(beta*degrad)

z3=8.0

c=c1+z3

int_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.0

int_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.1

xz2=xz+0.1

xy1=b-0.1

xy2=b+0.1

zx1=z1_5-0.1

zx2=z1_5+0.1

;;;;;;;;;;;;

n1=b/2.0

n2=long

n3=rad

n3_1=c1/2

n4=3

n5=3.0*kd

end

assi_parameter ;;;运行函数assi_parameter

gen 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,-10

gen 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 2

gen 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 1

gen 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 2

group 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 1

gen 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 1

gen 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