Flac3D中一个锚索例子
new
；建立网格（倾角1：1.5，边坡前缘距坡脚7.5m，边坡后缘距坡顶10m，前、后缘相距40m，深5m，坡体左右截面间距10m，下卧基岩坡度为1：2.5）
gen zone brick p0 0 0 0 p1 37.5 0 0 p2 0 3 0 p3 0 0 5 size 20 3 4 ratio 1 1 1
gen zone brick p0 37.5 0 0 p1 45 0 0 p2 37.5 3 0 p3 37.5 0 5 size 4 3 4 ratio 1 1 1
gen zone brick p0 0 0 5 p1 37.5 0 5 p2 0 3 5 p3 0 0 20 p4 37.5 3 5 p5 0 3 20 p6 15 0 20
p7 15 3 20 size 20 3 12 ratio 1 1 1
；plo block model id on
；分组
group hnt red range plane above ori 5 0 16 nor 2 0 5
group bedrock green range group hnt not
;;;;;;设置边界条件
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 44.9 45.1
fix y range y -0.1 0.1
fix y range y 2.9 3.1
;;;;;;;生成初始应力
m e定义弹性模型
prop bulk 3e7 shear 1e7
prop density 1970
set grav 0 0 -9.81
solve
save mg1.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;;;;;;;;安全系数求解（采用关联流动法则:fi=dila）
m m
ini dens 1970 range group hnt
prop bulk 1.15E7 shear 5.31E6 coh 20000 fri 12.1 ten 1000 dila 0 range group hnt
ini dens 2000 range group bedrock
prop bulk 22.6E9 shear 11.1E9 coh 6.72E6 fri 42 ten 1.58E6 dila 10.0 range group bedrock
sel cable id=1 begin=(20.63,1.5,16.24) end=(10.11,1.5,13.03) nseg=12
sel cable id=1 begin=(10.11,1.5,13.03) end=(5.33,1.5,11.57) nseg=4
sel cable id=1 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 1,12
sel cable id=1 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 13,16
sel delete link range id=1
sel link id=1 1 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 1
sel cable id=1 pretension 1e5 range cid 1 12
sel cable id=2 begin=(26.25,1.5,12.5) end=(15.73,1.5,9.28) nseg=10
sel cable id=2 begin=(15.73,1.5,9.28) end=(10.95,1.5,7.82) nseg=4
sel cable id=2 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 17,26
sel cable id=2 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 27,30
sel delete link range id=19
sel link id=19 18 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 19
sel cable id=2 pretension 1e5 range cid 17 26
sel cable id=3 begin=(31.88,1.5,8.74) end=(21.36,1.5,5.53) nseg=8
sel cable id=3 begin=(21.36,1.5,5.53) end=(16.58,1.5,4.07) nseg=3
sel cable id=3 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 0 gr_coh = 0 gr_fric 0 gr_per = 0.314 range cid 31,38
sel cable id=3 prop xcarea = 1.256e-3 emod = 45e9 ytension = 60e4 gr_k = 2.75e7 gr_coh = 2.2e4 gr_fric = 20 gr_per = 0.314 range cid 39,41
sel delete link range id=35
sel link id=35 33 target zone
sel link attach xdir = rigid ydir = rigid zdir = rigid xrdir = rigid yrdir = rigid zrdir = rigid range id 35
sel cable id=3 pretension 1e5 range cid 31 38
plo sk sel geo id on cid on node on scale 0.005
solve fos file mg2.sav associated ;;;;;;Fs=1.37
plot fos con ssi outline on vel red ;;;;;;绘制安全系数，剪应变增量及速度矢量图
;;;;;;;;;;监测采样
hist id=6 gp disp 20.63,1.5,16.24
hist id=7 gp disp 26.25,1.5,12.5
hist id=8 gp disp 31.88,1.5,8.74
plo hist 6 7 8
solve
save mg3.sav