关于最简单的摩擦滑动问题

zhly87

如上图所示，做了个简单的滑动例子。我想知道在摩擦系数和粘聚力分别为多大时接触面发生滑动，但是无论摩擦系数和粘聚力降低为多少，模型计算始终收敛，这是不是意味着接触面没滑动？
不知道是我对摩擦单元的定义有问题，还是对于是否滑动的判断方法有误，望高手指点，谢谢！
以下是这个例子的命令流~

1. 
   
2. 
   
3. /prep7
   
4. block,0,100,0,100,0,100
   
5. k,,0,50,100
   
6. k,,100,50,100
   
7. k,,100,100,50
   
8. k,,0,100,50
   
9. a,9,10,11,12
   
10. vsba,1,7
    
11. vdele,2
    
12. a,9,10,11,12
    
13. va,2,8,9,10,11
    
14. 
    
15. 
    
16. et,1,solid45
    
17. et,2,170
    
18. et,3,173
    
19. mp,ex,1,2.6e10
    
20. mp,prxy,1,0.167
    
21. mp,dens,1,2000
    
22. mp,mu,3,0 !摩擦系数为0
    
23. R,1,,,, ,,
    
24. RMORE,,,, ,,
    
25. RMORE,0,,, ,, !粘聚力为0
    
26. 
    
27. 
    
28. 
    
29. TYPE,1
    
30. MAT,1
    
31. REAL,1
    
32. ESYS,0
    
33. MSHAPE,1,3D
    
34. MSHKEY,0
    
35. CM,_Y,VOLU
    
36. VSEL,s,,,1
    
37. vsel,a,,,3
    
38. CM,_Y1,VOLU
    
39. CHKMSH,'VOLU'
    
40. CMSEL,S,_Y
    
41. VMESH,_Y1
    
42. CMDELE,_Y
    
43. CMDELE,_Y1
    
44. CMDELE,_Y2
    
45. 
    
46. 
    
47. 
    
48. ！定义接触对
    
49. asel,s,,,7
    
50. cm,target,area
    
51. type,2
    
52. nsla,s,1
    
53. esurf,all
    
54. esel,s,type,,2
    
55. eplot
    
56. allsel
    
57. 
    
58. asel,s,,,2
    
59. aplot
    
60. cm,contact,area
    
61. type,3
    
62. nsla,s,1
    
63. esurf,all
    
64. esel,s,type,,3
    
65. eplot
    
66. allsel
    
67. 
    
68. 
    
69. NSEL,S,,,140
    
70. NSEL,A,,,141
    
71. NSEL,A,,,103
    
72. NSEL,A,,,93
    
73. D,all, , , , , ,UX,UY,UZ, , ,
    
74. acel,,9.8,
    
75. finish
    
76. 
    
77. /solu
    
78. NSUBST,5,20,1
    
79. nlgeom,1
    
80. 
    
81. allsel
    
82. solve
    

复制代码

zhly87

期待高手解答~~

jacky_cai

这个问题解决了吗

satana123

楼主解决了吗