如何用fish语言编写巷道锚杆支护?
刚学了点FLAC,只会一根一根的加锚杆,太麻烦了。有哪个兄弟有fish语言编的锚杆支护程序?能给个实例看看吗? 我也是一根一根的加的,应该可以编写,好像手册上有这方面的例子,我也是勉强能计算,大家慢慢的学习!! 就是个loop循环而已,自己研究一下嘛 我还是个门外汉,自己研究进步慢啊,大牛们指导一下吧 http://www.simwe.com/forum/viewthread.php?tid=828765&extra=page%3D2%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D36这是我写的一个循环,希望对你有帮助! 学习一下,呵呵 学习一下,呵呵 呵呵呵我在看书时也是看到一段段开挖,一根一根的支护啊 学习学习 学习学习 学习一下 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 请 zhuimeng317 帮忙解答 请 zhuimeng317 帮忙解答
页:
[1]
2