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
save  dyn.sav
set log off

东黑子

我遇到过这个问题，是你的计算时间的问题，你试着调一下你的计算时间，你所输入的加速度波的持续时间可能不是30s吧

pisces15

我碰到同样的错误，是有一部分材料赋值赋掉了，检查一下材料吧。

oldfish2006

我也遇到了类似的问题，但是好像把计算时间调一少后就解决了

cuidazhen

我也出现了同样的问题
也是一个超大型模型的动力分析问题
错误一样是
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 sta  0
ini xvel 0   yvel 0  zvel 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文件）
前面进行了静力的地应力作用和开挖的计算
如果真的是材料赋值没有赋上的话，静力计算的时候为什么都都没有发现有任何问题呢？

还有请问楼上，把计算时间“调一少”是什么意思呢，怎么解决的呢

mihang

我也遇到了这样的问题。高手请指教一下呀。

mihang

我找到了。就静态边；界的时候不能用输入加速度。应该输入应力时程，把加速度积分成速度。再用公式转成应力。

cuidazhen

我找到了。就静态边；界的时候不能用输入加速度。应该输入应力时程，把加速度积分成速度。再用公式转成应力。
mihang 发表于 2009-8-22 14:51

你产生这个错误的原因应该和我及楼主的都不同吧
楼主用了自由边界，所以加的是加速度，我用了静态边界，加的是力。这一点上应该是没有错的啊。

mihang

LZ的底面明明是用的静态边界..四周用的才是自由场.加外载当然要加底面上了.

simweflac

我遇到的问题和你的一样
但我的只是静力计算
也是math error overflow
计算出来的位移和应力都超大，e141

计算的是一个井筒的某一段的应力和位移情况，分步开挖 位于深部软岩层 地应力较大，27MPa，静水压力场 井筒为圆形，模型仅为1/4圆， 模型沿重力方向(y轴正方向)10米厚，水平方向（x轴和z轴）则是80米*80米的正方形 最初的模型单元数约3.9万个，math errorverflow 后听取别人建议，减少单元数，到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

mihang

这个问题就没有高手能解决的吗？
我的程序简单。也是有这个问题。
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
谢谢高手门指点一二

mihang

能是网格划分有问题吗？还是材料参数赋值有问题？这个问题解决了的前辈们指点一下后生们。

simweflac

waiting for answer。。。

simweflac

我的math errorverflow问题解决了
是材料参数的问题

初始平衡和应力释放之后
给混凝土衬砌和卸压层赋材料参数的时候，也同时赋了密度的参数
但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

yinyunteng

我也遇到了同样的问题，请问楼主的问题解决了没有？

yinyunteng

我是用ansys建的模型，然后导入flac进行计算的，一开始静力计算都没有问题，但到了动力计算的时候就出现错误提示：math error : invalid
; @: U# i1 g, Y File:  \Solve.cpp
+ b; a5 K. F# v# b- t, ]Line 1276
一开始以为是由于开挖造成的，后来把模型直接建成开挖后的模型进行计算还是出现上述错误，另外开挖的那部分还出现了单元畸变。在下实在没有办法解决了，还望高手能够不吝赐教！

Euler

把此帖再顶起来。。

luoyecolin

我也遇到了这个问题，请高手指点一下。

luoyecolin

我的是在进行初始平衡时，  solve 后出现的。  请指教，谢谢！！

Davidzk89

围观 啊