CN112182808A - 闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法 - Google Patents
闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法 Download PDFInfo
- 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
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 70
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 57
- 238000013461 design Methods 0.000 title claims abstract description 53
- 238000005461 lubrication Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000000694 effects Effects 0.000 claims abstract description 31
- 230000001808 coupling effect Effects 0.000 claims abstract description 7
- 239000007791 liquid phase Substances 0.000 claims description 46
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000012546 transfer Methods 0.000 claims description 14
- 230000001050 lubricating effect Effects 0.000 claims description 9
- 239000012071 phase Substances 0.000 claims description 7
- 239000011800 void material Substances 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 6
- 238000007789 sealing Methods 0.000 claims description 6
- 238000010206 sensitivity analysis Methods 0.000 claims description 5
- 210000003934 vacuole Anatomy 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000012827 research and development Methods 0.000 abstract description 2
- 239000012530 fluid Substances 0.000 description 8
- 239000007788 liquid Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005514 two-phase flow Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16C—SHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
- F16C32/00—Bearings not otherwise provided for
- F16C32/06—Bearings 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/0629—Bearings 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16C—SHAFTS; FLEXIBLE SHAFTS; ELEMENTS OR CRANKSHAFT MECHANISMS; ROTARY BODIES OTHER THAN GEARING ELEMENTS; BEARINGS
- F16C32/00—Bearings not otherwise provided for
- F16C32/06—Bearings 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/0629—Bearings 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/064—Bearings 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force 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，θ)和计算平面(ξ，η)之间的转换关系为：
式(1)中：β为所述螺旋槽的螺旋角；rs为密封半径；r为轴承的半径方向坐标；θ为圆周方向坐标；ξ为螺旋槽方向坐标；η为圆周方向坐标；
s12：液相静动态广义雷诺方程的建立：
液相静态广义雷诺方程为：
式(2)中：下标中“1”、“2”、“3”、“4”表示控制体积中的区域编号；
p0为稳态水膜压力；
和
为瞬时质量流量的稳态分量，其可表示为：
式(3)中：“A”、“B”、“C”、“D”、“E”代表流量系数；αl为液相体积分数；ρl为液相密度；t为时间；h0为轴承处于稳态时的水膜厚度；
液相动态广义雷诺方程为：
式(4)中：pγ为扰动水膜压力；
和
为瞬时质量流量的扰动分量，均为复数，其可表示为：
式(5)中：“A”、“B”、“C”、“D”、“E”代表流量系数；υ为小扰动频率；
为扰动膜厚，其可表示为：
式(6)中：对于推力盘的下表面，α取+1，对于推力盘的上表面，α取-1；
s13：液相能量方程的建立：
式(7)中：h为水膜厚度；
为通过膜厚的液相平均温度；Tp和Td分别为推力瓦和推力盘温度；cv为液相定容比热；μl为液相动力粘度；Reh为当地雷诺数；kθ和kr分别为周向和径向紊流系数；Mi为周向动量传递项；Ei为两相之间的能量传递项；kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数；
s14：空泡体积分数输运方程的建立：
式(8)中：αb为空泡体积分数；
和
分别为空泡相的周向和径向平均速度；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是轴承静态性能随设计变量的变化曲线-抗倾覆力矩随倾斜角
的变化曲线(偏心率ε＝0.6时)。
图9是轴承静态性能随设计变量的变化曲线-体积流量随转速的变化曲线。
图10是轴承静态性能随设计变量的变化曲线-体积流量随偏心率的变化曲线。
图11是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随半径比的变化曲线。
图12是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随螺旋角的变化曲线。
图13是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随径向槽长比的变化曲线。
图14是轴承静态性能随设计变量的变化曲线-轴向承载力和摩擦功耗随槽深的变化曲线。
图15是轴承动态性能随设计变量的变化曲线-轴向刚度系数随倾斜角
的变化曲线。
图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)所示：
式(1)中：β为所述螺旋槽的螺旋角；rs为密封半径；r为轴承的半径方向坐标；θ为圆周方向坐标，逆时针为正；ξ为螺旋槽方向坐标；η为圆周方向坐标，逆时针为正；
物理平面(r，θ)和计算平面(ξ，η)均随推力盘一起作旋转运动，而坐标系(r，Φ，z)则静止；
s12：液相静动态广义雷诺方程的建立：
考虑流体离心效应和界面效应，按式(2)列写液相静态广义雷诺方程，按式(4)列写液相动态广义雷诺方程：
液相静态广义雷诺方程为：
式(2)中：下标中“1”、“2”、“3”、“4”表示如图5所示的控制体积中的区域编号，图中控制体积中的箭头表示通过控制体积边界流量的正方向；
稳态水压膜力和扰动水膜压力
在节点pj，i处取值；
空泡体积分数αb、液相平均温度
液相平均速度(
和
)和空泡平均速度(
和
)均在节点Nj，i处取值；
p0为稳态水膜压力；
式(3)中：αl为液相体积分数；ρl为液相密度；t为时间；h0为轴承处于稳态时的水膜厚度；
式(3)中
和
的各流量系数如下式(I.1)所示：
式(I.1)中：kθ和kr为紊流系数；α′、β′和γ′为惯性系数；Mi为周向动量传递项；μl为液相动力粘度；
液相动态广义雷诺方程为：
式(4)中：pγ为扰动水膜压力；
和
为瞬时质量流量的扰动分量，均为复数，其可表示为：
式(5)中：υ为小扰动频率；
式(6)中：对于推力盘的下表面，α取+1，对于推力盘的上表面，α取-1；
式(5)中
和
的各流量系数如下式(I.2)示：
式(I.2)中：kθ和kr为紊流系数；α′、β′和γ′为惯性系数；Mi为周向动量传递项；μl为液相动力粘度；
s13：液相能量方程的建立：
考虑紊流效应、气-液界面能量传递和液固间的热量传递，按式(7)列写液相能量方程如下：
式(7)中：h为水膜厚度；
为通过膜厚的液相平均温度；Tp和Td分别为推力瓦和推力盘温度；cv为液相定容比热；Reh为当地雷诺数；kθ和kr分别为周向和径向紊流系数；Mi为周向动量传递项；Ei为两相之间的能量传递项；kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数；
s14：空泡体积分数输运方程的建立：
采用由压力驱动的空化模型，按式(8)列写空泡体积分数输运方程如下：
式(8)中：αb为空泡体积分数；
和
分别为空泡相的周向和径向平均速度；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所示；
其中，轴承承载力表示如下：
式(10)中：下标“lower”和“upper”分别表示推力盘的下表面和上表面；Wz为轴向承载力；Mx、My分别为绕x轴和y轴的抗倾覆力矩。
轴承液相沿径向的质量流量表示如下：
式(11)中：r0和ri分别为轴承的外圆半径和内圆半径；
轴承摩擦功耗表示如下：
将扰动压力
在推力盘的上、下表面积分，可得由推力轴承刚度和阻尼系数组成的复阻抗表示如下：
s2：闭式高速水润滑动压螺旋槽推力轴承静动态性能分析计算：
闭式高速水润滑动压螺旋槽推力轴承的几何和工况参数列于表1；
表1
在确定某一设计变量对轴承静动态性能的影响时，其他变量均保持不变。
基于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，θ)和计算平面(ξ，η)之间的转换关系为：
式(1)中：β为所述螺旋槽的螺旋角；rs为密封半径；r为轴承的半径方向坐标；θ为圆周方向坐标；ξ为螺旋槽方向坐标；η为圆周方向坐标；
s12：液相静动态广义雷诺方程的建立：
液相静态广义雷诺方程为：
式(2)中：下标中“1”、“2”、“3”、“4”表示控制体积中的区域编号；
p0为稳态水膜压力；
和
为瞬时质量流量的稳态分量，其可表示为：
式(3)中：“A”、“B”、“C”、“D”、“E”代表流量系数；
αl为液相体积分数；ρl为液相密度；t为时间；h0为轴承处于稳态时的水膜厚度；
液相动态广义雷诺方程为：
式(4)中：pγ为扰动水膜压力；
和
为瞬时质量流量的扰动分量，均为复数，其可表示为：
式(5)中：“A”、“B”、“C”、“D”、“E”代表流量系数；υ为小扰动频率；
为扰动膜厚，其可表示为：
式(6)中：对于推力盘的下表面，α取+1，对于推力盘的上表面，α取-1；
s13：液相能量方程的建立：
式(7)中：h为水膜厚度；
为通过膜厚的液相平均温度；Tp和Td分别为推力瓦和推力盘温度；cv为液相定容比热；μl为液相动力粘度；Reh为当地雷诺数；kθ和kr分别为周向和径向紊流系数；Mi为周向动量传递项；Ei为两相之间的能量传递项；kp和kd分别为液相与推力瓦、液相与推力盘之间的传热系数；
s14：空泡体积分数输运方程的建立：
式(8)中：αb为空泡体积分数；
和
分别为空泡相的周向和径向平均速度；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所述的闭式高速水润滑动压螺旋槽推力轴承静动态性能设计方法，其特征在于，轴承静态性能包括承载力、抗倾覆力矩、体积流量和摩擦功耗，轴承动态性能包括刚度和阻尼系数；所述设计变量包括倾斜角、转速、偏心率、半径比、螺旋角、径向槽长比和槽深。
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)
| 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)
| 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自由度耦合刚性转子系统动态设计方法 |
-
2020
- 2020-10-26 CN CN202011159296.3A patent/CN112182808B/zh active Active
Patent Citations (4)
| 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)
| 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)
| 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 |