souce node not near enough to an existing non-null zone
newrest无支护开挖.sav
title '支护方案’
plot set title text 'fangan'
set log on
set logfile wuzhihufangan.txt
;赋值新的参数
;define aing
;command
;泥岩
prop bulk=9.1e9 shear=4.7e9coh=1.1e7tens=1.5e7 fric=33.7 range group 1
prop bulk=10.6e9 shear=5.5e9coh=1.1e7tens=1.6e7 fric=33.7 range group 3
prop bulk=13.9e9 shear=5.7e9coh=1.1e7tens=1.3e7 fric=33.7 range group 4
prop bulk=13.9e9 shear=5.7e9coh=1.1e7tens=1.3e7 fric=33.7 range group 6
prop bulk=29.5e9 shear=15.2e9coh=1.1e7tens=2.38e8 fric=33.7 range group 8
prop bulk=22.7e9 shear=11.7e9coh=1.1e7tens=1.8e8 fric=33.7 range group 9
prop bulk=19.7e9 shear=10.2e9coh=1.1e7tens=2.0e6 fric=33.7 range group 10
;9#煤层
prop bulk=2.2e9shear=0.7e9coh=5e6tens=8e6 fric=35 range group 2
;8#煤层下
prop bulk=2.2e9shear=0.8e9coh=5e6tens=9e6 fric=35 range group 5
;8#煤层上
prop bulk=2.2e9shear=0.7e9coh=5e6tens=10e6 fric=35 range group 7
;初始条件
initial density 2203 range group 1
initial density 1274 range group 2
initial density 2589 range group 3
initial density 2476 range group 4
initial density 1321 range group 5
initial density 2321 range group 6
initial density 1321 range group 7
initial density 2521 range group 8
initial density 2421 range group 9
initial density 2321 range group 10
;end_command
;end
;fuzhi
;重新设置初始位移
initial xdisp=0 ydisp=0 zdisp=0
;开挖
model null range x 14 18.2 y 0 10 z 14 16.8
;加锚杆
define ins_cab
loop i(0,6)
yi=0.2+i*1.6
command
;左帮锚杆
sel cable id 1 begin 13.99 yi 14.6 end 11.99 yi 14.6 nseg=10;左帮第一根锚杆
sel cable id 1 begin 13.99 yi 16.2 end 11.99 yi 16.2 nseg=10;左帮第二根锚杆
;右帮锚杆
sel cable id 1 begin 18.21 yi 14.6 end 20.21 yi 14.6 nseg=10;右帮第一根锚杆
sel cable id 1 begin 18.21 yi 16.2 end 20.21 yi 16.2 nseg=10;右帮第二根锚杆
;顶锚杆
sel cable id 1 begin 14.6 yi 16.81 end 14.6 yi 18.81 nseg=10;顶板第一根锚杆
sel cable id 1 begin 16.1 yi 16.81 end 16.1 yi 18.81 nseg=10;顶板第二根锚杆
sel cable id 1 begin 17.6 yi 16.81 end 17.6 yi 18.81 nseg=10;顶板第三根锚杆
end_command
end_loop
loop n(0,1)
yn=1+n*4.8
command
;顶板锚索
sel cable id 2 begin 16.1 yn 16.81 end 16.1 yn 22.11 nseg=10
end_command
end_loop
end
ins_cab
;给锚杆、锚索赋参数
define maogan
loop k(0,48)
id0=1+k*11
id1=5+k*11
id2=10+k*11
command
sel cable id 1 prop emod=200e9 ytens=100e7 xcarea=2.54e-4 &
gr_k=0 gr_coh=0 gr_per=0.088 range cid id0 id1
sel cable id 1 prop emod=200e9 ytens=100e7 xcarea=2.54e-4 &
gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.088 range cid id1 id2
end_command
end_loop
loop n(0,1)
id3=1+k*11
id4=8+k*11
id5=11+k*11
command
sel cable id 2 prop emod=200e9 ytens=186e7 xcarea=1.82e-4 &
gr_k=0 gr_coh=0 gr_per=0.088 range cid id3 id4
sel cable id 2 prop emod=200e9 ytens=186e7 xcarea=1.82e-4 &
gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.088 range cid id4 id5
end_command
end_loop
end
maogan
;建立锚杆连接link
def camd
new_id=10000
loop k(0,48)
id_node=1+k*11
new_id=new_id+1
command
sel delete link range id id_node
sel link id=new_id id_node target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id new_id
sel cable id=1 pretension 50e3
end_command
end_loop
end
camd
;建立锚索连接link
def cdma
new_id=20000
loop k(0,1)
id_node=540+k*11
new_id=new_id+1
command
sel delete link range id id_node
sel link id=new_id id_node target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id new_id
sel cable id=2 pretension 120e3
end_command
end_loop
end
cdma
;计算求解
solve
save 支护方案.sav
set log off
这是命令？还请各位指教，新手
页:
[1]