Math error：
overflow
invalid
File:./solve.cpp
Line:1276
怎么回事 谢谢
下面是我的代码
rest 加桩后初始应力.sav
ini state 0
set dyn on mech on fluid on
set dyn mult on
set fluid time 0
set mech rat 1e-100
set fluid pcut on ;土体不产生负孔压
model finn ran group soilb
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7 f 30 c 0 perm 1e-8 ran group soilb
prop ff_latency=50 ran group soilb
prop ff_switch = 0 ran group soilb
prop ff_c1=0.8 ff_c2=0.79 ran group soilb
prop ff_c3=0.45 ff_c4=0.73 ran group soilb
ini fmod 2e8 ran gro soila any group soilb any group soilc any group toushui any
fix pp range group toushui
ini damp local 0.1571 range group toushui
ini damp local 0.3142 range group soila a group soilc a
ini damp local 0.4713 range group soilb
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
free x range z -10.1 -9.9
;底部速度施加
def setup
freq=5.0
ampl=2
omega = 2.0 * pi * freq
end
setup
def sine_wave
vv = 9.36e-2 * sin(omega*dytime)
if dytime < 2.0
sine_wave = dytime / 2.0 * vv
else
if dytime < 20.0
sine_wave = vv
else
if dytime <= 30.
sine_wave = (30.0 - dytime) / 10.0 * vv
endif
endif
endif
if dytime > 30.0
sine_wave = 0.0
endif
end
apply xvel = 1.0 hist sine_wave ran z -9.9 -10.1
apply ff
call ppr.dat
;历史变量监测
his reset
hist id 1 unbal
hist id 2 dytime
hist id 3 sine_wave
;加速度(底部输入加速度时程)
hist id 111 gp xacc 0 -2 -10
hist id 112 gp xacc 0.2 -2 -10
hist id 113 gp xacc 0.8 -2 -10
hist id 114 gp xacc 1 -2 -10
hist id 115 gp xacc 2 -2 -10
hist id 116 gp xacc 2.5 -2 -10
hist id 117 gp xacc 3 -2 -10
hist id 118 gp xacc 4.5 -2 -10
hist id 119 gp xacc -0.2 -2 -10
hist id 120 gp xacc -0.8 -2 -10
hist id 121 gp xacc -1 -2 -10
hist id 122 gp xacc -2 -2 -10
hist id 123 gp xacc -2.5 -2 -10
hist id 124 gp xacc -3 -2 -10
hist id 125 gp xacc -3.5 -2 -10
hist id 126 gp xacc -4 -2 -10
hist id 127 gp xacc -5 -2 -10
hist id 128 gp xacc -6.5 -2 -10
def get_ppr_zone
whilestepping
p_z43 = z_near(-1 , -2, 0)
ppr_431 = z_extra(p_z43,2)
ppr_432 = z_extra(p_z43,3)
ppr_433 = z_extra(p_z43,6)
ppr_434 = z_extra(p_z43,7)
p_z44 = z_near(-1, -2, -2.5)
ppr_441 = z_extra(p_z44,2)
ppr_442 = z_extra(p_z44,3)
ppr_443 = z_extra(p_z44,6)
ppr_444 = z_extra(p_z44,7)
p_z45 = z_near(-1, -2, -5)
ppr_451 = z_extra(p_z45,2)
ppr_452 = z_extra(p_z45,3)
ppr_453 = z_extra(p_z45,6)
ppr_454 = z_extra(p_z45,7)
p_z46 = z_near(-1, -2, -7.5)
ppr_461 = z_extra(p_z46,2)
ppr_462 = z_extra(p_z46,3)
ppr_463 = z_extra(p_z46,6)
ppr_464 = z_extra(p_z46,7)
p_z47 = z_near(-1, -2, -10)
ppr_471 = z_extra(p_z47,2)
ppr_472 = z_extra(p_z47,3)
ppr_473 = z_extra(p_z47,6)
ppr_474 = z_extra(p_z47,7)
p_z48 = z_near(-1, -3, 0)
ppr_481 = z_extra(p_z48,2)
ppr_482 = z_extra(p_z48,3)
ppr_483 = z_extra(p_z48,6)
ppr_484 = z_extra(p_z48,7)
p_z49 = z_near(-1, -3, -2.5)
ppr_491 = z_extra(p_z49,2)
ppr_492 = z_extra(p_z49,3)
ppr_493 = z_extra(p_z49,6)
ppr_494 = z_extra(p_z49,7)
p_z50 = z_near(-1, -3, -5)
ppr_501 = z_extra(p_z50,2)
ppr_502 = z_extra(p_z50,3)
ppr_503 = z_extra(p_z50,6)
ppr_504 = z_extra(p_z50,7)
p_z51 = z_near(-1, -3, -7.5)
ppr_511 = z_extra(p_z51,2)
ppr_512 = z_extra(p_z51,3)
ppr_513 = z_extra(p_z51,6)
ppr_514 = z_extra(p_z51,7)
p_z52 = z_near(-1, -3, -10)
ppr_521 = z_extra(p_z52,2)
ppr_522 = z_extra(p_z52,3)
ppr_523 = z_extra(p_z52,6)
ppr_524 = z_extra(p_z52,7)
p_z53 = z_near(-2, -3, 0)
ppr_531 = z_extra(p_z53,2)
ppr_532 = z_extra(p_z53,3)
ppr_533 = z_extra(p_z53,6)
ppr_534 = z_extra(p_z53,7)
p_z54 = z_near(-2, -3, -2.5)
ppr_541 = z_extra(p_z54,2)
ppr_542 = z_extra(p_z54,3)
ppr_543 = z_extra(p_z54,6)
ppr_544 = z_extra(p_z54,7)
p_z55 = z_near(-2, -3, -5)
ppr_551 = z_extra(p_z55,2)
ppr_552 = z_extra(p_z55,3)
ppr_553 = z_extra(p_z55,6)
ppr_554 = z_extra(p_z55,7)
p_z56 = z_near(-2, -3, -7.5)
ppr_561 = z_extra(p_z56,2)
ppr_562 = z_extra(p_z56,3)
ppr_563 = z_extra(p_z56,6)
ppr_564 = z_extra(p_z56,7)
p_z57 = z_near(-2, -3, -10)
ppr_571 = z_extra(p_z57,2)
ppr_572 = z_extra(p_z57,3)
ppr_573 = z_extra(p_z57,6)
ppr_574 = z_extra(p_z57,7)
end
set hist_rep 300
set dyn dt 3e-4
def solve_ages
loop n (1,20)
time_n = n
save_file = 'Dyn-' + string(n) + 'x1s.sav'
command
sol age time_n
save save_file
endcommand
endloop
end
solve_ages
看帖子有人说是材料参数 本构模型 的问题 或者是网格的问题 我看了 都没有什么问题啊
我最后把set dyn dt 3e-4改成 set dyn muti on 可以运行下去 但是特别慢 不知怎么回事
多谢各位大神
