初始应力的生成？
n;***************************************************************
;导入模型
impgrid Wudongde.flac3d
;***************************************************************
plot bl gr
plot add axes red
;通过旋转坐标
def xiugai_11
p_gp=gp_head
loop while p_gp # null
w_temp1=gp_xpos(p_gp)
w_temp2=gp_ypos(p_gp)
w_temp3=gp_zpos(p_gp)
gp_xpos(p_gp)=w_temp1
gp_ypos(p_gp)=-w_temp3
gp_zpos(p_gp)=w_temp2
p_gp = gp_next(p_gp)
endloop xiugai_11
end
xiugai_11
ini x add 187.65
ini y add10
ini z add243.05
fix x y z ran z 0
fix x ran x 0
fix x ran x 689.96
fix y
model elas ran group Jiyan
model mohr ran group Fgc1
model mohr ran group Fgc2
model mohr ran group Fgc31
model mohr ran group Fgc32
model mohr ran group Kw1
model mohr ran group Kw2
model mohr ran group Kw3
model mohr ran group Kw4
model mohr ran group Fsq3
model null ran z 0 111 not
prop bulk 1.11e10 shear 0.833e10 ran group Jiyan
ini dens 2500 ran group Jiyan
prop bulk 3e8 shear 1e8 coh 30e3 fric 40 ran group fgc1 any group kw4 a gr fsq3 a
ini dens 2230 ran group fgc1 any group kw4 a gr fsq3 a
prop bulk 3e8 shear 1e8 coh 30e3 fric 40 ran group fgc2 any group kw3 a
ini dens 2340 ran group fgc2 any group kw3 a
prop bulk 3e8 shear 1e8 coh 30e3 fric 40 ran group fgc31 any group kw2 a
ini dens 2400 ran group fgc31 any group kw2 a
prop bulk 3e8 shear 1e8 coh 30e3 fric 40 ran group fgc32 any group kw1 a
ini dens 2400 ran group fgc32 any group kw1 a
set grav 0 0 -9.8
hist id=1 unbal
solve
save WDD1.sav
这是我算一个土石围堰的初始应力的命令流，为什么停不下来呢？？大家帮忙看看哈！ 附件的那个东西是命令流及附件 谢谢大家了 没人知道吗？？？？
页:
[1]