案例一：路堤填筑

model new
zone import "case1.f3grid"
zone cmodel assign null range group 'fill-1' or  'fill-2'
zone cmodel assign mohr-coulomb range group 'soil-1'
zone property young 4e6 poisson 0.35 cohesion 5000 tension 40 friction 10 density 1800 range group 'soil-1'
zone cmodel assign mohr-coulomb range group 'soil-2'
zone property young 15e6 poisson 0.35 cohesion 10 tension 30 friction 25 density 1800  range group 'soil-2'
;shear 234.9e6 bulk 5e6 
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 cmodel assign mohr-coulomb range group 'fill-1'
zone property young 8e6 poisson 0.35 cohesion 10000 tension 40 friction 10 density 1800 range group 'fill-1'
model solve
model save 'fill-1.sav'
;
zone cmodel assign mohr-coulomb range group 'fill-2'
zone property young 8e6 poisson 0.35 cohesion 10000 tension 40 friction 10 density 1800 range group 'fill-2'
model solve
model save 'fill-2.sav'