岩质边坡稳定分析 [复制链接]

;岩质边坡稳定分析 =8F]cW'1`  
;---------------------------------------------------------------------- oAgO3x   
; rocksl.dat rock slope stability analysis /"M7YPX;  
; comparison to solution in Hoek & Bray, 1981 |P"p/iY  
;---------------------------------------------------------------------- KvrcO#-sL  
config ats ch%-Cg~%  
grid 60,40 qs]7S^yw  
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 ?u9JRXj%  
model elastic &_EjP hZ  
; Add top surface ,\BGxGNAmV  
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & Uxq9H  
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & L[Dr[  
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) xN +j]LC  
gen table 1 }I\hOL
 
model null region 30 40 P'4j
z&4  
; B>CG/]  
fix x y j=1 k2xOu9ncEj  
fix x i=1 "Oh(&N:U  
fix x i=61 Lit@ m2{\  
; prC1<rm  
group 'rock_mass' reg tab 1 'o#ve72z1  
model mohr group 'rock_mass' notnull + W ? /A]  
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & rN'')n/F  
dilation=0.0 tension=12000.0 group 'rock_mass' notnull X~O2!F  
; Ifm|_  
set gravity=9.81 4q`$nI Bi  
history 999 unbalanced 2U-F}Z  
solve elastic 1R1DK$^c  
save rocsl_1.sav eJB !|  
; Add water table =?} t7}#  
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & d5zF9;[  
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & "YUyM5X  
(140.0 46.00) +@Oo)#V|.  
def wet_den Kzw)Q  
 loop i (1,izones) \k{UqU+s  
   loop j (1,jzones) d85\GEF9i  
   if model(i,j)>1 then 2\xv Yf-  
    xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) *;~*S4/P   
    xc=0.25*xa &J)q_Z8  
    ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) ?G[<~J3-E  
    yc=0.25*ya |=7ouFl  
    if yc < table(2,xc) then 'Jb6CRn  
     density(i,j) = 2800 !H`Q^Xf}  
    endif qYHAXc}$  
   endif =|lK
B;  
  endloop OIK14D:  
 endloop "JLKO${ Y  
end $td=h)S^`  
wet_den >g%^hjJ  
water density=1000.0 ~:Jw2 P2z
 
water table=2 E$-u:Z<-  
solve "#H@d+u  
save rocsl_2.sav 47R4gs#W  
;... State: rocksl fos.fsv .... I*/?*p/I  
solve fos  file rocksl_fos.fsv ]&za^%q0&  
;*** plot commands **** bwjjwu&  
;plot name: grid !c`KzqP  
plot hold grid water iwhite 7Vf2Qx1_  
;plot name: fail Mohr i2ap]  
set pltc 140000.0 pltf 45.0 pltt 12000.0 fUh7PF%
 
plot hold fail normal xkIRI1*!  
;plot name: fos - ssi D|W^PR:@h  
plot hold fos tension block inv ssi iwhite bound 8SvPDGu`]  
;plot name: fos - vel. vectors g">E it*[  
plot hold fos group water velocity plasticity no_past IL?mt2IQ>  

  

    

什么呀，看不懂

这是源程序代码吗？

看看，这样搞是否合理？ w'Y(doY,  

         %fT%,( w}t  
可惜是二维的，不是三维的

好资料  需要这方面的资料谢谢分享

FLAC编的吗