本帖最后由 zhly87 于 2012-4-27 17:21 编辑
如上图所示,做了个简单的滑动例子。我想知道在摩擦系数和粘聚力分别为多大时接触面发生滑动,但是无论摩擦系数和粘聚力降低为多少,模型计算始终收敛,这是不是意味着接触面没滑动?
不知道是我对摩擦单元的定义有问题,还是对于是否滑动的判断方法有误,望高手指点,谢谢!
以下是这个例子的命令流~
- /prep7
- block,0,100,0,100,0,100
- k,,0,50,100
- k,,100,50,100
- k,,100,100,50
- k,,0,100,50
- a,9,10,11,12
- vsba,1,7
- vdele,2
- a,9,10,11,12
- va,2,8,9,10,11
- et,1,solid45
- et,2,170
- et,3,173
- mp,ex,1,2.6e10
- mp,prxy,1,0.167
- mp,dens,1,2000
- mp,mu,3,0 !摩擦系数为0
- R,1,,,, ,,
- RMORE,,,, ,,
- RMORE,0,,, ,, !粘聚力为0
- TYPE,1
- MAT,1
- REAL,1
- ESYS,0
- MSHAPE,1,3D
- MSHKEY,0
- CM,_Y,VOLU
- VSEL,s,,,1
- vsel,a,,,3
- CM,_Y1,VOLU
- CHKMSH,'VOLU'
- CMSEL,S,_Y
- VMESH,_Y1
- CMDELE,_Y
- CMDELE,_Y1
- CMDELE,_Y2
- !定义接触对
- asel,s,,,7
- cm,target,area
- type,2
- nsla,s,1
- esurf,all
- esel,s,type,,2
- eplot
- allsel
- asel,s,,,2
- aplot
- cm,contact,area
- type,3
- nsla,s,1
- esurf,all
- esel,s,type,,3
- eplot
- allsel
- NSEL,S,,,140
- NSEL,A,,,141
- NSEL,A,,,103
- NSEL,A,,,93
- D,all, , , , , ,UX,UY,UZ, , ,
- acel,,9.8,
- finish
- /solu
- NSUBST,5,20,1
- nlgeom,1
- allsel
- solve
复制代码
|