FLAC3D滑坡稳定性分析，静力平衡和强度折减法实例

ajianme

那位前辈曾做过FLAC3D模拟边坡稳定，能否给我发一份建模及分析计算的全过程（邮箱34437751@qq.com），wo 看了几本书和一些相关资料,.自己也做了实例，可还是不清楚具体要怎么弄，望高手给予指点，十分感激。
ef _geoparm
size1=150      
size2=200     
filecal0=string(size1)+'_'+string(size2)+'_cal'+'.sav'  

bf=30      
af=bf*pi/180
h1=20        


gen zone brick size 150,200, 20  p0 0 ,0 ,0 &
p1 34.5,100,0 p2 0,2 ,0 &
p3 0 0 25 p4 34.6,100,0 &
p5 0,2,25 p6 34.6,-100,5 p7 34.6,100,5 size size1 1 size2
plot surf

def SSR


ait1=0.02         
k11=1.00         
k12=2.50         
ks=(k11+k12)/2

loop while (k12-k11)>ait1


flag=0
  E1=2.39e8  
  E2=7.38e9

  poi1=0.35
  poi2=0.22  

  coh1=1.7e4/ks
coh1=3.65e5/ks

  weight1=1.53e3
weight1=2.72e3

  dila1=0.0

  
  fri1=24.2
fri2=22.7
  
  
  ten1=0.05e6
  grav0=-9.80           
  dens1=-weight1/grav0   
  K1=E1/(3*(1-2*poi1))
  G1=E1/(2*(1+poi1))  
K2=E2/(3*(1-2*poi2))
  G2=E2/(2*(1+poi2))

command

model null      
model mohr  
pro bulk K1 she G1 dens dens1 coh coh1 &
friction fri1 dilation 0. tens ten1

fix x range x -0.1 0.1
fix x range x 104.9 105.1
fix y
fix x y z range z -0.1 0.1
set grav 0 0 grav0       

set mech ratio 9.8e-6
solve step 30000
endcommand


if mech_ratio<1.0e-5
  k11=ks
                k12=k12

else
  k12=ks
                k11=k11
endif
ks=(k11+k12)/2

endloop
fosfile0=string(size1)+'_'+string(size2)+'_fos'+'.sav'

command
save fosfile0
endcommand
end
SSR

pr ks
基于一个滑坡，长150m，宽200m，厚度5米，坡度45度， 上部土层，下部为弱风化基岩。
想用强度折减法和静力平衡法来做。
                 滑坡岩、土体物理力学指标
土体类型        ρ(kg/m3)        K(GPa)        G(GPa)        C(kPa)        Φ(°)        υ        ES(GPa)
土体        1530        2.66        0.88        17        24.2        0.35        2.39
岩体        2720        43.9        30.2        36.5        22.7        0.22        73.8

ajianme

请高手给予指正，谢谢。最好是能把在下的编程改完善，加上简单的解释，具有可读性，对于提高将会很有用，谢谢。有心者可以加我QQ，交个朋友