;岩质边坡稳定分析 ? { F{;r
;---------------------------------------------------------------------- pY^9l3y^
; rocksl.dat rock slope stability analysis N <_Ko+VF
; comparison to solution in Hoek & Bray, 1981 y9;#1:i c
;---------------------------------------------------------------------- f.pkQe(
config ats j%*7feSNC
grid 60,40 Y &r]lD
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 :PnSQjV:
model elastic 4R ) |->"
; Add top surface S HqyvF
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & ;ggy5? >Qu
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & M\ +* P,i
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) I/COq U 7~
gen table 1 Y5!b)vke
model null region 30 40 [ij,RE7,T
; g "!\\:M
fix x y j=1 jV7&Y.$zF]
fix x i=1 qM S}t3X
fix x i=61 qG >DTKIU
; D\T!4q'Q
group 'rock_mass' reg tab 1 'CS.p!Z\
model mohr group 'rock_mass' notnull HrqF![_
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & ]Bb7(JX
dilation=0.0 tension=12000.0 group 'rock_mass' notnull 3On
JWuVfZ
; 5>ST"l_ca
set gravity=9.81 76#.F
history 999 unbalanced 7-9;PkGG.A
solve elastic VKLU0*2R
save rocsl_1.sav . )8
; Add water table ]~a_d)
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & ^^$vR[7
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & O
rk
(140.0 46.00) bEr.nF
def wet_den s5Fr)q// !
loop i (1,izones) ?: |YGLaB
loop j (1,jzones) ":d*dl
if model(i,j)>1 then jgvh[@uB?
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) hD 46@
xc=0.25*xa G79C {|c\
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) K a r~I
yc=0.25*ya Wm6dQQ;Bj
if yc < table(2,xc) then A:Rw@ B$
density(i,j) = 2800 ~Y/z=^
endif ,p, Du
F
endif dB|Te "6
endloop \kwe51MQ
endloop %9c|%#3
end ^)cM&Bx t%
wet_den =
?N^>zie
water density=1000.0 -Q`C q |s
water table=2 q?VVYZXP
solve B {i&~k
save rocsl_2.sav iiwpSGFl]
;... State: rocksl fos.fsv .... ye r>
x
solve fos file rocksl_fos.fsv NFoZ4R1gy
;*** plot commands **** FH n,]Tfx
;plot name: grid owMuT^x?
plot hold grid water iwhite W:ixzpQ
;plot name: fail Mohr wd`R4CKhP]
set pltc 140000.0 pltf 45.0 pltt 12000.0 BpXEK.Xw
plot hold fail normal S.fXHt Sx
;plot name: fos - ssi 1XnBK$`
plot hold fos tension block inv ssi iwhite bound S_!R^^ySG9
;plot name: fos - vel. vectors Kpu<rKP`
plot hold fos group water velocity plasticity no_past tq E>Zx=X
