new
config dyn
set dyn off
gen zone brick p0 0 0 0 p1 30 0 0 p2 0 1 0 p3 0 0 10 size 30 1 10
gen zone brick p0 0 0 10 p1 20 0 10 p2 0 1 10 p3 0 0 20 p4 20 1 10 p5 0 1 20 p6 10 0 20 p7 10 1 20 size 20 1 10
plot zone
fix x y z range z -0.1 0.1
fix x range x 29.9 30.1
fix x range x -0.1 0.1
fix y
model elas
prop density 2000 bulk 3e9 shear 1e9
set gravity 0 0 -10.0
solve
ini xvel 0 yvel 0 zvel 0
ini xdis 0 ydis 0 zdis 0
model mech mohr
pro density 2000 bulk 1e8 shear 3e7 coh 1e3 dila 20 tens 1e10 fric 20
set dyn on
def wave
local per=0.1571
wave=1e5*(1.0-cos(2*pi*dytime/per))
end
free x y z ran x 19.9 30.1
apply nquiet squiet dquiet ran z -0.1 0.1
apply nstress 10.0 hist @wave ran x 24.9 25.1
apply ff
group ff_corner
group ff_side ran x 0 30
group ff_side ran y 0 1
group main_grid ran x 0 30 y 0 1
set dyn time=30
hist reset
hist unbal
hist dytime
hist gp xvel 5 0 20
hist gp xvel 5 1 20
hist gp xvel 15 0 12.5
hist gp xvel 15 1 12.5
hist gp xvel 25 0 10
hist gp xvel 25 1 10
solve age 0.015
