不可压缩、层流状态下动载荷作用无量纲雷诺方程为：

为润滑油粘度；h为油膜厚度； 为轴颈转速； 为从x周正向量起的油膜角； 为轴向坐标； 为径宽比。

油膜厚度为

其中：除ω外，其它量均为无量纲量H，x，y的量纲比为轴承半径间隙C， ， 的量纲比为 。

将黏度视为常数，Reynolds方程等价泛函L(p)为：

由变分原理知，Reynolds方程的解p应该满足

即

式中 为第i个节点压力，e为任意四边形单元，E为单元总数。

将所有单元都转换成局部坐标系下的正方形标准单元。则可用插值函数表达各量为：

式中， ， ，L为整个解区域的节点总数； 为第e个单元内第n个节点的插值函数；ξ、η为局部坐标，h_i、φ_i、λ_i分别为第i个节点的相应值。上标e表示对应单元内节点的值。

记： ， ，

逐单元积分后，对流量矩阵进行组装，得到对应的总刚矩阵K和载荷矩阵B，则压力P：

图1为利用有限元方法，求解的静态二维油膜力。其参数如表1所示，所采用的是40 × 10的九节点四边形单元。

五瓦可倾瓦轴承示意图如图2所示，可倾瓦轴承每块瓦块可以绕瓦块支点转动，各瓦相互独立，各瓦能够根据受力情况调整姿态，这也是可倾瓦轴承的优势所在。图3为轴承坐标系及其瓦块展开图， 表示油膜正压力区及油膜破裂区， 为油膜域 的边界，及 的交界线，轴颈中心为 ，轴承中心为 ，润滑域在动载情况下具有不确定性。

在计算单瓦油膜力时需要考虑瓦块的摆动，要求解瓦块压力，需要先确定瓦块的姿态和轴心坐标的相对位置。也就是求偏位角 值。此处为了简化计算，利用相对运动，把寻找β值转换为改变轴心的相对位置，试探性的微小改变轴心的x坐标，并在当前x坐标小计算得到压强值对枢轴的力矩

如果计算得到的力矩M不为零，则继续改变轴颈中心的x坐标，并计算相应的压强值，利用方程(5)计算力矩，直到M为零，其中 为瓦块第j个单元所对应的角度。同时还要考虑轴颈位置的改变是有极限位置的，无量纲化后，即是 。此时，按照当前x坐标下计算得到的压强值就是第i块瓦块的压强分布。求得相对位置后油膜力的计算方法与固定瓦相。

表1. 轴承参数

图1. 油膜压力分布

图2. 可倾瓦轴承结构示意图

图3. 轴承坐标系及其瓦块展开图

根据方程(4)可计算得到油膜压力 ，再通过积分得到油膜在x, y轴方向的合力 ，假设轴颈在运动过程中不与轴瓦接触，则轴颈在运动过程中所受到的力有轴瓦的支反力和不平衡质量引起的离心力。根据牛顿运动定律，轴颈的运动方程无量纲形式可以写为如下形式：

上式中，各量都是无量纲量，a的量纲比为 ；b的量纲比为 ；C为轴承的无量纲阻尼；K为轴承的无量纲刚度； ， 为第i块瓦在X和Y方向的分量，量纲比为 ；右边最后一项为存在质量偏心时的不平衡离心力在X和Y方向的分量； ， 量纲比为 。

由于考虑了轴的惯性力和油膜力 ， 的非线性，所以这种轨迹计算方法属于轴心轨迹的非线性分析，其计算步骤如下：

1) 设定初始的轴颈状态， ，计算压力p，油膜力 ， 。但应考虑到轴心位置X，Y，存在 。

2) 将计算得到的 ， ，代入公式(6-1)，(6-2)，计算得到加速度 ， 。

3) 用欧拉法计算轴下一时刻的状态，即 时间后，轴心的速度和位置 和 。

4) 计算得到 和 在回到1)，反复计算1)~3)，直到循环收敛或者小范围波动。

为了研究了阶跃载荷作用下的轴心轨迹运动规律，本文取不平衡力的变化规律如公式(8)所示，

图4是在阶跃载荷的作用下，轴颈中心的运动轨迹，在开始阶段q_x为0，轴心在自重和油膜力的作用下达到一个平衡点，即图示的1号平衡点。在16π的时间开始不平衡力x方向上的分力突然增加，而y方向上的分力不变为0，此时轴心必须找到一个新的平衡点，即2号平衡点，经过一段时间后轴颈开始平衡于2号平衡点。图5是轴颈中心分别在X和Y方向上轴颈中心的位移，从图5能够直观的看出，两个方向上的位移都出现了两个波动过程，波动出现的原因子图4中已经分析，从该图中能够更加清晰的知道稳定后的轴颈处于一个十分平稳的状态。

为了验证结果正确性，参考文献 [ 5 ] 所得的固定瓦滑动轴承在阶跃载荷作用下的轴颈中心运动规律，为了保持验证的可靠性，除固定瓦和可倾瓦的区别外，其他参数如载荷，径宽比，轴质量，转速等设置均一致，图6和图7为文献 [ 5 ] 得到的轴心轨迹图和轴心位移图，从两幅图看出在阶跃载荷作用下固定瓦

表2. 可倾瓦轴承参数

图4. 阶跃载荷作用时的轴心轨迹

图5. X和Y方向的位移

轴承收敛于两个平衡点且当载荷发生变化的时候会出现短暂的波动，和有限差分法所得的轴颈中心运动规律完全一致。

改变不平衡力取 ； ，图8为轴颈的轴心轨迹图，该图可以直观的看出轴心的运动轨迹，从图中可知，轴心初始阶段偏离平衡位置较远，但是经过一段时间后，逐渐趋于稳定，且稳定后的成8字型。图9分别为轴颈轴心在x轴和y轴方向的位移波形图，从中可以分别看出轴心位移在x和y方向分量的变化关系。结合三图可见，轴在一开始会经历一段时间的振荡，此时处于不稳定运转时期，轴心轨迹偏差较大，大概在经历650个时间步长后稳定下来，此时轴处于稳定运转时期，轴心轨迹趋于一致。图10为各瓦块1在油膜力对称轴处的压力，也即最大压力处随时间的变化关系图(图中

图6. 阶跃载荷作用时的轴心轨迹

图7. X和Y方向的位移

图8. 轴心轨迹图

取稳定后的100个时间步长)，该图可以直观的看出各瓦块油膜力相互调整的关系及各瓦块油膜力随时间的变化关系，结果显示各瓦块最大油膜力分布在时间轴均呈比较规则的周期变化的关系，这与轴心位移

图9. x，y方向位移

图10. 瓦块1最大压强随时间变化关系

在稳定后呈周期变化关系相符。