new
gen zone brick size 30 44 3 p0=(0,0,0) p1=(60,0,0) p2=(0,220,0) p3=(0,0,6) group meiceng
gen zone brick size 30 44 2 p0=(0,0,6) p1=(60,0,6) p2=(0,220,6) p3=(0,0,9) group niyan
gen zone brick size 30 44 3 p0=(0,0,9) p1=(60,0,9) p2=(0,220,9) p3=(0,0,15) group zhongxisha
gen zone brick size 30 44 3 p0=(0,0,15) p1=(60,0,15) p2=(0,220,15) p3=(0,0,21) group niyan
gen zone brick size 30 44 5 p0=(0,0,21) p1=(60,0,21) p2=(0,220,21) p3=(0,0,36) group shazhiniyan
gen zone brick size 30 44 10 p0=(0,0,36) p1=(60,0,36) p2=(0,220,36) p3=(0,0,66) group fenshayan
gen zone brick size 30 44 3 p0=(0,0,-5) p1=(60,0,-5) p2=(0,220,-5) p3=(0,0,0) group xishayan
gen zone brick size 30 44 10 p0=(0,0,-35) p1=(60,0,-35) p2=(0,220,-35) p3=(0,0,-5) group fenshayan
pl add surf yellow
pl show
model m
call prop.txt
fix z range z -35.1 -34.9
fix x range x -0.1 0.1
fix x range x 59.9 60.1
fix y range y -0.1 0.1
fix y range y 219.9 200.1
ini szz -6.75e6 grad 0 0 2.5e4
;rh=25000*270=6.75e6
ini sxx -2.02e6 grad 0,0,1.25e4
;(u/(1-u))*rh=2.02e6
ini syy -2.02e6 grad 0,0,1.25e4
apply szz -6.75e6 range z 65.9 66.1
apply szz -9.15e6 range z -35.1 -34.9
;rh=25000*366=9.15e6
apply sxx -2.02e6 grad 0 0 1.25e4 range x -0.1 0.1
apply sxx -2.02e6 grad 0 0 1.25e4 range x 59.9 60.1
apply syy -2.02e6 grad 0 0 1.25e4 range y -0.1 0.1
apply syy -2.02e6 grad 0 0 1.25e4 range y 219.9 220.1
set gravity 0 0 -10
set large
hist id=1 unbal
solve
ini xdis 0 ydis 0 zdis 0 xv 0 yv 0 zv 0
model null range x 0 45 y 0 220 z 0 6
step 3000