math error : invalid 是什么原因呀？
我做的大模型在动力计算时出错：math error : invalid
File: \Solve.cpp
Line 1276
是什么原因呀？
我的命令流：
set log on
set logfile dynhis.log
config dyn
set dyn on
table 100 read ns.txt
table 200 read ew.txt
table 300 read up.txt
free x y z range z -120.1 -119.9 ;去掉模型底部原有的静力条件
apply nquiet squiet dquiet ran z -120.1 -119.9
apply xaccel 1.0 his table 100 range z -120.1 -119.9
apply yaccel 1.0 his table 200 range z -120.1 -119.9
apply zaccel 1.0 his table 300 range z -120.1 -119.9
apply ff
group ff
set dyn damp local 0.05
his reset
his gp zdis 185,3.1,2.5
his gp xvel 0,0,-120
his gp yvel 0,0,-120
his gp zvel 0,0,-120
his gp xvel 0,0,2.5
his gp yvel 0,0,2.5
his gp zvel 0,0,2.5
his unbal
his dytime
set dyn multi on
pl ad his 1 vs 10
solve age 30
savedyn.sav
set log off 我遇到过这个问题，是你的计算时间的问题，你试着调一下你的计算时间，你所输入的加速度波的持续时间可能不是30s吧 我碰到同样的错误，是有一部分材料赋值赋掉了，检查一下材料吧。 我也遇到了类似的问题，但是好像把计算时间调一少后就解决了 我也出现了同样的问题
也是一个超大型模型的动力分析问题
错误一样是
math error :
overflow
invalid
File:\Solve.cpp
Line 1276
命令流如下：
new
set log on
rest ex2.sav
;读入3向地震波
table 1 read 50-5-1.str
table 2 read 50-5-2.str
table 3 read 50-5-3.str
config dyn
ini sta0
ini xvel 0 yvel 0zvel 0
ini xdisp 0 ydisp 0 zdisp 0
free x y z ;range z 499.9 500.1
app nquiet squiet dquiet rang x -800.1 -799.9
app nquiet squiet dquiet rang x 399.9 400.1
app nquiet squiet dquiet rang y -700.1 -699.9
app nquiet squiet dquiet rang y 299.9 300.1
app nquiet squiet dquiet rang z 499.9 500.1
app sxz 0.01 hist table 1 syz 0.01 hist table 2 szz 0.01 hist table 3 range z 499.9 500.1
set dyn time=0
set dyn multi on
set dyn damp local 0.0628
solve age 40.0
save dyn.sav
set log off
因为这个计算是在前面的静力计算的基础上进行的（见引用了ex2.sav文件）
前面进行了静力的地应力作用和开挖的计算
如果真的是材料赋值没有赋上的话，静力计算的时候为什么都都没有发现有任何问题呢？
还有请问楼上，把计算时间“调一少”是什么意思呢，怎么解决的呢 我也遇到了这样的问题。高手请指教一下呀。 我找到了。就静态边；界的时候不能用输入加速度。应该输入应力时程，把加速度积分成速度。再用公式转成应力。 本帖最后由 cuidazhen 于 2009-8-22 16:42 编辑
我找到了。就静态边；界的时候不能用输入加速度。应该输入应力时程，把加速度积分成速度。再用公式转成应力。
mihang 发表于 2009-8-22 14:51 http://forum.simwe.com/images/common/back.gif
你产生这个错误的原因应该和我及楼主的都不同吧
楼主用了自由边界，所以加的是加速度，我用了静态边界，加的是力。这一点上应该是没有错的啊。 LZ的底面明明是用的静态边界..四周用的才是自由场.加外载当然要加底面上了. 本帖最后由 simweflac 于 2009-9-14 20:18 编辑
我遇到的问题和你的一样
但我的只是静力计算
也是math error overflow
计算出来的位移和应力都超大，e141
计算的是一个井筒的某一段的应力和位移情况，分步开挖
位于深部软岩层
地应力较大，27MPa，静水压力场
井筒为圆形，模型仅为1/4圆，
模型沿重力方向(y轴正方向)10米厚，水平方向（x轴和z轴）则是80米*80米的正方形
最初的模型单元数约3.9万个，math error:overflow
后听取别人建议，减少单元数，到1.6万个左右
同时地应力减小到10MPa，依旧静水压力
可以计算完，但和预计的模型不同（网格比较粗糙，且地应力远小于原数值）
又在此基础上修改模型，加大地应力
却无论如何不能计算
现将命令流贴在这里
恳请大家指导:)不胜感激
因为是分步开挖，所以分组比较复杂
麻烦大家了，如有更好分组方法，也请大家指教:)
n
;井筒开挖前的岩心，5组
gen zon cyl p0 0 0 0 p1 2.74 0 0 p2 0 2 0 p3 0 0 2.74 size 5 3 20 gro c1
gen zon cyl p0 0 2 0 p1 2.74 2 0 p2 0 4 0 p3 0 2 2.74 size 5 3 20 gro c2
gen zon cyl p0 0 4 0 p1 2.74 4 0 p2 0 6 0 p3 0 4 2.74 size 5 3 20 gro c3
gen zon cyl p0 0 6 0 p1 2.74 6 0 p2 0 8 0 p3 0 6 2.74 size 5 3 20 gro c4
gen zon cyl p0 0 8 0 p1 2.74 8 0 p2 0 10 0 p3 0 8 2.74 size 5 3 20 gro c5
; 混凝土衬砌，5组
gen zon cshell p0 0 0 0 p1 3.49 0 0 p2 0 2 0 p3 0 0 3.49 dim 2.74 2.74 2.74 2.74 size 2 3 20 gro l1
gen zon cshell p0 0 2 0 p1 3.49 2 0 p2 0 4 0 p3 0 2 3.49 dim 2.74 2.74 2.74 2.74 size 2 3 20 gro l2
gen zon cshell p0 0 4 0 p1 3.49 4 0 p2 0 6 0 p3 0 4 3.49 dim 2.74 2.74 2.74 2.74 size 2 3 20 gro l3
gen zon cshell p0 0 6 0 p1 3.49 6 0 p2 0 8 0 p3 0 6 3.49 dim 2.74 2.74 2.74 2.74 size 2 3 20 gro l4
gen zon cshell p0 0 8 0 p1 3.49 8 0 p2 0 10 0 p3 0 8 3.49 dim 2.74 2.74 2.74 2.74 size 2 3 20 gro l5
; PU, 卸压层，类似于泡沫材料，非岩石，5组
gen zon cshell p0 0 0 0 p1 3.94 0 0 p2 0 2 0 p3 0 0 3.94 dim 3.49 3.49 3.49 3.49 size 1 3 20 gro u1
gen zon cshell p0 0 2 0 p1 3.94 2 0 p2 0 4 0 p3 0 2 3.94 dim 3.49 3.49 3.49 3.49 size 1 3 20 gro u2
gen zon cshell p0 0 4 0 p1 3.94 4 0 p2 0 6 0 p3 0 4 3.94 dim 3.49 3.49 3.49 3.49 size 1 3 20 gro u3
gen zon cshell p0 0 6 0 p1 3.94 6 0 p2 0 8 0 p3 0 6 3.94 dim 3.49 3.49 3.49 3.49 size 1 3 20 gro u4
gen zon cshell p0 0 8 0 p1 3.94 8 0 p2 0 10 0 p3 0 8 3.94 dim 3.49 3.49 3.49 3.49 size 1 3 20 gro u5
;风化围岩1
gen zon cshell p0 0 0 0 p1 19.94 0 0 p2 0 10 0 p3 0 0 19.94 dim 3.94 3.94 3.94 3.94 size 27 15 20 gro softmarl
;更远处的围岩2
gen zon cshell p0 0 0 0 p1 35 0 0 p2 0 10 0 p3 0 0 35 dim 19.94 19.94 19.94 19.94 size 10 15 20 ratio 1.16 1 1
gen zon radc p0 0 0 0 p1 60 0 0 p2 0 10 0 p3 0 0 60 p6 60 0 60 dim 35 35 35 35 size 0 15 20 6 ratio 1 1 1 1.15
gro marl range gro c1 not gro c2 not gro c3 not gro c4 not gro c5 not gro l1 not gro l2 not gro l3 not gro l4 not gro l5 not gro u1 not gro u2 not gro u3 not gro u4 not gro u5 not gro softmarl not
;初始平衡材料赋值，弹性
mo elas
pro bulk 0.5333e9 she 0.177778e9
; 边界条件，应力边界
; sxx=syy=szz=1100 m * 9.81 N/kg * 2500 kg/m3 = 27 MPa
; 10m*9.81*2300 =0.23 MPa
fix z range z -.1 .1
fix x range x -.1 .1
fix y range y 9.9 10.1
app szz -27e6 grad 0 -0.023e6 0 range z 59.9 60.1 y 0 10
app sxx -27e6 grad 0 -0.023e6 0 range x 59.9 60.1 y 0 10
app syy -27e6 range y -.1 .1
; 初始应力赋值，加快初始平衡计算
ini sxx=-27e6
ini syy=-27e6
ini szz=-27e6
set grav 0 9.81 0
ini density 2300.0
; 初始平衡计算
solve
sav liner0-0.sav
;围岩摩尔库仑模型，分步开挖第一组，井筒开挖前岩心弹模减半，应力释放模拟
ini xd=0 yd=0 zd=0 xv=0 yv=0 zv=0
m mohr
prop b=0.53333e9 s=0.177778e9 c=1.2e6 fric=25
prop b=0.26667E9 s=0.088889E9 c=0.4e6 fric=15 range gro c1
solve
; 第一步开挖，y方向2米深
m null range gro c1
ini xd=0 yd=0 zd=0 xv=0 yv=0 zv=0
; 混凝土衬砌，卸压层，风化围岩参数
prop b=0.104167E9 s=0.08474E9 c=0.4e6 fric=15 range gro softmarl
prop d=55 b=1666666.6 s=2.5E6 range gro u1
prop d=2400 b=1.50462E10 s=11.2847E9 c=3.34e6 fric=50 tension=2.8e6 range gro l1
hist gp zdisp 0 0 3.94
hist gp xdisp 3.94 0 0
hist gp zdisp 0 0 2.74
hist gp xdisp 2.74 0 0
hist gp ydisp 0 0 3.94
hist gp ydisp 0 0 2.74
solve
sav liner0-1.sav
我也出现了同样的问题
也是一个超大型模型的动力分析问题
错误一样是
math error :
overflow
invalid
File:\Solve.cpp
Line 1276
命令流如下：
new
set log on
rest ex2.sav
;读入3向地震波
t ...
cuidazhen 发表于 2009-8-21 15:33 http://forum.simwe.com/images/common/back.gif 这个问题就没有高手能解决的吗？
我的程序简单。也是有这个问题。
rest 11.sav
set dyn on
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
ini state 0
free x y z range z -0.1 0.1
table 1 read xacc.dat
apply xacc 1 hist table 1 yacc 0.0 zacc 0.0 range z -0.1 0.1
apply ff
set dyn multi on
solve age 1
谢谢高手门指点一二 能是网格划分有问题吗？还是材料参数赋值有问题？这个问题解决了的前辈们指点一下后生们。 waiting for answer。。。 我的math error:overflow问题解决了
是材料参数的问题
初始平衡和应力释放之后
给混凝土衬砌和卸压层赋材料参数的时候，也同时赋了密度的参数
但plot发现，此时密度赋值无效
后又单独用了ini density来赋值密度
就可以计算了
原来的名令：
prop d=2400 b=1666666.6 s=2.5E6 range gro u1
prop d=2400 b=1.50462E10 s=11.2847E9 c=3.34e6 fric=50 tension=2.8e6 range gro l1
改后的命令：
prop b=1666666.6 s=2.5E6 range gro u1
prop b=1.50462E10 s=11.2847E9 c=3.34e6 fric=50 tension=2.8e6 range gro l1
ini density 55 range gro u1
ini density 2400 range gro l1 我也遇到了同样的问题，请问楼主的问题解决了没有？ 我是用ansys建的模型，然后导入flac进行计算的，一开始静力计算都没有问题，但到了动力计算的时候就出现错误提示：math error : invalid
; @: U# i1 g, Y File:\Solve.cpp
+ b; a5 K. F# v# b- t, ]Line 1276
一开始以为是由于开挖造成的，后来把模型直接建成开挖后的模型进行计算还是出现上述错误，另外开挖的那部分还出现了单元畸变。在下实在没有办法解决了，还望高手能够不吝赐教！ 把此帖再顶起来。。 我也遇到了这个问题，请高手指点一下。 我的是在进行初始平衡时，solve 后出现的。请指教，谢谢！！ 围观 啊
页:
[1]
2