;岩质边坡稳定分析 Yjh02wo
;---------------------------------------------------------------------- |)b:@q3k+n
; rocksl.dat rock slope stability analysis &4jc3_UKV
; comparison to solution in Hoek & Bray, 1981 ]O+W+h{]
;---------------------------------------------------------------------- 1 'BC
R
config ats *QP+p,L *
grid 60,40 s =83a{#K
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 Vae=Yg=fw
model elastic )wfqGkr=m!
; Add top surface iJ!p9E*(
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & e>!=)6[*
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & 3s$vaV~(a
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) p [7?0 (
gen table 1 9<-7AN}Z
model null region 30 40 ,*d<hBGbh
; ExZ|_7^<
fix x y j=1 a0
8Wt
fix x i=1 +`'>
fix x i=61 \jHIjFwQ
; C=<PYkt,L
group 'rock_mass' reg tab 1 w)xfP^M#
model mohr group 'rock_mass' notnull W&;,7T8@
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & i
3i
dilation=0.0 tension=12000.0 group 'rock_mass' notnull 9?4:},FRmE
; noQS bI
@
set gravity=9.81 ,w$:=;i
history 999 unbalanced 4ZrRgx2MD
solve elastic e&ci\x%
save rocsl_1.sav rtZEK:.#
; Add water table ^#)]ICV
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & V D.T=(
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & }MW+K&sIh
(140.0 46.00) fW3NH7aUG
def wet_den xw~ 3x*{
loop i (1,izones) BFn4H%1
loop j (1,jzones) D>
E N:_v
if model(i,j)>1 then b!c2j
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) +$
0wBU
xc=0.25*xa I9O%/^5^[w
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) 4LkW`Sbm
yc=0.25*ya +9R@c Ur
if yc < table(2,xc) then zL/r V<
density(i,j) = 2800 bDT@E,cSi
endif <@J0
770
endif y.Y;<UGu
endloop H CZVvsG
endloop Q>
J9M` a
end $`6Q\=*R/
wet_den #Mrc!pT]xy
water density=1000.0 cOvdC4
water table=2 W?R@ eq.9
solve &^(4yw (~
save rocsl_2.sav i{4J$KT
;... State: rocksl fos.fsv .... X@H/"B%u2
solve fos file rocksl_fos.fsv 6QG"~>v7'(
;*** plot commands **** O7# 8g$ZIv
;plot name: grid
Ul a
h!s
plot hold grid water iwhite I7Xm~w!{qk
;plot name: fail Mohr *8I &|)x
set pltc 140000.0 pltf 45.0 pltt 12000.0 bSj-xxB]e
plot hold fail normal ISi^BFU
;plot name: fos - ssi Vl%UT@D|
plot hold fos tension block inv ssi iwhite bound ] Wx?k7T
;plot name: fos - vel. vectors (u-eL# @
plot hold fos group water velocity plasticity no_past 0artR~*}
