接触面计算结果显示的问题

yjw_413

本人最近遇到接触面的建模和计算问题，参照陈育民的《FLAC/FLAC3D基础与工程实例》一书P141，例9-4料仓下料的例子，命令流完全参照教材，建好模型后，计算，但是使用pl con zdisp ou on magf 100，却不能显示书上的接触面滑移的效果，问什么呢？还有请问接触面模拟的分离是指什么样的情况呢？以下是我计算显示的结果，另一幅是书上的！

1. 
   
2. 
   
3. new
   
4. ;创建模型
   
5. ;Creat Material Zones
   
6. gen zone brick size 5 5 5 &
   
7. p0 0,0,0 p1 3 0 0 p2 0 3 0 p3 0 0 5 &
   
8. p4 3 3 0 p5 0 5 5 p6 5 0 5 p7 5 5 5
   
9. gen zone brick size 5 5 5 p0 0 0 5 edge 5.0
   
10. group Material
    
11. ;Create Bin Zones
    
12. gen zone brick size 1 5 5 &
    
13. p0 3 0 0 p1 add 3 0 0 p2 add 0 3 0 & ;add means the relative cooridinates
    
14. p3 add 2 0 5 p4 add 3 6 0 p5 add 2 5 5 &
    
15. p6 add 3 0 5 p7 add 3 6 5
    
16. gen zone brick size 1 5 5 &
    
17. p0 5 0 5 p1 add 1 0 0 p2 add 0 5 0 &
    
18. p3 add 0 0 5 p4 add 1 6 0 p5 add 0 5 5 &
    
19. p6 add 1 0 5 p7 add 1 6 5
    
20. gen zone brick size 5 1 5 &
    
21. p0 0 3 0 p1 add 3 0 0 p2 add 0 3 0 &
    
22. p3 add 0 2 5 p4 add 6 3 0 p5 add 0 3 5 &
    
23. p6 add 5 2 5 p7 add 6 3 5
    
24. gen zone brick size 5 1 5 &
    
25. p0 0 5 5 p1 add 5 0 0 p2 add 0 1 0 &
    
26. p3 add 0 0 5 p4 add 6 1 0 p5 add 0 1 5 &
    
27. p6 add 5 0 5 p7 add 6 1 5
    
28. group Bin range group Material not
    
29. ;Create named range synonyms
    
30. range name Bin group Bin
    
31. range name Material group Material
    
32. ;Assign models to groups
    
33. model mohr range Material
    
34. model elas range Bin
    
35. gen separate  Material
    
36. interface 1 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 above
    
37. interface 2 wrap Material Bin range plane ori 0 0 0 normal 1 -1 0 below
    
38. int 1 maxedge 0.55
    
39. int 2 maxedge 0.55
    
40. ;Assign properties
    
41. prop shear 1e8 bulk 2e8 fric 30 range Material
    
42. prop shear 1e8 bulk 2e8 range Bin
    
43. ini den 2000
    
44. int 1 prop ks 2e9 kn 2e9 fric 15
    
45. int 2 prop ks 2e9 kn 2e9 fric 15
    
46. ;Assign Boundary Conditions
    
47. fix x range x -0.1 0.1 any x 5.9 6.1 any
    
48. fix y range y -0.1 0.1 any y 5.9 6.1 any
    
49. fix z range z -0.1 0.1
    
50. ;Monitor histories
    
51. hist unbal
    
52. hist gp zdisp 6 6 10
    
53. hist gp zdisp 0 0 10
    
54. hist gp zdisp 0 0 0
    
55. ;Setting
    
56. set large
    
57. set grav 0 0 -10
    
58. ;Cycling
    
59. step 4000
    
60. save bin.sav
    
61. pl con zdisp ou on magf 100

复制代码