设置接触，为什么两个面会粘住？

lilulilu

单元属性：
ET,4,TARGE170
ET,5,CONTA174
KEYOPT,5,1,1
KEYOPT,5,5,1
!r,1,0,0,0.01,0,0,0,
!rmore,0,0,1e+20,0,1,0,
!rmore,0,1000,.5,0,0,0

材料：
!steel
MP,EX,1,1.18e11,           !弹性模量
MP,DENS,1,7000,        !密度
MP,ALPX,1,8.7E-6,      !热膨胀系数
MP,NUXY,1,0.3,           !泊松比
MP,KXX,1,39.2,            !导热系数(加大导热系数会增加温度)
MP,C,1,470,              !比热
MP,MU,1,0.3,             !摩擦系数 ()

!pad(fenmoyejin)
MP,EX,2,1e8,           !弹性模量
MP,DENS,2,4000,        !密度
MP,ALPX,2,10E-6,      !热膨胀系数
MP,NUXY,2,0.3,           !泊松比
MP,KXX,2,40,            !导热系数
MP,C,2,335,              !比热
MP,MU,2,0.3,             !摩擦系数()

建单元：
ASEL,S,,,7!建接触面上圆柱下表面
NSLA,S,1
R,1
TYPE,4
ESURF
allsel

ASEL,S,,,2!下圆柱上表面
NSLA,S,1
R,1
TYPE,5
ESURF
allsel
设置是这样的，问题是为什么上下接触面会粘住不放，造成网格畸形？
到底是什么原因？谢谢了～

lilulilu

下面的盘旋转，上面的盘不动。
为什么下面的盘就是转不利索？
网格变形如图。哪位帮我看看，谢谢~~

zhlinzh

我也想做这样的模型，主要是转动的问题！也是转不利索的问题，请问你现在解决了吗？

imlini

太晚了，没有精力仔细看。不知道是不是两个同心圆环，里面是主动环，需要由它带动外面的。
如果是这样的话，里面那个面应该设为rigid。不过我忘记为什么你这样设置会网格畸变了。等我星期一回去看了笔记告诉你吧。

newton

你想让盘通过接触转起来？你应该修改接触单元的关键字，改为绑定接触或粗糙接触试试．

imlini

bonding接触怎么会转起来那？我查了一下以前笔记，看看是不是用penalty 算法？这种情况不能用这种算法的。因为这个算法会产生penetration，想想每个单元有一点点伸出的penetration，最后不就变成齿轮了马，所以就转不动了。当时老师给的解决方法是将里面作为rigid，设置一个pilot来控制转动。
我没有试过，不过我见过别人的例子，应该是可行的。

lilulilu

1.我粘住的情况其实是设置的热单元的问题
两个体的热单元一起建的，所以总是粘住，转不开。
2.转起来我用beam4单元在旋转中心与转动体内径连接，给中心点转速来做的
但是结果还是不太好