abaqus接触状态仿真结果与文献中不同，这是为什么？ 110

自己仿真的结果为图1，而文献中的结果如图2，接触部位的滑移区为什么这么小。这是根据文献中的要求进行设置的。具体设置内容如下：单元类型均为C3D8。采用接触对的形式建立轴和套环过盈配合面间的接触，轴环面为从面，套环内环面为主面。接触面间的摩擦行为遵循库伦定律，即 q=μp。过盈配合面间的摩擦系数取 0.6，切向接触行为由罚函数接触算法进行控制，其中为了更好地表征粘着和滑移状态，罚函数接触算法中的弹性滑移极限设为 0.0001。仿真分析时采用了线性随动强化本构模型对材料属性进行定义。为模拟四点弯曲疲劳试验，采用模型不动，在模型的B部位施加水平和垂直方向的简谐波加载力，在A处限制y和z方向的自由度，如图3所示。根据使轴与轴套配合部位的弯曲应力为220MPa，计算得到加载力幅值为1582N。附件：inp文件（重启动，共2个）、本构模型for（线性随机本构）、cae文件（abaqus 2020版）、参考文献。

20210914test.rar