;岩质边坡稳定分析 yuA+YZ
;---------------------------------------------------------------------- eF=cMC
; rocksl.dat rock slope stability analysis Ww)p&don
; comparison to solution in Hoek & Bray, 1981 s. ]<r5v7
;---------------------------------------------------------------------- DITo.PU
config ats l+HF+v $
grid 60,40 sjIUW$
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 9?J
3G,&
model elastic Y5,[udF:O
; Add top surface &Ay[mZQ 7
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & g*]/HS>e<G
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & J VxGS{Z
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) b;k3B7<
gen table 1 OPe3p {]
model null region 30 40 i&_& 4
; =R6IW, *
fix x y j=1 2*)2c[/0F
fix x i=1 vnr{Ekg
fix x i=61 bbe$6x wi
; ta+MH,
group 'rock_mass' reg tab 1 0+dc
model mohr group 'rock_mass' notnull =&~7Q"
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & [ ~&yLccN
dilation=0.0 tension=12000.0 group 'rock_mass' notnull OAw- -rl
; "kS!rJ[
set gravity=9.81 T0TgV
history 999 unbalanced 3[F9qDAy
solve elastic 1O ,:fTG<
save rocsl_1.sav l.nd Wv
; Add water table CyXFuk!R
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & E>I\m!ue
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & hb ="J349
(140.0 46.00) '
Ih f|;r
def wet_den 0Fc^c[
loop i (1,izones) }N,$4h9D j
loop j (1,jzones) {'@`: p&3r
if model(i,j)>1 then $s]vZ(H
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) scQnL '\
xc=0.25*xa ?Dr A@;IB
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) :tf'Gw6v
yc=0.25*ya ?[VL
2dP0
if yc < table(2,xc) then YM r2|VEU[
density(i,j) = 2800 +S 6(Fvp
endif !fmbm4 !a
endif A~) #
endloop !*G%vOa
endloop ;1AX u/
end r7^oqEp@B
wet_den QCAoL.v
water density=1000.0 Mj<T+Ohz
water table=2 'B\7P*L"p
solve riv8qg
save rocsl_2.sav fvBL? x
;... State: rocksl fos.fsv .... G0~6A@>
solve fos file rocksl_fos.fsv 38E
%]*5F
;*** plot commands **** O;lGh1.
;plot name: grid Rl{e<>O\^
plot hold grid water iwhite Lvd es.0|
;plot name: fail Mohr cXod43
set pltc 140000.0 pltf 45.0 pltt 12000.0 Us2I eR
plot hold fail normal 7SJR_G6,{
;plot name: fos - ssi lW|`8ykp
plot hold fos tension block inv ssi iwhite bound C C
B'
;plot name: fos - vel. vectors %urd;h D
plot hold fos group water velocity plasticity no_past Zk ] / m
