请教锚杆问题

li_2001

[命令/FISH] 锚杆支护锚杆, 支护锚杆, 支护

高手 您好 您编写的是全锚索支护的情况 一根锚索不需分段  现在我正在做一个有托盘的锚索  可是在运行过程中老是出现soure node not near enough to an exsiting non-null zone 现在将程序传上来望您指点 谢谢 new gen zone brick size 10 10 10 & p0 0 0 0 p1 10 0 0 p3 0 0 10 model mohr prop bulk 1e8 shear 0.3e8 fric 35 coh 1e10 ten 1e10 set grav 0 0 -9.81 ini dens 1000 fix x range x 0 fix x range x 10 fix y range y 0 fix y range y 10 fix z range z 0 hist unbal set hist_rep 10 solve save 初始平衡.sav ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 model null range x 0 2 y 0 2 z 0 2 solve save 1次开挖.sav restore 初始平衡.sav ini xdisp 0 ydisp 0 zdisp 0 ini xvel 0 yvel 0 zvel 0 model null range x 0 2 y 0 2 z 0 2 def place_cables     loop m(1,2)     caseof m     case 1     x_d1=0.5     x_d2=0.5     x_d3=0.5     y_d1=0.5     y_d2=0.5     y_d3=0.5     z_d1=1.5     z_d2=2.0     z_d3=3.5        a=1        b=2        c=3        d=6        id_11=1        id_111=1     case 2     x_d1=0.5     x_d2=0.5     x_d3=0.5     y_d1=1.5     y_d2=1.5     y_d3=1.5     z_d1=1.5     z_d2=2.0     z_d3=3.5        a=15        b=16        c=17        d=20        id_11=3        id_111=15    endcase    loop n(1,2)    x_d1=x_d1+1    x_d2=x_d2+1    x_d3=x_d3+1    y_d1=y_d1    y_d2=y_d2    y_d3=y_d3    z_d1=z_d1    z_d2=z_d2    z_d3=z_d3        a=a+(n-1)*7        b=b+(n-1)*7        c=c+(n-1)*7        d=d+(n-1)*7         id_11=id_11+1        id_111=id_111+(n-1)*7   command   sel cable id=id_11 beg x_d1 y_d1 z_d1 end x_d2 y_d2 z_d2 nseg 2   sel cable id=id_11 beg x_d2 y_d2 z_d2 end x_d3 y_d3 z_d3 nseg 4   sel cable id=id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &    gr_coh 1 gr_k 1 gr_per 0.0785 range cid a,b   sel cable id= id_11 prop emod 2e10 ytension 3.1e5 xcarea 0.0004906 &    gr_coh 10e5 gr_k 2e7 gr_per 0.0785 range cid c,d   sel delete link range id id_111   sel link id=id_111 id_111 target zone   sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id id_111   sel cable id=id_11 pretension 60e3 range cid a,d   endcommand   end_loop   end_loop   end   place_cables