桩-土水动力计算总是math error，程序如下求指点

saiyin

RT
模拟的是振动台试验，土层分两部分，中间有钢管桩，水位高1m。
计算到2万步左右   出现“Math errorverflow invalid File:.\Solve.cpp Line:1276”
一直都是自己研究，进展缓慢，上来求各位前辈指点，问题出在哪里？
程序如下：
new
title
试验测试一
config dynamic fluid
;
gen zone brick size 3 24 10 &
p0 (0,0,0) p1 (0.25,0,0) p2 (0,2.4,0) p3 (0,0,0.81) p4 (0.25,2.4,0) &
p5 (0,2.4,0.81) p6 (0.25,0,0.81) p7 (0.25,2.4,0.81)
gen zone brick size 21 24 10 &
p0 (0.25,0,0) p1 (2.35,0,0) p2 (0.25,2.4,0) p3 (0.25,0,0.81) p4 (2.35,2.4,0) &
p5 (0.25,2.4,0.81) p6 (2.35,0,1) p7 (2.35,2.4,1)
gen zone brick size 1 24 10 &
p0 (2.35,0,0) p1 (2.4,0,0) p2 (2.35,2.4,0) p3 (2.35,0,1) p4 (2.4,2.4,0) &
p5 (2.35,2.4,1) p6 (2.4,0,1) p7 (2.4,2.4,1)
group LooseSandLayer
gen zone brick size 21 24 6 &
p0 (0.25,0,1.81) p1 (2.35,0,2) p2 (0.25,2.4,1.81) p3 (0.25,0,2.11) p4 (2.35,2.4,2) &
p5 (0.25,2.4,2.11) p6 (2.35,0,2.3) p7 (2.35,2.4,2.3)
group FrozenSiltLayer range group LooseSandLayer not
range name=LooseSandLayer group LooseSandLayer
range name=FrozenSiltLayer group FrozenSiltLayer
;
model mohr range LooseSandLayer
model elastic range FrozenSiltLayer
interface 1 face range plane normal (-0.19,0,2.1) origin (0.25,0,1.81) dist 0.001 &
x (0.25,2.35) y (0,2.4)
ini z add -1 range z 1.5 2.5
prop bulk 1.5e7 shear 0.7e7 coh 1.7e4 fric 30 density 1520 range LooseSandLayer
prop bulk 1.435e8 shear 1.36e8 density 2000 range FrozenSiltLayer
interface 1 prop kn 4.9e9 ks 4.9e9 fric 20
model fl_iso range LooseSandLayer
prop poro 0.41 perm 1.73e-8
ini fmod 2e9 ftens -1e-3 fdens 1000
ini pp -10e3 grad 0 0 10e3 range z 0 1
model fl_null range FrozenSiltLayer
;
fix z range z (-0.001,0.001)
fix x range x (-0.001,0.001)
fix x range x (2.399,2.401)
fix y
;
set grav 0 0 -10
set fluid off dyn off
ini fmod 0
def ini_conf
_k0=0.4
pnt=zone_head
loop while pnt # null
val=_k0*z_szz(pnt)+(_k0-1.)*z_pp(pnt)
z_sxx(pnt)=val
z_syy(pnt)=val
pnt=z_next(pnt)
end_loop
end
ini_conf
solve;初始应力平衡
;
set dyn on fluid on
ini fmod 2e9
set fluid pcut on
model finn range LooseSandLayer
prop bulk 1.5e7 shear 0.7e7 coh 1.7e4 fric 30 range LooseSandLayer
ini density 1520
prop ff_latency 50 ff_switch 0 ff_c1 0.8 ff_c2 0.79 ff_c3 0.45 ff_c4 0.73 range LooseSandLayer
;
sel pile id=1 begin (1.2,1.2,0) end (1.2,1.2,0.9) nseg 18
sel pile id=1 begin (1.2,1.2,0.9) end (1.2,1.2,1.2) nseg 6
sel pile id=1 begin (1.2,1.2,1.2) end (1.2,1.2,1.4) nseg 4
sel pile id=1 prop dens 7850 emod 2e11 nu 0.3 xcarea 0.0002724  per 0.157 &
xcj 9.9e-9 xciy 4.95e-9 xciz 4.95e-9
sel pile id=1 prop cs_sk 1.3e11 cs_scoh 1.0e10 cs_sfric 30 &
cs_nk 1.3e9 cs_ncoh 1.0e4 cs_nfric 30 range z (0,0.9)
sel pile id=1 prop cs_sk 1.3e11 cs_scoh 1.0e10 cs_sfric 30 &
cs_nk 1.3e11 cs_ncoh 1.0e10 cs_nfric 30 range z (0.9,1.2)
;
sel node fix x y z xr yr zr range z (-0.001,0.001);端承桩，底部固定，对错未知
sel node apply force (0,0,-2500) range z (1.399,1.401);桩顶配重，不知这样表示合适不
;
table 1 read C:\Users\Administrator\Desktop\新建文件夹\第一个加载——正弦波f.txt
free x y z range z -0.001 0.001
apply nquiet dquiet squiet range z -0.001 0.001
apply dstress 1.0 hist table 1 range z -0.001 0.001
apply ff
group ff_corner
group ff_side range x 0 2.4
group ff_side range y 0 2.4
group main_grid range x 0 2.4 y 0 2.4
;
sel set damp combined
set dyn damp local 0.314
set dyn multi on
set small;若set large,则会出现错误zero volume tet in zone 4436
;hist id=1-24 xdisp pp moment xvel fx等
solve age 20

saiyin

没有人么。。求大神指点

saiyin

或者有人一起讨论也行。。共同提高嘛

高霞

你好，我无法加你好友，我qq号是850592360.我也做振动台，我是初学flac，不嫌弃的话加我好友吧，qq。也许咱们可以一起探讨

weijunzi006

高霞 发表于 2012-10-29 14:56
你好，我无法加你好友，我qq号是850592360.我也做振动台，我是初学flac，不嫌弃的话加我好友吧，qq。也许咱 ...

问题解决了么，我也遇到这样的问题

yhx_ustb

我也是这个问题，不知道是为什么，

yhx_ustb

我觉得可能是模型太小的原因，但是不确定