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

;岩质边坡稳定分析 VY "i>Ae  
;---------------------------------------------------------------------- CI@qT}Y_  
; rocksl.dat rock slope stability analysis m"mU:-jk`  
; comparison to solution in Hoek & Bray, 1981 {_]<mwd  
;---------------------------------------------------------------------- ,%,.c^-  
config ats o`\@Yq$.  
grid 60,40 %s[ n2w  
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 axdRV1+s  
model elastic ]f3R;d  
; Add top surface 2-gI@8NPI  
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & 1)ZdkTF@H  
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & TRQH{O\O  
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) jLreN#:9  
gen table 1 X8(Ws
N  
model null region 30 40 "$|ne[b2  
; @'FOM  
fix x y j=1 CL"q"  
fix x i=1 /7Ft1f  
fix x i=61 (W_U<~`t  
; be~'}`>  
group 'rock_mass' reg tab 1 &(rR
)cG  
model mohr group 'rock_mass' notnull Bc51 0I$c  
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & z229:L6"  
dilation=0.0 tension=12000.0 group 'rock_mass' notnull <84d Vg  
; w&LL-~KI+  
set gravity=9.81 1#^r5E4  
history 999 unbalanced HH'5kE0;d  
solve elastic n}4Lq^$  
save rocsl_1.sav 3+iQct[  
; Add water table _u8d`7$*%  
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & S$i3/t  
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & O:W4W=K  
(140.0 46.00) ,98`tB0  
def wet_den d#
q8-  
 loop i (1,izones) 
C9E@$4*  
   loop j (1,jzones) &BQ%df<y\  
   if model(i,j)>1 then Ozs&YZ  
    xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) IsP!ZcV;  
    xc=0.25*xa }B2qtb3  
    ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) ph=U<D4  
    yc=0.25*ya |BA<> WE  
    if yc < table(2,xc) then jW_FaPW(p  
     density(i,j) = 2800 Dhp|%_>  
    endif Mt+ggF.  
   endif pc/]t^]p  
  endloop \FjY;rqfKe  
 endloop |%5nV=&\  
end ;.b^A  
wet_den %1e{"_$O9  
water density=1000.0 (Kaunp5_`  
water table=2 :faB7wduW;  
solve rqYx\i?  
save rocsl_2.sav -LEpT$v|  
;... State: rocksl fos.fsv .... !!UQ,yU  
solve fos  file rocksl_fos.fsv 2f.4P]s`T  
;*** plot commands **** w"BMJ+  
;plot name: grid u YJL^I8M'  
plot hold grid water iwhite ah0  
;plot name: fail Mohr +v"%@lC};  
set pltc 140000.0 pltf 45.0 pltt 12000.0 Ixv/xI  
plot hold fail normal q<wQ/m  
;plot name: fos - ssi -gb'DN1BG  
plot hold fos tension block inv ssi iwhite bound qn~:B7f  
;plot name: fos - vel. vectors [Xo}CU  
plot hold fos group water velocity plasticity no_past 5`[B:<E4  

  

    

什么呀，看不懂

这是源程序代码吗？

看看，这样搞是否合理？ @O%d2bgEWV  

         37QXML  
可惜是二维的，不是三维的

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

FLAC编的吗