new
gen zone brick p0 0 0 0 p1 50 0 0 p2 0 25 0 p3 0 0 50 size 50 25 50
mo mohr range z 0 50
prop bulk 2615384615 shear 2903225806 fri 30 coh 1.00E+06 ten 1.00E+06 range x 0 50 y 0 50 z 0 50 ; 38-砂质泥岩
set g 0 0 -10
ini dens 2630
fix x range x -0.1 0.1
fix x range x 49.9 50.1
fix y range y -0.1 0.1
fix y range y 24.9 25.1
fix x y z range z -0.1 0.1
apply szz -6350000 range z 49.9 50.1
ini sxx -3832500 syy -3832500 szz -7665000 grad 0 0 26300 range x 0 50 y 0 50 z 0 50 ; 38-砂质泥岩
hist unbal
plot con szz gr
solve
save 初始平衡.sav
;开挖
model null range x 10 40 y 0 50 z 18 20
define place_cable
loop k(1,10)
y_a=2.5*float(k)+1
x_a=12*float(k-1)
command
sel cable beg 25 y_a 20 end 25 y_a 21.29 nseg 2 ;建立自由段
sel cable beg 25 y_a 21.29 end 25 y_a 22.4 nseg 2 ;建立锚固段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+00 gr_k 1 gr_per 0.0785 range cid=(x_a+2) ;自由段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.000491 gr_coh 1.00E+06 gr_k 2.00E+07 range cid x_a+3 x_a+4 ;锚固段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+09 gr_k 2.00E+10 range cid 1 1 ;模拟托盘
sel cable pretension 6.00E+04 range cid x_a+1 x_a+2
sel cable beg 15 y_a 20 end 15 y_a 21.29 nseg 2
sel cable beg 15 y_a 21.29 end 15 y_a 22.4 nseg 2
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+00 gr_k 1 gr_per 0.0785 range cid x_a+6 ;自由段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.000491 gr_coh 1.00E+06 gr_k 2.00E+07 range cid x_a+7 x_a+8 ;锚固段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+09 gr_k 2.00E+10 range cid x_a+5 x_a+5 ;模拟托盘
sel cable pretension 6.00E+04 range cid x_a+5 x_a+6
sel cable beg 35 y_a 20 end 35 y_a 21.29 nseg 2
sel cable beg 35 y_a 21.29 end 35 y_a 22.4 nseg 2
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+00 gr_k 1 gr_per 0.0785 range cid x_a+10 ;自由段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.000491 gr_coh 1.00E+06 gr_k 2.00E+07 range cid x_a+11 x_a+12 ;锚固段
sel cable prop emod 2.00E+10 ytension 3.10E+05 xcarea 0.0004906 gr_coh 1.00E+09 gr_k 2.00E+10 range cid x_a+9 ;模拟托盘
sel cable pretension 6.00E+04 range cid x_a+9 x_a+10
end_command
end_loop
end
place_cable
plot add sel
plot sk sel cable force red
pause
solve
save 打一列锚杆平衡.sav
