求助！！！！谢谢！计算边坡安全系数！ [复制链接]

使用FLAC3d计算某一边坡的安全系数，但是再计算出初始应力状态以后，在计算安全系数时出现计算错误

https://bbs.yantuchina.com/read.php?tid=323878

程序

new

gen zone brick p0 0 0 0 p1 140 0 0 p2 0 20 0 p3 104 0 120 p4 140 20 0 p5 104 20 120 p6 140 0 120 p7 140 20 120 size 10 5 12

gen zone brick p0 -14 0 -20 p1 140 0 -20 p2 -14 20 -20 p3 -14 0 0 size 11 5 2

fix x y z range z -20

fix x range x 140

fix x range x -14

fix y

;初始条件设置

model e

prop density 2800 bulk 4.72e12 shear 3.54e12

set gravity 0 0 -10

solve

save solpe1-chu-1.sav

;位移、速度清零

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

;安全系数求解

model mohr

prop density 2800 bulk 4.72e10 shear 3.54e10 coh 4.5e7 tens 1e6 fric 35 dila 35

solve fos file slopepoumian1fos2201.sav associated

请各位大神帮忙提一下修改意见，谢谢