CN112182808A - 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法 - Google Patents

闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法 Download PDF

Info

Publication number
CN112182808A
CN112182808A CN202011159296.3A CN202011159296A CN112182808A CN 112182808 A CN112182808 A CN 112182808A CN 202011159296 A CN202011159296 A CN 202011159296A CN 112182808 A CN112182808 A CN 112182808A
Authority
CN
China
Prior art keywords
bearing
static
dynamic
spiral groove
closed high
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011159296.3A
Other languages
English (en)
Other versions
CN112182808B (zh
Inventor
蒋书运
张少文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202011159296.3A priority Critical patent/CN112182808B/zh
Publication of CN112182808A publication Critical patent/CN112182808A/zh
Application granted granted Critical
Publication of CN112182808B publication Critical patent/CN112182808B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C32/00Bearings not otherwise provided for
    • F16C32/06Bearings not otherwise provided for with moving member supported by a fluid cushion formed, at least to a large extent, otherwise than by movement of the shaft, e.g. hydrostatic air-cushion bearings
    • F16C32/0629Bearings not otherwise provided for with moving member supported by a fluid cushion formed, at least to a large extent, otherwise than by movement of the shaft, e.g. hydrostatic air-cushion bearings supported by a liquid cushion, e.g. oil cushion
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16CSHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
    • F16C32/00Bearings not otherwise provided for
    • F16C32/06Bearings not otherwise provided for with moving member supported by a fluid cushion formed, at least to a large extent, otherwise than by movement of the shaft, e.g. hydrostatic air-cushion bearings
    • F16C32/0629Bearings not otherwise provided for with moving member supported by a fluid cushion formed, at least to a large extent, otherwise than by movement of the shaft, e.g. hydrostatic air-cushion bearings supported by a liquid cushion, e.g. oil cushion
    • F16C32/064Bearings not otherwise provided for with moving member supported by a fluid cushion formed, at least to a large extent, otherwise than by movement of the shaft, e.g. hydrostatic air-cushion bearings supported by a liquid cushion, e.g. oil cushion the liquid being supplied under pressure
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mechanical Engineering (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Sliding-Contact Bearings (AREA)

Abstract

本发明提供了一种闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其包括以下步骤:考虑空化效应、离心效应、紊流效应和热效应等多因素的耦合作用,建立闭式高速水润滑动压螺旋槽推力轴承润滑模型;闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算;闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,以获得尽可能大的承载力和轴向刚度,并始终保持正向流量。采用本发明提供的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,能够大幅提高该类轴承静动态性能设计精度,缩短研发周期,为此类轴承性能分析和工程设计提供有效途径。

Description

闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法
技术领域
本发明涉及滑动轴承性能设计技术领域,尤其是一种闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法。
背景技术
闭式水润滑动压螺旋槽推力轴承依靠主轴的转动自主地形成动压流体膜,既能承受双向轴向载荷,又能承受倾覆力矩,具有一系列非常突出的技术优势,如:结构简单,成本低,承载力大,摩擦功耗低,还具有自动泵送作用,可望作为高速精密机床主轴的一种新型支承轴承,具有广阔的工程应用前景。
高速工况下,水润滑动压螺旋槽推力轴承会出现一系列的现象。首先,水的初生空化数高,水流中的微气核在高速下更容易生长成空化泡,轴承往往处于气-液两相流润滑状态。空化现象是一个非常复杂的动力学过程,气相和液相在空化区域会交替出现,气-液界面上还存在质量、动量和能量传递。其次,水的运动粘度低,当轴承边缘速度超过60m/s时,雷诺数可达到2000,水流更趋于紊态流动;同时,流体离心惯性力凸显,从而严重影响螺旋槽的自动泵送作用。此外,水膜厚度非常薄,约为油膜厚度的1/10,随着轴承转速和尺寸的不断增大,水膜会受到强烈的剪切作用,摩擦功耗大,粘性生热效应已不容忽视。
现有的动压螺旋槽推力轴承润滑理论主要针对轻载小尺寸轴承的油(或气体)润滑而建立的,不能准确地描述高速水润滑动压螺旋槽推力轴承的动力学行为,其计算精度低,不适用于此类轴承静动态性能设计。
发明内容
本发明基于二相流体动力学理论,提出一种基于雷诺方程的异质二相流空化模型,并在此基础上,考虑流体离心效应、紊流效应和热效应等多种因素的耦合作用,建立一种适用的闭式高速水润滑动压螺旋槽推力轴承润滑模型,继而提出此类轴承静动态性能设计方法,旨在大幅提高轴承性能设计精度,降低设计风险。
本发明采用的技术方案如下:
一种闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,包括以下步骤:
S1:考虑空化效应、离心效应、紊流效应和热效应等多因素的耦合作用,建立闭式高速水润滑动压螺旋槽推力轴承润滑模型:建立液相静动态广义雷诺方程、液相能量方程、空泡体积分数输运方程、单个空泡运动方程,解得稳态压力和扰动压力,获得轴承的静动态性能参数;
s2:基于s1所建立的轴承润滑模型,开展闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算;
基于s2的分析结果,开展闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,以获得尽可能大的承载力和轴向刚度,并始终保持正向流量。
所述s1具体通过以下s11-s16六个步骤实现:
s11:坐标系的建立与转换:
所述推力轴承由两个相对的固定设置的推力瓦和一个顺时针旋转的推力盘组成;所述推力盘的上工作面、下工作面均刻有周期分布的螺旋槽,两个所述推力瓦的工作表面则完全光滑;
所建立的圆柱坐标系(r,θ,z)与以角速度ω旋转的推力盘固连,而所建立的惯性坐标系(r,Φ,z)与静止的推力瓦固连;此外,还建立笛卡尔坐标系(x,y,z);
运用边界拟合坐标系对轴承几何形状进行转换,物理平面(r,θ)和计算平面(ξ,η)之间的转换关系为:
Figure BDA0002741673940000021
式(1)中:β为所述螺旋槽的螺旋角;rs为密封半径;r为轴承的半径方向坐标;θ为圆周方向坐标;ξ为螺旋槽方向坐标;η为圆周方向坐标;
s12:液相静动态广义雷诺方程的建立:
液相静态广义雷诺方程为:
Figure BDA0002741673940000022
式(2)中:下标中“1”、“2”、“3”、“4”表示控制体积中的区域编号;
p0为稳态水膜压力;
Figure BDA0002741673940000023
Figure BDA0002741673940000024
为瞬时质量流量的稳态分量,其可表示为:
Figure BDA0002741673940000025
式(3)中:“A”、“B”、“C”、“D”、“E”代表流量系数;αl为液相体积分数;ρl为液相密度;t为时间;h0为轴承处于稳态时的水膜厚度;
液相动态广义雷诺方程为:
Figure BDA0002741673940000026
式(4)中:pγ为扰动水膜压力;
Figure BDA0002741673940000031
Figure BDA0002741673940000032
为瞬时质量流量的扰动分量,均为复数,其可表示为:
Figure BDA0002741673940000033
式(5)中:“A”、“B”、“C”、“D”、“E”代表流量系数;υ为小扰动频率;
Figure BDA0002741673940000034
为扰动膜厚,其可表示为:
Figure BDA0002741673940000035
式(6)中:对于推力盘的下表面,α取+1,对于推力盘的上表面,α取-1;
s13:液相能量方程的建立:
Figure BDA0002741673940000036
式(7)中:h为水膜厚度;
Figure BDA0002741673940000037
为通过膜厚的液相平均温度;Tp和Td分别为推力瓦和推力盘温度;cv为液相定容比热;μl为液相动力粘度;Reh为当地雷诺数;kθ和kr分别为周向和径向紊流系数;Mi为周向动量传递项;Ei为两相之间的能量传递项;kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数;
s14:空泡体积分数输运方程的建立:
Figure BDA0002741673940000041
式(8)中:αb为空泡体积分数;
Figure BDA0002741673940000042
Figure BDA0002741673940000043
分别为空泡相的周向和径向平均速度;Rv和Rc分别为空泡生长和溃灭的源项;
s15:单个空泡运动方程的建立:
FP+FD+FA=0 (9)
式(9)中:Fp为压差力;FD为阻力;FA为由空泡生长和溃灭引起的附加质量力。
s16:联立方程(2)、(4)、(7)~(9),设置边界条件,运用有限差分法离散方程组,采用超松弛迭代法求解方程(2)和(4),采用低松弛迭代法求解方程(7)和(8),解得稳态压力和扰动压力,从而获得轴承的静动态性能参数。
所述s2中,闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算,基于s1所建立的闭式高速水润滑动压螺旋槽推力轴承润滑模型和数值计算方法,系统地开展轴承静态性能和动态性能对设计变量的灵敏度分析,并对比研究泵入、泵出两类典型的闭式高速水润滑动压螺旋槽推力轴承静动态性能。
所述s2中,所述闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,基于s2建立的轴承静动态性能对设计变量的灵敏度分析结果和泵入、泵出两类典型轴承的静动态性能对比结果,以最大承载力和轴向刚度系数所对应的设计变量作为最佳值,以保持正向流量的结构作为合理的轴承选型。
轴承静态性能包括承载力、抗倾覆力矩、体积流量和摩擦功耗,轴承动态性能包括刚度和阻尼系数;所述设计变量包括倾斜角、转速、偏心率、半径比、螺旋角、径向槽长比和槽深。
本发明的有益效果如下:
本发明提供的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,提出一种基于雷诺方程的异质二相流空化模型,并在此基础上考虑流体离心效应、紊流效应和热效应等多种因素的耦合作用,建立一种适用的闭式高速水润滑动压螺旋槽推力轴承润滑模型,为闭式高速水润滑动压螺旋槽推力轴承静动态性能分析与设计,提供了行之有效的手段,同时能够大幅提高该类轴承性能设计的精度,降低设计风险、缩短研发周期。
本发明实现了系统地开展轴承静态(承载力、抗倾覆力矩、体积流量和摩擦功耗)和动态(刚度和阻尼系数)性能对设计变量(倾斜角、转速、偏心率、半径比、螺旋角、径向槽长比和槽深)的灵敏度分析,并可全面地对比研究泵入、泵出两类典型的闭式高速水润滑动压螺旋槽推力轴承静动态性能,能够以最大承载力和轴向刚度系数所对应的设计变量作为最佳值,以保持正向流量的结构作为合理的轴承选型,从而为此类轴承的选型和设计提供理论基础。
附图说明
图1是本发明所针对的推力轴承的几何结构与坐标系示意图。
图2是泵入式轴承槽型的上工作面、下工作面的结构与及坐标系示意图。
图3是泵出式轴承槽型的上工作面、下工作面的结构与及坐标系示意图。
图4是本发明基于边界拟合坐标系的坐标转换过程的示意图。
图5是本发明计算过程设定的控制体积内的流量平衡示意图。
图6是本发明轴承静动态性能具体计算流程图。
图7是本发明专利的设计方法的主要步骤和内容。
图8是轴承静态性能随设计变量的变化曲线-抗倾覆力矩随倾斜角
Figure BDA0002741673940000051
的变化曲线(偏心率ε=0.6时)。
图9是轴承静态性能随设计变量的变化曲线-体积流量随转速的变化曲线。
图10是轴承静态性能随设计变量的变化曲线-体积流量随偏心率的变化曲线。
图11是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随半径比的变化曲线。
图12是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随螺旋角的变化曲线。
图13是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随径向槽长比的变化曲线。
图14是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随槽深的变化曲线。
图15是轴承动态性能随设计变量的变化曲线-轴向刚度系数随倾斜角
Figure BDA0002741673940000052
的变化曲线。
图16是轴承动态性能随设计变量的变化曲线-轴向刚度系数随转速的变化曲线。
图17是轴承动态性能随设计变量的变化曲线-轴向刚度系数随偏心率的变化曲线。
图18是轴承动态性能随设计变量的变化曲线-轴向刚度系数随半径比的变化曲线。
图19是轴承动态性能随设计变量的变化曲线-轴向刚度系数随螺旋角的变化曲线。
图20是轴承动态性能随设计变量的变化曲线-轴向刚度系数随径向槽长比的变化曲线。
图21是轴承动态性能随设计变量的变化曲线-轴向刚度系数随槽深的变化曲线。
具体实施方式
以下结合附图说明本发明的具体实施方式。
如图7所示,给出了本实施例的高速水润滑动压螺旋槽推力轴承静动态性能设计方法的内容,主要包括以下步骤:
s1:考虑空化效应、离心效应、紊流效应和热效应等多因素的耦合作用,建立闭式高速水润滑动压螺旋槽推力轴承润滑模型:
闭式高速水润滑动压螺旋槽推力轴承的几何结构如图1所示,推力轴承由两个相对的固定设置的推力瓦和一个顺时针旋转的推力盘组成;推力盘的上工作面、下工作面均刻有周期分布的螺旋槽,而两个推力瓦的工作表面则完全光滑;
考虑空化效应、离心效应、紊流效应和热效应等多因素的耦合作用,建立闭式高速水润滑动压螺旋槽推力轴承润滑模型的过程可分为s11-s16六步:
S11:坐标系的建立与转换:
建立如图1-图3所示的坐标系,所建立的圆柱坐标系(r,θ,z)与以角速度ω旋转的推力盘固连,而所建立的惯性坐标系(r,Φ,z)与静止的推力瓦固连;此外,还建立笛卡尔坐标系(x,y,z),其中,x轴的正方向指向Φ=0,y轴的正方向指向Φ=π/2;
图中显示的轴承几何参数包括轴承间隙c、螺旋角β、槽深cg、槽宽bg、台宽br、外圆半径r0、内圆半径ri和密封半径rs;槽宽比κθ和径向槽长比κr是两个通常用来描述螺旋槽推力轴承的无量纲几何参数;ε为闭式推力轴承的偏心率。
比较图2和图3,泵入式轴承和泵出式轴承的主要区别在于:
1)螺旋角的方向相反;2)密封坝的位置不同。对于泵入式轴承而言,密封坝靠近轴承的内圆周;正常情况下,润滑水在螺旋槽的带动下,由轴承外侧流进轴承间隙,然后由轴承内侧流出。然而,对于泵出式轴承而言,密封坝靠近轴承的外圆周;润滑水在螺旋槽的带动下,会由轴承内侧流进轴承间隙,然后由轴承外侧流出。
如图4所示,运用边界拟合坐标系对轴承几何形状进行转换,物理平面(r,θ)和计算平面(ξ,η)之间的转换关系如式(1)所示:
Figure BDA0002741673940000061
式(1)中:β为所述螺旋槽的螺旋角;rs为密封半径;r为轴承的半径方向坐标;θ为圆周方向坐标,逆时针为正;ξ为螺旋槽方向坐标;η为圆周方向坐标,逆时针为正;
物理平面(r,θ)和计算平面(ξ,η)均随推力盘一起作旋转运动,而坐标系(r,Φ,z)则静止;
s12:液相静动态广义雷诺方程的建立:
考虑流体离心效应和界面效应,按式(2)列写液相静态广义雷诺方程,按式(4)列写液相动态广义雷诺方程:
液相静态广义雷诺方程为:
Figure BDA0002741673940000062
式(2)中:下标中“1”、“2”、“3”、“4”表示如图5所示的控制体积中的区域编号,图中控制体积中的箭头表示通过控制体积边界流量的正方向;
稳态水压膜力和扰动水膜压力
Figure BDA0002741673940000063
在节点pj,i处取值;
空泡体积分数αb、液相平均温度
Figure BDA00027416739400000610
液相平均速度(
Figure BDA0002741673940000064
Figure BDA0002741673940000065
)和空泡平均速度(
Figure BDA0002741673940000066
Figure BDA0002741673940000067
)均在节点Nj,i处取值;
p0为稳态水膜压力;
Figure BDA0002741673940000068
Figure BDA0002741673940000069
为瞬时质量流量的稳态分量,其可表示为:
Figure BDA0002741673940000071
式(3)中:αl为液相体积分数;ρl为液相密度;t为时间;h0为轴承处于稳态时的水膜厚度;
式(3)中
Figure BDA0002741673940000072
Figure BDA0002741673940000073
的各流量系数如下式(I.1)所示:
Figure BDA0002741673940000074
Figure BDA0002741673940000075
Figure BDA0002741673940000076
Figure BDA0002741673940000077
Figure BDA0002741673940000078
式(I.1)中:kθ和kr为紊流系数;α′、β′和γ′为惯性系数;Mi为周向动量传递项;μl为液相动力粘度;
Figure BDA0002741673940000079
J=rξ(rθη)-rη(rθξ);
液相动态广义雷诺方程为:
Figure BDA00027416739400000710
式(4)中:pγ为扰动水膜压力;
Figure BDA00027416739400000711
Figure BDA00027416739400000712
为瞬时质量流量的扰动分量,均为复数,其可表示为:
Figure BDA0002741673940000081
式(5)中:υ为小扰动频率;
Figure BDA0002741673940000082
Figure BDA0002741673940000083
为扰动膜厚,其可表示为:
Figure BDA0002741673940000084
式(6)中:对于推力盘的下表面,α取+1,对于推力盘的上表面,α取-1;
式(5)中
Figure BDA0002741673940000085
Figure BDA0002741673940000086
的各流量系数如下式(I.2)示:
Figure BDA0002741673940000087
Figure BDA0002741673940000088
Figure BDA0002741673940000089
Figure BDA00027416739400000810
Figure BDA00027416739400000811
式(I.2)中:kθ和kr为紊流系数;α′、β′和γ′为惯性系数;Mi为周向动量传递项;μl为液相动力粘度;
Figure BDA00027416739400000812
J=rξ(rθη)-rη(rθξ);
s13:液相能量方程的建立:
考虑紊流效应、气-液界面能量传递和液固间的热量传递,按式(7)列写液相能量方程如下:
Figure BDA0002741673940000091
式(7)中:h为水膜厚度;
Figure BDA0002741673940000092
为通过膜厚的液相平均温度;Tp和Td分别为推力瓦和推力盘温度;cv为液相定容比热;Reh为当地雷诺数;kθ和kr分别为周向和径向紊流系数;Mi为周向动量传递项;Ei为两相之间的能量传递项;kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数;
s14:空泡体积分数输运方程的建立:
采用由压力驱动的空化模型,按式(8)列写空泡体积分数输运方程如下:
Figure BDA0002741673940000093
式(8)中:αb为空泡体积分数;
Figure BDA0002741673940000095
Figure BDA0002741673940000096
分别为空泡相的周向和径向平均速度;Rv和Rc分别为空泡生长和溃灭的源项;
s15:单个空泡运动方程的建立:
考虑压差力、阻力和附加质量力之间的平衡,按式(9)列写单个空泡运动方程如下:
FP+FD+FA=0 (9)
式(9)中:Fp为压差力;FD为阻力;FA为由空泡生长和溃灭引起的附加质量力。
s16:联立方程(2)、(4)、(7)~(9),设置边界条件,运用有限差分法离散方程组,采用超松弛迭代法求解方程(2)和(4),采用低松弛迭代法求解方程(7)和(8),解得稳态压力和扰动压力,从而获得轴承的静动态性能参数,包括如式(10)-(13)所示意的轴承承载力、轴承液相沿径向的质量流量、轴承摩擦功耗、轴承刚度和阻尼系数组成的复阻抗,计算流程如图6所示;
其中,轴承承载力表示如下:
Figure BDA0002741673940000094
式(10)中:下标“lower”和“upper”分别表示推力盘的下表面和上表面;Wz为轴向承载力;Mx、My分别为绕x轴和y轴的抗倾覆力矩。
轴承液相沿径向的质量流量表示如下:
Figure BDA0002741673940000101
式(11)中:r0和ri分别为轴承的外圆半径和内圆半径;
轴承摩擦功耗表示如下:
Figure BDA0002741673940000102
将扰动压力
Figure BDA0002741673940000103
在推力盘的上、下表面积分,可得由推力轴承刚度和阻尼系数组成的复阻抗表示如下:
Figure BDA0002741673940000104
s2:闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算:
闭式高速水润滑动压螺旋槽推力轴承的几何和工况参数列于表1;
表1
Figure BDA0002741673940000105
在确定某一设计变量对轴承静动态性能的影响时,其他变量均保持不变。
基于s1所建立的轴承润滑模型和数值计算方法,系统地开展倾斜角、转速、偏心率、半径比、螺旋角、径向槽长比和槽深等设计变量对闭式高速水润滑动压螺旋槽推力轴承静动态性能的影响分析,计算结果如图8-图14、图15-图21所示。
计算结果反映出:流体离心效应对轴承流量具有十分显著的影响,在高转速和大偏心条件下,流体离心惯性力会导致泵入式轴承发生反流现象,进而引起轴承润滑性能的失效,而泵出式轴承却能始终保持正向流量;半径比、螺旋角、径向槽长比和槽深等设计变量都是影响轴承静动态性能的重要参数。
基于s2的分析结果,开展闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,以获得尽可能大的承载力和轴向刚度,并始终保持正向流量;
基于s2的计算结果,得出如下设计准则:
(1)在轴承选型方面:高速工况下,尽管泵出式轴承的摩擦功耗大于泵入式轴承,但是当半径比大于0.7时,泵出式轴承却能表现出更加优异的工作性能,如:具有较大的轴向承载力、更强的抗倾覆能力、更高的刚度和阻尼系数和始终保持正向的径向流量。
(2)在轴承结构参数方面:两种轴承最大轴向承载力所对应的结构参数相同,最佳的螺旋角、径向槽长比和槽台膜厚比分别约为25°、0.7和5.29;但两者最大轴向刚度系数所对应的结构参数却不相同,泵入式轴承最佳的螺旋角、径向槽长比和槽台区膜厚比分别为20°,0.7~0.75和4.21~4.93,而泵出式轴承最佳的螺旋角、径向槽长比和槽台区膜厚比分别为25°,0.65和4.93~5.29。

Claims (5)

1.一种闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其特征在于,包括以下步骤:
S1:考虑空化效应、离心效应、紊流效应和热效应等多因素的耦合作用,建立闭式高速水润滑动压螺旋槽推力轴承润滑模型:建立液相静动态广义雷诺方程、液相能量方程、空泡体积分数输运方程、单个空泡运动方程,解得稳态压力和扰动压力,获得轴承的静动态性能参数;
s2:基于s1所建立的轴承润滑模型,开展闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算;
基于s2的分析结果,开展闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,以获得尽可能大的承载力和轴向刚度,并始终保持正向流量。
2.根据权利要求1所述的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其特征在于,所述s1具体通过以下s11-s16六个步骤实现:
s11:坐标系的建立与转换:
所述推力轴承由两个相对的固定设置的推力瓦和一个顺时针旋转的推力盘组成;所述推力盘的上工作面、下工作面均刻有周期分布的螺旋槽,两个所述推力瓦的工作表面则完全光滑;
所建立的圆柱坐标系(r,θ,z)与以角速度ω旋转的推力盘固连,而所建立的惯性坐标系(r,Φ,z)与静止的推力瓦固连;此外,还建立笛卡尔坐标系(x,y,z);
运用边界拟合坐标系对轴承几何形状进行转换,物理平面(r,θ)和计算平面(ξ,η)之间的转换关系为:
Figure FDA0002741673930000011
式(1)中:β为所述螺旋槽的螺旋角;rs为密封半径;r为轴承的半径方向坐标;θ为圆周方向坐标;ξ为螺旋槽方向坐标;η为圆周方向坐标;
s12:液相静动态广义雷诺方程的建立:
液相静态广义雷诺方程为:
Figure FDA0002741673930000012
式(2)中:下标中“1”、“2”、“3”、“4”表示控制体积中的区域编号;
p0为稳态水膜压力;
Figure FDA0002741673930000013
Figure FDA0002741673930000014
为瞬时质量流量的稳态分量,其可表示为:
Figure FDA0002741673930000021
式(3)中:“A”、“B”、“C”、“D”、“E”代表流量系数;
αl为液相体积分数;ρl为液相密度;t为时间;h0为轴承处于稳态时的水膜厚度;
液相动态广义雷诺方程为:
Figure FDA0002741673930000022
式(4)中:pγ为扰动水膜压力;
Figure FDA0002741673930000023
Figure FDA0002741673930000024
为瞬时质量流量的扰动分量,均为复数,其可表示为:
Figure FDA0002741673930000025
式(5)中:“A”、“B”、“C”、“D”、“E”代表流量系数;υ为小扰动频率;
Figure FDA0002741673930000026
为扰动膜厚,其可表示为:
Figure FDA0002741673930000027
式(6)中:对于推力盘的下表面,α取+1,对于推力盘的上表面,α取-1;
s13:液相能量方程的建立:
Figure FDA0002741673930000031
式(7)中:h为水膜厚度;
Figure FDA0002741673930000032
为通过膜厚的液相平均温度;Tp和Td分别为推力瓦和推力盘温度;cv为液相定容比热;μl为液相动力粘度;Reh为当地雷诺数;kθ和kr分别为周向和径向紊流系数;Mi为周向动量传递项;Ei为两相之间的能量传递项;kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数;
s14:空泡体积分数输运方程的建立:
Figure FDA0002741673930000033
式(8)中:αb为空泡体积分数;
Figure FDA0002741673930000034
Figure FDA0002741673930000035
分别为空泡相的周向和径向平均速度;Rv和Rc分别为空泡生长和溃灭的源项;
s15:单个空泡运动方程的建立:
FP+FD+FA=0 (9)
式(9)中:Fp为压差力;FD为阻力;FA为由空泡生长和溃灭引起的附加质量力。
s16:联立方程(2)、(4)、(7)~(9),设置边界条件,运用有限差分法离散方程组,采用超松弛迭代法求解方程(2)和(4),采用低松弛迭代法求解方程(7)和(8),解得稳态压力和扰动压力,从而获得轴承的静动态性能参数。
3.根据权利要求2所述的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其特征在于,所述s2中,闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算,基于s1所建立的闭式高速水润滑动压螺旋槽推力轴承润滑模型和数值计算方法,系统地开展轴承静态性能和动态性能对设计变量的灵敏度分析,并对比研究泵入、泵出两类典型的闭式高速水润滑动压螺旋槽推力轴承静动态性能。
4.根据权利要求3所述的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其特征在于,所述s2中,所述闭式高速水润滑动压螺旋槽推力轴承静动态性能设计,基于s2建立的轴承静动态性能对设计变量的灵敏度分析结果和泵入、泵出两类典型轴承的静动态性能对比结果,以最大承载力和轴向刚度系数所对应的设计变量作为最佳值,以保持正向流量的结构作为合理的轴承选型。
5.根据权利要求3所述的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法,其特征在于,轴承静态性能包括承载力、抗倾覆力矩、体积流量和摩擦功耗,轴承动态性能包括刚度和阻尼系数;所述设计变量包括倾斜角、转速、偏心率、半径比、螺旋角、径向槽长比和槽深。
CN202011159296.3A 2020-10-26 2020-10-26 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法 Active CN112182808B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011159296.3A CN112182808B (zh) 2020-10-26 2020-10-26 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011159296.3A CN112182808B (zh) 2020-10-26 2020-10-26 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法

Publications (2)

Publication Number Publication Date
CN112182808A true CN112182808A (zh) 2021-01-05
CN112182808B CN112182808B (zh) 2024-05-03

Family

ID=73923282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011159296.3A Active CN112182808B (zh) 2020-10-26 2020-10-26 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法

Country Status (1)

Country Link
CN (1) CN112182808B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113027921A (zh) * 2021-02-09 2021-06-25 太原重工股份有限公司 获取静动压轴承油膜压力分布的方法和装置
CN113779712A (zh) * 2021-08-16 2021-12-10 东南大学 高速水润滑径向动压织构轴承静特性建模与设计方法
CN115855356A (zh) * 2022-11-24 2023-03-28 中国核动力研究设计院 一种用于支承圆柱挤压膜力测量的实验装置及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008038989A (ja) * 2006-08-03 2008-02-21 Tokai Univ 溝付き動圧スラスト気体軸受およびその製造方法
CN102269654A (zh) * 2011-05-10 2011-12-07 重庆大学 水润滑轴承及传动系统综合性能实验平台
CN107255636A (zh) * 2017-04-28 2017-10-17 东南大学 一种水润滑动压螺旋槽推力轴承的空化观测装置及方法
CN109359318A (zh) * 2018-08-16 2019-02-19 东南大学 空气轴承电主轴5自由度耦合刚性转子系统动态设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008038989A (ja) * 2006-08-03 2008-02-21 Tokai Univ 溝付き動圧スラスト気体軸受およびその製造方法
CN102269654A (zh) * 2011-05-10 2011-12-07 重庆大学 水润滑轴承及传动系统综合性能实验平台
CN107255636A (zh) * 2017-04-28 2017-10-17 东南大学 一种水润滑动压螺旋槽推力轴承的空化观测装置及方法
CN109359318A (zh) * 2018-08-16 2019-02-19 东南大学 空气轴承电主轴5自由度耦合刚性转子系统动态设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIN XIAOHUI 等: "Thermohydrodynamic model of cavitating flow and dynamic characteristic calculation for high-speed water-lubricated pump-out spiral groove bearing", TRIBOLOGY TRANSACTIONS, 19 May 2020 (2020-05-19) *
蒋书运 等: "高速精密水润滑电主轴关键技术研究进展", 机械设计与制造工程, vol. 45, no. 5, 31 May 2016 (2016-05-31) *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113027921A (zh) * 2021-02-09 2021-06-25 太原重工股份有限公司 获取静动压轴承油膜压力分布的方法和装置
CN113779712A (zh) * 2021-08-16 2021-12-10 东南大学 高速水润滑径向动压织构轴承静特性建模与设计方法
CN113779712B (zh) * 2021-08-16 2022-11-04 东南大学 高速水润滑径向动压织构轴承静特性建模与设计方法
CN115855356A (zh) * 2022-11-24 2023-03-28 中国核动力研究设计院 一种用于支承圆柱挤压膜力测量的实验装置及方法
CN115855356B (zh) * 2022-11-24 2024-04-02 中国核动力研究设计院 一种用于支承圆柱挤压膜力测量的实验装置及方法

Also Published As

Publication number Publication date
CN112182808B (zh) 2024-05-03

Similar Documents

Publication Publication Date Title
CN112182808B (zh) 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法
CN112580158B (zh) 高速水润滑动压螺旋槽径向轴承静动态性能设计方法
Guo et al. Application of CFD analysis for rotating machinery—Part I: Hydrodynamic, hydrostatic bearings and squeeze film damper
Srinivasan Analysis of static and dynamic load on hydrostatic bearing with variable viscosity and pressure
Zhang et al. Dynamic model and load sharing performance of planetary gear system with journal bearing
CN106354987A (zh) 一种重载静压转台承载力与油垫温度场分布规律关系计算方法
Du et al. Effect of different turbulent lubrication models on the lubrication characteristics of water-lubricated rubber bearings at a high Reynolds number
Zuo et al. Comparative performance analysis of conical hydrostatic bearings compensated by variable slot and fixed slot
Nagai et al. Numerical and experimental analyses of dynamic characteristics for liquid annular seals with helical grooves in seal stator
Nicholson et al. Fluid dynamics of the slip boundary condition for isothermal rimming flow with moderate inertial effects
Li et al. Manufacturing error and misalignment effect on the transient lubrication behavior of dynamically loaded journal bearing with micro-groove
Temis et al. Influence of centrifugal forces on oil flow in journal bearing of planetary gear
Zhang et al. Comparative study on static characteristics of double-pad inwardly and outwardly pumping water-lubricated spiral-groove thrust bearings
Voloshina et al. Improvement of the angular arrangement of distribution system windows when designing planetary hydraulic machines
Wei et al. Numerical calculation of fluid film force on journal bearings based on a biconjugate gradient-stabilized algorithm
Martin et al. Theoretical studies of a continuously adjustable hydrodynamic fluid film bearing
Yang et al. Analysis and Optimization of an Internal Feedback Hydrostatic Turntable Oil Pad Power Consumption Based on Finite Difference Method
Zuo et al. Performance analysis of multirecess angled-surface slot-compensated conical hydrostatic bearing
CN113779712B (zh) 高速水润滑径向动压织构轴承静特性建模与设计方法
Kim et al. Analyses of thrust bearings in scroll compressors considering keyways
Jianhua et al. Effect of film thickness on load-carrying property of seawater dynamic-hydrostatic hybrid thrust bearing
Feng et al. Numerical design method for water-lubricated hybrid sliding bearings
Soulas et al. Performance of damaged hydrostatic bearings: predictions versus experiments
Lu et al. Modeling and analysis for the transmission characteristics of a proposed dual-drive hydrostatic lead screw system
Martin Extended expansion of the Reynolds equation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant