以下是命令流和结果图片，

model mohr

prop dens 1998 young 37.5e6 poisson 0.3 coh 67e3 fric 15 ten 1 range group 1 any group 7 any group 8 any group 9 any

prop dens 1700 young 12.9e6 poisson 0.4 coh 10e3 fric 8 ten 1 range group 2

fix x range x -0.1 0.1

fix x range x 49.9 50.1

fix y range y -0.1 0.1

fix y range y 39.9 40.1

fix x y z range z -0.1 0.1

set gravity 0 0 -9.81

solve elas

save gravitysolve.sav

ini xdis 0 ydis 0 zdis 0 

ini xvel 0 yvel 0 zvel 0

ini state 0

model null range group 7 x 0 5

model null range group 8 x 0 5

model null range group 9 x 0 5

apply nstress -6.715e4 range x 4.9 5.1 group 7 any group 8 any group 9 any

solve

save gongkuang1.sav

model elas range group 8 x 0 5

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 0 5

model mohr range group 9 x 0 5

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 0 5

model null range group 7 x 5 10

model null range group 8 x 5 10

model null range group 9 x 5 10

apply nstress -6.715e4 range x 9.9 25.1 group 7 any group 8 any group 9 any

solve

save gongkuang2.sav

model elas range group 8 x 5 10

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 5 10

model mohr range group 9 x 5 10

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 5 10

model null range group 7 x 10 15

model null range group 8 x 10 15

model null range group 9 x 10 15

apply nstress -6.715e4 range x 14.9 15.1 group 7 any group 8 any group 9 any

solve

save gongkuang3.sav

model elas range group 8 x 10 15

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 10 15

model mohr range group 9 x 10 15

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 10 15

model null range group 7 x 15 20

model null range group 8 x 15 20

model null range group 9 x 15 20

apply nstress -6.715e4 range x 19.9 20.1 group 7 any group 8 any group 9 any

solve

save gongkuang4.sav

model elas range group 8 x 15 20

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 15 20

model mohr range group 9 x 15 20

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 15 20

model null range group 7 x 20 25

model null range group 8 x 20 25

model null range group 9 x 20 25

apply nstress -6.715e4 range x 24.9 25.1 group 7 any group 8 any group 9 any

solve

save gongkuang5.sav

model elas range group 8 x 20 25

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 20 25

model mohr range group 9 x 20 25

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 20 25

model null range group 7 x 25 30

model null range group 8 x 25 30

model null range group 9 x 25 30

apply nstress -6.715e4 range x 29.9 30.1 group 7 any group 8 any group 9 any

solve

save gongkuang6.sav

model elas range group 8 x 25 30

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 25 30

model mohr range group 9 x 25 30

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 25 30

model null range group 7 x 30 35

model null range group 8 x 30 35

model null range group 9 x 30 35

apply nstress -6.715e4 range x 34.9 35.1 group 7 any group 8 any group 9 any

solve

save gongkuang7.sav

model elas range group 8 x 30 35

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 30 35

model mohr range group 9 x 30 35

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 30 35

model null range group 7 x 35 40

model null range group 8 x 35 40

model null range group 9 x 35 40

apply nstress -6.715e4 range x 39.9 40.1 group 7 any group 8 any group 9 any

solve

save gongkuang8.sav

model elas range group 8 x 35 40

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 35 40

model mohr range group 9 x 35 40

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 35 40

model null range group 7 x 40 45

model null range group 8 x 40 45

model null range group 9 x 40 45

apply nstress -6.715e4 range x 44.9 45.1 group 7 any group 8 any group 9 any

solve

save gongkuang9.sav

model elas range group 8 x 40 45

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 40 45

model mohr range group 9 x 40 45

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 40 45

model null range group 7 x 45 50

model null range group 7 x 45 50

model null range group 7 x 45 50

solve

save gongkuang10.sav

model elas range group 8 x 45 50

prop dens 2500 young 30e9 poisson 0.2 range group 8 x 45 50

model mohr range group 9 x 45 50

prop dens 2400 young 4e6 poisson 0.3 coh 15e3 fric 25 ten 1 range group 9 x 45 50

solve

save gongkuang11.sav