我把问题蜕化为非常简单的一个地基模型，，是自由场动力问题计算。
底部给一个水平正弦波；
现在的问题是： 模型对参数比较挑剔，有时算几秒程序就退出，下面是一些材料参数和 可以计算的地震时长的关系：
; for prop dens 200 sh 1.4e6 bu 4.3e6 i=1,20 j=1,10 calculated time length = 100s about
; for prop dens 200 sh 1.4e8 bu 4.3e8 i=1,20 j=1,10 calculated time length = 8s about
; for prop dens 2000 sh 1.4e8 bu 4.3e8 i=1,20 j=1,10 calculated time length = 50s less
因为材料参数一般根据试验获得，不太可能随便调整。我是在模型计算中总达不到预定地震时长，才搞了这么一个简单模型来研究参数的敏感性的。
另外，也发现降低网格密度，可以延长计算地震时长。但通过不断调整网格大小来达到能够计算的地震时长，注定是不偿失、不可取的办法。
我检查了好多遍，都没有发现明显的错误。 请大牛指点如何解决在给定的参数下，计算达不到预定时长的问题。
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; this program is to demonstrate parameter sensitivity for dynamic response calcuation by Flac2D 5.0/4.0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
config dyn
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; define the configure and material
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
grid 20 10
m e
gen 0,0 0,10 20,10 20,0 i=1,21 j=1=11
prop dens 2000 sh 1.4e8 bu 4.3e8 i=1,20 j=1,10
fix x y j=1
fix i=1
fix i=21
set dyn off
set gravity=10.0
hist dytime
hist xvel i=11 j=2
hist xvel yvel i=11 j=11
solve ; solve the static equilibrium
;return
;save equili.sav
; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; continue to solve dynamic response
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
set dyn on
set large
; define a sine wave
;
def sine_wave
sine_wave = 1.0*sin(2.0*pi*freq*dytime)
end
;
; define rayleigh damping
;
set dy_damp=rayleigh 0.05 15
set dytime 0.0 freq=10
;
; assign initial boundary
;
ini xvel=0.0 yvel=0.0 xdisp=0.0 ydisp=0.0
;
;
apply ff ;apply free-field dynamic condition
;
; apply the wave load to the bottom
;
apply yvel=0 xvel 1.0 hist sine_wave j=1
hist reset
hist dytime
hist xvel i=11 j=2
hist xvel ydis i=11 j=11
solve dyt 100
; result:
;
; for prop dens 200 sh 1.4e6 bu 4.3e6 i=1,20 j=1,10 calculated time length = 100 s about
; for prop dens 200 sh 1.4e8 bu 4.3e8 i=1,20 j=1,10 calculated time length = 8 s about
; for prop dens 2000 sh 1.4e8 bu 4.3e8 i=1,20 j=1,10 calculated time length = 50 s less
