n
rest 111.sav
;====================================
;设置边界条件
fix x y z range z 542.3 542.5
fix x range x 1813.7 1813.9
fix x range x 608.8 609.0
fix y
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e10 shear 1e10
set gravity 0 0 -10.0
hist unb
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;===================================
;求解
model mohr
prop density 2200 bulk 8.5e8 shear 2.6e8 coh 1.5e6 tens 0.5e6 fric 35 ran group 1
prop density 2100 bulk 5.8e8 shear 1.9e8 coh 1.e6 tens 0.3e6 fric 27 ran group 3
prop density 1900 bulk 1.5e8 shear 0.5e8 coh 0.2e5 tens 0.3e5 fric 20 ran group 5
hist unb
hist gp xdisp 789.7 0 618
solve
采用如上计算边坡稳定状态，一直不收敛，不过选取参数就是这么大，看到书上有些例子参数还没用这么大都收敛的。检查模型没有错误的，因为模型本身也很小。。问过有些说边坡计算岩土体参数都偏小收敛都少，不太清楚，不过我计算几万步后ratio还是在 e-3左右，取消后查看位移云图这些根本不对。还望大神们多指教
