Reynolds方程的求解.pdfVIP

第6 章 滑动轴承中油膜力及其Jacobi矩阵的确定 在非线性转子动力系统中，滑动轴承中的油膜力是一主要的非线性因素。在 多自由度的转子动力系统中，非线性因素只表现在某些局部自由度上，但它却影 响着整个动力系统的特性，即它的准确与否决定着动力系统特性计算的正确性。 对于轴承油膜力的研究已有很长的历史，在一些特殊条件下：如长轴承、短 轴承等，可以采用解析方法确定其表达式，而对于有限长轴承，其油膜力的求解， 需要借助数值方法确定，常用的数值方法有有限差分法、有限元方法、边界元方 法等。 为了精确的求解油膜力及其相应的 Jacobi 矩阵，本文采用变分不等方程描 述具有 Reynolds 边界条件的流体润滑问题，对于这类连续的变分不等方程，应 用有限元方法将其离散化，从而形成有限维的代数不等方程组。 任一时刻油膜力的 Jacobi 矩阵的准确性，影响着判断周期解稳定性的 Floquet乘子的求解。本章中，将油膜力视为此时刻轴颈位移及速度的函数，然 后采用小扰动方法，将其展开，并只保留线性项，通过这种处理可以得到一组微 分方程，根据该方程组的特点，在求出油膜力的同时，可以很快求得其 Jacobi 矩阵。 2. 1 滑动轴承中油膜力及其 Jacobi 矩阵的数学模型(采用部分瓦块方式) 在下面的推导过程中，假设流体的温度、密度为常数，并且为牛顿流体。对 于第 i 块瓦建立坐标系：设瓦块 i 的中心位置为o (x , y ) ，轴颈中心位置为 b b b  c  o (x ,y ) 。间隙函数在允许误差o 的限度内可表示为： j j j R   j    h(,t) c e(t) cos  (t) 经过推导，得到下列形式的 Reynolds 方程： 1   3 P    3 P  h h h  h       6(  ) 12 （2.1.1） R 2    z  z  b j  t 在分析平衡点和周期解的稳定性时, 均需要求解油膜力的 Jacobi 矩阵，所 以，将任一时刻的压力分布，表示成轴颈及瓦块中心的位移和速度的函数，故有 系列方程。 第i 瓦块上的油膜对轴颈的作用力通过下列积分式求得： F  L / 2 dz a2 R d(P P ) cos x L / 2 a j a   1 F  L / 2 dz a2 R d(P P ) sin （2.1.2） y L / 2 a j a   1 其中P 为大气压，a a a 。 a 2 1 其Jacobi矩阵中各元素通过下列积分式求得： F F x L / 2 a2