;岩质边坡稳定分析 vp.?$(L^@/
;---------------------------------------------------------------------- ah_ >:x
; rocksl.dat rock slope stability analysis ciiI{T[Z
; comparison to solution in Hoek & Bray, 1981 tq3Rc}
;---------------------------------------------------------------------- '21gUYm
config ats %>_6&A{K,d
grid 60,40 V>,=%r4f
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 %=Z/Frd
model elastic 'P" i 9j
; Add top surface k4hk*
0Jq
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & 9=3DYCk/
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & + xU({/
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) ;b~\ [
gen table 1 A_9J ~3
model null region 30 40 I[v`)T'_{
; ^3S&LC
1;|
fix x y j=1 W]7 /
e
fix x i=1 \q4r/SbgW
fix x i=61 .-/IV^lGv
; '
|B 3@9 <
group 'rock_mass' reg tab 1 .UNV &R0
model mohr group 'rock_mass' notnull <F(2D<d{;)
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & !U>WAD9
dilation=0.0 tension=12000.0 group 'rock_mass' notnull {>9ED.t
; 3l41"5Fy&
set gravity=9.81 FKz5,PeL
history 999 unbalanced R LMn&j|?e
solve elastic wT6zeEV~*
save rocsl_1.sav e0(aRN{W
; Add water table p r7lm5
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & Cl9nmyf
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & #v xq|$e
(140.0 46.00) ..+#~3es#y
def wet_den m%apGp'=1
loop i (1,izones) Uc'}y!R
loop j (1,jzones) KR%WBvv
if model(i,j)>1 then )RvX}y-
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) x.wDA3ys
xc=0.25*xa u:2Ll[ eo
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) zxCx2.7
yc=0.25*ya $7c,<=
if yc < table(2,xc) then Xa.8-a"hz
density(i,j) = 2800 !|i #g$
endif {, +c
endif ;H.V-~:P)
endloop Ez0zk 9
endloop
Owi/e
end ~IjID
wet_den ujS oWs
water density=1000.0 _p+E(i 9
water table=2 h=: /9O{H
solve D9A%8[Yo
save rocsl_2.sav PnaiSt9p?r
;... State: rocksl fos.fsv .... jVQ89 vf
~
solve fos file rocksl_fos.fsv kaB4[u
;*** plot commands **** w4Df? )Z
;plot name: grid #Xo x2{~
plot hold grid water iwhite G $MEVfd"
;plot name: fail Mohr h_S>Q
set pltc 140000.0 pltf 45.0 pltt 12000.0 ^ bfZd
plot hold fail normal i!e8-gVMP&
;plot name: fos - ssi : S_]!'H
plot hold fos tension block inv ssi iwhite bound vr'cR2
;plot name: fos - vel. vectors &JqaIJh
plot hold fos group water velocity plasticity no_past xk%
62W
