巷道开挖支护.docVIP

巷道开挖支护 new rest 上区段开挖.sav ;巷道开挖支护 ============================================================================================================ plot set rot 0 0 0 plot set mag 10 plot set center 0 0 0 pl bl state plot show set large def tunnel_excavate cut=0 loop while cut=95 cut1=cut+5 i=cut1 j=cut command m n range x -1.75 1.75 y cut cut1 z -1.3 1.3 end_command loop while ji n1=13*j+1 n2=13*j+2 n3=13*j+3 n4=13*j+4 n5=13*j+5 n6=13*j+6 n7=13*j+7 n8=13*j+8 n9=13*j+9 n10=13*j+10 n11=13*j+11 n12=13*j+12 n13=13*j+13 a1=n1*10-9 a2=n1*10-8 a3=n1*10-4 a4=n1*10-3 a5=n1*10 b1=n2*10-9 b2=n2*10-8 b3=n2*10-4 b4=n2*10-3 b5=n2*10 c1=n3*10-9 c2=n3*10-8 c3=n3*10-4 c4=n3*10-3 c5=n3*10 d1=n4*10-9 d2=n4*10-8 d3=n4*10-4 d4=n4*10-3 d5=n4*10 e1=n5*10-9 e2=n5*10-8 e3=n5*10-4 e4=n5*10-3 e5=n5*10 f1=n6*10-9 f2=n6*10-8 f3=n6*10-4 f4=n6*10-3 f5=n6*10 g1=n7*10-9 g2=n7*10-8 g3=n7*10-4 g4=n7*10-3 g5=n7*10 h1=n8*10-9 h2=n8*10-8 h3=n8*10-4 h4=n8*10-3 h5=n8*10 i1=n9*10-9 i2=n9*10-8 i3=n9*10-4 i4=n9*10-3 i5=n9*10 j1=n10*10-9 j2=n10*10-8 j3=n10*10-4 j4=n10*10-3 j5=n10*10 k1=n11*10-9 k2=n11*10-8 k3=n11*10-4 k4=n11*10-3 k5=n11*10 l1=n12*10-9 l2=n12*10-8 l3=n12*10-4 l4=n12*10-3 l5=n12*10 m1=n13*10-9 m2=n13*10-8 m3=n13*10-4 m4=n13*10-3 m5=n13*10 command sel cable id=n1 beg -1.7,j,-1.05 end -4.2,j,-1.05 nseg 10 sel cable id=n1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on ;自由段 gr_coh 1 gr_k 1 range cid a2,a3 sel cable id=n1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on ;锚固段 gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid a4,a5 sel cable id=n1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on ;托盘 gr_coh 10e8 gr_k 2e10 gr_per 0.0785 range cid a1,a1 sel cable id=n1 pretension 50e3 range cid a2,a3 sel cable id=n2 beg -1.7, j, -0.35 end -4.2,j,-0.35 nseg 10 sel cable id=n2 prop emod 2e10 ytension 310e3 xcarea 0.0004906 slide=on ;自由段 gr_coh 1