请教锚杆问题
[命令/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
页:
[1]