FLAC3D隧道开挖模拟程序的锚杆部分

;第一种锚杆（左洞）
def fen_jiao
h1=float(164.5452/14) ;拱部锚杆
end
fen_jiao

def cab1_left ;左洞第一种锚杆
array x1_b(1,20),z1_b(1,20),x1_e(1,20),z1_e(1,20)
loop n(1,8)
x1_b(1,n)=5.85*cos(degrad*(90-(8-n)*h1))-11.85
z1_b(1,n)=5.85*sin(degrad*(90-(8-n)*h1))
x1_e(1,n)=8.85*cos(degrad*(90-(8-n)*h1))-11.85
z1_e(1,n)=8.85*sin(degrad*(90-(8-n)*h1))
end_loop

loop n(9,15)
x1_b(1,n)=-5.85*sin(degrad*(n-8)*h1)-11.85
z1_b(1,n)=5.85*cos(degrad*(n-8)*h1)
x1_e(1,n)=-8.85*sin(degrad*(n-8)*h1)-11.85
z1_e(1,n)=8.85*cos(degrad*(n-8)*h1)
end_loop

loop n(16,16)
x1_b(1,n)=-8.35*cos(degrad*2.9574)-9.35
z1_b(1,n)=-8.35*sin(degrad*2.9574)
x1_e(1,n)=-11.35*cos(degrad*2.9574)-9.35
z1_e(1,n)=-11.35*sin(degrad*2.9574)
end_loop

loop n(17,17)
x1_b(1,n)=-8.35*cos(degrad*11.1915)-9.35
z1_b(1,n)=-8.35*sin(degrad*11.1915)
x1_e(1,n)=-11.35*cos(degrad*11.1915)-9.35
z1_e(1,n)=-11.35*sin(degrad*11.1915)
end_loop

loop n(18,18)
x1_b(1,n)=-14.35+8.35*cos(degrad*11.1915)
z1_b(1,n)=-8.35*sin(degrad*11.1915)
x1_e(1,n)=-14.35+11.35*cos(degrad*11.1915)
z1_e(1,n)=-11.35*sin(degrad*11.1915)
end_loop

loop n(19,19)
x1_b(1,n)=-14.35+8.35*cos(degrad*2.9574)
z1_b(1,n)=-8.35*sin(degrad*2.9574)
x1_e(1,n)=-14.35+11.35*cos(degrad*2.9574)
z1_e(1,n)=-11.35*sin(degrad*2.9574)
end_loop

end

;左洞第二种锚杆

def fen_jiao2
h2=float(176.2934/15) ;左洞第二种锚杆拱部
end

fen_jiao2

def cab2_left
;array x2_b(1,20),z2_b(1,20),x2_e(1,20),z2_e(1,20)
loop n(1,6)
x1_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h2))
x1_e(1,n)=9.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_e(1,n)=9.85*sin(degrad*(1.8533+(n-1)*h2))
end_loop

loop n(7,10)
x1_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h2))
x1_e(1,n)=8.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_e(1,n)=8.85*sin(degrad*(1.8533+(n-1)*h2))

loop n(11,16)
x1_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h2))
x1_e(1,n)=9.85*cos(degrad*(1.8533+(n-1)*h2))-11.85
z1_e(1,n)=9.85*sin(degrad*(1.8533+(n-1)*h2))
end_loop

loop n(17,17)
x1_b(1,n)=-8.35*cos(degrad*6.9372)-9.35
z1_b(1,n)=-8.35*sin(degrad*6.9372)
x1_e(1,n)=-12.35*cos(degrad*6.9372)-9.35
z1_e(1,n)=-12.35*sin(degrad*6.9372)
end_loop

loop n(18,18)
x1_b(1,n)=-1.40*cos(degrad*22.7303)-16.1037
z1_b(1,n)=-1.40*sin(degrad*22.7303)-1.64
x1_e(1,n)=-5.4*cos(degrad*22.7303)-16.1037
z1_e(1,n)=-5.4*sin(degrad*22.7303)-1.64
end_loop

loop n(19,19)
x1_b(1,n)=1.4*cos(degrad*22.7303)-7.9563
z1_b(1,n)=-1.4*sin(degrad*22.7303)-1.64
x1_e(1,n)=5.4*cos(degrad*22.7303)-7.5963
z1_e(1,n)=-5.4*sin(degrad*22.7303)-1.64
end_loop

loop n(20,20)
x1_b(1,n)=8.35*cos

(degrad*6.9372)-11.85
z1_b(1,n)=-8.35*sin(degrad*6.9372)
x1_e(1,n)=12.35*cos(degrad*6.9372)-11.85
z1_e(1,n)=-12.35*sin(degrad*6.9372)
end_loop

end

;右洞第一种锚杆

def fen_jiao3
h3=float(164.5452/14) ;拱部锚杆段数
end
fen_jiao3

def cab1_right
array x2_b(1,20),z2_b(1,20),x2_e(1,20),z2_e(1,20)

loop n(1,8)
x2_b(1,n)=5.85*cos(degrad*(90-(8-n)*h3))+11.85
z2_b(1,n)=5.85*sin(degrad*(90-(8-n)*h3))
x2_e(1,n)=8.85*cos(degrad*(90-(8-n)*h3))+11.85
z2_e(1,n)=8.85*sin(degrad*(90-(8-n)*h3))
end_loop

loop n(9,15)
x2_b(1,n)=-5.85*sin(degrad*(n-8)*h3)+11.85
z2_b(1,n)=5.85*cos(degrad*(n-8)*h3)
x2_e(1,n)=-8.85*sin(degrad*(n-8)*h3)+11.85
z2_e(1,n)=-8.85*cos(degrad*(n-8)*h3)
end_loop

loop n(16,16)
x2_b(1,n)=14.35-8.35*cos(2.9574*degrad)
z2_b(1,n)=-8.35*sin(2.9574*degrad)
x2_e(1,n)=14.35-11.35*cos(2.9574*degrad)
z2_e(1,n)=-11.35*sin(2.9574*degrad)
end_loop

loop n(17,17)
x2_b(1,n)=14.35-8.35*cos(11.1915*degrad)
z2_b(1,n)=-8.35*sin(11.1915*degrad)
x2_e(1,n)=14.35-11.35*cos(11.1915*degrad)
z2_e(1,n)=-11.35*sin(11.1915*degrad)
end_loop

loop n(18,18)
x2_b(1,n)=9.35+8.35*cos(11.1915*degrad)
z2_b(1,n)=-8.35*sin(11.1915*degrad)
x2_e(1,n)=9.35+11.35*cos(11.1915*degrad)
z2_e(1,n)=-8.35*sin(11.1915*degrad)
end_loop

loop n(19,19)
x2_b(1,n)=9.35+8.35*cos(2.9574*degrad)
z2_b(1,n)=-8.35*sin(2.9574*degrad)
x2_e(1,n)=9.35+11.35*cos(2.9574*degrad)
z2_e(1,n)=-8.35*sin(2.9574*degrad)
end_loop

end

;右洞第二种锚杆

def fen_jiao4
h4=float(176.2934/15)
end
fen_jiao4

def cab2_right
;array x2_b(1,20),z2_b(1,20),x2_e(1,20),z2_e(1,20)
loop n(1,6)
x2_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h4))
x2_e(1,n)=9.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_e(1,n)=9.85*sin(degrad*(1.8533+(n-1)*h4))
end_loop

loop n(7,10)
x2_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h4))
x2_e(1,n)=8.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_e(1,n)=8.85*sin(degrad*(1.8533+(n-1)*h4))

loop n(11,16)
x2_b(1,n)=5.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_b(1,n)=5.85*sin(degrad*(1.8533+(n-1)*h4))
x2_e(1,n)=9.85*cos(degrad*(1.8533+(n-1)*h4))+11.85
z2_e(1,n)=9.85*sin(degrad*(1.8533+(n-1)*h4))
end_loop

loop n(17,17)
x2_b(1,n)=-8.35*cos(degrad*6.9372)+11.85
z2_b(1,n)=-8.35*sin(degrad*6.9372)
x2_e(1,n)=-12.35*cos(degrad*6.9372)+11.85
z2_e(1,n)=-12.35*sin(degrad*6.9372)
end_loop

loop n(18,18)
x2_b(1,n)=-1.40*cos(degrad*22.7303)+7.5963
z2_b(1,n)=-1.40*sin(degrad*22.7303)-1.64
x2_e(1,n)=-5.4*cos(degrad*22.7303)+7.5963
z2_e(1,n)=-5.4*sin(degrad*22.7303)-1.64
end_loop

loop n(19,19)
x2_b(1,n)=1.4*cos(degrad*22.7303)+16.1037
z2_b(1,n)=-1.4*sin(degrad*22.7303)-1.64
x2_e(1,n)=5.4*cos(degrad

*22.7303)+16.1037
z2_e(1,n)=-5.4*sin(degrad*22.7303)-1.64
end_loop

loop n(20,20)
x2_b(1,n)=8.35*cos(degrad*6.9372)+9.35
z2_b(1,n)=-8.35*sin(degrad*6.9372)
x2_e(1,n)=12.35*cos(degrad*6.9372)+9.35
z2_e(1,n)=-12.35*sin(degrad*6.9372)
end_loop

end

def chumao_left
x1_left=x1_b(1,n)
z1_left=z1_b(1,n)
x2_left=x1_e(1,n)
z2_left=z1_e(1,n)
end

def chumao_right
x1_right=x2_b(1,n)
z1_right=z2_b(1,n)
x2_right=x2_e(1,n)
z2_right=z2_e(1,n)
end

def number_1
a=u/2.0
b=int(u/2)
if a = b then
cab1_left
else
cab2_left
endif
end

def number_2
a=j/2.0
b=int(j/2)
if a = b then
cab1_left
else
cab2_left
endif
end

def number_3
a=k/2.0
b=int(k/2)
if a = b then
cab1_left
else
cab2_left
endif
end

def number_4
a=x/2.0
b=int(x/2)
if a = b then
cab1_right
else
cab2_right
endif
end


number_1
number_2
number_3
number_4


;左洞上台阶
def left_1
if u=2.0*int(u/2) then
loop n(4,12)
chumao_left
command
sel cable id n begin x1_left,chumaod1_y,z1_left end x2_left,chumaod1_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop
else
loop n(4,13)
chumao_left
command
sel cable id n begin x1_left,chumaod1_y,z1_left end x2_left,chumaod1_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop
end_if
end


;左洞两侧
def left_2
if j=2.0*int(j/2) then
loop n(1,3)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(13,17)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(18,19)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

else

loop n(1,3)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(13,18)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e

9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(19,20)
chumao_left
command
sel cable id n begin x1_left,chumaod2_y,z1_left end x2_left,chumaod2_y,z2_left nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop
end_if
end

;左洞下台阶(无锚杆)


;右洞上台阶

def right_1
if k=2.0*int(k/2) then
loop n(4,12)
chumao_right
command
sel cable id n begin x1_right,chumaod3_y,z1_right end x2_right,chumaod3_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop
else
loop n(4,13)
chumao_right
command
sel cable id n begin x1_right,chumaod3_y,z1_right end x2_right,chumaod3_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop
end_if
end

;右洞两侧

def right_2
if x=2.0*int(x/2) then
loop n(1,3)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(13,17)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(18,19)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

else

loop n(1,3)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(13,18)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_loop

loop n(19,20)
chumao_right
command
sel cable id n begin x1_right,chumaod4_y,z1_right end x2_right,chumaod4_y,z2_right nseg 10
sel cable pro emod 45e9 xcarea 0.3801e-3 gr_per 3.3 yten 25e5 gr_k 17.5e7 gr_c 20e5 range id n
step 1
end_command
end_

loop
end_if
end