gen zone brick size 12 12 4 p0 0.05 0 0.15 p1 0.35 0 0.15 p2 0.05 0.3 0.15 p3 0.05 0 0.25
gen zone brick size 16 12 4 p1 0.40 0 0 p2 0 0.30 0 p3 0 0 0.10
range name bot z 0 0.1
range name top z 0.15 0.25
interface 1 face range z 0.1
int 1 prop ks 64.8e9 kn 64.8e9 fric 20 dil 6 tension 0 bslip=on
ini z add -0.05 range top
plot bl group interface white axes black
pause
model mohr
prop bulk 2.3e8 sh 2.46e7 fric 20 cohesion 1e10 tension 0
set gravity 0,0,-9.81
initial density 1600
fix x y z range z 0
fix x range x 0
fix x range x 0.4
fix y range y 0
fix y range y 0.3
fix x range x 0.09 0.11
fix x range x 0.29 0.31
apply nstress -25e3 range z 0.2
step 0
plot contour szz interface white axes black
solve
save dsta.sav
pause
prop cohesion=0 tension=0
free x range x 0.09 0.11
free x range x 0.29 0.31
ini xvel 2.5e-6 range top
fix xvel range top
ini xdis 0 ydis 0 zdis 0
step 2000