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

;岩质边坡稳定分析 (Qk&g"I  
;---------------------------------------------------------------------- -DP8NTl"  
; rocksl.dat rock slope stability analysis As@ihB+(\  
; comparison to solution in Hoek & Bray, 1981 IXZ(]&we
 
;---------------------------------------------------------------------- d7"U WY^  
config ats Dac ^*k=D  
grid 60,40 bQwdgc),s{  
gen -40.0,-40.0 -40.0,80.0 140.0,80.0 140.0,-40.0 i=1,61 j=1,41 1C_'H.q<=  
model elastic +{xMIl_  
; Add top surface :tMre^oP  
table 1 (-130.0 -40.00) (-130.0 0) (-40.00 0) (0 0) (3.880 15.00) (11.88 & G{kj}>kS_  
15.00) (15.76 30.00) (23.76 30.00) (27.63 45.00) (43.15 60.50) (80.00 & vpm ]9>1[  
76.00) (140.0 80.00) (230.0 80.00) (230.0 -40.00) ;TL(w7vK  
gen table 1 Y
H[XRUa  
model null region 30 40 'e.q 7Jpd  
; 50LHF%  
fix x y j=1 4}Q O!(  
fix x i=1 J>=1dCK  
fix x i=61 Opv1B2  
; k42b:W5%  
group 'rock_mass' reg tab 1 +_qh)HX  
model mohr group 'rock_mass' notnull 8 7RHA $?  
prop density=2500.0 bulk=1E8 shear=3E7 cohesion=140000.0 friction=45.0 & pi`;I*f/  
dilation=0.0 tension=12000.0 group 'rock_mass' notnull 7qP4B9S  
; ~`t%M?l  
set gravity=9.81 >|a\>UgC  
history 999 unbalanced qyg*n>nt  
solve elastic 3ppuQQ  
save rocsl_1.sav
VQ`,#`wV  
; Add water table yS[z2:!  
table 2 (-40.00 0) (0 0) (3.880 4.500) (11.88 12.00) (15.76 16.00) & &/](HLdF  
(23.76 24.00) (27.64 27.00) (43.15 41.00) (52.00 43.00) (80.00 45.00) & ]3L@$`ys  
(140.0 46.00) ,IIZXl@  
def wet_den Hp\Ddx >Jd  
 loop i (1,izones) i8Fs0
U4"  
   loop j (1,jzones) V@vhj R4r\  
   if model(i,j)>1 then <ROpuY\!l  
    xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1)) Abj97S  
    xc=0.25*xa 9&jQ 35  
    ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)) jHAWK9fa  
    yc=0.25*ya $sEB'>:  
    if yc < table(2,xc) then /M3y)K`^  
     density(i,j) = 2800 ?"{QK:`  
    endif Z#-k.|}  
   endif PZys u  
  endloop `n 3FT=  
 endloop L5Urg*GNL  
end Z)m
X,=p  
wet_den >0PUWr$8  
water density=1000.0 w`L~#yu  
water table=2 f.||PH  
solve ,b4):{  
save rocsl_2.sav LthGZ|>  
;... State: rocksl fos.fsv .... S:ls[9G[3  
solve fos  file rocksl_fos.fsv }{ n\tzR  
;*** plot commands **** CI,lkO|C  
;plot name: grid i?F[||O"$  
plot hold grid water iwhite 0iYP  
;plot name: fail Mohr =~J
"kC  
set pltc 140000.0 pltf 45.0 pltt 12000.0 u4:\UC'  
plot hold fail normal 1"}B]5!  
;plot name: fos - ssi $ !v}xY  
plot hold fos tension block inv ssi iwhite bound br0u
@G  
;plot name: fos - vel. vectors 8+"10q-  
plot hold fos group water velocity plasticity no_past p?Ed- S  

  

    

什么呀，看不懂

这是源程序代码吗？

看看，这样搞是否合理？ qe$33f*  

         gWOt]D&#/  
可惜是二维的，不是三维的

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

FLAC编的吗