在做数值模拟的时候遇到这样一个问题,我在调用命令的时候,命令的内容是一样的,但是得到的结果不一样。
命令流:
model mohr
prop bulk = 1.0E8 shear = 1.2E8 coh 8e6 tens 2.0e6 dens 2400.0 friction = 37
pl cont pp
fix x y z ran z 0
fix x ran x 0
fix x ran x 600
fix y
set grav 0 0 -10
config fluid
set fluid on
model fl_iso
prop por 0.5 perm 1e-8
ini fden 1000 ften 0.0
ini dens 2000
ini sat 0.0 ran z 0 130
ini fmod 2e3
ini pp 1.3e6 grad 0 0 -1e4 range x -0.1 0.1 z 0 130
fix pp ran x -0.1 0.1
app pp 0.2 ran x 599.9 600.1
set mech ratio 1e-4
solve
save 233.sav