求助动力加载计算中出现Math error :zero divide File:.\Solve.cpp Line:

MMCHUL

我要做的模型其实挺简单的，就是一个地铁双洞隧道的地面动力响应问题
网格是用anasys生成的，经过了平衡→开挖→衬砌道床参数→平衡→地面明沟开挖→平衡→设置动力参数→动力加载

部分网格：
问题出在动力计算 的时候，是不是道床网格质量不好？可是更好的我画不出来了。。。求各位大神和路过的亲们帮我看看问题究竟出在哪里，帮我改一改，老师最近就要成果实在是急。。。上传不了那么大的附件，所以下面是部分代码

;-----------------静力参数----------------

plot set rot 90 0 0
plot set mag 1.3
set gravity 0 -10 0
fix y range y -39.9 -40.1
fix x range x -50.1 -49.9
fix x range x 49.9 50.1
fix z range z -0.1 0.1
fix z range z -0.9 -1.1
model  mohr  
ini density 2200
prop b 1.923e9 s 4.167e9 fric 35 coh 0.5e11 tension 0.5e11
hist unbal
plot con syy ou on
plot add hist 1
solve
prop b 1.923e9 s 4.167e9 fric 35 coh 0.5e6 tension 0.5e6
solve
;-----------------开挖和平衡-----------------------
def derive
;E_mod=Es_mod*(1+p_ratio)(1-2*p_ratio)/(1-p_ratio)
s_mod=E_mod/(2.0*(1.0+p_ratio))
b_mod=E_mod/(3.0*(1.0-2.0*p_ratio))
end
;开挖隧道部分
model null range group 1
model null range group 3
;衬砌参数
set E_mod=34.5e9
set p_ratio=0.25
derive

ini density 2600 range group 5
ini density 2600 range group 6
model elas range group 5
model elas range group 6
prop bulk b_mod s s_mod range group 5
prop bulk b_mod s s_mod range group 6
;道床参数
set E_mod=31.5e9
set p_ratio=0.25
derive

ini density 2500 range group 2
ini density 2500 range group 4
model elas range group 2
model elas range group 4
prop bulk b_mod s s_mod range group 2
prop bulk b_mod s s_mod range group 4
solve
;---------------------------------加载-------------

;加载
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
call hist.txt
call table2.txt
config dyn
set dyn=on
set dyn damp local 0.125
set dynamic time = 0.0
set dyn dt=5.0e-5
def the_jiedian
   jiedian=sqrt(2.7^2-2.056^2)
end
the_jiedian
def get_jiedian
   x1=-(2+jiedian)-0.768
   x2=-(2+jiedian)+0.768
   x3=(2+jiedian)-0.768
   x4=(2+jiedian)+0.768
end
get_jiedian

sel node id 444441 x1 16.944 0
sel node id 444442 x2 16.944 0
sel node id 444443 x3 16.944 0
sel node id 444444 x4 16.944 0
sel link 444441 target zone
sel link 444442 target zone
sel link 444443 target zone
sel link 444444 target zone

call boundary.txt
save 加载前1.sav
apply yf 1.0 hist table 2  range id 444441
apply yf 1.0 hist table 2  range id 444442
hist nstep=100
hist dytime
solve age=2.0
save con1.sav


求各位大神帮我看看啊！

MMCHUL

搞不好是静态边界问题？如下：
free x y z range x -49.9 -50.1
free x y z range x 49.9 50.1
free x y z range z -0.1 0.1
free x y z range z -0.9 -1.1
free x y z range y -40.1 -39.9

apply nquiet squiet range x -49.9 -50.1
apply nquiet squiet range x  49.9 50.1
apply nquiet squiet range z -0.1 0.1
apply nquiet squiet range z -0.9 -1.1
apply nquiet squiet range y -40.1 -39.9

MMCHUL

自己顶顶 求看到啊！！

weijunzi006

MMCHUL 发表于 2013-4-5 22:13
自己顶顶 求看到啊！！

我也遇到了这个问题，真实捉鸡啊！！

xiwanshan

楼主这个问题解决了吗，我也遇到了，能解答一下吗，谢谢

gaoquan2009

版本问题，换一个3.0版本试一下，我在网上重新下了一版本试了就没有出现这个问题了。

xiwanshan

gaoquan2009 发表于 2013-10-31 02:58
版本问题，换一个3.0版本试一下，我在网上重新下了一版本试了就没有出现这个问题了。 ...

你好我也出现这种问题，但是找不到其他版本的flac3d，你用的什么版本的，能不能给我传一份，不胜感激！

漫步云端1990

这到底是为什么啊？有大神能解决一下吗？很困恼啊