案例二：竖直开挖

model new
zone create brick size 10 10 10
zone cmodel assign elastic
zone property bulk 3e6 shear 1e6 density 1500

zone face skin ;break-angle 45

zone face apply velocity-normal 0.0 range group 'West' or 'East'
zone face apply velocity-normal 0.0 range group 'South' or 'North'
zone face apply velocity-normal 0.0 range group 'Bottom'

model gravity 0 0 -10

;zone initialize-stresses ratio 0.3

history interval 5
model history mechanical unbalanced-maximum

model solve ratio-local 1e-4
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0
model save 'initial'
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
model restore 'initial'
zone group "Space" range position-x 0 5 position-y 4 7 position-z 5 10
zone relax excavate range group "Space" ;将属于space的土体开挖掉
model solve  ;
model save "kaiwa.f3sav"  ;保存文件