;岩质边坡稳定分析 0-dh Gh?.
;---------------------------------------------------------------------- m
.2)P~a
; rocksl.dat rock slope stability analysis }<Y3 jQnl
; comparison to solution in Hoek & Bray, 1981 V%~u8 b
;---------------------------------------------------------------------- kTQ`$V(>&
config ats f#xqu +)Z
grid 60,40 'ad|@Bh
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 xe)< )y
model elastic h%kB>E~
; Add top surface wzAp`Zs2Dm
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 &
5" U8|
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & 7S<Z&1(
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) ^0t81,`
gen table 1 h#bpog
model null region 30 40 E.Hw|y0_(|
; 1a {~B#
fix x y j=1 C3`2{1
fix x i=1 C._I\:G^
fix x i=61 -CW$p=y}
; h~$Q\WCm#
group 'rock_mass' reg tab 1 X/,4hjg
model mohr group 'rock_mass' notnull @vf{_g<
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & %v++AcE
dilation=0.0 tension=12000.0 group 'rock_mass' notnull 7Kx3G{5ja
; xBGSj[1`i
set gravity=9.81 GiBq1U-Q
history 999 unbalanced ,wmPK;j
solve elastic hk"^3d!
save rocsl_1.sav GXaCH))TO
; Add water table i,/|H]Mzr
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & B^(0>Da\
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & KZV$rJ%G
(140.0 46.00) (tGK~!cAv
def wet_den i>bFQ1Rdx
loop i (1,izones) 7`3he8@ze
loop j (1,jzones) $jb3#Rj4
if model(i,j)>1 then BaIh, iu
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) A0*u(15%
xc=0.25*xa ["N>Po
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) ]2Aqqy
yc=0.25*ya ilv6A9/
if yc < table(2,xc) then ;F@dN,Y
density(i,j) = 2800 Vxif0Bx&/d
endif , cm;A'4]
endif bHcb.;<
endloop DBi3 j
endloop r]1|I6:&)
end * wN+Ak q
wet_den g<~[k?~J
water density=1000.0 UP:+1Sp9
water table=2 hB:R8Y^?H
solve
/A|cO
save rocsl_2.sav Fs:l"5~>1
;... State: rocksl fos.fsv .... tq9t(0EL
solve fos file rocksl_fos.fsv kKVq,41'
;*** plot commands **** >u5}5OP7
;plot name: grid 2S^xqvh
plot hold grid water iwhite ~ S Js2- 2
;plot name: fail Mohr fU~>A-P
set pltc 140000.0 pltf 45.0 pltt 12000.0 di6A.N5A
plot hold fail normal 5v8&C2Jy@
;plot name: fos - ssi %j=E}J<H5*
plot hold fos tension block inv ssi iwhite bound Ch
` Omq
;plot name: fos - vel. vectors c Xcn}gKV
plot hold fos group water velocity plasticity no_past <_""4
