自己仿真的結果為圖1，而文獻中的結果如圖2，接觸部位的滑移區為什么這么小。這是根據文獻中的要求進行設置的。具體設置內容如下：單元類型均為C3D8。采用接觸對的形式建立軸和套環過盈配合面間的接觸，軸環面為從面，套環內環面為主面。接觸面間的摩擦行為遵循庫倫定律，即 q=μp。過盈配合面間的摩擦系數取 0.6，切向接觸行為由罰函數接觸算法進行控制，其中為了更好地表征粘著和滑移狀態，罰函數接觸算法中的彈性滑移極限設為 0.0001。仿真分析時采用了線性隨動強化本構模型對材料屬性進行定義。為模擬四點彎曲疲勞試驗，采用模型不動，在模型的B部位施加水平和垂直方向的簡諧波加載力，在A處限制y和z方向的自由度，如圖3所示。根據使軸與軸套配合部位的彎曲應力為220MPa，計算得到加載力幅值為1582N。附件：inp文件（重啟動，共2個）、本構模型for（線性隨機本構）、cae文件（abaqus 2020版）、參考文獻。

20210914test.rar