FLAC3D滑坡稳定性分析,静力平衡和强度折减法实例
那位前辈曾做过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, 20p0 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
请高手给予指正,谢谢。最好是能把在下的编程改完善,加上简单的解释,具有可读性,对于提高将会很有用,谢谢。有心者可以加我QQ,交个朋友
页:
[1]