找回密码
 注册
Simdroid-非首页
查看: 137|回复: 3

[命令/FISH] FLAC3D自带的例子(EXCAV.DAT)计算出错的疑惑

[复制链接]
发表于 2006-10-31 22:59:00 | 显示全部楼层 |阅读模式 来自 新加坡
初次使用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
group  wall    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 x  0.0 4.0 y  3.9 4.1 z  0.0 5.0
apply nstress 0 grad 0 0 -1.2e4 range x  3.9 4.1 y  0.0 4.0 z  0.0 5.0
apply nstress            -7.5e4 range x  0.0 4.0 y  0.0 4.0 z  4.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 x  0.0 4.0 y  3.9 4.1 z  0.0 5.0
apply nstress 0 grad 0 0 -1.e4 range x  3.9 4.1 y  0.0 4.0 z  0.0 5.0
apply nstress            -5.e4 range x  0.0 4.0 y  0.0 4.0 z  4.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 x  0.0 4.0 y  3.9 4.1 z  0.0 5.0
apply nstress 0 grad 0 0 -1.e4 hist relax &
            range x  3.9 4.1 y  0.0 4.0 z  0.0 5.0
apply nstress            -5.e4 hist relax &
            range x  0.0 4.0 y  0.0 4.0 z  4.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 编辑 ]
发表于 2006-11-1 00:10:16 | 显示全部楼层 来自 湖北武汉
Simdroid开发平台
有个问题是
FLAC3D 3.00-257
这个版本FLAC3D的pj有点问题
我遇到过

如果是直接使用软件自带示例文件,应该不是数据输入的问题。

供参考
发表于 2006-11-3 12:06:13 | 显示全部楼层 来自 四川成都
这个版本的是有问题。很多计算中都会出现错误,包括interface计算等。。。
不建议使用,后头貌似还有个261的是没问题的。
发表于 2011-8-23 11:01:37 | 显示全部楼层 来自 江苏南京
回复 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 x  0.0 4.0 y  3.9 4.1 z  0.0 5.0  
apply nstress 0 grad 0 0 -1.2e4 range x  3.9 4.1 y  0.0 4.0 z  0.0 5.0  
apply nstress       -7.5e4 range x  0.0 4.0 y  0.0 4.0 z  4.9 5.1  这几个力是为了平衡墙体而设置的吗?
(3)再一个是函数relaxsetup,它在整个程序运行过程中共执行了多少次啊?还是只执行了一次?在函数relax执行过程中,step0的值变吗?step的值是随着程序执行逐渐累加的吧?
很期待前辈的解答,谢谢!
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 02:28 , Processed in 0.040450 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表