new
gen zone brick p0 0 0 0 p1 300 0 0 p2 0 50 0 p3 0 0 20 size 30 15 15 group 1
gen zone brick p0 0 0 20 p1 300 0 20 p2 0 50 20 p3 0 0 25 size 30 15 10 group 2
gen zone brick p0 0 0 25 p1 300 0 25 p2 0 50 25 p3 0 0 28 size 30 15 15 group 3
gen zone brick p0 0 0 28 p1 300 0 28 p2 0 50 28 p3 0 0 30 size 30 15 10 group 4
gen zone brick p0 0 0 30 p1 300 0 30 p2 0 50 30 p3 0 0 34 size 30 15 10 group 5
gen zone brick p0 0 0 34 p1 300 0 34 p2 0 50 34 p3 0 0 38 size 30 15 10 group 6
gen zone brick p0 0 0 38 p1 300 0 38 p2 0 50 38 p3 0 0 39 size 30 15 10 group 7
gen zone brick p0 0 0 39 p1 300 0 39 p2 0 50 39 p3 0 0 59 size 30 15 20 group 8
plot add surface red range group 1
plot add surface black range group 2
plot add surface blue range group 3
plot add surface black range group 4
plot add surface Cyan range group 5
plot add surface blue range group 6
plot add surface green range group 7
plot add surface YELLOW range group 8
model mohr
prop bulk 4.8e10 shear 9.6e9 coh 5e6 fric 29 tens 6e6 ran group 1
prop bulk 3.57e9 shear 7.14e8 coh 8e6 fric 30 tens 2e6 ran group 2
prop bulk 3.45e10 shear 6.85e9 coh 5e6 fric 24.4 tens 5e6 ran group 3
prop bulk 1.97e10 shear 1.14e10 coh 4e6 fric 20.9 tens 4e6 ran group 4
prop bulk 1.65e10 shear 1.25e10 coh 8e6 fric 30.5 tens 25.7e6 ran group 5
prop bulk 2.66e10 shear 1.22e10 coh 6e6 fric 25.6 tens 9.4e6 ran group 6
prop bulk 3.57e9 shear 7.14e8 coh 8e6 fric 30 tens 2e6 ran group 7
prop bulk 4.8e10 shear 9.6e9 coh 5e6 fric 29 tens 6e6 ran group 8
ini dens 2400 ran group 1
ini dens 1400 ran group 2
ini dens 2500 ran group 3
ini dens 2400 ran group 4
ini dens 2600 ran group 5
ini dens 2560 ran group 6
ini dens 1400 ran group 7
ini dens 2400 ran group 8
fix x ran x -.1 .1
fix x ran x 299.9 300 .1
fix y ran y -.1 .1
fix y ran y 49.9 50.1
fix z ran z -.1 .1
set grav 0 0 -9.8
app nstress -8e6 ran x 0 300 y 0 50 z 59
hist unbal
;solve
step 1000
plot con szz ou