液体动静压轴承油膜的压力场和温度场分析

液体动静压轴承油膜的压力场和温度场分析

刘蕾;刘保国;王攀;申会鹏;丁浩;郑金勇

【摘要】针对深浅腔液体动静压轴承的承载特性等问题,对液体动静压轴承的油膜压力场和温度场进行了仿真分析.以超高速磨削电主轴系统中常用的深浅腔液体动静压轴承为研究对象,建立了液体动静压轴承油膜的三维有限元模型,对油膜进行了网格划分,并对划分后的网格进行了质量评定;采用动网格技术实现了对油膜偏心率的变更,在不同主轴转速、偏心率的工作条件下,计算了深浅腔动静压轴承油膜压力和温度的分布情况,分析了其油膜压力分布和温度分布的变化规律;研究了转速、偏心率对动静压轴承的承载力和油膜温升的影响规律.研究结果表明:在深浅腔液体动静压轴承运转过程中,随着转速和偏心率的提高,油膜承载力和温升也随之提高,且转速对油膜温升的影响要比偏心率大.

【期刊名称】《机电工程》

【年(卷),期】2019(036)009

【总页数】7页(P900-906)

【关键词】油膜压力场;温度场;转速;偏心率;计算流体动力学

【作者】刘蕾;刘保国;王攀;申会鹏;丁浩;郑金勇

【作者单位】河南工业大学机电工程学院,河南郑州450001;河南工业大学河南省超硬磨料磨削装备重点实验室,河南郑州450001;河南工业大学机电工程学院,河南郑州450001;河南工业大学河南省超硬磨料磨削装备重点实验室,河南郑州450001;河南工业大学机电工程学院,河南郑州450001;河南工业大学机电工程学院,河南郑州450001;河南工业大学河南省超硬磨料磨削装备重点实验室,河南郑

州450001;河南工业大学机电工程学院,河南郑州450001;河南工业大学河南省超硬磨料磨削装备重点实验室,河南郑州450001;河南工业大学机电工程学院,河南郑州450001;河南工业大学河南省超硬磨料磨削装备重点实验室,河南郑州450001

【正文语种】中文

【中图分类】TH133.3

0 引言

超高速磨削已经广泛应用于汽车制造、航空航天等重要领域，成为磨削工艺未来的重要方向，电主轴是超高速磨削机床的重要组成部分，轴承直接影响电主轴寿命和超高速磨床的加工精度和可靠性[1]。在超高速磨削电主轴中，液体动静压轴承由于其回转精度高、动态刚性好和承载性能优等特点，具有广泛应用前景。超高速磨削电主轴系统中常用的动静压轴承由于经常在超高速状态下运转，油膜内部压力、温升等变化使得油膜稳定性下降，进而影响轴承的承载特性，严重时甚至会影响超高速磨削电主轴系统的稳定性。

针对液体动静压轴承的承载特性问题，国内外学者大多以雷诺方程为基础建立轴承的数学模型，从而对液体动静压轴承的承载特性进行了研究，而雷诺方程是N—S 方程的简化形式，忽略了惯性项、油膜曲率、径向流场变化等因素，对于液体动静压轴承的承载特性研究具有一定的局限性[2-5]。随着计算机技术的发展，基于N—S方程的CFD方法在动静压轴承的承载特性计算中得到了广泛应用。郭胜安等[6]运用CFD软件对轴承性能进行了研究，得到了轴承结构参数和工作参数对刚度、流量和温升的影响机制；刘豪杰[7-8]运用Fluent软件对数控铣床上动静压轴承的压力场和温度场进行了分析，得到了不同工况下油膜压力场和温度场的变化规

律，并搭建了实验平台以验证分析结果；陈超等[9]基于CFD软件对水润滑动静压轴承的承载力进行了研究，发现了水润滑动静压轴承适用于高速轻载的工况；王攀[10]采用基于6DOF模型及动网格的计算方法，对动静压轴承的油膜刚度和阻尼

进行了计算，得到了转速对油膜刚度和阻尼的影响规律，这种计算方法可以反映在外载荷作用下油膜的运动变化过程，更加准确地计算轴颈静平衡位置，使得计算结果与实际工况更为吻合。

本文将在前期研究基础上，以超高速磨削电主轴液体动静压轴承为研究对象，采用动网格技术通过数值计算，得到不同转速和偏心率下轴承的压力分布和温度分布情况，进而研究转速、偏心率对油膜承载力和温升的影响机制，为后续优化设计做准备。

1 CFD控制方程及边界条件

计算流体动力学(CFD)以计算机为载体，通过离散求解流动方程[11-12]，从而进

行数值分析，以解决工程中遇到的各种流体问题。

1.1 质量守恒方程

即流体连续方程。在流体流动的空间中，任取一六面体空间，将该六面体空间称为控制体，根据质量守恒定律，控制体内质量的减少率等于流出控制体的质量净流量。其表达式为：

(1)

式中：ρ—密度；t—时间；ν—速度矢量。

1.2 动量守恒方程

即N—S方程。根据牛顿第二定律，作用在流体微团上力的总和等于微团的质量乘以微团的加速度。其表达式为：

(2)

(3)

式中：p—压力；F—外部体积力；τ—应力张量；I—单位张量。

1.3 能量守恒方程

即热力学第一定律。微元体中能量的变化率等于进入微元体的净热流量加上体积力与表面力对微元体所做的功率之和。其表达式为：

(4)

式中：Cp—比热容；K—热传动系数；ST—耗散项。

1.4 边界条件

采用Fluent进行求解时，油膜压力边界条件采用雷诺边界条件，设定进油孔为压力入口，压力值为供油压力，封油面间隙为出油口，压力值为大气压力，进油孔温度为环境温度。

2 计算模型的建立

2.1 物理模型及网格划分

本研究以超高速磨削电主轴系统常用的深浅腔动静压轴承为研究对象，其结构参数如表1所示。

表1 结构参数模型参数数值模型参数数值轴承直径/mm80油腔轴向宽度/mm64轴承宽度/mm80浅腔包角/(°)46油膜厚度/mm0.025深腔包角/(°)8浅油腔深度/mm0.050节流小孔直径/mm0.6深油腔深度/mm0.100腔体个数4

本研究将轴承的实体模型导入到ICEM软件中，建立轴承油膜的三维有限元模型，并进行网格划分。划分过程中，考虑到网格生成质量和后续计算速度等问题，采用结构化网格划分。同时，考虑到油膜的腔体形状和进油孔的存在，在分块时采用