n
title
土钉（用锚索）
gen zone b size 11 1 11
model mohr
prop bulk 5e9 shear 1e9 coh 4e4 fric 30
ini dens 2000
fix x z range x -0.1 0.1
fix x range x -0.1 0.1
fix x range x 10.9 11.1
fix y
;
ini szz=-2.2e5 grad 0,0,2e4
ini sxx=-1.32e5 grad 0,0,1.2e4
ini syy=-0.88e5 grad 0,0,0.8e4
set grav 0,0,-10
;
hist unbal
hist gp xdisp 0,0,11
;
solve
;
free x range x -0.1 0.1
;
sel cable id=1 begin 0.0 0.5 3.5 end 8 0.5 3.5 nseg=8
sel cable id=2 begin 0.0 0.5 6.5 end 8 0.5 6.5 nseg=8
sel cable id=3 begin 0.0 0.5 9.5 end 8 0.5 9.5 nseg=8
sel cable prop xcarea=8.5e-3 emod=200e9 ytens=1e10 &
gr_k=7e6 gr_coh=1e2
;sel cable prop gr_per=0.314 gr_fric=25
;
plot create snail
plot set cap size 25
plot add sketch
plot add disp max=2e-2
plot add sel cable force yellow red
plot add sel cable grout slip
plot add sel geom node=off black black
plot show
;
step 5000
;save cable-sn.save
plot show
到 solve 就计算不收敛了，不知道后面的结构他是怎么算出来的
