动力分析-3

1.边界：抵消或者吸收振动波的反射响应； →  选择并施加合适的边界条件

特殊边界减小反射：①静态/粘性边界 ②自由场边界

2.软土地基和刚性地基--地震荷载的施加

①对于模型底座地层为软土时，外部施加的动荷载程只能选择应力或者点时程，

因为底部边界需要粘性条件，四周采用自由场边界；

②刚性基础一般选择速度或者加时程，自由场边界即可；内部动荷载无此约定。

3.静态边界注意事项-关于荷载源（经验）

①对于震动荷载源于模型内部（ 如地铁、隧道中的列车振动、内部爆破等问题）， 可以将震动荷载直接施加在节点上，这种情况下使用静态边界可以有效减小人工上的反射现象， 并且不需要再施加自由场边界。

②如果动力源在顶部或者底，则用自由场更能抵消反射对结果的影响；

注意：关于荷载施加过程中的计算 ：动力学分析过 程中应避免静力荷载的变化。

4.静态（粘性）边界条件的一个案例

zone face apply quiet-dip range group 'Bottom'
zone face apply quiet-strike range group 'Bottom'

model new
model largestrain off
model config dynamic
zone create brick size 1,1,50
zone face skin ; Name boundaries
zone cmodel assign elastic
zone property shear 1e7 bulk 2e7 density 1000
;
zone gridpoint fix velocity-z ;固定z方向的位移
zone face apply quiet-dip range group 'Bottom'
zone face apply quiet-strike range group 'Bottom'

fish define wave
    wave = 0.5 * (1.0 - math.cos(8.0*math.pi*math.min(zone.dynamic.time.total,0.25)))
end

zone face apply stress-xz -2e5 fish @wave range group 'Bottom'

zone history name='bot' velocity-x position (0,0, 0)
zone history name='mid' velocity-x position (0,0,25)
zone history name='top' velocity-x position (0,0,50)
model history name='time' dynamic time-total
;
model solve time-total 1.8
model save "40.f3sav"

5.自由场边界：通过在模型边界侧面四周生成二维网格，角落生成一维网格的方法来实现；

zone dynamic free-field on

model new
model largestrain off
fish automatic-create off
model title "Shear wave loading of a model mech with free-field boundaries"
model config dynamic
zone create brick size 6 3 2
zone create brick size 2 3 2 point 0 (0,0,2)
zone create brick size 2 3 2 point 0 (4,0,2)
zone create wedge size 1 3 2 point 0 (2,0,2)
zone create wedge size 1 3 2 point 0 (4,3,2) point 1 (3,3,2) point 2 (4,0,2) point 3 (4,3,4)  point 4 (3,0,2) point 5 (4,0,4)
zone face skin
;; Assign model and properties
zone cmodel assign elastic
zone property bulk 66667 shear 40000 density 0.0025
; Set Initial Conditions
model gravity 10
zone initialize-stresses ratio 0.25
; Set Boundary Conditions - Roller boundaries
zone face apply velocity-normal 0 range group 'West2' or 'East2'
zone face apply velocity-normal 0 range group 'North' or 'South'
zone face apply velocity-normal 0 range group 'Bottom'
; Static equilibrium
model dynamic active off
model solve convergence 1
model save "FreeFieldInit"

; --- dynamic run A: applied velocity at base ---
zone face apply-remove range group 'Bottom'  ; Remove velocity boundary
zone face apply quiet range group 'Bottom' ; Add quiet boundary

fish define wave
wave = 0.5 * (1.0 - math.cos(math.pi*zone.dynamic.time.total*200.0))
end

zone face apply stress-dip 1.0 fish @wave range group 'Bottom' ; Time varying stress boundary
zone dynamic free-field on ; Free field boundaries
; Solve to time 0.015
model dynamic active on
model solve time-total 0.015
model save 'FreeField'