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