n
config dyn
set dyn off
;===================================
;建立网格模型
gen zone brick &
p0 0 0 0 p1 2 0 0 p2 0 2 0 p3 0 0 3 size 3 10 3
gen zone brick &
p0 2 0 0 p1 20 0 0 p2 2 2 0 p3 2 0 3 &
size 17 10 3 ratio 1.03 1 1
gen zone brick &
p0 2 0 3 p1 20 0 3 p2 2 2 3 p3 12 0 13 &
p4 20 2 3 p5 12 2 13 p6 20 0 13 &
p7 20 2 13 size 17 10 17 ratio 1.03 1 1
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 19.9 20.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 9.9 10.1
;====================================
;初始地应力的生成
model elas
prop density 2000 bulk 3e12 shear 1e12
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
;===================================
;安全系数求解
model mohr
prop density 2000 bulk 1.0e8 shear 3.0e7 &
coh 12380.0 tens 1.0e6 fric 20 dila 20
hist unbal
solve
save 1.sav
这是对边坡进行静力分析的步骤，可是我的程序算了250多万步了还没有算完，不知道哪里有问题？求教
