在动力分析中的例题:
; --- Free-field test ---
;
config dyn
;
def iniwave
per = 0.01
end
def wave
wave = 0.5 * (1.0 - cos (2*pi*dytime/per))
end
iniwave
;
gen zone brick size 6 3 2
gen zone brick size 2 3 2 p0 0 0 2
gen zone brick size 2 3 2 p0 4 0 2
gen zone wedge size 1 3 2 p0 2 0 2
gen zone wedge size 1 3 2 p0 4 3 2 p1 3 3 2 p2 4 0 2 p3 4 3 4 &
p4 3 0 2 p5 4 0 4
;
model elastic
prop bulk 66667 shear 40000
ini dens 0.0025
;
; --- static state ---
;
set grav 0 0 -10
;
fix x range x -0.01 0.01
fix x range x 5.99 6.01
fix y range y -0.01 0.01
fix y range y 2.99 3.01
fix z range z -0.1 0.1
;
set dyn off
;
hist unbal
hist dytime
hist gp zdis 2 1 0
hist gp zdis 2 1 2.0
hist gp zdis 2 1 5.0
;
solve
;
save ff0.sav
; --- dynamic run A: applied velocity at base ---
;
set dyn on
; --- set dyn b.c. at base ---
free x y z ran z -0.1 0.1
apply nquiet squiet dquiet ran z -0.1 0.1
apply dstress 1.0 hist wave ran z -0.1 0.1
; --- FF ---
;apply ff
;
group ff_corner 这四个group命令不是很清楚。
group ff_side ran x 0 6
group ff_side ran y 0 3
group main_grid ran x 0 6 y 0 3 ;
;
set dyn time = 0
hist reset
;
hist unbal
hist dytime
;
; main grid
hist gp xvel 2 1 0
hist gp xvel 2 1 5.0
;
; corner FF
hist gp xvel -1 -1 0
hist gp xvel -1 -1 5.0
;
; side FF parallel to y
hist gp xvel -1 0 0
hist gp xvel -1 0 5.0
;
; side FF parallel to x
hist gp xvel 2 -1 0
hist gp xvel 2 -1 5.0
;
solve age 0.015
save noff1.sav