动力分析-2

时程转换与地震波滤波

1.滤波机理

过滤最高频率（波长变小）→增大网格尺寸→减少网格数量→减少计算时间

2.基线校正

输入加速度时程曲线，进行积分得到的最终速度和最终位移不为0，则在动力计算结束时，

模型底部出现继续的速度和残余的位移，此时需要对加速度时程进行基线校正；

如图所示-----滤波以及基线校正后

3.动态多步（分区使用时间步计算）

zone dynamic multi-step on

计算模型中存在刚度差异较大、模型网格尺寸不均匀的情况时，flac3d可以采用“动态多步”的过程来有效减小计算所需要的时间。

model title 'Shear wave applied to a stiff wall in a soft soil with dyn. multistep'
;program threads 8
model configure dynamic
zone create brick size 40 10 20 ...
             point 0 (-5,0,0) point 1 (15,0,0) ...
             point 2 (-5,5,0) point 3 (-5,0,10)
zone delete range position-x -5 5 position-z 5 10
zone group "soil"
zone group "wall" range position-x 5 6 position-z 5 10
;
; import gridpoint file
;zone import '40.f3grid'
zone face skin ; Label model boundaries
;
; Assign model and properties
zone cmodel assign elastic

zone property young 20e4 poisson 0.2 density 2000 ;range group 'soil'
zone property young 30e6 poisson 0.2 density 2000 range group 'wall'
;
; Boundary conditions
zone face apply velocity-z 0 range group 'East' or 'West2'
zone face apply velocity-normal 0 range group 'North' or 'South'


fish define wave
    wave = math.sin(math.pi * dynamic.time.total / 16.0)
end
zone face apply velocity-x = 1 fish @wave range group 'Bottom'
zone face apply velocity-normal 0 range group 'Bottom'


;
; Histories
zone history name='xvelb' velocity-x position (5,2, 0)
zone history name='xvelt' velocity-x position (5,2,10)
zone history name='zvelt' velocity-z position (5,2,10)
model history name='time' dynamic time-total
;set the dynamic timestep
;model dynamic timestep fix 9e-5
; Set multi-step status, and cycle to 32, timing how long that takes

zone dynamic multi-step on
[global start = time.clock] 
model solve time-total 32
[global time = 0.01 * (time.clock - start)]
[time]