new
gen zone brick size 90 50 60 &
p0=(0 0 0 ) p1=(90 0 0) p2=(0 50 0) p3=(0 0 60)
plot
show
add surface yellow
add axes black
model mohr
property bulk 8e9 shear 5e9 friction 44
proterty cohesion 2.72e5 tension 2e5
set gravity 0,0,-10
ini density=2500
fix z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 89.9 90.1
fix y range y -0.1 0.1
fix y range y 49.9 50.1
;grz=10*2500*500=12500000
ini szz -1.4e7 grad 0,0,2.5e4
ini syy -7e6 grad 0,0,1.25e4
ini sxx -7e6 grad 0,0,1.25e4
apply szz -1.25e7 range z 59.9 60.1
apply szz -1.4e7 range z -0.1 0.1
apply sxx -7e6 grad 0 0 1.25e4 range x -0.1 0.1
apply sxx -7e6 grad 0 0 1.25e4 range x 89.9 90.1
apply syy -7e6 grad 0 0 1.25e4 range y -0.1 0.1
apply syy -7e6 grad 0 0 1.25e4 range y 49.9 50.1
history nstep=5
history unb
set mechanical force 50
solve
plot
history
print history
