new
config dynamic
set dynamic dt 1.0e-5
set dynamic time 0.0
set dynamic multi on
set dynamic damp rayleigh 0.02 30
gen zone brick p0 0 0 0 p1 10 0 0 p2 0 1 0 p3 0 0 1 size 50 5 5
gen zone brick p0 10.1 0 0 p1 20.1 0 0 p2 10.1 1 0 p3 10.1 0 1 size 50 5 5
group 1
gen zone brick p0 10 0 0 p1 10.1 0 0 p2 10 1 0 p3 10 0 1 size 4 10 10
group 2 range group 1 not
;plot block group
model elastic
prop bulk 3e10 shear 1.8e10 range group 1
prop bulk 2.5e10 shear 1.5e10 range group 2
ini dens 2500 range group 1
ini dens 2200 range group 2
fix x y z range x -.1 .1
def wave
wave=-1e6*sin(2*pi*1000.0*dytime)
if dytime>1.0/2000.0
wave=0.0
endif
end
apply sxx 1.0 hist wave range x 20.09 20.11
hist gp xdisp 0 0.5 1 ;1
hist gp xdisp 10.1 0.5 1 ;2
hist gp xdisp 20.1 0.5 1 ;3
hist gp xdisp 10 0.5 1 ;4
hist dytime ;5
plot create hhh
plot add hist 1 2 3 4 v 5
plot show
solve age 0.1
hist write 1 2 3 4 v 5 file 位移.txt
