FLAC3D 施加锚杆命令流

gen zone radtun  p0 0,0,0   p1 25,0,0  p2 0,50, 0  p3  0,0,25  size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1  fill

gen  zone  reflect  normal  1 0 0  ori  0 0 0  

gen  zone  reflect  normal  0 0 1  ori  0 0 0  

mo mohr

pro   bulk  2.2e9  she  1.3e9  fric  30    coh  1.3e6  ten 1.5e5

ini dens 2000

fix x range x -25.1,-24.9

fix x range x 24.9 25.1

fix y range y 49.9 50.1

fix z range z -25.1 -24.9

fix z range z 24.9 25.1

PLOT SURF

PAUSE K

sel cable id=1 beg  0, 0, 0   end 0 ,29, 0  nseg 10 

sel cable id=1 beg  0,29,0  end 0,35,0 nseg 6

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

                    gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10

sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &

                    gr_coh 10e5 gr_k 2e7 range cid 11,17

sel delete link range id 1

sel link id=100 1 target zone

sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100

sel cable id=1 pretension 60e3 range cid 1,10

set mech ratio 7.0e-3

solve 

SAVE 2CALE.SAV

plot set plane ori 0 5 0 norm 0 0 1

plot cont syy plane ou on

plot add axe 

plot set rot 90 0 0

plot set mag 2