请教版主及各位大峡，为何动力计算只计算2秒？

lsy

请教版主及各位大峡！
最近在学FLAC2D，在计算一动力分析时，为何只计算2秒钟？敬请各位大侠指出错误！
我的数据文件：
config gw ats dyn extra 10
grid 30,20
;region-1-底部土体
gen (-40,-40) (-40,-3) (-11.1,-3) (-11.1,-40) ratio 1.0,1.0 i=1,11 j=1,11
gen (-11.1,-40) (-11.1,-3) (-7.2325,-3) (-7.2325,-40) ratio 1.0,1.0 &
       i=11,13 j=1,11
gen (-7.2325,-40) (-7.2325,-3) (7.2325,-3) (7.2325,-40) ratio 1.0,1.0 &
       i=13,19 j=1,11
gen (7.2325,-40) (7.2325,-3) (11.1,-3) (11.1,-40) ratio 1.0,1.0  &
       i=19,21 j=1,11
gen (11.1,-40) (11.1,-3) (40,-3) (40,-40) ratio 1.0,1.0  &
       i=21,31 j=1,11
;region-1-material-底部土体
model mohr i=1,30 j=1,10
;region-2-底部以上，天然以下土体
gen (-40,-3) (-40,0) (-11.1,0) (-11.1,-3) ratio 1.0,1.0 i=1,11 j=11,13
gen (-11.1,-3) (-11.1,0) (-7.2325,0) (-7.2325,-3) ratio 1.0,1.0 &
       i=11,13 j=11,13
gen (-7.2325,-3) (-7.2325,0) (7.2325,0) (7.2325,-3) ratio 1.0,1.0 &
       i=13,19 j=11,13
gen (7.2325,-3) (7.2325,0) (11.1,0) (11.1,-3) ratio 1.0,1.0  &
       i=19,21 j=11,13
gen (11.1,-3) (11.1,0) (40,0) (40,-3) ratio 1.0,1.0  &
       i=21,31 j=11,13
;region-2-material-底部以上，天然以下土体
model mohr i=1,30 j=11,12
;region-3-路基填土
gen (-11.1,0) (-3.6,5) (-2.4,5) (-7.2325,0) ratio 1.0,1.0 i=11,13 j=13,19
gen (-7.2325,0) (-2.4,5) (2.4,5) (7.2325,0) ratio 1.0,1.0 &
       i=13,19 j=13,19
gen (7.2325,0) (2.4,5) (3.6,5) (11.1,0) ratio 1.0,1.0  &
       i=19,21 j=13,19
;region-3-material-路基填土
model mohr i=11,20 j=13,18
;region-2-material-底部以上，天然以下土体
model mohr i=1,30 j=11,12
;region-4-道碴
gen (-2.4,5.0) (-1.7,5.4) (1.7,5.4) (2.4,5.0) ratio 1.0,1.0 &
       i=13,19 j=19,21
;region-4-material-道碴
model mohr i=13,18 j=19,20


; 定义材料
group 'foundation:soil 2'  i=1,30 j=1,10
model mohr group 'foundation:soil 2'
prop density=1.92 bulk=7.5e5 shear=3.46e5 cohesion=500.0 &
   tension=0.0 friction=45.0 dilation=0.0 group 'foundation:soil 2'
group 'foundation:soil 1'  i=1,30 j=11,12
model mohr group 'foundation:soil 1'
prop density=1.92 bulk=4.17e4 shear=1.92e4 cohesion=60.0 &
   tension=0.0 friction=25.0 dilation=0.0 group 'foundation:soil 1'
group 'sugrade:soil'  i=11,20 j=13,18
model mohr group 'sugrade:soil'
prop density=2.1 bulk=4.17e4 shear=1.92e4 cohesion=60.0 &
   tension=0.0 friction=25.0 dilation=0.0 group 'sugrade:soil'
group 'ballast'  i=13,18 j=19,20
model elastic group 'ballast'
prop density=2.1 bulk=2.46e7 shear=1.13e7      group 'ballast'
; 设置地下水
prop por=0.3 perm=3.06e-8 group 'foundation:soil 2'
prop por=0.3 perm=3.06e-8 group 'foundation:soil 1'
prop por=0.3 perm=9.18e-6 group 'sugrade:soil'
prop por=0.7 perm=3e-1 group 'ballast'
; 定义边界和初始条件
fix x y j 1
fix x i 1
fix x i 31
set gravity=9.8
set flow=off
water density=1.0
water bulk=0.0   ;2e6
set dyn=off
set=large
fix pp i 1 11 j 13
fix pp i 21 31 j 13
fix pp i 11 j 13 19
fix pp i 21 j 13 19
fix pp i 11 21 j 19
fix pp i 1
fix pp i 31
hist 1 unbal
hist 2 ydisp i=16 j=21
set mechanical=on
set sratio=0.0001
solve
save subgrade1.sav


; *** 输入地震波***
water bulk=2e6
set dyn=on
set=large
hist 100 read acc.his
hist write 100 table 100
call filter.fis
set filter_in 100
set filter_out 101
set fc=5.0
filter
call int.fis
set int_in 101
set int_out 102
integrate   
initial xdisp 0 ydisp 0
initial xvel 0 yvel 0
history 3 pp i=16, j=4
history 4 pp i=16, j=11
history 5 pp i=16, j=13
history 6 syy i=16, j=11
history 7 syy i=16, j=13
history 8 ydisp  i=16, j=11
history 9 ydisp  i=16, j=13
history 10 ydisp i=16, j=21   
history 11 xvel i=16, j=21
history 12 dytime
save subgrade2.sav

;***加入动力条件***
apply ffield
apply sxy -5154.88 hist table 102 from 1,1 to 31,1
apply xquiet yquiet from 1,1 to 31,1
save subgrade3.sav

;***无阻尼的情况***
solve dytime 30.02
hist write 11 vs 12 table 200
call tab_ind.fis
set fft_inp1=200 fft_inp2=201
tab_ind
call fft.fis
fftransform
save subgrade4.sav

sch

你的代码中错误很多啊！首先释放约束边界，在施加粘性边界！然后施加输入，最后应用自由场，顺序不要搞错，另外solve dytime 30.02 应改为solve age 30.02 。

lsy

老兄，你说的这个我试过了，好象不是这个原因！！！

sch

是不是这个原因倒是不知道，但是你的那些错误却是肯定的。

lsy

我的问题已解决！谢谢楼上的仁兄参与讨论！

benjackxu

LZ可否把你解决该问题的心得体会写下来，以便后来者少走弯路？

lsy

就像SCH所说的一样将solve dytime 30.02 应改为solve age 30.02 。

ransheng77

关于solve age能不能有更详细的讲解呢？