;岩质边坡稳定分析 _+PiaJ&'
;---------------------------------------------------------------------- Klfg:q:j+b
; rocksl.dat rock slope stability analysis rD=8O#m
g
; comparison to solution in Hoek & Bray, 1981 YoyJnl.?u
;---------------------------------------------------------------------- %B?@le+%
config ats Qt` hUyL
grid 60,40 !>QS746S@
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 JPQWRK^
model elastic f%vJmpg
; Add top surface _:@~bHd
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & \7rAQ[\#V
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & d 90
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) HVO
mM17
gen table 1 {gh41G;n
model null region 30 40 X`#,*HkK
; ^X-3YhJ4U
fix x y j=1 $A~UA
fix x i=1 e_'/4
n
fix x i=61 o6`4y^Q{/
; J{>9ctN
group 'rock_mass' reg tab 1 HYVSi3[
model mohr group 'rock_mass' notnull ^u'hl$`^
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & (NfP2E|B
dilation=0.0 tension=12000.0 group 'rock_mass' notnull BGM5pc (ei
; _vOV(#q2a
set gravity=9.81 @1+C*
history 999 unbalanced )Z6bMAb0'N
solve elastic aD)$aK
save rocsl_1.sav 7$3R}=Z`\q
; Add water table C]S~DK1
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & $t*>A+J
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & *R9s0;&:
(140.0 46.00) Bmx+QO
def wet_den Mdk(FG(
loop i (1,izones) U'(Exr[
loop j (1,jzones) vW6
a=j8
if model(i,j)>1 then U@t"o3E
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) D5]AL5=Xt2
xc=0.25*xa b"D? @dGB,
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) k^<s|8Y
yc=0.25*ya pe^hOzVv
if yc < table(2,xc) then 5L ]TV\\
density(i,j) = 2800 `Fn"%P!
endif H~Z$ pk%
endif j
/d?c5
endloop (<M^C>pldf
endloop E}Ir<\
end (pBPf
wet_den '8%pEl^
water density=1000.0 0*e)_l!
water table=2 G.j R
solve @|I:A
save rocsl_2.sav $H0diwl9R
;... State: rocksl fos.fsv .... PK&3nXF%4
solve fos file rocksl_fos.fsv By3y.}'Ub9
;*** plot commands **** OGl>i
;plot name: grid mT]+wi&
plot hold grid water iwhite 3k%fY
;plot name: fail Mohr Mn)>G36(
set pltc 140000.0 pltf 45.0 pltt 12000.0 B #o/3
plot hold fail normal x"N{5
;plot name: fos - ssi ~*R:UTBtw
plot hold fos tension block inv ssi iwhite bound >Q(+H-w
;plot name: fos - vel. vectors aJ!(c}N~97
plot hold fos group water velocity plasticity no_past eA4dDKX+