apply和fix应该都可以固定边界吧！
我将书上的例子有apply方式换成fix之后，就不能绘制速率图了不知道为什么？
model mohr
prop bulk 5e9 shear 3e9 fric 35
ini density 2500
set gravity 0,0,-10
fix x y z range z -0.1 0.1
ini szz -5.0e6 grad 0,0,2.5e4
ini syy -2.5e6 grad 0,0,1.25e4
ini sxx -2.5e6 grad 0,0,1.25e4
apply szz -4.5e6 range z 19.9 20.1
;apply szz -5.0e6 range z -0.1 0.1
;apply sxx -2.5e6 grad 0,0,1.25e4 range x -0.1 0.1
;apply sxx -2.5e6 grad 0,0,1.25e4 range x 19.9 20.1
;apply syy -2.5e6 grad 0,0,1.25e4 range y -0.1 0.1
;apply syy -2.5e6 grad 0,0,1.25e4 range y 19.9 20.1
fix x ra x .1 -.1
fix y ra y .1 -.1
fix z ra z .1 -.1
fix x ra x 19.9 20.1                         s
fix y ra y 19.9 20.1

pl creat yingliyuntu
pl add con szz

pl creat sudutu
pl add grid red vel blue

pl show

step 1
刚刚开始学习，不懂。还有apply和ini施加应力时候有什么不同~~

这会儿自己弄好了，自己改动的程序有点错误。