卷板时主动辊对板材的摩擦力出现正负波动？

leojcf

做的卷板成形如下图所示，板材为壳体，三个辊子为离散刚体。
上辊为主动辊，绕自身的轴线旋转，转速为3rad/s，赋于参考点质量为50，绕轴线的转动惯量为3200
两下辊为从动轴，同样赋于参考点质量为50，绕轴线的转动惯量为3200，约束参考点除了UR3外五个自由度。
成形分两步，也给上辊一定的下压量。接着给上辊一转速，带动板材连续成形，两从动下辊也发生转动。

先说一下模拟结果吧。从应力应变图看，板材的卷板成形效果很好。现在主要是研究辊子与板材间作用的法向正压力及摩擦力。
先输出辊子与板材作用的正压力CFNM，得到的结果与理论计算值较接近。
再输出辊子对板材的摩擦力。先输出主动辊的，CFS1，1是沿X轴方向的分力。可能实际摩擦力与X轴方向有一夹角，但由于角度很小，所以在X轴方向的分力接近实际摩擦力大小，因此输出CFS1更能反映摩擦力的正负。结果输出的摩擦力的结果很令人意外。主动辊，即上辊对板材的摩擦力理论上应该沿板材的进给方向，即沿X轴正方向，为正值。但现在输出的结果在正负之间波动。
上辊对板材作用的摩擦力曲线
再来看看两从动下辊对板材的摩擦力曲线，同样是沿X轴方向的分力，同样也是出现了在正负之间波动。理论上两从动辊是靠板材的摩擦力发生转动。摩擦力应该是负值才对。
两下辊对板材的摩擦力曲线
在前0.50s，属于上下辊不转动，上辊下移的过程。此时摩擦力的情况方向大小都比较合理。也符合库伦摩擦定律。

先来说说自己接触的设置吧。采用库伦摩擦模型，即罚函数摩擦公式，摩擦系数0.15；法向采用默认的硬接触；面面显式接触，有限滑移。然后就出现上述问题了。
把模型传上来，各位大侠们帮忙看看怎么回事。一直没能解决这个问题。

leojcf

在旋转阶段，取主动辊及从动辊上的任意一点的线速度，发现它们的线速度相等，也就是从动辊在主动辊的带动下已经发生了旋转，并且转速相同。当然，为了解决这个问题，自己也做过很多尝试。其中有，把三个辊子都设为主动辊，网格细分，降低转速。把辊子由离散刚体改为解板刚体，这样一来就不用划分网格。给以两个从动下辊很大的转动惯量，目的是使两从动的转动变得相对困难。最近怀疑是库伦摩擦定律和动态显式算法导致的，理由是：取辊子与板材上任意一点的线速度曲线，如下图所示：

由图可见，板材与辊子的线速度非常接近，但是板材的线速度有波动，时而比辊子大，时而小，当然变化很微小。实际上板材线速度不可能比主动辊大。这应该是动态算法引起的。而正是这种线速度的微小变化，导致板材与辊子间的相对运动方向相反，经典库伦摩擦对相对运动方向的敏感性立刻表现在摩擦力的正负上。后来看论文看到修正的库伦摩擦定律也许能解决问题，但一直没弄明白如何使用修正的库伦摩擦定律，似乎需要编程。

leojcf

模型稍后传上，之前传的inp自己试了一下导入后发现有点小问题，还是传cae文件吧。先吃饭

leojcf

3# leojcf