pnt = zone_head loop while pnt # null
val = _k0 * z_szz(pnt) + (_k0-1.) * z_pp(pnt) z_sxx(pnt) = val z_syy(pnt) = val pnt=z_next(pnt) end_loop end ini_conf solve save 11-8.sav
2. 动力分析
动力计算的命令文件见例11.8所示。在本例中对计算结果采取了分阶段保存的方法，定义了一个FISH函数，其功能是在动力计算完成特定的时刻即保存一个文件，这样便于计算结束后对中间结果的调用。
例11.8：振动台模拟的动力计算
rest 11-7.sav set dyn on fluid on ini fmod 2e8 set fluid pcut on model finn ran gro sand
prop bulk=3e7 shear=1e7 fric=35 ran gro sand ini dens 2000 ran gro sand
prop ff_latency=50 ff_switch = 0 ff_c1=0.8 ff_c2=0.79 ff_c3=0.45 ff_c4=0.73 ran gro sand ;扭剪试验结果 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 free x
apply xvel = 1.0 hist sine_wave ran z -9.9 -10.1 apply xvel = 1.0 hist sine_wave ran x -.1 .1 apply xvel = 1.0 hist sine_wave ran x 29.9 30.1 set dyn damp local .314 call ppr.dat4 set hist_rep 100 set large ;set dyn dt 3e-4 set mech rat 1e-20 def solve_ages loop n(1,39)
save_file = '11-8_'+string(n)+'s.sav' command sol age n ;save save_file endcommand endloop end solve_ages save save_file
hist write 20 21 22 23 24 30 31 32 33 34 vs 2 file 10-8_Outfile_pp.txt hist write 40 41 42 43 44 vs 2 file 10-8_Outfile_xdis.txt
11.10.3 计算结果与分析
将Finn模型应用不同的阻尼形式得到的结果进行了汇总，见图11-21所示。可以发现，利用不同阻尼形式得到的计算存在一定的差异，但超孔压时程和平均应力时程曲线的变化趋势基本一致。
本书中仅提供了不耦合模式下的计算结果，读者可以自行开展耦合模式下的动力计算，可以发现，在耦合模式下，在动荷载作用到20s后，随着荷载幅值的减小，超孔隙水压力会逐渐消散，而在不耦合模式下，积累的超孔压将不会消散。因此在实际分析中，如果需要考虑动孔压在震动过程中的积累和消散效应，则在流固耦合的同时也需要将渗流计算打开。
ppr..dat文件的作用是监测计算过程中一些变量的值，由于监测的内容较多，为了使命令的可读性更强，所以将历史变量的监测单独形成了一个文件。
4
200p / kPa150100 局部阻尼 瑞利阻尼 滞后阻尼0?m' / kPa-50-100051015202530 局部阻尼 瑞利阻尼 滞后阻尼 3540时间 / s
图11-21 阻尼形式对计算结果的影响 (不耦合)
11.11 本章小结
本章以FLAC3D为例，对动力分析中边界条件的设置、阻尼形式的选择、网格尺寸的要求、动荷载的调整、土体的液化等内容进行了介绍。FLAC / FLAC3D软件在动力分析方面具有强大的功能，同时也具有一定的难度。读者在进行动力分析前要对静力计算非常熟悉，因为初始应力计算的正确才能保证后续动力计算的顺利开展。
FLAC3D的动力计算往往要耗费很多的时间，主要原因是由于有限差分法本身需要较小的时间步，因此建议读者在进行动力分析前首先进行简单模型的分析，熟悉边界条件、荷载施加、阻尼选择等方面的基本内容，虽然在这些“小例子”中会花费读者一定的时间，但是这样会大大提高分析的效率。