n
;===================================
;建立网格模型
gen zone brick &
p0 0 0 0 p1 100 0 0 p2 0 1 0 p3 0 0 20 size 25 1 5
gen zone brick &
p0 20 0 20 p1 100 0 20 p2 20 1 20 p3 60 0 60 p4 100 1 20 p5 60 1 60 p6 100 0 60 p7 100
1 60 size 20 1 10
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 99.9 101.1
fix x range x -0.1 0.1
fix y
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;===================================
;安全系数求解
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 &
coh 11000.0 tens 1.0e6 fric 20 dila 20
sel cable id=1 beg 42.7 0.5 42.7 end 52 0.5 40 nseg 10 ;自由端坐标
sel cable id=1 beg 52 0.5 40 end 67.7 0.5 35.2 nseg 6 ;锚固端
sel cable prop emod 4.5e10 ytension 60e4 xcarea 1.256e-3 &
gr_coh 1 gr_k 1 gr_per 0.314 range cid 2,10 ;自由端参数
sel cable prop emod 4.5e10 ytension 60e4 xcarea 1.256e-3 &
gr_coh 10e7 gr_k 2e9 range cid 11,17 ;锚固端参数
sel cable prop emod 4.5e14 ytension 60e4 xcarea 1.256e-3 &
gr_coh 10e10 gr_k 2e12 range cid 1,1 ;端头参数
sel cable id=1 pretension 30e5 range cid 1,10 ;自由端施加预紧力
plo con ssi out on
step 63000
solve fos file slope3dfosB.sav associated
plo fos con ssi outli on vel red
