基于分离变量法的四瓦块可倾瓦轴承瓦块摆动特性分析.docxVIP

基于分离变量法的四瓦块可倾瓦轴承瓦块摆动特性分析 随着旋转轴承系统的大型化和高速化，旋转轴承的旋转元件的稳定性问题日益突出。在实际设计中，可倾斜辊轴承载的旋转系统的油膜失稳偶尔发生。可倾斜辊轴的滑动轴的油膜力模型是分析可倾斜辊轴的非线性动态性能的基础和关键。因此，有必要建立可靠的可倾斜辊轴的线性油膜力模型。目前，可倾斜辊轴的线性油膜力解法主要包括三种类型：直接法、数据库法和非线性油膜力公式。虽然直接的解算方法具有很高的精度，但计算效率低，基于数据库法的效率低。虽然建模过程中存在复杂且不相关的问题，但ola采用短圆柱形瓦臂轴承模型，进一步导出了可倾斜辊轴的非线性油膜力模型。对该模型与有限差分法的计算结果进行了比较和分析，但未提供瓦块的样本角公式。王立平等人基于瓦块和轴颈之间的运动耦合关系，建立了瓦块的局部坐标系统。通过变量变换，将瓦块和轴颈运动的整体位移向量转换为局部动静态下的轴颈部位移向量。利用固定瓦轴动态特性方法获得的局部动态坐标系下的油膜力。然后，局部坐标系下的油膜力应精确地转换为全局广义坐标系下的表达形式。 分析计算可倾瓦滑动轴承油膜动力特性时,给定轴颈在可倾瓦轴承中的扰动位移和扰动速度,如果能够快速获得瓦块摆角,从而快速获得油膜厚度分布规律,高效求解可倾瓦轴承油膜动力特性,对于可倾瓦滑动轴承动力特性的研究,以及可倾瓦轴承支承转子系统非线性动力学特性分析来说是具有重要意义的. 本文推导可倾瓦轴承瓦块油膜特性几何关系,并基于分离变量法建立有限长可倾瓦块油膜压力近似解析表达式,根据静平衡时瓦块受油膜压力力矩平衡,推导可倾瓦块摆角与可倾瓦轴承几何参数、间隙比、及原始偏心率间关系式,最后研究摆角随间隙比和原始偏心率变化规律. 1 可倾瓦块油膜特征角与油膜原始几何角的关系 图1(a)为一静止不动的轴颈和其支撑可倾瓦瓦块,轴颈半径为R,可倾瓦块与轴颈之间半径间隙为C,瓦块轴向长度为L,图中Φ1、Φ2、Φ3为瓦块原始几何角,确定了瓦块的结构参数.可倾瓦块与轴颈之间充满了粘度为μ的润滑介质.由于轴颈是没有运动的,瓦块圆弧中心O、轴颈中心Oj以及瓦块支点Op处于同一直线上,轴颈中心和瓦块中心间距离定义为原始偏心距e. 轴颈以角速度Ω旋转,假定轴颈旋转过程中,轴颈中心Oj不动,根据流体润滑理论知,可倾瓦块必定绕支点转动一定角度φ,如图1(b)所示. 这时,可倾瓦轴承瓦块油膜特性就转化为轴承中心O′、轴颈中心为Oj、偏心距为E、偏位角为α的有限弧度圆瓦轴承油膜特性,进而利用圆瓦轴承油膜相关理论分析可倾瓦块特性.图1(b)中θ为油膜周向坐标,起始位置为名义最大油膜厚度处.θ1、θ2、θ3定义为油膜特征角,分别为瓦块出油边边界,瓦块支点角,瓦块进油边边界. 由图1(b)可知,轴颈开始旋转后,油膜特征角与瓦块原始几何角有如下关系: {θ1={π,whileπ+α≤Φ1;Φ1-α,whileπ+αΦ1.θ2=Φ2-α.θ3=Φ3-α.(1) 图1(c)为图1(b)相关几何关系放大图,根据三角关系可得 E2=(R-e)2+R2-2R(R-e)cosφ, (2) A2=2R2(1-cosφ)=e2+E2-2Eecosβ. (3) 将式(2)带入式(3)整理后可得 β=arccose-R+Rcosφ√(R-e)2+R2-2R(R-e)cosφ.(4) 定义ε0=e/C为原始偏心率,ψ=C/R为半径间隙比.对式(4)进行无量纲化,得 β=arccosε0ψ-1+cosφ√(1-ε0ψ)2+1-2(1-ε0ψ)cosφ. 又由图1(c)可得知瓦块摆动后实际偏位角表达式为 α=arccosε0ψ-1+cosφ√(1-ε0ψ)2+1-2(1-ε0ψ)cosφ-φ,(5) 而瓦块摆动后实际偏心率为 ε=EC=√(1-ε0ψ)2+1-2(1-ε0ψ)cosφψ.(6) 2 量法求解压力轴向分布 本文研究瓦块摆动达到平衡后的摆角,即轴颈在轴承中无涡动,油膜只有动压效应,没有挤压效应.因此忽略流体润滑雷诺方程中所有动态效应项,在此条件下适合图1(b)中油膜润滑Reynolds方程形式如下: ??θ(h3?p?θ)+R2??z(h3?p?z)=-6μR2CεΩsinθ.(7) 式中:p为瓦块与轴颈之间油膜的压力;z为油膜轴向坐标,其零点在瓦块中分面处;h为油膜厚度,h=C(1+εcosθ). 轴向和周向边界条件为 p(θ1,z)=p(θ3,z)=p(θ,L/2)=p(θ,-L/2)=0. 利用分离变量法对雷诺方程进行求解,以获得压力近似解析解. 假设方程(7)解的形式如下: p(θ,z)=Μ(θ)Ν(z).(8) 其中,M和N分别为变量θ和z的二次可微函数,将式(8)带入式(7)整理后可得 Ν??θ(h3?Μ?θ)/h3+ΜR2?2Ν?z2=-6μ