sos~~恳请各位大侠给小菜鸟的指点~~sos

ziguayida

我是一名刚学习FLAC 3d的菜鸟，最近在模拟一个库区。我模拟的是一个长宽为530X400的库区，四周围有8m高的大坝，我取1/4作为模拟的对象。库区内填的是污泥。我初始化应力时，是将污泥设置为空网格。模型是建好了，初始化应力后，重新赋值，却发现计算步只有10，不知道是怎么回事，恳请各位大侠能够给点指点！

ziguayida

这是我建模的命令
new
generate zone brick p0=(35 35 12) p1=(255 35 12) p2=(69.6 69.6 12) p3=(55 55 20) p4=(255 69.6 12) &
                    p5=(65 65 20) p6=(255 55 20) p7=(255 65 20) size 100 17 4 group daba;大坝x方向
pl set center auto
pl set rot 100 0 0
pl bl gr sh on
generate zone brick p0=(35 320 12) p1=(35 35 12) p2=(69.6 320 12) p3=(55 320 20) p4=(69.6 69.6 12) &
                    p5=(65 320 20) p6=(55 55 20) p7=(65 65 20) size 180 17 4 group daba;大坝y方向

generate zone brick p0=(255 69.6 12) p1=(255 320 12) p2=(69.6 69.6 12) p3=(255 65 20) p4=(69.6 320 12) &
                    p5=(65 65 20) p6=(255 320 20) p7=(65 320 20) size 180 100 4 group wuni;污泥

save shangbu.sav


gen zone brick  p0=(0 0 8) p1=(255 0 8) p2=(35 35 8) p3=(0 0 12) p4=(255,35,8) p5=(35,35,12) &
                p6=(255,0,12) p7=(255,35,12)  size 100 17 2 group badi1;地基一层

gen zone brick  p0=(35 35 8) p1=(35,320,8) p2=(0,0,8) p3=(35,35,12) p4=(0,320,8) p5=(0,0,12) &
                p6=(35,320,12) p7=(0,320,12) size 180 17 2 group badi1;地基一层

gen zone brick  p0=(69.6 69.6 8) p1=(69.6 320 8) p2=(35 35 8) p3=(69.6 69.6 12) p4=(35 320 8) &
                    p5=(35 35 12) p6=(69.6 320 12) p7=(35 320 12) size 180 17 2 group badi1;

gen zone brick  p0=(35 35 8) p1=(255 35 8) p2=(69.6 69.6 8) p3=(35 35 12) p4=(255 69.6 8) &
                    p5=(69.6 69.6 12) p6=(255 35 12) p7=(255 69.6 12) size 100 17 2 group badi1;

gen zone brick  p0=(69.6 69.6 8) p1=(255 69.6 8) p2=(69.6 320 8) p3=(69.6 69.6 12) size 100 180 2 group badi1;



save shangbu-yiceng



gen zone brick  p0=(0 0 0) p1=(255 0 0) p2=(35 35 0) p3=(0 0 8) p4=(255,35,0) p5=(35,35,8) &
                p6=(255,0,8) p7=(255,35,8)  size 100 17 4 group badi2;地基二层
gen zone brick  p0=(35 35 0) p1=(35,320,0) p2=(0,0,0) p3=(35,35,8) p4=(0,320,0) p5=(0,0,8) &
                p6=(35,320,8) p7=(0,320,8) size 180 17 4 group badi2;地基二层

gen zone brick  p0=(69.6 69.6 0) p1=(69.6 320 0) p2=(35 35 0) p3=(69.6 69.6 8) p4=(35 320 0) &
                    p5=(35 35 8) p6=(69.6 320 8) p7=(35 320 8) size 180 17 4 group badi2;地基二层

gen zone brick  p0=(35 35 0) p1=(255 35 0) p2=(69.6 69.6 0) p3=(35 35 8) p4=(255 69.6 0) &
                    p5=(69.6 69.6 8) p6=(255 35 8) p7=(255 69.6 8) size 100 17 4 group badi2;地基二层

gen zone brick  p0=(69.6 69.6 0) p1=(255 69.6 0) p2=(69.6 320 0) p3=(69.6 69.6 8) size 100 180 4 group badi2;地基二层




pl set center 1.293e2,1.465e2,-3.923e2
pl set rot 30,0,310
pl bl gr sh on


save shangbu-yiceng-erceng


range name=daba group daba
range name=wuni group wuni
range name=badi1 group badi1
range name=badi2 group badi2


model mohr


save moxingwangge.sav

ziguayida

;初始化应力
prop bulk=21.7e6 shear=10e6 cohesion=1e10 friction=29.81 tension=1e10 range daba
ini density=1551 range daba

model null range wuni

prop bulk=21.7e6 shear=10e6 cohesion=1e10 friction=29.81 tension=1e10 range badi1
ini density=1551 range badi1
prop bulk=7.8e6 shear=5.8e6 cohesion=1e10 friction=15 tension=1e10 range badi2
ini density=1800 range badi2
set grav 0 0 -9.81

；边界条件
fix x range x -0.1 0.1
fix x range x 254.9 255.1
fix y range y -0.1 0.1
fix y range y 319.9 320.1
fix x range z -0.1 0.1
fix y range z -0.1 0.1
fix z range z -0.1 0.1


history nstep=10
hist unbal
hist gp xdisp (55 160 20)
hist gp xdisp (65 160 20)
hist gp ydisp (165 55 20)
hist gp ydisp (165 65 20)
hist gp zdisp (165 60 20)
set mechanical force 50
solve
plot
history 1
history 6
print history


;等值线
plot
creat dengzhixian
add contour disp
add axes
show


;应力等值线
clear
add bcontour szz
add axes


;任意剖面上的等值线
creat poumian
set plane dip=90 dd=0 origin=165 55 20
add boundary behind
add bcontour szz plane
add axes
show
save dengzhixian.sav
plot current dengzhixian
plot show


save chushihuayingli





;赋实际值
prop cohesion=19.1e3 tension=19.1e3 range daba

model moh range wuni

prop bulk=4.7e6 shear=5.1e6 cohesion=2.7e3 friction=7.16 tension=2.7e3 range wuni
ini density=1195 range wuni

prop cohesion=19.1e3 tension=19.1e3 range badi1

prop cohesion=10e3 tension=10e3 range badi2


;初始化位移
initial xdisp=0 ydisp=0 zdisp=0
step 10
plot creat weiyi
plot copy poumian weiyi settings
plot add contour disp plane behind shade on
plot add axes
plot show

ziguayida

大侠们，我很急着要计算结果，小弟先感谢了您们了！！！！

ziguayida

大侠们，请帮帮忙吧！！