外文翻译--流体动压轴承-挠性转子系统的非线性动态特性中文版.doc

流体动压轴承-挠性转子系统的非线性动态特性 (西安交通大学润滑理论及轴承研究所 中国西安 710049) 摘要：分析了带液压轴承座的挠性转子系统的非线性动力性。1所示的转子-轴承系统是一个典型的非线性动力系统。该系统由线性部分(柔性转轴)和局部非线性部分(径向轴承)组成。图2所示的2节点具有8自由度的Timoshenko梁轴单元模型，由于其可以计及转动惯量与剪切变形的影响，故更接近实际运行的转子系统。因此采用有限元方法建立如下形式的柔性转轴横向振动方程 （1） 式中分别为转轴的质量矩阵、陀螺矩阵、刚度矩阵、周期为 的外力向量(包括重力和不平衡力)和轴承施加于转轴的非线性力向量。对于具有P个节点的转轴，节点位移向量可表示为 式中 分别为第j个节点沿水平和铅垂方向的横向位移与弯曲转角。非线性力向量可表示为 式中Fxj，Fyj分别为轴承作用在轴第j个节点上的水平和铅锤方向的油膜力。由于轴承的非线性油膜力孤立地作用于转子的个别节点上，因此对具有m个轴承支承的转子系统，轴承力 具有如下局部性质 式中xsb∈R4m, Fsb(xsb,xsb)∈R4m可被写成： 图1 转子-轴承系统示意图 图2 转轴有限单元模型 为了简化符号，将式(1)中各元素重新排序且表示为如下分块形式 由于需要花很多的时间计算多自由度的转子系统，在维持系统响应准确的情况下，减少系统的自由度是非常重要的。由于系统是局部非线性的，仅非线性自由度受控于非线性方程，而线性自由度又依赖于非线性自由度，因此可对线性自由度进行减缩，以使该系统降阶。为了避免缩减自由度时，坐标转换给系统的非线性因素带来的数值误差，仅将线性自由度转换为模态坐标，而非线性自由度和决定系统动力特性的非线性力仍保留在物理空间中，使降阶后的系统仍具有局部非线性特征。 为了降低线性组合的自由度，将XS表示为列的线性组合 式中 因此，矩阵的保留弹性特征模态的列是ωk∈(0，ωcut)的无阻尼特征值问题(ks-ω2jms)ψj=0(j=1,…，nk)的质量正则解。矩阵的剩余柔性模态的列可表示为 其中对角矩阵是角频率小于或等于时的谱矩阵。因此 从（11）式开始，可被写成： 这样就有如下整个变换 在以上等式中，矩阵变换T=T1T2T2，运用式(13)，可得缩减的系统方程 通过缩减把n(n=nb +nc)阶方程组减缩为s(s=nb +nc)阶方程组，由式(11)- (I4)，可见转轴的不平衡力及非线性项的影响全部保留在缩减的方程组(14)中。考虑到圆盘的不平衡力的影响，可得动力系统的运动方程 有 式中分md，Gd, Kd, Fdex别为圆盘的质量矩阵、阻尼矩阵、刚度矩阵、不平衡力向量。 当引入状态变量后，其相应的系统方程在状态空间中为 2.非线性力的计算及其流体动压轴承的Jacobians 矩阵 对于实际轴承，不具有解析形式的油膜力，而在计算系统非线性响应时，每一时刻动力积分均需要非线性力及其Jacobians矩阵的求解。非线性油膜力及Jacobians矩阵协调一致的精度不仅影响到求周期解的PNF是否收敛，而且对周期解的稳定性及分岔的分析有着极其重要的影响。同时任一时刻油膜力的Jacobians矩阵的准确性又影响着判断周期解稳定性的Floquet乘法的求解。基于以上问题，运用有限元法求解具有变分不等式形式的流体润滑Reynolds边值条件(油膜区域下游边值条件)问题。将油膜力视为某时刻轴颈中心位移及速度的函数， 由此可以得到一组微分方程，根据该方程组的特点，在求出油膜力的同时，可很快求得Jacobians矩阵。 对于有限长轴承流体润滑的Reynolds边值问题: 式中 — —油膜压力函数(表压) — — 润滑油动力粘度 ——轴承长径比的倒数 — — 油膜厚度 — — 从Y轴负方向到油膜位置的角度 — — 偏位角 — — 偏位线与轴承中心连线至油膜位置的角度，如图3和图4所示 式(18)可等价于如下离散的椭圆型变分不等式[18] 式中Φ(p,q)= 是H(Ω)× H(Ω) 上的强制、对称，椭圆双线性泛函,ψ(q)=,k={p∈H(Ω)，p≥0 in Ω}是线性泛函，H(Ω) )是Sobolev空间。是唯一的一层油，是的微分。是偏位线与轴承中心连线至油膜完整区和油膜破裂区交界线（此交界线是随位移及速度扰动变化的曲线)的夹角。由以上等式，油膜厚度h和变角可写成如下方程式： 图3 单块瓦计算坐标 图4 椭圆轴承及其计算坐标 由于非线性油膜力是轴颈中心动态位置的函数，所以油膜力可被写成如下等式： 油膜力Fx ,Fy过对动态压力分布分别积分后得到: 运用8节点等参有限元法可求出油膜区域各节点的压力分布Pi。压力