请问哪位做过转子-滑动轴承系统的扭转振动

health

本人目前采用每个节点具有5个自由度的转子有限元模型，采用有限元方法，将转子-轮盘系统划分节点，形成单元的刚度矩阵、质量矩阵和陀螺矩阵，然后将单元的矩阵叠加形成系统的整体刚度矩阵、陀螺矩阵和质量矩阵，再根据刚度矩阵和质量矩阵形成瑞利阻尼矩阵，对一根转子-轮盘-滑动轴承系统的运动微分方程进行数值仿真，激励为轮盘的不平衡质量，在两个油膜支承位置，考虑了非线性油膜力的作用，通过newmark方法求解了系统的动态响应，绘制时间历程图和频谱图，但是扭转方程中并没有加激励，扭转方程应该是有阻尼自由振动方程，计算的结果应该能稳定到一定的值，但是输出部分节点的扭转角位移时，发现扭转角位移并不是稳定在一个值，而是有10的负数12此方的变化，即小数点后面最后的一位数始终在逐渐的加1，我想了很久也没有想出原因来。扭转振动的时间历程图和频谱图