abaqus圆柱滚子轴承隐式动力学分析，载荷步设置要点

圆柱滚子轴承隐式动力学跑了一晚上，滚子刚转起来就崩，msg里全是接触状态反复开合的警告。静力学算法算到20%就卡死，换隐式动力学才勉强能动，载荷步的设置是核心。

分析步用Dynamic, Implicit，Hilber-Hughes-Taylor积分器，α默认取-0.05，提供少量数值阻尼抑制高频振荡，同时保持二阶精度。时间周期按实际物理时间给，内圈从静止加速到工作转速的时间、稳定运转的时间，加起来就是总时间。时间增量不要一上来就设太大，初始增量取总时间的1/1000到1/500，最小增量设1e-8，最大增量取总时间的1/100。Application选Transient fidelity，瞬态保真度模式会在接触状态变化时自动缩小时间增量，让增量结束点恰好落在接触状态改变的时刻。如果算的是准静态问题比如缓慢加载到额定载荷，Application选Quasi-static，时间增量控制逻辑会偏向收敛而非精度。

圆柱滚子轴承的接触定义比球轴承麻烦，滚子跟滚道是线接触，网格稍微粗一点就是点接触，接触力忽大忽小。Interaction模块用Surface-to-Surface接触，主面选滚道面，从面选滚子外表面，滚道网格比滚子网格粗一些，主从面刚度分配合理。法向用Hard Contact，切向用Penalty摩擦，摩擦系数取0.05到0.1，滚子跟滚道之间的摩擦小，取值大了滚子转不起来。接触对的数量是滚子数乘以2，内圈和保持架另算。初始分析步里加Adjust Tolerance自动消除初始穿透，几何模型建的时候滚子和滚道之间如果已经过盈或者有间隙，这一步能调平。

内圈和外圈用Kinematic Coupling耦合到各自轴心的参考点，滚子如果建了保持架，保持架也耦合到中心参考点，滚子本身不耦合，让接触力驱动它自转和公转。参考点上施加边界条件和载荷。外圈参考点约束六个自由度，模拟轴承座固定。内圈参考点约束径向平动和三个转动，释放轴向平动和绕轴转动，轴向平动上加径向载荷，绕轴转动上加转速。转速的施加方式有两种，一种是在Predefined Field里给初始角速度，另一种是在边界条件里给转动自由度指定幅值曲线。初始角速度用INITIAL CONDITIONS, TYPE=ROTATING VELOCITY，Load模块里Predefined Field → Create → Mechanical → Velocity → Rotating，Step选Initial，轴线和角速度值填进去。如果转速是随时间变化的，用边界条件加Amplitude曲线，在分析步里逐步加速到目标转速。

隐式动力学做轴承模型，收敛困难是常态。滚子数量多的时候，同时有十几个接触对处于开合边缘，牛顿迭代的雅可比矩阵容易奇异。Automatic Stabilization打开，初始比例因子取1e-4到1e-3， dissipated energy fraction保持默认的2e-4，让求解器在收敛困难的时候施加少量阻尼，算完之后检查ALLSD和ALLIE的比值，超过5%说明人工阻尼介入过多，结果不可信。Line Search打开，牛顿法在每次迭代后沿搜索方向走一步，放宽步长加速收敛。Maximum Increments默认1000不够用，改成10000，让求解器有更多尝试机会。

有一次做圆柱滚子轴承的启动过程，内圈从0加速到3000rpm，用了五个滚子先试。周五晚上提交，机房里空调关了，窗户开着一条缝，外面的蝉叫得很响，显示器亮度调到最低省电。跑了三个小时，滚子刚脱离静止状态就崩了，msg里全是“contact status changed”和“excessive rotation”。把外圈约束从六个自由度改成只约束径向和轴向、释放切向，让外圈可以微小转动，接触力分布均匀了很多。滚子网格从1mm加密到0.5mm，接触区域至少三层单元参与，穿透量降下来了。时间增量从1e-5调到1e-6，收敛了，但计算时间翻了十倍，跑了两天。后来又试着把滚子数量从五个减到三个，验证阶段用，确认载荷步设置没问题之后再加滚子，节省调试时间。

后处理关注三个量。滚子接触力用CPRESS在History Output里输出，每个滚子的接触力随时间的变化能看出载荷分布是否均匀，理想状态下同一时刻各滚子接触力应该按正弦规律分布。内圈和保持架的转速用UR3或VR3输出，看内圈加速过程中保持架有没有跟上，如果保持架转速明显滞后，可能是滚子跟滚道之间有打滑。能量输出ALLIE、ALLKE、ALLFD看数值合理性，ALLKE在稳定运转阶段应该接近零，如果一直震荡说明时间增量还是太大。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删