自己写的边坡稳定性分析的命令，求高人指点下。
newgen zone brick size 20 20 10 p0 0 0 0 p1 30 0 0 &
p2 0 40 0 p3 0 0 10 p4 30 40 0 p5 0 40 10 &
p6 30 0 5 p7 30 40 5 group ceng1
gen zone brick size 20 20 10 p0 0 0 10 p1 30 0 5 &
p2 0 40 10 p3 0 0 20 p4 30 40 5 p5 0 40 20 &
p6 24.5 0 14 p7 24.5 40 14 group ceng2
gen zone brick size 20 20 5 p0 0 0 20 p1 24.5 0 14 &
p2 0 40 20 p3 0 0 27 p4 24.5 40 14 p5 0 40 27 &
p6 14 0 25.3 p7 14 40 25.3 group ceng3
gen zone brick size 20 20 3 p0 0 0 27 p1 14 0 25.3 &
p2 0 40 27 p3 0 0 30 p4 14 40 25.3 p5 0 40 30 &
p6 8 0 29.5 p7 8 40 29.5 group ceng4
plot block group
pause
model mohr range group ceng1
model mohr range group ceng2
model mohr range group ceng3
model mohr range group ceng4
prop bulk 3.51e6shear 1.17e6coh 31000 fric 25.5range group ceng1
prop bulk 3.06e6shear 0.94e6coh 28000 fric 23.0range group ceng2
prop bulk 2.64e6shear 0.88e6coh 22000 fric 17.0range group ceng3
prop bulk 3.76e6shear 0.8e6coh 19500 fric 16.8range group ceng4
fix x range x -0.1 0.1
fix x range x 29.9 30.1
fix y range y -0.1 0.1
fix y range y 39.9 40.1
fix z range z -0.1 0.1
fix z range z 29.9 30.1
set gravity 0 0 -10
initial density=1870 range group ceng1
initial density=1850 range group ceng2
initial density=1860 range group ceng3
initial density=1950 range group ceng4
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
plot add gpfix scal 0.01
pause
plot cont pp
plot add sk
pause
plot con szz
plot add sk
pause
hist unbal
hist gp xdisp 8 0 29.5
hist gp zdisp 8 0 29.5
hist gp xdisp 30 0 5
hist gp ydisp 30 0 5
hist gp zdisp 30 0 5
hist gp xdisp 30 40 5
hist gp ydisp 30 40 5
hist gp zdisp 30 40 5
pause
solve fos
请高手给指点下，为什么我运行了后，半天还没有出来结果啊。是不是哪里出现问题了？？谢谢。 强度折减法本来就很慢，请耐心等待 2# zlxiaoxi3 算了一下午，终于出来了。fos为0.88.。 强度折减法求解安全系数很慢，慢慢等，别着急....
页:
[1]