命令如下：
new
gen z radt p0 (0,0,0) p1 20,0,0 p2 0 0.5 0 p3 0 0 20 &
size 8 4 8 18 ratio 1 1 1 1.1 dim 2 2 2 2 group 1
gen z b p0 2 0 0 p1 2 0 2 p2 2 0.5 0 p3 0 0 0 size 8 4 8 group 2
pause
gen z reflect norm (1,0 0) origin = (0 0 0 )
gen z reflect norm (0,0,1) origin =(0 0 0 )
attach face all
save 模型m.sav
rest 模型m.sav
model mohr
prop bulk 26.8e9 shear 7.0e9 coh 37.2e6 fric 27.8 ten 2.6e7
apply szz=-20e6 range z= 19.9 20.1
init sxx=-20e6
;fix x range x= -20.1 -19.9
;fix x range x= 19.9 20.1
fix y range y= -0.01 0.01
fix y range y= 0.499 0.501
fix z rang z=-20.1 -19.9
set large
solve
save 原岩m.sav ;原岩
new
rest 原岩m.sav
;开挖
ini xdis=0 ydis=0 zdis=0
prop bulk 26.8e9 shear 7.8e9 coh 27.2e6 fric 27.8 ten 2.17e6
mod null range group 2
pause
;加锚杆
;定义顶锚杆
sel cable id 1 begin -2 0.25 -1 end -4.5 0.25 -1 nseg=10;左帮第一根锚杆
sel cable id 1 begin -2 0.25 0 end -4.5 0.25 0 nseg=10;左帮第二根锚杆
sel cable id 1 begin -2 0.25 1 end -4.5 0.25 1 nseg=10;左帮第三根锚杆
sel cable id 1 begin -1 0.25 2 end -1.0 0.25 4.5 nseg=10;顶板第四根锚杆
sel cable id 1 begin 0 0.25 2 end 0 0.25 4.5 nseg=10;顶板第五根锚杆
sel cable id 1 begin 1 0.25 2 end 1 0.25 4.5 nseg=10;顶板第六根锚杆
sel cable id 1 begin 2 0.25 1 end 4.5 0.25 1 nseg=10;右帮第七根锚杆
sel cable id 1 begin 2 0.25 0 end 4.5 0.25 0 nseg=10;右板第八根锚杆
sel cable id 1 begin 2 0.25 -1 end 4.5 0.25 -1 nseg=10;右板第九根锚杆
;给锚杆、锚索赋参数===========================================
sel cable id 1 prop emod 200e9 ytens 190e7 xcarea 3.14e-4 &
gr_k 5.35e9 gr_coh 4.2e5 gr_per 0.75
;建立锚杆连接link
def camd
new_id=150
loop k(0,8)
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 100e3
end_command
end_loop
end
camd
pause
set large
solve
save 锚杆支护m.sav
总是提示：source node not near enough to a existing non-null zone !求高手指教。
