FLAC3D自带的例子(EXCAV.DAT)计算出错的疑惑
初次使用FLAC3D，在跑其自带的例子(EXCAV.DAT)的时候，出现以下错误：***Math Error: invalid number
不知是何原因，在跑其他自带例子的时候，也会出现类似的错误。
我用的版本是：FLAC3D 3.00-257
;------------------------------------------------------
; Excavation in a saturated soil (EXCAV.DAT)
;------------------------------------------------------
config fluid
; --- geometrical model ---
gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1
group soil
group excavate range x 0 4 y 0 4 z 0 5
group wal1 range x 4 5 y 0 5 z 0 7
group wal2 range x 0 4 y 4 5 z 0 7
groupwall range group wal1 any group wal2 any
; --- fluid flow model ---
model fl_iso
prop perm 1e-12 poro 0.3
ini fdensity 1e3
ini fmod 2.0e9 ftens -1e-3
model fl_null range group excavate
model fl_null range group wall
ini pp 0 grad 0 0 1e4
fix pp range z -0.1 0.1
fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1
; --- mechanical model ---
model elas
prop bul 3.9e6 shea 2.8e6
model null range group excavate
ini density 1.2e3
ini density 1.5e3 range group wall
fix x range x -.1 .1
fix x range x 11.9 12.1
fix y range y -.1 .1
fix y range y 11.9 12.1
fix z range z 11.9 12.1
; initial total stresses
ini szz 0 grad 0 0 -1.5e4
ini sxx 0 grad 0 0 -1.2e4
ini syy 0 grad 0 0 -1.2e4
apply nstress 0 grad 0 0 -1.2e4 range x0.0 4.0 y3.9 4.1 z0.0 5.0
apply nstress 0 grad 0 0 -1.2e4 range x3.9 4.1 y0.0 4.0 z0.0 5.0
apply nstress -7.5e4 range x0.0 4.0 y0.0 4.0 z4.9 5.1
; --- setting ---
set gravity 0 0 10
; --- initial state ---
solve force 1;check initial equilibrium
; --- histories ---
set hist_rep 40
hist fltime
hist gp pp 0 0 6
hist gp xdis 4 0 0
hist gp xdis 4 0 2
hist gp xdis 4 2 0
hist gp zdis 0 0 5
hist gp zdis 2 0 5
hist gp zdis 4 0 5
hist gp zdis 2 2 5
hist gp zdis 4 2 5
hist gp zdis 4 4 5
hist gp zdis 10 0 1
hist gp zdis 10 0 2
;
; --- excavation ---
set fluid off
;apply pore pressure at walls
apply nstress 0 grad 0 0 -1.e4 range x0.0 4.0 y3.9 4.1 z0.0 5.0
apply nstress 0 grad 0 0 -1.e4 range x3.9 4.1 y0.0 4.0 z0.0 5.0
apply nstress -5.e4 range x0.0 4.0 y0.0 4.0 z4.9 5.1
solve ;force 1
save exc1.sav
;
; --- drainage ---
apply remove nstress
def relaxsetup
step0 = step
end
relaxsetup
def relax
rstep = step - step0
if rstep < ncyc then
relax=1.0-(float(rstep)/float(ncyc))
else
relax = 0.0
endif
end
set ncyc = 1000
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x0.0 4.0 y3.9 4.1 z0.0 5.0
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x3.9 4.1 y0.0 4.0 z0.0 5.0
apply nstress -5.e4 hist relax &
range x0.0 4.0 y0.0 4.0 z4.9 5.1
cyc ncyc
solve
save exc2.sav
;
; --- percolation ---
set fluid on
fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1
cyc 9000
save exc3.sav
plot create excav
plot set rot 200 0 195
plot set cent 6 6 6
plot set dist 39.18
plot set magn 0.8
plot set plane ori 0 0 0 normal 0 -1 0
plot add cont pp out on
plot add block group lgra range group wall
plot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7
plot add flow plane
plot show
ret
[ 本帖最后由 benjackxu 于 2006-12-31 00:34 编辑 ] 有个问题是
FLAC3D 3.00-257
这个版本FLAC3D的pj有点问题
我遇到过
如果是直接使用软件自带示例文件，应该不是数据输入的问题。
供参考 这个版本的是有问题。很多计算中都会出现错误，包括interface计算等。。。
不建议使用，后头貌似还有个261的是没问题的。 回复 1# mygeotech
前辈对于这个例子我有一些不明白的地方，希望你能指导一下？
（1）开挖完成后有这样几句命令：
; initial total stresses
ini szz 0 grad 0 0 -1.5e4
ini sxx 0 grad 0 0 -1.2e4
ini syy 0 grad 0 0 -1.2e4，我的理解是，这是对整个土体进行应力设置但为什么sxx和syy取成1.2e4呢？按道理水平方向应该是竖向的一半吧（根据公式v/1-v）。
（2）接下来又设置了如下命令：
apply nstress 0 grad 0 0 -1.2e4 range x0.0 4.0 y3.9 4.1 z0.0 5.0
apply nstress 0 grad 0 0 -1.2e4 range x3.9 4.1 y0.0 4.0 z0.0 5.0
apply nstress -7.5e4 range x0.0 4.0 y0.0 4.0 z4.9 5.1这几个力是为了平衡墙体而设置的吗？
（3）再一个是函数relaxsetup，它在整个程序运行过程中共执行了多少次啊？还是只执行了一次？在函数relax执行过程中，step0的值变吗？step的值是随着程序执行逐渐累加的吧？
很期待前辈的解答，谢谢！
页:
[1]