15-1边坡稳定性分析

gen zone radcyl &

p0 0 0 25 p1 add 80 0 0 p2 add 0 0 -25 p3 add 0 80 0 &

dim 24.5 24.5 12 12 ratio 1 1 1 1.1

gen zone radcyl &

p0 0 0 0 p1 add 80 0 0 p2 add 0 0 -15 p3 add 0 80 0 &

dim 12 12 12 12 ratio 1 1 1 1.1 fill size 5 7 10 10

gen zone b &

p0 0 -40 -15 p1 add 12 0 0 p2 add 0 40 0 p3 add 0 0 15 &

size 5 12 7

gen zone b &

p0 12 -40 -15 p1 add 68 0 0 p2 add 0 40 0 p3 add 0 0 15 &

size 10 12 7 ratio 1.1 1 1

gen zone b &

p0 12 -40 0 p1 add 68 0 0 p2 add 0 40 0 p3 add 12.5 0 25 &

p4 add 68 40 25 p5 add 12.5 40 25 p6 add 68 0 25 p7 &

add 68 40 25 size 10 12 10 rat 1.1 1 1

ini pp 0 grad 0 0 -1e4 range z -15.1 0

call water.fis

;设置mohr coulomb 模型及材料参数

mod mo

pro bulk 2e8 she 1e8 frick 45 coh 1e5 ten 1e5

;边界条件

fix x range x -.1 1

fix x range x 79.9 80.1

fix y range y -40.1 -39.9

fix y range y 79.9 80.1

fix x range z -15.1 -14.9

fix y range z -15.1 -14.9

fix z range z -15.1 -14.9

;初始化密度：地下水位之上为2500，之下为2600

def ini_dens

pnt = zone_head

loop while pnt # null

if z_pp(pnt)#0.0 then

z_denisty(pnt) = 2600

else

z_denisty(pnt) = 2500

endif

pnt = z_next(pnt)

endloop

end

ini_dens

set g 0 0 -10

;初始化重应力

ini syy -6.25e5 grad 0 0 2.5e4

ini sxx -6.25e5 grad 0 0 2.5e4

ini szz -6.25e5 grad 0 0 2.5e4

;监测

hist unbal

hist gp xdisp 24.5 -40 25

hist gp zdisp 24.5 -40 25

hist gp xdisp 24.5 0 25

hist gp ydisp 24.5 0 25

hist gp zdisp 24.5 0 25

hist gp ydisp 0 24.5 25

hist gp zdisp 0 24.5 25

;自动安全因子结算

solve fos

ret